[論文レビュー] Embedded error estimation and adaptive step-size control for optimal explicit strong stability preserving Runge--Kutta methods

Sidafa Conde, Imre Fekete+1 more|arXiv (Cornell University)|Jun 22, 2018

Numerical methods for differential equations参考文献 22被引用数 5

ひとこと要約

本稿では、2次から4次までの最適な明示的強安定性保存（SSP）ルンゲ＝クッタ法の埋め込み対を導入し、信頼性の高い局所誤差推定と適応的ステップサイズ制御を可能にするとともに、非線形安定性を維持する。提案されたSSPERK対は、空間的に離散化された双曲型PDEの解法において優れた効率性と頑健性を示し、複数のテスト問題で性能が検証されたが、1つの過剰に慎重な対は、ステップサイズの過剰な縮小により性能が劣る。

ABSTRACT

We construct a family of embedded pairs for optimal strong stability preserving explicit Runge-Kutta methods of order $2 \leq p \leq 4$ to be used to obtain numerical solution of spatially discretized hyperbolic PDEs. In this construction, the goals include non-defective methods, large region of absolute stability, and optimal error measurement as defined in [5,19]. The new family of embedded pairs offer the ability for strong stability preserving (SSP) methods to adapt by varying the step-size based on the local error estimation while maintaining their inherent nonlinear stability properties. Through several numerical experiments, we assess the overall effectiveness in terms of precision versus work while also taking into consideration accuracy and stability.

研究の動機と目的

• 適応的ステップサイズ制御を可能にしつつ強安定性の性質を保持する最適な明示的SSPルンゲ＝クッタ法の埋め込み対の開発を目的とする。
• 大きな絶対安定性領域と非退化しない方法を確保し、双曲型保存則の解法における頑健性を向上させることを目的とする。
• 埋め込み式による正確な局所誤差推定を提供し、空間的に離散化されたPDEの効率的かつ信頼性の高い時間積分を可能とすることを目的とする。
• 新規のSSPERK対の精度、作業量、安定性の観点から、さまざまな双曲型問題における性能を評価することを目的とする。
• 数値実験に基づき、段数、安定性、計算効率のバランスを考慮した実用的な対を提案することを目的とする。

提案手法

• 同じブチャールテーブル係数を持つ高次数のSSPルンゲ＝クッタ法と低次数の埋め込み法を組み合わせることで、埋め込み対を構築する。
• 2つの解の差分を用いて、各ステップにおける局所切り捨て誤差を推定する。
• 誤差推定に基づき適応的ステップサイズ制御を適用し、誤差の均等分布を維持するとともに計算作業量を最小化する。
• 先行研究[5, 19]で定義された最適な誤差測定と非退化する方法であることを保証する。
• SSP法の安定性特性を活用し、双曲型PDEの時間積分中に非線形安定性を維持する。
• 異なる許容誤差を設定した剛性および非剛性の双曲型問題に対して、数値実験を用いて手法の有効性を検証する。

実験結果

リサーチクエスチョン

• RQ1 最適な明示的SSPルンゲ＝クッタ法に、強安定性を損なわずに適応的ステップサイズ制御を可能にする埋め込み対を構築できるか？
• RQ2 Bogacki–Shampine法やDormand–Prince法といった既存手法と比較して、新規のSSPERK埋め込み対は効率性および正確性においてどのように異なるか？
• RQ3 段数とSSP係数の影響は、埋め込み対の性能および頑健性にどのように現れるか？
• RQ4 理論的に最適であるにもかかわらず、実用的にはSSPERK(10,4)–b̃₂対が失敗するのはなぜか？
• RQ5 双曲型PDEを解く際に、安定性、正確性、計算効率のバランスが最も良い埋め込み対は何か？

主な発見

• 2次問題には、段数が少なく、頑健性に優れるSSPERK(4,2)–b̃₂、SSPERK(6,2)–b̃₂、およびSSPERK(3,2)–w対が推奨される。
• 3次問題においては、SSPERK(4,3)–b̃₂対がBogacki–Shampine 3(2)対と同程度の効率性と正確性を達成している。
• SSPERK(10,4)–b̃₂対は、局所誤差を一貫して過剰に推定し、ステップサイズが極端に小さくなり、計算コストが高くなるため、理論的最適性にもかかわらず実用的ではない。
• 残りのSSPERK(10,4)対（例：b̃₁, b̃₃, b̃₅, b̃₈）は類似した性能を示し、頑健性と効率性に優れるため、4次問題に適している。
• 安定性領域とSSP係数が大きい高段数のSSPERK法は、総作業量を削減するため、長時間シミュレーションでは段数の増加が全体の効率を向上させることを示している。
• 非常に厳しい許容誤差（≤1e-5）では、SSP性がもはや重要ではなくなる。この場合、非SSP法も同程度の性能を示す可能性がある。