[論文レビュー] Runko: Modern multi-physics toolbox for simulating plasma

Joonas Nättilä|arXiv (Cornell University)|Jun 14, 2019

Magnetic confinement fusion research参考文献 1被引用数 7

ひとこと要約

Runkoは、マルチフィジックスおよびマルチスケールモデリングをサポートする、オープンソースでハイブリッドなC++/Pythonフレームワークであり、天体物理学的プラズマのシミュレーションを可能にします。高性能を実現する現代的なC++14と、ハイレベルな制御を可能にするPythonを組み合わせ、相対論的PICモジュールを通じて、磁場支配のペアプラズマにおける相対論的乱流の崩壊が、熱的乱流カスケードと別個の非熱的粒子集団を生じさせ、エネルギーの吸収源として機能することを示しています。

ABSTRACT

Runko is a new open-source plasma simulation framework implemented in C++ and Python. It is designed to function as an easy-to-extend general toolbox for simulating astrophysical plasmas with different theoretical and numerical models. Computationally intensive low-level kernels are written in modern C++14 taking advantage of polymorphic classes, multiple inheritance, and template metaprogramming. High-level functionality is operated with Python3 scripts. This hybrid program design ensures fast code and ease of use. The framework has a modular object-oriented design that allow the user to easily add new numerical algorithms to the system. The code can be run on various computing platforms ranging from laptops (shared-memory systems) to massively parallel supercomputer architectures (distributed-memory systems). The framework also supports heterogeneous multi-physics simulations in which different physical solvers can be combined and run simultaneously. Here we report on the first results from the framework's relativistic particle-in-cell (PIC) module. Using the PIC module, we simulate decaying relativistic kinetic turbulence in suddenly stirred magnetically-dominated pair plasma. We show that the resulting particle distribution can be separated into a thermal part that forms the turbulent cascade and into a separate decoupled non-thermal particle population that acts as an energy sink for the system.

研究の動機と目的

• 複雑な天体物理学的プラズマ系をモデリングするための柔軟で拡張可能かつ高性能なシミュレーションフレームワークの開発。
• 異なる物理的ソルバーが単一の計算環境内で共存・結合可能であることを通じて、マルチフィジックスプラズマ現象のシミュレーションという課題に取り組む。
• 共有メモリおよび分散メモリアーキテクチャの両方をサポートし、ラップトップからスーパーコンピュータまでにおけるポータビリティとスケーラビリティを確保。
• 磁場支配のペアプラズマにおける相対論的キネティック乱流の研究を可能にし、特に非熱的粒子集団の出現を明らかにすること。

提案手法

• フレームワークはハイブリッドアーキテクチャを採用：低レベルの計算カーネルは、パフォーマンスと拡張性を実現するためのポリモービッククラス、多重継承、テンプレートメタプログラミングを用いた現代的なC++14で実装。
• ハイレベルな制御とワークフロー管理はPython3スクリプトで処理され、使いやすさと迅速なプロトタイピングを実現。
• モジュラーでオブジェクト指向の設計により、ユーザーは簡単に新しい数値アルゴリズムや物理モデルをシミュレーションパイプラインに統合可能。
• 異なるソルバー（例：相対論的PICモジュール）を併用して並列実行できる仕組みにより、異種マルチフィジックスシミュレーションをサポート。
• 相対論的PICモジュールは、粒子-場法を用いて磁場がかかるプラズマ内での相対論的電子および陽電子のダイナミクスをシミュレート。
• ラップトップからスーパーコンピュータまで多様なコンピューティングプラットフォームに展開され、共有メモリおよび分散メモリ並列処理を活用。

実験結果

リサーチクエスチョン

• RQ1 急激に攪拌された磁場支配のペアプラズマにおいて、相対論的キネティック乱流はどのように発展するか？
• RQ2 このような乱流系に現れる、明確に分離された粒子集団の成分は何か？
• RQ3 乱流カスケードの過程で、エネルギーが熱的および非熱的粒子集団の間でどのように分配されるか？
• RQ4 非熱的粒子集団は、相対論的プラズマ系において安定なエネルギー吸収源として機能できるか？
• RQ5 Runkoフレームワークは、複雑なプラズマ現象のスケーラブルなマルチフィジックスシミュレーションをどの程度サポートできるか？

主な発見

• シミュレーションにより、乱流エネルギーカスケードに参加する熱的粒子集団と、別個の非熱的集団との明確な分離が確認された。
• 非熱的粒子集団は独立に形成され、エネルギーの主な吸収源として機能し、乱流カスケードからのエネルギーを吸収している。
• 熱的成分は、発達した乱流カスケードを示す、コルモゴロフに類似したエネルギースペクトルを示している。
• 非熱的集団は空間的および運動論的に熱的成分から分離されており、別個の加速機構が存在すると示唆される。
• フレームワークは、ラップトップからスーパーコンピュータまで多様なハードウェアプラットフォームでスケーリングに成功し、パフォーマンスと使いやすさを維持している。