[論文レビュー] The effects of numerical resolution on hydrodynamical surface convection simulations and spectral line formation

M. Asplund, H.‐G. Ludwig+2 more|arXiv (Cornell University)|May 15, 2000

Solar and Space Plasma Dynamics参考文献 2被引用数 30

ひとこと要約

本研究では、太陽表面対流の流体力学的シミュレーションにおい、数値的分解能と次元性（2次元対3次元）の影響、およびそれらが得られるスペクトル線プロファイルに与える影響を調査する。異なる分解能の3次元および2次元シミュレーションを用いて、約100³の分解能における3次元シミュレーションが観測された線の形状や非対称性に収束することを明らかにした。一方、2次元シミュレーションは系統的に広く浅い線を生成し、誤ったドップラー速度シフトを引き起こすため、正確な元素含有量分析やスペクトル診断には3次元シミュレーションが不可欠である。

ABSTRACT

The computationally demanding nature of radiative-hydrodynamical simulations of stellar surface convection warrants an investigation of the sensitivity of the convective structure and spectral synthesis to the numerical resolution and dimension of the simulations, which is presented here. With too coarse a resolution the predicted spectral lines tend to be too narrow, reflecting insufficient Doppler broadening from the convective motions, while at the currently highest affordable resolution the line shapes have converged essentially perfectly to the observed profiles. Similar conclusions are drawn from the line asymmetries and shifts. In terms of abundances, weak FeI and FeII lines show a very small dependence (~0.02 dex) while for intermediate strong lines with significant non-thermal broadening the sensitivity increases (~0.10 dex). Problems arise when using 2D convection simulations to describe an inherent 3D phenomenon, which translates to inaccurate atmospheric velocity fields and temperature and pressure structures. In 2D the theoretical line profiles tend to be too shallow and broad compared with the 3D calculations and observations, in particular for intermediate strong lines. In terms of abundances, the 2D results are systematically about 0.1 dex lower than for the 3D case for FeI lines. Furthermore, the predicted line asymmetries and shifts are much inferior in 2D. Given these shortcomings and computing time considerations it is better to use 3D simulations of even modest resolution than high-resolution 2D simulations.

研究の動機と目的

• 3次元流体力学的対流シミュレーションにおけるスペクトル線プロファイルおよび非対称性予測の正確性を確保するための最小数値的分解能を特定すること。
• 特にスペクトル線形成に関して、3次元星の対流現象をモデル化するにあたり、2次元シミュレーションの信頼性を評価すること。
• モデル大気における分解能および次元性の変化が、元素含有量決定に与える感受性を評価すること。
• 観測された太陽スペクトルと比較して、3次元と2次元シミュレーションの間で線プロファイル、非対称性、およびシフトの収束性を比較すること。
• 計算制約を考慮した場合、低分解能の3次元シミュレーションが高分解能の2次元シミュレーションに代わってスペクトル解析に使用可能かどうかを検討すること。

提案手法

• 異なるグリッド分解能（50³から200³まで）を用いた3次元コードを用いた、太陽表面対流の時間発展的かつ圧縮性・放射線を含む流体力学的シミュレーションと、それに対応する2次元シミュレーションを実施する。
• 分解能と次元性の影響を明確に分離するために、すべてのシミュレーションで同一の数値的パラメータを用い、差分比較を可能にする。
• 統計的有意性を確保するため、1つの線あたり10⁵個以上の1次元的計算を要する、空間的・時間的平均スペクトル線プロファイルの計算を、局所熱平衡（LTE）スペクトル合成法を用いて行う。
• 形状、深さ、非対称性、およびシフトの観点から線プロファイルを分析し、観測された太陽スペクトルと比較する。
• 合成されたプロファイルを用いてFe I、Fe II、およびHの線から元素含有量を算出し、分解能および次元性に依存する度合いを定量化する。
• スペクトル線形成に寄与するモデル大気内の温度構造および速度構造を評価し、その役割を理解する。

実験結果

リサーチクエスチョン

• RQ1 3次元シミュレーションにおける数値的分解能が、観測結果と比較して予測されるスペクトル線プロファイルおよび非対称性の正確性にどのように影響するか？
• RQ2 2次元シミュレーションは、3次元シミュレーションと比較して、太陽スペクトルで観測された線プロファイル、シフト、非対称性をどの程度正しく再現できないか？
• RQ3 線プロファイルおよび非対称性予測の収束を達成するための3次元シミュレーションにおける最小分解能は何か？
• RQ4 Fe I、Fe II、およびHの線から得られる元素含有量は、対流シミュレーションの分解能および次元性にどのように依存するか？
• RQ5 低分解能でも、減衰翼を有する強い線やH線を用いてT_effおよびlog gの補正に信頼性を持って使用可能か？その理由は何か？

主な発見

• 約100³グリッドポイントの分解能において、3次元シミュレーションは観測値にほぼ完全に収束する線プロファイルおよび非対称性を生成しており、現在の観測不確実性を考慮しても十分な精度であることが示された。
• 2次元シミュレーションは、特に中程度強度の線に対して、系統的に浅く広がった線プロファイルを生成しており、これは誤った対流速度構造に起因する。
• 2次元シミュレーションにおける線シフトおよび非対称性は、観測結果と最大で約200 m s⁻¹の差を示しており、対流ダイナミクスの不適切な表現を示している。
• 2次元シミュレーションにおけるFe I線からの元素含有量は、3次元シミュレーションと比較して系統的に約0.1 dex低くなるため、顕著な系統的誤差が生じている。
• 弱いFe IおよびFe II線は分解能依存性がわずか（約0.02 dex）である一方、中程度強度の線はより大きな感受性（約<0.10 dex）を示す。
• 強いFe IおよびH線は、分解能にほとんど依存せず、温度および圧力構造が低分解能で収束するため、低分解能でもT_effおよびlog gの補正に信頼性を持って使用可能である。