[论文解读] Onsager reciprocity principle for kinetic models and kinetic schemes

A.K. Mahendra, Ramesh Singh|arXiv (Cornell University)|Aug 19, 2013

Gas Dynamics and Kinetic Theory参考文献 47被引用 8

一句话总结

该论文提出了一种基于昂萨格互易原理和非平衡态热力学的新型动力学模型，以确保动力学模拟中传输系数和熵产生率的正确性。通过引入碰撞概率函数 $\mathbb{P}(\bar{y})$，该模型实现了对稀薄流与连续流的统一单体求解器，其动力学壁面边界条件可准确捕捉跨尺度的滑移流，而无需昂贵的多尺度耦合计算。

ABSTRACT

Boltzmann equation requires some alternative simpler kinetic model like BGK to replace the collision term. Such a kinetic model which replaces the Boltzmann collision integral should preserve the basic properties and characteristics of the Boltzmann equation and comply with the requirements of non equilibrium thermodynamics. Most of the research in development of kinetic theory based methods have focused more on entropy conditions, stability and ignored the crucial aspect of non equilibrium thermodynamics. The paper presents a new kinetic model formulated based on the principles of non equilibrium thermodynamics. The new kinetic model yields correct transport coefficients and satisfies Onsager's reciprocity relationship. The present work also describes a novel kinetic particle method and gas kinetic scheme based on this linkage of non-equilibrium thermodynamics and kinetic theory. The work also presents derivation of kinetic theory based wall boundary condition which complies with the principles of non-equilibrium thermodynamics, and can simulate both continuum and rarefied slip flow in order to avoid extremely costly multi-scale simulation.

研究动机与目标

• 开发一种动力学模型，使其在保持玻尔兹曼方程基本性质的同时，符合非平衡态热力学，特别是昂萨格互易关系与熵产生原理。
• 解决现有动力学格式中普遍存在的热力学不一致性问题，这些格式往往优先考虑熵产生，而非正确分配熵在共轭热力学力之间的分布。
• 构建一种统一的动力学粒子方法与气体动力学格式，使其在从稀薄流到连续流的全流态范围内无缝运行，通过将非平衡态热力学嵌入碰撞模型实现。
• 推导一种基于动力学理论的壁面边界条件，使其满足非平衡态热力学要求，并能准确模拟连续流与稀薄流中的滑移流，而无需多尺度建模。
• 在涵盖连续流与稀薄流状态的测试案例上验证所提框架的准确性与鲁棒性。

提出的方法

• 基于昂萨格变分原理构建新动力学模型，确保传输系数的对称性并满足昂萨格互易关系。
• 引入碰撞概率函数 $\mathbb{P}(\bar{y})$，将其作为Knudsen层中法向距离 $\bar{y} = y/\lambda_e$ 的函数，调节弛豫时间 $\tilde{t}_{R(j)} = \mathbb{P}(\bar{y})t_{R(j)}$。
• 基于新模型推导动力学通量矢量分裂与粒子方法，其中粘性通量由 $\mathbb{P}(\bar{y})$ 缩放，以反映非平衡效应。
• 基于Maxwell气体相互作用模型构建壁面边界条件，将入射、反射与漫反射分布按accommodation系数 $\sigma$ 加权，且非平衡部分在壁面处消失（$\mathbb{P}(0) = 0$）。
• 应用扩展不可逆热力学（EIT）框架，对Knudsen层中非线性通量-力关系进行建模，通过 $\boldsymbol{J}_i \approx \mathbb{P}(\bar{y}) \sum_j \boldsymbol{L}_{ij} \odot \boldsymbol{X}_j$ 近似EIT通量。
• 在Knudsen层中采用Chapman-Enskog分布，将非平衡部分表示为 $f_1(\bar{y}) = f_0 - \mathbb{P}(\bar{y}) \sum_j \boldsymbol{\Upsilon}_j \odot \boldsymbol{X}_j$，确保热力学一致性。

实验结果

研究问题

• RQ1 如何构建一种动力学模型，使其满足昂萨格互易原理，同时正确保持如普朗特数等传输系数？
• RQ2 能否设计一种单一动力学格式，通过将非平衡态热力学嵌入碰撞模型，实现对稀薄流与连续流的精确模拟？
• RQ3 碰撞概率函数 $\mathbb{P}(\bar{y})$ 在实现Knudsen层中粘性通量的统一处理中起到何种作用？
• RQ4 如何推导一种满足非平衡态热力学的动力学壁面边界条件，以正确模拟稀薄流与连续流状态下的滑移流？
• RQ5 所提模型在多大程度上能确保熵产生率在不同热力学力（如应力与温度梯度）之间的正确分配？

主要发现

• 所提动力学模型满足昂萨格互易关系，并正确生成普朗特数，确保了传输性质的热力学一致性。
• 动力学格式中的粘性通量分量由碰撞概率函数 $\mathbb{P}(\bar{y})$ 缩放，从而能够准确捕捉Knudsen层中的非平衡效应。
• 基于Maxwell气体相互作用模型推导的壁面边界条件能正确反映镜面反射与漫反射行为，且由于 $\mathbb{P}(0) = 0$，非平衡部分在壁面处消失。
• Knudsen层中的状态更新方程引入了依赖于 $\mathbb{P}(\bar{y})$ 的弛豫时间，使模型能平滑过渡至壁面附近的平衡行为。
• 基于新模型的动力学粒子方法与通量矢量分裂格式成功模拟了连续流与稀薄流滑移流测试案例，验证了该框架的鲁棒性与准确性。
• 该模型通过将非平衡态热力学直接嵌入动力学框架，避免了昂贵的多尺度模拟，实现了全流态范围内的单体求解器。