Issue 
A&A
Volume 649, May 2021



Article Number  A62  
Number of page(s)  7  
Section  Stellar structure and evolution  
DOI  https://doi.org/10.1051/00046361/202039281  
Published online  11 May 2021 
Anisotropic turbulent transport with horizontal shear in stellar radiative zones
AIM, CEA, CNRS, Université ParisSaclay, Université Paris Diderot, Sorbonne Paris Cité, 91191 GifsurYvette, France
email: vincent.prat@cea.fr
Received:
28
August
2020
Accepted:
3
March
2021
Context. Turbulent transport in stellar radiative zones is a key ingredient of stellar evolution theory, but the anisotropy of the transport due to the stable stratification and the rotation of these regions is poorly understood. The assumption of shellular rotation, which is a cornerstone of the socalled rotational mixing, relies on an efficient horizontal transport. However, this transport is included in many stellar evolution codes through phenomenological models that have never been tested.
Aims. We investigate the impact of horizontal shear on the anisotropy of turbulent transport.
Methods. We used a relaxation approximation (also known as τ approximation) to describe the anisotropising effect of stratification, rotation, and shear on a background turbulent flow by computing velocity correlations.
Results. We obtain new theoretical scalings for velocity correlations that include the effect of horizontal shear. These scalings show an enhancement of turbulent motions, which would lead to a more efficient transport of chemicals and angular momentum, in better agreement with helio and asteroseismic observations of rotation in the whole HertzsprungRussell diagram. Moreover, we propose a new choice for the nonlinear time used in the relaxation approximation, which characterises the source of the turbulence.
Conclusions. For the first time, we describe the effect of stratification, rotation, and vertical and horizontal shear on the anisotropy of turbulent transport in stellar radiative zones. The new prescriptions need to be implemented in stellar evolution calculations. To do so, it may be necessary to implement nondiffusive transport.
Key words: hydrodynamics / turbulence / methods: analytical / stars: evolution / stars: interiors / stars: rotation
© V. Prat and S. Mathis 2021
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1. Introduction
The transport of angular momentum and chemical species plays a key role in stellar evolution. In particular, rotation induces transport both through largescale motions such as meridional circulation and through turbulent, smallscale motions generated by hydrodynamical instabilities (Zahn et al. 1983; Zahn 1992; Maeder & Zahn 1998; Mathis & Zahn 2004; Maeder 2009). These multidimensional smallscale transport processes cannot be resolved in onedimensional (1D) stellar evolution codes, so their secular effects have to be modelled in those codes.
In the last decades, helio and asteroseismology have provided us with many constraints on the internal rotation of the Sun and distant stars. In the case of the Sun, helioseismic data bring to light an almost flat rotation profile in the solar radiative zone down to 0.2 R_{⊙}, and a possibly faster core (Brown et al. 1989; Thompson et al. 2003; García et al. 2007; Fossat et al. 2017). Benomar et al. (2015) detected only weak differential rotation on a few solartype stars. For a large number of subgiant and red giant stars, asteroseismic data show a weak coresurface rotation contrast (Beck et al. 2012; Mosser et al. 2012; Deheuvels et al. 2012, 2014, 2015; Triana et al. 2017; Gehan et al. 2018). Weak differential rotation has also been detected in intermediatemass and massive stars (Kurtz et al. 2014; Saio et al. 2015; Triana et al. 2015; Murphy et al. 2016; Van Reeth et al. 2016; Aerts et al. 2017; Van Reeth et al. 2018; Ouazzani et al. 2019). Finally, a strong transport of angular momentum is needed to explain data obtained for white dwarves (Suijs et al. 2008; Hermes et al. 2017) and neutron stars (Heger et al. 2005; Hirschi & Maeder 2010).
Zahn (1992) introduced a formalism for turbulent transport accounting for both advection by the meridional circulation and the anisotropic turbulent transport induced by differential rotation, which is modelled thanks to vertical and horizontal turbulent diffusion coefficients. The anisotropy of the transport is due to the combined effect of rotation and stable stratification, which limit horizontal and vertical motions through the buoyancy force and the Coriolis acceleration, respectively. Those coefficients are based on phenomenological arguments and thus introduce large uncertainties in stellar evolution models. A lot of work has been done to improve the physics of the vertical diffusion coefficient by adding physical ingredients such as gradients of mean molecular weight (Maeder & Meynet 1996; Talon & Zahn 1997), and local numerical simulations have been performed to test existing models (Prat & Lignières 2013, 2014; Garaud & Kulenthirarajah 2016) and investigate the effect of viscosity (Prat et al. 2016; Garaud et al. 2017; Gagnier & Garaud 2018). Current models predict less transport than what is needed to explain observations of internal rotation (Eggenberger et al. 2012; Ceillier et al. 2013; Marques et al. 2013; Cantiello et al. 2014; Ouazzani et al. 2019).
Concerning the horizontal diffusion coefficient, since the inital model of Zahn (1992) that was based on phenomenological arguments, other models based on energetic arguments (Maeder 2003) or results of an unstratified TaylorCouette experiment (Mathis et al. 2004) have been proposed, but this coefficient is still poorly constrained, and its impact on stellar evolution is clearly overlooked. The first numerical simulations of the turbulent transport triggered by the instabilities of a horizontal shear in stellar conditions have also been performed (Cope et al. 2020; Garaud 2020). Finally, Park et al. (2020, 2021) studied the linear instabilities of the horizontal differential rotation in stellar radiation zones as a function of stratification, rotation, and thermal diffusivity.
However, the widely used formalism of Zahn (1992) assumes that the turbulent transport of angular momentum is viscous, which is not necessarily the case. Indeed, anisotropic turbulent motions can generate nonzero turbulent fluxes of angular momentum even in the absence of differential rotation, which is a characteristic of nondiffusive transport. The existence of such nondiffusive turbulent transport associated with the Coriolis force, known as the Λ effect (Rüdiger 1989), has been confirmed by direct numerical simulations (see e.g. Käpylä 2019). Originally, this effect was introduced in the context of solar and stellar convective zones to explain the persistance of differential rotation despite the presence of a magnetic field, which was supposed to erase it. Kitchatinov & Brandenburg (2012, hereafter referred to as KB12) later investigated the effect of the stable stratification of stellar radiative zones on the Reynolds stress.
Based on this work, Mathis et al. (2018, hereafter referred to as M+18) proposed a new model of the horizontal transport induced by the vertical shear instability. In particular, M+18 added a vertical shear (due to radial differential rotation) to the formalism of KB12 and show that it has no significant effect on the anisotropy of turbulence. Besides, they used estimates of mean velocity correlations to propose new prescriptions for the horizontal turbulent diffusion coefficient of chemical elements. The implementation of these prescriptions in stellar evolution computations shows a slightly enhanced transport of angular momentum throughout the main sequence, but not enough to fit helio and asteroseismic observations. Furthermore, a major caveat of this implementation is that it assumed that the transport of angular momentum had the same (diffusive) nature as the transport of chemical elements, and thus completely ignored the Λ effect.
In the present work, we perform a new generalisation of the formalism of KB12 in the presence of a general (both vertical and horizontal) shear. It allows us to propose new prescriptions for the transport of both angular momentum and chemical elements that include the Λ effect.
In Sect. 2, we derive spectral properties of the flow anisotropised by stratification, rotation, and shear. Then, in Sect. 3, we deduce scalings for velocity correlations as a function of stratification, rotation, and shear. We propose possible interpretations of the new scalings in Sect. 4. We explain how nonviscous turbulent transport of angular momentum should be implemented in stellar evolution codes to account for the Λ effect in Sect. 5. Finally, we conclude on the impact of these results on the modelling of turbulence in stellar evolution codes in Sect. 6.
2. Turbulent spectrum
In the Boussinesq approximation, which is justified here because turbulent motions are at length scales much smaller than the pressure scale height and velocity scales much smaller than the sound speed, the equations governing the flow are
where V is the velocity, P′ and ρ′ are pressure and density fluctuations, respectively, ρ is the background density, g is the gravity vector, ν is the kinematic viscosity, f is a forcing term, s is the specific entropy, and κ is the thermal diffusivity. If one neglects the largescale meridional circulation, velocity can be split into a mean rotation part and a fluctuating part:
where r, θ, and φ are the standard spherical coordinates, (e_{r}, e_{θ}, e_{φ}) is the associated orthonormal basis, Ω is the rotation rate, and u is the velocity fluctuation. The continuity Eq. (1) implies that
The lefthand side of the momentum Eq. (2) can be developed into
where s is the distance to the rotation axis and e_{s} the associated unit vector. The second term, which is the centrifugal acceleration, will be neglected in the following. Then, in the frame that locally rotates at the same rate as the fluid, Eq. (2) becomes
where S = r sin θ∇Ω is the shear rate.
Similarly, the entropy can be split into a mean part ⟨s⟩ and a fluctuating one s′:
Equation (3) yields
Nonlinear and dissipative terms are usually difficult to model analytically in turbulent flows. To simplify the problem, we use here a relaxation approximation, also known as τ approximation, which assumes that in a stationary steady state, the main effect of these terms is to tend the flow to relax with a given time constant τ. Equations (7) and (9) can then be approximated by
For a complete treatment of thermal diffusion, we refer the reader to Park et al. (2020, 2021).
These two equations can be combined into a single one by expressing entropy fluctuations as a function of density fluctuations:
where c_{P} is the specific heat capacity at constant pressure. This yields
where the BruntVäisälä frequency N is defined by
Using the Fourier decomposition
the continuity Eq. (5) becomes
and pressure fluctuations can be eliminated from Eq. (13):
where , , , , , and . This way of writing the forcing term means that we assume that it comes from a preexisting background turbulent flow. The previous equation can be written as a matrix relation
where ℳ is the matrix
and .
When ℳ is inversible, Eq. (18) implies
Using the fact that also verifies the continuity equation, the previous equation reduces to
where
is the determinant of ℳ, , 𝒟 is the matrix given by
and δ_{ij} and ε_{ijk} are the Kronecker and LeviCivita tensors, respectively.
The spectrum of the background turbulence is assumed to be characterised by the correlation tensor , where ⟨⟩ denotes a statistical average. In the continuity of Mathis et al. (2018), we assume a purely horizontal background turbulence, for which the Fourier transform of is
where E(k) is the kinetic energy spectrum (see e.g. Davidson 2004) such that
The turbulent spectrum of the sheared turbulence is deduced from that of the background turbulence using the relation
After some algebra, this leads to
where
3. Mean correlations
Mean velocity correlations can be computed using the relation
Stellar radiative zones are usually strongly stratified. It implies that (for example, the typical value of the ratio N/(2Ω) is around 500 for the Sun; see also André et al. 2018, 2019, for the evolution of this ratio along stellar evolution for low and intermediatemass stars), , and . As a first approximation,
where μ = cos α, and (with these notations, σ = cos θ cos α − sin θ sin α cos β). For the same reason,
Using parity arguments and integrating over k (see Appendix A), this leads, when , to
This term corresponds to the Reynolds stress in the meridional plane. Similarly,
The first two terms correspond to the vertical and horizontal transport of angular momentum, respectively. The other two describe the horizontal turbulent velocity scales.
We note that when tends towards zero, the scalings for ⟨u_{r}u_{θ}⟩ and ⟨u_{θ}u_{φ}⟩ (in Eqs. (34) and (36), respectively) vanish, while the others tend towards finite values. These values, as well as the zero limit of ⟨u_{θ}u_{φ}⟩, are consistent with the results of KB12. In contrast, KB12 found a nonzero value for ⟨u_{r}u_{θ}⟩:
which is of lower order in than the new scaling of Eq. (34).
The case of is more complicated, since it contains three terms that come from
and a priori have the same order of magnitude. The first term leads to
which is consistent with the scaling derived in KB12 and M+18. The second term vanishes for parity reasons. The third term leads to
When , this reduces to
which is small compared to . In total,
4. Interpretation
Interestingly, none of the terms computed in Sect. 3 explicitely depend on the radial shear , although it was included in the governing equations. In contrast, all share a common feature: they have a denominator in . This implies that when gets close to 2, the turbulence is enhanced. Physically, this might be related to the limit of the RayleighTaylor instability. The RayleighTaylor stability criterion is
Assuming that Ω > 0, this criterion can be rewritten
When 2Ω + sin θS_{r} > 0, the last inequality is verified for all positive values of cos θS_{θ}, but only for negative values verifying
By comparison with , this suggests that
This expression is similar to τ = (2Ω + S_{r})^{−1}, proposed by M+18. However, the presence of cos θ at the numerator means that τ would vanish at the equator, in which case the scalings derived in Sect. 3 are not valid there. Besides, the comparison with the RayleighTaylor instability stands only in the special case where 2Ω + sin θS_{r} > 0 and cos S_{θ} < 0. Thus, the criterion could be linked to some other hydrodynamical instability instead, such as the inertial or the inflectional instability studied by Park et al. (2020, 2021).
The form of the expression of ⟨u_{θ}u_{φ}⟩ in Eq. (36), which describes the horizontal transport of angular momentum, suggests that this flux term is viscous, with a turbulent viscosity coefficient equal to
We note that this coefficient has no explicit dependence on the BruntVäisälä frequency or on the vertical shear. However, τ most likely depends on S_{r}, and not necessarily on N. Thus, the viscosity coefficient probably depends on both S_{r} and S_{θ}. The dependence on ⟨u^{2}⟩^{(0)} is problematic for the use in stellar evolution codes, since there is no easy way to reliably prescribe this quantity. One way to roughly estimate it is to assume that horizontal motions generated by the background turbulence are global. This can be written
where R is the radial extent of the radiative zone. Introducing this expression in Eq. (48) yields
In the case where the dynamics is dominated by the horizontal shear, assuming τ ∼ 1/S_{θ} further leads to a horizontal viscosity coefficient that is proportional to the horizontal shear rate, which is consistent with the viscosity coefficient derived by Park et al. (2021). We note that, in contrast, Garaud (2001) found that the horizontal shear instability can lead to an antidiffusive horizontal transport of angular momentum.
Although the expression in Eq. (48) is not directly usable in stellar evolution codes, it can be used to express other flux terms as a function of ν_{h}. Thus,
In the present work, we consider the case of a horizontally isotropic turbulence (i.e. with identical statistical properties in the latitudinal and azimuthal directions). This is rigourously justified for the radiative zone of slowly rotating stars, in which turbulent structures would have the shape of a horizontal pancake (see e.g. Fig. 1 of M+18). Because of that, there is no horizontal Λ effect. In rapidly rotating stars, turbulent structures are elongated along the rotation axis, so the horizontal isotropy is broken (see e.g. Rüdiger 1980) and a nonzero horizontal flux of angular momentum due to the Λ effect is present. However, in stellar radiative zones, provided that the BruntVäisälä frequency is still much larger than the rotation rate, the horizontal flux due to the Λ effect remains negligle compared to the leading diffusive term, which reduces to Eq. (48) in the horizontally symmetric case.
In contrast with the horizontal transport, the vertical transport may not be viscous, since the term ⟨u_{r}u_{φ}⟩ given in Eq. (35) is not explicitly proportional to S_{r}. Instead, it is proportional to the rotation rate Ω, which is a characteristic of the Λ effect. However, when the dynamics is dominated by the vertical shear, one can assume that τ ∼ 1/S_{r}, in which case Eq. (35) becomes
This then leads to a vertical turbulent viscosity coefficient
which is negative. This corresponds to an antidiffusive transport of angular momentum. Cope et al. (2020) and Garaud (2020) estimated the vertical turbulent diffusion coefficient associated with the horizontal shear instability using scaling laws but they did not provide any proof that the transport of chemical elements is indeed diffusive. In addition, they did not discuss the transport of angular momentum. Using the approximation of Eq. (49) leads to
The same approximation can be used in the nonviscous case to express the flux from Eq. (35):
Following KB12 and M+18, we estimate the turbulent diffusion coefficients using a mixinglength approximation and , where . This allows us to estimate the ratio between vertical and horizontal diffusion coefficients:
which has the same scaling with N as found by M+18. Besides, in the limit where , Eq. (57) reduces to
which also has the same scaling with τ as found by M+18.
Moreover, in contrast with M+18, where ⟨u_{θ}u_{φ}⟩ was zero and no horizontal turbulent viscosity coefficient could be directly derived, we can now compare horizontal viscosity and diffusion coefficients. We thus find that D_{h} ≃ 4ν_{h}. This suggests that the common assumption that they are equal may be qualitatively valid, but given the uncertainties introduced by the approximation we used, we cannot really conclude on the quantitative validity. In the vertical direction, we cannot directly compute the turbulent viscosity coefficient, so there is no evidence of such an equality.
5. Nonviscous transport in stellar evolution codes
If a prescription of ⟨u^{2}⟩^{(0)} is found, our expressions for the flux terms in Eqs. (35) and (36) can be used in stellar evolution calculations. For nonviscous transport, the evolution of angular momentum is governed by
where u(r, θ) is the largescale velocity field of the meridional circulation in stellar radiative zones. 1D stellar evolution codes that take rotation into account usually are based on the shellular approximation, which assumes that , where
and (see e.g. Zahn 1992). In this approximation, the meridional circulation reads
where P_{2} is the Legendre polynomial of degree 2. Averaging Eq. (59) over isobars yields
where
We note that the term of horizontal transport is no longer present since it has a zero average. We showed in Sect. 4 that the horizontal transport can be seen as viscous. Therefore, the equation for the fluctuations is the same as in Mathis et al. (2004), for example:
where and
6. Conclusion
In the current paper, we investigated the impact of horizontal shear on the anisotropy of the transport in stably stratified, rotating stellar radiative zones. We were able to derive new scalings for mean velocity correlations that confirmed that the main effect is to significantly enhance turbulent motions. Thus, horizontal shear can be considered as an additional source of transport of angular momentum and chemical elements both in the horizontal and vertical directions.
In the regime where the stratification has a much stronger effect than rotation, we found that the horizontal transport of angular momentum is diffusive, while the vertical transport is dominated by the Λ effect. Additionally, the horizontal turbulent viscosity and diffusion coefficients are proportional to each other, which partly justifies the approximation that they are equal, wellspread in stellar evolution codes. The expressions of all the mean Reynolds stresses computed in the current paper depend on a property of the background turbulent flow, which is unknown a priori. Nevertheless, if one of the components is known by another way (for example using a phenomenological model for the horizontal turbulent viscosity coefficient), this allows one to compute fluxes of angular momentum. We finally show how the nonviscous transport of angular momentum predicted by our model should be implemented in stellar evolution codes to go beyond the current diffusive and viscous formalism.
Although the relaxation approximation has been validated in special cases by numerical simulations (see e.g. KB12), its validity in the configuration considered in the present work is uncertain. In particular, the dependence of the nonlinear time τ on the physical parameters is undetermined. In M+18, three different physically motivated but phenomenological expressions were proposed. Here, we derived from more rigourous arguments a new expression for the nonlinear turbulent time scale which is close to one of those proposed in M+18, τ = (2Ω + S_{r})^{−1}, which corresponds to a turbulence dominated by the effects of rotation and vertical shear. Ultimately, direct numerical simulations should be performed to test the formalism and the new prescriptions in the presence of a general shear.
The results presented in this work have been obtained in the Boussinesq approximation, which neglects density fluctuations except in the buoyancy term. This treatment is valid as long as turbulent structures are small compared to the density scale height. When this condition is not satisfied, for example in regions with strong density gradients, similar calculations should be performed in the more general framework of the anelastic approximation, which accounts for density gradients, but filters pressure waves out.
Acknowledgments
V.P. and S.M. acknowledge support from the European Research Council through ERC grant SPIRE 647383 and from the CNES PLATO/GOLF grant at CEASaclay. The authors thank the anonymous referee for their constructive comments, which have helped us improve the manuscript.
References
 Aerts, C., Van Reeth, T., & Tkachenko, A. 2017, ApJ, 847, L7 [NASA ADS] [CrossRef] [Google Scholar]
 André, Q., Mathis, S., & Amard, L. 2018, SF2A2018: Proceedings of the Annual meeting of the French Society of Astronomy and Astrophysics, Di [Google Scholar]
 André, Q., Mathis, S., & Amard, L. 2019, https://doi.org/10.5281/zenodo.2545306 [Google Scholar]
 Beck, P. G., Montalban, J., Kallinger, T., et al. 2012, Nature, 481, 55 [NASA ADS] [CrossRef] [Google Scholar]
 Benomar, O., Takata, M., Shibahashi, H., Ceillier, T., & García, R. A. 2015, MNRAS, 452, 2654 [NASA ADS] [CrossRef] [Google Scholar]
 Brown, T. M., ChristensenDalsgaard, J., Dziembowski, W. A., et al. 1989, ApJ, 343, 526 [NASA ADS] [CrossRef] [Google Scholar]
 Cantiello, M., Mankovich, C., Bildsten, L., ChristensenDalsgaard, J., & Paxton, B. 2014, ApJ, 788, 93 [Google Scholar]
 Ceillier, T., Eggenberger, P., García, R. A., & Mathis, S. 2013, A&A, 555, A54 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Cope, L., Garaud, P., & Caulfield, C. P. 2020, J. Fluid Mech., 903, A1 [Google Scholar]
 Davidson, P. A. 2004, Turbulence: an Introduction for Scientists and Engineers (Oxford University Press) [Google Scholar]
 Deheuvels, S., García, R. A., Chaplin, W. J., et al. 2012, ApJ, 756, 19 [NASA ADS] [CrossRef] [Google Scholar]
 Deheuvels, S., Doğan, G., Goupil, M. J., et al. 2014, A&A, 564, A27 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Deheuvels, S., Ballot, J., Beck, P. G., et al. 2015, A&A, 580, A96 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Eggenberger, P., Montalbán, J., & Miglio, A. 2012, A&A, 544, L4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Fossat, E., Boumier, P., Corbard, T., et al. 2017, A&A, 604, A40 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Gagnier, D., & Garaud, P. 2018, ApJ, 862, 36 [NASA ADS] [CrossRef] [Google Scholar]
 Garaud, P. 2001, MNRAS, 324, 68 [Google Scholar]
 Garaud, P. 2020, ApJ, 901, 146 [Google Scholar]
 Garaud, P., & Kulenthirarajah, L. 2016, ApJ, 821, 49 [Google Scholar]
 Garaud, P., Gagnier, D., & Verhoeven, J. 2017, ApJ, 837, 133 [NASA ADS] [CrossRef] [Google Scholar]
 García, R. A., TurckChièze, S., JiménezReyes, S. J., et al. 2007, Science, 316, 1591 [Google Scholar]
 Gehan, C., Mosser, B., Michel, E., Samadi, R., & Kallinger, T. 2018, A&A, 616, A24 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Heger, A., Woosley, S. E., & Spruit, H. C. 2005, ApJ, 626, 350 [NASA ADS] [CrossRef] [Google Scholar]
 Hermes, J. J., Gänsicke, B. T., Kawaler, S. D., et al. 2017, ApJS, 232, 23 [Google Scholar]
 Hirschi, R., & Maeder, A. 2010, A&A, 519, A16 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Käpylä, P. J. 2019, A&A, 622, A195 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Kitchatinov, L. L., & Brandenburg, A. 2012, Astron. Nachr., 333, 230 [NASA ADS] [CrossRef] [Google Scholar]
 Kurtz, D. W., Saio, H., Takata, M., et al. 2014, MNRAS, 444, 102 [Google Scholar]
 Maeder, A. 2003, A&A, 399, 263 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Maeder, A. 2009, Physics, Formation and Evolution of Rotating Stars (Berlin Heidelberg: Springer) [Google Scholar]
 Maeder, A., & Meynet, G. 1996, A&A, 313, 140 [NASA ADS] [Google Scholar]
 Maeder, A., & Zahn, J.P. 1998, A&A, 334, 1000 [NASA ADS] [Google Scholar]
 Marques, J. P., Goupil, M. J., Lebreton, Y., et al. 2013, A&A, 549, A74 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mathis, S., & Zahn, J. P. 2004, A&A, 425, 229 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mathis, S., Palacios, A., & Zahn, J.P. 2004, A&A, 425, 243 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mathis, S., Prat, V., Amard, L., et al. 2018, A&A, 620, A22 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mosser, B., Goupil, M. J., Belkacem, K., et al. 2012, A&A, 548, A10 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Murphy, S. J., Fossati, L., Bedding, T. R., et al. 2016, MNRAS, 459, 1201 [Google Scholar]
 Ouazzani, R. M., Marques, J. P., Goupil, M. J., et al. 2019, A&A, 626, A121 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Park, J., Prat, V., & Mathis, S. 2020, A&A, 635, A133 [EDP Sciences] [Google Scholar]
 Park, J., Prat, V., Mathis, S., & Bugnet, L. 2021, A&A, 646, A64 [EDP Sciences] [Google Scholar]
 Prat, V., & Lignières, F. 2013, A&A, 551, L3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Prat, V., & Lignières, F. 2014, A&A, 566, A110 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Prat, V., Guilet, J., Viallet, M., & Müller, E. 2016, A&A, 592, A59 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Rüdiger, G. 1980, Geophys. Astro. Fluid, 16, 239 [Google Scholar]
 Rüdiger, G. 1989, Differential Rotation and Stellar Convection (Gordon & Breach Science Publishers) [Google Scholar]
 Saio, H., Kurtz, D. W., Takata, M., et al. 2015, MNRAS, 447, 3264 [Google Scholar]
 Suijs, M. P. L., Langer, N., Poelarends, A.J., et al. 2008, A&A, 481, L87 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Talon, S., & Zahn, J.P. 1997, A&A, 317, 749 [Google Scholar]
 Thompson, M. J., ChristensenDalsgaard, J., Miesch, M. S., & Toomre, J. 2003, ARA&A, 41, 599 [NASA ADS] [CrossRef] [Google Scholar]
 Triana, S. A., Moravveji, E., Pápics, P. I., et al. 2015, ApJ, 810, 16 [NASA ADS] [CrossRef] [Google Scholar]
 Triana, S. A., Corsaro, E., De Ridder, J., et al. 2017, A&A, 602, A62 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Van Reeth, T., Tkachenko, A., & Aerts, C. 2016, A&A, 593, A120 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Van Reeth, T., Mombarg, J. S. G., Mathis, S., et al. 2018, A&A, 618, A24 [NASA ADS] [EDP Sciences] [Google Scholar]
 Zahn, J. P. 1983, in SaasFee Advanced Course 13: Astrophysical Processes in Upper Main Sequence Stars, eds. A. N. Cox, S. Vauclair, & J. P. Zahn, 253 [Google Scholar]
 Zahn, J.P. 1992, A&A, 265, 115 [NASA ADS] [Google Scholar]
Appendix A: Integrals
The derivation of Eq. (34) from Eq. (33) requires the computation of the double integral
The first and third terms of the integrand are antisymmetric with respect to α = π/2, so the corresponding contribution to the integral is zero. The remaining term is then proportional to
Using the periodicity of the integrand and the variables t = tanβ and , one obtains
These integrals are convergent only when . Using , the integral now reads
Similarly, the derivation of Eq. (35) requires the computation of
and the derivation of Eq. (37) requires the computation of
It follows that
which is required for the derivation of Eq. (40). Finally, the derivation of Eq. (41) requires the computation of
Using the partial fraction decomposition
one obtains
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.