We investigate the linear properties of the steady stress-driven spin-down flow of a viscous fluid inside a spherical shell, both within the incompressible and anelastic approximations, and in the asymptotic limit of small viscosities. Neglecting buoyancy force, we find a steady-state to settle in a viscous time for the incompressible flow. From boundary layer analysis, we derive an analytical geostrophic solution for the 3D incompressible steady flow, inside and outside a cylinder $\mathcal{C}$ that is tangent to the inner shell. The meridional circulation associated with the spin-down flow in these two regions scales with the Ekman number $E$. A power expansion of the fields with the small parameter $E^{1/3}$ corresponding to the width of the inner Stewartson shear layer laying along the tangent cylinder, reveals it to be the site of an upwelling flow of amplitude $O(E^{1/3})$. Despite its narrowness, this shear layer concentrates most of the global meridional kinetic energy of the system. Furthermore, a stable stratification does not perturb the spin-down flow provided the Prandtl number is small enough. If this is not the case, the Stewartson layer disappears and meridional circulation is confined within the thermal layers. The scalings for the amplitude of the anelastic secondary flow have been found to be the same as for the incompressible flow in all three regions, at the lowest order. However, because the velocity no longer conforms the Taylor-Proudman theorem, its shape differs outside the tangent cylinder $\mathcal{C}$, that is, where differential rotation takes place. Finally, we find the settling of the steady-state to be reached on a viscous time for the weakly, strongly and thermally unstratified incompressible flows. Large density variations relevant to astro- and geophysical systems, tend to slightly shorten the transient.

The understanding of the evolution of early-type stars is tightly related to that of the effects of rapid rotation. For massive stars, rapid rotation combines with their strong radiation-driven wind. The aim of this paper is to investigate two questions that are prerequisite to the study of the evolution of massive rapidly rotating stars: (i) What is the critical angular velocity of a star when radiative acceleration is significant in its atmosphere? (ii) How do mass and angular momentum loss depend on the rotation rate? To investigate fast rotation, which makes stars oblate, we used the 2D ESTER models and a simplified approach, the $\omega$-model, which gives the latitudinal dependence of the radiative flux in a centrifugally flattened radiative envelope. We find that radiative acceleration only mildly influences the critical angular velocity, at least for stars with masses lower than 40 Msun. We explain this mild reduction of the critical angular velocity compared to the classical Keplerian angular velocity by the combined effects of gravity darkening and a reduced equatorial opacity that is due to the centrifugal acceleration. To answer the second question, we first devised a model of the local surface mass flux, which we calibrated with previously developed 1D models. The discontinuity (the so-called bi-stability jump) included in the $\dot{M}-T_{\rm eff}$ relation of 1D models means that the mass flux of a fast-rotating star is controlled by either a single wind or a two-wind regime. Mass and angular momentum losses are strong around the equator if the star is in the two-wind regime. We also show that the difficulty of selecting massive stars that are viewed pole-on makes detecting the discontinuity in the relation between mass loss and effective temperature also quite challenging.

The understanding of the rotational evolution of early-type stars is deeply related to that of anisotropic mass and angular momentum loss. In this paper, we aim to clarify the rotational evolution of rapidly rotating early-type stars along the main sequence (MS). We have used the 2D ESTER code to compute and evolve isolated rapidly rotating early-type stellar models along the MS, with and without anisotropic mass loss. We show that stars with $Z=0.02$ and masses between $5$ and $7~M_\odot$ reach criticality during the main sequence provided their initial angular velocity is larger than 50% of the Keplerian one. More massive stars are subject to radiation-driven winds and to an associated loss of mass and angular momentum. We find that this angular momentum extraction from the outer layers can prevent massive stars from reaching critical rotation and greatly reduce the degree of criticality at the end of the MS. Our model includes the so-called bi-stability jump of the $\dot{M}-T_{\rm eff}$ relation of 1D-models. This discontinuity now shows up in the latitude variations of the mass-flux surface density, endowing rotating massive stars with either a single-wind regime (no discontinuity) or a two-wind regime (a discontinuity). In the two-winds-regime, mass loss and angular momentum loss are strongly increased at low latitudes inducing a faster slow-down of the rotation. However, predicting the rotational fate of a massive star is difficult, mainly because of the non-linearity of the phenomena involved and their strong dependence on uncertain prescriptions. Moreover, the very existence of the bi-stability jump in mass-loss rate remains to be substantiated by observations.

This paper continues the systematic investigation of diffusive shear instabilities initiated in Part I of this series. In this work, we primarily focus on quantifying the impact of non-local mixing, which is not taken into account in Zahn's mixing model \citep{Zahn92}. We present the results of direct numerical simulations in a new model setup designed to contain coexisting laminar and turbulent shear layers. As in Part I, we use the Low P\'eclet Number approximation of \citet{Lign1999} to model the evolution of the perturbations. Our main findings are twofold. First, turbulence is not necessarily generated whenever Zahn's nonlinear criterion \citep{Zahn1974} $J{\rm Pr} < (J{\rm Pr})_c$ is satisfied, where $J=N^2/S^2$ is the local gradient Richardson number, ${\rm Pr} = \nu/ \kappa_T$ is the Prandtl number, and $(J{\rm Pr})_c \simeq 0.007$. We have demonstrated that the presence or absence of turbulent mixing in this limit hysteretically depends on the history of the shear layer. Second, Zahn's nonlinear instability criterion only approximately locates the edge of the turbulent layer, and mixing beyond the region where $J{\rm Pr} < (J{\rm Pr})_c$ can also take place in a manner analogous to convective overshoot. We found that the turbulent kinetic energy decays roughly exponentially beyond the edge of the shear-unstable region, on a lengthscale $\delta$ that is directly proportional to the scale of the turbulent eddies, which are themselves of the order of the Zahn scale (see Part I). Our results suggest that mixing by diffusive shear instabilities should be modeled with more care than is currently standard in stellar evolution codes.

Shear-induced turbulence could play a significant role in mixing momentum and chemical species in stellar radiation zones, as discussed by Zahn (1974). In this paper we analyze the results of direct numerical simulations of stratified plane Couette flows, in the limit of rapid thermal diffusion, to measure the turbulent diffusivity and turbulent viscosity as a function of the local shear and the local stratification. We find that the stability criterion proposed by Zahn (1974), namely that the product of the gradient Richardson number and the Prandtl number must be smaller than a critical values $(J\Pr)_c$ for instability, adequately accounts for the transition to turbulence in the flow, with $(J\Pr)_c \simeq 0.007$. This result recovers and confirms the prior findings of Prat et al. (2016). Zahn's model for the turbulent diffusivity and viscosity (Zahn 1992), namely that the mixing coefficient should be proportional to the ratio of the thermal diffusivity to the gradient Richardson number, does not satisfactorily match our numerical data when applied as is. It fails (as expected) in the limit of large stratification where the Richardson number exceeds the aforementioned threshold for instability, but it also fails in the limit of low stratification where the turbulent eddy scale becomes limited by the computational domain size. We propose a revised model for turbulent mixing by diffusive stratified shear instabilities, that now properly accounts for both limits, fits our data satisfactorily, and recovers Zahn's 1992 model in the limit of large Reynolds numbers.