[論文レビュー] Explicit non-canonical symplectic algorithms for charged particle dynamics

He Yang, Yajuan Sun+3 more|arXiv (Cornell University)|Sep 25, 2015

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

ひとこと要約

本稿では、ハミルトニアンを4つの解ける部分系に分解することで、荷電粒子運動における明示的K-シンプレクティック数値解法を提示する。各部分系は非正則シンプレクティック構造を保存する。得られた解法は2次精度を示し、長期的安定性とエネルギー保存性を有しており、ガロア運動やガイドセンター拡散の長期シミュレーションにおいて、標準的なルンゲ・クッタ法を上回る性能を発揮する。

ABSTRACT

We study the non-canonical symplectic structure, or K-symplectic structure inherited by the charged particle dynamics. Based on the splitting technique, we construct non-canonical symplectic methods which is explicit and stable for the long-term simulation. The key point of splitting is to decompose the Hamiltonian as four parts, so that the resulting four subsystems have the same structure and can be solved exactly. This guarantees the K-symplectic preservation of the numerical methods constructed by composing the exact solutions of the subsystems. The error convergency of numerical solutions is analyzed by means of the Darboux transformation. The numerical experiment display the long-term stability and efficiency for these methods.

研究の動機と目的

• 電磁場内での荷電粒子運動に対する長期的安定な数値解法の開発、特にトカマクのようなマルチタイムスケール問題に特化して。
• 従来の正則シンプレクティック法に起因する制限、すなわち非正則系に対しては陰的かつ計算コストが高くなるという点を解決する。
• ローレンツ力系が持つ内在的な非正則シンプレクティック（K-シンプレクティック）構造を活用し、明示的かつ構造保存型積分子を構築する。
• ハミルトニアンの分割と正確な解の合成により、K-シンプレクティック構造を保存することで、長期的な精度と安定性を確保する。
• ガロア運動およびガイドセンター拡散の長期的シミュレーションにおいて、標準的なルンゲ・クッタ法と比較して本手法の優位性を検証する。

提案手法

• スケューマトリック行列 $ K({\bf z}) $ を用いてローレンツ力系を非正則ハミルトニアン形式に再定式化し、$ ({\bf x}, {\bf v}) $ 座標系におけるK-シンプレクティック構造を明らかにする。
• ハミルトニアンを3つの速度依存項 $ H_{v_i} $ と1つのポテンシャル項 $ H_{\rho} $ に分割し、それぞれが同一のK-シンプレクティック構造を持つ部分系に対応する。
• 磁場および電場の構造に基づく解析的解を用いて各部分系を正確に解き、明示的時間ステッピングを可能にする。
• 4つの部分系の正確なフローを合成することで全数値解法を構築し、合成により離散的K-シンプレクティック保存が保証される。
• 誤差解析には、K-シンプレクティック構造を正則構造に写像するダーブォー変換を用い、2次精度解法の $ h^2 $ に基づく誤差境界を導出する。
• 2次元対称電磁場内での数値実験により、長期的挙動、エネルギー誤差、解誤差について4次ルンゲ・クッタ法と比較する。

実験結果

リサーチクエスチョン

• RQ1 荷電粒子運動の非正則シンプレクティック（K-シンプレクティック）構造を保存する明示的数値解法を構築可能か？
• RQ2 ハミルトニアン分割法をどのようにしてローレンツ力系を同一のK-シンプレクティック構造を持つ解ける部分系に分解できるか？
• RQ3 得られたK-シンプレクティック解法の長期的安定性および精度は、標準的なルンゲ・クッタ法と比較してどの程度優れているか？
• RQ4 K-シンプレクティック構造の保存が、長時間積分におけるエネルギー誤差および解誤差をどの程度制限するか？
• RQ5 固定誤差許容範囲下で、本手法の計算コストは古典的手法と比較してどの程度か？

主な発見

• 構築された2次精度K-シンプレクティック解法は、$ 2 \times 10^4 $ ステップの間、エネルギー誤差を $ O(h^2) $ の範囲に維持するが、4次ルンゲ・クッタ法は顕著なエネルギー減衰を示す。
• 190,000ステップ経過後、K-シンプレクティック解法は安定したガイドセンター拡散を伴う正しいらせん軌道を保存するが、ルンゲ・クッタ法はガロア運動の数値的減衰により失敗する。
• K-シンプレクティック解法の解誤差は時間経過に伴い緩やかに増加し、特に大きなステップサイズにおいて、長時間積分後もルンゲ・クッタ法の誤差を下回る。
• 解誤差が1%未満である場合、2次精度K-シンプレクティック解法は1次精度バージョンよりも計算時間を短縮するため、より高い効率性を示す。
• ダーブォー変換により理論的誤差解析が可能となり、K-シンプレクティック解法を正則シンプレクティック解法に結びつけることで、$ O(h^2) $ 収束率の正当化がなされた。
• 数値結果は、K-シンプレクティック解法が構造保存性のおかげで、特にマルチタイムスケール系において長期的シミュレーションで顕著に優れていることを確認した。