[论文解读] Stabilized integrating factor Runge-Kutta method and unconditional preservation of maximum bound principle

Jingwei Li, Xiao Li+2 more|arXiv (Cornell University)|Apr 20, 2021

Advanced Numerical Methods in Computational Mathematics参考文献 44被引用 5

一句话总结

本文提出了一种稳定化积分因子Runge-Kutta（sIFRK）方法，可无条件保持半线性抛物方程（如Allen-Cahn方程）的最大界原理（MBP）。通过将线性稳定化与积分因子Runge-Kutta框架相结合，作者推导出MBP保持的充分条件，并构建了达到三阶精度的sIFRK格式，满足这些条件，从而在离散设置下实现无条件MBP保持，这与许多现有高阶格式（包括超过一阶的SSP-sIFRK变体）形成对比。

ABSTRACT

Maximum bound principle (MBP) is an important property for a large class of semilinear parabolic equations, in the sense that the time-dependent solution of the equation with appropriate initial and boundary conditions and nonlinear operator preserves for all time a uniform pointwise bound in absolute value. It has been a challenging problem on how to design unconditionally MBP-preserving high-order accurate time-stepping schemes for these equations. In this paper, we combine the integrating factor Runge-Kutta (IFRK) method with the linear stabilization technique to develop a stabilized IFRK (sIFRK) method, and successfully derive sufficient conditions for the proposed method to preserve MBP unconditionally in the discrete setting. We then elaborate some sIFRK schemes with up to the third-order accuracy, which are proven to be unconditionally MBP-preserving by verifying these conditions. In addition, it is shown that many classic strong stability-preserving sIFRK schemes do not satisfy these conditions except the first-order one. Extensive numerical experiments are also carried out to demonstrate the performance of the proposed method.

研究动机与目标

• 开发高阶时间积分格式，无条件保持半线性抛物方程的最大界原理（MBP）。
• 克服现有高阶指数积分格式（如ETD和IFRK）在二阶以上无法无条件保持MBP的局限性。
• 建立稳定化积分因子Runge-Kutta（sIFRK）格式中无条件MBP保持的充分条件。
• 识别并验证阶数最高达三阶的sIFRK格式，使其满足这些条件并实现无条件MBP保持。
• 通过数值实验表明，许多经典的强稳定保形（SSP）sIFRK格式在非一阶时无法保持MBP。

提出的方法

• 该方法将积分因子Runge-Kutta（IFRK）方法与线性稳定化技术相结合，为半线性抛物方程构建sIFRK格式。
• 以Butcher表系数和稳定化参数的形式推导出无条件MBP保持的充分条件。
• 通过解析验证，确认sIFRK格式在阶数最高达三阶时满足这些条件，确保无论时间步长大小，离散MBP均被保持。
• 该方法应用于Allen-Cahn方程及其他相场模型，空间离散化采用有限差分或有限元方法。
• 对非线性项施加稳定化以控制刚性，实现无条件稳定且MBP保持的时间积分。
• 通过在不同时间步长及二维和三维区域中的数值验证，评估精度、稳定性与MBP保持性能。

Fig. 1: Evolutions of energies (left) and supremum norms (right) calculated by using sIFRK(1,1) (top), sIFRK( $s$ ,2) (middle) and Heun-sIFRK(3,3) (bottom) schemes with $\tau=0.01$ .

实验结果

研究问题

• RQ1 高阶（超过二阶）时间积分格式能否无条件保持半线性抛物方程的最大界原理？
• RQ2 sIFRK格式需满足何种条件，才能在离散设置中实现无条件MBP保持？
• RQ3 为何许多现有强稳定保形（SSP）sIFRK格式即使在二阶时也无法无条件保持MBP？
• RQ4 所提出的sIFRK框架是否与标准空间离散化方法（如有限差分和有限元）兼容？
• RQ5 sIFRK方法能否在大时间步长下保持精度与稳定性，同时保持MBP？

主要发现

• 通过合理选择稳定化参数，sIFRK方法可无条件保持半线性抛物方程的最大界原理，且不受时间步长影响。
• 成功构建了三阶sIFRK格式，并通过满足推导出的充分条件，证明其可无条件保持MBP。
• 数值实验表明，SSP-sIFRK(2,2)违反了MBP保持条件，在t=5.5后无法保持边界，而标准sIFRK(2,2)则完美保持边界。
• sIFRK(2,2)在不同时间步长（τ=0.1, 0.05, 0.01）下性能一致，解分布与能量演化差异可忽略。
• 在Allen-Cahn方程的三维模拟中，sIFRK(2,2)保持了MBP，并对所有测试时间步长（包括τ=0.1）确保了能量单调衰减。
• 研究发现，除一阶外，大多数现有SSP-sIFRK格式不满足推导出的条件，因此无法保证无条件MBP保持。

Fig. 2: Initial configuration of the bubble. Left: the surface-project view; right: the cross-section view at $y=0$ .