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

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.

arXiv comment14 pages, 6 figures, submitted to MNRAS; Code available for download at https://github.com/quokka-astro/quokka

References in corpus (25)

• Radiation Pressure Supported Starburst Disks and AGN Fueling
• Theoretical Challenges in Galaxy Formation
• The Formation of Massive Star Systems by Accretion
• Multi--dimensional Cosmological Radiative Transfer with a Variable Eddington Tensor Formalism
• Super-Eddington Accretion Disks around Supermassive black Holes
• A scheme for radiation pressure and photon diffusion with the M1 closure in RAMSES-RT
• Fornax: a Flexible Code for Multiphysics Astrophysical Simulations
• Equations and Algorithms for Mixed Frame Flux-Limited Diffusion Radiation Hydrodynamics
• Magnetohydrodynamic Simulations of Active Galactic Nucleus Disks and Jets
• An Unstable Truth: How Massive Stars get their Mass
• A Godunov Method for Multidimensional Radiation Magnetohydrodynamics based on a variable Eddington tensor
• A Two-moment Radiation Hydrodynamics Module in Athena Using a Time-explicit Godunov Method
• An Algorithm for Radiation Magnetohydrodynamics Based on Solving the Time-dependent Transfer Equation
• Atmospheric Regimes and Trends on Exoplanets and Brown Dwarfs
• CASTRO: A New Compressible Astrophysical Solver. II. Gray Radiation Hydrodynamics
• Should One Use the Ray-by-Ray Approximation in Core-Collapse Supernova Simulations?
• Turbulence in Core-Collapse Supernovae
• Simulating Star Clusters Across Cosmic Time: I. Initial Mass Function, Star Formation Rates and Efficiencies
• Outflows Driven by Direct and Reprocessed Radiation Pressure in Massive Star Clusters
• Saturation of the MRI in Strongly Radiation Dominated Accretion Disks
• Quokka: A code for two-moment AMR radiation hydrodynamics on GPUs
• VETTAM: A scheme for radiation hydrodynamics with adaptive mesh refinement using the variable Eddington tensor method
• Realizability-Preserving DG-IMEX Method for the Two-Moment Model of Fermion Transport
• A radiation-hydrodynamics scheme valid from the transport to the diffusion limit
• Time Dependent Radiation Hydrodynamics on a Moving Mesh

Cited by in corpus (4)

• Quokka-based understanding of outflows (QED) -- II. X-ray metallicity gradients as a signature of galactic wind metal loading
• A novel numerical method for mixed-frame multigroup radiation-hydrodynamics with GPU acceleration implemented in the QUOKKA code
• QED V: Variations in metal loading of galactic winds with element nucleosynthetic origin
• Quokka-based understanding of outflows (QED) -- IV. Limitations of H$\alpha$ as an outflow diagnostic