[論文レビュー] Runge-Kutta Characterization of the Generalized Summation-by-Parts Approach in Time

Pieter D. Boom, David W. Zingg|arXiv (Cornell University)|Oct 1, 2014

Numerical methods for differential equations参考文献 23被引用数 6

ひとこと要約

本稿では、一般化された和分型（GSBP）時間積分法が、陰的ルンゲ＝クッタ法と同等であることを確立し、SBPおよびGSBP法に既存のルンゲ＝クッタの次数条件および安定性理論を適用可能にする。L-安定かつ代数的安定性を満たす新しいGSBP法、特に4段7次完全陰的法を導出し、密度ノルム条件下での非線形安定性を証明する。

ABSTRACT

This article extends the theory of dual-consistent summation-by-parts (SBP) and generalized SBP (GSBP) time-marching methods by showing that they are implicit Runge-Kutta schemes. Through this connection, the accuracy theory for the pointwise solution, as well as the solution projected to the end of each time step, is extended for nonlinear problems. Furthermore, it is shown that these minimum guaranteed order results can be superseded by leveraging the full nonlinear order conditions of Runge-Kutta methods. The connection to Runge-Kutta methods is also exploited to derive conditions under which SBP and GSBP time-marching methods associated with dense norms are nonlinearly stable. A few known and novel Runge-Kutta methods with associated GSBP operators are presented. The novel methods, all of which are L-stable and algebraically-stable, include a four-stage seventh-order fully-implicit method, a three-stage third-order diagonally-implicit method, and a fourth-order four-stage diagonally-implicit method.

研究の動機と目的

• SBPおよびGSBP時間積分法の精度および安定性理論を、線形問題を超えて非線形問題へと拡張すること。
• GSBP時間積分と陰的ルンゲ＝クッタ法の間の正式な関係を確立すること。
• ルンゲ＝クッタの次数条件および安定性理論を活用して、保証された安定性を有する新しい高次GSBP法を導出すること。
• 密度ノルムを有するGSBP法が非線形安定性を達成するための条件を同定すること。
• L-安定および代数的安定性を有する新規のGSBP埋め込みルンゲ＝クッタ法を提示すること。

提案手法

• 理論的分析により、GSBP時間積分スキームがそのブッチャー表の構造を通じて、陰的ルンゲ＝クッタ法と同等であることが示される。
• 古典的ルンゲ＝クッタ次数条件を適用し、非線形問題におけるGSBPスキームの最小保証次数を導出する。
• 非線形安定性は、特にL-安定性および代数的安定性条件を活用したルンゲ＝クッタ理論との関係を介して分析される。
• 既知のルンゲ＝クッタ法（例えば、対角陰的および完全陰的スキーム）をGSBPフレームワークに埋め込むことで、新しいGSBP演算子が構築される。
• 密度ノルムを有するGSBPスキームの安定性条件の導出は、基礎となるルンゲ＝クッタ法の代数的構造に依存する。
• 新規の方法が構築され、検証され、特に4段7次完全陰的GSBP法および3種の対角陰的法（3段3次および4段4次）が含まれる。

実験結果

リサーチクエスチョン

• RQ1 ルンゲ＝クッタの次数条件を用いて、SBPおよびGSBP時間積分法の精度理論を非線形問題へと拡張可能か？
• RQ2 密度ノルムを有するGSBPスキームが非線形安定性を達成するための条件は何か？
• RQ3 既知のどのルンゲ＝クッタ法をGSBPフレームワークに埋め込むことで、新しいGSBP時間積分法が得られるか？
• RQ4 ルンゲ＝クッタ理論から、系統的に新しい高次、L-安定、代数的安定GSBP法を導出可能か？
• RQ5 ルンゲ＝クッタ特徴付けは、GSBP時間積分の次数および安定性特性をどのように向上させるか？

主な発見

• GSBP時間積分スキームは、正式に陰的ルンゲ＝クッタ法と同等であり、既存のルンゲ＝クッタの次数および安定性理論の適用が可能になる。
• 非線形問題におけるGSBPスキームの最小保証精度は、ルンゲ＝クッタ法の完全非線形次数条件を用いて拡張される。
• 密度ノルムを有するGSBP法の非線形安定性は、基礎となるルンゲ＝クッタ法がL-安定性または代数的安定性を満たす場合に達成される。
• L-安定かつ代数的安定性を満たす4段7次完全陰的GSBP法が新たに導出された。
• 3段3次および4段4次の2種類の新規対角陰的GSBP法が提示された。両者ともL-安定かつ代数的安定性を満たす。
• 理論的枠組みにより、既知のルンゲ＝クッタスキームを活用して高次かつ安定なGSBP時間積分スキームを構築可能となる。