Axisymmetric Magnetorotational Instability in Viscous Accretion Disks

TL;DR

This paper investigates the axisymmetric magnetorotational instability (MRI) in viscous accretion disks through linear analysis and nonlinear simulations, revealing how viscosity influences MRI growth and energy conversion.

Contribution

It provides new insights into the nonlinear behavior of viscous MRI and its dependence on Reynolds and Lundquist numbers, extending understanding to regimes with both viscosity and resistivity.

Findings

Linear growth rate suppressed at low Reynolds number

Nonlinear evolution shows continued growth of two-channel flow

Magnetic energy dominates energy conversion processes

Abstract

Axisymmetric magnetorotational instability (MRI) in viscous accretion disks is investigated by linear analysis and two-dimensional nonlinear simulations. The linear growth of the viscous MRI is characterized by the Reynolds number defined as $R_{\rm MRI} \equiv v_A^2/\nu\Omega $, where $v_A$ is the Alfv{\'e}n velocity, $\nu$ is the kinematic viscosity, and $\Omega$ is the angular velocity of the disk. Although the linear growth rate is suppressed considerably as the Reynolds number decreases, the nonlinear behavior is found to be almost independent of $R_{\rm MRI}$. At the nonlinear evolutionary stage, a two-channel flow continues growing and the Maxwell stress increases until the end of calculations even though the Reynolds number is much smaller than unity. A large portion of the injected energy to the system is converted to the magnetic energy. The gain rate of the thermal energy, on the other hand, is found to be much larger than the viscous heating rate. Nonlinear behavior of the MRI in the viscous regime and its difference from that in the highly resistive regime can be explained schematically by using the characteristics of the linear dispersion relation. Applying our results to the case with both the viscosity and resistivity, it is anticipated that the critical value of the Lundquist number $S_{\rm MRI} \equiv v_A^2/\eta\Omega$ for active turbulence depends on the magnetic Prandtl number $S_{{\rm MRI},c} \propto Pm^{1/2}$ in the regime of $Pm \gg 1$ and remains constant when $Pm \ll 1$, where $Pm \equiv S_{\rm MRI}/R_{\rm MRI} = \nu/\eta$ and $\eta$ is the magnetic diffusivity.