𝔖 Bobbio Scriptorium
✦   LIBER   ✦

Semi-implicit time integration for thermal radiative transfer

✍ Scribed by Ryan G. McClarren; Thomas M. Evans; Robert B. Lowrie; Jeffery D. Densmore


Publisher
Elsevier Science
Year
2008
Tongue
English
Weight
406 KB
Volume
227
Category
Article
ISSN
0021-9991

No coin nor oath required. For personal study only.

✦ Synopsis


Implicit time integration involving the solution of large systems of equations is the current paradigm for time-dependent radiative transfer. In this paper we present a semi-implicit, linear discontinuous Galerkin method for the spherical harmonics ðP N Þ equations for thermal radiative transfer in planar geometry. Our method is novel in that the material coupling terms are treated implicitly (via linearizing the emission source) and the streaming operator is treated explicitly using a second-order accurate Runge-Kutta method. The benefit of this approach is that each time step only involves the solution of equations that are local to each cell. This benefit comes at the cost of having the time step limited by a CFL condition based on the speed of light. To guarantee positivity and avoid artificial oscillations, we use a slope-limiting technique. We present analysis and numerical results that show the method is robust in the diffusion limit when the photon mean-free path is not resolved by the spatial mesh. Also, in the diffusion limit the time step restriction relaxes to a less restrictive explicit diffusion CFL condition. We demonstrate with numerical results that away from the diffusion limit our method demonstrates second-order error convergence as the spatial mesh is refined with a fixed CFL number.


πŸ“œ SIMILAR VOLUMES


Implicit Monte Carlo Diffusionβ€”An Accele
✍ N.A. Gentile πŸ“‚ Article πŸ“… 2001 πŸ› Elsevier Science 🌐 English βš– 189 KB

Implicit Monte Carlo (IMC) is often employed to numerically simulate radiative transfer. In problems with regions that are characterized by a small mean free path, IMC can take a prohibitive amount of time, because many particle steps must be simulated to advance the particle through the time step.