[論文レビュー] A class of generalized additive Runge-Kutta methods

Adrian Sandu, Michael Günther|arXiv (Cornell University)|Oct 21, 2013

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

ひとこと要約

本稿では、右辺の異なる成分に異なる段値を引数として使用できるように一般化された、一般化加法的ルンゲ＝クッタ（GARK）法を導入する。この枠組みにより、高次精度、安定性の向上、マルチスケールまたは複数の物理的プロセスを含む系に対する柔軟性が向上し、特に安定性特性が分離された、implicit-explicit（IMEX）およびimplicit-implicit形式が可能になる。

ABSTRACT

This work generalizes the additively partitioned Runge-Kutta methods by allowing for different stage values as arguments of different components of the right hand side. An order conditions theory is developed for the new family of generalized additive methods, and stability and monotonicity investigations are carried out. The paper discusses the construction and properties of implicit-explicit and implicit-implicit,methods in the new framework. The new family, named GARK, introduces additional flexibility when compared to traditional partitioned Runge-Kutta methods, and therefore offers additional opportunities for the development of flexible solvers for systems with multiple scales, or driven by multiple physical processes.

研究の動機と目的

• 右辺の異なる成分に異なる段値を入力として使用できるように一般化された加法的ルンゲ＝クッタ法のフレームワークを構築すること。
• この新しいクラスの方法について、次数条件理論、安定性、単調性解析を拡張すること。
• マルチスケールまたはマルチフィジックス問題に適した、安定性と柔軟性が向上したimplicit-explicitおよびimplicit-implicit GARKスキームの構築を可能にすること。
• 複数の剛性および非剛性成分を含む系に対する、高次、安定、かつ効率的なスケーラーの設計の基盤を提供すること。

提案手法

• GARK法は、従来のARKを一般化し、各成分 $ f^{\{q\}} $ に対して別個の段値 $ Y_i^{\{q\}} $ を導入することで、右辺評価における異なる引数を可能にする。
• 次数条件はNB級数理論を用いて導出され、古典的な次数条件フレームワークが一般化された設定に拡張される。
• 安定性は代数的安定性および絶対単調性を用いて分析され、GARK構造に適合した新しい定義が導入される。
• 結合係数 $ A^{\{q,m\}} $ を用いて、段値の依存関係を柔軟に制御できるため、implicit-explicit（IMEX）および完全にimplicit（implicit-implicit）形式の両方が可能である。
• 各成分間の相互作用を表すために、別個の係数行列を備えた一般化されたブーター表が導入される。
• 安定性特性が分離されたスキームの構築が可能であり、たとえば $ P^{\{i,j\}} = 0 $（$ i \neq j $）となるように設計できるため、安定性の制御が向上する。

実験結果

リサーチクエスチョン

• RQ1 右辺の異なる成分に異なる段値を引数として使用できるように一般化された加法的ルンゲ＝クッタフレームワークを構築することは可能か？
• RQ2 NB級数理論を用いて、この新しいクラスの方法の次数条件をどのように導出できるか？
• RQ3 一般化されたGARKフレームワークにおける代数的安定性および絶対単調性の条件は何か？
• RQ4 安定性と次数精度が向上したimplicit-explicitおよびimplicit-implicit GARKスキームを構築することは可能か？
• RQ5 マルチフィジックス系に対して、高次かつ安定性が分離されたGARK法を設計することは可能か？

主な発見

• 本稿では、NB級数を用いた完全な次数条件理論をGARK法に確立し、高次スキームの構築を可能にした。
• 代数的安定性の新しい定義を導入し、一般条件のもとで安定性を保証するGARKスキームへの応用を実現した。
• 段値 $ P^{\{i,j\}} = 0 $（$ i \neq j $）となるように設計された安定性分離スキームの構築が可能となり、安定性挙動の制御が向上した。
• $ P^{\{1,2\}} = P^{\{2,1\}} = 0 $ を満たす2次で代数的安定なGARKスキームが構築され、実用的応用の有効性が示された。
• 転置的古典的IMEXスキームの開発が可能となり、剛性精度の概念が一般化された設定に拡張された。
• 理論的分析により、ステップサイズ制限のもとで単調性が保証され、今後の研究のためのマルチレートおよびシンプレクティック拡張が可能である。