[論文レビュー] A general symplectic integrator for canonical Hamiltonian systems

Yonghui Bo, Wenjun Cai+1 more|arXiv (Cornell University)|Dec 2, 2019

Numerical methods for differential equations参考文献 22被引用数 1

ひとこと要約

本稿では、実数パラメータを用いて、正則ハミルトニアン系の一般化されたシンプレクティック積分法を提案する。この手法は、シンプレクティックEuler法と陰的中点法を統合する。生成関数と対称的合成を活用し、高次のシンプレクティックスキームを構築し、パラメータの調整によりエネルギー保存型スキームを実現する。厳密な可解性の証明と数値的検証を実施している。

ABSTRACT

The focus of this paper is to recommend a novel symplectic scheme for canonical Hamiltonian systems. The new scheme contains a real parameter which makes the symplectic Euler methods and implicit midpoint rule as its special cases. The validity of the symplecticity of this scheme is well explained from the perspectives of the generating function and partitioned Runge-Kutta methods. The generating function of the new symplectic scheme with new coordinates is studied, and these coordinates include the three typical coordinates. Employing the generating function method and symmetric composition methods, two new classes of symplectic schemes of any high order are devised, respectively. Furthermore, based on these symplectic schemes, two energy-preserving schemes and the feasibility of constructing higher order energy-preserving schemes are presented by the parameter serving for clever tuning. The solvability of all the schemes mentioned is proved, and the numerical performances of these schemes are demonstrated with numerical experiments.

研究の動機と目的

• 正則ハミルトニアン系のための統一的なシンプレクティック積分法を構築し、既存の手法を一般化すること。
• 生成関数と分割Runge-Kutta理論を用いて、新規スキームのシンプレクティック性を確立すること。
• 生成関数と対称的合成法を用いて、高次のシンプレクティックスキームを構築すること。
• スキームのパラメータを巧みに調整することで、エネルギー保存型スキームを設計すること。
• 提案されたすべてのスキームの可解性を証明し、数値的性能を検証すること。

提案手法

• 実数パラメータを有する一般化されたシンプレクティックスキームを導入し、特殊な場合としてシンプレクティックEuler法と陰的中点法に還元されることを示す。
• 生成関数と分割Runge-Kutta形式を用いてシンプレクティック性を分析し、幾何的構造の保存を保証する。
• 生成関数から新しい座標を導出し、3つの代表的な座標タイプを含むことで、スキームの構築を容易にする。
• 対称的合成法を用いて、基本スキームから高次のシンプレクティックスキームを生成する。
• 生成関数と対称的合成技術を用いて、2つの新しい高次シンプレクティックスキームのクラスを構築する。
• スキームのパラメータを調整することでエネルギー保存を達成し、エネルギー保存型積分法の設計を可能にする。

実験結果

リサーチクエスチョン

• RQ1 1つのシンプレクティック積分法フレームワークが、パラメータを用いてシンプレクティックEuler法と陰的中点法を統合できるか？
• RQ2 生成関数と分割Runge-Kutta理論を用いて、一般化されたスキームのシンプレクティック性をどのように証明できるか？
• RQ3 対称的合成は、基本積分法から高次のシンプレクティックスキームを構築する際に果たす役割は何か？
• RQ4 導入されたパラメータを調整することで、ハミルトニアン系におけるエネルギー保存を達成できるか？
• RQ5 提案されたシンプレクティックスキームおよびエネルギー保存型スキームの可解性と数値的性能はいかがなものか？

主な発見

• 提案された一般化されたシンプレクティックスキームは、シンプレクティックEuler法と陰的中点法を特別な場合として含み、統一的なフレームワークを確立する。
• 生成関数と分割Runge-Kutta理論を用いた厳密な証明により、スキームのシンプレクティック性が裏付けられる。
• 生成関数と対称的合成法を用いて、2つの新しい高次シンプレクティックスキームのクラスが成功裏に設計された。
• スキームのパラメータを巧みに調整することでエネルギー保存型スキームが構築され、高次エネルギー保存型積分法の実現可能性が示された。
• すべての提案スキームの可解性が数学的に証明され、数値的適用可能性が保証された。
• 数値実験により、幾何的および物理的構造を良好に保存するスキームの優れた数値的性能が確認された。