[論文レビュー] An asynchronous leapfrog method II

Ulrich Mutze|arXiv (Cornell University)|Nov 26, 2013

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

ひとこと要約

本稿では、時間刻みの固定入力に依存する古典的ルンゲ・クッタ法の欠点を解消し、速度を状態変数として用いることで、時間刻み依存性を排除した非同期フロント法を提案する。この手法により、効果的な自動時間刻み制御が可能となり、2次精度と安定性を維持する。平均化を施したバージョンは、剛性問題に対して優れた耐性を示し、適応的時間刻み制御下でルンゲ・クッタ法と比較して最大4倍の精度を達成する。

ABSTRACT

A second order explicit one-step numerical method for the initial value problem of the general ordinary differential equation is proposed. It is obtained by natural modifications of the well-known leapfrog method, which is a second order, two-step, explicit method. According to the latter method, the input data for an integration step are two system states, which refer to different times. The usage of two states instead of a single one can be seen as the reason for the robustness of the method. Since the time step size thus is part of the step input data, it is complicated to change this size during the computation of a discrete trajectory. This is a serious drawback when one needs to implement automatic time step control. The proposed modification transforms one of the two input states into a velocity and thus gets rid of the time step dependency in the step input data. For these new step input data, the leapfrog method gives a unique prescription how to evolve them stepwise. The stability properties of this modified method are the same as for the original one: the set of absolute stability is the interval [-i,+i] on the imaginary axis. This implies exponential growth of trajectories in situations where the exact trajectory has an asymptote. By considering new evolution steps that are composed of two consecutive old evolution steps we can average over the velocities of the sub-steps and get an integrator with a much larger set of absolute stability, which is immune to the asymptote problem. The method is exemplified with the equation of motion of a one-dimensional non-linear oscillator describing the radial motion in the Kepler problem.

研究の動機と目的

• 固定時間刻み入力に依存する古典的フロント法が、適応的時間刻み制御を扱えないという制限を解消すること。
• 変動する剛性を示す現実世界の力学系に適した一般用途用の、安定的で2次精度の明示的積分法を開発すること。
• 特に漸近的挙動や剛性を示す問題に対して、フロント法の絶対安定領域を改善すること。
• 安定性を損なわず精度を維持する新たな平均化技術の有効性を示すこと。
• 積分法の精度を軌道の指紋として可視化するための「数値的相互作用像」を導入・検証すること。

提案手法

• 時間刻みサイズに依存しない入力として、2番目の状態の代わりに速度を用いることで、フロント法を再定式化し、ステップ入力を時間刻みから分離する。
• 部分ステップにおける2段階の平均化処理を適用し、絶対安定領域が著しく拡大された新しい積分法を構築する。
• 初期および最終の速度値から導出される曲率に基づく誤差推定子（κ）を用い、自動時間刻み調整を支援する。
• 推定された局所誤差に基づいて、しきい値（kinkCrit）とステップサイズ修正係数（f1, f2）を用いたフィードバック制御方式を採用し、動的時間刻み調整を実現する。
• 正確な軌道と数値的軌道を比較する「数値的相互作用像」を導入し、位相空間における積分法固有の指紋を明らかにする。
• 適応的制御下で公平な比較が行えるよう、標準的な2次ルンゲ・クッタ法およびストルマー・ヴェルレ法を同じ入力形式に適合させる。

実験結果

リサーチクエスチョン

• RQ1 フロント法を、時間刻み依存性を失い、安定性や精度を損なわず自動時間刻み制御を可能にするように変更可能か？
• RQ2 2番目の状態の代わりに速度を入力に用いる場合、フロント法の安定領域はどのように変化するか？
• RQ3 連続するフロントステップを平均化することで、剛性問題に対する積分法の安定性と精度にどのような影響を与えるか？
• RQ4 適応的時間刻み制御下で、提案手法は既存の2次ルンゲ・クッタ法およびストルマー・ヴェルレ法と比較して、精度と効率の面で優れているか？
• RQ5 数値的相互作用像は、位相空間における異なる積分法の特徴的な診断的特徴を明らかにできるか？

主な発見

• 自動時間刻み制御下で、非同期フロント法は、関数評価回数が約1/4に抑えられるにもかかわらず、2次ルンゲ・クッタ法と比較して最大4倍の精度を達成する。
• 非同期フロント法の平均化バージョンは、わずかな減衰と広い絶対安定領域のおかげで、剛性初期条件に対して優れた耐性を示す。
• 非剛性問題においても、非同期フロント法は、ルンゲ・クッタ法やストルマー・ヴェルレ法を含む全テスト手法の中で最も高い精度を示す。
• 数値的相互作用像は、明確に区別できる軌道の指紋を明らかにし、非同期フロント法は位相空間内での滑らかでより正確な進化を示す。
• 曲率推定子κに基づく自動時間刻み制御は、剛性の広い範囲にわたり精度を維持することができ、適応的実行では時間刻みが125倍の変動を示す一方、固定刻みではたった3倍に留まる。
• 適応的制御下での本手法の性能は、理想の定常時間刻み性能の約3倍低効率であるため、ステップサイズ選択論理の最適化の余地がある。