Standard Magnetorotational Instability Research Articles

Taylor-Couette flow for astrophysical purposes.

A concise review is given of astrophysically motivated experimental and theoretical research on Taylor-Couette flow. The flows of interest rotate differentially with the inner cylinder faster than the outer, but are linearly stable against Rayleigh's inviscid centrifugal instability. At shear Reynolds numbers as large as [Formula: see text], hydrodynamic flows of this type (quasi-Keplerian) appear to be nonlinearly stable: no turbulence is seen that cannot be attributed to interaction with the axial boundaries, rather than the radial shear itself. Direct numerical simulations agree, although they cannot yet reach such high Reynolds numbers. This result indicates that accretion-disc turbulence is not purely hydrodynamic in origin, at least insofar as it is driven by radial shear. Theory, however, predicts linear magnetohydrodynamic (MHD) instabilities in astrophysical discs: in particular, the standard magnetorotational instability (SMRI). MHD Taylor-Couette experiments aimed at SMRI are challenged by the low magnetic Prandtl numbers of liquid metals. High fluid Reynolds numbers and careful control of the axial boundaries are required. The quest for laboratory SMRI has been rewarded with the discovery of some interesting inductionless cousins of SMRI, and with the recently reported success in demonstrating SMRI itself using conducting axial boundaries. Some outstanding questions and near-future prospects are discussed, especially in connection with astrophysics. This article is part of the theme issue 'Taylor-Couette and related flows on the centennial of Taylor's seminal Philosophical Transactions paper (Part 2)'.

Observation of Axisymmetric Standard Magnetorotational Instability in the Laboratory.

We report the first direct evidence for the axisymmetric standard magnetorotational instability (SMRI) from a combined experimental and numerical study of a magnetized liquid-metal shear flow in a Taylor-Couette cell with independently rotating and electrically conducting end caps. When a uniform vertical magnetic field B_{i} is applied along the rotation axis, the measured radial magnetic field B_{r} on the inner cylinder increases linearly with a small magnetic Reynolds number Rm due to the magnetization of the residue Ekman circulation. Onset of the axisymmetric SMRI is identified from the nonlinear increase of B_{r} beyond a critical Rm in both experiments and nonlinear numerical simulations. The axisymmetric SMRI exists only at sufficiently large Rm and intermediate B_{i}, a feature consistent with theoretical predictions. Our simulations further show that the axisymmetric SMRI causes the velocity and magnetic fields to contribute an outward flux of axial angular momentum in the bulk region, just as it should in accretion disks.

A Nonlocal Magneto-curvature Instability in a Differentially Rotating Disk

A global mode is shown to be unstable to nonaxisymmetric perturbations in a differentially rotating Keplerian disk containing either vertical or azimuthal magnetic fields. In an unstratified cylindrical disk model, using both global eigenvalue stability analysis and linear global initial-value simulations, it is demonstrated that this instability dominates at strong magnetic fields where local standard magnetorotational instability (MRI) becomes stable. Unlike the standard MRI mode, which is concentrated in the high flow shear region, these distinct global modes (with low azimuthal mode numbers) are extended in the global domain and are Alfvén-continuum-driven unstable modes. As its mode structure and relative dominance over MRI are inherently determined by the global spatial curvature as well as the flow shear in the presence of a magnetic field, we call it the magneto-curvature (magneto-spatial-curvature) instability. Consistent with the linear analysis, as the field strength is increased in the nonlinear simulations, a transition from MRI-driven turbulence to a state dominated by global nonaxisymmetric modes is obtained. This global instability could therefore be a source of nonlinear transport in accretion disks at a higher magnetic field than predicted by local models.

Identification of a non-axisymmetric mode in laboratory experiments searching for standard magnetorotational instability

