Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userTwo-Dimensional Hydrodynamic Simulations of Convection in Radiation-Dominated Accretion Disks
109
0
0.0
(
0
)
Authors
Eric Agol
Ask ChatGPT about the research
No Arabic abstract
The standard equilibrium for radiation-dominated accretion disks has long been known to be viscously, thermally, and convectively unstable, but the nonlinear development of these instabilities---hence the actual state of such disks---has not yet been identified. By performing local two-dimensional hydrodynamic simulations of disks, we demonstrate that convective motions can release heat sufficiently rapidly as to substantially alter the vertical structure of the disk. If the dissipation rate within a vertical column is proportional to its mass, the disk settles into a new configuration thinner by a factor of two than the standard radiation-supported equilibrium. If, on the other hand, the vertically-integrated dissipation rate is proportional to the vertically-integrated total pressure, the disk is subject to the well-known thermal instability. Convection, however, biases the development of this instability toward collapse. The end result of such a collapse is a gas pressure-dominated equilibrium at the original column density.
rate research
Read More
We use global three dimensional radiation magneto-hydrodynamic simulations to study the properties of inner regions of accretion disks around a 5times 10^8 solar mass black hole with mass accretion rates reaching 7% and 20% of the Eddington value. This region of the disk is supported by magnetic pressure with surface density significantly smaller than the values predicted by the standard thin disk model but with a much larger disk scale height. The disks do not show any sign of thermal instability over many thermal time scales. More than half of the accretion is driven by radiation viscosity in the optically thin corona region for the lower accretion rate case, while accretion in the optically thick part of the disk is driven by the Maxwell and Reynolds stresses from MRI turbulence. Coronae with gas temperatures > 10^8 K are generated only in the inner approx 10 gravitational radii in both simulations, being more compact in the higher accretion rate case. In contrast to the thin disk model, surface density increases with increasing mass accretion rate, which causes less dissipation in the optically thin region and a relatively weaker corona. The simulation results may explain the formation of X-ray coronae in Active Galactic Nuclei (AGNs), the compact size of such coronae, and the observed trend of optical to X-ray luminosity with Eddington ratio for many AGNs.
We examine the effects of photon bubble instability in radiation-dominated accretion disks such as those found around black holes in active galactic nuclei and X-ray binary star systems. Two- and 3-D numerical radiation MHD calculations of small patches of disk are used. Modes with wavelengths shorter than the gas pressure scale height grow faster than the orbital frequency in the surface layers. The fastest growth rate observed is five times the orbital frequency and occurs on nearly-vertical magnetic fields. The spectrum of linear modes agrees with a WKB analysis indicating still faster growth at unresolved scales, with a maximum proportional to the gravity and inversely proportional to the gas sound speed. Disturbances reaching non-linear amplitudes steepen into trains of shocks similar to a 1-D periodic non-linear analytic solution. Variations in propagation speed result in merging of adjacent fronts, and over time the shock spacing and amplitude increase. Growth is limited by the strength of the field, and the structure is disrupted when the ram pressure exceeds the magnetic pressure. The largest horizontal density variations are similar to the ratio of magnetic to gas pressure, and in our calculations are more than 100. Under the conditions considered, radiation diffuses through the inhomogeneneous flow 5x faster than through the initial hydrostatic equilibrium, and the net cooling rate is several times greater than in a similar calculation without magnetic fields showing the effects of convection. These results indicate photon bubbles may be important in cooling radiation-dominated disks.
(Abridged) Stars more massive than $20-30M_{odot}$ are so luminous that the radiation force on the cooler, more opaque outer layers can balance or exceed the force of gravity. These near or super-Eddington outer envelopes represent a long standing challenge for calculating the evolution of massive stars in one dimension, a situation that limits our understanding of the stellar progenitors of some of the most exciting and energetic explosions in the universe. In particular, the proximity to the Eddington limit has been the suspected cause for the variability, large mass loss rate and giant eruptions of an enigmatic class of massive stars: the luminous blue variables (LBVs). When in quiescence, LBVs are usually found on the hot ($T_{eff} approx 2 - 4 times 10^4$ K) S Dor instability strip. While in outburst, most LBVs stay on the cold S Dor instability strip with a $T_{eff} approx 9000$ K. Here we show that physically realistic three dimensional global radiation hydrodynamic simulations of radiation dominated massive stars with the largest supercomputers in the world naturally reproduce many observed properties of LBVs, specifically their location in the Hertzsprung-Russell (HR) diagram and their episodic mass loss with rates of $10^{-7}-10^{-5} M_{odot}/yr$. The helium opacity peak is found to play an important role to determine these properties, which is not realized in the traditional one dimensional models of massive stars. The simulations also predict that convection causes irregular envelope oscillations and 10-30% brightness variations on a typical timescale of a few days. The variability is more prominent in our models that are on the cool part of the S Dor instability. These calculations pave the way to a quantitative understanding of the structure, stability and the dominant mode of mass loss of massive stars.
Recently, the vertical shear instability (VSI) has become an attractive purely hydrodynamic candidate for the anomalous angular momentum transport required for weakly ionized accretion disks. In direct three-dimensional numerical simulations of VSI turbulence in disks, a meridional circulation pattern was observed that is opposite to the usual viscous flow behavior. Here, we investigate whether this feature can possibly be explained by an anisotropy of the VSI turbulence. Using three-dimensional hydrodynamical simulations, we calculate the turbulent Reynolds stresses relevant for angular momentum transport for a representative section of a disk. We find that the vertical stress is significantly stronger than the radial stress. Using our results in viscous disk simulations with different viscosity coefficients for the radial and vertical direction, we find good agreement with the VSI turbulence for the stresses and meridional flow; this provides additional evidence for the anisotropy. The results are important with respect to the transport of small embedded particles in disks.
We perform detailed variability analysis of two-dimensional viscous, radiation hydrodynamic numerical simulations of Shakura-Sunyaev thin disks around a stellar mass black hole. Disk models are initialized on both the gas-, as well as radiation-, pressure-dominated branches of the thermal equilibrium curve, with mass accretion rates spanning the range from $dot{M} = 0.01 L_mathrm{Edd}/c^2$ to $10 L_mathrm{Edd}/c^2$. An analysis of temporal variations of the numerically simulated disk reveals multiple robust, coherent oscillations. Considering the local mass flux variability, we find an oscillation occurring at the maximum radial epicyclic frequency, $3.5times 10^{-3},t_mathrm{g}^{-1}$, a possible signature of a trapped fundamental ${it g}$-mode. Although present in each of our simulated models, the trapped ${it g}$-mode feature is most prominent in the gas-pressure-dominated case. The total pressure fluctuations in the disk suggest strong evidence for standing-wave ${it p}$-modes, some trapped in the inner disk close to the ISCO, others present in the middle/outer parts of the disk. Knowing that the trapped ${it g}$-mode frequency and maximum radial epicyclic frequency differ by only $0.01%$ in the case of a non-rotating black hole, we simulated an additional initially gas-pressure-dominated disk with a dimensionless black hole spin parameter $a_* = 0.5$. The oscillation frequency in the spinning black hole case confirms that this oscillation is indeed a trapped ${it g}$-mode. All the numerical models we report here also show a set of high frequency oscillations at the vertical epicyclic and breathing mode frequencies. The vertical oscillations show a 3:2 frequency ratio with oscillations occurring approximately at the radial epicyclic frequency, which could be of astrophysical importance in observed twin peak, high-frequency quasi-periodic oscillations.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
