[論文レビュー] A unified formulation of splitting-based implicit time integration schemes

S. González‐Pinto, D. Hernández‐Abreu+4 more|arXiv (Cornell University)|Mar 1, 2021

Numerical methods for differential equations参考文献 71被引用数 7

ひとこと要約

本稿では、implicit-implicit一般構造加法的 Runge–Kutta (IMIM-GARK) 法を用いた、分割に基づく陰的時間積分スキームの統一的枠組みを提案する。ADI、作用素分割、分離ステップスキームといった古典的手法をこの枠組みに組み込み、次数条件を導出し、安定性および誤差特性が最適化された高次（3次および4次）の IMIM-GARK スキームを構築した。放物型 PDE に対する数値実験により、最適な収束性と効率性が確認された。

ABSTRACT

Splitting-based time integration approaches such as fractional steps, alternating direction implicit, operator splitting, and locally one-dimensional methods partition the system of interest into components and solve individual components implicitly in a cost-effective way. This work proposes a unified formulation of splitting time integration schemes in the framework of general-structure additive Runge-Kutta (GARK) methods. Specifically, we develop implicit-implicit (IMIM) GARK schemes, provide the order conditions and stability analysis for this class, and explain their application to partitioned systems of ordinary differential equations. We show that classical splitting methods belong to the IMIM GARK family, and therefore can be studied in this unified framework. New IMIM-GARK splitting methods are developed and tested using parabolic systems.

研究の動機と目的

• 分離ステップ、ADI、作用素分割、LOD などの多様な分割に基づく陰的時間積分法を、一つの数学的枠組みに統一すること。
• 分割系 ODE システムに適用可能な、implicit-implicit GARK (IMIM-GARK) スキームの体系的次数条件理論を構築すること。
• 精度と効率の向上を図るため、安定性定数および誤差定数が最適化された高次（3次および4次）IMIM-GARK スキームを構築すること。
• 正確な解が得られる放鳥型 PDE に対して、提案手法を検証し、最適な収束性および計算性能を示すこと。

提案手法

• 一般構造加法的 Runge–Kutta (GARK) フレームワーク内に、分割に基づく時間積分スキームを定式化し、成分ごとの陰的積分を可能にする。
• Butcher 表が下三角行列に並び替え可能である GARK 法として IMIM-GARK スキームを定義し、段階を逐次的かつ分離的に陰的に解けるようにする。
• 4次までの一貫した次数条件を導出し、高次スキームの体系的構築を可能にする。
• 古典的手法（ADI、作用素分割、分離ステップ Runge–Kutta）に本フレームワークを適用し、これらが IMIM-GARK の特殊な場合であることを示す。
• 導出した次数条件を用いて、安定性および誤差定数が最適化された3次および4次 IMIM-GARK スキームを構築する。
• 空間方向の半離散化に2次中央差分法を用い、時間誤差を正確な解との ℓ2 範囲で測定する。

実験結果

リサーチクエスチョン

• RQ1 ADI、作用素分割、分離ステップ法などの古典的手法が、一つの数学的枠組みに統合可能かどうか。
• RQ2 implicit-implicit GARK (IMIM-GARK) スキームの次数条件は何か。また、それらを用いて高次スキームを構築する方法は。
• RQ3 放鳥型 PDE に対して、分割に基づく IMIM-GARK スキームの安定性および誤差特性をどのように最適化できるか。
• RQ4 新規 IMIM-GARK スキームは、剛性が増す条件下でも、その名目上の次数の精度を実際の計算で達成できるか。
• RQ5 ベンチマーク放鳥型問題において、新規 IMIM-GARK スキームは、既存の GLM-ADI および IMEX-RK 手法と比べて、精度および効率の面で優れているか。

主な発見

• 提案された IMIM-GARK フレームワークは、AD I、作用素分割、分離ステップ法といった古典的手法を、一つの数学的定式化に統一的に統合することに成功した。
• 安定性定数および誤差定数が最適化された3次および4次 IMIM-GARK スキームが構築され、滑らかな問題において期待される古典的収束次数を達成した。
• 2次元および3次元の放鳥型テスト問題において、粗いメッシュでは名目上の収束次数が達成されたが、細かいメッシュでは剛性および時間依存境界条件の影響により、観測された収束次数が低下した。
• 並列 ADI-GARK 変種は、古典的収束次数を達成し、計算効率も高いことが示された。4次スキームは、一般の IMEX-RK 4 よりも誤差-CPU 時間のトレードオフにおいて優れていた。
• 効率性比較において、ADI-GARK 3 は IMEX-RK 4 に近い性能を示した。一方、ADI-GARK 4 は、最良の GLM-ADI 4 とほぼ同等の効率性を達成した。
• 数値結果から、IMIM-GARK フレームワークが、特に剛性および多次元 PDE に対して、高次、安定的、かつ効率的な時間積分スキームの設計を可能にすることが確認された。