[论文解读] Preconditioned Implicit-Exponential (IMEXP) Time Integrators for Stiff Differential Equations

Vũ Thái Luân, Mayya Tokman+1 more|arXiv (Cornell University)|May 10, 2016

Numerical methods for differential equations参考文献 34被引用 3

一句话总结

本文提出了两种新型时间积分格式——隐式-指数（IMEXP）和混合指数（HImExp）——用于处理刚性微分方程，其中线性（L）和非线性（N）算子均为刚性。通过将L隐式处理并与N的指数积分相结合，并支持任意预条件器的使用，该方法在高度刚性问题中实现了更优的稳定性与效率，相较于预条件化的2-sBDF方法，计算时间节省最高达87%。

ABSTRACT

We propose two new classes of time integrators for stiff DEs: the implicit-explicit exponential (IMEXP) and the hybrid exponential methods. In contrast to the existing exponential schemes, the new methods offer significant computational advantages when used with preconditioners. Any preconditioner can be used with any of these new schemes. This leads to a broader applicability of exponential methods. The proof of stability and convergence of these integrators and numerical demonstration of their efficiency are presented.

研究动机与目标

• 为解决现有IMEX与指数积分格式在刚性ODE中线性和非线性算子均为刚性时计算效率低下的问题。
• 开发新型时间积分格式，能够有效利用现有预条件器处理两个刚性算子，克服标准指数方法的局限性。
• 为具有刚性非线性的问题提供一种稳定、收敛且计算高效的IMEX格式替代方案。
• 通过理论分析与数值实验表明，新方法在高刚性及使用预条件器时，显著优于标准2-sBDF与IMEX格式。

提出的方法

• 提出两类新型时间积分格式：隐式-指数（IMEXP）与混合指数（HImExp），专为分别处理刚性线性和非线性项而设计。
• 对非线性项N(u(t))采用指数积分，同时对线性项L(u(t))采用隐式处理，从而提升对刚性非线性的稳定性。
• 提出一种新公式，允许将任意现有预条件器应用于L或N，增强灵活性与计算效率。
• 推导了刚性精确格式，并建立了新方法的收敛性与稳定性理论，确保其在长时间积分中的鲁棒性。
• 采用Krylov子空间方法计算矩阵指数，降低指数积分格式的计算成本。
• 通过二维Allen-Cahn与Schnakenberg反应-扩散方程的数值验证，采用不同刚性参数与预条件处理方式。

实验结果

研究问题

• RQ1 当标准IMEX方法因显式处理N而失效时，指数积分格式能否有效扩展以处理刚性非线性算子？
• RQ2 如何系统性地将预条件器整合到指数积分格式中，以提升对线性和非线性刚性项的效率？
• RQ3 当L与N均为刚性时，新IMEXP与HImExp格式的稳定性与收敛性行为如何？
• RQ4 新方法能否在高刚性与预条件处理下，显著节省2-sBDF及其他IMEX格式的计算时间？
• RQ5 新格式在非线性项N(u(t))刚性程度增加时，性能如何变化？

主要发现

• HImExp2J格式在所有测试问题中均优于2-sBDF与其它IMEXP方法，在最刚性情况下CPU时间最高可减少87%。
• 对于γ = 5×10⁴的二维Schnakenberg系统，预条件化的HImExp2J方法在10⁻³精度下将CPU时间降至2-sBDF的20%，而非预条件化时仅为70%。
• HImExp2J的计算优势随刚性增加而增强：在10⁻³容差下，预条件化处理的最刚性Schnakenberg案例中，节省高达87%。
• 对于Allen-Cahn方程，在ε = 0.02且10⁻²精度下，HImExp2J实现55%的CPU时间减少，即使未使用预条件器。
• 在ε = 0.005且10⁻³精度下，非预条件化的HImExp2J比2-sBDF更快，表明其改进的稳定性允许采用更大的时间步长。
• 数值结果验证了理论预测：HImExp2J在不同刚性水平与预条件处理下均保持优异性能，证实其鲁棒性与可扩展性。