[論文レビュー] Isospectral Property of Hamiltonian Boundary Value Methods (HBVMs) and their blended implementation

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

Numerical methods for differential equations参考文献 20被引用数 4

ひとこと要約

本稿では、ハミルトニアン境界値法（HBVMs）の等スペクトル性を確立し、そのバーサテーブル行列が対応するガウス＝レジェンドル法と同一の非ゼロ固有値を持つことを証明することで、A安定性を保証するとともに、効率的なブレンド反復ソルバーの構築を可能にする。主な貢献は、最適収束を実現するパラメータ化されたブレンド反復法であり、線形安定性解析およびさまざまな s に対する数値結果によって検証されている。

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 Runge-Kutta methods, a matrix of the Butcher tableau with the same spectrum (apart the zero eigenvalues) as that of the corresponding Gauss-Legendre method, independently of the considered abscissae. Consequently, HBVMs are always perfectly A-stable methods. Moreover, this allows their efficient "blended" implementation, for solving the generated discrete problems.

研究の動機と目的

• HBVMsをルンゲ＝クッタ法として見た場合のスペクトル構造を分析すること。
• HBVMsがゼロ固有値を除き、ガウス＝レジェンドル法と等スペクトルであることを確立し、A安定性を保証すること。
• HBVMsに由来する離散非線形系に対する効率的なブレンド反復ソルバーを開発すること。
• 線形収束解析を用いて、ブレンド反復における自由パラメータ γ の最適選択を行うこと。
• ステップ数 k を増加させても、離散問題の条件数が有界のままであることを示すこと。

提案手法

• HBVMsはブロック境界値法（BVMs）として再定式化され、その後、ブレンド陰解法に再構成される。
• 行列 C 及びその逆行列を用いて、離散問題を2通りの同等な形で表現し、対角分割戦略を可能にする。
• 自由パラメータ γ を用いて、2つの表現を凸結合により組み合わせたブレンド反復を構築する。
• 反復の収束性は、スカラーテスト方程式 y′ = λy と q = hλ を用いて分析される。
• 反復行列 Z(q) が導出され、そのスペクトル半径 ρ(q) が解析され、A収束性および L 収束性が評価される。
• 最適な γ は γ = |μ_min| として選ばれ、最大増幅係数 ρ* = 1 − cos(φ_min) を最小化する。ここで μ_min は C の最小絶対値固有値である。

実験結果

リサーチクエスチョン

• RQ1 HBVMsのバーサテーブル行列は、使用される区間の選択に関わらず、対応するガウス＝レジェンドル法と同一の非ゼロ固有値を持つのか？
• RQ2 等スペクトル性を活用して、HBVMsのための効率的かつ収束性を有する反復ソルバーを設計できるか？
• RQ3 ブレンド反復における自由パラメータ γ の最適選択は何か、L 収束性を保証するには？
• RQ4 ステップ数 k を増加させても、離散問題の条件数はどのように変化するか？
• RQ5 s ≥ 1 および k ≥ s に対して、ブレンド反復は常に A 安定かつ L 収束するか？

主な発見

• HBVMsはガウス＝レジェンドル法と等スペクトルである：バーサテーブル行列は同一の非ゼロ固有値を持つため、A 安定性が保証される。
• 自由パラメータ γ を C の最小固有値の絶対値に設定した場合、HBVMsのブレンド反復は L 収束する。
• 最大増幅係数 ρ* は 1 − cos(φ_min) < 1 で有界であり、L 収束性が確認される。
• s = 2 から 10 の範囲で、最適な γ と対応する ρ* が計算され、表1に示されている。ρ* は 0.1340 から 0.6467 に増加している。
• 離散問題の条件数は s にのみ依存し、k → ∞ のときにも有界のままであるため、安定性が保証される。
• 1ステップあたり1回の行列分解で実装可能であり、計算効率が著しく向上する。