The standard magnetorotational instability (SMRI) is a promising mechanism for turbulence and rapid accretion in astrophysical disks. It is a magnetohydrodynamic (MHD) instability that destabilizes otherwise hydrodynamically stable disk flow. Due to its microscopic nature at astronomical distances and stringent requirements in laboratory experiments, SMRI has remained unconfirmed since its proposal, despite its astrophysical importance. Here we report a nonaxisymmetric MHD instability in a modified Taylor-Couette experiment. To search for SMRI, a uniform magnetic field is imposed along the rotation axis of a swirling liquid-metal flow. The instability initially grows exponentially, becoming prominent only for sufficient flow shear and moderate magnetic field. These conditions for instability are qualitatively consistent with SMRI, but at magnetic Reynolds numbers below the predictions of linear analyses with periodic axial boundaries. Three-dimensional numerical simulations, however, reproduce the observed instability, indicating that it grows linearly from the primary axisymmetric flow modified by the applied magnetic field.

From helical to standard magnetorotational instability: Predictions for upcoming liquid sodium experiments

This work is a first preparatory theoretical study for upcoming magnetorotational instability (MRI) experiments with liquid sodium planned within the DRESDYN project. Therefore, we focus on the characteristics of MRI for typical values of the main parameters of the basic magnetized Taylor-Couette (TC) flow achievable in these experiments. In contrast to previous attempts at finding MRI in the laboratory, our results clearly indicate that standard MRI and its helically modified version can be detected for the first time in the DRESDYN-TC device, including the astrophysically important Keplerian rotation, despite the very small magnetic Prandtl number of liquid sodium.

High-speed standard magneto-rotational instability

The large Reynolds number asymptotic approximations of the neutral curve of Taylor–Couette flow subject to an axial uniform magnetic field are analysed. The flow has been extensively studied since the early 1990s as the magneto-rotational instability (MRI) occurring in the flow may explain the origin of the instability observed in some astrophysical objects. Elsewhere, the ideal approximation has been used to study high-speed flows, which sometimes produces paradoxical results. For example, ideal flows must be completely stabilised for a sufficiently strong applied magnetic field. On the other hand, the vanishing magnetic Prandtl number limit of the stability should be purely hydrodynamic, so instability must occur when Rayleigh’s stability condition is violated. Our first discovery is that this apparent contradiction can be resolved by showing the abrupt appearance of the hydrodynamic instability at a certain critical value of the magnetic Prandtl number. This is found using the asymptotically large Reynolds number limit but with a sufficiently long wavelength to retain some diffusive effects. Our second finding concerns the so-called Velikhov–Chandrasekhar paradox, namely the mismatch of the zero external magnetic field limit of the Velikhov–Chandrasekhar stability criterion and Rayleigh’s stability criterion. We show for fully wide-gap cases that the high Reynolds number asymptotic analysis of the MRI naturally yields the simple stability condition that describes smooth transition from Rayleigh to Velikhov–Chandrasekhar stability criteria with increasing Lundquist number.

Experimental confirmation of the standard magnetorotational instability mechanism with a spring-mass analogue

The magnetorotational instability (MRI) has long been considered a plausibly ubiquitous mechanism to destabilize otherwise stable Keplerian flows to support radially outward transport of angular momentum. Such an efficient transport process would allow fast accretion in astrophysical objects such as stars and black holes to release copious kinetic energy that powers many of the most luminous sources in the universe. But the standard MRI under a purely vertical magnetic field has heretofore never been directly measured despite numerous efforts over more than a decade. Here we report an unambiguous laboratory demonstration of the spring-mass analogue to the standard MRI by comparing motion of a spring-tethered ball within different rotating flows. The experiment corroborates the theory: efficient outward angular momentum transport manifests only for cases with a weak spring in quasi-Keperian flow. Our experimental method accomplishes this in a new way, thereby connecting solid and fluid mechanics to plasma astrophysics.

Azimuthal magnetorotational instability with super-rotation

