A modified implicit Monte Carlo method for time-dependent radiative transfer with adaptive material coupling
- Institute for Applied Mathematics and Computational Science, Texas A and M University, College Station, TX 77843-3133 (United States), E-mail: rgm@tamu.edu
- Computational Physics Group (CCS-2), Los Alamos National Laboratory, P.O. Box 1663, MS D409, Los Alamos, NM 87545 (United States), E-mail: tmonster@lanl.gov
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, Vol. 228, Issue 16; Other Information: DOI: 10.1016/j.jcp.2009.04.028; PII: S0021-9991(09)00209-5; Copyright (c) 2009 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
Similar Records
An adaptive alpha for the implicit Monte Carlo equations
Implicit Monte Carlo with a linear discontinuous finite element material solution and piecewise non-constant opacity