Helical Magnetorotational Instability Research Articles

Axisymmetric helical magnetorotational instability in a Taylor-Couette flow with sufficiently large positive shear

Axisymmetric helical magnetorotational instability in a Taylor-Couette flow with sufficiently large positive shear

Destabilization of super-rotating Taylor–Couette flows by current-free helical magnetic fields

In an earlier paper we showed that the combination of azimuthal magnetic fields and super-rotation in Taylor–Couette flows of conducting fluids can be unstable against non-axisymmetric perturbations if the magnetic Prandtl number of the fluid is $\textrm {Pm}\neq 1$ . Here we demonstrate that the addition of a weak axial field component allows axisymmetric perturbation patterns for $\textrm {Pm}$ of order unity depending on the boundary conditions. The axisymmetric modes only occur for magnetic Mach numbers (of the azimuthal field) of order unity, while higher values are necessary for the non-axisymmetric modes. The typical growth time of the instability and the characteristic time scale of the axial migration of the axisymmetric mode are long compared with the rotation period, but short compared with the magnetic diffusion time. The modes travel in the positive or negative $z$ direction along the rotation axis depending on the sign of $B_\phi B_z$ . We also demonstrate that the azimuthal components of flow and field perturbations travel in phase if $|B_\phi |\gg |B_z|$ , independent of the form of the rotation law. Within a short-wave approximation for thin gaps it is also shown (in an appendix) that for ideal fluids the considered helical magnetorotational instability only exists for rotation laws with negative shear.

Two types of axisymmetric helical magnetorotational instability in rotating flows with positive shear

We reveal and investigate a new type of linear axisymmetric helical magnetorotational instability which is capable of destabilizing viscous and resistive rotational flows with radially increasing angular velocity, or positive shear. This instability is double-diffusive by nature and is different from the more familiar helical magnetorotational instability, operating at positive shear above the Liu limit, in that it works instead for a wide range of the positive shear when ${\rm (i)}$ a combination of axial/poloidal and azimuthal/toroidal magnetic fields is applied and ${\rm (ii)}$ the magnetic Prandtl number is not too close to unity. We study this instability first with radially local WKB analysis and then confirm its existence using a global stability analysis of the magnetized flow between two rotating cylinders with conducting or insulating boundaries. From an experimental point of view, we also demonstrate the presence of the new instability in a magnetized viscous and resistive Taylor-Couette flow with positive shear for such values of the flow parameters, which can be realized in upcoming experiments at the DRESDYN facility. Finally, this instability might have implications for the dynamics of the equatorial parts of the solar tachocline and dynamo action there, since the above two necessary conditions for the instability to take place are satisfied in this region. Our global stability calculations for the tachocline-like configuration, representing a thin rotating cylindrical layer with the appropriate boundary conditions -- conducting inner and insulating outer cylinders -- and the values of the flow parameters, indicate that it can indeed arise in this case with a characteristic growth time comparable to the solar cycle period.

Local study of helical magnetorotational instability in viscous Keplerian disks

In this paper, regarding the recent detection of significant azimuthal magnetic field in some accretion disks such as protostellar (Donati et al. in Nature 438:466, 2005), the multi-fluid model has been employed to analysis the stability of Keplerian rotational viscous dusty plasma system in a current-free helical magnetic field structure. Using the fluid-Maxwell equations, the general dispersion relation of the excited modes in the system has been obtained by applying the local approximation method in the linear perturbation theory. The typical numerical analysis of the obtained dispersion relation in the high-frequency regime shows that the presence of azimuthal magnetic field component in Keplerian flow has a considerable role in the stability conditions of the system. It also shows that the magnetic field helicity has a stabilization role against the magnetorotational instability (MRI) in the system due to contraction of the unstable wavelength region and decreasing the maximum growth rate of the instability. In this sense, the stabilization role of the viscosity term is more considerable for HMRI (instability in the presence of azimuthal magnetic field component) than the corresponding MRI (instability in the absence of azimuthal magnetic field component). Moreover, considering the discovered azimuthal magnetic field in these systems, the MRI can be arisen in the over-all range of dust grains construction values in contract with traditional MRI. This investigation can greatly contribute to better understanding the physics of some astrophysical phenomena, such as the main source of turbulence and angular momentum transport in protostellar and the other sufficiently ionized astrophysical disks, where the azimuthal magnetic field component in these systems can play a significant role.

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.

