A modified implicit Monte Carlo method for time-dependent radiative transfer with adaptive material coupling

Journal Article · Tue Sep 01 04:00:00 EDT 2009 · Journal of Computational Physics

DOI:https://doi.org/10.1016/j.jcp.2009.04.028· OSTI ID:21308105

In this paper we develop a robust implicit Monte Carlo (IMC) algorithm based on more accurately updating the linearized equilibrium radiation energy density. The method does not introduce oscillations in the solution and has the same limit as {delta}t{yields}{infinity} as the standard Fleck and Cummings IMC method. Moreover, the approach we introduce can be trivially added to current implementations of IMC by changing the definition of the Fleck factor. Using this new method we develop an adaptive scheme that uses either standard IMC or the modified method basing the adaptation on a zero-dimensional problem solved in each cell. Numerical results demonstrate that the new method can avoid the nonphysical overheating that occurs in standard IMC when the time step is large. The method also leads to decreased noise in the material temperature at the cost of a potential increase in the radiation temperature noise.

• OSTI ID: 21308105
• Journal Information: Journal of Computational Physics, Journal Name: Journal of Computational Physics Journal Issue: 16 Vol. 228; ISSN JCTPAH; ISSN 0021-9991
• Country of Publication: United States
• Language: English

Related Subjects

71 CLASSICAL AND QUANTUM MECHANICS
GENERAL PHYSICS
ALGORITHMS
ENERGY DENSITY
EQUILIBRIUM
MATHEMATICAL SOLUTIONS
MONTE CARLO METHOD
OSCILLATIONS
RADIANT HEAT TRANSFER
TEMPERATURE NOISE
TIME DEPENDENCE