Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userResolution Dependence of Magnetorotational Turbulence in the Isothermal Stratified Shearing Box
94
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Magnetohydrodynamic (MHD) turbulence driven by the magnetorotational instability can provide diffusive transport of angular momentum in astrophysical disks, and a widely studied computational model for this process is the ideal, stratified, isothermal shearing box. Here we report results of a convergence study of such boxes up to a resolution of $N = 256$ zones per scale height, performed on blue waters at NCSA with ramses-gpu. We find that the time and vertically integrated dimensionless shear stress $overline{alpha} sim N^{-1/3}$, i.e. the shear stress is resolution dependent. We also find that the magnetic field correlation length decreases with resolution, $lambda sim N^{-1/2}$. This variation is strongest at the disk midplane. We show that our measurements of $alpha$ are consistent with earlier studies. We discuss possible reasons for the lack of convergence.
rate research
Read More
We develop a framework for magnetohydrodynamical (MHD) simulations in a local cylindrical shearing box by extending the formulation of the Cartesian shearing box. We construct shearing-periodic conditions at the radial boundaries of a simulation box from the conservation relations of the basic MHD equations, taking into account the explicit radial dependence of physical quantities. We demonstrate quasi-steady mass accretion, which cannot be handled by the standard Cartesian shearing box model, with an ideal MHD simulation in a vertically unstratified cylindrical shearing box up to 200 rotations. In this demonstrative run we set up (i) net vertical magnetic flux, (ii) a locally isothermal equation of state, and (iii) a sub-Keplerian equilibrium rotation, whereas the sound velocity and the initial Alfven velocity have the same radial dependence as that of the Keplerian velocity. Inward mass accretion is induced to balance with the outward angular momentum flux of the MHD turbulence triggered by the magnetorotational instability in a self-consistent manner. We discuss detailed physical properties of the saturated magnetic field, in comparison to the results of a Cartesian shearing box simulation.
The launching process of a magnetically driven outflow from an accretion disk is investigated in a local, shearing box model which allows a study of the feedback between accretion and angular momentum loss. The mass-flux instability found in previous linear analyses of this problem is recovered in a series of 2D (axisymmetric) simulations in the MRI-stable (high magnetic field strength) regime. At low field strengths that are still sufficient to suppress MRI, the instability develops on a short radial length scale and saturates at a modest amplitude. At high field strengths, a long-wavelength clump instability of large amplitude is observed, with growth times of a few orbits. As speculated before, the unstable connection between disk and outflow may be relevant for the time dependence observed in jet-producing disks. The success of the simulations is due in a large part to the implementation of an effective wave-transmitting upper boundary condition.
The phenomenological Disc Instability Model has been successful in reproducing the observed light curves of dwarf nova outbursts by invoking an enhanced Shakura-Sunyaev $alpha$ parameter $sim0.1-0.2$ in outburst compared to a low value $sim0.01$ in quiescence. Recent thermodynamically consistent simulations of magnetorotational (MRI) turbulence with appropriate opacities and equation of state for dwarf nova accretion discs have found that thermal convection enhances $alpha$ in discs in outburst, but only near the hydrogen ionization transition. At higher temperatures, convection no longer exists and $alpha$ returns to the low value comparable to that in quiescence. In order to check whether this enhancement near the hydrogen ionization transition is sufficient to reproduce observed light curves, we incorporate this MRI-based variation in $alpha$ into the Disc Instability Model, as well as simulation-based models of turbulent dissipation and convective transport. These MRI-based models can successfully reproduce observed outburst and quiescence durations, as well as outburst amplitudes, albeit with different parameters from the standard Disc Instability Models. The MRI-based model lightcurves exhibit reflares in the decay from outburst, which are not generally observed in dwarf novae. However, we highlight the problematic aspects of the quiescence physics in the Disc Instability Model and MRI simulations that are responsible for this behavior.
We present results from the first 3D kinetic numerical simulation of magnetorotational turbulence and dynamo, using the local shearing-box model of a collisionless accretion disc. The kinetic magnetorotational instability grows from a subthermal magnetic field having zero net flux over the computational domain to generate self-sustained turbulence and outward angular-momentum transport. Significant Maxwell and Reynolds stresses are accompanied by comparable viscous stresses produced by field-aligned ion pressure anisotropy, which is regulated primarily by the mirror and ion-cyclotron instabilities through particle trapping and pitch-angle scattering. The latter endow the plasma with an effective viscosity that is biased with respect to the magnetic-field direction and spatio-temporally variable. Energy spectra suggest an Alfven-wave cascade at large scales and a kinetic-Alfven-wave cascade at small scales, with strong small-scale density fluctuations and weak non-axisymmetric density waves. Ions undergo non-thermal particle acceleration, their distribution accurately described by a kappa distribution. These results have implications for the properties of low-collisionality accretion flows, such as that near the black hole at the Galactic center.
Hot accretion flows contain collisionless plasmas that are believed to be capable of accelerating particles to very high energies, as a result of turbulence generated by the magnetorotational instability (MRI). We conduct unstratified shearing-box simulations of the MRI turbulence in ideal magnetohydrodynamic (MHD), and inject energetic (relativistic) test particles in simulation snapshots to conduct a detailed investigation on particle diffusion and stochastic acceleration. We consider different amount of net vertical magnetic flux to achieve different disk magnetizations levels at saturated states, with sufficiently high resolution to resolve the gyro-radii ($R_g$) of most particles. Particles with large $R_g$ ($gtrsim0.03$ disk scale height $H$) show spatial diffusion coefficients of $sim30$ and $sim5$ times Bohm values in the azimuthal and poloidal directions, respectively. We further measure particle momentum diffusion coefficient $D(p)$ by applying the Fokker-Planck equation to particle momentum evolution. For these particles, contribution from turbulent fluctuations scales as $D(p)propto p$, and shear acceleration takes over when $R_ggtrsim0.1H$, characterized by $D(p)propto p^3$. For particles with smaller $R_g$ ($lesssim0.03H$), their spatial diffusion coefficients roughly scale as $sim p^{-1}$, and show evidence of $D(p)propto p^2$ scaling in momentum diffusion but with large uncertainties. We find that multiple effects contribute to stochastic acceleration/deceleration, and the process is also likely affected by intermittency in the MRI turbulence. We also discuss the potential of accelerating PeV cosmic-rays in hot accretion flows around supermassive black holes.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
