[論文レビュー] Symplectic integration with Jacobi polynomials

Wensheng Tang|arXiv (Cornell University)|Jun 8, 2018

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

ひとこと要約

本稿では、基底関数としてジャコビ多項式を用いることで、ハミルトニアン系に対するシンプレクティック連続段階ルンゲ＝クッタ法を提案する。連続段階フレームワーク内での直交多項式展開を活用することで、系の幾何的構造を保存する高次シンプレクティック積分法を構築し、数値実験において標準的な Radau IIB 法と比較して優れた長期的エネルギー保存性と解の精度を示した。

ABSTRACT

In this paper, we study symplectic integration of canonical Hamiltonian systems with Jacobi polynomials. The relevant theoretical results of continuous-stage Runge-Kutta methods are revisited firstly and then symplectic methods with Jacobi polynomials will be established. A few numerical experiments are well performed to verify the efficiency of our new methods.

研究の動機と目的

• 正準ハミルトニアン系の構造を保存する数値積分法の開発を目的とする。
• ルジャンドルおよびチェビシェフ多項式に限定されない一般のジャコビ多項式を用いた直交多項式に基づくシンプレクティック法の拡張を目的とする。
• ジャコビ多項式展開を用いた連続段階シンプレクティック法の体系的構築フレームワークを確立することを目的とする。
• ケプラー問題に対する数値実験を通じて、提案手法の効率性および幾何的精度を検証することを目的とする。

提案手法

• 本手法は、区間 [0,1] 上の関数 Aτ,σ, Bτ, Cτ でパrameter化された連続段階ルンゲ＝クッタ（csRK）スキームに基づく。
• 係数 Aτ,σ および Bτ の展開にジャコビ多項式を基底関数として用いることで、方法構築におけるスペクトル的精度を実現する。
• 得られる積分法の順序条件およびシンプレクティック性を保証するために、簡略化条件 (B̌(ξ), Č(η), Ď(ζ)) が適用される。
• 係数は、ジャコビ重み関数との内積を含む方程式系を解くことで導出され、シンプレクティック性および高次精度が保証される。
• 連続段階法を実用的で高次な陽的ルンゲ＝クッタスキームに離散化するために、5点ガウス＝ルジャンドル則が適用される。
• エネルギー保存性および解の誤差増大度を評価するために、解析解が既知のケプラー問題を用いて手法の妥当性が検証された。

実験結果

リサーチクエスチョン

• RQ1 ジャコビ多項式は、ハミルトニアン系に対する高次シンプレクティック積分法を構築するために効果的に用いることができるか？
• RQ2 連続段階フレームワークにおいて、ルジャンドルおよびチェビシェフ多項式からジャコビ多項式への直交多項式展開技術の一般化はどのように行われるか？
• RQ3 異なるジャコビ多項式パラメータの使用が、シンプレクティック積分法の安定性および精度に与える影響は何か？
• RQ4 得られたシンプレクティック法は、既存の Radau IIB 法と比較して、長期的エネルギー保存性および解の誤差においてどのように差を示すか？

主な発見

• パラメータ (α,β) = (0,1) のジャコビ多項式を用いた3段3次シンプレクティックルンゲ＝クッタ法が導出され、高い精度と安定性を示した。
• パラメータ (α,β) = (1,2) のジャコビ多項式を用いた5段5次シンプレクティック法が構築され、同程度の次数の Radau IIB 法と比較して優れたエネルギー保存性と解の精度を示した。
• ケプラー問題における数値実験では、提案手法が同等次数の Radau IIB 法と比較してわずかに優れたエネルギー保存性および解の誤差増大度を示した。
• 本手法フレームワークは一般性を有し、最小限の修正で任意の重み付き直交多項式系へ拡張可能である。
• 5段法のバーサ・テーブルが明示的に提示されており、コサイン関数の三角関数的表現から導かれた複雑だが構造的な係数が示されている。
• 理論的解析により、構築された手法がシンプレクティック条件 φ′h(z)TJφ′h(z) = J を満たすことが確認され、ハミルトニアン系の幾何的構造が保存されることを裏付けた。