Multi-Frequency General Relativistic Radiation-Hydrodynamics with ${\mathbf{M}}_1$ Closure

Abstract

We report on recent upgrades to our general relativistic radiation-magnetohydrodynamics code, Cosmos++, which expands the two-moment, ${\mathbf{M}}_1$, radiation treatment from grey to multi-frequency transport, including Doppler and gravitational frequency shifts. The solver accommodates either photon (Bose-Einstein) or neutrino (Fermi-Dirac) statistical distribution functions with absorption, emission, and elastic scattering processes. An implicit scheme is implemented to simultaneously solve the primitive inversion problem together with the radiation-matter coupling source terms, providing stability over a broad range of opacities and optical depths where the interactions terms can be stiff. We discuss our formulations and numerical methods, and validate our methods against a wide variety of test problems spanning optically thin to thick regimes in flat, weakly curved, and strongly curved spacetimes.

22 pages, 11 figures, submitted to Astrophysical Journal

References in corpus (13)

• Multi-messenger Observations of a Binary Neutron Star Merger
• The Zwicky Transient Facility: Science Objectives
• Truncated Moment Formalism for Radiation Hydrodynamics in Numerical Relativity
• 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
• Impact of an improved neutrino energy estimate on outflows in neutron star merger simulations
• The explosion mechanism of core-collapse supernovae: progress in supernova theory and experiments
• Relativistic Radiation Magnetohydrodynamics in Dynamical Spacetimes: Numerical Methods and Tests
• High-Energy Neutrino Flares From X-Ray Bright and Dark Tidal Disruptions Events
• Numerical Simulations of Optically Thick Accretion onto a Black Hole - II. Rotating Flow
• Multigroup radiation hydrodynamics with flux-limited diffusion and adaptive mesh refinement
• CosmosDG: An hp-Adaptive Discontinuous Galerkin Code for Hyper-Resolved Relativistic MHD
• Covariant Radiative Transfer for Black Hole Spacetimes

Cited by in corpus (5)

• An Implicit Finite Volume Scheme to Solve the Time Dependent Radiation Transport Equation Based on Discrete Ordinates
• An Extension of the Athena++ Code Framework for Radiation-Magnetohydrodynamics in General Relativity Using a Finite-Solid-Angle Discretization
• The Guided Moments formalism: a new efficient full-neutrino treatment for astrophysical simulations
• Thermal aspects of neutron star mergers
• Incorporating a radiative hydrodynamics scheme in the numerical-relativity code BAM