[論文レビュー] Achieving Extreme Resolution in Numerical Cosmology Using Adaptive Mesh Refinement: Resolving Primordial Star Formation

Greg L. Bryan, Tom Abel+1 more|arXiv (Cornell University)|Dec 4, 2001

Computer Graphics and Visualization Techniques参考文献 12被引用数 6

ひとこと要約

本論文は、128ビット拡張精度算術と34段階のリファインメントレベルを持つ8,000個のサブグリッドを有する階層的グリッドを用いて、標準的宇宙論的物理法則に基づく初代星形成を、空間的および時間的動的範囲が未曾有の高さに達する3次元適応メッシュリファインメント（AMR）宇宙論的シミュレーションを提示する。Enzoコードを用いることで、初代星の形成を第一原理的にモデル化することが可能になった。

ABSTRACT

As an entry for the 2001 Gordon Bell Award in the "special" category, we describe our 3-d, hybrid, adaptive mesh refinement (AMR) code, Enzo, designed for high-resolution, multiphysics, cosmological structure formation simulations. Our parallel implementation places no limit on the depth or complexity of the adaptive grid hierarchy, allowing us to achieve unprecedented spatial and temporal dynamic range. We report on a simulation of primordial star formation which develops over 8000 subgrids at 34 levels of refinement to achieve a local refinement of a factor of 10^12 in space and time. This allows us to resolve the properties of the first stars which form in the universe assuming standard physics and a standard cosmological model. Achieving extreme resolution requires the use of 128-bit extended precision arithmetic (EPA) to accurately specify the subgrid positions. We describe our EPA AMR implementation on the IBM SP2 Blue Horizon system at the San Diego Supercomputer Center.

研究の動機と目的

• 適応メッシュリファインメントを用いて、宇宙における最初の星の形成を極めて高分解能でシミュレートすること。
• 静的グリッドや従来のN体コードの限界を克服し、宇宙論的スケールから星間スケールまで局所的な解像度を確保すること。
• 空間的および時間的両方の極めて広い動的範囲を実現し、初期星形成ガスの完全な進化を解明すること。
• 複雑で結合された物理現象を扱える高性能で並列処理可能なAMRコードの実装とベンチマーク化。
• 標準的物理法則と宇宙論的モデルに基づき、初期状態からの正確な第一原理的星形成シミュレーションを可能にすること。

提案手法

• 物理的状態に応じて計算グリッドを動的にリファインする3次元ハイブリッド適応メッシュリファインメント（AMR）コード「Enzo」を採用。
• 最大34段階のリファインメントを持つ構造的グリッド階層を用い、8,000以上のサブグリッドを生成することで、星間密度にまで到達する局所スケールを解像。
• 複数のソルバーを統合：Euler形式の流体力学、ダークマターのN体重力、Poisson重力ソルバー、初期星形成ガスの12種類の反応機構を有する剛性化学ソルバー。
• 極めて高い解像度下でも数値的精度を維持できるよう、128ビット拡張精度算術（EPA）を用いて小スケール位置を正確に表現。
• IBM SP2 Blue Horizonシステムの64プロセッサを用いて並列処理を実行し、最適化された通信と負荷分散により、約13 Gflop/sの継続的性能を達成。
• GUIベースの可視化ツール「Jacques」を用いて、階層的AMRデータをインタラクティブに探索・レンダリング可能にし、細分化されたスケールでの分析が可能（最大$10^{10}$倍のズーム機能を備える）。

実験結果

リサーチクエスチョン

• RQ1 適応メッシュリファインメントは、3次元宇宙論的シミュレーションにおいて空間的および時間的動的範囲$10^{12}$を達成できるか？
• RQ2 128ビット拡張精度算術は、極めて高解像度のシミュレーションにおいて、数値的精度を維持するために効果的に利用可能か？
• RQ3 流体力学、重力、化学反応の完全な結合を伴う初期星形成のシミュレーションにおける計算コストと性能はいかほどか？
• RQ4 34段階のリファインメントレベルにわたる、結合された偏微分方程式および常微分方程式の複雑さを処理できる並列AMRコードは実現可能か？
• RQ5 メモリ制限があるシステムにおいて、多段階かつ高動的範囲を持つデータのインタラクティブ可視化はどのように実現できるか？

主な発見

• 本シミュレーションは、3次元天体物理学的シミュレーションで報告された中で最高水準の空間的および時間的動的範囲$10^{12}$を達成した。
• コードは、約100太陽半径の原始星の形成を正確に解像し、初期宇宙スケールに比べて$10^{12}$倍の解像度向上を実現した。
• IBM SP2 Blue Horizonシステムの64プロセッサ上で、シミュレーションの継続的性能は約13 Gflop/sであった。
• 流体力学ソルバーが合計計算時間の36％を占め、次いでPoisson重力ソルバー（17％）、境界条件処理（15％）が続く。
• 静的グリッドの同等シミュレーションでは、仮想的なフロップレートが$10^{44}$ flop/sに達するため、AMRが均一グリッドに比べて極めて効率的であることが示された。
• 可視化ツール「Jacques」により、AMR階層のインタラクティブな探索が可能となり、$10^{10}$倍のズーム機能を活用して希少で高解像度の特徴を特定・分析できるようになった。