[论文解读] Local models of stellar convection: Reynolds stresses and turbulent heat transport

P. J. Käpylä, M. J. Korpi+1 more|CERN Bulletin|Dec 15, 2003

Solar and Space Plasma Dynamics参考文献 33被引用 46

一句话总结

本研究利用三维局部MHD模拟，探究恒星对流层中雷诺应力与湍流热输运的特性，重点关注自转与磁场如何影响角动量与热输运。结果表明，自转通过负的水平雷诺应力（Qθφ）诱导赤道方向的角动量输运，且输运方向（向内或向外）取决于自转速率；磁场虽增强应力但未导致抑制效应，径向热通量在赤道处达到峰值，中纬度处出现最小值。

ABSTRACT

We study stellar convection using a local three-dimensional MHD model, with which we investigate the influence of rotation and large-scale magnetic fields on the turbulent momentum and heat transport. The former is studied by computing the Reynolds stresses, the latter by calculating the correlation of velocity and temperature fluctuations, both as functions of rotation and latitude. We find that the horisontal correlation, Q_(theta phi), capable of generating horisontal differential rotation, is mostly negative in the southern hemisphere for Coriolis numbers exceeding unity, corresponding to equatorward flux of angular momentum in accordance with solar observations. The radial component Q_(r phi) is negative for slow and intermediate rotation indicating inward transport of angular momentum, while for rapid rotation, the transport occurs outwards. Parametrisation in terms of the mean-field Lambda-effect shows qualitative agreement with the turbulence model of Kichatinov & Rüdiger (1993) for the horisontal part H \propto Q_(theta phi)/cos(theta), whereas for the vertical part, V \propto Q_(r phi)/sin(theta), agreement only for intermediate rotation exists. The Lambda-coefficients become suppressed in the limit of rapid rotation, this rotational quenching being stronger for the V component than for H. We find that the stresses are enhanced by the presence of the magnetic field for field strengths up to and above the equipartition value, without significant quenching. Concerning the turbulent heat transport, our calculations show that the transport in the radial direction is most efficient at the equatorial regions, obtains a minimum at midlatitudes, and shows a slight increase towards the poles. The latitudinal heat transport does not show a systematic trend as function of latitude or rotation.

研究动机与目标

• 理解自转与磁场如何通过湍流动量与热输运在恒星对流层中生成大尺度流。
• 量化雷诺应力与湍流热通量相关性作为自转速率与纬度的函数，基于局部三维MHD模型。
• 检验平均场参数化（Λ-效应）对直接模拟结果的有效性，尤其在不同自转与磁场强度下。
• 评估湍流输运的自转抑制效应，以及大尺度方位磁场对应力生成与热输运的影响。
• 通过量化局部湍流输运特性，为未来恒星自转与热结构的全球平均场建模提供数据支持。

提出的方法

• 采用具有周期性边界条件的三维局部非等熵MHD模型，模拟旋转稳定分层层中的对流。
• 通过速度涨落计算雷诺应力张量 Qij = ⟨u′i u′j⟩，以量化湍流动量输运。
• 计算湍流热通量相关性 ⟨u′i T′⟩，以评估纬向与径向热输运效率。
• 施加平均方位磁场（最高达能量平衡强度），研究其对湍流应力与麦克斯韦应力的影响。
• 采用Λ-效应形式化方法，将雷诺应力与平均场输运关联：Qij = Λijk Ωk + 高阶项。
• 对不同科里奥利数与磁场强度下的应力分量（Qθφ, Qrφ, Qrθ）进行体积平均分析。

实验结果

研究问题

• RQ1 自转如何影响水平雷诺应力Qθφ的符号与大小？是否导致赤道方向的角动量输运？
• RQ2 角动量输运的径向依赖性（Qrφ）如何随自转速率（慢速与快速自转）变化？
• RQ3 大尺度方位磁场如何影响湍流应力的生成与大小？是否存在抑制或增强的证据？
• RQ4 湍流热通量如何随纬度与半径变化？是否表现出系统性模式，如赤道最大值或中纬度最小值？
• RQ5 模拟得到的Λ-效应（H与V分量）在多大程度上与KR93等分析模型一致，尤其是在自转抑制条件下？

主要发现

• 当科里奥利数 >1 时，南半球的水平雷诺应力Qθφ为负，表明存在赤道方向的角动量输运，与太阳观测结果一致。
• 在慢速与中等自转条件下，径向应力Qrφ为负，意味着角动量向内输运；在快速自转条件下，Qrφ变为正，表明角动量向外输运。
• 垂直Λ-效应（V ∝ Qrφ / sinθ）仅在中等自转条件下与KR93模型有良好一致性；在快速自转区域出现明显偏差。
• 水平Λ-效应（H ∝ Qθφ / cosθ）在极区附近与KR93模型定性一致，但模拟值更大，且在赤道处抑制效应更弱。
• 自转抑制效应同时抑制两种Λ-效应，其中垂直分量V的抑制更强，且在更低自转速率下即发生，而水平分量H的抑制较弱。
• 方位磁场的存在使总湍流应力单调增强，直至并超过能量平衡强度，未出现显著抑制；在强磁场下，水平Λ-效应符号反转为负，表明角动量向两极输运。