[論文レビュー] JPIC & How to make a PIC code

H. Wu|arXiv (Cornell University)|Apr 15, 2011

Ionosphere and magnetosphere dynamics被引用数 4

ひとこと要約

この論文は、数値分散フリーな電磁場ソルバーを含む高度な数値アルゴリズムを用いて、レーザー・プラズマ相互作用をシミュレートすることを目的とした、並列処理可能で完全な運動論的粒子-場（PIC）コード、JPICを提示する。本論文は、粒子プッシュ、場の計算、電流および電荷密度の計算、MPIを用いた並列化を含む、カスタムPICコードの開発に関する包括的で段階的なガイドを提供する。応用範囲は相対論的ソリトン形成から、ブーストフレームシミュレーションを用いたkeVオーダーのアトセカンドX線パルス生成まで広がる。

ABSTRACT

Author developed the parallel fully kinetic particle-in-cell (PIC) code JPIC based on updated and advanced algorithms (e.g. numerical-dispersion-free electromagnetic field solver) for simulating laser plasma interactions. Basic technical points and hints of PIC programming and parallel programming by message passing interface (MPI) are reviewed. Most of contents come from Author's notes when writing up JPIC and experiences when using the code to solve different problems. Enough "how-to-do-it" information should help a new beginner to effectively build up his/her own PIC code. General advices on how to use a PIC code are also given.

研究の動機と目的

• 研究者が最初から完全な運動論的PICコードを実装できる包括的で実用的なガイドを提供すること。
• PICシミュレーションにおける電荷保存スキームや数値分散フリーな場のソルバーなどの高度な数値技術を紹介し、実証すること。
• 最適化されたコード設計と並列化を通じて、相対論的自己焦点化、ウェーブフィールド加速、およびコherentなトンプソン散乱といった複雑なレーザー・プラズマ現象の効率的シミュレーションを可能にすること。
• ブーストフレーム法などのフレーム変換技術を用いて、極端な状態（例：a₀ > 100、keV X線生成）へのPICシミュレーションの適用範囲を拡大すること。
• 数値精度とアルゴリズムの安定性に配慮しつつ、将来の3次元PICコードおよびGPUアクセcelerationの開発を支援すること。

提案手法

• コードの簡略化と異なる物理的スケールへの普遍性を確保するため、レーザー関連の単位で正規化されたマクスウェル方程式とローレンツ力方程式を用いる。
• 高k数でのシミュレーションにおける位相速度誤差を排除するため、数値分散フリー（NDF）電磁場ソルバーを採用し、高周波数および短波長現象の精度を向上させる。
• ポアソン方程式を解かずに全電荷保存を維持する電荷保存スキーム（CCS）を実装し、効率的な並列化を可能にするとともに、グローバル通信のボトルネックを回避する。
• 粒子形状関数を用いた粒子-場（PIC）アルゴリズムを採用し、場の補間と電流のデポジションを実行する。二重精度浮動小数点演算を用いることで、数値的精度を維持する。
• 複数プロセッサ間でのドメイン分割と並列化にMessage Passing Interface（MPI）を用い、移動ウィンドウおよび斜め入射に対応する。
• 長距離かつ高相対論的プロセス（例：keV X線生成）のシミュレーションを可能にするために、共動する慣性系に変換するブーストフレーム技術を適用し、シミュレーションコストを数個のオーダー低下させる。

実験結果

リサーチクエスチョン

• RQ1 初心者が、場の計算、粒子プッシュ、電流デポジションを含め、最初から完全な運動論的PICコードを効果的に実装する方法は何か？
• RQ2 高強度レーザー・プラズマ領域において、高精度および安定性を達成するために不可欠な数値アルゴリズム（例：NDF、CCS）は何か？
• RQ3 1次元および2次元PICコードにおけるMPIを用いた並列化を、大規模システムおよび長時間シミュレーションにスケーリングするために効果的に実装する方法は何か？
• RQ4 ブーストPIC技術は、相対論的および超短パルス現象のシミュレーションにおいて、計算コストをどのように大幅に削減できるか？
• RQ5 高精度および高動的レンジ問題を対象とした3次元シミュレーションおよびGPUアクセcelerationへのPICコードの拡張において、主な課題とその解決策は何か？

主な発見

• JPICに数値分散フリー（NDF）電磁場ソルバーを適用することで、位相速度誤差が完全に排除され、高kモードおよび短波長現象の正確なシミュレーションが可能になった。
• 電荷保存スキーム（CCS）により、ポアソン方程式を解かずに全電荷保存が維持され、効率的な並列化が可能になり、グローバル通信のボトルネックを回避できるようになった。
• ブーストPIC技術により、コherentなトンプソン散乱（CTS）におけるkeV X線生成のシミュレーション時間が最大2桁短縮され、従来では不可能とされていたシミュレーションが可能になった。
• ブーストフレームにおける非線形CTSのシミュレーション結果は、解析理論と良好に一致し、非線形ドップラー効果、出力の飽和、パルスの鋭利化といった新規効果が明らかになった。
• コードは相対論的ソリトン形成（a₀=1, ne=0.5nc）、ビーム駆動ウェーブフィールドにおけるマイクロバンチング、およびマルチ粒子ウェイトを用いた一様な相対論的電子シート生成といった複雑な現象を成功裏にモデル化した。
• 今後の計画として、NDFXYZソルバーおよび3次元ジグザグCCSを用いた3次元実装、およびGPUアクセcelerationの導入を予定しており、1次元シミュレーションにおける低メモリ要件を考慮すると、1つのGPUで10倍の速度向上が見込まれる。