2001/02/08 by N. J. Turner, J. M. Stone · 3 citations
Engineering · Mathematics · Physics and Astronomy · #Computational Fluid Dynamics and Aerodynamics #Diffusion #Gas Dynamics and Kinetic Theory #Heavy traffic approximation #Laser-Plasma Interactions and Diagnostics #Opacity #Radiant energy #Radiation #Radiation flux #Radiative transfer #Shock (circulatory) #Source function #astro-ph
paper · pdf · doi:10.1086/321779
published as Astrophys.J.Suppl. 135 (2001) 95-108 · Accepted by ApJS; 30 pages, 11 figures
arxiv created 2001/02/08 · openalex publication_date 2001/07/01 · arxiv updated 2009/12/01 · openalex created_date 2016/06/24 · openalex updated_date 2026/08/05
A module for the ZEUS-2D code is described that may be used to solve the equations of radiation hydrodynamics to order unity in v / c , in the flux-limited diffusion (FLD) approximation. In this approximation, the factor Eddington tensor f , which closes the radiation moment equations, is chosen to be an empirical function of the radiation energy density. This is easier to implement and faster than full-transport techniques, in which f is computed by solving the transfer equation. However, FLD is less accurate when the flux has a component perpendicular to the gradient in radiation energy density and in optically thin regions when the radiation field depends strongly on angle. The material component of the fluid is here assumed to be in local thermodynamic equilibrium. The energy equations are operator split, with transport terms, radiation diffusion term, and other source terms evolved separately. Transport terms are applied using the same consistent transport algorithm as in ZEUS-2D. The radiation diffusion term is updated using an alternating direction-implicit method with convergence checking. Remaining source terms are advanced together implicitly using numerical root finding. However, when absorption opacity is zero, accuracy is improved by instead treating the compression and expansion source terms using a time-centered differencing scheme. Results are discussed for test problems including radiation-damped linear waves, radiation fronts propagating in optically thin media, subcritical and supercritical radiating shocks, and an optically thick shock in which radiation dominates downstream pressure.