[論文レビュー] Introducing ZEUS-MP: A 3D, Parallel, Multiphysics Code for Astrophysical Fluid Dynamics

Michael L. Norman|arXiv (Cornell University)|May 5, 2000

Geophysics and Gravity Measurements被引用数 6

ひとこと要約

ZEUS-MP は、ドメイン分割と MPI ベースのメッセージパッシングを用いて、理想流体力学、MHD、自己重力の実装を含む、アストロフィジカル流体力学のための巨大並列的3次元マルチフィジックスコードである。256プロセッサまででほぼ理想に近い強化スケーリングを達成し、超音速乱流が存在するにもかかわらず、磁化された自己重力分子雲における重力収縮を示しており、磁化が過剰な環境における星形成の研究にその有用性を裏付けている。

ABSTRACT

We describe ZEUS-MP: a Multi-Physics, Massively-Parallel, Message-Passing code for astrophysical fluid dynamics simulations in 3 dimensions. ZEUS-MP is a follow-on to the sequential ZEUS-2D and ZEUS-3D codes developed and disseminated by the Laboratory for Computational Astrophysics (lca.ncsa.uiuc.edu) at NCSA. V1.0 released 1/1/2000 includes the following physics modules: ideal hydrodynamics, ideal MHD, and self-gravity. Future releases will include flux-limited radiation diffusion, thermal heat conduction, two-temperature plasma, and heating and cooling functions. The covariant equations are cast on a moving Eulerian grid with Cartesian, cylindrical, and spherical polar coordinates currently supported. Parallelization is done by domain decomposition and implemented in F77 and MPI. The code is portable across a wide range of platforms from networks of workstations to massively parallel processors. Some parallel performance results are presented as well as an application to turbulent star formation.

研究の動機と目的

• 3次元アストロフィジカル流体力学シミュレーションのためのポータブルで、巨大並列的かつ拡張可能なコードの開発。
• 順次処理に限界のある ZEUS-2D および ZEUS-3D の課題を克服し、分散メモリシステムにおけるスケーラブルなシミュレーションを可能にすること。
• 流体力学、MHD、自己重力に加え、将来の追加として放射拡散や熱伝導を含む、複数の物理モジュールのサポート。
• 磁化され、自己重力を持つ雲における乱流星形成のような、複雑な天体物理学的現象の高精度なシミュレーションの実現。
• 天体物理学コミュニティが広く利用可能で、良好にドキュメント化され、拡張可能なシミュレーションフレームワークの提供。

提案手法

• プロセッサ間通信には MPI を用いた単一プログラム・マルチデータ（SPMD）モデルを採用し、さまざまな HPC プラットフォームへのポータビリティを実現。
• 計算グリッドをプロセッサ間で分割するドメイン分割を用い、性能最適化のため計算と通信をオーバラップさせる。
• 直交座標、円柱座標、球座標を用い、幾何学的柔軟性を確保するための共変表現を用いて、移動するEuler的グリッド上で方程式を解く。
• 流体力学モジュールは、ZEUS-2D で採用された2次精度の単調性保証の輸送スキームを3次元に拡張し、衝撃の捕捉における安定性と精度を確保。
• MHD モジュールは、磁場の発散がゼロとなるように保証するための「制約輸送」と「特徴線法」を組み合わせた MOC-CT アルゴリズムを実装し、波の分解能を向上。
• 自己重力モジュールは、多重グリッド法と FFT を用いたソルバーでポアソン方程式を解き、後者は分散メモリシステムにおける性能最適化を図っている。

実験結果

リサーチクエスチョン

• RQ1 巨大並列的マルチフィジックスコードは、3次元アストロフィジカル流体力学シミュレーションにおいて、大規模 HPC システムでほぼ理想に近い強化スケーリングを達成できるか？
• RQ2 質量・磁束比が臨界値を超える磁化された分子雲において、超音速乱流は重力収縮をどの程度妨げるか？
• RQ3 ZEUS-MP は、乱流、磁場、自己重力の相互作用を、高密度で重力的に束縛されたコア形成の文脈でどの程度正確にシミュレートできるか？
• RQ4 コードのモジュラー設計は、フラックス制限付き放射拡散や熱伝導の今後の拡張をサポートできるか？
• RQ5 実際の天体物理学ワークロードにおいて、異なるネットワークおよびプロセッサ構成での ZEUS-MP の性能スケーリングはどの程度か？

主な発見

• ZEUS-MP は、SGI/Cray Origin2000 システム上で、256プロセッサまでで、$256^3$ および $512^3$ グリッド解像度の流体力学および MHD シミュレーションにおいて、ほぼ理想に近い強化スケーリングを達成した。
• MIPS R10000 プロセッサ搭載の単一ノードでは約 100 MFlop/s の性能を示し、ネットワークおよびシステムの特性に強く依存していた。
• $M/M_{cr} = 1.1$ の乱流的で磁化され、自己重力を持つ分子雲シミュレーションにおいて、超音速乱流（マッハ数 $M_{turb} = 5$）が存在するにもかかわらず重力収縮が進行した。これは、過剰な磁化雲では乱流そのものが収縮を防げないことを示している。
• フィラメント状および平板状の凝縮が、平均磁場に平行および垂直に形成され、高密度コアはフィラメントの中心部に形成された。これは観測された星形成の形態と整合的である。
• 臨界質量・磁束比が、乱流環境における収縮の予測に十分ではないことがシミュレーションで示された。$M/M_{cr} = 1.1 > 1$ であっても収縮が進行したため、従来の仮定に疑問を呈するものである。
• コードのモジュラーアーキテクチャは、複数の物理モジュールの統合を成功裏にサポートしており、放射拡散および熱伝導の今後のリリースに向けた計画が進行中である。