No Arabic abstract
(abriged) MRI turbulence is a leading mechanism for the generation of an efficient turbulent transport of angular momentum in an accretion disk through a turbulent viscosity effect. It is believed that the same process could also transport large-scale magnetic fields in disks, reshaping the magnetic structures in these objects. This process, known as turbulent resistivity, has been suggested and used in several accretion-ejection models and simulations to produce jets. Still, the efficiency of MRI-driven turbulence to transport large-scale magnetic fields is largely unknown. We investigate this problem both analytically and numerically. We introduce a linear calculation of the MRI in the presence of a spatially inhomogeneous mean magnetic field. We show that, in this configuration, MRI modes lead to an efficient magnetic field transport, on the order of the angular momentum transport. We next use fully non linear simulations of MRI turbulence to compute the turbulent resistivity in several magnetic configurations. We find that the turbulent resistivity is on the order of the turbulent viscosity in all our simulations, although somewhat lower. The turbulent resistivity tensor is found to be highly anisotropic with a diffusion coefficient 3 times greater in the radial direction than in the vertical direction. These results support the possibility of driving jets from turbulent disks; the resulting jets may not be steady.
We investigate sustenance and dependence on magnetic Prandtl number (${rm Pm}$) for magnetorotational instability (MRI)-driven turbulence in astrophysical Keplerian disks with zero net magnetic flux using standard shearing box simulations. We focus on the turbulence dynamics in Fourier space, capturing specific/noncanonical anisotropy of nonlinear processes due to disk flow shear. This is a new type of nonlinear redistribution of modes over wavevector orientations in Fourier space -- the nonlinear transverse cascade -- which is generic to shear flows and fundamentally different from usual direct/inverse cascade. The zero flux MRI has no exponentially growing modes, so its growth is transient, or nonmodal. Turbulence self-sustenance is governed by constructive cooperation of the transient growth of MRI and the nonlinear transverse cascade. This cooperation takes place at small wavenumbers (on the flow size scales) referred to as the vital area in Fourier space. The direct cascade transfers mode energy from the vital area to larger wavenumbers. At large ${rm Pm}$, the transverse cascade prevails over the direct one, keeping most of modes energy contained in small wavenumbers. With decreasing ${rm Pm}$, however, the action of the transverse cascade weakens and can no longer oppose the action of direct cascade which more efficiently transfers energy to higher wavenumbers, leading to increased resistive dissipation. This undermines the sustenance scheme, resulting in the turbulence decay. Thus, the decay of zero net flux MRI-turbulence with decreasing ${rm Pm}$ is attributed to topological rearrangement of the nonlinear processes when the direct cascade begins to prevail over the transverse cascade.
We studied dynamical balances in magnetorotational instability (MRI) turbulence with a net vertical field in the shearing box model of disks. Analyzing the turbulence dynamics in Fourier (${bf k}$-)space, we identified three types of active modes that define turbulence characteristics. These modes have lengths similar to the box size, i.e., lie in the small wavenumber region in Fourier space labeled the vital area and are: (i) the channel mode - uniform in the disk plane with the smallest vertical wavenumber,(ii) the zonal flow mode - azimuthally and vertically uniform with the smallest radial wavenumber and (iii) the rest modes. The rest modes comprise those harmonics in the vital area whose energies reach more than $50 %$ of the maximum spectral energy. The rest modes individually are not so significant compared to the channel and zonal flow modes, however, the combined action of their multitude is dominant over these two modes. These three mode types are governed by interplay of the linear and nonlinear processes, leading to their interdependent dynamics. The linear processes consist in disk flow nonmodality-modified classical MRI with a net vertical field. The main nonlinear process is transfer of modes over wavevector angles in Fourier space - the transverse cascade. The channel mode exhibits episodic bursts supplied by linear MRI growth, while the nonlinear processes mostly oppose this, draining the channel energy and redistributing it to the rest modes. As for the zonal flow, it does not have a linear source and is fed by nonlinear interactions of the rest modes.
Magnetorotational turbulence draws its energy from gravity and ultimately releases it via dissipation. However, the quantitative details of this energy flow have not been assessed for global disk models. In this work we examine the energetics of a well-resolved, three-dimensional, global magnetohydrodynamic accretion disk simulation by evaluating statistically-averaged mean-field equations for magnetic, kinetic, and internal energy using simulation data. The results reveal that turbulent magnetic (kinetic) energy is primarily injected by the correlation between Maxwell (Reynolds) stresses and shear in the (almost Keplerian) mean flow, and removed by dissipation. This finding differs from previous work using local (shearing-box) models, which indicated that turbulent kinetic energy was primarily sourced from the magnetic energy reservoir. Lorentz forces provide the bridge between the magnetic and kinetic energy reservoirs, converting ~ 1/5 of the total turbulent magnetic power input into turbulent kinetic energy. The turbulent energies (both magnetic and kinetic) are mainly driven by terms associated with the turbulent fields, with only a minor influence from mean magnetic fields. The interaction between mean and turbulent fields is most evident in the induction equation, with the mean radial magnetic field being strongly influenced by the turbulent electromotive force (EMF). During the quasi-steady turbulent state roughly 2/3 of the Poynting flux travels into the corona, with the remainder transporting magnetic energy in the radial direction. In contrast to previous studies, the stress-related part of the Poynting flux is found to dominate, which may have important implications for reflection models of Seyfert galaxy coronae that typically invoke a picture of buoyant rising of magnetic flux tubes via advection.
Accretion disks are likely threaded by external vertical magnetic flux, which enhances the level of turbulence via the magnetorotational instability (MRI). Using shearing-box simulations, we find that such external magnetic flux also strongly enhances the amplitude of banded radial density variations known as zonal flows. Moreover, we report that vertical magnetic flux is strongly concentrated toward low-density regions of the zonal flow. Mean vertical magnetic field can be more than doubled in low-density regions, and reduced to nearly zero in high density regions in some cases. In ideal MHD, the scale on which magnetic flux concentrates can reach a few disk scale heights. In the non-ideal MHD regime with strong ambipolar diffusion, magnetic flux is concentrated into thin axisymmetric shells at some enhanced level, whose size is typically less than half a scale height. We show that magnetic flux concentration is closely related to the fact that the magnetic diffusivity of the MRI turbulence is anisotropic. In addition to a conventional Ohmic-like turbulent resistivity, we find that there is a correlation between the vertical velocity and horizontal magnetic field fluctuations that produces a mean electric field that acts to anti-diffuse the vertical magnetic flux. The anisotropic turbulent diffusivity has analogies to the Hall effect, and may have important implications for magnetic flux transport in accretion disks. The physical origin of magnetic flux concentration may be related to the development of channel flows followed by magnetic reconnection, which acts to decrease the mass-to-flux ratio in localized regions. The association of enhanced zonal flows with magnetic flux concentration may lead to global pressure bumps in protoplanetary disks that helps trap dust particles and facilitates planet formation.
Global three dimensional magnetohydrodynamic (MHD) simulations of turbulent accretion disks are presented which start from fully equilibrium initial conditions in which the magnetic forces are accounted for and the induction equation is satisfied. The local linear theory of the magnetorotational instability (MRI) is used as a predictor of the growth of magnetic field perturbations in the global simulations. The linear growth estimates and global simulations diverge when non-linear motions - perhaps triggered by the onset of turbulence - upset the velocity perturbations used to excite the MRI. The saturated state is found to be independent of the initially excited MRI mode, showing that once the disk has expelled the initially net flux field and settled into quasi-periodic oscillations in the toroidal magnetic flux, the dynamo cycle regulates the global saturation stress level. Furthermore, time-averaged measures of converged turbulence, such as the ratio of magnetic energies, are found to be in agreement with previous works. In particular, the globally averaged stress normalized to the gas pressure, <alpha_{rm P}> = 0.034, with notably higher values achieved for simulations with higher azimuthal resolution. Supplementary tests are performed using different numerical algorithms and resolutions. Convergence with resolution during the initial linear MRI growth phase is found for 23-35 cells per scaleheight (in the vertical direction).