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 $\varepsilon=(L/L_{\rm Edd})/(\dot{M}/\dot{M}_{\rm 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.

Nonequilibrium–diffusion limit of the compressible Euler radiation model in R3

A Three-Dimensional Finite Volume Method for Conservation Laws in Conjunction with Modified Solution for Nonlinear Coupled Constitutive Relations

A radiation-hydrodynamics scheme valid from the transport to the diffusion limit

The Combined Non-Equilibrium Diffusion and Low Mach Number Limits of a Model Arising in Radiation Magnetohydrodynamics

Non-equilibrium diffusion limit of the compressible Navier-Stokes-Fourier-P1 approximation model at low Mach number: general initial data case

Diffusion asymptotics of a coupled model for radiative transfer: general initial data

Global Smooth Radially Symmetric Solutions to a Multidimensional Radiation Hydrodynamics Model

High order conservative Lagrangian schemes for two-dimensional radiation hydrodynamics equations in the equilibrium-diffusion limit

Compressible Euler Limit from Boltzmann Equation with Maxwell Reflection Boundary Condition in Half-Space

Equilibrium and Non-Equilibrium diffusion approximation for the radiative transfer equation

Radiative hydrodynamic equations with nonequilibrium radiative transfer

Non-relativistic Limit of the Euler-HMPN Approximation Model Arising in Radiation Hydrodynamics

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

3D hydrodynamic simulations of massive main-sequence stars. III. The effect of radiation pressure and diffusion leading to a 1D equilibrium model

Nonrelativistic Limit of the Compressible Navier--Stokes--Fourier--P1 Approximation Model Arising in Radiation Hydrodynamics

Non-Relativistic and Low Mach Number Limits of Two P1 Approximation Model Arising in Radiation Hydrodynamics

Implicit-explicit Runge-Kutta for radiation hydrodynamics I: gray diffusion

Kinetic modelling of rarefied gas flows with radiation

An asymptotic-preserving method for the three-temperature radiative transfer model

Quasineutral limit for the compressible two-fluid Euler–Maxwell equations for well-prepared initial data

DISSIPATIVE ENTROPY AND GLOBAL SMOOTH SOLUTIONS IN RADIATION HYDRODYNAMICS AND MAGNETOHYDRODYNAMICS