[论文解读] Radiative hydrodynamics simulations of red supergiant stars. IV gray versus non-gray opacities

A. Chiavassa, B. Freytag+2 more|arXiv (Cornell University)|Sep 16, 2011

Stellar, planetary, and galactic studies参考文献 48被引用 63

一句话总结

本研究利用非灰度消光系数的CO5BOLD代码对红超巨星开展三维辐射流体动力学模拟，结果表明非灰度处理使温度梯度更陡峭，温度涨落减弱，相较于灰度模型表现更优。非灰度处理显著提升了光谱能量分布和干涉测量半径的精度，而通过经验校准发现，1D模型中的微湍流与巨湍流参数必须调整至与3D结果一致，才能准确匹配有效重力和谱线轮廓。

ABSTRACT

Red supergiants are massive evolved stars that contribute extensively to the chemical enrichment of our Galaxy. It has been shown that convection in those stars gives rise to large granules that cause surface inhomogeneities and shock waves in the photosphere. The understanding of their dynamics is crucial to unveil the unknown mass-loss mechanism, their chemical composition and stellar parameters. We present a new generation of red supergiants simulations with a more sophisticated opacity treatment done with 3D radiative- hydrodynamics CO5BOLD. In the code, the coupled equations of compressible hydrodynamics and non-local radiation transport are solved in the presence of a spherical potential. The stellar core is replaced by a special spherical inner boundary condition, where the gravitational potential is smoothed and the energy production by fusion is mimicked by a simply producing heat corresponding to the stellar luminosity. The post-processing radiative transfer code OPTIM3D is used to extract spectroscopic and interferometric observables. We show that the relaxation of the assumption of frequency-independent opacities shows a steeper mean thermal gradient in the optical thin region that affect strongly the atomic strengths and the spectral energy distribution. Moreover, the weaker temperature fluctuations reduce the incertitude on the radius determination with interferometry. We show that 1D models of red supergiants must include a turbulent velocity calibrated on 3D simulations to obtain the effective surface gravity that mimic the effect of turbulent pressure on the stellar atmosphere. We provide an empirical calibration of the ad-hoc micro- and macroturbulence parameters for 1D models using the 3D simulations: we find that there is not a clear distinction between the different macroturbulent profiles needed in 1D models to fit 3D synthetic lines.

研究动机与目标

• 研究非灰度消光系数对红超巨星（RSG）大气热结构与动力学的影响。
• 评估频率依赖性辐射转移相较于灰度近似对谱线形成与干涉观测量的影响。
• 校准1D模型中的微湍流与巨湍流参数，使其与3D模拟所得湍流速度场一致。
• 评估对流越激与边界条件对RSG表面非均质性与能量传输的影响。
• 确定非灰度模型是否能更准确再现红超巨星的观测属性，如有效温度、光度与半径。

提出的方法

• 使用具有球形引力势与模拟恒星光度的热源的CO5BOLD代码执行三维辐射流体动力学模拟。
• 采用多组分频率依赖性消光系数（5个波段）替代传统灰度近似以实现辐射转移。
• 利用OPTIM3D后处理代码从3D模拟输出中计算合成光谱与干涉可见度。
• 比较3D非灰度与灰度模拟在温度结构、速度场与表面非均质性方面的差异。
• 通过匹配3D模拟的合成谱线轮廓，校准1D模型的微湍流与巨湍流参数。
• 分析数值分辨率（255³与401³网格点）对对流结构与收敛性的影响。

实验结果

研究问题

• RQ1 非灰度消光系数如何影响红超巨星大气外层的平均热梯度与温度涨落？
• RQ2 相较于灰度模型，非灰度模拟在多大程度上提升了与观测光谱能量分布及干涉测量结果的一致性？
• RQ3 为使1D模型重现3D非灰度模拟所得的有效表面重力与谱线轮廓，所需的微湍流与巨湍流参数值是多少？
• RQ4 非灰度与灰度模拟中激波与对流运动有何差异，其对表面强度对比度有何影响？
• RQ5 提高3D模拟的数值分辨率是否能实现对流颗粒数量与尺寸的收敛？

主要发现

• 非灰度模拟在光学稀薄区域表现出比灰度模型更陡的平均热梯度，从而改变分子与谱线强度。
• 3D非灰度模型产生的光谱能量分布比3D灰度模拟更接近1D模型，表明非灰度处理对准确光谱合成至关重要。
• 由于增强了辐射传热，非灰度模型中温度涨落更弱，表面强度对比度降低，从而提升了干涉测量半径的精度。
• 灰度与非灰度3D模型中的微湍流速度相近，支持在1D模型中采用与深度无关的微湍流作为合理近似。
• 1D模型中的巨湍流剖面在拟合3D合成谱线时无明显区分，但非灰度模型表现出更高的速度场标准差，表明谱线轮廓更复杂。
• 非灰度模型因压力与密度标高更小，在外层表现出更强的激波，导致更复杂的速度场与谱线轮廓。