[論文レビュー] A Hybrid Nodal-Staggered Pseudo-Spectral Electromagnetic Particle-In-Cell Method with Finite-Order Centering

Edoardo Zoni, Remi Lehé+6 more|arXiv (Cornell University)|Jun 23, 2021

Laser-Plasma Interactions and Diagnostics参考文献 61被引用数 4

ひとこと要約

本論文は、有限次数中心化を用いてノード格子とステaggerド格子を組み合わせる、新しいハイブリッド擬スペクトル電磁PIC法を提案する。この手法は、マクスウェル方程式をステaggerドYee格子で解き、高次補間を用いて電流をノード格子にデポジットし、力の収集をノード格子で行うことで、数値的チェレンコフ不安定性（NCI）を低減し、数値的安定性と精度を向上させる。この方法は、レーザー・プラズマ加速や電子-陽電子対生成を含む多様なプラズマ物理学的状況において、純粋なノード格子またはステaggerド格子の手法を凌駕する性能を示す。

ABSTRACT

Electromagnetic particle-in-cell (PIC) codes are widely used to perform computer simulations of a variety of physical systems, including fusion plasmas, astrophysical plasmas, plasma wakefield particle accelerators, and secondary photon sources driven by ultra-intense lasers. In a PIC code, Maxwell's equations are solved on a grid with a numerical method of choice. This article focuses on pseudo-spectral analytical time-domain (PSATD) algorithms and presents a novel hybrid PSATD PIC scheme that combines the respective advantages of standard nodal and staggered methods. The novelty of the hybrid scheme consists in using finite-order centering of grid quantities between nodal and staggered grids, in order to combine the solution of Maxwell's equations on a staggered grid with the deposition of charges and currents and the gathering of electromagnetic forces on a nodal grid. The correctness and performance of the novel hybrid scheme are assessed by means of numerical tests that employ different classes of PSATD equations in a variety of physical scenarios, ranging from the modeling of electron-positron pair creation in vacuum to the simulation of laser-driven and particle beam-driven plasma wakefield acceleration. It is shown that the novel hybrid scheme offers significant numerical and computational advantages, compared to purely nodal or staggered methods, for all the test cases presented.

研究の動機と目的

• 相対論的および低密度プラズマにおいて顕著な数値的チェレンコフ不安定性（NCI）を低減すること。
• 格子のステagger化に起因する高分散性や不安定性に苦しむ純粋なノード格子またはステaggerド格子PIC手法の限界を克服すること。
• 場の計算にステaggerド格子の安定性を活用し、粒子プッシャーの力収集にノード格子の高精度を活かすハイブリッド手法を開発すること。
• Fornberg係数を用いた有限次数補間を導入し、ノード格子とステaggerド格子間で場と電流データを高精度に中心化すること。
• レーザー駆動およびビーム駆動プラズマウェーブフィールド加速、電子-陽電子対生成を含む多様な物理的領域において、性能の向上を実証すること。

提案手法

• 電荷および電流のデポジットをノード格子で行い、マクスウェル方程式をステaggerドYee格子で解くハイブリッドPICサイクルを採用する。
• Fornberg係数に基づく有限次数補間を用い、ノード格子とステaggerド格子間で場および電流データを高精度に転送する。
• 標準的、ガリレオ変換型、平均化ガリレオ変換型の3種類のPSATDスキームをステaggerド格子に適応させ、時間領域における精度を向上させる。
• マクスウェル方程式をフーリエ空間で解析的に積分することで、数値的分散を最小限に抑え、安定性を維持する。
• Fornberg補間係数を用いて、グリッド量のフーリエ変換をモodulateする転送演算子を定義し、補間係数ζFが1から逸脱する度合いを効果的に制御する。
• 高次補間により、ζFが1から逸脱する領域（0 < ζF < 1）を最小限に抑え、運動量保存を確保し、NCIを低減する。

実験結果

リサーチクエスチョン

• RQ1 ノード格子とステaggerド格子を組み合わせたハイブリッドPSATD-PIC法は、純粋なノード格子またはステaggerド格子手法よりも、数値的チェレンコフ不安定性（NCI）をより効果的に低減できるか？
• RQ2 ノード格子とステaggerド格子間の有限次数中心化が、相対論的プラズマシミュレーションにおける電磁場解の安定性および精度にどのように影響するか？
• RQ3 標準的、ガリレオ変換型、平均化ガリレオ変換型のPSATDバリアントの選択が、異なる物理的領域におけるハイブリッドスキームの性能に及ぼす影響はどの程度か？
• RQ4 補間係数ζFがNCI低減に果たす役割は何か？ また、その1からの逸脱度合いは、有限次数中心化の次数にどのように依存するか？
• RQ5 電子-陽電子対生成や高強度レーザー・プラズマ相互作用のような極限状態でも、ハイブリッドスキームは精度と安定性を維持できるか？

主な発見

• ハイブリッドスキームは、純粋なノード格子またはステaggerド格子PIC手法と比較して、特に低密度および相対論的領域において、数値的チェレンコフ不安定性（NCI）を顕著に低減する。
• 高次補間（最大2m = 8）により、補間係数ζFの1からの逸脱が小さくなり、不安定性を引き起こしやすい領域（0 < ζF < 1）が最小限に抑えられる。
• レーザー駆動および粒子ビーム駆動プラズマウェーブフィールド加速を含む多様な物理的状況において、高い精度と安定性を維持する。
• 数値的テストの結果、ハイブリッドスキームは、標準的なノード格子およびステaggerド格子PSATD手法に比べ、場の精度と長期的なシミュレーション安定性において優れていることが示された。
• Fornbergに基づく補間により、ノード格子とステaggerド格子間で高次精度のデータ転送が可能となり、精度と計算効率のバランスを効果的に実現できる。
• 分散解析により、ζFが1に近い場合にNCIが低減されることを確認した。これは、特にハイブリッド構成において、高次補間によりより効果的に達成される。