[論文レビュー] An Energy Conserving Parallel Hybrid Plasma Solver

Mats Holmström|arXiv (Cornell University)|Oct 15, 2010

Magnetic confinement fusion research参考文献 5被引用数 4

ひとこと要約

本論文は、FLASHフレームワーク内において、セル中心有限差分法と予測子補正型ルンゲ・クッタ法（サイクル付きルンゲ・クッタ法）を用いた並列的でエネルギー保存型のハイブリッドプラズマソルバを提案する。50回のガロア周期にわたり0.1%未満のエネルギー誤差を示し、大規模クラスタ上で優れた弱スケーリング性能を発揮する。フィールド平滑化を施さないまま、従来のソルバーよりも安定性と精度に優れている。

ABSTRACT

We investigate the performance of a hybrid plasma solver on the test problem of an ion beam. The parallel solver is based on cell centered finite differences in space, and a predictor-corrector leapfrog scheme in time. The implementation is done in the FLASH software framework. It is shown that the solver conserves energy well over time, and that the parallelization is efficient (it exhibits weak scaling).

研究の動機と目的

• 衝突なしプラズマシミュレーション用の数値的に安定でエネルギー保存型のハイブリッドプラズマソルバの開発。
• 大規模シミュレーションを対象として、並列処理フレームワークFLASH内での効率的実装。
• フィールド平滑化を施さずに、長時間スケールでのエネルギー保存性と数値的安定性の検証。
• 高性能コンピューティングクラスタ上での弱スケーリング性能の実証。
• イオンビームおよび宇宙プラズマ応用に適した強固で再利用可能なソルバの提供。

提案手法

• 空間微分に2次精度を持つセル中心有限差分法を用いる。
• 時間積分に予測子補正型ルンゲ・クッタ法（サイクル付きルンゲ・クッタ法）を採用し、フィールド更新を部分的にサイクリングする。
• マクロ粒子の電荷をグリッド上にCIC（Cloud-in-Cell）線形重み付け法で配置する。
• 力計算のため、線形補間によりグリッド上の場を粒子位置に補間する。
• ハイブリッドモデルとして、イオンを粒子として、電子を無限小質量の流体として扱い、ローレンツ力とファラデーの法則を解く。
• MPIを用いた並列化とHDF5を用いた並列I/Oにより、FLASHフレームワーク上で実行する。

実験結果

リサーチクエスチョン

• RQ1 有限差分ハイブリッドソルバは、50回のガロア周期にわたり0.1%未満のエネルギー保存性を達成できるか？
• RQ2 フィールド平滑化を施さずに、長時間積分においても数値的安定性を維持できるか？
• RQ3 大規模HPCクラスタ上での弱スケーリングモードにおいて、どの程度の性能を示すか？
• RQ4 これまでに発表されたソルバーよりもエネルギー保存性の性能が優れているか？
• RQ5 FLASHフレームワークは、粒子および場の計算を伴うハイブリッドプラズマシミュレーションを、効率的にサポートできるか？

主な発見

• t = 86 s（約50 Ω⁻¹）で、合計エネルギー誤差が0.1%未満に抑えられ、Matthews（1994）が報告した11%の誤差よりも著しく優れている。
• 2次元でもエネルギー保存性が優れており、t = 173 s（約100 Ω⁻¹）で合計相対エネルギー増加が0.004%未満にとどまる。
• 強い弱スケーリングを示す：672コア上で異なる問題サイズにおいて、1ステップあたりのウォールクロック時間が約70秒でほぼ一定である。
• フィールド平滑化を必要とせず、長時間のシミュレーションにおいても数値的安定性が維持され、従来の実装とは異なり安定している。
• FLASHにおける実装により、グリッド演算、ゴーストセル更新、I/Oを含む既存の並列インfraストラクチャを効率的に活用できる。
• t = 34.6 sにおける速度空間分布は、Winske & Leroy（1984）の既存結果とよく一致しており、物理的正確性が裏付けられている。