Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userBreathing Oscillations in a Global Simulation of a Thin Accretion Disk
75
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We study the oscillations of an axisymmetric, viscous, radiative, general relativistic hydrodynamical simulation of a geometrically thin disk around a non-rotating, $6.62,M_odot$ black hole. The numerical setup is initialized with a Novikov-Thorne, gas-pressure-dominated accretion disk, with an initial mass accretion rate of $dot{m} = 0.01,L_mathrm{Edd}/c^2$ (where $L_mathrm{Edd}$ is the Eddington luminosity and $c$ is the speed of light). Viscosity is treated with the $alpha$-prescription. The simulation was evolved for about $1000$ Keplerian orbital periods at three Schwarzschild radii (ISCO radius). Power density spectra of the radial and vertical fluid velocity components, the total (gas $+$ radiation) midplane pressure, and the vertical component of radiative flux from the photosphere, all reveal strong power at the local breathing oscillation frequency. The first, second and third harmonics of the breathing oscillation are also clearly seen in the data. We quantify the properties of these oscillations by extracting eigenfunctions of the radial and vertical velocity components and total pressure. This confirms that these oscillations are associated with breathing motion.
rate research
Read More
We present the results of a 3D global magnetohydrodynamic (MHD) simulation of an AM CVn system that was aimed at exploring eccentricity growth in the accretion disc self-consistently from a first principles treatment of the MHD turbulence. No significant eccentricity growth occurs in the simulation. In order to investigate the reasons why, we ran 2D alpha disc simulations with alpha values of 0.01, 0.1, and 0.2, and found that only the latter two exhibit significant eccentricity growth. We present an equation expressing global eccentricity evolution in terms of contributing forces and use it to analyze the simulations. As expected, we find that the dominant term contributing to the growth of eccentricity is the tidal gravity of the companion star. In the 2D simulations, the alpha viscosity directly contributes to eccentricity growth. In contrast, the overall magnetic forces in the 3D simulation damp eccentricity. We also analyzed the mode-coupling mechanism of Lubow, and confirmed that the spiral wave excited by the 3:1 resonance was the dominant contributor to eccentricity growth in the 2D $alpha=0.1$ simulations, but other waves also contribute significantly. We found that the $alpha=0.1$ and 0.2 simulations had more relative mass at larger radii compared to the $alpha=0.01$ and 3D MHD simulation, which also had an effective $alpha$ of 0.01. This suggests that in 3D MHD simulations without sufficient poloidal magnetic flux, MRI turbulence does not saturate at a high enough $alpha$ to spread the disc to large enough radii to reproduce the superhumps observed in real systems.
141 -
Nagendra Kumar
, Banibrata Mukhopadhyay (Department of Physics
, n Indian Institute of Science
2021
Many astrophysical sources, e.g., cataclysmic variables, X-ray binaries, active galactic nuclei, exhibit a wind outflow, when they reveal a multicolor blackbody spectrum, hence harboring a geometrically thin Keplerian accretion disk. Unlike an advective disk, in the thin disk, the physical environment, like, emission line, external heating, is expected to play a key role to drive the wind outflow. We show the wind outflow in a thin disk attributing a disk irradiation effect, probably from the inner to outer disks. We solve the set of steady, axisymmetric disk model equations in cylindrical coordinates along the vertical direction for a given launching radius $(r)$ from the midplane, introducing irradiation as a parameter. We obtain an acceleration solution, for a finite irradiation in the presence of a fixed but tiny initial vertical velocity (hence thin disk properties practically do not alter) at the midplane, upto a maximum height ($z^{max}$). We find that wind outflow mainly occurs from the outer region of the disk and its density decreases with increasing launching radius, and for a given launching radius with increasing ejection height. Wind power decreases with increasing ejection height. For $z^{max} < 2r$, wind outflow is ejected tangentially (or parallel to the disk midplane) in all directions with the fluid speed same as the azimuthal speed. This confirms mainly, for low mass X-ray binaries, (a) wind outflow should be preferentially observed in high-inclination sources, (b) the expectation of red and blue shifted absorption lines.
We investigated a semi-analytic and numerical model to study the geometrically thin and optically thick accretion disk around Maclaurin spheroid (MS). The main interest is in the inner region of the so called {alpha}-disk, {alpha} being the viscosity parameter. Analytical calculations are done assuming radiation pressure and gas pressure dominated for close to Eddington mass accretion rate and $dot{M}lesssim 0.1dot{M_{Edd}}$ respectively. We found that the change in eccentricity of MS gives a change at high frequency region in the emitted spectra. We found that disk parameters are dependent on eccentricity of MS. Our semi-analytic results show that qualitatively an increase in eccentricity of MS has same behavior as decrease in mass accretion rate. Numerical work has been carried out to see the viscous time evolution of the accretion disk around MS. In numerical model we showed that if the eccentricity of the object is high the matter will diffuse slowly during its viscous evolution. This gives a clue that how spin-up or spin-down can change the time evolution of the accretion disk using a simple Newtonian approach. The change in spectra can be used to determine the eccentricity of MS and thus period of the MS.
We present the results of a global, three-dimensional magnetohydrodynamics simulation of an accretion disk with a rotating, weakly magnetized central star. The disk is threaded by a weak, large-scale poloidal magnetic field, and the central star has no strong stellar magnetosphere initially. Our simulation investigates the structure of the accretion flows from a turbulent accretion disk onto the star. The simulation reveals that fast accretion onto the star at high latitudes occurs even without a stellar magnetosphere. We find that the failed disk wind becomes the fast, high-latitude accretion as a result of angular momentum exchange mediated by magnetic fields well above the disk, where the Lorentz force that decelerates the rotational motion of gas can be comparable to the centrifugal force. Unlike the classical magnetospheric accretion scenario, fast accretion streams are not guided by magnetic fields of the stellar magnetosphere. Nevertheless, the accretion velocity reaches the free-fall velocity at the stellar surface due to the efficient angular momentum loss at a distant place from the star. This study provides a possible explanation why Herbig Ae/Be stars whose magnetic fields are generally not strong enough to form magnetospheres also show indications of fast accretion. A magnetically driven jet is not formed from the disk in our model. The differential rotation cannot generate sufficiently strong magnetic fields for the jet acceleration because the Parker instability interrupts the field amplification.
We examine the properties of strongly magnetized accretion discs in a global framework, with particular focus on the evolution of magnetohydrodynamic instabilities such as the magnetorotational instability (MRI). Work by Pessah and Psaltis showed that MRI is stabilized beyond a critical toroidal field in compressible, differentially rotating flows and, also, reported the appearance of two new instabilities beyond this field. Their results stemmed from considering geometric curvature effects due to the suprathermal background toroidal field, which had been previously ignored in weak-field studies. However, their calculations were performed under the local approximation, which poses the danger of introducing spurious behavior due to the introduction of global geometric terms in an otherwise local framework. In order to avoid this, we perform a global eigenvalue analysis of the linearized MHD equations in cylindrical geometry. We confirm that MRI indeed tends to be highly suppressed when the background toroidal field attains the Pessah-Psaltis limit. We also observe the appearance of two new instabilities that emerge in the presence of highly suprathermal toroidal fields. These results were additionally verified using numerical simulations in PLUTO. There are, however, certain differences between the the local and global results, especially in the vertical wavenumber occupancies of the various instabilities, which we discuss in detail. We also study the global eigenfunctions of the most unstable modes in the suprathermal regime, which are inaccessible in the local analysis. Overall, our findings emphasize the necessity of a global treatment for accurately modeling strongly magnetized accretion discs.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
