[論文レビュー] Locally exact modifications of numerical integrators

Jan L. Cieśliński|arXiv (Cornell University)|Jan 3, 2011

Numerical methods for differential equations被引用数 2

ひとこと要約

本稿では、任意の点で常微分方程式の線形化を正確に保つ局所的正確な数値積分スキームの修正を導入し、より高い精度とA安定性を実現する。既存のスキーム、特に離散勾配法を修正することで、正準座標における任意の多次元ハミルトニアン系に対してエネルギー保存型でA安定な積分スキームを構築し、特に1次元系では精度向上が最大8桁に達するなど、数桁の精度向上を達成する。

ABSTRACT

We present a new class of exponential integrators for ordinary differential equations. They are locally exact, i.e., they preserve the linearization of the original system at every point. Their construction consists in modifying existing numerical schemes in order to make them locally exact. The resulting schemes preserve all fixed points and are A-stable. The most promising results concern the discrete gradient method (modified implicit midpoint rule) where we succeeded to preserve essential geometric properties and the final results have a relatively simple form. In the case of one-dimensional Hamiltonian systems numerical experiments show that our modifications can increase the accuracy by several orders of magnitude. The main result of this paper is the construction of energy-preserving locally exact discrete gradient schemes for arbitrary multidimensional Hamiltonian systems in canonical coordinates.

研究の動機と目的

• すべての点でODEの線形化を正確に保つ、局所的正確な数値積分スキームの開発。
• 特に離散勾配法を含む既存の数値スキームを、エネルギー保存などの幾何的性質を維持する形で修正すること。
• 局所的正確な手法を、正準座標における任意の多次元ハミルトニアン系に拡張すること。
• A安定性と幾何的構造を保ちながら、著しい数値精度の向上を達成すること。
• 従来のシンプレクティック法に制限されるが、エネルギー保存型積分スキームに可変時間ステップを組み込むこと。

提案手法

• 本手法は、各点まわりの線形化系に対する正確な離散化を用いて、既存の数値スキームを修正することで局所的正確な積分スキームを構築する。
• ヤコビ行列の行列指数を用いて線形化ODEの正確な解を適用し、δₙ = 2Ω⁻¹ₙ tan(hₙΩₙ/2) により修正された時間ステップを定義する。
• 本手法は離散勾配法を活用し、エネルギー保存性を持つ2段階スキームに再定式化する。
• 重要な式は δₙ⁻¹(xₙ₊₁ − xₙ) = ∇ₛᵀ(pₙ, pₙ₊₁) であり、任意のδₙに対して正確なエネルギー保存を保証する。
• 元の系の線形化ダイナミクスと線形化ODEの正確解を一致させることで、ヤコビ行列の行列指数を用いて導出する。
• 平衡点におけるTₚₚVₓₓを用いてΩₙ² = TₚₚVₓₓを定義することで、多次元ハミルトニアン系に一般化する。

実験結果

リサーチクエスチョン

• RQ1 エネルギー保存などの幾何的性質を保ちながら、数値積分スキームを局所的正確に修正することは可能か？
• RQ2 一般の多次元ハミルトニアン系に対して、局所的正確な修正を体系的に構築する方法は何か？
• RQ3 標準的な離散勾配法と比較して、局所的正確なスキームはどの程度精度を向上させられるか？
• RQ4 エネルギー保存型積分スキームに可変時間ステップを組み込むことは可能か？幾何的構造を損なわないか？
• RQ5 離散勾配スキームに対して、正確なエネルギー保存を維持する局所的正確な修正を構築することは可能か？

主な発見

• 提案された離散勾配法の局所的正確な修正は、正準座標における任意の多次元ハミルトニアン系に対しても、エネルギー積分を正確に保存する。
• 1次元ハミルトニアン系では、標準的な離散勾配スキームと比較して、精度が最大8桁向上する。
• 得られたスキームはA安定であり、元の系のすべての固定点を保存する。
• 本手法により可変時間ステップが可能となり、従来のシンプレクティック法では実現不可能な適応的統合戦略が可能になる。
• 線形化系の正確な離散化が基盤となり、線形系では既知の正確な調和振動子解に還元される。
• 本手法の構築は一般であり、平衡点におけるヤコビ行列を評価し、行列指数を用いることで、任意の自己準同型ODE系に適用可能である。