Suppressing nonphysical overheating with a modified implicit Monte Carlo method for time-dependent radiative transfer
- Los Alamos National Laboratory
In this note 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 alleviates both the nonphysical overheating that occurs in standard IMC when the time step is large and significantly diminishes the statistical noise in the solution.
- Research Organization:
- Los Alamos National Laboratory (LANL), Los Alamos, NM (United States)
- Sponsoring Organization:
- USDOE
- DOE Contract Number:
- AC52-06NA25396
- OSTI ID:
- 956608
- Report Number(s):
- LA-UR-08-07561; LA-UR-08-7561; JCTPAH; TRN: US201016%%2293
- Journal Information:
- Journal of Computational Physics, Journal Name: Journal of Computational Physics; ISSN 0021-9991
- Country of Publication:
- United States
- Language:
- English
Similar Records
An adaptive alpha for the implicit Monte Carlo equations
Four decades of implicit Monte Carlo