[論文レビュー] SWIFT: A modern highly-parallel gravity and smoothed particle hydrodynamics solver for astrophysical and cosmological applications

Matthieu Schaller, Josh Borrow+19 more|arXiv (Cornell University)|May 22, 2023

Galaxies: Formation, Evolution, Phenomena参考文献 220被引用数 23

ひとこと要約

SWIFT は、共有メモリおよび分散メモリ並列処理を組み合わせた、完全にオープンソースで高並列なタスクベースの天体物理学シミュレーションコードであり、重力、滑らかめの粒子流体力学（SPH）、および宇宙論的物理学を統合しています。高速多重極法（FMM）および粒子メッシュ重力ソルバーを用いて、約3000億個の粒子を含む宇宙論的シミュレーションにおいて優れた弱スケーリング性能を達成しており、複数のSPHバージョンと、銀河形成および惑星物理学のためのモジュラーな小スケールモデルを備えています。

ABSTRACT

Numerical simulations have become one of the key tools used by theorists in all the fields of astrophysics and cosmology. The development of modern tools that target the largest existing computing systems and exploit state-of-the-art numerical methods and algorithms is thus crucial. In this paper, we introduce the fully open-source highly-parallel, versatile, and modular coupled hydrodynamics, gravity, cosmology, and galaxy-formation code SWIFT. The software package exploits hybrid shared- and distributed-memory task-based parallelism, asynchronous communications, and domain-decomposition algorithms based on balancing the workload, rather than the data, to efficiently exploit modern high-performance computing cluster architectures. Gravity is solved for using a fast-multipole-method, optionally coupled to a particle mesh solver in Fourier space to handle periodic volumes. For gas evolution, multiple modern flavours of Smoothed Particle Hydrodynamics are implemented. SWIFT also evolves neutrinos using a state-of-the-art particle-based method. Two complementary networks of sub-grid models for galaxy formation as well as extensions to simulate planetary physics are also released as part of the code. An extensive set of output options, including snapshots, light-cones, power spectra, and a coupling to structure finders are also included. We describe the overall code architecture, summarise the consistency and accuracy tests that were performed, and demonstrate the excellent weak-scaling performance of the code using a representative cosmological hydrodynamical problem with $\approx$$300$ billion particles. The code is released to the community alongside extensive documentation for both users and developers, a large selection of example test problems, and a suite of tools to aid in the analysis of large simulations run with SWIFT.

研究の動機と目的

• 大規模な天体物理学的および宇宙論的問題を扱うためのスケーラブルでモジュラーかつオープンソースのシミュレーションコードの開発。
• タスクベース並列処理とワークロードバランスの取れたドメイン分割を活用し、現代のハイパフォーマンスコンピューティングアーキテクチャを効率的に活用すること。
• 拡張可能な小スケール物理モデルを備えた、重力、流体力学、ニュートリノ、銀河形成の正確で柔軟なモデル化を支援すること。
• 公開リリース、包括的なドキュメンテーション、共有の例題問題を通じて、再現可能性とコミュニティ連携を促進すること。

提案手法

• 現代のHPCクラスタで効率的にスケーリングできるように、非同期通信を伴うハイブリッド共有メモリおよび分散メモリのタスクベース並列処理を採用。
• 重力計算には高速多重極法（FMM）を用い、周期的境界条件の場合は、フーリエ空間の粒子メッシュソルバーと optionally 組み合わせる。
• Sphenix を含む複数の最新のSPHバージョンを実装し、天体環境における正確なガス力学を再現。
• 最新の粒子ベース手法を用いてニュートリノの進化をモデル化し、大規模構造シミュレーションにおけるその宇宙論的影響を反映。
• 銀河形成のための2つの補完的な小スケールモデルを統合し、惑星系物理学への拡張をサポート。
• スナップショット、光円錐、パワースペクトルなど、多様な出力形式をサポートし、構造探索ツールとの統合により解析を可能に。

実験結果

リサーチクエスチョン

• RQ1 現代の宇宙論的および天体物理学的問題の複雑さを扱える、高スケーラブルでモジュラーかつオープンソースのシミュレーションコードをどのように設計できるか？
• RQ2 数百億個の粒子を含むシミュレーションにおいて、現代のCPUベースHPCシステムで最適なパフォーマンスを達成するための並列化戦略は何か？
• RQ3 重力、SPH、ニュートリノ動力学、銀河形成モデルといった複数の数値手法を、1つの拡張可能なフレームワーク内で効率的に統合する方法は何か？
• RQ4 このコードを用いた大規模な宇宙論的流体力学において、達成可能なパフォーマンスと正確性のベンチマークは何か？
• RQ5 計算天文学分野における再現可能性とコミュニティ主導の開発を支えるために、コードはどのように支援できるか？

主な発見

• SWIFT は、約3000億個の粒子を含む宇宙論的流体力学的問題において、大規模なコンピュータクラスタを横断して優れた弱スケーリング性能を示しており、スケーラビリティを実証。
• Sphenixスキームを含む複数のSPHバージョンが正常に統合されており、周期的体積におけるFMMと粒子メッシュ法を組み合わせた高度な重力ソルバーもサポート。
• 粒子ベースのアプローチによりニュートリノの進化が正確にモデル化され、大規模構造シミュレーションにおけるその宇宙論的効果の組み込みが可能に。
• モジュラー設計により、銀河形成の小スケールモデルがスムーズに統合可能で、コアコードとは別に2つの補完的ネットワークが公開。
• 包括的なドキュメンテーション、例題問題、解析ツールとともに、コードは公開されており、再現可能性とコミュニティ貢献を支援。
• 広範な一貫性および正確性のテストにより、さまざまな天体物理学的テスト問題において数値手法の信頼性が確認された。