Absolute versus convective helical magnetorotational instabilities in Taylor–Couette flows

We study magnetic Taylor–Couette flow in a system having nondimensional radii ri = 1 and ro = 2, and periodic in the axial direction with wavelengths . The rotation ratio of the inner and outer cylinders is adjusted to be slightly in the Rayleigh-stable regime, where magnetic fields are required to destabilize the flow, in this case triggering the axisymmetric helical magnetorotational instability (HMRI). Two choices of imposed magnetic field are considered, both having the same azimuthal component , but differing axial components. The first choice has Bz = 0.1, and yields the familiar HMRI, consisting of unidirectionally traveling waves. The second choice has , and yields HMRI waves that travel in opposite directions depending on the sign of Bz. The first configuration corresponds to a convective instability, the second to an absolute instability. The two variants behave very similarly regarding both linear onset as well as nonlinear equilibration.

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.

Nonmodel analysis of helical and azimuthal magnetorotational instabilities

Helical and azimuthal magnetorotational instabilities operate in rotating magnetized flows with relatively steep negative or extremely steep positive shear. The corresponding lower and upper Liu limits of the shear, which determine the threshold of modal growth of these instabilities, are continuously connected when some axial electrical current is allowed to pass through the rotating fluid. We investigate the nonmodal dynamics of these instabilities arising from the non-normality of shear flow in the local approximation, generalizing the results of the modal approach. It is demonstrated that moderate transient/nonmodal amplification of both types of magnetorotational instability occurs within the Liu limits, where the system is stable according to modal analysis. We show that for the helical magnetorotational instability this magnetohydrodynamic behavior is closely connected with the nonmodal growth of the underlying purely hydrodynamic problem.

Linking dissipation-induced instabilities with nonmodal growth: The case of helical magnetorotational instability.

The helical magnetorotational instability is known to work for resistive rotational flows with comparably steep negative or extremely steep positive shear. The corresponding lower and upper Liu limits of the shear are continuously connected when some axial electrical current is allowed to flow through the rotating fluid. Using a local approximation we demonstrate that the magnetohydrodynamic behavior of this dissipation-induced instability is intimately connected with the nonmodal growth and the pseudospectrum of the underlying purely hydrodynamic problem.

Generalised quasilinear approximation of the helical magnetorotational instability

Motivated by recent advances in direct statistical simulation (DSS) of astrophysical phenomena such as out-of-equilibrium jets, we perform a direct numerical simulation (DNS) of the helical magnetorotational instability (HMRI) under the generalised quasilinear approximation (GQL). This approximation generalises the quasilinear approximation (QL) to include the self-consistent interaction of large-scale modes, interpolating between fully nonlinear DNS and QL DNS whilst still remaining formally linear in the small scales. In this paper we address whether GQL can more accurately describe low-order statistics of axisymmetric HMRI when compared with QL by performing DNS under various degrees of GQL approximation. We utilise various diagnostics, such as energy spectra in addition to first and second cumulants, for calculations performed for a range of Reynolds and Hartmann numbers (describing rotation and imposed magnetic field strength respectively). We find that GQL performs significantly better than QL in describing the statistics of the HMRI even when relatively few large-scale modes are kept in the formalism. We conclude that DSS based on GQL (GCE2) will be significantly more accurate than that based on QL (CE2).

Nonaxisymmetric linear instability of cylindrical magnetohydrodynamic Taylor-Couette flow.

