[論文レビュー] Geometry of Numerical Complex Time Integration

Thorsten Orendt, Jürgen Richter-Gebert+1 more|ArXiv.org|Mar 9, 2009

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

ひとこと要約

本稿は、線形常微分方程式における超収束を達成するために、複素数時間グリッドを用いたルンゲ・クッタ法の幾何的枠組みを提案する。再帰的に定義された基本ルンゲ・クッタ法をフラクタルに似た複素数の経路に沿って繰り返し適用することで、収束次数を1または2上昇させることができ、複素係数を有する組み合わせ法の新しい導出を可能にするとともに、実時間を超えた数値積分の構造的洞察を明らかにする。

ABSTRACT

We are studying Runge-Kutta methods along complex paths of integration from a geometric point of view. Thereby we derive special complex time grids, which applied to the problem of integrating a linear autonomous system of ordinary differential equations, can be used to achieve a classical superconvergence effect. The approach is also adapted for arbitrary ODEs. Furthermore we draw a connection from our geometric reasoning to the class of composition methods with complex coefficients. Thereby, our main goal is to introduce a new point of view on these methods.

研究の動機と目的

• 常微分方程式の解法において、数値的不安定性や特異点を克服するため、複素数時間統合経路を用いることの可能性を探ること。
• 複素数時間グリッドにおける幾何的選択が、数値積分における超収束にどのように寄与するかを調査すること。
• 代数的に順序条件を解くのとは異なる、複素係数を有する組み合わせ法の新しい幾何的導出を提供すること。
• 複素数常微分方程式統合における経路依存解とモノドロミーの影響を分析すること。
• 制限付き三体問題などのベンチマーク問題において、複素数時間統合の実用的・理論的利点を示すこと。

提案手法

• 基本法Ψから導出される再帰的に定義された複素数時間グリッドΓに沿って適用されるルンゲ・クッタ法を形式化する。
• 各段階rで、法の次数がg（明示的ルンゲ・クッタ法ではg=1、対称的ルンゲ・クッタ法ではg=2）だけ上昇する再帰的構成を用いる。
• 時間グリッドΓを複素平面内の経路として定義し、2項組み合わせ法における係数に対応するステップサイズσr,jを導入する。
• 解析的接続が単一分岐解であるため自明である線形自己同次系ẋ = Axに適用する。
• 統合経路が解析的接続を模倣する幾何的解釈を採用し、特異点を避けるために実数の経路を回避する迂回路を形成する。
• 等間隔グリッドを用いて、明示的オイラー法とΥ1組み合わせ法を用いて制限付き三体問題に対して妥当性を検証する。

実験結果

リサーチクエスチョン

• RQ1 線形常微分方程式において、ルンゲ・クッタ法が超収束を達成するような複素数時間グリッドを設計可能か？
• RQ2 複素数時間経路の幾何的構造は、複素係数を有する組み合わせ法の順序条件とどのように関係するか？
• RQ3 複素曲線に沿って統合する場合、経路依存性とモノドロミーが常微分方程式の解に与える影響は何か？
• RQ4 誤差低減と計算コストの観点から、複素数時間統合はどのような状況で標準的な実数値法を上回るか？
• RQ5 複素数時間グリッドの再帰的構造を用いて、高次組み合わせ法を体系的に導出可能か？

主な発見

• 線形自己同次常微分方程式では、適切に選択された複素数時間グリッドにより、収束次数が1（g=1）または2（g=2）上昇し、超収束が達成される。
• 法Υh_rの次数はp(r) = p + r·gで与えられ、gは基本ルンゲ・クッタ法の対称性および構造に依存する。
• 等間隔実数ステップを用いた明示的オイラー法は、同じステップ数でΥ1組み合わせ法に比べて誤差が55倍以上大きい結果となった。
• 複素数時間グリッドΓはフラクタルに似た構造を示し、そのステップサイズは先行研究における2項組み合わせ法の係数と完全に一致する。
• 幾何的アプローチにより、代数的に順序条件を解くのとは独立した、直感的で新しい複素係数組み合わせ法の導出が可能となった。
• 実数値法よりも計算コストが約6倍高いものの、剛性や特異性が顕著な状況において顕著な誤差低減効果を示した。