An Implicit Finite Volume Scheme to Solve the Timedependent Radiation Transport Equation Based on Discrete Ordinates
Abstract
We describe a new algorithm to implicitly solve the timedependent, frequencyintegrated radiation transport (RT) equation, which is coupled to an explicit solver for equations of magnetohydrodynamics (MHD) using Athena++. The radiation field 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 terms for spatial transport of photons and interactions between gas and radiation are calculated implicitly together. An efficient Jacobilike iteration scheme is used to solve the implicit equations. This removes any timestep constraint 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 comoving 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 CourantFriedrichsLewy 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 pressuredominated flows to demonstrate its accuracy and efficiency.
 Publication:

The Astrophysical Journal Supplement Series
 Pub Date:
 April 2021
 DOI:
 10.3847/15384365/abe303
 arXiv:
 arXiv:2102.02212
 Bibcode:
 2021ApJS..253...49J
 Keywords:

 Computational methods;
 Radiative transfer;
 1965;
 1335;
 Astrophysics  Instrumentation and Methods for Astrophysics
 EPrint:
 24 pages, 16 Figures, accepted by ApJS