[论文解读] Transport and mixing in the radiation zones of rotating stars: I-Hydrodynamical processes

S. Mathis, Jean-Paul Zahn|ArXiv.org|Jun 18, 2004

Stellar, planetary, and galactic studies参考文献 20被引用 104

一句话总结

该论文通过将流体动力学处理扩展至更高阶球谐函数并保留全部时间导数，改进了恒星辐射区中旋转混合的建模，实现了对较差对流层环流和快速演化阶段的精确模拟。通过捕捉纬向速度差和非定常动力学，该模型解决了先前模型的局限性，尤其适用于具有强剪切驱动湍流和经向环流的大质量恒星。

ABSTRACT

The purpose of this paper is to improve the modelization of the rotational mixing which occurs in stellar radiation zones, through the combined action of the thermally driven meridional circulation and of the turbulence generated by the shear of differential rotation. The turbulence is assumed to be anisotropic, due to the stratification, with stronger transport in the horizontal directions than in the vertical. The main difference with the former treatments by Zahn (1992) and Maeder & Zahn (1998) is that we expand here the departures from spherical symmetry to higher order, and include explicitly the differential rotation in latitude, to first order. This allows us to treat simultaneously the bulk of a radiation zone and its tachocline(s). Moreover, we take fully into account the non-stationarity of the problem, which will enable us to tackle the rapid phases of evolution. The system of partial differential equations, which govern the transport of angular momentum, heat and chemical elements, is written in a form which makes it ready to implement in a stellar evolution code. Here the effect of a magnetic field is deliberately ignored; it will be included in forthcoming papers.

研究动机与目标

• 通过将偏离球对称性的处理扩展至标准1.5D近似之外，改进恒星辐射区中旋转混合的建模。
• 在线性近似下显式包含纬向速度差的影响，从而实现对较差对流层环流的处理。
• 保留所有时间导数（动力松弛除外），以更好地解析旋转和成分中出现陡峭梯度的快速演化阶段。
• 提供一套适合在恒星演化代码中实现的偏微分方程系统，重点聚焦于角动量、热量和化学元素的输运。
• 为未来在旋转混合过程中引入磁场（目前被忽略）奠定基础。

提出的方法

• 将所有流体动力学变量展开为更高阶球谐函数（超越标准1.5D方法），从而显式处理纬向变化。
• 采用松弛法处理温度涨落（Ψₗ），在成分和旋转梯度陡峭区域稳定数值解。
• 由于稳定层结，将湍流输运建模为各向异性（νᵥ ≪ νₕ, Dᵥ ≪ Dₕ），偏好水平方向而非垂直方向的混合。
• 包含由热不平衡驱动的经向环流（源自von Zeipel的工作）及其通过角动量输运对旋转的反馈作用。
• 施加边界条件，确保在恒星表面（r = R）与势场连续，包括Ψₗ的径向导数为零，以及Λₗ的混合条件。
• 采用线性化、轴对称的公式描述通过Ψₗ耦合的差速旋转与通过Λₗ耦合的成分不均匀性之间的相互作用，避免引入使数值处理复杂化的非线性项。

实验结果

研究问题

• RQ1 如何改进辐射区中经向环流与剪切驱动湍流的处理，以更准确地捕捉较差对流层动力学？
• RQ2 包含更高阶球谐项对纬向速度差和环流建模有何影响？
• RQ3 保留全部时间导数如何提升对旋转和成分中出现强梯度的快速演化阶段的模拟能力？
• RQ4 温度涨落Ψₗ的松弛公式在恒星内部陡峭梯度附近如何提升数值稳定性？
• RQ5 如何构建方程系统，使其可直接在恒星演化代码中实现，同时保持物理解释的一致性？

主要发现

• 该模型成功捕捉了由对流区与辐射区之间差速旋转引发的较差对流层中的八极环流，该现象此前被人为处理。
• 通过保留全部时间导数（动力松弛除外），该模型能够精确模拟角速度和成分梯度显著增大的快速演化阶段。
• 温度涨落Ψₗ的松弛公式降低了对陡峭成分梯度的数值敏感性，提升了恒星演化代码中的稳定性。
• 引入更高阶球谐函数使得旋转和环流的纬向变化得以更准确地表示，突破了1.5D壳层近似。
• 该模型提供了一个自洽的框架，用于角动量、热量和化学元素输运，适合集成到恒星演化代码中。
• 该框架明确设计为未来可纳入磁场，而磁场被识别为解释类太阳恒星中平坦旋转剖面的关键因素，与当前模型预测的快速旋转核心相反。