[論文レビュー] Systematic Improvement of Splitting Methods for the Hamilton Equations

Asif Mushtaq, Anne Kværnø+1 more|arXiv (Cornell University)|Apr 18, 2012

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

ひとこと要約

本稿では、標準的な Störmer-Verlet スキームを上回る精度を達成するための体系的な手法を提示する。組み合わせ技術を用いずに、$\tau^8$ の精度を達成する。時間積分の正確なフローのテイラー展開から導かれる時間依存係数を用いて、分割スキームにおける運動エネルギー項とポテンシャルエネルギー項を修正することで、シンプレクティック性を保ちつつ、長期的なエネルギー保存性と位相空間の精度を著しく向上させる。

ABSTRACT

We show how the standard (St{ö}rmer-Verlet) splitting method for differential equations of Hamiltonian mechanics (with accuracy of order $τ^2$ for a timestep of length $τ$) can be improved in a systematic manner without using the composition method. We give the explicit expressions which increase the accuracy to order $τ^8$, and demonstrate that the method work on a simple anharmonic oscillator.

研究の動機と目的

• ハミルトニアン系のシンプレクティック分割スキームの精度を、組み合わせを用いない体系的で非組み合わせ的な手法により向上させること。
• 標準的な分割法の限界に対処すること。標準的な手法は第二位数精度に留まり、長期的なエネルギー漂移を生じる。
• シンプレクティック構造を維持しながら、$\tau^8$ の精度を達成するための明示的補正項を導出すること。
• 非線形系、例えば非調和振動子に対する本手法の有効性を示すこと。エネルギー保存性と周期の保存性に優れていることを示す。

提案手法

• 本手法は、正確なフローのテイラー展開から導かれる時間依存係数を用いて、分割スキームにおける運動エネルギー項とポテンシャルエネルギー項を修正する。
• ハミルトニアン成分のネストされた交換子から導かれる補正項 $G_2, G_4, \dots, G_8$ を、フローの生成子に加える。
• 離散的時間発展演算子が調和振動子の正確な解と $\tau^8$ まで一致するように、係数 $m(\tau)$ と $k(\tau)$ を選ぶ。
• 複数ステップの組み合わせを避ける代わりに、直接的に分割パラメータを調整することで、組み合わせを用いない。
• 非線形系では、各時間ステップで補正パラメータを求めるために非線形代数方程式を解く必要がある。
• シンプレクティック性は構成上保たれるため、エネルギーと位相空間軌道の長期的安定性が保証される。

実験結果

リサーチクエスチョン

• RQ1 組み合わせ技術を用いずに、ハミルトニアン系の標準的分割法の精度を $\tau^2$ を超えて体系的に向上させることは可能か？
• RQ2 シンプレクティック性を保ちながら $\tau^8$ の精度を達成するための明示的補正項は何か？
• RQ3 運動エネルギー項およびポテンシャルエネルギー項に対する時間依存係数の修正が、エネルギー保存性および位相空間構造の長期的保存性にどのように影響するか？
• RQ4 非線形系、例えば非調和振動子において、本手法は周期の精度とエネルギー保存性をどの程度向上させるか？
• RQ5 高次補正を適用する際の計算コストは、標準的手法と比較してどの程度か？

主な発見

• 本手法は、正確な調和振動子のフロー展開を用いて分割パラメータを体系的に補正することで、ハミルトニアン系の解に対して $\tau^8$ の精度を達成する。
• エネルギー保存性が著しく向上する：長時間積分において、$(H - \frac{1}{2})/\tau^8$ の値が小さく保たれる。これは高い精度を示している。
• ハミルトニアン $H = \frac{1}{2}p^2 + \frac{1}{4}q^4$ の非調和振動子に対して、$\tau^6$ 補正スキームは260,000周期以上にわたりエネルギー保存性を維持する。
• シンプレクティック性が保たれるため、大きな時間刻みでも長期的安定性と正確な位相空間構造が保たれる。
• $\tau^6$ 補正スキームは、正確な解と比較して振動周期のずれが最小限に抑えられ、低次のスキームを上回る性能を示す。
• ステップごとの計算コストは非線形方程式の解法により増加するが、各ステップの高い精度が、高精度シミュレーションにおいてコストを上回る利点をもたらす。