اشترك بالحزمة الذهبية واحصل على وصول غير محدود شمرا أكاديميا
تسجيل مستخدم جديد$ utexttt{bhlight}$: Radiation GRMHD for Neutrino-Driven Accretion Flows
211
0
0.0
(
0
)
اسأل ChatGPT حول البحث
ﻻ يوجد ملخص باللغة العربية
The 2017 detection of the in-spiral and merger of two neutron stars was a landmark discovery in astrophysics. We now know that such mergers are central engines of short gamma ray bursts and sites of r-process nucleosynthesis, where the heaviest elements in our universe are formed. In the coming years, we expect many more such mergers. Modeling such systems presents a significant computational challenge along with the observational one. To meet this challenge, we present $ utexttt{bhlight}$, a scheme for solving general relativistic magnetohydrodynamics with energy-dependent neutrino transport in full (3+1)-dimensions, facilitated by Monte Carlo methods. We present a suite of tests demonstrating the accuracy, efficacy, and necessity of our scheme. We demonstrate the potential of our scheme by running a sample calculation in a domain of interest---the dynamics and composition of the accretion disk formed by a binary neutron star merger.
قيم البحث
اقرأ أيضاً
We present the implementation of an implicit-explicit (IMEX) Runge-Kutta numerical scheme for general relativistic hydrodynamics coupled to an optically thick radiation field in two existing GR-hydrodynamics codes. We argue that the necessity of such
an improvement arises naturally in astrophysically relevant regimes where the optical thickness is high as the equations become stiff. By performing several 1D tests we verify the codes new ability to deal with this stiffness and show consistency. Then, still in 1D, we compute a luminosity versus accretion rate diagram for the setup of spherical accretion onto a Schwarzschild black hole and find good agreement with previous work. Lastly, we revisit the supersonic Bondi Hoyle Lyttleton (BHL) accretion in 2D where we can now present simulations of realistic temperatures, down to T~10^6 K. Here we find that radiation pressure plays an important role, but also that these highly dynamical set-ups push our approximate treatment towards the limit of physical applicability. The main features of radiation hydrodynamics BHL flows manifest as (i) an effective adiabatic index approaching gamma_effective ~ 4/3; (ii) accretion rates two orders of magnitude lower than without radiation pressure; (iii) luminosity estimates around the Eddington limit, hence with an overall radiative efficiency as small as eta ~ 10^{-2}; (iv) strong departures from thermal equilibrium in shocked regions; (v) no appearance of the flip-flop instability. We conclude that the current optically thick approximation to the radiation transfer does give physically substantial improvements over the pure hydro also in set-ups departing from equilibrium, and, once accompanied by an optically thin treatment, is likely to provide a fundamental tool for investigating accretion flows in a large variety of astrophysical systems.
We present 3D general relativistic magnetohydrodynamic(GRMHD) simulations of zero angular momentum accretion around a rapidly rotating black hole, modified by the presence of initially uniform magnetic fields. We consider serveral angles between the
magnetic field direction and the black hole spin. In the resulting flows, the midplane dynamics are governed by magnetic reconnection-driven turbulence in a magnetically arrested (or a nearly arrested) state. Electromagnetic jets with outflow efficiencies ~10-200% occupy the polar regions, reaching several hundred gravitational radii before they dissipate due to the kink instability. The jet directions fluctuate in time and can be tilted by as much as ~30 degrees with respect to black hole spin, but this tilt does not depend strongly on the tilt of the initial magnetic field. A jet forms even when there is no initial net vertical magnetic flux since turbulent, horizon-scale fluctuations can generate a net vertical field locally. Peak jet power is obtained for an initial magnetic field tilted by 40-80 degrees with respect to the black hole spin because this maximizes the amount of magnetic flux that can reach the black hole. These simulations may be a reasonable model for low luminosity black hole accretion flows such as Sgr A* or M87.
The underlying hypothesis of this work is that the active galactic nuclei (AGNs) are wormhole mouths rather than supermassive black holes (SMBHs). Under some - quite general - assumptions such wormholes may emit gamma radiation as a result of a colli
sion of accreting flows inside the wormholes. This radiation has a distinctive spectrum much different from those of jets or accretion disks of AGNs. An observation of such radiation would serve as evidence of the existence of wormholes.
Monte Carlo approaches to radiation transport have several attractive properties such as simplicity of implementation, high accuracy, and good parallel scaling. Moreover, Monte Carlo methods can handle complicated geometries and are relatively easy t
o extend to multiple spatial dimensions, which makes them potentially interesting in modeling complex multi-dimensional astrophysical phenomena such as core-collapse supernovae. The aim of this paper is to explore Monte Carlo methods for modeling neutrino transport in core-collapse supernovae. We generalize the Implicit Monte Carlo photon transport scheme of Fleck & Cummings and gray discrete-diffusion scheme of Densmore et al. to energy-, time-, and velocity-dependent neutrino transport. Using our 1D spherically-symmetric implementation, we show that, similar to the photon transport case, the implicit scheme enables significantly larger timesteps compared with explicit time discretization, without sacrificing accuracy, while the discrete-diffusion method leads to significant speed-ups at high optical depth. Our results suggest that a combination of spectral, velocity-dependent, Implicit Monte Carlo and discrete-diffusion Monte Carlo methods represents a robust approach for use in neutrino transport calculations in core-collapse supernovae. Our velocity-dependent scheme can easily be adapted to photon transport.
We present bhlight, a numerical scheme for solving the equations of general relativistic radiation magnetohydrodynamics (GRRMHD) using a direct Monte Carlo solution of the frequency-dependent radiative transport equation. bhlight is designed to evolv
e black hole accretion flows at intermediate accretion rate, in the regime between the classical radiatively efficient disk and the radiatively inefficient accretion flow (RIAF), in which global radiative effects play a sub-dominant but non-negligible role in disk dynamics. We describe the governing equations, numerical method, idiosyncrasies of our implementation, and a suite of test and convergence results. We also describe example applications to radiative Bondi accretion and to a slowly accreting Kerr black hole in axisymmetry.
سجل دخول لتتمكن من نشر تعليقات
التعليقات
جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