It is demonstrated that the azimuthal magnetorotational instability (AMRI) also works with radially increasing rotation rates contrary to the standard magnetorotational instability for axial fields which requires negative shear. The stability against non-axisymmetric perturbations of a conducting Taylor–Couette flow with positive shear under the influence of a toroidal magnetic field is considered if the background field between the cylinders is current free. For small magnetic Prandtl number $Pm\rightarrow 0$ the curves of neutral stability converge in the (Hartmann number,Reynolds number) plane approximating the stability curve obtained in the inductionless limit $Pm=0$. The numerical solutions for $Pm=0$ indicate the existence of a lower limit of the shear rate. For large $Pm$ the curves scale with the magnetic Reynolds number of the outer cylinder but the flow is always stable for magnetic Prandtl number unity as is typical for double-diffusive instabilities. We are particularly interested to know the minimum Hartmann number for neutral stability. For models with resting or almost resting inner cylinder and with perfectly conducting cylinder material the minimum Hartmann number occurs for a radius ratio of $r_{\text{in}}=0.9$. The corresponding critical Reynolds numbers are smaller than $10^{4}$.

Quasi-two-dimensional nonlinear evolution of helical magnetorotational instability in a magnetized Taylor–Couette flow

Magnetorotational instability (MRI) is one of the fundamental processes in astrophysics, driving angular momentum transport and mass accretion in a wide variety of cosmic objects. Despite much theoretical/numerical and experimental efforts over the last decades, its saturation mechanism and amplitude, which sets the angular momentum transport rate, remains not well understood, especially in the limit of high resistivity, or small magnetic Prandtl numbers typical to interiors (dead zones) of protoplanetary disks, liquid cores of planets and liquid metals in laboratory. Using direct numerical simulations, in this paper we investigate the nonlinear development and saturation properties of the helical magnetorotational instability (HMRI)—a relative of the standard MRI—in a magnetized Taylor–Couette flow at very low magnetic Prandtl number (correspondingly at low magnetic Reynolds number) relevant to liquid metals. For simplicity, the ratio of azimuthal field to axial field is kept fixed. From the linear theory of HMRI, it is known that the Elsasser number, or interaction parameter determines its growth rate and plays a special role in the dynamics. We show that this parameter is also important in the nonlinear problem. By increasing its value, a sudden transition from weakly nonlinear, where the system is slightly above the linear stability threshold, to strongly nonlinear, or turbulent regime occurs. We calculate the azimuthal and axial energy spectra corresponding to these two regimes and show that they differ qualitatively. Remarkably, the nonlinear state remains in all cases nearly axisymmetric suggesting that this HMRI-driven turbulence is quasi two-dimensional in nature. Although the contribution of non-axisymmetric modes increases moderately with the Elsasser number, their total energy remains much smaller than that of the axisymmetric ones.

Pressure-anisotropy-induced nonlinearities in the kinetic magnetorotational instability.

In collisionless and weakly collisional plasmas, such as hot accretion flows onto compact objects, the magnetorotational instability (MRI) can differ significantly from the standard (collisional) MRI. In particular, pressure anisotropy with respect to the local magnetic-field direction can both change the linear MRI dispersion relation and cause nonlinear modifications to the mode structure and growth rate, even when the field and flow perturbations are very small. This work studies these pressure-anisotropy-induced nonlinearities in the weakly nonlinear, high-ion-beta regime, before the MRI saturates into strong turbulence. Our goal is to better understand how the saturation of the MRI in a low-collisionality plasma might differ from that in the collisional regime. We focus on two key effects: (i) the direct impact of self-induced pressure-anisotropy nonlinearities on the evolution of an MRI mode, and (ii) the influence of pressure anisotropy on the 'parasitic instabilities' that are suspected to cause the mode to break up into turbulence. Our main conclusions are: (i) The mirror instability regulates the pressure anisotropy in such a way that the linear MRI in a collisionless plasma is an approximate nonlinear solution once the mode amplitude becomes larger than the background field (just as in magnetohyrodynamics). This implies that differences between the collisionless and collisional MRI become unimportant at large amplitudes. (ii) The break up of large-amplitude MRI modes into turbulence via parasitic instabilities is similar in collisionless and collisional plasmas. Together, these conclusions suggest that the route to magnetorotational turbulence in a collisionless plasma may well be similar to that in a collisional plasma, as suggested by recent kinetic simulations. As a supplement to these findings, we offer guidance for the design of future kinetic simulations of magnetorotational turbulence.

