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.

An asymptotically correct implicit–explicit time integration scheme for finite volume radiation-hydrodynamics

An Asymptotically-Correct Implicit-Explicit Time Integration Scheme for Finite Volume Radiation-Hydrodynamics

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

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

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

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

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

A Spatial-Temporal Asymptotic Preserving Scheme for Radiation Magnetohydrodynamics in the Equilibrium and Non-Equilibrium Diffusion Limit

Implicit-Explicit Finite-Difference Lattice Boltzmann Method For Compressible Flows

Scalable explicit implementation of anisotropic diffusion with Runge-Kutta-Legendre super-time-stepping

A conservative implicit scheme for three-dimensional steady flows of diatomic gases in all flow regimes using unstructured meshes in the physical and velocity spaces

Semi-implicit Front Capturing Schemes for the Degenerate Nonlinear Radiative Diffusion Equation.

A direct Eulerian GRP scheme for radiation hydrodynamical equations in diffusion limit

Efficient solver of relativistic hydrodynamics with implicit Runge-Kutta method

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

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

A Fourth-Order Accurate Finite-Volume Method with Structured Adaptive Mesh Refinement for Solving the Advection-Diffusion Equation

A coupled high-order implicit-explicit flux reconstruction lattice Boltzmann method for nearly incompressible thermal flows

Implicit high-order gas-kinetic schemes for compressible flows on three-dimensional unstructured meshes I: Steady flows

High order implicit time integration schemes on multiresolution adaptive grids for stiff PDEs