On the application of Jacobian-free Riemann solvers for relativistic radiation magnetohydrodynamics under M1 closure
arXiv:2212.00370 · doi:10.1016/j.cpc.2022.108630
Abstract
Radiative transfer plays a major role in high-energy astrophysics. In multiple scenarios and in a broad range of energy scales, the coupling between matter and radiation is essential to understand the interplay between theory, observations and numerical simulations. In this paper, we present a novel scheme for solving the equations of radiation relativistic magnetohydrodynamics within the parallel code Lóstrego. These equations, which are formulated taking successive moments of the Boltzmann radiative transfer equation, are solved under the gray-body approximation and the M1 closure using an IMEX time integration scheme. The main novelty of our scheme is that we introduce for the first time in the context of radiation magnetohydrodynamics a family of Jacobian-free Riemann solvers based on internal approximations to the Polynomial Viscosity Matrix, which were demonstrated to be robust and accurate for non-radiative applications. The robustness and the limitations of the new algorithms are tested by solving a collection of one-dimensional and multi-dimensional test problems, both in the free-streaming and in the diffusion radiation transport limits. Due to its stable performance, the applicability of the scheme presented in this paper to real astrophysical scenarios in high-energy astrophysics is promising. In future simulations, we expect to be able to explore the dynamical relevance of photon-matter interactions in the context of relativistic jets and accretion discs, from microquasars and AGN to gamma-ray bursts.
21 pages, 13 figures. Accepted for publication in Computer Physics Communications
References in corpus (19)
- ECHO: an Eulerian Conservative High Order scheme for general relativistic magnetohydrodynamics and magnetodynamics
- The Nature of SS433 and the Ultraluminous X-ray Sources
- Semi-implicit scheme for treating radiation under M1 closure in general relativistic conservative fluid dynamics codes
- Post-merger evolution of a neutron star-black hole binary with neutrino transport
- An HLLC Solver for Relativistic Flows -- II. Magnetohydrodynamics
- Three-Dimensional Supernova Explosion Simulations of 9-, 10-, 11-, 12-, and 13-M Stars
- Relativistic Radiation Magnetohydrodynamics in Dynamical Spacetimes: Numerical Methods and Tests
- A powerful hydrodynamic booster for relativistic jets
- General-relativistic neutrino-radiation magnetohydrodynamic simulation of seconds-long black hole-neutron star mergers
- Numerical Simulations of Optically Thick Accretion onto a Black Hole - II. Rotating Flow
- A radiative transfer module for relativistic magnetohydrodynamics in the PLUTO code
- Two-moment scheme for general-relativistic radiation hydrodynamics: a systematic description and new applications
- Equations of General Relativistic Radiation Hydrodynamics from Tensor Formalism
- 3D RMHD simulations of jet-wind interactions in High Mass X-ray Binaries
- CAFE-R a code that solves the special relativistic radiation hydrodynamics equations
- On the correction of conserved variables for numerical RMHD with staggered constrained transport
- Simulations of radiation driven winds from Keplerian discs
- Classifying initial conditions of long GRBs modeled with relativistic radiation hydrodynamics
- Future Simulations of Tidal Disruption Events