Numerical Simulations of Optically Thick Accretion onto a Black Hole - II. Rotating Flow

Abstract

In this paper we report on recent upgrades to our general relativistic radiation magnetohydrodynamics code, Cosmos++, including the development of a new primitive inversion scheme and a hybrid implicit-explicit solver with a more general closure relation for the radiation equations. The new hybrid solver helps stabilize the treatment of the radiation source terms, while the new closure allows for a much broader range of optical depths to be considered. These changes allow us to expand by orders of magnitude the range of temperatures, opacities, and mass accretion rates, and move a step closer toward our goal of performing global simulations of radiation-pressure-dominated black hole accretion disks. In this work we test and validate the new method against an array of problems. We also demonstrate its ability to handle super-Eddington, quasi-spherical accretion. Even with just a single proof-of-principle simulation, we already see tantalizing hints of the interesting phenomenology associated with the coupling of radiation and gas in super-Eddington accretion flows.

18 pages, 9 figures

References in corpus (1)

• Semi-implicit scheme for treating radiation under M1 closure in general relativistic conservative fluid dynamics codes

Cited by in corpus (27)

• First Sagittarius A* Event Horizon Telescope Results. V. Testing Astrophysical Models of the Galactic Center Black Hole
• Iharm3D: Vectorized General Relativistic Magnetohydrodynamics
• Beyond moments: relativistic Lattice-Boltzmann methods for radiative transport in computational astrophysics
• An Extension of the Athena++ Code Framework for Radiation-Magnetohydrodynamics in General Relativity Using a Finite-Solid-Angle Discretization
• The Role of Strong Magnetic Fields in Stabilizing Highly Luminous, Thin Disks
• Common-envelope Dynamics of a Stellar-mass Black Hole: General Relativistic Simulations
• Large-scale outflow structure and radiation properties of super-Eddington flow: Dependence on the accretion rates
• Multi-Frequency General Relativistic Radiation-Hydrodynamics with ${\mathbf{M}}_1$ Closure
• Global Three-Dimensional Radiation Magnetohydrodynamic Simulations of Accretion onto a Stellar Mass Black Hole at Sub- and Near-critical Accretion Rates
• General-relativistic radiation transport scheme in $\text{Gmunu}$ I: Implementation of two-moment based multi-frequency radiative transfer and code tests
• Looking for the underlying cause of black hole X-ray variability in GRMHD simulations
• CosmosDG: An hp-Adaptive Discontinuous Galerkin Code for Hyper-Resolved Relativistic MHD
• Estimating Photometric Redshifts for X-ray sources in the X-ATLAS field, using machine-learning techniques
• Covariant Radiative Transfer for Black Hole Spacetimes
• General Relativistic Implicit Monte Carlo Radiation-Hydrodynamics
• Local models of two-temperature accretion disc coronae. I. Structure, outflows, and energetics
• Performance-Portable Binary Neutron Star Mergers with AthenaK
• CAFE-R a code that solves the special relativistic radiation hydrodynamics equations
• Two-dimensional Inflow-outflow Solution of Supercritical Accretion Flow
• Nuclear Ignition of White Dwarf Stars by Relativistic Encounters with Rotating Intermediate Mass Black Holes
• Simulations of Tidal Disruption Events
• On the application of Jacobian-free Riemann solvers for relativistic radiation magnetohydrodynamics under M1 closure
• Dispersion relations of relativistic radiation hydrodynamics
• Polarization Signatures of Quasi-Periodic Oscillations in Simulated Tilted, Truncated Disks
• Radiative GRMHD simulations of accretion disks
• Resolution Study of Thermonuclear Initiation in White Dwarf Tidal Disruption Events
• Non-Newtonian corrections to radiative viscosity: Israel-Stewart theory as a viscosity limiter