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

Abstract

We describe a new algorithm to solve the time dependent, frequency integrated radiation transport (RT) equation implicitly, which is coupled to an explicit solver for equations of magnetohydrodynamics (MHD) using {\sf Athena++}. The radiation filed is represented by specific intensities along discrete rays, which are evolved using a conservative finite volume approach for both cartesian and curvilinear coordinate systems. All the terms for spatial transport of photons and interactions between gas and radiation are calculated implicitly together. An efficient Jacobi-like iteration scheme is used to solve the implicit equations. This removes any time step constrain due to the speed of light in RT. We evolve the specific intensities in the lab frame to simplify the transport step. The lab-frame specific intensities are transformed to the co-moving frame via Lorentz transformation when the source term is calculated. Therefore, the scheme does not need any expansion in terms of $v/c$. The radiation energy and momentum source terms for the gas are calculated via direct quadrature in the angular space. The time step for the whole scheme is determined by the normal Courant -- Friedrichs -- Lewy condition in the MHD module. We provide a variety of test problems for this algorithm including both optically thick and thin regimes, and for both gas and radiation pressure dominated flows to demonstrate its accuracy and efficiency.

24 pages, 16 Figures, accepted by ApJS

References in corpus (13)

• The Athena++ Adaptive Mesh Refinement Framework: Design and Magnetohydrodynamic Solvers
• Semi-implicit scheme for treating radiation under M1 closure in general relativistic conservative fluid dynamics codes
• Equations and Algorithms for Mixed Frame Flux-Limited Diffusion Radiation Hydrodynamics
• An Algorithm for Radiation Magnetohydrodynamics Based on Solving the Time-dependent Transfer Equation
• Three-dimensional Boltzmann-Hydro code for core-collapse in massive stars I. special relativistic treatments
• A New Spherical Harmonics Scheme for Multi-Dimensional Radiation Transport I: Static Matter Configurations
• Radiation Feedback in ULIRGS: Are Photons Movers and Shakers?
• Monte Carlo Radiation Hydrodynamics with Implicit Methods
• Monte-Carlo closure for moment-based transport schemes in general relativistic radiation hydrodynamics simulations
• Saturation of the MRI in Strongly Radiation Dominated Accretion Disks
• Multi-Frequency General Relativistic Radiation-Hydrodynamics with ${\mathbf{M}}_1$ Closure
• Covariant Radiative Transfer for Black Hole Spacetimes
• Numerical non-LTE 3D radiative transfer using a multigrid method

Cited by in corpus (4)

• Global 3D Radiation Hydrodynamic Simulations of Proto-Jupiter's Convective Envelope
• VETTAM: A scheme for radiation hydrodynamics with adaptive mesh refinement using the variable Eddington tensor method
• No self-shadowing instability in 2D radiation-hydrodynamical models of irradiated protoplanetary disks
• Thermal instability in radiation hydrodynamics: instability mechanisms, position-dependent S-curves, and attenuation curves