An Extension of the Athena++ Code Framework for Radiation-Magnetohydrodynamics in General Relativity Using a Finite-Solid-Angle Discretization

Abstract

We extend the general-relativistic magnetohydrodynamics (GRMHD) capabilities of Athena++ to incorporate radiation. The intensity field in each finite-volume cell is discretized in angle, with explicit transport in both space and angle properly accounting for the effects of gravity on null geodesics, and with matter and radiation coupled in a locally implicit fashion. Here we describe the numerical procedure in detail, verifying its correctness with a suite of tests. Motivated in particular by black hole accretion in the high-accretion-rate, thin-disk regime, we demonstrate the application of the method to this problem. With excellent scaling on flagship computing clusters, the port of the algorithm to the GPU-enabled AthenaK code now allows the simulation of many previously intractable radiation-GRMHD systems.

arXiv commentTo be published in ApJS

References in corpus (40)

• Modelling the behaviour of accretion flows in X-ray binaries
• Explosion Mechanisms of Core-Collapse Supernovae
• HARM: A Numerical Scheme for General Relativistic Magnetohydrodynamics
• The Athena++ Adaptive Mesh Refinement Framework: Design and Magnetohydrodynamic Solvers
• Binary neutron-star mergers: a review of Einstein's richest laboratory
• Kilonovae
• Three-Dimensional General Relativistic Radiation Magnetohydrodynamical Simulation of Super-Eddington Accretion, using a new code HARMRAD with M1 Closure
• A Finite Difference Representation of Neutrino Radiation Hydrodynamics in Spherically Symmetric General Relativistic Space-Time
• Binary Neutron Star Mergers
• Dimension as a Key to the Neutrino Mechanism of Core-Collapse Supernova Explosions
• Truncated Moment Formalism for Radiation Hydrodynamics in Numerical Relativity
• An Extension of the Athena++ Code Framework for GRMHD Based on Advanced Riemann Solvers and Staggered-Mesh Constrained Transport
• The Status of Multi-Dimensional Core-Collapse Supernova Models
• An Open-Source Neutrino Radiation Hydrodynamics Code for Core-Collapse Supernovae
• Semi-implicit scheme for treating radiation under M1 closure in general relativistic conservative fluid dynamics codes
• A new multi-dimensional general relativistic neutrino hydrodynamics code for core-collapse supernovae. I. Method and code tests in spherical symmetry
• Post-merger evolution of a neutron star-black hole binary with neutrino transport
• A New Multi-Energy Neutrino Radiation-Hydrodynamics Code in Full General Relativity and Its Application to Gravitational Collapse of Massive Stars
• General relativistic radiation hydrodynamics of accretion flows. I: Bondi-Hoyle accretion
• An Algorithm for Radiation Magnetohydrodynamics Based on Solving the Time-dependent Transfer Equation
• Relativistic Radiation Magnetohydrodynamics in Dynamical Spacetimes: Numerical Methods and Tests
• A New Spherical Harmonics Scheme for Multi-Dimensional Radiation Transport I: Static Matter Configurations
• Formation of Overheated Regions and Truncated Disks around Black Holes; Three-dimensional General Relativistic Radiation-magnetohydrodynamics Simulations
• A New Code for Proto-Neutron Star Evolution
• Numerical Simulations of Optically Thick Accretion onto a Black Hole - II. Rotating Flow
• An Implicit Finite Volume Scheme to Solve the Time Dependent Radiation Transport Equation Based on Discrete Ordinates
• Monte-Carlo closure for moment-based transport schemes in general relativistic radiation hydrodynamics simulations
• Three-dimensional Boltzmann-Hydro code for core-collapse in massive stars II. The Implementation of moving-mesh for neutron star kicks
• General relativistic radiation hydrodynamics of accretion flows: II. Treating stiff source terms and exploring physical limitations
• GRAthena++: puncture evolutions on vertex-centered oct-tree AMR
• Two-moment scheme for general-relativistic radiation hydrodynamics: a systematic description and new applications
• Multi-Frequency General Relativistic Radiation-Hydrodynamics with ${\mathbf{M}}_1$ Closure
• MHD modeling on geodesic grids
• Shock wave structure in astrophysical flows with an account of photon transfer
• A Novel Multi-Dimensional Boltzmann Neutrino Transport Scheme for Core-Collapse Supernovae
• MOCMC: Method of Characteristics Moment Closure, a Numerical Method for Covariant Radiation Magnetohydrodynamics
• NADA-FLD: A General Relativistic, Multi-dimensional Neutrino-hydrodynamics Code Employing Flux-limited Diffusion
• Covariant Radiative Transfer for Black Hole Spacetimes
• A Monte-Carlo based relativistic radiation hydrodynamics code with a higher-order scheme
• Radiative MHD Simulations of Photon Bubbles in Radiation-Supported Magnetized Atmospheres of Neutron Stars with Isotropic Thomson Scattering

Cited by in corpus (6)

• General-relativistic radiation transport scheme in $\text{Gmunu}$ I: Implementation of two-moment based multi-frequency radiative transfer and code tests
• The plunging region of a thin accretion disc around a Schwarzschild black hole
• General Relativistic Radiation Magnetohydrodynamics Simulations of Precessing Tilted Super-Eddington Disks
• KHARMA: Flexible, Portable Performance for GRMHD
• A User-Friendly Python Interface for the Numerical Relativity Code AMSS-NCKU
• Convective stability analysis of massive neutron stars formed in binary mergers