No Arabic abstract
Rotational mixing, the key process in stellar evolution, transports angular momentum and chemical elements in stellar radiative zones. In the past two decades, an emphasis has been placed on the turbulent transport induced by the vertical shear instability. However, instabilities arising from horizontal shear and the strength of the anisotropic turbulent transport that they may trigger remain relatively unexplored. This paper investigates the combined effects of stable stratification, rotation, and thermal diffusion on the horizontal shear instabilities in the context of stellar radiative zones. The eigenvalue problem describing linear instabilities of a flow with a hyperbolic-tangent horizontal shear profile was solved numerically for a wide range of parameters. As a first step, we consider a polar $f$-plane where the gravity and rotation vector are aligned. Two types of instabilities are identified: the inflectional and inertial instabilities. The inflectional instability that arises from the inflection point is the most unstable when at a zero vertical wavenumber and a finite wavenumber in the streamwise direction along the imposed-flow direction. The three-dimensional inflectional instability is destabilized by stratification, while it is stabilized by thermal diffusion. The inertial instability is rotationally driven, and a WKBJ analysis reveals that its growth rate reaches the maximum $sqrt{f(1-f)}$ in the inviscid limit as the vertical wavenumber goes to infinity, where $f$ is the dimensionless Coriolis parameter. The inertial instability for a finite vertical wavenumber is stabilized as the stratification increases, whereas it is destabilized by the thermal diffusion. Furthermore, we found a self-similarity in both the inflectional and inertial instabilities based on the rescaled parameter $PeN^2$ with the P{e}clet number $Pe$ and the Brunt-V{a}is{a}l{a} frequency $N$.
Stellar interiors are the seat of efficient transport of angular momentum all along their evolution. Understanding the dependence of the turbulent transport triggered by the shear instabilities due to the differential rotation in stellar radiation zones is mandatory. Indeed, it constitutes one of the cornerstones of the rotational transport and mixing theory which is implemented in stellar evolution codes to predict the rotational and chemical evolutions of stars. We investigate horizontal shear instabilities in stellar radiation zones by considering the full Coriolis acceleration with both the dimensionless horizontal component $tilde{f}$ and the vertical component $f$. We performed a linear stability analysis for a horizontal shear flow with a hyperbolic tangent profile, both numerically and asymptotically using the WKBJ approximation. As in the traditional approximation, we identified the inflectional and inertial instabilities. The inflectional instability is destabilized as $tilde{f}$ increases and its maximum growth rate increases significantly, while the thermal diffusivity stabilizes the inflectional instability similarly to the traditional case. The inertial instability is also strongly affected; for instance, the inertially unstable regime is extended in the non-diffusive limit as $0<f<1+tilde{f}^{2}/N^{2}$, where $N$ is the dimensionless Brunt-Vaisala frequency. More strikingly, in the high-thermal-diffusivity limit, it is always inertially unstable at any colatitude $theta$ except at the poles (i.e., $0^{circ}<theta<180^{circ}$). Using the asymptotic and numerical results, we propose a prescription for the effective turbulent viscosities induced by the instabilities to be possibly used in stellar evolution models. The characteristic time of this turbulence is short enough to redistribute efficiently angular momentum and mix chemicals in the radiation zones.
We examine the MHD instabilities arising in the radiation zone of a differentially rotating star, in which a poloidal field of fossil origin is sheared into a toroidal field. We focus on the non-axisymmetric instability that affects the toroidal magnetic field in a rotating star, which was first studied by Pitts and Tayler in the non-dissipative limit. According to Spruit, it could also drive a dynamo. The Pitts & Tayler instability is manifestly present in our simulations, with its conspicuous m=1 dependence in azimuth. But its analytic treatment used so far is too simplified to be applied to the real stellar situation. Although the instability generated field reaches an energy comparable to that of the mean poloidal field, that field seems unaffected by the instability: it undergoes Ohmic decline, and is neither eroded nor regenerated by the instability. The toroidal field is produced by shearing the poloidal field and it draws its energy from the differential rotation. The small scale motions behave as Alfven waves; they cause negligible eddy-diffusivity and contribute little to the net transport of angular momentum. In our simulations we observe no sign of dynamo action, of either mean field or fluctuation type, up to a magnetic Reynolds number of 10^5. However the Pitts & Tayler instability is sustained as long as the differential rotation acting on the poloidal field is able to generate a toroidal field of sufficient strength.
We present results from two-dimensional, general relativistic, viscous, radiation hydrodynamic numerical simulations of Shakura-Sunyaev thin disks accreting onto stellar mass Schwarzschild black holes. We consider cases on both the gas- and radiation-pressure-dominated branches of the thermal equilibrium curve, with mass accretion rates spanning the range from $dot{M} = 0.01 L_mathrm{Edd}/c^2$ to $10 L_mathrm{Edd}/c^2$. The simulations directly test the stability of this standard disk model on the different branches. We find clear evidence of thermal instability for all radiation-pressure-dominated disks, resulting universally in the vertical collapse of the disks, which in some cases then settle onto the stable, gas-pressure-dominated branch. Although these results are consistent with decades-old theoretical predictions, they appear to be in conflict with available observational data from black hole X-ray binaries. We also find evidence for a radiation-pressure-driven instability that breaks the unstable disks up into alternating rings of high and low surface density on a timescale comparable to the thermal collapse. Since radiation is included self-consistently in the simulations, we are able to calculate lightcurves and power density spectra (PDS). For the most part, we measure radiative efficiencies (ratio of luminosity to mass accretion rate) close to 6%, as expected for a non-rotating black hole. The PDS appear as broken power laws, with a break typically around 100 Hz. There is no evidence of significant excess power at any frequencies, i.e. no quasi-periodic oscillations are observed.
Radiation-dust driven outflows, where radiation pressure on dust grains accelerates gas, occur in many astrophysical environments. Almost all previous numerical studies of these systems have assumed that the dust was perfectly-coupled to the gas. However, it has recently been shown that the dust in these systems is unstable to a large class of resonant drag instabilities (RDIs) which de-couple the dust and gas dynamics and could qualitatively change the nonlinear outcome of these outflows. We present the first simulations of radiation-dust driven outflows in stratified, inhomogeneous media, including explicit grain dynamics and a realistic spectrum of grain sizes and charge, magnetic fields and Lorentz forces on grains (which dramatically enhance the RDIs), Coulomb and Epstein drag forces, and explicit radiation transport allowing for different grain absorption and scattering properties. In this paper we consider conditions resembling giant molecular clouds (GMCs), HII regions, and distributed starbursts, where optical depths are modest ($lesssim 1$), single-scattering effects dominate radiation-dust coupling, Lorentz forces dominate over drag on grains, and the fastest-growing RDIs are similar, such as magnetosonic and fast-gyro RDIs. These RDIs generically produce strong size-dependent dust clustering, growing nonlinear on timescales that are much shorter than the characteristic times of the outflow. The instabilities produce filamentary and plume-like or horsehead nebular morphologies that are remarkably similar to observed dust structures in GMCs and HII regions. Additionally, in some cases they strongly alter the magnetic field structure and topology relative to filaments. Despite driving strong micro-scale dust clumping which leaves some gas behind, an order-unity fraction of the gas is always efficiently entrained by dust.
Helicity and alpha effect driven by the nonaxisymmetric Tayler instability of toroidal magnetic fields in stellar radiation zones are computed. In the linear approximation a purely toroidal field always excites pairs of modes with identical growth rates but with opposite helicity so that the net helicity vanishes. If the magnetic background field has a helical structure by an extra (weak) poloidal component then one of the modes dominates producing a net kinetic helicity anticorrelated to the current helicity of the background field. The mean electromotive force is computed with the result that the alpha effect by the most rapidly growing mode has the same sign as the current helicity of the background field. The alpha effect is found as too small to drive an alpha^{2} dynamo but the excitation conditions for an alphaOmega dynamo can be fulfilled for weak poloidal fields. Moreover, if the dynamo produces its own alpha effect by the magnetic instability then problems with its sign do not arise. For all cases, however, the alpha effect shows an extremely strong concentration to the poles so that a possible alphaOmega dynamo might only work at the polar regions. Hence, the results of our linear theory lead to a new topological problem for the existence of large-scale dynamos in stellar radiation zones on the basis of the current-driven instability of toroidal fields.