[論文レビュー] A New Public Release of the GIZMO Code

Philip F. Hopkins|arXiv (Cornell University)|Dec 4, 2017

Astrophysics and Star Formation Studies参考文献 2被引用数 16

ひとこと要約

この論文は、GIZMOコードの大幅な公開リリースを発表している。GIZMOは、適応メッシュフリー有限体積流体動力学、SPH、MHD、異方性熱伝導/粘性、星およびブラックホール形成、ダスト-ガスダイナミクス、弾性/塑性固体、および任意の状態方程式をサポートする、大規模並列でマルチフィジックスを扱う天体物理学的シミュレーションコードである。このリリースにより、100万スレッドを超える規模でほぼ理想に近い強拡張および弱拡張が実現可能となり、広範なドキュメンテーション、テスト問題、およびGADGETツールとの互換性が提供される。

ABSTRACT

We describe a major update to the public GIZMO code. GIZMO has been used in simulations of cosmology; galaxy and star formation and evolution; black hole accretion and feedback; proto-stellar disk dynamics and planet formation; fluid dynamics and plasma physics; dust-gas dynamics; giant impacts and solid-body interactions; collisionless gravitational dynamics; and more. This release of the public code supports: hydrodynamics (using various mesh-free finite-volume Godunov methods or SPH), ideal and non-ideal MHD, anisotropic conduction and viscosity, radiative cooling and chemistry, star and black hole formation and feedback, sink particles, dust-gas (aero)-dynamics (with or without magnetic fields), elastic/plastic dynamics, arbitrary (gas, stellar, degenerate, solid/liquid material) equations of state, passive scalar/turbulent diffusion, large-eddy and shearing boxes, self-gravity with fully-adaptive force softenings, arbitrary cosmological expansion, and on-the-fly group-finding. It is massively-parallel with hybrid MPI+OpenMP scaling verified up to >1 million threads. The code is extensively documented, with test problems and tutorials provided for these different physics modules.

研究の動機と目的

• 宇宙論、銀河形成、星・惑星形成、プラズマ物理学をカバーする多様な天体物理学的問題に適した包括的で公開可で、かつ非常にスケーラブルなシミュレーションコードを提供すること。
• 自己重力、フィードバック、多相流体ダイナミクスを含む高度な物理モジュールを用いて、大規模シミュレーションを柔軟でモジュラーなフレームワークで実行可能にすること。
• 広範なドキュメンテーション、チュートリアル、テスト問題、および既存のGADGETワークフローとデータ形式との互換性を提供することで、研究者による広範なアクセス性と再現可能性を確保すること。
• エクサスケールクラスのスーパーコンピュータ上で、10億個を超える解像度要素を含むシミュレーションをサポートする、ほぼ理想に近い強拡張および弱拡張を達成すること。
• パubblicおよびプライベート開発ブランチを分離することで、コードの拡張性と将来対応性を維持し、特許または不安定な物理モジュールの責任ある共有を可能にすること。

提案手法

• GIZMOは、適応的ドメイン分割と階層的タイムステップを用いた、ハイブリッドMPI+OpenMP並列化戦略を採用しており、大規模スーパーコンピュータでの効率的スケーリングを実現している。
• コードは複数の数値法をサポートする：メッシュフリー有限体積ゴドノフ法、SPH、固定直交座標グリッドで、ユーザー定義のメッシュ運動（例：せん断、拡張、回転）が可能である。
• 任意のラグランジュ・オイラー流体動力学を実装し、完全に適応可能な重力ソフトニングと、任意の状態方程式（固体用のTillotson型、縮重物質用のHelmholtz型を含む）をサポートしている。
• 主な物理モジュールには、オンザフライでのグループ検出、シンク粒子形成、自己一貫的な星およびブラックホールフィードバック、異方性熱伝導および粘性、Smagorinskiモデルによる乱流拡散が含まれる。
• ダスト-ガスダイナミクスは、柔軟な抵抗則（Epstein、Stokes、Coulomb）とローレンツ力に対応し、ガスへの反作用が運動量方程式に含まれている。
• 宇宙論的シミュレーションでは、標準でない膨張歴史、動的ダークエネルギー、任意の境界条件（周期的、開口、反射）をサポートし、オンザフライでのパワースペクトルとホールファインディングが可能である。

実験結果

リサーチクエスチョン

• RQ1 1つのシミュレーションコードが、宇宙論から惑星形成まで広範な天体物理学的問題を効率的にサポートできるか？
• RQ2 流体動力学、重力、磁場、フィードバックを含むマルチフィジックスシミュレーションで、どの程度の並列拡張が達成可能か？
• RQ3 適応的メッシュフリー法とラグランジュ粒子追跡を組み合わせることで、複雑な流体および固体ダイナミクスをどのようにモデル化できるか？
• RQ4 GIZMOは、強い磁場乱流、異方性熱伝導、多相性星間媒体のシミュレーションにおいて、性能と正確さをどの程度達成できるか？
• RQ5 モジュラーで拡張可能な物理モジュールを、1つのスケーラブルで十分にドキュメント化されたシミュレーションフレームワークに統合する方法は何か？

主な発見

• GIZMOは、宇宙論的シミュレーションにおいて100万スレッドを超える規模でほぼ理想に近い弱拡張を達成しており、2〜10,000 Mpc³の体積範囲で性能が安定している。
• 1億5,000万個の粒子を用いた銀河系質量および矮星銀河ズームインシミュレーションにおいて、強拡張は25万CPUおよび1万6千ノードにわたり、ほぼ理想に近い効率を維持している。
• コードは、Orszag-Tang渦、磁気熱不安定性、ケルビン＝ヘルムホルツ不安定性といった複雑な現象を、高い正確性と解像度で正確に再現している。
• GIZMOは、乱流的で磁場のある雲におけるダスト-ガスダイナミクスを、Epsteinおよびローレンツ力下での粒子運動を含めて正確にモデル化できており、ガスへの反作用も含んでいる。
• 重力的捕獲とフィードバック駆動型噴出により、ブラックホール降着を解像度を保って再現でき、正確にエネルギーと運動量をISMに供給している。
• Tillotson型状態方程式を用いた弾性および塑性ダイナミクスの導入により、物質混合と降伏基準を伴う巨大衝突や惑星体衝突の現実的なシミュレーションが可能になった。