[論文レビュー] Strong Stability Preserving Multistep Runge-Kutta Methods

Christopher Bresten, Sigal Gottlieb+4 more|arXiv (Cornell University)|Jul 30, 2013

Computational Fluid Dynamics and Aerodynamics参考文献 25被引用数 12

ひとこと要約

本稿では、単調性および安定性の性質を保つ高次の明示的強安定性保存（SSP）多段ルンゲ＝クッタ法を、双曲型PDEを解くために開発する。次数および単調性条件を導出し、2次以上の方法についてSSP係数の上限を証明し、数値最適化を用いて10次まで、8段まで、5段までにわたる最適化された方法を構築し、標準的手法よりも顕著に大きな有効SSP係数を達成する。

ABSTRACT

High-order spatial discretizations with strong stability properties (such as monotonicity) are desirable for the solution of hyperbolic PDEs. Methods may be compared in terms of the strong stability preserving (SSP) time-step. We prove an upper bound on the SSP coefficient of explicit multistep Runge--Kutta methods of order two and above. Order conditions and monotonicity conditions for such methods are worked out in terms of the method coefficients. Numerical optimization is used to find optimized explicit methods of up to five steps, eight stages, and tenth order. These methods are tested on the advection and Buckley-Leverett equations, and the results for the observed total variation diminishing and positivity preserving time-step are presented.

研究の動機と目的

• 全変動を減衰させる（TVD）および正性を保つといった強安定性特性を保持する高次明示的多段ルンゲ＝クッタ法の開発。
• 多段ルンゲ＝クッタ法の係数に基づいて、必要な次数および単調性条件を導出すること。
• 2次以上の明示的多段ルンゲ＝クッタ法のSSP係数の上限を確立すること。
• 数値アルゴリズムを用いて、次数2から10、段数5まで、8段までにわたる方法を最適化し、有効SSP係数を最大化すること。
• 移流方程式およびBuckley-Leverett方程式に対する数値実験を通じて、理論的SSP性質の妥当性を検証すること。

提案手法

• 方法の係数を用いて、明示的多段ルンゲ＝クッタ法の次数条件および単調性条件を導出する。
• 強安定性保存（SSP）性質を前進オイラー手順の凸結合として定式化し、時間刻みを小さくした場合の安定性を保証する。
• SSP係数Cを前進オイラー手順に対する許容可能な最大時間刻みとして定義し、有効SSP係数をC/s（sは段数）として定義する。
• 次数2から10、段数5まで、8段までにわたる方法について、有効SSP係数を最大化するように数値最適化を実施する。
• Buckley-Leverett方程式にはKorenリミッターを用いた保存的空間離散化を、移流方程式には半ラグランジュスキームを用いて方法を実装する。
• TVDおよび正性保存の観察時間刻みを理論的予測と比較し、SSP性質の妥当性を検証する。

実験結果

リサーチクエスチョン

• RQ1 明示的多段ルンゲ＝クッタ法の2次以上のSSP係数の理論的上限は何か？
• RQ2 多段ルンゲ＝クッタ法に対して、次数および単調性条件を体系的に導出する方法は何か？
• RQ3 段数および段数に制限のある高次明示的多段ルンゲ＝クッタ法で達成可能な最大の有効SSP係数は何か？
• RQ4 数値実験において、TVDおよび正性保存の観察時間刻みは理論的予測とどのように異なるか？
• RQ5 最適化された多段ルンゲ＝クッタ法は、双曲型PDEの安定性および効率性において、標準的手法を上回ることができるか？

主な発見

• 2次以上の明示的多段ルンゲ＝クッタ法について、凸結合フレームワーク下での最大安定性を制限するSSP係数の上限が証明された。
• 5段8段の10次最適化法は、有効SSP係数0.603を達成し、非SSP RK4,4法の0.250よりも顕著に高い値を示した。
• 移流方程式において、最適化された(7,3,3)法のTVD時間刻みは理論値の1.113倍に達し、理論的境界を超える強い安定性保存を示した。
• Buckley-Leverett問題において、観察された正性保存時間刻みは常にTVD時間刻みよりも大きく、安定性性質の選択が観察される安定性限界に影響することを示した。
• (20,3,10)法は、TVD時間刻み0.356Δxを観察し、有効SSP係数0.110を達成した。これは段数を増やすことで高次法の安定性が向上することを示している。
• TVDおよび正性保存の観察時間刻みは段数の増加に伴い増加し、より多くの段数が安定性特性を向上させることを理論的予想が裏付けた。