Residual Monte Carlo high-order solver for Moment-Based Accelerated Thermal Radiative Transfer equations

Journal Article · Sat Nov 01 00:00:00 EDT 2014 · Journal of Computational Physics

DOI:https://doi.org/10.1016/J.JCP.2014.07.039· OSTI ID:22382136

Willert, Jeffrey; Park, H.

In this article we explore the possibility of replacing Standard Monte Carlo (SMC) transport sweeps within a Moment-Based Accelerated Thermal Radiative Transfer (TRT) algorithm with a Residual Monte Carlo (RMC) formulation. Previous Moment-Based Accelerated TRT implementations have encountered trouble when stochastic noise from SMC transport sweeps accumulates over several iterations and pollutes the low-order system. With RMC we hope to significantly lower the build-up of statistical error at a much lower cost. First, we display encouraging results for a zero-dimensional test problem. Then, we demonstrate that we can achieve a lower degree of error in two one-dimensional test problems by employing an RMC transport sweep with multiple orders of magnitude fewer particles per sweep. We find that by reformulating the high-order problem, we can compute more accurate solutions at a fraction of the cost.

• OSTI ID: 22382136
• Journal Information: Journal of Computational Physics, Vol. 276; Other Information: Copyright (c) 2014 Elsevier Science B.V., Amsterdam, The Netherlands, All rights reserved.; Country of input: International Atomic Energy Agency (IAEA); ISSN 0021-9991
• Country of Publication: United States
• Language: English

Related Subjects

71 CLASSICAL AND QUANTUM MECHANICS
GENERAL PHYSICS
ACCELERATION
ALGORITHMS
DIFFERENTIAL EQUATIONS
IMPLEMENTATION
MATHEMATICAL SOLUTIONS
MONTE CARLO METHOD
NOISE
ONE-DIMENSIONAL CALCULATIONS
RADIANT HEAT TRANSFER
STOCHASTIC PROCESSES
THERMAL CONDUCTION