[論文レビュー] Asymptotically stable particle-in-cell methods for the Vlasov-Poisson system with a strong external magnetic field

Francis Filbet, L. Miguel Rodrigues|arXiv (Cornell University)|Nov 23, 2015

Gas Dynamics and Kinetic Theory参考文献 43被引用数 63

ひとこと要約

本稿では、強い外部磁場下におけるVlasov-Poisson系の漸近的保存性を持つ半陽的粒子法（PIC）スキームを提示する。この手法により、制限的な時間ステップ制約がなく、安定かつ高精度なシミュレーションが可能となる。Larmor半径が小さい極限において、ガイドセンター運動を漸近的に捉えることができ、渦度形式の非圧縮性Euler方程式を一貫して離散化する。

ABSTRACT

This paper deals with the numerical resolution of the Vlasov-Poissonsystem with a strong external magnetic field by Particle-In-Cell(PIC) methods. In this regime, classical PIC methods are subject tostability constraints on the time and space steps related to the smallLarmor radius and plasma frequency. Here, we propose anasymptotic-preserving PIC scheme which is not subjected to theselimitations. Our approach is based on first and higher order semi-implicit numericalschemes already validated on dissipative systems. Additionally, when the magnitude of the external magneticfield becomes large, this method provides a consistent PICdiscretization of the guiding-center equation, that is, incompressibleEuler equation in vorticity form. We propose several numerical experiments which provide a solid validation of the method and its underlying concepts.

研究の動機と目的

• 強い外部磁場下でプラズマをシミュレートする際、古典的PIC法に見られる不安定性および厳しい時間ステップ制約を解消すること。
• 磁場強度が増加する（Larmor半径が小さくなる）状況でも、安定かつ高精度に保てる数値スキームの開発。
• Larmor半径がゼロに近づく極限において、ガイドセンター近似を漸近的に回復することの保証。
• 強い磁場下において、渦度形式の非圧縮性Euler方程式に対する一貫した粒子ベースの離散化を提供すること。
• 単粒子運動およびdiocotron不安定性の数値シミュレーションを通じて、長時間にわたる安定性と精度を検証すること。

提案手法

• 強い磁場下におけるVlasov-Poisson系に対して、1次、2次、3次の半陽的Runge-Kutta時間積分スキームを提案する。
• 高速な回転運動と遅いガイドセンター運動を分離する時間分割法を導入し、より大きな時間ステップを可能にする。
• 分布関数を粒子ベースで離散化し、Vlasov方程式の特性曲線に従って粒子軌道を進める。
• 比較のため、ガイドセンターモデルにおける電場ポテンシャルを解くために有限差分法を適用する。
• ε → 0 の極限において、初期の精度の順序が保たれることを保証する（εはLarmor周波数の逆数に比例）。
• 形式的に、半陽的スキームがε → 0 の極限において、ガイドセンターモデルと安定かつ一貫性を持つことを示す。

実験結果

リサーチクエスチョン

• RQ1 半陽的時間積分スキームは、強い外部磁場下におけるVlasov-Poisson系のPICシミュレーションを安定化できるか？
• RQ2 磁場強度が増加するに従い、提案スキームはガイドセンターモデルに漸近的に収束するか？
• RQ3 時間ステップ制約がなくとも、diocotron不安定性のような長時間ダイナミクスを正確に捉えることができるか？
• RQ4 εが小さい極限において、エネルギー保存則および電場ノルム‖E(t)‖∞はどのように変化するか？また、このスキームはこれらの性質を保持するか？
• RQ5 非圧縮性Euler方程式の渦度形式をシミュレートする際、粒子法はどの程度の精度と安定性を保つのか？

主な発見

• 1次半陽的スキームは、無条件安定であり、滑らかな電磁場下でガイドセンターモデルと漸近的に一貫している。
• ε = 10⁻² の場合、長時間にわたり相対的エネルギー変動が10⁻³未満に保たれ、優れたエネルギー保存性を示す。
• Vlasov-Poisson系とガイドセンターモデルとの間で、電場ノルム‖E(t)‖∞の時間発展が良好に一致しており、一貫性が確認された。
• 大時間ステップおよび多数の粒子を用いても、t = 120における密度プロファイルは、Vlasov-Poisson系とガイドセンターモデルでほぼ一致した。
• ε → 0 の領域においても、粒子法は過剰な数値的ゆらぎを伴わず、diocotron不安定性の発展を的確に捉えた。
• 3次半陽的スキームは、長時間シミュレーションにおいても高い精度と安定性を維持し、観察された振動や発散は一切なかった。