[論文レビュー] The Hamiltonian BVMs (HBVMs) Homepage

Luigi Brugnano, Felice Iavernaro+1 more|arXiv (Cornell University)|Feb 14, 2010

Numerical methods for differential equations参考文献 30被引用数 12

ひとこと要約

本稿は、離散線積分と直交多項式コロケーションを活用することで、任意の次数の多項式ハミルトニアンに対してハミルトニアン関数を正確に保存する、新しい数値積分法であるハミルトニアン境界値法（HBVMs）を導入する。主な貢献は、高次精度と優れた安定性を備えたエネルギー保存型スキームの開発であり、収束解析およびBiM/BiMDコードにおける計算実装によって検証されている。

ABSTRACT

Hamiltonian Boundary Value Methods (in short, HBVMs) is a new class of numerical methods for the efficient numerical solution of canonical Hamiltonian systems. In particular, their main feature is that of exactly preserving, for the numerical solution, the value of the Hamiltonian function, when the latter is a polynomial of arbitrarily high degree. Clearly, this fact implies a practical conservation of any analytical Hamiltonian function. In this notes, we collect the introductory material on HBVMs contained in the HBVMs Homepage, available at http://web.math.unifi.it/users/brugnano/HBVM/index.html

研究の動機と目的

• 任意の次数の多項式ハミルトニアンに対してハミルトニアン関数を正確に保存する数値解法の開発。
• ハミルトニアン系における標準的数値積分法が示すエネルギー保存の喪失を是正すること。
• 長期シミュレーションに適した高次精度および強い安定性を有するスキームのクラスの設計。
• 直交多項式およびコロケーション技術を用いたHBVMの理論的基盤の構築。
• ブレンド型式と収束解析を通じた効率的計算実装の実現。

提案手法

• HBVMsは、時間ステップ内での解を次数sの多項式σ(t)で近似し、初期条件と整合性を保つ。
• ハミルトニアンの時間ステップ内での変化を計算するために、離散線積分形式を採用し、正確な台形積分を実現するため、直交多項式（シフトされたルジャンドル基底）を用いる。
• 2段階の台形則を用いる：解の段階を表すabscissae {c_i} と、高次多項式被積分関数の正確な積分を保証するための別セット {ĉ_i}。
• 係数は直交多項式の射影および重み付き積分から導出されるコロケーション・ルンゲ＝クッタスキームとして定式化される。
• 収束性と安定性を向上させるために、パrameter化された反復行列を用いた予測子・修正子戦略を有するブレンド型実装が導入される。
• 反復解法の収束性は、固有値半径を用いて分析され、すべてのHBVM(k,s)スキーム（s ≥ 1 かつ k ≥ s）に対してL収束性が確立される。

実験結果

リサーチクエスチョン

• RQ1 任意の高次多項式ハミルトニアンに対して、正確にハミルトニアンを保存する数値解法を構築可能か？
• RQ2 ハミルトニアン系の数値解法において、高次精度とエネルギー保存を同時に達成できるか？
• RQ3 直交多項式は、最適な安定性を有するエネルギー保存型積分法を構築する上で果たす役割は何か？
• RQ4 HBVMから生じる非線形方程式系の反復的解法を、計算的に効率的かつ収束可能にできるか？
• RQ5 HBVMのブレンド型実装の収束性はどのようなものか？また、標準的陰的スキームと比較してどうなるか？

主な発見

• HBVMsは、任意の次数の多項式ハミルトニアンに対してハミルトニアンを正確に保存し、長期的なエネルギー保存を保証する。
• 解の段階を越える拡張されたabscissaeセットを用いることで、線積分の正確な台形積分が実現され、高次精度が達成される。
• HBVMのブレンド型実装はL収束性が証明されており、固有値半径ρ* < 1を満たすため、すべてのs ≥ 1およびk ≥ sに対して収束性が保証される。
• ガウス＝ルジャンドルに基づくHBVMsでは、最適パrameter γ = |μ_min|が最大増幅係数ρ*を最小化し、ρ* = 1 − cos(φ_min) < 1を達成する。
• 数値結果により、HBVMアプローチが、エネルギー保存および幾何的構造の維持において、標準的対称的・シンプレクティック法を上回ることが確認された。
• BiMおよびBiMDコードへの実装は、計算的有効性を示しており、理論的収束性および安定性結果の妥当性が裏付けられた。