[论文解读] Projective and telescopic projective integration for the nonlinear BGK and Boltzmann equations

Ward Melis, Thomas Rey+1 more|arXiv (Cornell University)|Dec 18, 2017

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

一句话总结

本文提出了一种针对非线性BGK方程和玻尔兹曼方程的高阶、完全显式的望远镜投影积分（TPI）格式，利用内层时间步长来抑制刚性分量，外层Runge-Kutta阶段来推进解。该方法的计算成本几乎与刚性程度无关，2D模拟中BGK方程的加速比超过100，玻尔兹曼方程的加速比达5.9，这是由于投影层级随刚性程度对数缩放所致。

ABSTRACT

We present high-order, fully explicit projective integration schemes for nonlinear collisional kinetic equations such as the BGK and Boltzmann equation. The methods first take a few small (inner) steps with a simple, explicit method (such as direct forward Euler) to damp out the stiff components of the solution. Then, the time derivative is estimated and used in an (outer) Runge-Kutta method of arbitrary order. The procedure can be recursively repeated on a hierarchy of projective levels to construct telescopic projective integration methods. Based on the spectrum of the linearized collision operator, we deduce that the computational cost of the method is essentially independent of the stiffness of the problem: with an appropriate choice of inner step size, the time step restriction on the outer time step, as well as the number of inner time steps, is independent of the stiffness of the (collisional) source term. In some cases, the number of levels in the telescopic hierarchy depends logarithmically on the stiffness. We illustrate the method with numerical results in one and two spatial dimensions.

研究动机与目标

• 解决由于平均自由程较短而引起的非线性动力学方程（如BGK和玻尔兹曼方程）中的刚性问题。
• 开发完全显式的时间积分格式，在避免隐式求解器的同时保持高阶精度。
• 克服在多维相空间中刚性碰撞算子带来的计算瓶颈。
• 即使在刚性程度高（ε较小）的流体区域中，也能保证稳定性和效率。
• 通过分层投影时间推进实现可扩展、自适应的时间积分。

提出的方法

• 采用分层投影层级：内层使用显式欧拉法来抑制刚性分量，外层使用高阶Runge-Kutta法来推进解。
• 从内层步骤中估计时间导数，以构建显式外层时间积分，从而实现高阶精度。
• 通过在多级之间递归应用投影步骤，实现望远镜投影积分（TPI），以减少对外层时间步长的限制。
• 利用线性化碰撞算子的谱理论来确保稳定性，并指导时间步长的选择。
• 使用快速谱方法计算碰撞算子，采用WENO2空间离散化以精确捕捉强梯度。
• 选择内层时间步长 h₀ = ε，以实现外层时间步长与刚性程度的解耦，确保计算成本与ε无关。

实验结果

研究问题

• RQ1 能否在不使用隐式求解器的情况下，为刚性非线性BGK和玻尔兹曼方程实现完全显式、高阶的时间积分？
• RQ2 在望远镜投影积分中，投影层级的数量如何随刚性程度（ε）变化？
• RQ3 与标准显式格式相比，TPI方法的计算成本是多少，且是否保持与刚性程度无关？
• RQ4 TPI格式能否在具有强梯度和多尺度动力学的问题（如激波管）中保持精度？
• RQ5 在刚性区域中，内层时间步长和外层时间步长的选择如何影响稳定性与效率？

主要发现

• 望远镜投影积分方法的计算成本几乎与刚性程度ε无关，这是由于采用了最优的内层时间步长 h₀ = ε。
• 在BGK方程中实现了超过100倍的加速比，在2D玻尔兹曼方程（伪麦克斯韦粒子）中实现了5.9倍的加速比。
• 所需投影层级的数量随刚性程度呈对数缩放，使得该方法在ε极小时依然高效。
• 通过在外层使用PRK4、内层最底层使用FE作为积分器，并结合WENO2空间离散化，成功保持了高阶精度，其中WENO2对捕捉涡旋结构至关重要。
• 即使在流体区域（ε = 5·10⁻⁵）中，该方法仍保持稳定和准确，而标准显式格式会因CFL限制而失效。
• 该框架可自然推广至自适应、空间依赖的投影层级结构，为未来进一步提升效率提供了可能。