We consider the nonaxisymmetric modes of instability present in Taylor-Couette flow under the application of helical magnetic fields, mainly for magnetic Prandtl numbers close to the inductionless limit, and conduct a full examination of marginal stability in the resulting parameter space. We allow for the azimuthal magnetic field to be generated by a combination of currents in the inner cylinder and fluid itself and introduce a parameter governing the relation between the strength of these currents. A set of governing eigenvalue equations for the nonaxisymmetric modes of instability are derived and solved by spectral collocation with Chebyshev polynomials over the relevant parameter space, with the resulting instabilities examined in detail. We find that by altering the azimuthal magnetic field profiles the azimuthal magnetorotational instability, nonaxisymmetric helical magnetorotational instability, and Tayler instability yield interesting dynamics, such as different preferred mode types and modes with azimuthal wave number m>1. Finally, a comparison is given to the recent WKB analysis performed by Kirillov etal. [Kirillov, Stefani, and Fukumoto, J. Fluid Mech. 760, 591 (2014)JFLSA70022-112010.1017/jfm.2014.614] and its validity in the linear regime.

Metamorphosis of helical magnetorotational instability in the presence of axial electric current.

This paper presents numerical linear stability analysis of a cylindrical Taylor-Couette flow of liquid metal carrying axial electric current in a generally helical external magnetic field. Axially symmetric disturbances are considered in the inductionless approximation corresponding to zero magnetic Prandtl number. Axial symmetry allows us to reveal an entirely new electromagnetic instability. First, we show that the electric current passing through the liquid can extend the range of helical magnetorotational instability (HMRI) indefinitely by transforming it into a purely electromagnetic instability. Two different electromagnetic instability mechanisms are identified. The first is an internal pinch-type instability, which is due to the interaction of the electric current with its own magnetic field. Axisymmetric mode of this instability requires a free-space component of the azimuthal magnetic field. When the azimuthal component of the magnetic field is purely rotational and the axial component is nonzero, a new kind of electromagnetic instability emerges. The latter, driven by the interaction of electric current with a weak collinear magnetic field in a quiescent fluid, gives rise to a steady meridional circulation coupled with azimuthal rotation.

Extending the Range of the Inductionless Magnetorotational Instability

The magnetorotational instability (MRI) can destabilize hydrodynamically stable rotational flows, thereby allowing angular momentum transport in accretion disks. A notorious problem for the MRI is its questionable applicability in regions with low magnetic Reynolds number. Using the WKB method, we extend the range of applicability of the MRI by showing that the inductionless versions of the MRI, such as the helical MRI and the azimuthal MRI, can easily destabilize Keplerian profiles ∝r(-3/2) if the radial profile of the azimuthal magnetic field is only slightly modified from the current-free profile ∝r(-1). This way we further show how the formerly known lower Liu limit of the critical Rossby number Ro≈-0.828 connects naturally with the upper Liu limit Ro≈+4.828.

Standard and Helical Magnetorotational Instability

The magnetorotational instability (MRI) triggers turbulence and enables outward transport of angular momentum in hydrodynamically stable rotating shear flows, e.g., in accretion disks. What laws of differential rotation are susceptible to the destabilization by axial, azimuthal, or helical magnetic field? The answer to this question, which is vital for astrophysical and experimental applications, inevitably leads to the study of spectral and geometrical singularities on the instability threshold. The singularities provide a connection between seemingly discontinuous stability criteria and thus explain several paradoxes in the theory of MRI that were poorly understood since the 1950s.

Inviscid helical magnetorotational instability in cylindrical Taylor-Couette flow

This paper presents the analysis of axisymmetric helical magnetorotational instability (HMRI) in the inviscid limit, which is relevant for astrophysical conditions. The inductionless approximation defined by zero magnetic Prandtl number is adopted to distinguish the HMRI from the standard MRI in the cylindrical Taylor-Couette flow subject to a helical magnetic field. Using a Chebyshev collocation method convective and absolute instability thresholds are computed in terms of the Elsasser number for a fixed ratio of inner and outer radii λ = 2 and various ratios of rotation rates and helicities of the magnetic field. It is found that the extension of self-sustained HMRI modes beyond the Rayleigh limit does not reach the astrophysically relevant Keplerian rotation profile not only in the narrow- but also in the finite-gap approximation. The Keppler limit can be attained only by the convective HMRI mode provided that the boundaries are perfectly conducting. However, this mode requires not only a permanent external excitation to be observable but also has a long axial wavelength, which is not compatible with limited thickness of astrophysical accretion disks.

