Two-moment scheme for general-relativistic radiation hydrodynamics: a systematic description and new applications

Abstract

We provide a systematic description of the steps necessary -- and of the potential pitfalls to be encountered -- when implementing a two-moment scheme within an Implicit-Explicit (IMEX) scheme to include radiative-transfer contributions in numerical simulations of general-relativistic (magneto-)hydrodynamics. We make use of the M1 closure, which provides an exact solution for the optically thin and thick limit, and an interpolation between these limits. Special attention is paid to the efficient solution of the emerging set of implicit conservation equations. In particular, we present an efficient method for solving these equations via the inversion of a $4\times 4$-matrix within an IMEX scheme. While this method relies on a few approximations, it offers a very good compromise between accuracy and computational efficiency. After a large number of tests in special relativity, we couple our new radiation code, \texttt{FRAC}, with the general-relativistic magnetohydrodynamics code \texttt{BHAC} to investigate the radiative Michel solution, namely, the problem of spherical accretion onto a black hole in the presence of a radiative field. By performing the most extensive exploration of the parameter space for this problem, we find that the accretion's efficiency can be expressed in terms of physical quantities such as temperature, $T$, luminosity, $L$, and black-hole mass, $M$, via the expression $\epsilon =(L/L_{\mathrm{Edd}})/(\dot{M}/{\dot{M}}_{\mathrm{Edd}})=7.41\times 10^{-7}{\left(T/10^6\mathrm{K}\right)}^{0.22}{\left(L/L_{\odot }\right)}^{0.48}{\left(M/M_{\odot }\right)}^{0.48}$, where $L_{\mathrm{Edd}}$ and ${\dot{M}}_{\mathrm{Edd}}$ are the Eddington luminosity and accretion rate, respectively. Finally, we also consider the accretion problem away from spherical symmetry, finding that the solution is stable under perturbations in the radiation field.

arXiv comment22 pages, 15 figures, matches version accepted to MNRAS

References in corpus (76)

• GW170817: Observation of Gravitational Waves from a Binary Neutron Star Inspiral
• First M87 Event Horizon Telescope Results. I. The Shadow of the Supermassive Black Hole
• Multi-messenger Observations of a Binary Neutron Star Merger
• Gravitational Waves and Gamma-rays from a Binary Neutron Star Merger: GW170817 and GRB 170817A
• First M87 Event Horizon Telescope Results. VI. The Shadow and Mass of the Central Black Hole
• First M87 Event Horizon Telescope Results. V. Physical Origin of the Asymmetric Ring
• Gravitational-wave constraints on the neutron-star-matter Equation of State
• Constraining the Maximum Mass of Neutron Stars From Multi-Messenger Observations of GW170817
• Using gravitational-wave observations and quasi-universal relations to constrain the maximum mass of neutron stars
• New constraints on radii and tidal deformabilities of neutron stars from GW170817
• Foundations of Black Hole Accretion Disk Theory
• The Einstein Toolkit: A Community Computational Infrastructure for Relativistic Astrophysics
• GW170817: Joint Constraint on the Neutron Star Equation of State from Multimessenger Observations
• Neutron-star radius constraints from GW170817 and future detections
• Binary neutron-star mergers: a review of Einstein's richest laboratory
• Neutrino-driven winds from neutron star merger remnants
• GW170817, General Relativistic Magnetohydrodynamic Simulations, and the Neutron Star Maximum Mass
• Signatures of quark-hadron phase transitions in general-relativistic neutron-star mergers
• A critical examination of constraints on the equation of state of dense matter obtained from GW170817
• Three-Dimensional General Relativistic Radiation Magnetohydrodynamical Simulation of Super-Eddington Accretion, using a new code HARMRAD with M1 Closure
• Three-dimensional GRMHD simulations of the remnant accretion disks from neutron star mergers: outflows and r-process nucleosynthesis
• Dynamical mass ejection from the merger of asymmetric binary neutron stars: Radiation-hydrodynamics study in general relativity
• Constraint on the maximum mass of neutron stars using GW170817 event
• The Current Ability to Test Theories of Gravity with Black Hole Shadows
• High Resolution Calculations of Merging Neutron Stars II: Neutrino Emission
• The long-term evolution of neutron star merger remnants - I. The impact of r-process nucleosynthesis
• Long-term GRMHD Simulations of Neutron Star Merger Accretion Disks: Implications for Electromagnetic Counterparts
• The Black Hole Accretion Code
• Multi-messenger Bayesian parameter inference of a binary neutron-star merger
• AT2017gfo: an anisotropic and three-component kilonova counterpart of GW170817
• Mass Ejection from the Remnant of a Binary Neutron Star Merger: Viscous-Radiation Hydrodynamics Study
• Truncated Moment Formalism for Radiation Hydrodynamics in Numerical Relativity
• An Open-Source Neutrino Radiation Hydrodynamics Code for Core-Collapse Supernovae
• Low mass binary neutron star mergers : gravitational waves and neutrino emission
• Beyond ideal MHD: towards a more realistic modeling of relativistic astrophysical plasmas
• Semi-implicit scheme for treating radiation under M1 closure in general relativistic conservative fluid dynamics codes
• A Module for Radiation Hydrodynamic Calculations With ZEUS-2D Using Flux-Limited Diffusion
• General-relativistic resistive magnetohydrodynamics in three dimensions: Formulation and tests
• Full Transport Model of GW170817-Like Disk Produces a Blue Kilonova
• Post-merger evolution of a neutron star-black hole binary with neutrino transport
• Impact of an improved neutrino energy estimate on outflows in neutron star merger simulations
• On r-process nucleosynthesis from matter ejected in binary neutron starmergers
• Two-temperature, Magnetically Arrested Disc simulations of the jet from the supermassive black hole in M87
• Modeling dynamical ejecta from binary neutron star mergers and implications for electromagnetic counterparts
• General relativistic simulations of compact binary mergers as engines of short gamma-ray bursts
• A new multidimensional, energy-dependent two-moment transport code for neutrino-hydrodynamics
• A General-relativistic Determination of the Threshold Mass to Prompt Collapse in Binary Neutron Star Mergers
• Implementation of a simplified approach to radiative transfer in general relativity
• A New Multi-Energy Neutrino Radiation-Hydrodynamics Code in Full General Relativity and Its Application to Gravitational Collapse of Massive Stars
• Constrained transport and adaptive mesh refinement in the Black Hole Accretion Code
• General relativistic radiation hydrodynamics of accretion flows. I: Bondi-Hoyle accretion
• Beyond second-order convergence in simulations of magnetised binary neutron stars with realistic microphysics
• Are small radii of compact stars ruled out by GW170817/AT2017gfo?
• Accurate Simulations of Binary Black-Hole Mergers in Force-Free Electrodynamics
• Pair Production in Low Luminosity Galactic Nuclei
• Modeling non-thermal emission from the jet-launching region of M 87 with adaptive mesh refinement
• A New Spherical Harmonics Scheme for Multi-Dimensional Radiation Transport I: Static Matter Configurations
• GRMHD simulations of accretion onto Sgr A*: How important are radiative losses?
• General relativistic resistive magnetohydrodynamics with robust primitive variable recovery for accretion disk simulations
• Electromagnetic emission from long-lived binary neutron star merger remnants I: formulation of the problem
• Conservative 3+1 General Relativistic Variable Eddington Tensor Radiation Transport Equations
• Numerical Simulations of Optically Thick Accretion onto a Black Hole - I. Spherical Case
• Numerical Simulations of Optically Thick Accretion onto a Black Hole - II. Rotating Flow
• Monte-Carlo closure for moment-based transport schemes in general relativistic radiation hydrodynamics simulations
• Analytic Closures for M1 Neutrino Transport
• Evaluating radiation transport errors in merger simulations using a Monte-Carlo algorithm
• General relativistic radiation hydrodynamics of accretion flows: II. Treating stiff source terms and exploring physical limitations
• Relativistic Radiative Transfer for Spherical Flows
• Hyperbolic character of the angular moment equations of radiative transfer and numerical methods
• $\nu \text{bhlight}$: Radiation GRMHD for Neutrino-Driven Accretion Flows
• A radiative transfer module for relativistic magnetohydrodynamics in the PLUTO code
• Time--dependent analysis of spherical accretion onto black holes
• NADA-FLD: A General Relativistic, Multi-dimensional Neutrino-hydrodynamics Code Employing Flux-limited Diffusion
• Choked accretion onto a Schwarzschild black hole: A hydrodynamical jet-launching mechanism
• Outflows from inflows: the nature of Bondi-like accretion
• CAFE-R a code that solves the special relativistic radiation hydrodynamics equations

