Monte Carlo simulation methods in momentbased scalebridging algorithms for thermal radiativetransfer problems
Abstract
We present a momentbased acceleration algorithm applied to Monte Carlo simulation of thermal radiativetransfer problems. Our acceleration algorithm employs a continuum system of moments to accelerate convergence of stiff absorption–emission physics. The combination of energyconserving tallies and the use of an asymptotic approximation in optically thick regions remedy the difficulties of local energy conservation and mitigation of statistical noise in such regions. We demonstrate the efficiency and accuracy of the developed method. We also compare directly to the standard linearizationbased method of Fleck and Cummings [1]. A factor of 40 reduction in total computational time is achieved with the new algorithm for an equivalent (or more accurate) solution as compared with the Fleck–Cummings algorithm.
 Authors:
 Bettis Atomic Power Laboratory, P.O. Box 79, West Mifflin, PA 15122 (United States)
 Fluid Dynamics and Solid Mechanics Group, Los Alamos National Laboratory, P.O. Box 1663, MS B216, Los Alamos, NM 87545 (United States)
 Computational Physics and Methods Group, Los Alamos National Laboratory, P.O. Box 1663, MS D409, Los Alamos, NM 87545 (United States)
 Publication Date:
 OSTI Identifier:
 22465600
 Resource Type:
 Journal Article
 Resource Relation:
 Journal Name: Journal of Computational Physics; Journal Volume: 284; Other Information: Copyright (c) 2014 Elsevier Science B.V., Amsterdam, The Netherlands, All rights reserved.; Country of input: International Atomic Energy Agency (IAEA)
 Country of Publication:
 United States
 Language:
 English
 Subject:
 71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; 97 MATHEMATICAL METHODS AND COMPUTING; ABSORPTION; ACCURACY; ALGORITHMS; ASYMPTOTIC SOLUTIONS; COMPUTERIZED SIMULATION; CONVERGENCE; EFFICIENCY; EMISSION; ENERGY CONSERVATION; MONTE CARLO METHOD; RADIANT HEAT TRANSFER; STATISTICS
Citation Formats
Densmore, J.D., Email: jeffery.densmore@unnpp.gov, Park, H., Email: hkpark@lanl.gov, Wollaber, A.B., Email: wollaber@lanl.gov, Rauenzahn, R.M., Email: rick@lanl.gov, and Knoll, D.A., Email: nol@lanl.gov. Monte Carlo simulation methods in momentbased scalebridging algorithms for thermal radiativetransfer problems. United States: N. p., 2015.
Web. doi:10.1016/J.JCP.2014.12.020.
Densmore, J.D., Email: jeffery.densmore@unnpp.gov, Park, H., Email: hkpark@lanl.gov, Wollaber, A.B., Email: wollaber@lanl.gov, Rauenzahn, R.M., Email: rick@lanl.gov, & Knoll, D.A., Email: nol@lanl.gov. Monte Carlo simulation methods in momentbased scalebridging algorithms for thermal radiativetransfer problems. United States. doi:10.1016/J.JCP.2014.12.020.
Densmore, J.D., Email: jeffery.densmore@unnpp.gov, Park, H., Email: hkpark@lanl.gov, Wollaber, A.B., Email: wollaber@lanl.gov, Rauenzahn, R.M., Email: rick@lanl.gov, and Knoll, D.A., Email: nol@lanl.gov. 2015.
"Monte Carlo simulation methods in momentbased scalebridging algorithms for thermal radiativetransfer problems". United States.
doi:10.1016/J.JCP.2014.12.020.
@article{osti_22465600,
title = {Monte Carlo simulation methods in momentbased scalebridging algorithms for thermal radiativetransfer problems},
author = {Densmore, J.D., Email: jeffery.densmore@unnpp.gov and Park, H., Email: hkpark@lanl.gov and Wollaber, A.B., Email: wollaber@lanl.gov and Rauenzahn, R.M., Email: rick@lanl.gov and Knoll, D.A., Email: nol@lanl.gov},
abstractNote = {We present a momentbased acceleration algorithm applied to Monte Carlo simulation of thermal radiativetransfer problems. Our acceleration algorithm employs a continuum system of moments to accelerate convergence of stiff absorption–emission physics. The combination of energyconserving tallies and the use of an asymptotic approximation in optically thick regions remedy the difficulties of local energy conservation and mitigation of statistical noise in such regions. We demonstrate the efficiency and accuracy of the developed method. We also compare directly to the standard linearizationbased method of Fleck and Cummings [1]. A factor of 40 reduction in total computational time is achieved with the new algorithm for an equivalent (or more accurate) solution as compared with the Fleck–Cummings algorithm.},
doi = {10.1016/J.JCP.2014.12.020},
journal = {Journal of Computational Physics},
number = ,
volume = 284,
place = {United States},
year = 2015,
month = 3
}

Monte Carlo solution methods in a momentbased scalebridging algorithm for thermal radiative transfer problems: Comparison with Fleck and Cummings
We have developed a momentbased scalebridging algorithm for thermal radiative transfer problems. The algorithm takes the form of wellknown nonlineardiffusion acceleration which utilizes a loworder (LO) continuum problem to accelerate the solution of a highorder (HO) kinetic problem. The coupled nonlinear equations that form the LO problem are efficiently solved using a preconditioned Jacobianfree NewtonKrylov method. This work demonstrates the applicability of the scalebridging algorithm with a Monte Carlo HO solver and reports the computational efficiency of the algorithm in comparison to the wellknown FleckCummings algorithm. (authors) 
Residual Monte Carlo highorder solver for MomentBased Accelerated Thermal Radiative Transfer equations
In this article we explore the possibility of replacing Standard Monte Carlo (SMC) transport sweeps within a MomentBased Accelerated Thermal Radiative Transfer (TRT) algorithm with a Residual Monte Carlo (RMC) formulation. Previous MomentBased Accelerated TRT implementations have encountered trouble when stochastic noise from SMC transport sweeps accumulates over several iterations and pollutes the loworder system. With RMC we hope to significantly lower the buildup of statistical error at a much lower cost. First, we display encouraging results for a zerodimensional test problem. Then, we demonstrate that we can achieve a lower degree of error in two onedimensional test problems bymore » 
Vector and parallel Monte Carlo radiative heat transfer simulation
A fully vectorized version of a Monte Carlo algorithm of radiative heat transfer in twodimensional geometries is presented. This algorithm differs from previous applications in that its capabilities are more extensive, with arbitrary numbers of surfaces, arbitrary numbers of material properties, and surface characteristics that include transmission, specular reflection, and diffuse reflection (all of which may be functions of the angle of incidence). The algorithm is applied to an irregular, experimental geometry and implemented on a Cyber 205. A speedup factor of approximately 16, for this combination of geometry and material properties, is achieved for the vector version over themore »