[论文解读] High order semi-Lagrangian methods for the BGK equation

Maria Groppi, Giovanni Russo+1 more|arXiv (Cornell University)|Nov 28, 2014

Gas Dynamics and Kinetic Theory参考文献 15被引用 4

一句话总结

本论文提出了一种基于BDF和DIRK时间积分格式的BGK方程高阶半拉格朗日格式，沿特征线进行时间推进，实现了在稀薄流与流体流动两种状态下的渐近保持性与高精度。数值实验表明，特别是不采用插值的BDF3方法，在计算效率和鲁棒性方面表现优异，相较于DIRK格式在成本-性能权衡上更具优势。

ABSTRACT

A new class of high-order accuracy numerical methods for the BGK model of the Boltzmann equation is presented. The schemes are based on a semi-lagrangian formulation of the BGK equation; time integration is dealt with DIRK (Diagonally Implicit Runge Kutta) and BDF methods; the latter turn out to be accurate and computationally less expensive than the former. Numerical results and examples show that the schemes are reliable and efficient for the investigation of both rarefied and fluid regimes in gasdynamics.

研究动机与目标

• 开发适用于BGK方程的高阶精度数值格式，确保在稀薄流与流体流动两种状态下均保持稳定与高精度。
• 比较半拉格朗日格式下BDF与DIRK时间积分方法在BGK方程求解中的性能表现。
• 研究空间插值（尤其是速度空间中）对精度与计算成本的影响。
• 设计渐近保持格式，确保在流体极限（ε → 0）下保持精度，且无需解析克努森层。
• 通过楚约化方法将方法扩展至一维空间与三维速度空间，以适用于实际气体动力学应用。

提出的方法

• 采用半拉格朗日方法沿特征线重写BGK方程，将时间演化与空间网格约束解耦。
• 时间积分采用BDF（后向差分公式）与DIRK（对角隐式龙格-库塔）格式，其中BDF3展现出更优的计算效率。
• 通过WENO重构（WENO23与WENO35）实现高阶空间精度，适用于非均匀分布函数，有效抑制间断附近的振荡。
• 提出一种避免插值的格式变体，直接在相空间网格上运算，显著降低速度空间中的计算开销。
• 通过楚约化方法将方法扩展至一维空间与三维速度空间，保持精度的同时降低维度。
• 该格式保证了质量、动量与能量的守恒，并满足熵耗散特性，维持物理一致性。

实验结果

研究问题

• RQ1 高阶BDF与DIRK格式能否与半拉格朗日方法有效结合，以实现BGK方程的高精度求解？
• RQ2 时间积分格式的选择（BDF与DIRK）在稀薄流与流体流动两种状态下如何影响计算成本与精度？
• RQ3 在速度空间中避免空间插值在多大程度上可提升计算效率而不损失精度？
• RQ4 所提出的格式在流体极限（ε → 0）下是否保持渐近保持特性，尤其在低克努森数条件下？
• RQ5 WENO重构技术在多大程度上提升了分布函数间断附近稳定性与分辨率？

主要发现

• 不采用插值的BDF3方法在所有测试案例中均展现出最佳的整体性能，兼顾精度与计算成本。
• BDF格式在保持相近精度的同时，计算成本显著低于DIRK格式。
• 采用WENO重构的格式实现了高阶精度，并有效抑制了激波与间断附近的虚假振荡。
• 该方法保持了渐近保持特性，即使在ε极小的情况下也能正确捕捉流体极限（可压缩纳维-斯托克斯行为）。
• 在速度空间中避免插值显著降低了计算成本，尤其在速度网格较粗的问题中优势明显。
• 数值结果验证了该格式在高分辨率模拟稀薄流与近平衡态流体动力学中的鲁棒性。