A High-Order Low-Order Algorithm with Exponentially Convergent Monte Carlo for Thermal Radiative Transfer

Bolding, Simon R.; Cleveland, Mathew A.; Morel, Jim E.

doi:10.13182/nse16-36

Title: A High-Order Low-Order Algorithm with Exponentially Convergent Monte Carlo for Thermal Radiative Transfer

Journal Article · Wed Feb 08 00:00:00 EST 2017 · Nuclear Science and Engineering

DOI:https://doi.org/10.13182/nse16-36· OSTI ID:1673359

Bolding, Simon R. ^[1]; Cleveland, Mathew A. ^[2]; Morel, Jim E. ^[1]

Texas A & M Univ., College Station, TX (United States)
Los Alamos National Lab. (LANL), Los Alamos, NM (United States)

In this paper, we have implemented a new high-order low-order (HOLO) algorithm for solving thermal radiative transfer problems. The low-order (LO) system is based on the spatial and angular moments of the transport equation and a linear-discontinuous finite-element spatial representation, producing equations similar to the standard S₂ equations. The LO solver is fully implicit in time and efficiently resolves the nonlinear temperature dependence at each time step. The high-order (HO) solver utilizes exponentially convergent Monte Carlo (ECMC) to give a globally accurate solution for the angular intensity to a fixed-source pure-absorber transport problem. This global solution is used to compute consistency terms, which require the HO and LO solutions to converge toward the same solution. The use of ECMC allows for the efficient reduction of statistical noise in the Monte Carlo solution, reducing inaccuracies introduced through the LO consistency terms. Finally, we compare results with an implicit Monte Carlo code for one-dimensional gray test problems and demonstrate the efficiency of ECMC over standard Monte Carlo in this HOLO algorithm.

View Accepted Manuscript (DOE)

Cite

Export

Save

Research Organization:: Los Alamos National Laboratory (LANL), Los Alamos, NM (United States)

Sponsoring Organization:: USDOE Office of Nuclear Energy (NE). Nuclear Energy University Programs (NEUP); USDOE National Nuclear Security Administration (NNSA)

Grant/Contract Number:: 89233218CNA000001; AC52-06NA25396; NA0002376

OSTI ID:: 1673359

Alternate ID(s):: OSTI ID: 1338753

Report Number(s):: LA-UR-14-29653; LA-UR-16-20497

Journal Information:: Nuclear Science and Engineering, Vol. 185, Issue 1; ISSN 0029-5639

Publisher:: American Nuclear Society - Taylor & FrancisCopyright Statement

Country of Publication:: United States

Language:: English

Citation Metrics:

Cited by: 1 work

Citation information provided by
Web of Science

References (12)

Asymptotic analysis of radiative transfer problems Larsen, E. W.; Pomraning, G. C.; Badham, V. C. Journal of Quantitative Spectroscopy and Radiative Transfer, Vol. 29, Issue 4 https://doi.org/10.1016/0022-4073(83)90048-1	journal	April 1983
A Linear-Discontinuous Spatial Differencing Scheme forSnRadiative Transfer Calculations Morel, J. E.; Wareing, Todd A.; Smith, Kenneth Journal of Computational Physics, Vol. 128, Issue 2 https://doi.org/10.1006/jcph.1996.0223	journal	October 1996
Properties of the implicitly time-differenced equations of thermal radiation transport Larsen, Edward W.; Kumar, Akansha; Morel, Jim E. Journal of Computational Physics, Vol. 238 https://doi.org/10.1016/j.jcp.2012.11.034	journal	April 2013
Finite-Difference Approximations and Superconvergence for the Discrete-Ordinate Equations in Slab Geometry Larsen, Edward W.; Nelson, Paul SIAM Journal on Numerical Analysis, Vol. 19, Issue 2 https://doi.org/10.1137/0719020	journal	April 1982
Comparison of implicit and symbolic implicit Monte Carlo line transport with frequency weight vector extension McKinley, Michael Scott; Brooks III, Eugene D.; Szoke, Abraham Journal of Computational Physics, Vol. 189, Issue 1 https://doi.org/10.1016/S0021-9991(03)00213-4	journal	July 2003
Asymptotic equilibrium diffusion analysis of time-dependent Monte Carlo methods for grey radiative transfer Densmore, Jeffery D.; Larsen, Edward W. Journal of Computational Physics, Vol. 199, Issue 1 https://doi.org/10.1016/j.jcp.2004.02.004	journal	September 2004
Residual Monte Carlo high-order solver for Moment-Based Accelerated Thermal Radiative Transfer equations Willert, Jeffrey; Park, H. Journal of Computational Physics, Vol. 276 https://doi.org/10.1016/j.jcp.2014.07.039	journal	November 2014
An implicit Monte Carlo scheme for calculating time and frequency dependent nonlinear radiation transport Fleck, J. A.; Cummings, J. D. Journal of Computational Physics, Vol. 8, Issue 3 https://doi.org/10.1016/0021-9991(71)90015-5	journal	December 1971
A hybrid transport-diffusion Monte Carlo method for frequency-dependent radiative-transfer simulations Densmore, Jeffery D.; Thompson, Kelly G.; Urbatsch, Todd J. Journal of Computational Physics, Vol. 231, Issue 20 https://doi.org/10.1016/j.jcp.2012.06.020	journal	August 2012
A new proof of geometric convergence for general transport problems based on sequential correlated sampling methods Kong, Rong; Spanier, Jerome Journal of Computational Physics, Vol. 227, Issue 23 https://doi.org/10.1016/j.jcp.2008.07.016	journal	December 2008
Implicit Monte Carlo Diffusion—An Acceleration Method for Monte Carlo Time-Dependent Radiative Transfer Simulations Gentile, N. A. Journal of Computational Physics, Vol. 172, Issue 2 https://doi.org/10.1006/jcph.2001.6836	journal	September 2001
A Discrete Maximum Principle for the Implicit Monte Carlo Equations Wollaber, Allan B.; Larsen, Edward W.; Densmore, Jeffery D. Nuclear Science and Engineering, Vol. 173, Issue 3 https://doi.org/10.13182/NSE11-101	journal	March 2013

Cited By (1)

Physical, numerical, and computational challenges of modeling neutrino transport in core-collapse supernovae Mezzacappa, Anthony; Endeve, Eirik; Messer, O. E. Bronson arXiv https://doi.org/10.48550/arxiv.2010.09013	preprint	January 2020