[论文解读] A New Class of Conservative Large Time Step Methods for the BGK Models of the Boltzmann Equation

Pietro Santagati, Giovanni Russo|arXiv (Cornell University)|Mar 27, 2011

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

一句话总结

本文提出了一类新型的保守性、高阶、大时间步长的半拉格朗日方法，用于玻尔兹曼方程的BGK模型。通过将拉格朗日框架与对松弛项的对角隐式龙格-库塔（DIRK）时间积分相结合，该方法消除了CFL稳定性限制，从而在稀薄气体和流体区域中实现了稳定且精确的模拟，兼具高效率和极小的守恒误差。

ABSTRACT

This work is aimed to develop a new class of methods for the BGK model of the Boltzmann equation. This technique allows to get high order of accuracy both in space and time, theoretically without CFL stability limitation. It's based on a Lagrangian formulation of the problem: information is stored on a fixed grid in space and velocity, and the equation is integrated along the characteristics. The source term is treated implicitly by using a DIRK (Diagonally Implicit Runge Kutta) scheme in order to avoid the time step restriction due to the stiff relaxation. In particular some L-stable schemes are tested by smooth and Riemann problems, both in rarefied and fully fluid regimes. Numerical results show good accuracy and efficiency of the method.

研究动机与目标

• 开发一种用于玻尔兹曼方程BGK模型的数值方法，以克服传统方法中常见的CFL稳定性限制。
• 在保持大时间步长稳定性的同时，实现空间和时间上的高阶精度。
• 在离散格式中确保质量、动量和能量的守恒，尤其在流体动力学极限下至关重要。
• 实现对从光滑到不连续流动（如激波）的广泛Knudsen数范围内的精确模拟。
• 为未来在MEMS和微流体等实际问题中的多维扩展提供基础。

提出的方法

• 在半拉格朗日框架中表述BGK方程，沿特征线积分以避免CFL条件限制。
• 采用高阶WENO重构方法，从特征轨迹上精确重构网格点处的分布函数。
• 应用对角隐式龙格-库塔（DIRK）格式对刚性松弛项进行隐式处理，确保在小松弛时间下的稳定性。
• 在物理空间和速度空间中使用固定的结构化笛卡尔网格，并在网格节点处存储分布函数。
• 通过投影或校正步骤强制实现离散时刻的守恒（密度、动量、能量），从而实现完全保守的变体。
• 采用L-稳定DIRK格式，即使在流体极限（小Knudsen数）下也能保持精度与稳定性。

实验结果

研究问题

• RQ1 是否可通过隐式处理BGK松弛项的半拉格朗日方法，在无CFL限制下实现高阶精度与稳定性？
• RQ2 该方法在稀薄与流体动力学区域中的表现如何，特别是在激波等不连续解情况下的表现？
• RQ3 非保守变体中的守恒误差积累程度如何？在保守形式中能否有效抑制这些误差？
• RQ4 在不同时间步长和CFL数下，该方法的收敛行为与精度如何，尤其是在大CFL值下？
• RQ5 当松弛时间τ变得极小时（趋近于流体极限），该方法是否仍能保持高精度与稳定性？

主要发现

• 对于光滑流动，该方法在空间和时间上均达到二阶精度，网格细化下密度、速度和温度的观测收敛阶数约为2.7–2.8。
• 在Riemann问题测试中，即使在大CFL数（高达10.5）下，该方法仍能保持对激波结构的高分辨率，且当网格分辨率加倍时，相对L2误差降低一个数量级。
• 非保守版本表现出极小的守恒误差：当τ = 10⁻²时，密度、动量和能量的相对误差在Nx = 1280下仍低于10⁻⁴，且随网格细化而减小。
• 保守版本在τ = 10⁻⁶且Nx = 1280时，成功将质量、动量和能量守恒至机器精度，误差低于10⁻⁸。
• 熵演化保持单调且符合物理规律，证实了该方法在稀薄与流体区域中均满足热力学第二定律。
• 该方法在大CFL值（如CFL = 10.5）下仍保持稳定与精确，未观察到数值耗散或振荡，而传统显式格式则存在此类问题。