The Weakly Nonlinear Magnetorotational Instability in a Global, Cylindrical Taylor–Couette Flow

Abstract We conduct a global, weakly nonlinear analysis of the magnetorotational instability (MRI) in a Taylor–Couette flow. This is a multiscale, perturbative treatment of the nonideal, axisymmetric MRI near threshold, subject to realistic radial boundary conditions and cylindrical geometry. We analyze both the standard MRI, initialized by a constant vertical background magnetic field, and the helical MRI, with an azimuthal background field component. This is the first weakly nonlinear analysis of the MRI in a global Taylor–Couette geometry, as well as the first weakly nonlinear analysis of the helical MRI. We find that the evolution of the amplitude of the standard MRI is described by a real Ginzburg–Landau equation (GLE), whereas the amplitude of the helical MRI takes the form of a complex GLE. This suggests that the saturated state of the helical MRI may itself be unstable on long spatial and temporal scales.

Angular Momentum Transport in Accretion Disk Boundary Layers Around Weakly Magnetized Stars

The standard model for turbulent shear viscosity in accretion disks is based on the assumption that angular momentum transport is opposite to the radial angular frequency gradient of the disk. This implies that the turbulent stress must be negative and thus transport angular momentum inwards, in the boundary layer where the accretion disk meets the surface of a weakly magnetized star. However, this behavior is not supported by numerical simulations of turbulent magnetohydrodynamic (MHD) accretion disks, which show that angular momentum transport driven by the magnetorotational instability (MRI) is inefficient in disk regions where, as expected in boundary layers, the angular frequency increases with radius. Motivated by the need of a deeper understanding of the behavior of an MHD fluid in a differentially rotating background that deviates from a Keplerian profile, we study the dynamics of MHD waves in configurations that are stable to the standard MRI. Employing the shearing-sheet framework, we show that transient amplification of shearing MHD waves can generate magnetic energy without leading to a substantial generation of hydromagnetic stresses. While these results are in agreement with numerical simulations, they emphasize the need to better understand the mechanism for angular momentum transport in the inner disk regions on more solid grounds.

A UNIFYING PICTURE OF HELICAL AND AZIMUTHAL MAGNETOROTATIONAL INSTABILITY, AND THE UNIVERSAL SIGNIFICANCE OF THE LIU LIMIT

The magnetorotational instability (MRI) plays a key role for cosmic structure formation by triggering turbulence in the rotating flows of accretion disks that would be otherwise hydrodynamically stable. In the limit of small magnetic Prandtl number, the helical and the azimuthal version of MRI are known to be governed by a quite different scaling behaviour than the standard MRI with a vertical applied magnetic field. Using the short-wavelength approximation for an incompressible, resistive, and viscous rotating fluid we present a unified description of helical and azimuthal MRI, and we identify the universal character of the Liu limit Ro= -0.8284 for the critical Rossby number. From this universal behaviour we are also lead to the prediction of higher azimuthal wavenumber for rather small ratios of azimuthal to axial applied fields.

ON HYDROMAGNETIC STRESSES IN ACCRETION DISK BOUNDARY LAYERS

