A Monte Carlo Synthetic-Acceleration Method for Solving the Thermal Radiation Diffusion Equation

Journal Article · Wed Jan 01 00:00:00 EST 2014 · Journal of Computational Physics

DOI:https://doi.org/10.1016/j.jcp.2013.10.043· OSTI ID:1105918

Evans, Thomas M ^[1]; Mosher, Scott W ^[1]; Slattery, Stuart ^[2]

---

1. ORNL
2. University of Wisconsin, Madison

We present a novel synthetic-acceleration based Monte Carlo method for solving the equilibrium thermal radiation diusion equation in three dimensions. The algorithm performance is compared against traditional solution techniques using a Marshak benchmark problem and a more complex multiple material problem. Our results show that not only can our Monte Carlo method be an eective solver for sparse matrix systems, but also that it performs competitively with deterministic methods including preconditioned Conjugate Gradient while producing numerically identical results. We also discuss various aspects of preconditioning the method and its general applicability to broader classes of problems.

• Research Organization: Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States)
• Sponsoring Organization: Work for Others (WFO)
• DOE Contract Number: DE-AC05-00OR22725
• OSTI ID: 1105918
• Journal Information: Journal of Computational Physics, Vol. 258; ISSN 0021-9991
• Country of Publication: United States
• Language: English

Related Subjects

radiation diffusion
synthetic acceleration
Monte Carlo
sparse matrix systems