[论文解读] Beyond pressureless gas dynamics: Quadrature-based velocity moment models

Christophe Chalons, Damien Kah+1 more|arXiv (Cornell University)|Nov 12, 2010

Gas Dynamics and Kinetic Theory参考文献 18被引用 42

一句话总结

本文将基于积分的速率矩模型扩展至无限Knudsen数极限下的动力学方程，提出一种基于两节点积分的四矩模型，能够捕捉光滑解与奇异解。该研究建立了熵测度解的严格框架，通过动力学格式验证了方法的有效性，并展示了对粒子轨迹交叉及广义δ-激波的收敛性与合理处理能力。

ABSTRACT

Following the seminal work of F. Bouchut on zero pressure gas dynamics which has been extensively used for gas particle-flows, the present contribution investigates quadrature-based velocity moments models for kinetic equations in the framework of the infinite Knudsen number limit, that is, for dilute clouds of small particles where the collision or coalescence probability asymptotically approaches zero. Such models define a hierarchy based on the number of moments and associated quadrature nodes, the first level of which leads to pressureless gas dynamics. We focus in particular on the four moment model where the flux closure is provided by a two-node quadrature in the velocity phase space and provide the right framework for studying both smooth and singular solutions. The link with both the kinetic underlying equation as well as with zero pressure gas dynamics is provided and we define the notion of measure solutions as well as the mathematical structure of the resulting system of four PDEs. We exhibit a family of entropies and entropy fluxes and define the notion of entropic solution. We study the Riemann problem and provide a series of entropic solutions in particular cases. This leads to a rigorous link with the possibility of the system of macroscopic PDEs to allow particle trajectory crossing (PTC) in the framework of smooth solutions. Generalized $\\delta$-choc solutions resulting from Riemann problem are also investigated. Finally, using a kinetic scheme proposed in the literature without mathematical background in several areas, we validate such a numerical approach in the framework of both smooth and singular solutions.

研究动机与目标

• 将基于积分的矩模型扩展至压强为零的气体动力学之外，通过速度相空间中的两节点积分引入四矩系统。
• 在由动力学方程导出的弱双曲守恒律背景下，建立测度解与熵解的数学框架。
• 研究矩空间边界处的系统行为，特别是当速度合并时的情况，并确保数值可实现性。
• 通过动力学格式验证数值方法，展示其在捕捉光滑解与奇异解方面的收敛性与准确性。
• 为包含轨迹交叉与复杂相空间动力学的颗粒载流体建模提供理论与数值基础。

提出的方法

• 基于速度分布的两节点积分近似，推导出四矩系统，将分布表示为狄拉克δ函数之和。
• 将系统定义为弱双曲守恒律系统，并表征其数学结构，包括可实现性条件与矩空间的边界。
• 通过构造与系统动力学相容的一族熵与熵通量，引入熵测度解的概念。
• 应用基于积分节点与权重的动力学格式数值求解系统，确保可实现性并合理处理奇异性。
• 通过黎曼问题解与收敛性研究验证格式，采用结构化网格上的L₁误差范数。
• 通过监测q/e^{3/2}等量的变化，分析从矩空间内部到边界的过渡，并验证锥保持性。

实验结果

研究问题

• RQ1 如何将基于积分的矩模型扩展至压强为零的气体动力学之外，以捕捉稀疏颗粒流中的复杂相空间动力学？
• RQ2 四矩系统的数学结构是什么？它如何支持光滑解与奇异解（包括广义δ-激波）？
• RQ3 如何为该系统定义并表征熵测度解？它们在确保物理一致性方面发挥什么作用？
• RQ4 动力学格式能否保持可实现性，并准确模拟粒子轨迹交叉及在矩空间边界的过渡？
• RQ5 数值格式的收敛行为如何？其在捕捉黎曼问题解析解方面的表现如何？

主要发现

• 四矩基于积分的模型成功捕捉了光滑解与广义δ-激波，展示了其在模拟包含粒子轨迹交叉的复杂动力学方面的能力。
• 动力学格式在L₁范数下表现出约0.5的收敛率，这通过从400到3200个单元的网格细化过程中误差的降低得到验证。
• 该格式保持了可实现性锥，并在矩空间边界处维持了数值稳定性，表现为q/e^{3/2}的行为稳定且无虚假振荡。
• 该系统支持一族熵与熵通量，使得可定义与底层动力学方程一致的熵测度解。
• 在单流极限下，该模型退化为压强为零的气体动力学，确认了其与既有框架的一致性。
• 数值验证表明，该格式能准确复现黎曼问题的解析解，包括具有奇异性与轨迹交叉的情形。