Helical magnetorotational instability in a Taylor-Couette flow with strongly reduced Ekman pumping

The magnetorotational instability (MRI) is thought to play a key role in the formation of stars and black holes by sustaining the turbulence in hydrodynamically stable Keplerian accretion disks. In previous experiments the MRI was observed in a liquid metal Taylor-Couette flow at moderate Reynolds numbers by applying a helical magnetic field. The observation of this helical MRI (HMRI) was interfered with a significant Ekman pumping driven by solid end caps that confined the instability only to a part of the Taylor-Couette cell. This paper describes the observation of the HMRI in an improved Taylor-Couette setup with the Ekman pumping significantly reduced by using split end caps. The HMRI, which now spreads over the whole height of the cell, appears much sharper and in better agreement with numerical predictions. By analyzing various parameter dependencies we conclude that the observed HMRI represents a self-sustained global instability rather than a noise-sustained convective one.

Absolute versus convective helical magnetorotational instability in a Taylor-Couette flow

We analyze numerically the magnetorotational instability of a Taylor-Couette flow in a helical magnetic field [helical magnetorotational instability (HMRI)] using the inductionless approximation defined by a zero magnetic Prandtl number (Pr_{m}=0) . The Chebyshev collocation method is used to calculate the eigenvalue spectrum for small-amplitude perturbations. First, we carry out a detailed conventional linear stability analysis with respect to perturbations in the form of Fourier modes that corresponds to the convective instability which is not in general self-sustained. The helical magnetic field is found to extend the instability to a relatively narrow range beyond its purely hydrodynamic limit defined by the Rayleigh line. There is not only a lower critical threshold at which HMRI appears but also an upper one at which it disappears again. The latter distinguishes the HMRI from a magnetically modified Taylor vortex flow. Second, we find an absolute instability threshold as well. In the hydrodynamically unstable regime before the Rayleigh line, the threshold of absolute instability is just slightly above the convective one although the critical wavelength of the former is noticeably shorter than that of the latter. Beyond the Rayleigh line the lower threshold of absolute instability rises significantly above the corresponding convective one while the upper one descends significantly below its convective counterpart. As a result, the extension of the absolute HMRI beyond the Rayleigh line is considerably shorter than that of the convective instability. The absolute HMRI is supposed to be self-sustained and, thus, experimentally observable without any external excitation in a system of sufficiently large axial extension.

NOISE-SUSTAINED CONVECTIVE INSTABILITY IN A MAGNETIZED TAYLOR-COUETTE FLOW

The helical magnetorotational instability of the magnetized Taylor-Couette flow is studied numerically in a finite cylinder. A distant upstream insulating boundary is shown to stabilize the convective instability entirely while reducing the growth rate of the absolute instability. The reduction is less severe with larger height. After modeling the boundary conditions properly, the wave patterns observed in the experiment turn out to be a noise-sustained convective instability. After the source of the noise resulted from unstable Ekman and Stewartson layers is switched off, a slowly-decaying inertial oscillation is observed in the simulation. We reach the conclusion that the experiments completed to date have not yet reached the regime of absolute instability.

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)

Comment on “Helical magnetorotational instability in magnetized Taylor-Couette flow”

Liu [Phys. Rev. E 74, 056302 (2006)] have presented a WKB analysis of the helical magnetorotational instability (HMRI), and claim that it does not exist for Keplerian rotation profiles. We show that, if radial boundary conditions are included, the HMRI can exist even for rotation profiles as flat as Keplerian, provided only that at least one of the boundaries is sufficiently conducting.