Cited by in corpus (22)

• MPI-AMRVAC 3.0: updates to an open-source simulation framework
• Neutrino transport in general relativistic neutron star merger simulations
• General-relativistic hydrodynamics of non-perfect fluids: 3+1 conservative formulation and application to viscous black-hole accretion
• Separating Astrophysics and Geometry in Black Hole Images
• 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
• Multi-Frequency General Relativistic Radiation-Hydrodynamics with ${\mathbf{M}}_1$ Closure
• A two-moment radiation hydrodynamics scheme applicable to simulations of planet formation in circumstellar disks
• General-relativistic radiation transport scheme in $\text{Gmunu}$ I: Implementation of two-moment based multi-frequency radiative transfer and code tests
• An extension of Gmunu: General-relativistic resistive magnetohydrodynamics based on staggered-meshed constrained transport with elliptic cleaning
• Spherical accretion: Bondi, Michel and rotating black holes
• Energy-dependent and energy-integrated two-moment general-relativistic neutrino transport simulations of hypermassive neutron star
• The Guided Moments formalism: a new efficient full-neutrino treatment for astrophysical simulations
• Choked accretion onto a Kerr black hole
• Relativistic radiation hydrodynamics in a reference-metric formulation
• Global high-order numerical schemes for the time evolution of the general relativistic radiation magneto-hydrodynamics equations
• On the application of Jacobian-free Riemann solvers for relativistic radiation magnetohydrodynamics under M1 closure
• Multidimensional half-moment multigroup radiative transfer. Improving moment-based thermal models of circumstellar disks
• Simulating Binary Neutron Star Mergers with Finite-temperature Equations of State: The influences of the slope of the symmetry energy and artificial heating
• Fully general-relativistic simulations of isolated and binary strange quark stars
• A New Moment-Based General-Relativistic Neutrino-Radiation Transport Code: Methods and First Applications to Neutron Star Mergers
• Incorporating a radiative hydrodynamics scheme in the numerical-relativity code BAM