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

JOURNAL OF COMPUTATIONAL PHYSICS

卷:276

Article

Willert, Jeffrey¹  Park, H.¹ 

[1] Los Alamos Natl Lab, Div Theoret, MS B216, Los Alamos, NM 87545 USA

关键词: Thermal Radiative Transfer;    Moment-Based Acceleration;    Monte Carlo;    Residual Monte Carlo;   

DOI  :  10.1016/j.jcp.2014.07.039

来源: Elsevier

【 摘 要 】

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. (C) 2014 Elsevier Inc. All rights reserved.

【 授权许可】

Free