[論文レビュー] Overcoming timestep limitations in boosted-frame Particle-In-Cell simulations of plasma-based acceleration

Olga V. Shapoval, Remi Lehé+4 more|arXiv (Cornell University)|Apr 28, 2021

Laser-Plasma Interactions and Diagnostics参考文献 29被引用数 7

ひとこと要約

本論文は、プラズマベース加速器のブーストフレーム粒子-場（PIC）シミュレーションにおけるタイムステップ制限を克服するため、平均化されたガリレオ変換PSATDアルゴリズムを提案する。時間領域で電磁場をフーリエ空間で解析的に平均化することにより、c∆t ≫ ∆x, ∆y となる大きなタイムステップを安定に使用可能にし、数値的不安定性を引き起こさない。これにより、高分解能の横方向シミュレーションにおいて、ブーストフレーム手法の計算高速化の潜在的効果を完全に解き放つことが可能になる。

ABSTRACT

Explicit electromagnetic Particle-In-Cell (PIC) codes are typically limited by the Courant- Friedrichs-Lewy (CFL) condition, which implies that the timestep multiplied by the speed of light must be smaller than the smallest cell size. In the case of boosted-frame PIC simulations of plasma-based acceleration, this limitation can be a major hinderance as the cells are often very elongated along the longitudinal direction and the timestep is thus limited by the small, transverse cell size. This entails many small-timestep PIC iterations, and can limit the potential speed-up of the boosted-frame technique. Here, by using a CFL-free analytical spectral solver, and by mitigating additional numerical instabilities that arise at large timestep, we show that it is possible to overcome traditional limitations on the timestep and thereby realize the full potential of the boosted-frame technique over a much wider range of parameters.

研究の動機と目的

• 横方向セルサイズの制約によりタイムステップ制限が生じるが、縦方向分解能は緩和されているブーストフレームPICシミュレーションにおけるタイムステップ制限を解消すること。
• PSATDに基づくブーストフレームシミュレーションにおいて、大きなタイムステップで再発する数値チェレンコフ不安定性（NCI）を解消すること。
• PSATDシミュレーションにおいて、c∆t ≫ ∆x, ∆y となる大きなタイムステップを用いながらも、シミュレーションの正確性と安定性を維持すること。
• 高横方向分解能を要するプラズマ加速器構成の幅広い範囲にわたり、ブーストフレーム手法の計算高速化を拡張すること。

提案手法

• 時間間隔 [(n−1/2)∆t, (n+1/2)∆t] において、フーリエ空間で電磁場を解析的に平均化する、平均化されたガリレオ変換PSATDアルゴリズムを導入する。
• タイムステップ内での電流が一定であると仮定し、マクスウェル方程式の解析解を用いて、時間平均場 ⟨ˆE⟩n および ⟨ˆB⟩n をPSATDフレームワークに適応して計算する。
• ガリレオ変換のドリフトとスペクトル精度を反映する複素係数（Ψ1, Ψ2, A1, A2, Cρ）を含む、新たな場の更新式を導出する。
• 標準的なPICループを変更し、粒子の運搬に時間平均場を使用することで、場の収集と運動量更新の整合性を保証する。
• フーリエ空間におけるp次元離散化を採用し、グリッド解像度が変化しても高い正確性と安定性を維持する。
• CPUおよびGPUプラットフォームをサポートするオープンソースのWarpXコードにアルゴリズムを実装する。

実験結果

リサーチクエスチョン

• RQ1 c∆t ≫ ∆x, ∆y となる大きなタイムステップを、ブーストフレームPSATDシミュレーションで安全に使用できるか？（数値的不安定性が発生しないか？）
• RQ2 CFLフリーであるにもかかわらず、c∆t > ∆x, ∆y の場合に標準的なガリレオ変換PSATDアルゴリズムが不安定化するのはなぜか？
• RQ3 時間領域でフーリエ空間での時間平均化を用いることで、場の収集プロセスを変更することで不安定性を軽減できるか？
• RQ4 本新アルゴリズムは、現実のプラズマ加速器シナリオにおいて、物理的整合性と計算高速化の両方をどの程度維持できるか？
• RQ5 本手法は、高横方向分解能を要するが高時間分解能を必要としない他のPICシミュレーションにも一般化可能か？

主な発見

• 平均化されたガリレオ変換PSATDアルゴリズムは、c∆t ≫ ∆x, ∆y となる大きなタイムステップでも、再発する数値チェレンコフ不安定性（NCI）を効果的に抑制し、安定なシミュレーションを可能にする。
• c∆t = ∆z（横方向セルサイズよりもはるかに大きい）のシミュレーションでは、NCIモードの成長が認められず、標準的なガリレオ変換PSATDアルゴリズムとは対照的である。
• 2次元均一プラズマシミュレーションでは、新しいアルゴリズムはc∆tが∆zまで安定に保たれる一方、標準的手法はc∆t > ∆x で不安定化する。
• 2次元レーザー・ウェーブフィールド加速（LWFA）では、平均化PSATD法が大きなタイムステップでも正しいウェーブフィールド構造と電子ビームダイナミクスを再現し、参照シミュレーションと一致する。
• 3次元プラズマウェーブフィールド加速（PWFA）では、c∆t ≫ ∆x, ∆y となる安定で高速なシミュレーションが可能となり、顕著な計算高速化が達成された。
• アルゴリズムはWarpXコードに実装され、適応メッシュリファインメントおよびCPU/GPUハイパフォーマンスコンピューティング環境への対応が確認された。