[論文レビュー] Radiative transfer in moving media II. Solution of the radiative transfer equation in axial symmetry

D. Korčáková, J. Kubát|arXiv (Cornell University)|Jun 16, 2005

Astrophysics and Star Formation Studies参考文献 21被引用数 10

ひとこと要約

本稿では、軸対称な2次元放射移動方程式に対して、長距離特徴線と短距離特徴線を組み合わせ、局所ローレンツ変換を用いることで、星風や回転する大気のような任意の速度場を正確にモデル化する、新しい形式的解法を提示する。この手法によりドップラーシフトや周波数依存の縁減光を正確に扱えるようになり、ソボレフ近似に依存せずに、特に速度勾配が小さい領域でも正確な線分光プロファイルと放射束を得ることができる。

ABSTRACT

A new method for the formal solution of the 2D radiative transfer equation in axial symmetry in the presence of arbitrary velocity fields is presented. The combination of long and short characteristics methods is used to solve the radiative transfer equation. We include the velocity field in detail using the Local Lorentz Transformation. This allows us to obtain a significantly better description of the photospheric region, where the gradient of the global velocity is too small for the Sobolev approximation to be valid. Sample test calculations for the case of a stellar wind and a rotating atmosphere are presented.

研究の動機と目的

• 高温星の風と回転を伴う状況を想定し、任意の速度場を有する2次元軸対称における放射移動方程式を解くための堅牢な数値手法の開発。
• 速度勾配が小さい領域では崩壊するソボレフ近似の限界を克服すること。
• 非径方向の光線におけるドップラーシフトと周波数依存の縁減光を、特に回転する星の大気に正確にモデル化すること。
• 一貫した放射移動を用いて、星の光球層、星風、および降着円盤のような拡張構造を同時にモデル化できること。
• 多次元で動く媒体における合成分光合成のための、計算的に効率的かつ正確な解法を提供すること。

提案手法

• 長距離特徴線と短距離特徴線を、軸対称対象と交差する平面内で組み合わせたハイブリッド手法を用い、精度と計算コストのバランスをとる。
• 速度場は局所ローレンツ変換により組み込み、フレームの移動に伴う周波数と強度の正しい変換を実現し、ソボレフ近似の必要性を回避する。
• ドップラーシフトを施した周波数グリッドを用いて、共動フレームで移動媒体内における分光線形成を正確に解く。
• セル境界での強度補間はドップラーシフトを考慮し、速度勾配にわたって分光線構造を保持する。
• 適応的解像度を有する構造化グリッド上に実装されており、大きな光学的厚さの対比に対しても効率的に処理できる。
• 計算時間は角度および幾何グリッドポイント数に線形に比例し、複雑な星系の実用的シミュレーションが可能になる。

実験結果

リサーチクエスチョン

• RQ1 任意の速度勾配（高・低を含む）を持つ2次元軸対称における放射移動方程式を、形式的にどのように解くことができるか。
• RQ2 周波数依存の縁減光は、回転する星の大気における回転幅広がり線分光プロファイルにどのような影響を与えるか。
• RQ3 ソボレフ近似に依存せずに、星風の光球層から風領域への遷移を正確にモデル化できるか。
• RQ4 非径方向の光線に沿ったドップラーシフトの取り入れ方は、回転する星の顕在線分光プロファイルにどのように影響するか。
• RQ5 角度、幾何、周波数グリッド解像度の観点から、この手法の計算的効率性とスケーラビリティはいかほどか。

主な発見

• 静的かつ球対称モデルにおいて、正しい縁減光則を正確に再現でき、速度場が存在しない場合の手法の正確性を検証した。
• β速度則に従う星風に対して、コードはP Cygni型特徴に類似した線分光プロファイルの青方シフトと変形を再現したが、完全な動的風解を用いていなくても成立した。
• 周波数依存の縁減光を含めることは、正確な回転幅広がり線分光プロファイルのための不可欠な要素であり、これを無視すると顕著な放射束誤差が生じる。
• 薄い大気では標準的な畳み込み法と完全な移動方程式解法との間で放射束の差は小さいが、物理的正確性を確保するためには完全な手法が不可欠である。
• 計算時間は角度および幾何グリッドポイント数に線形に比例しており、実用的シミュレーションにおける良好なスケーラビリティを示している。
• この手法は星風および回転する星のモデル化に有効であり、干渉観測や降着円盤系への応用の可能性を有する。