Abstract:Mixed-frame formulations of radiation-hydrodynamics (RHD), where the radiation quantities are computed in an inertial frame but matter quantities are in a comoving frame, are advantageous because they admit algorithms that conserve energy and momentum to machine precision and combine more naturally with adaptive mesh techniques, since unlike pure comoving-frame methods they do not face the problem that radiation quantities must change frame every time a cell is refined or coarsened. However, implementing multigroup RHD in a mixed-frame formulation presents challenges due to the complexity of handling frequency-dependent interactions and the Doppler shift of radiation boundaries. In this paper, we introduce a novel method for multigroup RHD that integrates a mixed-frame formulation with a piecewise powerlaw approximation for frequency dependence within groups. This approach ensures the exact conservation of total energy and momentum while effectively managing the Lorentz transformation of group boundaries and evaluation of group-averaged opacities. Our method takes advantage of the locality of matter-radiation coupling, allowing the source term for $N_g$ frequency groups to be handled with simple equations with a sparse Jacobian matrix of size $N_g + 1$, which can be inverted with $O(N_g)$ complexity. This results in a computational complexity that scales linearly with $N_g$ and requires no more communication than a pure hydrodynamics update, making it highly efficient for massively parallel and GPU-based systems. We implement our method in the GPU-accelerated RHD code QUOKKA and demonstrate that it passes a wide range of numerical tests. We demonstrate that the piecewise powerlaw method shows significant advantages over traditional opacity averaging methods for handling rapidly variable opacities with modest frequency resolution.

Multi-group Radiation Magneto-hydrodynamics based on Discrete Ordinates including Compton Scattering

A numerical model for multigroup radiation hydrodynamics

An Implicit Finite Volume Scheme to Solve the Time Dependent Radiation Transport Equation Based on Discrete Ordinates

A Unified Gas-Kinetic Particle Method for Frequency-Dependent Radiative Transfer Equations with Isotropic Scattering Process on Unstructured Mesh

Multiscale Simulation for the System of Radiation Hydrodynamics

Fervent: Chemistry-coupled, ionising and non-ionising radiative feedback in magnetohydrodynamical simulations

Multiscale Radiative Transfer in Cylindrical Coordinates

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

A general hybrid radiation transport scheme for star formation simulations on an adaptive grid

An Adaptive High-Order Transient Algorithm to Solve Large-Scale Anisotropic Maxwell’s Equations

A novel numerical method for mixed-frame multigroup radiation-hydrodynamics with GPU acceleration implemented in the QUOKKA code

Unified Gas-Kinetic Particle Method for Frequency-dependent Radiation Transport

An Implicit Unified Gas-Kinetic Particle Method with Large Time Steps for Gray Radiation Transport

Unified gas-kinetic simulation for the system of multiscale radiation hydrodynamics

A fully asymptotic preserving decomposed multi-group method for the frequency-dependent radiative transfer equations

Frequency-dependent Discrete Implicit Monte-Carlo Scheme for the Radiative Transfer Equation

An Extension of the Athena++ Code Framework for Radiation-magnetohydrodynamics in General Relativity Using a Finite-solid-angle Discretization

General-relativistic radiation transport scheme in $\texttt{Gmunu}$ I: Implementation of two-moment based multi-frequency radiative transfer and code tests

A radiation hydrodynamics scheme on adaptive meshes using the Variable Eddington Tensor (VET) closure

A new MHD code with adaptive mesh refinement and parallelization for astrophysics

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