The iterative thermal emission method: A more implicit modification of IMC

For over 40 years, the Implicit Monte Carlo (IMC) method has been used to solve challenging problems in thermal radiative transfer. These problems typically contain regions that are optically thick and diffusive, as a consequence of the high degree of "pseudo-scattering" introduced to model the absorption and reemission of photons from a tightly-coupled, radiating material. IMC has several well-known features that could be improved: a) it can be prohibitively computationally expensive, b) it introduces statistical noise into the material and radiation temperatures, which may be problematic in multiphysics simulations, and c) under certain conditions, solutions can be nonphysical, in that they violate a maximum principle, where IMC-calculated temperatures can be greater than the maximum temperature used to drive the problem.We have developed a variant of IMC called iterative thermal emission IMC, which is designed to have a reduced parameter space in which the maximum principle is violated. ITE IMC is a more implicit version of IMC in that it uses the information obtained from a series of IMC photon histories to improve the estimate for the end of time step material temperature during a time step. A better estimate of the end of time step material temperature allows for a more implicit estimate of other temperature-dependent quantities: opacity, heat capacity, Fleck factor (probability that a photon absorbed during a time step is not reemitted) and the Planckian emission source.We have verified the ITE IMC method against 0-D and 1-D analytic solutions and problems from the literature. These results are compared with traditional IMC. We perform an infinite medium stability analysis of ITE IMC and show that it is slightly more numerically stable than traditional IMC. We find that significantly larger time steps can be used with ITE IMC without violating the maximum principle, especially in problems with non-linear material properties. The ITE IMC method does however yield solutions with larger variance because each sub-step uses a different Fleck factor (even at equilibrium).

[1]  Eugene D. Brooks,et al.  Symbolic implicit Monte Carlo , 1989 .

[2]  Edward W. Larsen,et al.  Analysis of a Monte Carlo method for nonlinear radiative transfer , 1987 .

[3]  H. H. Rosenbrock,et al.  Some general implicit processes for the numerical solution of differential equations , 1963, Comput. J..

[4]  J. A. Fleck,et al.  An implicit Monte Carlo scheme for calculating time and frequency dependent nonlinear radiation transport , 1971 .

[5]  Allan B. Wollaber,et al.  Advanced Monte Carlo methods for thermal radiation transport , 2008 .

[6]  John Giuliani,et al.  Multi-dimensional radiation transport for modeling axisymmetric Z pinches: Ray tracing compared to Monte Carlo solutions for a two-level atom , 2010 .

[7]  Edward W. Larsen,et al.  TOWARDS A FREQUENCY-DEPENDENT DISCRETE MAXIMUM PRINCIPLE FOR THE IMPLICIT MONTE CARLO EQUATIONS , 2010 .

[8]  N. D. Keen,et al.  The Kull IMC package , 1998 .

[9]  Jeffery D. Densmore,et al.  Stability and monotonicity conditions for linear, grey, 0-D implicit Monte Carlo calculations , 2005 .

[10]  Gordon L. Olson,et al.  An analytical benchmark for non-equilibrium radiative transfer in an isotropically scattering medium , 1997 .

[11]  William J. Hogan,et al.  The National Ignition Facility , 2001 .

[12]  Ryan G. McClarren,et al.  A modified implicit Monte Carlo method for time-dependent radiative transfer with adaptive material coupling , 2009, J. Comput. Phys..

[13]  William H. Press,et al.  Numerical recipes in C. The art of scientific computing , 1987 .

[14]  Nick A. Gentile.  Including the effects of temperature-dependent opacities in the implicit Monte Carlo algorithm , 2011, J. Comput. Phys..

[15]  W. Press,et al.  Numerical Recipes: The Art of Scientific Computing , 1987 .

[16]  William H. Press,et al.  Numerical Recipes 3rd Edition: The Art of Scientific Computing , 2007 .

[17]  Leema A. Al-Makhadmeh,et al.  Oxy-fuel coal combustion—A review of the current state-of-the-art , 2011 .

[18]  Jesse R. Cheatham,et al.  Truncation Analysis and Numerical Method Improvements for the Thermal Radiative Transfer Equations. , 2010 .

[19]  G. C. Pomraning.  The non-equilibrium Marshak wave problem , 1979 .

[20]  Scott W. Mosher.  Exact solution of a nonlinear, time-dependent, infinite-medium, grey radiative transfer problem , 2006 .

[21]  Edward W. Larsen,et al.  A linear stability analysis for nonlinear, grey, thermal radiative transfer problems , 2011, J. Comput. Phys..