[论文解读] Flexible Macroscopic Models for Dense-Fluid Shockwaves: Partitioning Heat and Work; Delaying Stress and Heat Flux; Two-Temperature Thermal Relaxation

Wm. G. Hoover, Carol G. Hoover+1 more|arXiv (Cornell University)|May 10, 2010

Gas Dynamics and Kinetic Theory参考文献 21被引用 7

一句话总结

该论文提出了一种用于致密流体激波的灵活宏观模型，区分纵向温度与横向温度，引入应力、热流和热能的三种弛豫时间，并对功和热量贡献进行分解。该模型通过带有麦克斯韦-卡塔内奥型和克鲁克型弛豫的双曲偏微分方程系统，再现了微观模拟中的关键特征——如应力与热流的滞后、张量温度以及双温热弛豫，实现了数值模拟中稳定且物理解释一致的激波特性的模拟。

ABSTRACT

Macroscopic models which distinguish the longitudinal and transverse temperatures can provide improved descriptions of the microscopic shock structures as revealed by molecular dynamics simulations. Additionally, we can include three relaxation times in the models, two based on Maxwell's viscoelasticity and its Cattaneo-equation analog for heat flow, and a third thermal, based on the Krook-Boltzmann equation. This approach can replicate the observed lags of stress (which lags behind the strain rate) and heat flux (which lags behind the temperature gradient), as well as the eventual equilibration of the two temperatures. For profile stability the time lags cannot be too large. By partitioning the longitudinal and transverse contributions of work and heat and including a tensor heat conductivity and bulk viscosity, all the qualitative microscopic features of strong simple-fluid shockwave structures can be reproduced.

研究动机与目标

• 开发一种宏观模型，以捕捉分子动力学模拟中观察到的致密流体激波的非平衡行为。
• 通过在应力和热流中引入非局部、延迟响应，克服经典纳维-斯托克斯方程和 Burnett 方程的局限性。
• 在强激波中建模纵向与横向温度不同的物理现实，与微观观测一致。
• 通过引入三种弛豫时间（应力与热流的麦克斯韦-卡塔内奥型模型，以及热弛豫的克鲁克-玻尔兹曼方程）确保数值稳定性和物理解释一致性。

提出的方法

• 为质量、动量和能量建立五个双曲偏微分方程系统，将应力和热流作为独立变量。
• 引入三种弛豫时间：τσ 用于应力（麦克斯韦型粘弹性），τq 用于热流（卡塔内奥方程类比），τT 用于热弛豫（克鲁克-玻尔兹曼模型）。
• 通过张量热导率和体胀粘度，将功和热量贡献划分为纵向与横向分量。
• 使用二阶空间精度的交错有限差分格式与四阶龙格-库塔时间积分方法求解系统。
• 通过 Lucy 型加权函数实现空间平均，以处理基于粒子的平均中的非局部性，确保与原子模拟的一致性。
• 基于 Hugoniot 关系施加边界条件，两端分别为冷流体和热流体状态，并确保边界处热流为零且剪切应力为零。

实验结果

研究问题

• RQ1 宏观模型如何捕捉致密流体激波中应力与热流相对于应变率和温度梯度的观测滞后？
• RQ2 张量温度（纵向与横向）在超越标量温度近似时，对准确描述激波特性的角色是什么？
• RQ3 一个包含三种不同弛豫时间的统一模型能否再现强激波分子动力学模拟中观察到的全部定性特征？
• RQ4 将功和热量划分为纵向与横向分量如何影响激波特性的稳定性和形状？
• RQ5 弛豫时间的稳定性约束是什么？它们如何限制物理有效解的参数空间？

主要发现

• 该模型成功再现了应力与热流的观测滞后，延迟量级与粒子碰撞时间一致，与分子动力学模拟一致。
• 应力与热流被证明滞后于其驱动力梯度（应变率与温度梯度），证实了强激波中因果关系的弛豫行为。
• 纵向与横向温度在激波前缘演化不同，横向温度更早上升，两者的平衡仅在弛豫后实现。
• 当弛豫时间不过大时，数值解才保持稳定，表明非平衡状态下流体记忆存在物理极限。
• 引入张量热导率与体胀粘度，可准确再现激波的完整定性结构，包括非平衡温度各向异性。
• 该模型在广泛参数范围内实现稳定、收敛的解，数值稳定性通过自适应时间步长和弛豫参数的人工时间依赖性进一步增强。