Subcritical dynamos in shear flows

Identifying generic physical mechanisms responsible for the generation of magnetic fields and turbulence in differentially rotating flows is fundamental to understand the dynamics of astrophysical objects such as accretion disks and stars. In this paper, we discuss the concept of subcritical dynamo action and its hydrodynamic analogue exemplified by the process of nonlinear transition to turbulence in non‐rotating wall‐bounded shear flows. To illustrate this idea, we describe some recent results on nonlinear hydrodynamic transition to turbulence and nonlinear dynamo action in rotating shear flows pertaining to the problem of turbulent angular momentum transport in accretion disks. We argue that this concept is very generic and should be applicable to many astrophysical problems involving a shear flow and non‐axisymmetric instabilities of shearinduced axisymmetric toroidal velocity or magnetic fields, such as Kelvin‐Helmholtz, magnetorotational, Tayler or global magnetoshear instabilities. In the light of several recent numerical results, we finally suggest that, similarly to a standard linear instability, subcritical MHD dynamo processes in high‐Reynolds number shear flows could act as a large‐scale driving mechanism of turbulent flows that would in turn generate an independent small‐scale dynamo. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Turbulence and its parameterization in accretion discs

Accretion disc turbulence is investigated in the framework of the shearing box approximation. The turbulence is either driven by the magneto‐rotational instability or, in the non‐magnetic case, by an explicit and artificial forcing term in the momentum equation. Unlike the magnetic case, where most of the dissipation occurs in the disc corona, in the forced hydrodynamic case most of the dissipation occurs near the midplane. In the hydrodynamic case evidence is presented for the stochastic excitation of epicycles. When the vertical and radial epicyclic frequencies are different (modeling the properties around rotating black holes), the beat frequency between these two frequencies appear to show up as a peak in the temporal power spectrum in some cases. Finally, the full turbulent resistivity tensor is determined and it is found that, if the turbulence is driven by a forcing term, the signs of its off‐diagonal components are such that this effect would not be capable of dynamo action by the shear–current effect. (© 2005 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Linear stability analysis of the Hall magnetorotational instability in a spherical domain

We investigate the stability of the Hall‐MHD system and determine its importance for neutron stars at their birth, when they still consist of differentially rotating plasma permeated by extremely strong magnetic fields. We solve the linearised HallMHD equations in a spherical shell threaded by a homogeneous magnetic field. With the fluid/flow coupling and the Hall effect included, the magnetorotational instability and the Hall effect are both acting together. Results differ for magnetic fields aligned with the rotation axis and anti‐parallel magnetic fields. For a positive alignment of the magnetic field the instability grows on a rotational time‐scale for any sufficiently large magnetic Reynolds number. Even the magnetic fields which are stable against the MRI due to the magnetic diffusion are now susceptible to the shear‐Hall instability. In contrast, the negative alignment places strong restrictions on the growth and the magnitude of the fields, hindering the effectiveness of the Hall‐MRI. While non‐axisymmetric modes of the MRI can be suppressed by strong enough rotation, there is no such restriction when the Hall effect is present. The implications for the magnitude and the topology of the magnetic field of a young neutron star may be significant (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

The response of a turbulent accretion disc to an imposed epicyclic shearing motion

We excite an epicyclic motion, the amplitude of which depends on the vertical position, z , in a simulation of a turbulent accretion disc. An epicyclic motion of this kind may be caused by a warping of the disc. By studying how the epicyclic motion decays, we can obtain information about the interaction between the warp and the disc turbulence. A high-amplitude epicyclic motion decays first by exciting inertial waves through a parametric instability, but its subsequent exponential damping may be reproduced by a turbulent viscosity. We estimate the effective viscosity parameter, α _v, pertaining to such a vertical shear. We also gain new information on the properties of the disc turbulence in general, and measure the usual viscosity parameter, α _h, pertaining to a horizontal (Keplerian) shear. We find that, as is often assumed in theoretical studies, α _v is approximately equal to α _h and both are much less than unity, for the field strengths achieved in our local box calculations of turbulence. In view of the smallness (∼0.01) of α _v and α _h we conclude that for β p _gas p _mag∼10 the time-scale for diffusion or damping of a warp is much shorter than the usual viscous time-scale. Finally, we review the astrophysical implications.     

Latitudinal shear instability in the solar tachocline

In an attempt to explain the observed rotation profile in the solar radiative zone and the tachocline, Spiegel & Zahn proposed a model based on anisotropic turbulent angular momentum transport. Although very successful in reproducing some of the features of the solar tachocline, their model assumes without verification that the origin of the turbulence could be caused by latitudinal shear instability. This paper studies the weakly non-linear evolution of two-dimensional shear instability, in which the interaction between the global rotation profile and the Reynolds stresses can be described self-consistently. Provided that the initial rotation profile is sufficiently close to marginal stability (which is the case of the solar tachocline), the instability is shown to saturate and to relax to a marginally stable state, which differs very little from the observed rotation profile. It is therefore likely that the tachocline is in a state of marginal stability with respect to latitudinal shear instability, and shows that angular momentum transport in the tachocline is unlikely to be caused by shear-induced turbulence.     

Turbulent stresses as a function of shear rate in a local disk model

We present local numerical models of accretion disk turbulence driven by the magnetorotational instability with varying shear rate. The resulting turbulent stresses are compared with predictions of a closure model in which triple correlations are modelled in terms of quadratic correlations. This local model uses five nondimensional parameters to describe the properties of the flow. We attempt to determine these closure parameters for our simulations and find that the model does produce qualitatively correct behaviour. In addition, we present results concerning the shear rate dependency of the magnetic to kinetic energy ratio. We find both the turbulent stress ratio and the total stress to be strongly dependent on the shear rate (© 2009 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Thermomagnetic and magnetorotational instabilities of thin Keplerian discs in poloidal magnetic fields

Stability of thin hot Keplerian discs is investigated asymptotically in small disc's aspect ratio, ε. The study is carried out in the local approximation for short vertical waves in the disc‐thickness scale. Besides the radial rotation shear and the vertical magnetic field, the background configuration is characterized by a vertically near‐constant temperature profile with a small vertical gradient. The temperature‐gradient term in Ohm's law, which characterizes the thermomagnetic transport is found to be of the order of ε. The effect of the thermomagnetic transport slightly modifies the conventional magnetorotational instability (MRI), while a new thermomagnetic instability (TMI) emerges in regions of the wavenumber space where MRI is absent. Explicit solutions are obtained for a wide range of values of plasma beta, β, and thermomagnetic transport coefficient, λ. In particular, it is shown for λ ≪ 1 that the MRI dominates in weak magnetic fields, β ≫ 1, while the TMI is exhibited in strong magnetic fields, β ∼ 1, also with the growth rate of the order of inverse rotation period (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Transport of angular momentum and chemical species by anisotropic mixing in stellar radiative interiors

Small levels of turbulence can be present in stellar radiative interiors due to, e.g., the instability of rotational shear. In this paper we estimate turbulent transport coefficients for stably stratified rotating stellar radiation zones. Stable stratification induces strong anisotropy with a very small ratio of radial‐to‐horizontal turbulence intensities. Angular momentum is transported mainly due to the correlation between azimuthal and radial turbulent motions induced by the Coriolis force. This non‐diffusive transport known as the Λ‐effect has outward direction in radius and is much more efficient compared to the effect of radial eddy viscosity. Chemical species are transported by small radial diffusion only. This result is confirmed using direct numerical simulations combined with the test‐scalar method. As a consequence of the non‐diffusive transport of angular momentum, the estimated characteristic time of rotational coupling (≲100 Myr) between radiative core and convective envelope in young solar‐type stars is much shorter compared to the time‐scale of Lithium depletion (∼1 Gyr) (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Saturation of the magnetorotational instability at large Elsasser number

The magnetorotational instability is investigated within the shearing box approximation in the large Elsasser number regime. In this regime, which is of fundamental importance to astrophysical accretion disk theory, shear is the dominant source of energy, but the instability itself requires the presence of a weaker vertical magnetic field. Dissipative effects are weaker still but not negligible. The regime explored retains the condition that (viscous and ohmic) dissipative forces do not play a role in the leading order linear instability mechanism. However, they are sufficiently large to permit a nonlinear feedback mechanism whereby the turbulent stresses generated by the MRI act on and modify the local background shear in the angular velocity profile. To date this response has been omitted in shearing box simulations and is captured by a reduced pde model derived here from the global MHD fluid equations using multiscale asymptotic perturbation theory. Results from numerical simulations of the reduced pde model indicate a linear phase of exponential growth followed by a nonlinear adjustment to algebraic growth and decay in the fluctuating quantities. Remarkably, the velocity and magnetic field correlations associated with these algebraic growth and decay laws conspire to achieve saturation of the angular momentum transport. The inclusion of subdominant ohmic dissipation arrests the algebraic growth of the fluctuations on a longer, dissipative time scale. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Magnetic and vertical shear instabilities in accretion discs

The stability properties of magnetized discs rotating with angular velocity Ω = Ω( s ,  z ), dependent on both the radial and the vertical coordinates s and z , are considered. Such a rotation law is adequate for many astrophysical discs (e.g., galactic and protoplanetary discs, as well as accretion discs in binaries). In general, the angular velocity depends on height, even in thin accretion discs. A linear stability analysis is performed in the Boussinesq approximation, and the dispersion relation is obtained for short-wavelength perturbations. Any dependence of Ω on z can destabilize the flow. This concerns primarily small-scale perturbations for which the stabilizing effect of buoyancy is strongly suppressed due to the energy exchange with the surrounding plasma. For a weak magnetic field, instability of discs is mainly associated with vertical shear, whilst for an intermediate magnetic field the magnetic shear instability, first considered by Chandrasekhar and Velikhov, is more efficient. This instability is caused by the radial shear which is typically much stronger than the vertical shear. Therefore the growth time for the magnetic shear instability is much shorter than for the vertical shear instability. A relatively strong magnetic field can suppress both these instabilities. The vertical shear instability could be the source of turbulence in protoplanetary discs, where the conductivity is low.     

Turbulent viscosity and Λ-effect from numerical turbulence models

Homogeneous anisotropic turbulence simulations are used to determine off-diagonal components of the Reynolds stress tensor and its parameterization in terms of turbulent viscosity and Λ-effect. The turbulence is forced in an anisotropic fashion by enhancing the strength of the forcing in the vertical direction. The Coriolis force is included with a rotation axis inclined relative to the vertical direction. The system studied here is significantly simpler than that of turbulent stratified convection which has often been used to study Reynolds stresses. Certain puzzling features of the results for convection, such as sign changes or highly concentrated latitude distributions, are not present in the simpler system considered here. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Dynamical aspects of the venus 4-day circulation

Evidence for the 4-day retrograde zonal circulation of the upper Venus atmosphere is summarized. The ‘moving flame’ phenomenon, convective instability to a mean shear and tidal forcing are discussed as possible dynamical explanations for the 4-day rotation. Tidal forcing seems feasible only if momentum diffusion is molecular in nature. Convective instability to a mean shear, although it can account for the magnitude of the 4-day circulation, must be supplemented by another mechanism, the ‘moving flame’ say, to explain the direction of the zonal motion. However, numerical computations indicate that the ‘moving flame’ by itself can account for both the magnitude and direction of the 4-day rotation. It appears that the stable stratification above altitudes of about 60 km is an essential factor in the ‘moving flame’ mechanism for generating the retrograde atmospheric rotation.     

Kinetic helicity decay in linearly forced turbulence

The decay of kinetic helicity is studied in numerical models of forced turbulence using either an externally imposed forcing function as an inhomogeneous term in the equations or, alternatively, a term linear in the velocity giving rise to a linear instability. The externally imposed forcing function injects energy at the largest scales, giving rise to a turbulent inertial range with nearly constant energy flux while for linearly forced turbulence the spectral energy is maximum near the dissipation wavenumber. Kinetic helicity is injected once a statistically steady state is reached, but it is found to decay on a turbulent time scale regardless of the nature of the forcing and the value of the Reynolds number (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Stability of accretion discs threaded by a strong magnetic field

We study the stability of poloidal magnetic fields anchored in a thin accretion disc. The two-dimensional hydrodynamics in the disc plane is followed by a grid-based numerical simulation including the vertically integrated magnetic forces. The three-dimensional magnetic field outside the disc is calculated in a potential field approximation from the magnetic flux density distribution in the disc. For uniformly rotating discs we confirm numerically the existence of the interchange instability as predicted by Spruit, Stehle & Papaloizou . In agreement with predictions from the shearing sheet model, discs with Keplerian rotation are found to be stabilized by the shear, as long as the contribution of magnetic forces to support against gravity is small. When this support becomes significant, we find a global instability which transports angular momentum outwardly and allows mass to accrete inwardly. The instability takes the form of a m =1 rotating 'crescent', reminiscent of the purely hydrodynamic non-linear instability previously found in pressure-supported discs. A model where the initial surface mass density Σ( r ) and B _z ( r ) decrease with radius as power laws shows transient mass accretion during about six orbital periods, and settles into a state with surface density and field strength decreasing approximately exponentially with radius. We argue that this instability is likely to be the main angular momentum transport mechanism in discs with a poloidal magnetic field sufficiently strong to suppress magnetic turbulence. It may be especially relevant in jet-producing discs.     

Can turbulent plane‐shear flows support large‐scale dynamos?

Turbulent plane‐shear flow is found to show same basic effects of mean‐fieldMHD as rotating turbulence. In particular, the mean electromotive force (EMF) includes highly anisotropic turbulent diffusion and alpha‐effect. Only magnetic diffusion remains for spatially‐uniform turbulence. The question is addressed whether in this case a self‐excitation of a magnetic field by so‐called sher‐current dynamo is possible and the quasilinear theory provides a negative answer. The streamaligned component of the EMF has the sign opposite to that required for dynamo. If, however, the turbulence is not uniform across the flow direction then a dynamo‐active α ‐effect emerges. The critical magnetic Reynolds number for the alpha‐shear dynamo is estimated to be slightly above ten. Possibilities for cross‐checking theoretical predictions with MHD experiments are discussed. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Growth of the MRI in accretion discs – the influence of radiation transport

In this paper, we investigate the influence of radiative transport on the growth of the magnetorotational instability (MRI) in accretion discs. The analysis is performed by the use of analytical and numerical means. We provide a general dispersion relation together with the corresponding eigenfunctions describing the growth rates of small disturbances on a homogeneous background shear flow. The dispersion relation includes compressibility and radiative effects in the flux-limited diffusion approximation. By introducing an effective speed of sound, all the effects of radiation transport can be subsumed into one single parameter. It can be shown that the growth rates of the vertical modes – which are the fastest growing ones – are reduced by radiative transport. For the case of non-vertical modes, the growth rates may instead be enhanced. We quantify the effects of compressibility and radiative diffusion on the growth rates for the gas-pressure dominated case. The analytical discussion is supplemented by numerical simulations, which are also used for a first investigation of the non-linear stage of the MRI in gas-pressure dominated accretion discs with radiation transport included.     

The dual role of shear in large‐scale dynamos

The role of shear in alleviating catastrophic quenching by shedding small‐scale magnetic helicity through fluxes along contours of constant shear is discussed. The level of quenching of the dynamo effect depends on the quenched value of the turbulent magnetic diffusivity. Earlier estimates that might have suffered from the force‐free degeneracy of Beltrami fields are now confirmed for shear flows where this degeneracy is lifted. For a dynamo that is saturated near equipartition field strength those estimates result in a 5‐fold decrease of the magnetic diffusivity as the magnetic Reynolds number based on the wavenumber of the energy‐carrying eddies is increased from 2 to 600. Finally, the role of shear in driving turbulence and large‐scale fields by the magneto‐rotational instability is emphasized. New simulations are presented and the 3π /4 phase shift between poloidal and toroidal fields is confirmed. It is suggested that this phase shift might be a useful diagnostic tool in identifying mean‐field dynamo action in simulations and to distinguish this from other scenarios invoking magnetic buoyancy as a means to explain migration away from the midplane. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Non-linear fluid dynamics of eccentric discs

A new theory of eccentric accretion discs is presented. Starting from the basic fluid-dynamical equations in three dimensions, I derive the fundamental set of one-dimensional equations that describe how the mass, angular momentum and eccentricity vector of a thin disc evolve as a result of internal stresses and external forcing. The analysis is asymptotically exact in the limit of a thin disc, and allows for slowly varying eccentricities of arbitrary magnitude. The theory is worked out in detail for a Maxwellian viscoelastic model of the turbulent stress in an accretion disc. This generalizes the conventional alpha viscosity model to account for the non-zero relaxation time of the turbulence, and is physically motivated by a consideration of the nature of magnetohydrodynamic turbulence. It is confirmed that circular discs are typically viscously unstable to eccentric perturbations, as found by Lyubarskij, Postnov & Prokhorov, if the conventional alpha viscosity model is adopted. However, the instability can usually be suppressed by introducing a sufficient relaxation time and/or bulk viscosity. It is then shown that an initially uniformly eccentric disc does not retain its eccentricity as had been suggested by previous analyses. The evolutionary equations should be useful in many applications, including understanding the origin of planetary eccentricities and testing theories of quasi-periodic oscillations in X-ray binaries.     

MRI channel flows and their parasites

Local simulations of the magnetorotational instability (MRI) in accretion discs can exhibit recurrent coherent structures called channel flows. The formation and destruction of these structures may play a role in the development and saturation of MRI-induced turbulence, and consequently help us understand the time-dependent accretion behaviour of certain astrophysical objects. Previous investigations have revealed that channel solutions are attacked by various parasitic modes, foremost of which is an analogue of the Kelvin–Helmholtz instability. We revisit these instabilities and show how they relate to the classical instabilities of plasma physics, the kink and pinch modes. However, we argue that in most cases channels emerge from developed turbulence and are eventually destroyed by turbulent mixing, not by the parasites. The exceptions are clean isolated channels, which appear in systems near criticality or which emerge from low amplitude initial conditions. These structures can achieve large amplitudes and are only then destroyed, giving rise to eruptive behaviour.