[论文解读] Asymptotic preserving Implicit-Explicit Runge-Kutta methods for non linear kinetic equations

Giacomo Dimarco, Lorenzo Pareschi|arXiv (Cornell University)|May 4, 2012

Gas Dynamics and Kinetic Theory参考文献 30被引用 25

一句话总结

本文提出了一种用于刚性非线性 kinetic 方程（特别是 Boltzmann 方程）的惩罚隐式-显式（IMEX）Runge-Kutta 格式，实现了无须昂贵的非线性碰撞算子隐式求逆的渐近保持性与精度。通过惩罚技术将碰撞项分解为平衡部分与非平衡部分，该方法在广泛变化的弛豫时间下实现了稳定性与高效性，经由三阶格式与数值实验验证。

ABSTRACT

We discuss Implicit-Explicit (IMEX) Runge Kutta methods which are particularly adapted to stiff kinetic equations of Boltzmann type. We consider both the case of easy invertible collision operators and the challenging case of Boltzmann collision operators. We give sufficient conditions in order that such methods are asymptotic preserving and asymptotically accurate. Their monotonicity properties are also studied. In the case of the Boltzmann operator, the methods are based on the introduction of a penalization technique for the collision integral. This reformulation of the collision operator permits to construct penalized IMEX schemes which work uniformly for a wide range of relaxation times avoiding the expensive implicit resolution of the collision operator. Finally we show some numerical results which confirm the theoretical analysis.

研究动机与目标

• 开发用于刚性非线性 kinetic 方程的高效数值格式，特别是在 Knudsen 数较小时的流体动力学极限情形。
• 通过避免直接求解非线性 Boltzmann 碰撞算子的隐式逆，克服刚性区域中隐式求解器带来的计算负担。
• 利用惩罚策略将 IMEX Runge-Kutta 方法扩展至完整的 Boltzmann 方程，同时保持渐近保持性与渐近精度特性。
• 分析并确保格式在同伦情形下的单调性与稳定性。
• 提供一个适用于具有弛豫算子与守恒量的大规模刚性 ODE 系统的通用框架。

提出的方法

• 将碰撞算子分解为 $ R(Y) = N(Y) + L(Y) $，其中 $ L(Y) $ 为线性化弛豫项，近似于平衡态处的雅可比矩阵。
• 引入惩罚技术以重构 Boltzmann 碰撞积分的生成项，从而实现对非线性部分的显式处理与对线性化弛豫部分的隐式处理。
• 构建 IMEX 格式，其中隐式阶段求解涉及惩罚算子的线性系统，避免了对完整非线性碰撞项求逆的需求。
• 通过底层 ODE 系统渐近分析导出的充分条件，确保渐近保持性。
• 设计兼具高阶精度与强稳定性的格式，包括广义阶段阶（GSA）与单调性条件。
• 使用 Butcher 表构造具体格式，包括一阶、二阶与三阶方法，通过参数调节实现对单调性区间的显式控制。

实验结果

研究问题

• RQ1 IMEX Runge-Kutta 格式能否在应用于完整非线性 Boltzmann 方程时仍保持渐近保持性与渐近精度？
• RQ2 在刚性区域中，如何在不牺牲稳定性或精度的前提下降低隐式碰撞求解器的计算成本？
• RQ3 何种惩罚策略可实现对 Boltzmann 碰撞项的高效隐式-显式积分，同时保持渐近行为？
• RQ4 在何种条件下，惩罚 IMEX 格式能保持单调性，特别是在同伦情形下？
• RQ5 能否构造出高阶（最高至三阶）的 IMEX 格式，使其同时具备渐近保持性与无条件稳定性？

主要发现

• 在 Butcher 表与惩罚参数满足充分条件时，惩罚 IMEX 格式对完整 Boltzmann 方程实现了渐近保持性与渐近精度。
• 惩罚技术使得避免对完整非线性碰撞算子进行昂贵的隐式求解成为可能，从而在广泛弛豫时间范围内实现高效计算。
• 一阶格式如 DP-A $(1,2,1)$ 与 DP-ARS $(1,2,1)$ 满足 AA（渐近精度）条件，且当 $ \rho \neq 1 $ 时可实现单调性。
• 二阶格式如 DP2-A $(2,4,2)$ 在 DIRK（对角隐式 Runge-Kutta）部分达到三阶精度，当 $ \rho \to \rho_{\text{min}} $ 时，若 $ \rho \neq 1 $，则保持单调性。
• 三阶格式如 DP1-A $(2,4,2)$ 采用四阶段构造，同时满足 GSA 与 AA 条件，可在刚性区域实现高阶精度。
• 数值结果验证了理论分析，显示在不同 $ \rho $ 值（包括 $ \rho = 10^{-6} $）下均表现出收敛性与稳定性，准确解析了流体动力学极限。