[論文レビュー] 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

ひとこと要約

本稿では、アレン＝チャン方程式を含む半線形放物型方程式に対して、最大有界原理（MBP）を無条件に保存する安定化積分因子ルンゲ＝クッタ（sIFRK）法を導入する。線形安定化と積分因子ルンゲ＝クッタフレームワークを組み合わせることで、MBP保存の十分条件を導出し、それらを満たす3次までのsIFRKスキームを構築した。これにより、時間刻みサイズに依存せずに離散的条件下での無条件MBP保存が保証される。これは、2次以降の精度を持つ多くの既存の高次スキーム、特に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 などの既存の高次指数積分スキームが2次を超えてMBPを無条件に保存できないという制限を克服すること。
• 安定化積分因子ルンゲ＝クッタ（sIFRK）スキームにおける無条件MBP保存のための十分条件を確立すること。
• これらの条件を満たし、MBPを無条件に保存する3次までのsIFRKスキームを特定・検証すること。
• 数値実験を通じて、多くの古典的強安定性保存（SSP）sIFRKスキームが1次でない限りMBPを無条件に保存しないことを示すこと。

提案手法

• 積分因子ルンゲ＝クッタ（IFRK）法と線形安定化技術を組み合わせることで、半線形放物型方程式のためのsIFRKスキームを構築する。
• MBPの無条件保存のための十分条件を、バーサールテーブル係数と安定化パラメータの観点から導出する。
• 3次までのsIFRKスキームについて、これらの条件を解析的に検証し、時間刻みサイズに依存せずに離散的MBP保存が保証されることを確認する。
• アレン＝チャン方程式やその他の相場モデルにこの手法を適用し、有限差分法や有限要素法による空間離散化を実施する。
• 非線形項に安定化を施すことで剛性を制御し、無条件に安定かつMBPを保存する時間積分を可能にする。
• さまざまな時間刻みサイズと2次元・3次元領域を用いた数値的検証により、精度、安定性、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 2次を超える高次（2次を超える）時間積分スキームは、半線形放物型方程式に対して最大有界原理を無条件に保存できるか？
• RQ2 sIFRKスキームが離散的条件下で無条件にMBPを保存するには、どのような条件を満たす必要があるか？
• RQ3 多くの既存の強安定性保存（SSP）sIFRKスキームが、2次でも無条件にMBPを保存できないのはなぜか？
• RQ4 提案されたsIFRKフレームワークは、有限差分法や有限要素法といった標準的な空間離散化と互換性があるか？
• RQ5 sIFRK法は、大きな時間刻みサイズを用いても、MBPを保存しながら精度と安定性を維持できるか？

主な発見

• 適切な安定化パラメータを備えたsIFRK法は、時間刻みサイズに依存せず、半線形放物型方程式に対して最大有界原理（MBP）を無条件に保存する。
• 3次までのsIFRKスキームが、導出された十分条件を満たすように明示的に構築され、MBPを無条件に保存することが証明された。
• 数値実験の結果、SSP-sIFRK(2,2)はMBP保存条件を満たさず、t=5.5以降に境界を逸脱するが、標準的sIFRK(2,2)は境界を完璧に維持する。
• sIFRK(2,2)は、時間刻みサイズτ=0.1、0.05、0.01のすべてにおいて、解のプロファイルとエネルギーの変化にほとんど差がなく、一貫性を示した。
• アレン＝チャン方程式の3次元シミュレーションにおいて、sIFRK(2,2)はMBPを保存し、全時間刻みサイズ（τ=0.1を含む）でエネルギーの単調減少を保証した。
• 本研究により、1次スキームを除き、大多数の既存のSSP-sIFRKスキームが導出された条件を満たさないことが判明し、したがって無条件にMBP保存を保証できないことが明らかになった。

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