[論文レビュー] A 3D radiative transfer framework: X. Arbitrary Velocity Fields in the Co-moving Frame

E. Baron, P. H. Hauschildt+2 more|Oct 24, 2012

Cosmology and Gravitation Theories参考文献 24被引用数 3

ひとこと要約

本論文は、光子の経路をアフィンパラメータ化することで、共動フレームにおける任意の速度場にまで拡張された3次元放射平衡移動フレームワークを提示する。準解析的解法とNg加速、ドメイン分割を組み合わせることで、複雑な波長結合に対しても高い計算効率と収束性を達成し、ジェットや衝撃波面などの天体的環境の正確なモデル化を可能にする。

ABSTRACT

3-D astrophysical atmospheres will have random velocity fields. We seek to combine the methods we have developed for solving the 1-D problem with arbitrary flows to those that we have developed for solving the fully 3-D relativistic radiative transfer problem in the case of monotonic flows. The methods developed in the case of 3-D atmospheres with monotonic flows, solving the fully relativistic problem along curves defined by an affine parameter, are very flexible and can be extended to the case of arbitrary velocity fields in 3-D. Simultaneously, the techniques we developed for treating the 1-D problem with arbitrary velocity fields are easily adapted to the 3-D problem. The algorithm we present allows the solution of 3-D radiative transfer problems that include arbitrary wavelength couplings. We use a quasi-analytic formal solution of the radiative transfer equation that significantly improves the overall computation speed. We show that the approximate lambda operator developed in previous work gives good convergence, even neglecting wavelength coupling. Ng acceleration also gives good results. We present tests that are of similar resolution to what has been presented using Monte-Carlo techniques, thus our methods will be applicable to problems outside of our test setup. Additional domain decomposition parallelization strategies will be explored in future work.

研究の動機と目的

• 共動フレームにおける任意の速度場を扱える一般化された3次元放射平衡移動フレームワークの開発。単調または球対称な流れにとどまらない。
• クェーサーのジェット、降着衝撃波、非単調な速度場を示す星風など、複雑な天体的系への相対論的放射平衡移動の適用範囲の拡張。
• 準解析的形式解法とNg加速を完全に3次元的・相対論的フレームワークに統合することで、計算効率の向上。
• チェックポイント技術を用いた再起動および解像度の適応を可能とする仕組みにより、高解像度のシミュレーションを実現。

提案手法

• 平坦な時空における光子経路のアフィンパラメータ化を用い、v/cのすべての次数において相対論的効果を正確に取り扱える。
• 光子の経路はアフィンパラメータξに沿って測地線として記述され、特定強度Iλは位置、時刻、方向、共動波長λに依存する。
• 移動方程式に準解析的形式解法を適用することで、標準的な数値積分法と比較して計算速度が著しく向上する。
• 全相対論的移動方程式に含まれる共動強度の時間発展項および源関数寄与項を含めた波長結合を組み込む。
• 大規模な計算リソースへのスケーリングを実現するため、ドメイン分割と並列化戦略を導入。
• 中間結果（Λ* と J）を保存するチェックポイント技術を用い、完全な再起動と反復回数の削減による解像度の向上を可能にする。

実験結果

リサーチクエスチョン

• RQ1 共動フレームにおける任意の非単調な速度場に対して、3次元的相対論的放射平衡移動をどのように正確に解けるか？
• RQ2 複雑な速度構造を有する3次元放射平衡移動において、準解析的形式解法が計算効率をどの程度向上させるか？
• RQ3 Ng加速とチェックポイント戦略は、任意の流れを有する高解像度3次元シミュレーションにおける収束時間をどの程度短縮できるか？
• RQ4 波長結合を無視した場合、近似ラムダ作用素は複雑な速度場においてどの程度の性能を示すか？
• RQ5 速度場の非対称性が、3次元大気における共動平均強度分布に及ぼす影響は何か？

主な発見

• アルゴリズムは、非球対称な流れ（例：ジェット構造）を含む、任意の速度場における3次元放射平衡移動問題を正常に解ける。
• 準解析的形式解法により計算時間が著しく短縮され、モンテカルロ法と同等の解像度のシミュレーションが可能になった。
• Ng加速は良好な収束率を示したが、再起動実行では初期履歴が欠落するため、効果の発現が遅れる傾向にあった。
• チェックポイント技術により完全な再起動が可能となり、低解像度結果を再利用することで、高解像度計算の反復回数を約半分に削減できた。
• 非対称な速度場では、極から赤道へと移行する間に共動平均強度が10倍近く変化し、連続体レベルに最大速度の影響が現れている。