skip to main content

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

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.
 [1] ;  [1] ;  [2]
  1. ORNL
  2. University of Wisconsin, Madison
Publication Date:
OSTI Identifier:
DOE Contract Number:
Resource Type:
Journal Article
Resource Relation:
Journal Name: Journal of Computational Physics; Journal Volume: 258
Research Org:
Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States)
Sponsoring Org:
Work for Others (WFO)
Country of Publication:
United States
radiation diffusion; synthetic acceleration; Monte Carlo; sparse matrix systems