A radiative transfer module for relativistic magnetohydrodynamics in the PLUTO code
arXiv:1903.10456 · doi:10.3847/1538-4365/ab18ff
Abstract
We present a numerical implementation for the solution of the relativistic radiation hydrodynamics and magnetohydrodynamics equations, designed as an independent module within the freely available code PLUTO. The radiation transfer equations are solved under the grey approximation and imposing the M1 closure, which allows the radiation transport to be handled in both the free-streaming and diffusion limits. Equations are integrated following an implicit-explicit scheme, where radiation-matter interaction terms are integrated implicitly, whereas transport and all of the remaining source terms are solved explicitly by means of the same Godunov-type solvers included in PLUTO. Among these, we introduce a new Harten-van Leer-contact (HLLC) solver for optically thin radiation transport. The code is suitable for multidimensional computations in Cartesian, spherical and cylindrical coordinates, using either a single processor or parallel architectures. Adaptive grid computations are also made possible, by means of the CHOMBO library. The algorithm performance is demonstrated through a series of numerical benchmarks by investigating various different configurations with a particular emphasis on the behavior of the solutions in the free-streaming and diffusion limits.
25 pages, 17 figures
References in corpus (26)
- PLUTO: a Numerical Code for Computational Astrophysics
- Gamma-Ray Bursts
- The PLUTO Code for Adaptive Mesh Computations in Astrophysical Fluid Dynamics
- Three-Dimensional General Relativistic Radiation Magnetohydrodynamical Simulation of Super-Eddington Accretion, using a new code HARMRAD with M1 Closure
- On the Divergence-Free Condition in Godunov-Type Schemes for Ideal Magnetohydrodynamics: the Upwind Constrained Transport Method
- An HLLC Riemann Solver for Relativistic Flows: I. Hydrodynamics
- Time-dependent 3D spectrum synthesis for type Ia supernovae
- Equation of State in Relativistic Magnetohydrodynamics: variable versus constant adiabatic index
- 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
- An HLLC Solver for Relativistic Flows -- II. Magnetohydrodynamics
- Radiation hydrodynamics with Adaptive Mesh Refinement and application to prestellar core collapse. I Methods
- High-order conservative finite difference GLM-MHD schemes for cell-centered MHD
- High-order conservative reconstruction schemes for finite volume methods in cylindrical and spherical coordinates
- Spectral evolution of superluminal components in parsec-scale jets
- General relativistic radiation hydrodynamics of accretion flows. I: Bondi-Hoyle accretion
- Beyond Flux-Limited Diffusion: Parallel Algorithms for Multidimensional Radiation Hydrodynamics
- A Two-moment Radiation Hydrodynamics Module in Athena Using a Time-explicit Godunov Method
- Relativistic Radiation Magnetohydrodynamics in Dynamical Spacetimes: Numerical Methods and Tests
- Numerical Simulations of Optically Thick Accretion onto a Black Hole - I. Spherical Case
- General relativistic radiation hydrodynamics of accretion flows: II. Treating stiff source terms and exploring physical limitations
- Numerical Treatment of Anisotropic Radiation Field Coupling with the Relativistic Resistive Magnetofluids
- Light Curve Calculations of Supernovae from Fallback Gamma-Ray Bursts
- Interface Conditions for Wave Propagation Through Mesh Refinement Boundaries
- Modeling long GRBs using a single shock with relativistic radiation hydrodynamics
Cited by in corpus (21)
- General relativistic resistive magnetohydrodynamics with robust primitive variable recovery for accretion disk simulations
- Beyond moments: relativistic Lattice-Boltzmann methods for radiative transport in computational astrophysics
- Two-moment scheme for general-relativistic radiation hydrodynamics: a systematic description and new applications
- A two-moment radiation hydrodynamics scheme applicable to simulations of planet formation in circumstellar disks
- Successive interacting coronal mass ejections: How to create a perfect storm?
- Three-temperature radiation hydrodynamics with PLUTO: Tests and applications to protoplanetary disks
- General-relativistic radiation transport scheme in I: Implementation of two-moment based multi-frequency radiative transfer and code tests
- No self-shadowing instability in 2D radiation-hydrodynamical models of irradiated protoplanetary disks
- Vertical shear instability in two-moment radiation-hydrodynamical simulations of irradiated protoplanetary disks I. Angular momentum transport and turbulent heating
- The Guided Moments formalism: a new efficient full-neutrino treatment for astrophysical simulations
- Vertical shear instability in two-moment radiation-hydrodynamical simulations of irradiated protoplanetary disks II. Secondary instabilities and stability regions
- Radiation transport methods in star formation simulations
- Radiative shock oscillation model for the long-term flares of Sgr A*
- Three-temperature radiation hydrodynamics with PLUTO: Thermal and kinematic signatures of accreting protoplanets
- Non-LTE radiation hydrodynamics in PLUTO
- Can gap-edge illumination excite spirals in protoplanetary disks? Three-temperature radiation hydrodynamics and NIR image modelling
- A new ionisation network and radiation transport module in PLUTO
- 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
- Time delays between radio and X-ray and between narrow radio bands of Sgr A* flares in the shock oscillation model
- MeV, GeV and TeV neutrinos from binary-driven hypernovae