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 HMRI 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.

The magnetorotational instability (MRI) is thought to be a powerful source of turbulence in Keplerian accretion disks. Motivated by recent laboratory experiments, we study the MRI driven by an azimuthal magnetic field in an electrically conducting fluid sheared between two concentric rotating cylinders. By adjusting the rotation rates of the cylinders, we approximate angular velocity profiles $\omega \propto r^{q}$. We perform direct numerical simulations of a steep profile close to the Rayleigh line $q \gtrsim -2 $ and a quasi-Keplerian profile $q \approx -3/2$ and cover wide ranges of Reynolds ($Re\le 4\cdot10^4$) and magnetic Prandtl numbers ($0\le Pm \le 1$). In the quasi-Keplerian case, the onset of instability depends on the magnetic Reynolds number, with $Rm_c \approx 50$, and angular momentum transport scales as $\sqrt{Pm} Re^2$ in the turbulent regime. The ratio of Maxwell to Reynolds stresses is set by $Rm$. At the onset of instability both stresses have similar magnitude, whereas the Reynolds stress vanishes or becomes even negative as $Rm$ increases. For the profile close to the Rayleigh line, the instability shares these properties as long as $Pm\gtrsim0.1$, but exhibits a markedly different character if $Pm\rightarrow 0$, where the onset of instability is governed by the Reynolds number, with $Re_c \approx 1250$, transport is via Reynolds stresses and scales as $Re^2$. At intermediate $Pm=0.01$ we observe a continuous transition from one regime to the other, with a crossover at $Rm=\mathcal{O}(100)$. Our results give a comprehensive picture of angular momentum transport of the MRI with an imposed azimuthal field.

We consider axially periodic Taylor-Couette geometry with insulating boundary conditions. The imposed basic states are so-called Chandrasekhar states, where the azimuthal flow $U_\phi$ and magnetic field $B_\phi$ have the same radial profiles. Mainly three particular profiles are considered: the Rayleigh limit, quasi-Keplerian, and solid-body rotation. In each case we begin by computing linear instability curves and their dependence on the magnetic Prandtl number Pm. For the azimuthal wavenumber m=1 modes, the instability curves always scale with the Reynolds number and the Hartmann number. For sufficiently small Pm these modes therefore only become unstable for magnetic Mach numbers less than unity, and are thus not relevant for most astrophysical applications. However, modes with m>10 can behave very differently. For sufficiently flat profiles, they scale with the magnetic Reynolds number and the Lundquist number, thereby allowing instability also for the large magnetic Mach numbers of astrophysical objects. We further compute fully nonlinear, three-dimensional equilibration of these instabilities, and investigate how the energy is distributed among the azimuthal (m) and axial (k) wavenumbers. In comparison spectra become steeper for large m, reflecting the smoothing action of shear. On the other hand kinetic and magnetic energy spectra exhibit similar behavior: if several azimuthal modes are already linearly unstable they are relatively flat, but for the rigidly rotating case where m=1 is the only unstable mode they are so steep that neither Kolmogorov nor Iroshnikov-Kraichnan spectra fit the results. The total magnetic energy exceeds the kinetic energy only for large magnetic Reynolds numbers Rm>100.