[論文レビュー] Energy-preserving continuous-stage partitioned Runge-Kutta methods

Wensheng Tang|arXiv (Cornell University)|Aug 7, 2018

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

ひとこと要約

本稿では、正規化されたシフト付きルジャンドル多項式と単純化仮定を用いることで、エネルギー保存を保証するための十分条件を導出し、ハミルトニアン系におけるエネルギー保存型連続段ルンゲ＝クッタ（csPRK）法を提案する。この手法は、両方のバーサー係数 $B_\tau$ と $\widehat{B}_\tau$ が 1 に等しいことを要求することで、正確なエネルギー保存を実現し、非線形系に対する数値実験により、機械精度に近いエネルギー保存性を確認した高次の積分法を構築する。

ABSTRACT

In this paper, we present continuous-stage partitioned Runge-Kutta (csPRK) methods for energy-preserving integration of Hamiltonian systems. A sufficient condition for the energy preservation of the csPRK methods is derived. It is shown that the presented condition contains the existing condition for energy-preserving continuous-stage Runge-Kutta methods as a special case. A noticeable and interesting result is that when we use the simplifying assumptions of order conditions and the normalized shifted Legendre polynomials for constructing high-order energy-preserving csPRK methods, both the Butcher "weight" coefficients $B_τ$ and $\widehat{B}_τ$ must be equal to $1$. As illustrative examples, new energy-preserving integrators are acquired by virtue of the presented condition, and for the sake of verifying our theoretical results, some numerical experiments are reported.

研究の動機と目的

• 長時間にわたるシミュレーションにおいても全エネルギーを正確に保存するエネルギー保存型数値積分法の開発。
• 連続段ルンゲ＝クッタ法の理論を、幾何的構造のより良い保存を実現するための分割型定式化に拡張すること。
• 連続段ルンゲ＝クッタ法に対する既存の結果を一般化する、csPRK法におけるエネルギー保存の十分条件の確立。
• 単純化仮定と正規化されたシフト付きルジャンドル多項式を用いて高次エネルギー保存型積分法を構築し、両方の $B_\tau = 1$ および $\widehat{B}_\tau = 1$ を保証すること。
• 非線形ハミルトニアン系における実験を通じて理論的結果の数値的妥当性を検証し、エネルギーおよび不変量の保存性をシンプレクティック法と比較すること。

提案手法

• 級数展開法を用いて、連続段分割ルンゲ＝クッタ（csPRK）法におけるエネルギー保存の十分条件を導出する。
• 次数条件の単純化仮定を適用し、正規化されたシフト付きルジャンドル多項式を用いて高次積分法を構築する。
• エネルギー保存を保証するため、両方のバーサー係数 $B_\tau$ と $\widehat{B}_\tau$ が 1 に等しいことを制約条件として課す。
• 数値的検証とテストを可能にするために、手法を非線形ハミルトニアン系に再定式化する。
• 数値積分に3点ガウス求積法を用い、2段4次精度のガウス＝ルジャンドルRK法（GLRK-4）と結果を比較する。
• 解のグローバル誤差および $H$、$I$、$L$ の不変量誤差の最大ノルム誤差を用いて性能を評価する。

実験結果

リサーチクエスチョン

• RQ1 連続段分割ルンゲ＝クッタ法におけるハミルトニアン系のエネルギー保存のための十分条件を導出可能か？
• RQ2 単純化仮定およびルジャンドル多項式に基づく構築において、バーサー係数 $B_\tau$ と $\widehat{B}_\tau$ はどのように振る舞うか？
• RQ3 提案されたcsPRK法は、GLRK-4などのシンプレクティック法と比較して、エネルギーおよび他の不変量をどの程度正確に保存するか？
• RQ4 導出されたエネルギー保存型csPRK法は、カオス的または剛性のある系において線形誤差成長と正確な長期的ダイナミクスを達成できるか？
• RQ5 提案された手法は、標準的なルンゲ＝クッタフレームワークを越えて、特にP系列積分法の文脈において適用可能か？

主な発見

• csPRK法におけるエネルギー保存の導出された十分条件は、既存の連続段ルンゲ＝クッタ法に対する条件を一般化する。
• 単純化仮定と正規化されたシフト付きルジャンドル多項式を用いる場合、エネルギー保存を達成するためには $B_\tau$ と $\widehat{B}_\tau$ がともに正確に 1 でなければならない。
• 数値実験の結果、提案されたすべてのエネルギー保存型手法が、ハミルトニアン $H$ を機械精度まで正確に保存している一方、シンプレクティック法のGLRK-4は二次的角運動量 $I$ を正確に保存している。
• エネルギー保存型手法は、グローバル解誤差および不変量誤差において線形誤差成長を示し、シンプレクティック法と同等の性能を示している。
• エネルギー保存型手法のすべての数値的軌道は、正確な楕円軌道をよく近似しており、高い精度と安定性を示している。
• エネルギー保存型手法は、シンプレクティックでないにもかかわらず、幾何的および動的性質の保存においてシンプレクティック法と同等の性能を発揮している。