[論文レビュー] Conservative Finite Element Time-Domain PIC Algorithm for Relativistic Maxwell-Vlasov Equations on Irregular Meshes

D. Y. Na, Haksu Moon+2 more|arXiv (Cornell University)|Jul 13, 2017

Laser-Plasma Interactions and Diagnostics被引用数 1

ひとこと要約

この論文は、非構造化メッシュ上で相対論的 Maxwell-Vlasov 方程式を解くための保存的有限要素時域粒子-場（FETD-PIC）アルゴリズムを提示する。Boris、Vay、Higuera-Cary の3つの相対論的粒子プッシュ方式を統合することで、電荷保存を厳密に保ちつつ、サイクロトロン運動、ブーストフレームにおける高調波振動、ペアプラズマにおける相対論的 Bernstein モードを含む、相対論的プラズマダイナミクスの高精度なシミュレーションが可能になる。

ABSTRACT

In many problems involving particle accelerators and relativistic plasmas, the accurate modeling of relativistic particle motion is essential for accurate physical predictions. Here, we extend a charge-conserving finite element time-domain (FETD) particle-in-cell (PIC) algorithm for the time-dependent Maxwell-Vlasov equations on irregular (unstructured) meshes to the relativistic regime by implementing and comparing three particle pushers: (relativistic) Boris, Vay, and Higuera-Cary. We illustrate the application of the proposed relativistic FETD-PIC algorithm for the analysis of particle cyclotron motion at relativistic speeds, harmonic particle oscillation in the Lorentz-boosted frame, and relativistic Bernstein modes in magnetized charge-neutral (pair) plasmas.

研究の動機と目的

• 粒子加速器や核融合プラズマに一般的に見られる複雑で不規則な幾何形状における相対論的粒子運動の高精度なモデリングを可能にすること。
• 電荷保存型有限要素法を相対論的領域に拡張し、電磁気的プラズマシミュレーションにおける精度を向上させること。
• 有限要素 PIC フレームワーク内での3つの相対論的粒子プッシュ方式（Boris、Vay、Higuera-Cary）の性能を評価・比較すること。
• サイクロトロン運動や Bernstein モードといった、アルゴリズムの能力を示す主要な相対論的プラズマ現象のシミュレーションを実証すること。

提案手法

• 相対論的運動量とローレンツ力項を組み込むことで、電荷保存型有限要素時域（FETD）フレームワークを相対論的粒子運動に適応させる。
• FETD-PIC フレームワーク内に、Boris、Vay、Higuera-Cary の3つの相対論的粒子プッシュ方式を実装し、時間発展に伴う粒子の位置と速度を進める。
• 複雑なプラズマおよび加速器構成をモデル化するための幾何的柔軟性を実現するため、非構造化（不規則）メッシュの使用。
• 有限要素の変分定式化から導かれる離散的連続の式を用いて、電荷保存を強制的に満たす。
• 粒子-場法を介して Maxwell 方程式と Vlasov 方程式を結合し、非構造化グリッド上での場の補間と電流のデポジションを実行する。
• 相対論的サイクロトロン運動、ブーストフレーム振動、ペアプラズマモードのベンチマークシミュレーションを通じて、アルゴリズムの妥当性を検証する。

実験結果

リサーチクエスチョン

• RQ1 FETD-PIC アルゴリズムは、非構造化メッシュ上での相対論的領域において、どのように電荷保存を維持しているか？
• RQ2 Boris、Vay、Higuera-Cary の3つのプッシュ方式は、相対論的プラズマシミュレーションにおいて、相対的な精度と安定性の性能でどのように差がつくか？
• RQ3 アルゴリズムは、ローレンツブーストフレームにおける相対論的サイクロトロン運動と高調波振動を正確に捉えられるか？
• RQ4 アルゴリズムは、磁場を有する電荷中性ペアプラズマにおける相対論的 Bernstein モードをどのようにシミュレートするか？
• RQ5 メッシュの不規則性は、相対論的 FETD-PIC スキームの精度および保存性にどの程度影響を与えるか？

主な発見

• 離散的変分定式化のおかげで、非構造化メッシュ上でもすべてのシミュレーションで厳密な電荷保存が達成された。
• 高相対論的領域では、標準的な Boris 法や Higuera-Cary 法に比べ、Vay プッシュ方式が優れた精度と安定性を示した。
• 相対論的ローレンツ力下での理論的予測と一致する粒子軌道が正確に再現された。
• ローレンツブーストフレームにおける高調波的粒子振動も正しくシミュレートされ、異なる参照枠での一貫性が裏付けられた。
• 磁場を有するペアプラズマにおける相対論的 Bernstein モードが正常に励起され、解像された。これにより、複雑な波動現象のシミュレーション能力が確認された。
• 非構造化メッシュの使用により、保存性や精度を損なうことなく、複雑な幾何形状の高精度なモデル化が可能になった。