Abstract:Numerical radiation-hydrodynamics (RHD) for non-relativistic flows is a challenging problem because it encompasses processes acting over a very broad range of timescales, and where the relative importance of these processes often varies by orders of magnitude across the computational domain. Here we present a new implicit-explicit (IMEX) method for numerical RHD that has a number of desirable properties that have not previously been combined in a single method. Our scheme is based on moments and allows machine-precision conservation of energy and momentum, making it highly suitable for adaptive mesh refinement applications; it requires no more communication than hydrodynamics and includes no non-local iterative steps, making it highly suitable for massively parallel and GPU-based systems where communication is a bottleneck; and we show that it is asymptotically-accurate in the streaming, static diffusion, and dynamic diffusion limits, including in the so-called asymptotic diffusion regime where the computational grid does not resolve the photon mean free path. We implement our method in the GPU-accelerated RHD code QUOKKA and show that it passes a wide range of numerical tests.

What problem does this paper attempt to address?

The problem that this paper attempts to solve is: in non - relativistic fluids, numerical radiation - hydrodynamics (RHD) problems become very complicated due to processes involving multiple time scales. Specifically, the paper aims to develop a new implicit - explicit (IMEX) time - integration method to address the following key challenges: 1. **Multi - time - scale problem**: The time scales of radiation and hydrodynamics processes differ greatly, where the radiation process is determined by the speed of light \( c \) and the gas flow velocity \( v \). This large difference in time scales makes it difficult for traditional numerical methods to handle both efficiently and accurately at the same time. 2. **Optical depth variation**: Actual RHD problems usually cover a wide range of opacity, so that the mean free path of photons in some regions may be comparable to the size of the entire simulation domain, while in other regions it may be so small that it cannot be resolved at a reasonable computational cost. 3. **Non - local implicit update**: To ensure stability, many RHD codes use implicit methods to handle the transport term, but this requires non - local implicit updates, which will lead to communication bottlenecks in modern large - scale parallel and GPU - accelerated architectures, thus affecting performance. 4. **Accuracy of asymptotic limits**: Existing operator - splitting methods often have difficulty maintaining high accuracy when dealing with different asymptotic limits (such as static diffusion, dynamic diffusion and fluid limits), especially in terms of the balance between these limits. To solve these problems, the author proposes a new IMEX method, which has the following advantages: - **Energy and momentum conservation at machine precision**: The moment - based method can achieve highly accurate energy and momentum conservation. - **Efficient parallelism and GPU compatibility**: No additional cross - domain communication is required, and there are no non - local iteration steps, which is very suitable for large - scale parallel and GPU systems. - **Asymptotic accuracy**: This method can maintain asymptotic accuracy in the fluid limit, static diffusion limit and dynamic diffusion limit, including in the so - called asymptotic diffusion region, even if the computational grid cannot resolve the mean free path of photons. Through these improvements, this method can maintain high accuracy in all RHD limits while avoiding the communication bottlenecks in traditional methods, thus improving computational efficiency and scalability.

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

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

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

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

A structure-preserving semi-implicit IMEX finite volume scheme for ideal magnetohydrodynamics at all Mach and Alfvén numbers

One-sweep moment-based semi-implicit-explicit integration for gray thermal radiation transport

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

Efficient solver of relativistic hydrodynamics with implicit Runge-Kutta method

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

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

A Spectral Radius Scaling Semi-Implicit Iterative Time Stepping Method for Reactive Flow Simulations with Detailed Chemistry

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

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

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

Fast Iterative Method with a Second-Order Implicit Difference Scheme for Time-Space Fractional Convection–Diffusion Equation

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

A New Discrete Implicit Monte Carlo Scheme for Simulating Radiative Transfer Problems