[論文レビュー] A 3D radiative transfer framework: II. line transfer problems

E. Baron, P. H. Hauschildt|ArXiv.org|Mar 16, 2007

Radiative Heat Transfer Studies参考文献 10被引用数 26

ひとこと要約

本稿では、波長および立体角方向の並列化と演算子分割を用いて、散乱を伴う線遷移問題を解くために長特徴線法を拡張した3次元放射移動フレームワークを提示する。この手法は、3次元空間分解能が低いにもかかわらず1次元ベンチマークと非常に良好な一致を示し、100,000プロセッサを超えて効率的なスケーリングを達成しており、現代の天体物理学的シミュレーションにおける高精度な3次元スペクトル合成を可能にする。

ABSTRACT

Higher resolution telescopes as well as 3D numerical simulations will require the development of detailed 3D radiative transfer calculations. Building upon our previous work we extend our method to include both continuum and line transfer. We present a general method to calculate radiative transfer including scattering in the continuum as well as in lines in 3D static atmospheres. The scattering problem for line transfer is solved via means of an operator splitting (OS) technique. The formal solution is based on a long-characteristics method. The approximate $Λ$ operator is constructed considering nearest neighbors {\em exactly}. The code is parallelized over both wavelength and solid angle using the MPI library. We present the results of several test cases with different values of the thermalization parameter and two choices for the temperature structure. The results are directly compared to 1D spherical tests. With our current grid setup the interior resolution is much lower in 3D than in 1D, nevertheless the 3D results agree very well with the well-tested 1D calculations. We show that with relatively simple parallelization that the code scales to very large number of processors which is mandatory for practical applications. Advances in modern computers will make realistic 3D radiative transfer calculations possible in the near future. Our current code scales to very large numbers of processors, but requires larger memory per processor at high spatial resolution.

研究の動機と目的

• 静的大気中における散乱を伴う線遷移を扱える3次元放射移動フレームワークの開発。
• 論文Iで提示された3次元連続体放射移動法を、スペクトル線および散乱支配の線遷移に拡張すること。
• 高分解能望遠鏡からの観測と比較可能な高分解能3次元スペクトル合成を可能にすること。
• 大規模並列アーキテクチャ上で効率的なスケーリングを確保し、実用的な天体物理学的応用を可能にすること。
• さまざまな熱化パラメータおよび温度構造の下で、3次元結果を検証し、テスト済みの1次元球対称計算と比較すること。

提案手法

• 放射移動方程式の形式的解法は、3次元空間および角度グリッド上で長特徴線法を用いる。
• 線遷移は演算子分割（OS）により解き、散乱問題は近似的なΛ演算子を最近接ノード寄与から構築することで扱う。
• 波長依存の平均強度JλおよびΛ*はプロファイル平均化され、源関数の反復的解法に使用される線統合量J̄およびΛ̄*を計算する。
• 源関数はS = (1−ε)J̄ + εBとして計算され、ここでεは線熱化パラメータ、Bは線プロファイルにわたるプランク関数の平均値である。
• 加速λ反復（ALI）法を用いてΛ*演算子を反復的に更新し、収束は反復式[1−Λ*(1−ε)]J̄_new = J̄_fs − Λ*(1−ε)J̄_oldにより監視する。
• コードは波長および立体角方向にMPIを用いて並列化され、形式的解法後に結果を集約するために集約操作が用いられる。

実験結果

リサーチクエスチョン

• RQ1 さまざまな熱化パラメータの下で、3次元長特徴線法は1次元球対称線遷移結果をどの程度正確に再現できるか？
• RQ2 3次元空間分解能の低さが、1次元と比較して線遷移計算の精度にどの程度影響を与えるか？
• RQ3 波長および立体角方向のMPI並列化を用いた場合、コードは大規模なプロセッサ数に対してどの程度効率的にスケーリングするか？
• RQ4 演算子分割およびΛ*近似戦略は、3次元散乱支配の線遷移問題においても収束性と正確性を維持できるか？
• RQ5 現実的な3次元グリッドサイズおよび高分解能スペクトルグリッドに適用した場合、フレームワークの計算コストとスケーラビリティはいかほどか？

主な発見

• 3次元グリッドの内部空間分解能が低いにもかかわらず、3次元結果はテスト済みの1次元球対称計算と非常に良好に一致する。
• コードは100,000プロセッサを超えて効率的にスケーリングされ、128CPUを用いた設定ではほぼ理想に近い負荷バランスが達成され、すべての並列構成で一貫した結果が得られた。
• Λ*演算子の構築は最初の反復でのみ必要であり、形式的解法とほぼ同等の時間がかかるが、以降の反復では事前に計算されたΛ*により高速化される。
• 波長クラスタと立体角のためのワーカープロセスを用いた並列化戦略により、通信オーバーヘッドを最小限に抑え、主要段階で集約MPI操作が効果的に使用された。
• 散乱支配の線（ε = 1）に対しても、温度構造や線熱化パラメータの違いに関わらず、1次元ベンチマークと一致する正確な結果が得られた。
• フレームワークは最大131,072プロセッサでの理論的実行が可能であり、今後の天体物理学的シミュレーションにおける大規模3次元スペクトル合成への強い可能性を示している。