[论文解读] Exponential Runge-Kutta schemes for inhomogeneous Boltzmann equations with high order of accuracy

Qin Li, Lorenzo Pareschi|arXiv (Cornell University)|Aug 13, 2012

Gas Dynamics and Kinetic Theory参考文献 26被引用 9

一句话总结

本文提出用于刚性非均匀玻尔兹曼方程的高阶指数 Runge-Kutta (ExpRK) 格式，采用时变麦克斯韦平衡态以实现渐近保持性与一致精度。该方法避免了非线性求解，可在所有克努森数下保持高阶精度，并通过 Shu-Osher 表示法保证非负性与强稳定性，数值实验验证了其在三阶精度范围内的有效性。

ABSTRACT

We consider the development of exponential methods for the robust time discretization of space inhomogeneous Boltzmann equations in stiff regimes. Compared to the space homogeneous case, or more in general to the case of splitting based methods, studied in Dimarco Pareschi (SIAM J. Num. Anal. 2011) a major difficulty is that the local Maxwellian equilibrium state is not constant in a time step and thus needs a proper numerical treatment. We show how to derive asymptotic preserving (AP) schemes of arbitrary order and in particular using the Shu-Osher representation of Runge-Kutta methods we explore the monotonicity properties of such schemes, like strong stability preserving (SSP) and positivity preserving. Several numerical results confirm our analysis.

研究动机与目标

• 开发适用于刚性非均匀玻尔兹曼方程的高阶时间格式，使其在所有克努森数范围内保持稳定与高精度。
• 通过避免非线性系统求解，克服基于分裂法与隐式格式的局限性，同时维持高阶精度。
• 利用指数积分格式，在刚性情形下实现渐近保持性（AP）与单调性（如非负性与强稳定性）。
• 将高阶指数方法从均匀情形推广至非均匀玻尔兹曼方程，无需域分解或分裂技术。
• 为带有源项的 kinetic 方程提供一种鲁棒、一致精度且保持非负性的数值框架。

提出的方法

• 提出两类指数 Runge-Kutta 格式：ExpRK-F（固定平衡态）与 ExpRK-V（时变平衡态），用于非均匀玻尔兹曼方程。
• 采用 Runge-Kutta 格式的 Shu-Osher 表示法，以保证强稳定性保持性（SSP）与非负性保持性。
• 通过将玻尔兹曼方程重新表述为指数积分格式，结合碰撞项的指数时间积分，推导出数值格式。
• 利用矩匹配与渐近分析，计算时变麦克斯韦平衡态 $M$ 及其时间导数 $\partial_t M$。
• 采用 $\varepsilon$-正则化技术处理刚性碰撞项，确保在流体极限下的统一精度。
• 通过矩保持性与压缩性论证验证方法，证明在每一阶段 $f^{(i)}$ 与 $M^{(i)}$ 具有相同的矩。

实验结果

研究问题

• RQ1 能否为非均匀玻尔兹曼方程构造出在所有克努森数下保持一致高阶精度的高阶指数 Runge-Kutta 格式？
• RQ2 在刚性 kinetic 方程中，如何通过指数积分格式实现渐近保持性与强稳定性（SSP）？
• RQ3 在 ExpRK 格式中，使用固定平衡态与时变平衡态对精度与单调性有何影响？
• RQ4 能否在高阶指数积分格式中保持玻尔兹曼方程分布函数的非负性？
• RQ5 在刚性区域中，所提格式与基于分裂法或 IMEX 方法相比，在精度与稳定性方面表现如何？

主要发现

• 所提出的 ExpRK-V 格式通过采用时变平衡态，在所有克努森数范围内实现了统一的高阶精度，尤其在中间区域优于固定平衡态格式。
• Shu-Osher 表示法实现了强稳定性保持性（SSP）与非负性保持性，确保了物理解释的合理性。
• 数值实验验证了 ExpRK-F 与 ExpRK-V 格式在时间方向均达到三阶收敛率，即使在刚性区域亦成立。
• 在 $\varepsilon \to 0$ 极限下，该方法正确保持了流体极限（Euler 方程），证实了渐近保持性。
• 分布函数 $f$ 与平衡态 $M$ 之间的差异被控制在 $O(\varepsilon)$ 范围内，确保了矩方程中校正项的小幅且稳定。
• 对于截断的麦克斯韦分子，证明了 $d_2$-范数的压缩性性质 $\mathrm{d}_2(P(f), P(g)) \leq \mathrm{d}_2(f,g)$，支持了傅里叶域中的稳定性。