[论文解读] Globally Hyperbolic Regularization of Grad's Moment System in One Dimensional Space

Zhenning Cai, Yuwei Fan+1 more|arXiv (Cornell University)|Nov 15, 2011

Gas Dynamics and Kinetic Theory参考文献 13被引用 10

一句话总结

本文通过修改最高阶矩方程以消除非双曲项，提出了一维Grad矩系统的全局双曲正则化方法，确保特征速度为实数并实现稳定的波动力学。该方法不依赖碰撞项即可实现全局双曲性，保留了关键波结构（稀疏波、接触间断、激波），并随着矩阶数增加实现收敛。

ABSTRACT

In this paper, we present a regularization to 1D Grad's moment system to achieve global hyperbolicity. The regularization is based on the observation that the characteristic polynomial of the Jacobian of the flux in Grad's moment system is independent of the intermediate moments. The method is not relied on the form of the collision at all, thus this regularization is applicable to the system without collision terms. Moreover, the proposed approach is proved to be the unique one if only the last moment equation is allowed to be alternated to match the condition of non-equilibrium independent characteristic speeds. The hyperbolic structure of the regularized system, including the signal speeds, Riemann invariants and the properties of the characteristic waves including the rarefaction wave, contact discontinuity and shock are provided in the perfect formations.

研究动机与目标

• 解决Grad矩系统缺乏全局双曲性的问题，该问题在远离平衡态时会导致数值发散。
• 开发一种正则化方法，确保系统在所有状态（而不仅在平衡附近）均保持双曲性。
• 保留特征波结构（稀疏波、接触间断、激波）的物理特性，与Euler方程相似。
• 提供一种不依赖碰撞项的正则化方法，适用于无碰撞和有碰撞两种情形。
• 确定对最后一阶矩方程的唯一修改方式，使其与高斯-埃尔米特插值点匹配并确保双曲性。

提出的方法

• 识别通量雅可比矩阵的特征多项式仅依赖于宏观速度、温度以及两个高阶埃尔米特系数。
• 观察到这两个系数引入了复特征值，破坏了双曲性。
• 仅修改最高阶矩方程，以消除特征多项式中的问题项。
• 在特征速度与高斯-埃尔米特求积点匹配的约束下，证明该修改的唯一性。
• 推导出信号传播速度、黎曼不变量以及波特性（真正非线性或线性退化场）的精确表达式。
• 通过分析稀疏波、接触间断和激波，验证了系统的双曲结构，表明其与Euler方程具有结构相似性。

实验结果

研究问题

• RQ1 能否在不依赖碰撞项的前提下，为一维Grad矩系统构造出全局双曲正则化？
• RQ2 对最高阶矩方程的唯一修改方式是什么，可确保特征速度为实数并匹配高斯-埃尔米特点？
• RQ3 正则化系统中特征波结构（稀疏波、接触间断、激波）与Euler方程中的结构相比如何？
• RQ4 增加矩的阶数是否会导致激波管问题中向动能解的收敛？
• RQ5 正则化系统在远离平衡态的状态下是否仍能保持双曲性和稳定性？

主要发现

• 通量雅可比矩阵的特征多项式仅依赖于速度、温度以及两个高阶埃尔米特系数，简化了双曲性分析。
• 原系统非双曲性的根源在于两个最高阶系数引入的复特征值。
• 所提出的正则化方法唯一地修改了最后一阶矩方程，将这些系数从特征多项式中消除，从而确保特征值为实数。
• 正则化系统中所有特征场均为真正非线性或线性退化，支持适定的黎曼问题。
• 系统支持所有三种基本波型——稀疏波、接触间断和激波，其动力学结构与Euler方程相似。
• 激波管问题的数值结果表明，随着矩阶数增加（M=2至M=10），系统实现收敛，证实了精度和稳定性的提升。