[論文レビュー] A 3D radiative transfer framework: VI. PHOENIX/3D example applications

Peter H. Hauschildt, E. Baron|arXiv (Cornell University)|Nov 17, 2009

Radiative Heat Transfer Studies被引用数 5

ひとこと要約

本論文は、PHOENIX/3Dコード内での3次元放射移動フレームワークの実装と検証を提示しており、散乱を含む完全な3次元スペクトル合成を可能にするとともに、1次元モデルと同一の微視的物理学を維持している。現代のスーパーコンputersリソースを用いることで、階層的ドメイン分割とMPI並列化により強力なスケーリングとメモリ効率を達成し、M型矮星、太陽型星、A型星を含む多様な星の対象で一貫した結果を得られた。

ABSTRACT

We demonstrate the application of our 3D radiative transfer framework in the model atmosphere code PHOENIX/3D for a number of spectrum synthesis calculations for very different conditions. The 3DRT framework discussed in the previous papers of this series was added to our general-purpose model atmosphere code PHOENIX/1D and an extended 3D version PHOENIX/3D was created. The \phxT code is parallelized via the MPI library using a hierarchical domain decomposition and displays very good strong scaling. We present the results of several test cases for widely different atmosphere conditions and compare the 3D calculations with equivalent 1D models to assess the internal accuracy of the 3D modeling. In addition, we show the results for a number of parameterized 3D structures. With presently available computational resources it is possible to solve the full 3D radiative transfer (including scattering) problem with the same micro-physics as included in 1D modeling.

研究の動機と目的

• PHOENIXモデル大気コードを拡張し、散乱を含む完全な3次元放射移動をサポートする。1次元物理学と整合性を保つ。
• 多様な星のパラメータにわたる1次元モデルとの直接比較を通じて、3次元放射移動計算の正確性を検証する。
• 現在の計算リソースを用いて、散乱を含む完全な3次元放射移動問題を解く可能性を実証する。
• MPIと階層的ドメイン分割を用いた効率的かつスケーラブルな並列化戦略を実装する。
• フラックスベクトルとスペクトル合成結果の定量的比較を通じて、3次元フレームワークの内部正確性を評価する。

提案手法

• PHOENIX/3Dコードは、先行論文の3次元放射移動フレームワークを統合し、線および連続体における散乱を処理するために非局所的オペレータ分割を用いている。
• コードはMPIを用いた階層的ドメイン分割戦略を採用しており、プロセスのクラスタが波長と空間ボクセルのサブセットをそれぞれ処理する。
• 方程式状態と光学厚さデータをプロセス間で分散させることで、CPUとメモリの効率を向上させ、1024コアで1プロセスあたり約0.5GBのメモリ使用量にまで削減した。
• 3次元放射移動ソルバーは、角度分解能として$n_\theta = 64$および$n_\phi = 64$の立体角点を用い、超新星モデルでは球座標系でフラックスベクトルを計算している。
• スペクトル合成には、PHOENIX/1Dと同一のLTEラインプロファイル（ガウス型およびボイチ型）、光学厚さ（b-f、f-f、ダスト）、および状態方程式を用いており、物理的整合性を保証している。
• フレームワークは直角座標系および球座標系の両方をサポートしており、後者は超新星モデルで使用され、径方向および角度方向のフラックス構造を捉えている。

実験結果

リサーチクエスチョン

• RQ1 PHOENIX/3Dの3次元放射移動フレームワークは、同一の物理的条件下で1次元モデルと整合する合成スペクトルを生成できるか？
• RQ2 構造的星の大気において、3次元放射移動解は1次元解と比べてフラックスベクトルの非等方性およびエネルギーの流れのパターンでどのように異なるか？
• RQ3 散乱を含む3次元放射移動で達成できる数値的正確性の程度は、フラックス成分のずれによってどのように測定できるか？
• RQ4 階層的ドメイン分割を用いた現代のスーパーコンピュータ上で、3次元放射移動問題のスケーリングはどの程度効率的に行えるか？
• RQ5 完全な散乱と微視的物理学を含む3次元放射移動の計算要件とパフォーマンス特性は何か？

主な発見

• PHOENIX/3Dの3次元放射移動フレームワークは、M型矮星、太陽型星、A型星において、1次元の合成スペクトルを良好に再現しており、ほとんどの場合でスペクトル差が1%未満である。
• 超新星モデルでは、外側のボクセルにおいて$F_\theta/F_r$および$F_\phi/F_r$のフラックス比が0%から1%未満のずれに留まり、高い内部正確性を示している。
• 1024プロセスで強力なスケーリングを達成し、方程式状態のためのコアあたりメモリ使用量は数MBにまで削減され、完全な光学厚さおよび$\Lambda^*$-演算子の保存には約0.5GBとなった。
• 3次元構造における水平エネルギーの流れは、複雑で非等方的なフラックスパターンを示し、流れ線から非球形領域における顕著な横方向エネルギー輸送が確認された。
• 階層的ドメイン分割の使用により、メモリ使用量の効率化とパフォーマンス最適化が可能となり、大規模システムで最大約275,000ボクセルのシミュレーションが可能になった。
• 超新星の3次元モデルでは、フラックス成分$F_\theta$および$F_\phi$がゼロでないことが確認され、1次元モデルには存在しない3次元構造における横方向エネルギー輸送が存在することが裏付けられた。