[論文レビュー] Numerical Solution of ODEs and the Columbus' Egg: Three Simple Ideas for Three Difficult Problems

Luigi Brugnano, Felice Iavernaro+1 more|arXiv (Cornell University)|Aug 27, 2010

Numerical methods for differential equations被引用数 6

ひとこと要約

本稿では、常微分方程式（ODE）における長年の課題を解決する3つの革新的な数値的手法を提示する：コロケーション技法の単純な修正により線形 multistep 法におけるA-安定性を達成する方法、 stiff さを明確で計算可能な基準として定義する方法、およびハミルトニアン境界値法（HBVMs）を用いてハミルトニアン系におけるエネルギー保存を実現する方法。主な貢献は、HBVMs が、従来の保存則を満たさないシンプレクティック法がエネルギー漂移により失敗する高離心率を示すケプラー問題のような、 stiff かつ保存則を満たす系についても、安定かつ高精度な可変ステップサイズ統合を可能にすることを示したことである。

ABSTRACT

On computers, discrete problems are solved instead of continuous ones. One must be sure that the solutions of the former problems, obtained in real time (i.e., when the stepsize h is not infinitesimal) are good approximations of the solutions of the latter ones. However, since the discrete world is much richer than the continuous one (the latter being a limit case of the former), the classical definitions and techniques, devised to analyze the behaviors of continuous problems, are often insufficient to handle the discrete case, and new specific tools are needed. Often, the insistence in following a path already traced in the continuous setting, has caused waste of time and efforts, whereas new specific tools have solved the problems both more easily and elegantly. In this paper we survey three of the main difficulties encountered in the numerical solutions of ODEs, along with the novel solutions proposed.

研究の動機と目的

• A-安定な線形 multistep 法（LMM）を構築する長年の難問に対処する。これは、かつてダールグイスの障壁により制限されていた。
• ODEにおける stiff さの明確で計算可能な定義を提供し、既存の文献における曖昧さを解消し、効果的な数値的取り扱いを可能にする。
• 可変ステップサイズ統合においてもエネルギー漂移を回避し、エネルギーを保存するハミルトニアン問題用の数値的手法を開発する。
• 従来のシンプレクティック法（例：ガウス＝レジェンドル法）が可変ステップサイズに適さない理由を示し、HBVMs がエネルギー保存を維持し、誤差の線形増加を達成することを示す。

提案手法

• 境界値問題の定式化に埋め込むことで、コロケーション手法を修正し、A-安定なLMMを構築する。これにより、次数を損なわずに安定性を保証する。
• ヤコビ行列の最大固有値と最小固有値の比に基づく計算可能な条件を用いて、stiff さを定義する。これにより、一般のODEにおいて信頼性の高いstiffさ検出が可能になる。
• ハミルトニアン境界値法（HBVMs）を考案し、定式化における積分を高精度に計算すれば、ハミルトニアンを正確に保存する。
• 高次精度（例：12点則）のガウス＝レジェンドル則を用いてHBVMの積分を近似し、高次精度とエネルギー保存を両立させる。
• HBVMを用いた可変ステップサイズ戦略を実装し、エネルギー保存性を維持し、従来のシンプレクティック法で見られるエネルギー漂移を回避する。
• 可変ステップサイズ制御下で、HBVM(12,3)と6次精度のガウス＝レジェンドル法を比較し、離心率 e = 0.99 のケプラー問題をベンチマークとして用いる。

実験結果

リサーチクエスチョン

• RQ1 ダールグイスの障壁に反するものなく、A-安定な線形 multistep 法を構築可能か？ もしそうなら、これを達成する最も単純な修正は何か？
• RQ2 現代の数値ソフトウェアで使用可能な、明確で計算可能なODEにおけるstiffさの定義は存在するか？ これにより、曖昧またはヒューリスティックな基準に取って代われるか？
• RQ3 可変ステップサイズ統合においても安定かつ高精度を維持するハミルトニアン系のエネルギー保存型数値手法は設計可能か？ これは従来のシンプレクティック法とは対照的である。
• RQ4 なぜ従来のシンプレクティック法（例：ガウス＝レジェンドル法）は可変ステップサイズ制御下でエネルギー漂移と2次誤差増加を示すのか？ そして、これを回避するにはどうすればよいか？
• RQ5 HBVMsは、stiff あるいは高離心率を示すハミルトニアン問題において、ステップ数と精度の面でシンプレクティック法をどの程度上回れるか？

主な発見

• HBVM(12,3)は、離心率 e = 0.99 のケプラー問題において1000周期にわたり、ハミルトニアン誤差の増加が無視できるほど小さい一方、ガウス＝レジェンドル法では顕著なエネルギー漂移が観察された。
• HBVM(12,3)では解の誤差が時間とともに線形に増加するが、同じ可変ステップサイズ戦略下でガウス＝レジェンドル法では2次的に増加する。
• HBVM(12,3)は1000周期経過後に解の誤差が約6.85e-3に達するまでに1周期あたり約153ステップで十分であるのに対し、ガウス＝レジェンドル法は同等の精度を得るためには1周期あたり約2×10⁵の定常ステップが必要となる。
• HBVM(12,3)は、定式化における積分が高次則（例：12点ガウス則）で計算されれば、エネルギー保存性を維持することが確認され、エネルギー保存性の本質が裏付けられた。
• 従来のメッシュ選択戦略は、エネルギー漂移と誤差増加のため、シンプレクティック法では効果が薄いが、HBVMでは非常に効果的かつ安定に機能する。
• HBVM(3,3)は数学的に6次精度のガウス＝レジェンドル法と同等であるが、HBVMsはこれを高次、エネルギー保存型、可変ステップサイズ対応のスキームへ一般化している。