Detailed calculations of the physical structure of accretion disk boundary layers, and thus their inferred observational properties, rely on the assumption that angular momentum transport is opposite to the radial angular frequency gradient of the disk. The standard model for turbulent shear viscosity satisfies this assumption by construction. However, this behavior is not supported by numerical simulations of turbulent magnetohydrodynamic (MHD) accretion disks, which show that angular momentum transport driven by the magnetorotational instability (MRI) is inefficient in disk regions where, as expected in boundary layers, the angular frequency increases with radius. In order to shed light into physically viable mechanisms for angular momentum transport in this inner disk region, we examine the generation of hydromagnetic stresses and energy density in differentially rotating backgrounds with angular frequencies that increase outward in the shearing-sheet framework. We isolate the modes that are unrelated to the standard MRI and provide analytic solutions for the long-term evolution of the resulting shearing MHD waves. We show that, although the energy density of these waves can be amplified significantly, their associated stresses oscillate around zero, rendering them an inefficient mechanism to transport significant angular momentum (inward). These findings are consistent with the results obtained in numerical simulations of MHD accretion disk boundary layers and challenge the standard assumption of efficient angular momentum transport in the inner disk regions. This suggests that the detailed structure of turbulent MHD accretion disk boundary layers could differ appreciably from those derived within the standard framework of turbulent shear viscosity.

Dissipative Taylor-Couette flows under the influence of helical magnetic fields

The linear stability of magnetohydrodynamic Taylor-Couette flows in axially unbounded cylinders is considered for magnetic Prandtl number unity. Magnetic background fields varying from purely axial to purely azimuthal are imposed, with a general helical field parametrized by β=B(ϕ)/B(z). We map out the transition from the standard magnetorotational instability (MRI) for β=0 to the nonaxisymmetric azimuthal magnetorotational instability for β→∞. For finite β, positive and negative wave numbers m , corresponding to right and left spirals, are no longer degenerate. For the nonaxisymmetric modes, the most unstable mode spirals in the opposite direction to the background field. The standard (β=0) MRI is axisymmetric for weak fields (including the instability with the lowest Reynolds number) but is nonaxisymmetric for stronger fields. If the azimuthal field is due in part to an axial current flowing through the fluid itself (and not just along the central axis), then it is also unstable to the nonaxisymmetric Tayler instability which is most effective without rotation. For purely toroidal fields the solutions for m=±1 are identical so that in this case no preferred helicity results. For large β the wave number m=-1 is preferred, whereas for β≲1 the mode with m=-2 is most unstable. The most unstable modes always spiral in the same direction as the background field. For background fields with positive and not too large β the kinetic helicity of the fluctuations proves to be negative for all the magnetic instabilities considered.

ON THE RELATION OF STANDARD AND HELICAL MAGNETOROTATIONAL INSTABILITY

The magnetorotational instability (MRI) plays a crucial role for cosmic structure formation by enabling turbulence in Keplerian disks which would be otherwise hydrodynamically stable. With particular focus on MRI experiments with liquid metals, which have small magnetic Prandtl numbers, it has been shown that the helical version of this instability (HMRI) has a scaling behavior that is quite different from that of the standard MRI (SMRI). We discuss the relation of HMRI to SMRI by exploring various parameter dependencies. We identify the mechanism of transfer of instability between modes through a spectral exceptional point that explains both the transition from a stationary instability (SMRI) to an unstable traveling wave (HMRI) and the excitation of HMRI in the inductionless limit. For certain parameter regions, we find new islands of the HMRI.

Accretion and Destruction of Planetesimals in Turbulent Disks

