[論文レビュー] Initializing and stabilizing variational multistep algorithms for modeling dynamical systems

C. Leland Ellison, J. W. Burby+3 more|arXiv (Cornell University)|Mar 4, 2014

Numerical methods for differential equations参考文献 24被引用数 3

ひとこと要約

本稿では、退化したラグランジュ関数を伴う力学系に対する変分多段積分法の安定化を目的として、後退誤差初期化および周期的再初期化を提案する。後退誤差解析によって得られる滑らかな修正系の解に沿って初期条件をサンプリングすることで、余剰モードの振幅を最小化し、標準的一段法初期化と比較して、保存量の滑らかさが著しく向上し、長期的な数値的整合性が向上する。

ABSTRACT

Backward error initialization and parasitic mode control are well-suited for use in algorithms that arise from a discrete variational principle on phase-space dynamics. Dynamical systems described by degenerate Lagrangians, such as those occurring in phase-space action principles, lead to variational multistep algorithms for the integration of first-order differential equations. As multistep algorithms, an initialization procedure must be chosen and the stability of parasitic modes assessed. The conventional selection of initial conditions using accurate one-step methods does not yield the best numerical performance for smoothness and stability. Instead, backward error initialization identifies a set of initial conditions that minimize the amplitude of undesirable parasitic modes. This issue is especially important in the context of structure-preserving multistep algorithms where numerical damping of the parasitic modes would violate the conservation properties. In the presence of growing parasitic modes, the algorithm may also be periodically re-initialized to prevent the undesired mode from reaching large amplitude. Numerical examples of variational multistep algorithms are presented in which the backward error initialized trajectories outperform those initialized using highly accurate approximations of the true solution.

研究の動機と目的

• 非物理的固有値に起因する余剰モードによって引き起こされる変分多段アルゴリズムの不安定性および滑らかさの欠如を解消する。
• 高精度であっても余剰モード振幅を最小化できない標準的一段法による初期化の限界を克服する。
• 保存性を維持しつつ数値アーチファクトを低減する実用的で構造保存型の初期化手法を開発する。
• 後退誤差解析を用いた周期的再初期化を導入し、保存則を破らずに増大する余剰モードを制御する。
• 非正準ハミルトニアン系の長期シミュレーションにおいて、従来手法に比べて後退誤差初期化が優れていることを実証する。

提案手法

• 多段変分積分法の滑らかな修正系を後退誤差解析により導出し、精度を向上させるために 𝒪(h²) 項まで保持する。
• 滑らかな修正系を用いて、余剰モードの振幅を最小化する初期条件を生成する。これは、望ましい物理的軌道と整合するように初期条件を設定することに相当する。
• 1段法を用いて滑らかな修正系を解くことで、周期的再初期化を実装し、余剰モード振幅をほぼゼロにリセットする。
• 修正系を元のODEの摂動として定式化する： ẋ̃ = f(x̃) + h²[−1/24 f''(x̃)f(x̃)f(x̃) + 1/12 f'(x̃)f'(x̃)f(x̃)]。この式は主な余剰モードの振る舞いを捉える。
• 記号的または自動微分を用いて、修正系に必要なfの高階微分を計算し、正確な初期化を可能にする。
• 後退誤差解析による初期化と高次の一様法（例：2次ルンゲ＝クッタ法）による初期化の軌道を比較し、保存量の滑らかさを評価する。

実験結果

リサーチクエスチョン

• RQ1 後退誤差初期化は、標準的一段法初期化に比べて、変分多段積分法における余剰モード振幅をより効果的に低減できるか？
• RQ2 後退誤差解析によって導かれる滑らかな修正系は、非正準ハミルトニアン系における多段アルゴリズムの長期的挙動をどのように改善するか？
• RQ3 滑らかな修正系を用いた周期的再初期化は、保存性を損なわず、増大する余剰モードをどの程度効果的に抑制できるか？
• RQ4 後退誤差初期化は、平面的点渦のような系におけるエネルギーなどの保存量の滑らかさに、定量的にどのような影響を与えるか？
• RQ5 後退誤差初期化の計算コストは標準初期化に比べてどの程度高く、その性能向上は妥当性を有するか？

主な発見

• 後退誤差初期化は、余剰モードの振動が検出不能に小さい軌道を生成する一方、標準的一段法初期化では高精度近似であっても明確な振動が残存する。
• 2段中点法は、後退誤差初期化を施した場合、2次ルンゲ＝クッタ法で初期化した場合と同等の初期精度を有しながらも、保存量の滑らかさが著しく向上する。
• 初期条件が後退誤差解析から導かれた滑らかな修正系の解に沿っているとき、余剰モード振幅が最小化される。
• 滑らかな修正系を用いた周期的再初期化は、余剰モードを低振幅にリセットする効果を発揮し、長時間にわたる積分時間における増大を防ぐ。
• 2段中点法の修正系は明示的に ẋ̃ = f(x̃) + h²[−1/24 f''(x̃)f(x̃)f(x̃) + 1/12 f'(x̃)f'(x̃)f(x̃)] として導出され、主な余剰モードダイナミクスを捉えている。
• 高階微分の計算と追加計算を要するが、標準的多段更新に比べて使用頻度が低いため、後退誤差初期化は計算的に実用可能である。