[論文レビュー] Partially implicit Runge-Kutta methods for wave-like equations in spherical-type coordinates

I. Cordero-Carrión, P. Cerdá–Durán|arXiv (Cornell University)|Nov 26, 2012

Numerical methods for differential equations参考文献 13被引用数 4

ひとこと要約

本稿では、球対称座標系における波動型偏微分方程式の時間発展を安定化するため、部分的に暗黙のルンゲ・クッタ法を提案する。これは、剛性項や座標系に起因する要因によって引き起こされる数値的不安定性に対処するものである。この手法は、明示的ルンゲ・クッタ法よりも安定性を向上させつつ、長時間シミュレーションにおいても高次精度を維持する。

ABSTRACT

Partially implicit Runge-Kutta methods are presented in this work in order to numerically evolve in time a set of partial differential equations. These methods are designed to overcome numerical instabilities appearing during the evolution of a system of equations due to potential numerical unstable terms in the sources, such as stiff terms or the presence of factors as a result of a partic- ular chosen system of coordinates. In this article, partially implicit Runge-Kutta methods for several convergence orders have been derived and stability properties have been analyzed. These methods are shown to be appropriated to avoid the de- velopment of numerical instabilities in the evolution in time of wave-like equations in spherical-type coordinates, in contrast to the explicit Runge-Kutta methods.

研究の動機と目的

• 剛性の源となる項や座標系の選択によって生じる時間発展における数値的不安定性を解消すること。
• 球対称座標系に特化した高次部分的暗黙ルンゲ・クッタスキームの開発。
• 明示的スキームが失敗するような長時間シミュレーションにおいても安定性と精度を確保すること。
• 提案手法の安定性特性を分析し、その頑健性を検証すること。

提案手法

• 波動型PDEの時間積分に適した、複数の収束次数を有する部分的暗黙ルンゲ・クッタ法を導出する。
• 源項に含まれる剛性または不安定性を引き起こす項を、暗黙的成分に分離することで時間発展を安定化する。
• 座標特異性や因子が不安定性を引き起こす可能性がある球対称座標系における偏微分方程式系に、本手法を適用する。
• ルンゲ・クッタスキームの係数をバーサー図表形式で構造化する。
• ヴォン・ノイマン安定性解析を実施し、導出された手法の安定性特性を評価する。
• 標準的な明示的ルンゲ・クッタ法と比較して、部分的暗黙スキームの性能および安定性を検証する。

実験結果

リサーチクエスチョン

• RQ1 剛性項や座標に起因する要因によって引き起こされる、波動型方程式の時間積分における数値的不安定性は、どのように軽減可能か？
• RQ2 球対称座標系における波動型PDEに部分的暗黙ルンゲ・クッタ法を適用した際の安定性挙動はいかなるものか？
• RQ3 明示的ルンゲ・クッタ法と比較して、安定性が向上した状態で高次精度を維持できるか？
• RQ4 剛性系において安定性を確保するための、部分的暗黙ルンゲ・クッタスキームの最適係数は何か？
• RQ5 球座標系の文脈において、提案手法の安定性領域は明示的ルンゲ・クッタ法のそれと比べてどのように異なるか？

主な発見

• 提案された部分的暗黙ルンゲ・クッタ法は、剛性または座標に起因する項によって生じる数値的不安定性を効果的に抑制する。
• 安定性解析により、明示的ルンゲ・クッタ法が同様の条件下で発散するのに対し、本手法は長時間にわたり安定性を維持することが確認された。
• 導出されたスキームにおいても高次収束が保持されており、高精度な長時間シミュレーションが可能である。
• 本手法は、球対称座標系が引き起こす数学的課題（特異点や径方向要因など）に対処するように特に設計されている。
• 部分的暗黙の定式化により、明示的スキームと比較してより大きな安定な時間ステップが可能となり、計算効率が向上する。
• 提案されたスキームの安定性領域は、明示的ルンゲ・クッタ法と比較して顕著に広く、特に源項内の剛性成分に対して顕著である。