Shear mixing in stellar radiative zones I. Effect of thermal diffusion and chemical stratification
Shear mixing in stellar radiative zones I. Effect of thermal diffusion and chemical stratification
复制标题
恒星辐射区的剪切混合 I. 热扩散和化学分层的影响
DOI:
10.1051/0004-6361/201423655
复制
发表时间:
2014
影响因子:
6.5
通讯作者:
F. Lignieres
中科院分区:
文献类型:
--
作者:
V. Prat;F. Lignieres
Context. Turbulent transport of chemical elements in radiative zones of stars is considered in current stellar evolution codes thanks to phenomenologically derived di usion coe cients. Recent local numerical simulations (Prat & Lignieres 2013, A&A, 551, L3) suggest that the coe cient for radial turbulent di usion due to radial di erential rotation satisfies Dt ’ 0:058= Ri, in qualitative agreement with the model of Zahn (1992, A&A, 265, 115). However, this model does not apply (i) when di erential rotation is strong with respect to stable thermal stratification or (ii) when chemical stratification has a significant dynamical e ect, a situation encountered at the outer boundary of nuclear-burning convective cores. Aims. We extend our numerical study to consider the e ects of chemical stratification and of strong shear, and compare the results with prescriptions used in stellar evolution codes. Methods. We performed local, direct numerical simulations of stably stratified, homogeneous, sheared turbulence in the Boussinesq approximation. The regime of high thermal di usivities, typical of stellar radiative zones, is reached thanks to the so-called smallPeclet-number approximation, which is an asymptotic development of the Boussinesq equations in this regime. The dependence of the di usion coe cient on chemical stratification was explored in this approximation. Results. Maeder’s extension of Zahn’s model in the strong-shear regime (Maeder 1995, A&A, 299, 84) is not supported by our results, which are better described by a model found in the geophysical literature. As regards the e ect of chemical stratification, our quantitative estimate of the di usion coe cient as a function of the mean gradient of mean molecular weight leads to the formula Dt ’ 0:45 (0:12 Ri )=Ri, which is compatible in the weak-shear regime with the model of Maeder & Meynet (1996, A&A, 313, 140) but not with Maeder’s (1997, A&A, 321, 134).