[論文レビュー] Isospectral Property of Hamiltonian Boundary Value Methods (HBVMs) and their connections with Runge-Kutta collocation methods

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

Numerical methods for differential equations参考文献 16被引用数 7

ひとこと要約

この論文は、ハミルトニアン境界値法（HBVMs）が、Runge-Kuttaスキームとして定式化された場合、ガウス＝ルジャンドル法と等スペクトル的であることを確立している。これは、両者のバーサテーブル行列が同一の非ゼロ固有値を持つことを意味する。この性質により、完全なA安定性が保証され、HBVMsと古典的Runge-Kuttaコロケーション法の間の深い構造的関係が明らかになる。特にガウス求積ノードを用いる場合に顕著である。

ABSTRACT

One main issue, when numerically integrating autonomous Hamiltonian systems, is the long-term conservation of some of its invariants, among which the Hamiltonian function itself. Recently, a new class of methods, named Hamiltonian Boundary Value Methods (HBVMs) has been introduced and analysed, which are able to exactly preserve polynomial Hamiltonians of arbitrarily high degree. We here study a further property of such methods, namely that of having, when cast as a Runge-Kutta method, a matrix of the Butcher tableau with the same spectrum (apart from the zero eigenvalues) as that of the corresponding Gauss-Legendre method, independently of the considered abscissae. Consequently, HBVMs are always perfectly A-stable methods. This, in turn, allows to elucidate the existing connections with classical Runge-Kutta collocation methods.

研究の動機と目的

• HBVMsをRunge-Kutta法として定式化した際のスペクトル的性質を調査すること。
• 特にエネルギー保存性の文脈において、HBVMsと古典的Runge-Kuttaコロケーション法との間の関係を明確にすること。
• HBVMsが等スペクトル性を通じてガウス＝ルジャンドル法のA安定性を継承することを示すこと。
• スペクトル的および求積に基づく解析を用いて、HBVMsの次数と安定性の理論的基盤を提供すること。
• ガウスノードに基づくHBVMsが、高次数のガウスコロケーション法と密接に関連していることを確立すること。

提案手法

• HBVMsは、特にルジャンドル多項式を用いた直交多項式に基づくコロケーションフレームワークを用いて、明示的Runge-Kutta法に再定式化される。
• この手法は、ガウス求積ノードから導かれる射影行列 $\mathcal{P}_s$ とコロケーション行列 $\mathcal{A}$ を用いた変換によりバーサテーブルを構築する。
• 等スペクトル性は、HBVMテーブルの行列 $A$ が $\mathcal{P}^{-1}A\mathcal{P} = \widetilde{X}_s \oplus O$ を満たすことで証明され、これはガウス法の構造と一致する。
• 解析は、多項式の次数が $2s-1$ まで正確に積分できるという単純化仮定 $B(2s)$ を活用する。
• HBVM$(k,s)$ がガウスノードを用いる場合、フィルタリング変換を通じて次数 $2k$ のガウスコロケーション法と関連していることが示され、関係性が明確化される。
• 次数 $2s$ の別証明が、$T(s,s)$ 条件と行列分解を用いて提示され、既知の次数理論と整合していることが確認される。

実験結果

リサーチクエスチョン

• RQ1 HBVMsがRunge-Kutta法として定式化された場合、ガウス＝ルジャンドル法と同一の非ゼロ固有値を持つのか？
• RQ2 ガウス求積ノードを用いる場合、HBVMsと古典的Runge-Kuttaコロケーション法との間の構造的関係は何か？
• RQ3 HBVMsの等スペクトル性が、A安定性および長期的エネルギー保存性にどのように影響するか？
• RQ4 求積および行列スペクトル解析を用いて、HBVMsの次数を厳密に導出できるか？
• RQ5 フィルタリング行列 $\mathcal{P}_s \mathcal{P}_s^T \Omega$ は、HBVMsとガウス法を結ぶ役割を果たすか？

主な発見

• HBVMsはガウス＝ルジャンドル法と等スペクトル的である：バーサテーブル行列の非ゼロ固有値は、使用されるAbszissenにかかわらず同一である。
• この等スペクトル性により、HBVMsは完全なA安定性を保証され、剛性のあるハミルトニアン問題において顕著な利点を有する。
• HBVMsがガウスノードに基づく場合、次数 $2k$ のガウスコロケーション法と構造的に関連しており、これはシンプレクティック性が保証されていることで知られている。
• HBVM$(k,s)$ の次数は $2s$ であり、$T(s,s)$ 条件と行列分解を用いて証明され、古典的次数理論と整合していることが確認された。
• 安定性および次数は、$B(2s)$ 条件を満たす限り、Abszissenの分布とは無関係に成立する。
• 変換 $A = \mathcal{A} (\mathcal{P}_s \mathcal{P}_s^T \Omega)$ は、HBVMsがガウス法のフィルタリング版として作用することを示し、重要なスペクトル的および安定性的性質を保持している。