[論文レビュー] Chebyshev symplectic methods based on continuous-stage Runge-Kutta methods

Wensheng Tang|arXiv (Cornell University)|May 29, 2018

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

ひとこと要約

本稿では、第一種および第二種チェビシェフ直交多項式に基づく、重み付き連続段Runge-Kutta法の新しい枠組みを用いて構築された、新しいシンプレクティックRunge-Kutta法の族を提案する。提案手法は、高次の精度を達成するとともに、シンプレクティック構造を正確に保存し、ハミルトニアン系における数値実験で優れた長期安定性とエネルギー保存を示す。古典的ガウス＝ルジャンドル法よりも優れた定性的な挙動を示すが、初期段階での誤差成長はわずかに大きい。

ABSTRACT

We develop Chebyshev symplectic methods based on Chebyshev orthogonal polynomials of the first and second kind separately in this paper. Such type of symplectic methods can be conveniently constructed with the newly-built theory of weighted continuous-stage Runge-Kutta methods. A few numerical experiments are well performed to verify the efficiency of our new methods.

研究の動機と目的

• 長期間にわたり幾何的構造を保存する新しいハミルトニアン系用シンプレクティック積分法の開発を目的とする。
• 連続段Runge-Kutta法の理論を拡張し、第一種および第二種チェビシェフ多項式を用いてシンプレクティック法を構築することを目的とする。
• 古典的W変換技術はルジャンドル多項式に限定されるが、これを任意の重み付き直交多項式を用いて構築可能な代替手法に拡張することを目的とする。
• 摂動付きケプラー問題に対する数値実験を通じて、提案手法の効率性と幾何的精度を検証することを目的とする。

提案手法

• 重み付き連続段Runge-Kutta（csRK）法の新たに開発された理論を活用してシンプレクティック積分法を構築する。
• csRK法の係数を定義するために、第一種および第二種チェビシェフ多項式を基底関数として用いる。
• 高次精度およびシンプレクティック性を保証するために、単純化仮定（B̌(ξ), Č(η), Ď(ζ)）を適用する。
• 離散フローにおける正確なシンプレクティック性を保証するために、チェビシェフの求積点を用いた補間型求積則を適用する。
• 直交多項式の根に基づいて、3段4次および5段6次法の具体的なバーサ型テーブルを導出する。
• 積分条件を満たす関数Aτ,σ、Bτ、Cτを用いた連続段形式により、これらの手法を実装する。

実験結果

リサーチクエスチョン

• RQ1 連続段Runge-Kutta枠組み内において、チェビシェフ多項式を用いてシンプレクティック積分法を体系的に構築できるか？
• RQ2 チェビシェフに基づくシンプレクティック法と古典的ガウス＝ルジャンドル法との間で、エネルギー保存性および誤差成長の観点から性能に差異はあるか？
• RQ3 この枠組みはルジャンドル多項式に限定されず、例えばチェビシェフ多項式を含む他の多項式へも一般化可能か？
• RQ4 剛性または摂動付きハミルトニアン系に適用した場合、得られるシンプレクティック法の次数および安定性挙動はいかなるものか？

主な発見

• 4次および6次チェビシェフシンプレクティック法は、長時間積分区間において線形の誤差成長と有界なエネルギー誤差を示し、幾何的精度が確認された。
• 数値結果から、チェビシェフI型およびII型4次法は4次ガウス＝ルジャンドル法と同等の性能を示し、初期解誤差はわずかに大きいがエネルギー保存性は類似していた。
• 6次チェビシェフシンプレクティック法（I型およびII型）は、解誤差およびエネルギー保存性の両面で6次ガウス＝ルジャンドル法とほぼ同一の性能を示した。
• 離散フローにおけるシンプレクティック構造の保存が確認されたことから、これらの手法は正確にシンプレクティックである。これはスペクトルコロケーション法が近似的にしか実現しないのとは対照的である。
• 本手法の構築フレームワークは汎用的であり、チェビシェフに限らず、ヤコビ多項式やルジャンドル多項式など、他の重み付き直交多項式へも拡張可能である。