[论文解读] On the uniform accuracy of implicit-explicit backward differentiation formulas (IMEX-BDF) for stiff hyperbolic relaxation systems and kinetic equations

Jingwei Hu, Ruiwen Shu|arXiv (Cornell University)|Dec 2, 2019

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

一句话总结

本文证明了隐式-显式后向差分公式（IMEX-BDF）在所有松弛参数 ε 的参数范围内，对刚性双曲弛豫系统和动力学方程具有统一精度。通过一种新颖的能量估计方法，结合自定义的乘子技术，作者建立了与 ε 无关的最优先验误差界，确保在标准格式常出现阶数降低的中间 regime 中仍保持高阶精度。

ABSTRACT

Many hyperbolic and kinetic equations contain a non-stiff convection/transport part and a stiff relaxation/collision part (characterized by the relaxation or mean free time $\varepsilon$). To solve this type of problems, implicit-explicit (IMEX) multistep methods have been widely used and their performance is understood well in the non-stiff regime ($\varepsilon=O(1)$) and limiting regime ($\varepsilon ightarrow 0$). However, in the intermediate regime (say, $\varepsilon=O(Δt)$), uniform accuracy has been reported numerically without a complete theoretical justification (except some asymptotic or stability analysis). In this work, we prove the uniform accuracy -- an optimal {\it a priori} error bound -- of a class of IMEX multistep methods, IMEX backward differentiation formulas (IMEX-BDF), for linear hyperbolic systems with stiff relaxation. The proof is based on the energy estimate with a new multiplier technique. For nonlinear hyperbolic and kinetic equations, we numerically verify the same property using a series of examples.

研究动机与目标

• 建立 IMEX-BDF 格式在中间 regime（ε = O(Δt)）中统一精度的严格理论基础，此前虽有数值证据但缺乏完整理论证明。
• 弥补对 IMEX 格式在刚性双曲方程和动力学方程中阶数降低现象的理论理解空白，特别是在 ε 既非极小也非 O(1) 的情况下。
• 开发并应用一种专为双曲弛豫系统的对流-弛豫结构量身定制的新型乘子技术，以证明稳定性和收敛性。
• 通过数值验证将分析扩展至非线性动力学方程，确认其在非线性情形下与线性情形中观察到的统一精度行为一致。

提出的方法

• 作者采用一种新设计的乘子技术，结合能量估计，专门适配刚性双曲弛豫系统的结构，克服了现有乘子技术的局限性。
• 该方法应用于具有刚性弛豫的线性双曲系统，如模型系统 (1.1)，并进一步推广至变系数情形。
• 对于动力学方程，通过 BGK 模型进行数值验证，其中 IMEX-BDF 格式通过基于矩的非线性系统求解实现。
• 采用谱空间离散化方法，实现空间上的高阶精度，从而支持严格的误差分析。
• 证明建立了形式为 ‖u(T,·)−U(T,·)‖L² + ‖v(T,·)−V(T,·)‖L² ≤ CΔt^q 的最优先验误差界，且常数 C 与 ε 无关。
• 使用符号计算推导 q=4 的 IMEX-BDF 格式的高阶系数，确保稳定性与相容性条件以高达 10^-31 的精度满足。

实验结果

研究问题

• RQ1 IMEX-BDF 格式在所有松弛参数 ε 取值下，对刚性双曲弛豫系统是否可严格证明其统一精度？
• RQ2 在中间 regime 中，由于标准乘子失效，为建立 IMEX-BDF 格式的稳定性，需要何种新颖的数学技术？
• RQ3 IMEX-BDF 格式的统一精度特性是否可推广至非线性动力学方程（如 BGK 模型）？
• RQ4 当 ε 与时间步长 Δt 同阶时，IMEX-BDF 格式的精度阶数如何表现？其理论基础是什么？

主要发现

• 本文证明了对线性双曲弛豫系统应用 IMEX-BDF 格式时，存在形式为 ‖u(T,·)−U(T,·)‖L² + ‖v(T,·)−V(T,·)‖L² ≤ CΔt^q 的最优先验误差界，且常数 C 与 ε 无关。
• 统一精度在 IMEX-BDF 格式中可保持至四阶，无论 ε 处于非刚性、刚性还是中间 regime，均能实现高阶收敛。
• 提出了一种新型乘子技术，成功捕捉了对流-弛豫结构，从而实现了统一稳定性和收敛性的证明。
• 数值结果表明，即使在非线性隐式格式下，该统一精度特性在非线性动力学方程（如 BGK 方程）中依然成立。
• 分析表明，一阶格式无法在 ε 较小时捕捉解中的 O(ε) 项，凸显了高阶格式对解析中间尺度动力学的必要性。
• 对于 q=4，作者将系数计算精度达到 10^-31，相关矩阵 G 和 A 的最小特征值分别约为 5.03×10^-6 和 0.0765，确认了数值稳定性。