We study the conditions for collisions between planetesimals to be accretional or disruptive in turbulent disks, through analytical arguments based on fluid dynamical simulations and orbital integrations. In turbulent disks, the velocity dispersion of planetesimals is pumped up by random gravitational perturbations from density fluctuations of the disk gas. When the velocity dispersion is larger than the planetesimals' surface escape velocity, collisions between planetesimals do not result in accretion and may even lead to their destruction. In disks with a surface density equal to that of the minimum-mass solar nebula and with nominal magnetorotational instability (MRI) turbulence, we find that accretion proceeds only for planetesimals with sizes above ~300 km at 1 AU and ~1000 km at 5 AU. We find that accretion is facilitated in disks with smaller masses. However, at 5 AU and for nominal turbulence strength, km-sized planetesimals are in a highly erosive regime even for a disk mass as small as a fraction of the mass of Jupiter. The existence of giant planets implies that either turbulence was weaker than calculated by standard MRI models or some mechanism was capable of producing Ceres-mass planetesimals on very short timescales. In any case, our results show that in the presence of turbulence planetesimal accretion is most difficult in massive disks and at large orbital distances.

Helical magnetorotational instability of Taylor‐Couette flows in the Rayleigh limit and for quasi‐Kepler rotation

AbstractThe magnetorotational instability (MRI) of differential rotation under the simultaneous presence of axial and azimuthal components of the (current‐free) magnetic field is considered. For rotation with uniform specific angular momentum the MHD equations for axisymmetric perturbations are solved in a local short‐wave approximation. All the solutions are overstable for Bz · Bϕ ≠ 0 with eigenfrequencies approaching the viscous frequency. For more flat rotation laws the results of the local approximation do not comply with the results of a global calculation of the MHD instability of Taylor‐Couette flows between rotating cylinders. – With Bϕ and Bz of the same order the traveling‐mode solutions are also prefered for flat rotation laws such as the quasi‐Kepler rotation. For magnetic Prandtl number Pm → 0 they scale with the Reynolds number of rotation rather than with the magnetic Reynolds number (as for standard MRI) so that they can easily be realized in MHD laboratory experiments. – Regarding the nonaxisymmetric modes one finds a remarkable influence of the ratio Bϕ/Bz only for the extrema. For Bϕ ≫ Bz and for not too small Pm the nonaxisymmetric modes dominate the traveling axisymmetric modes. For standard MRI with Bz ≫ Bϕ, however, the critical Reynolds numbers of the nonaxisymmetric modes exceed the values for the axisymmetric modes by many orders so that they are never prefered. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)

Global Axisymmetric Magnetorotational Instability with Density Gradients

We examine global incompressible axisymmetric perturbations of a differentially rotating MHD plasma with radial density gradients. It is shown that the standard magnetorotational instability (MRI) criterion drawn from the local dispersion relation is often misleading. If the equilibrium magnetic field is either purely axial or purely toroidal, the problem reduces to finding the global radial eigenvalues of an effective potential. The standard Keplerian profile including the origin is mathematically ill-posed, and thus any solution will depend strongly on the inner boundary. We find a class of unstable modes localized by the form of the rotation and density profiles, with reduced dependence on boundary conditions.

Helical magnetorotational instability in magnetized Taylor-Couette flow

Hollerbach and Rüdiger have reported a new type of magnetorotational instability (MRI) in magnetized Taylor-Couette flow in the presence of combined axial and azimuthal magnetic fields. The salient advantage of this "helical" MRI (HMRI) is that marginal instability occurs at arbitrarily low magnetic Reynolds and Lundquist numbers, suggesting that HMRI might be easier to realize than standard MRI (axial field only), and that it might be relevant to cooler astrophysical disks, especially those around protostars, which may be quite resistive. We confirm previous results for marginal stability and calculate HMRI growth rates. We show that in the resistive limit, HMRI is a weakly destabilized inertial oscillation propagating in a unique direction along the axis. But we report other features of HMRI that make it less attractive for experiments and for resistive astrophysical disks. Large axial currents are required. More fundamentally, instability of highly resistive flow is peculiar to infinitely long or periodic cylinders: finite cylinders with insulating endcaps are shown to be stable in this limit, at least if viscosity is neglected. Also, Keplerian rotation profiles are stable in the resistive limit regardless of axial boundary conditions. Nevertheless, the addition of a toroidal field lowers thresholds for instability even in finite cylinders.