Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userEquilibrium disks, MRI mode excitation, and steady state turbulence in global accretion disk simulations
117
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
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).
rate research
Read More
Magnetorotational turbulence provides a viable mechanism for angular momentum transport in accretion disks. We present global, three dimensional (3D), MHD accretion disk simulations that investigate the dependence of the turbulent stresses on resolution. Convergence in the time-and-volume-averaged stress-to-gas-pressure ratio, at a value of $sim0.04$, is found for a model with radial, vertical, and azimuthal resolution of 12-51, 27, and 12.5 cells per scale-height (the simulation mesh is such that cells per scale-height varies in the radial direction). A control volume analysis is performed on the main body of the disk (|z|<2H) to examine the production and removal of magnetic energy. Maxwell stresses in combination with the mean disk rotation are mainly responsible for magnetic energy production, whereas turbulent dissipation (facilitated by numerical resistivity) predominantly removes magnetic energy from the disk. Re-casting the magnetic energy equation in terms of the power injected by Maxwell stresses on the boundaries of, and by Lorentz forces within, the control volume highlights the importance of the boundary conditions (of the control volume). The different convergence properties of shearing-box and global accretion disk simulations can be readily understood on the basis of choice of boundary conditions and the magnetic field configuration. Periodic boundary conditions restrict the establishment of large-scale gradients in the magnetic field, limiting the power that can be delivered to the disk by Lorentz forces and by stresses at the surfaces. The factor of three lower resolution required for convergence in turbulent stresses for our global disk models compared to stratified shearing-boxes is explained by this finding. (Abridged)
We use global magnetohydrodynamic simulations to study the influence of net vertical magnetic fields on the structure of geometrically thin ($H/r approx 0.05$) accretion disks in the Newtonian limit. We consider initial mid-plane gas to magnetic pressure ratios $beta_0 = 1000,, 300$ and $100$, spanning the transition between weakly and strongly magnetized accretion regimes. We find that magnetic pressure is important for the disks vertical structure in all three cases, with accretion occurring at $z/Rapprox 0.2$ in the two most strongly magnetized models. The disk midplane shows outflow rather than accretion. Accretion through the surface layers is driven mainly by stress due to coherent large scale magnetic field rather than by turbulent stress. Equivalent viscosity parameters measured from our simulations show similar dependencies on initial $beta_0$ to those seen in shearing box simulations, though the disk midplane is not magnetic pressure dominated even for the strongest magnetic field case. Winds are present but are not the dominant driver of disk evolution. Over the (limited) duration of our simulations, we find evidence that the net flux attains a quasi-steady state at levels that can stably maintain a strongly magnetized disk. We suggest that geometrically thin accretion disks in observed systems may commonly exist in a magnetically elevated state, characterized by non-zero but modest vertical magnetic fluxes, with potentially important implications for disk phenomenology in X-ray binaries (XRBs) and active galactic nuclei (AGN).
In this paper we consider two outstanding intertwined problems in modern high-energy astrophysics: (1) the vertical thermal structure of an optically thick accretion disk heated by the dissipation of magnetohydrodynamic (MHD) turbulence driven by the magneto-rotational instability (MRI), and (2) determining the fraction of the accretion power released in the corona above the disk. For simplicity, we consider a gas-pressure-dominated disk and assume a constant opacity. We argue that the local turbulent dissipation rate due to the disruption of MRI channel flows by secondary parasitic instabilities should be uniform across most of the disk, almost up to the disk photosphere. We then obtain a self-consistent analytical solution for the vertical thermal structure of the disk, governed by the balance between the heating by MRI turbulence and the cooling by radiative diffusion. Next, we argue that the coronal power fraction is determined by the competition between the Parker instability, viewed as a parasitic instability feeding off of MRI channel flows, and other parasitic instabilities. We show that the Parker instability inevitably becomes important near the disk surface, leading to a certain lower limit on the coronal power. While most of the analysis in this paper focuses on the case of a disk threaded by an externally imposed vertical magnetic field, we also discuss the zero-net-flux case, in which the magnetic field is produced by the MRI dynamo itself, and show that most of our arguments and conclusions should be valid in this case as well.
It is widely accepted that quasars and other active galactic nuclei (AGN) are powered by accretion of matter onto a central supermassive black hole. While numerical simulations have demonstrated the importance of magnetic fields in generating the turbulence believed necessary for accretion, so far they have not produced the high mass accretion rates required to explain the most powerful sources. We describe new global 3D simulations we are developing to assess the importance of radiation and non-ideal MHD in generating magnetized outflows that can enhance the overall rates of angular momentum transport and mass accretion.
In the course of its evolution, a black hole (BH) accretes gas from a wide range of directions. Given a random accretion event, the typical angular momentum of an accretion disc would be tilted by $sim$60$^circ$ relative to the BH spin. Misalignment causes the disc to precess at a rate that increases with BH spin and tilt angle. We present the first general-relativistic magnetohydrodynamic (GRMHD) simulations spanning a full precession period of highly tilted (60$^circ$), moderately thin ($h/r=0.1$) accretion discs around a rapidly spinning ($asimeq0.9$) BH. While the disc and jets precess in phase, we find that the corona, sandwiched between the two, lags behind by $gtrsim 10^{circ}$. For spectral models of BH accretion, the implication is that hard non-thermal (corona) emission lags behind the softer (disc) emission, thus potentially explaining some properties of the hard energy lags seen in Type-C low frequency quasi-periodic oscillations in X-Ray binaries. While strong jets are unaffected by this disc-corona lag, weak jets stall when encountering the lagging corona at distances $r sim 100$ black hole radii. This interaction may quench large-scale jet formation.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
