[論文レビュー] Hamiltonian Boundary Value Methods (Energy Conserving Discrete Line Integral Methods)

Luigi Brugnano, Felice Iavernaro+1 more|arXiv (Cornell University)|Oct 19, 2009

Numerical methods for differential equations参考文献 17被引用数 171

ひとこと要約

本稿は、離散線積分とサイレントステージを用いて多項式ハミルトニアンを正確に保存するエネルギー保存型数値積分法であるハミルトニアン境界値法（HBVMs）を導入する。主な貢献は、サイレントステージの数が無限大に近づく際の理論的極限の導出であり、これにより最適な次数2sを達成する新しい方法—無限HBVMs—が得られ、エネルギー保存型コロケーション法を一般化する。

ABSTRACT

Recently, a new family of integrators (Hamiltonian Boundary ValueMethods) has been introduced, which is able to precisely conserve the energy function of polynomial Hamiltonian systems and to provide a practical conservation of the energy in the non-polynomial case. We settle the definition and the theory of such methods in a more general framework. Our aim is on the one hand to give account of their good behavior when applied to general Hamiltonian systems and, on the other hand, to find out what are the optimal formulae, in relation to the choice of the polynomial basis and of the distribution of the nodes. Such analysis is based upon the notion of extended collocation conditions and the definition of discrete line integral, and is carried out by looking at the limit of such family of methods as the number of the so called silent stages tends to infinity.

研究の動機と目的

• 離散線積分と拡張コロケーション条件に基づく、ハミルトニアン境界値法（HBVMs）の統一的枠組みを確立すること。
• サイレントステージ数が無限大に近づく際のHBVMsの極限を調査し、既存のエネルギー保存型手法との関連を明確にすること。
• 多項式基底とノード分布の影響を解析することで、最適なHBVM公式を同定すること。
• シフトされたルジャンドル多項式に基づく新しい極限法—無限HBVMs—を導入し、その特徴を明らかにすること。
• 極限法がノード分布に依存せず、同じ固有関数に収束することを示し、あらゆる定式化において一貫性のある振る舞いを示すこと。

提案手法

• 連続線積分の一般化として、離散線積分を用いてHBVMsを定式化し、多項式ハミルトニアンに対して正確なエネルギー保存を保証する。
• サイレントステージを基本ステージの線形結合として定義し、計算コストをsステージに増加させることなくエネルギー保存を実現する。
• 拡張コロケーション条件を用いてHBVM族を導出し、ハミルトニアン多項式の次数νに合わせてサイレントステージ数rを設定する。
• k = s + r → ∞ として極限をとることで極限法を導出し、ノード分布に依存せずに一意な極限に収束することを示す。
• ラグランジュ基底を用いたHBVMの極限とエネルギー保存型コロケーション法（EPCMs）の等価性を確立し、シフトされたルジャンドル基底を用いた新たな極限法を導出する。
• 得られた無限HBVMs（HBVM(∞,s)）が次数2sに達することを証明し、ガウス＝ルジャンドル法と同一の次数を達成するとともに、ノード選択に依存しないことを示す。

実験結果

リサーチクエスチョン

• RQ1 サイレントステージ数が無限大に近づく際のHBVMsの理論的極限は何か？
• RQ2 ラグランジュ基底とシフトされたルジャンドル基底の選択が、得られる極限法にどのように影響するか？
• RQ3 コロケーション過程で用いられるノード分布に依存せずに、極限法を独立に特徴づけられるか？
• RQ4 極限法は古典的ガウス＝ルジャンドル法と同一の次数を保つのか？
• RQ5 実用的観点から、既存のエネルギー保存型積分法と比較して、保存性および安定性の観点でどのように差がつくか？

主な発見

• サイレントステージ数k → ∞ におけるHBVMsの極限は、コロケーションノードの分布に依存しない一意な極限法に収束する。
• ラグランジュ基底を用いる場合、極限法はエネルギー保存型コロケーション法（EPCMs）と一致し、ハイメルの予想を裏付ける。
• シフトされたルジャンドル基底を用いることで、新たな極限法のクラス—無限HBVMs（HBVM(∞,s)）—が得られ、次数2sを達成する。
• 極限法は、ある特定の作用素の単位固有値に対応する固有関数によって特徴づけられ、ノード配置に依存しない。
• 数値実験では、HBVM(10,2)が長時間積分においてエネルギー保存性を10−12の精度で維持しており、Itoh-Abe法や4次公式に比べて優れた性能を示す。
• 多項式近似が機械精度に近い高次精度で行われるため、非多項式ハミルトニアンに対しても実用的なエネルギー保存性を示す。