Asymptotic analysis of the spatial discretization of radiation absorption and re-emission in Implicit Monte Carlo

We perform an asymptotic analysis of the spatial discretization of radiation absorption and re-emission in Implicit Monte Carlo (IMC), a Monte Carlo technique for simulating nonlinear radiative transfer. Specifically, we examine the approximation of absorption and re-emission by a spatially continuous artificial-scattering process and either a piecewise-constant or piecewise-linear emission source within each spatial cell. We consider three asymptotic scalings representing (i) a time step that resolves the mean-free time, (ii) a Courant limit on the time-step size, and (iii) a fixed time step that does not depend on any asymptotic scaling. For the piecewise-constant approximation, we show that only the third scaling results in a valid discretization of the proper diffusion equation, which implies that IMC may generate inaccurate solutions with optically large spatial cells if time steps are refined. However, we also demonstrate that, for a certain class of problems, the piecewise-linear approximation yields an appropriate discretized diffusion equation under all three scalings. We therefore expect IMC to produce accurate solutions for a wider range of time-step sizes when the piecewise-linear instead of piecewise-constant discretization is employed. We demonstrate the validity of our analysis with a set of numerical examples.

[1]  G. Habetler,et al.  Uniform asymptotic expansions in transport theory with small mean free paths, and the diffusion approximation , 1975 .

[2]  E. Larsen,et al.  Asymptotic solutions of numerical transport problems in optically thick, diffusive regimes II , 1989 .

[3]  Edward W. Larsen,et al.  Asymptotic analysis of radiative transfer problems , 1983 .

[4]  Edward W. Larsen.  Diffusion theory as an asymptotic limit of transport theory for nearly critical systems with small mean free paths , 1980 .

[5]  Marvin L. Adams,et al.  Asymptotic Analysis of a Computational Method for Time- and Frequency-Dependent Radiative Transfer , 1998 .

[6]  G. C. Pomraning.  The Equations of Radiation Hydrodynamics , 2005 .

[7]  G. Samba,et al.  Asymptotic diffusion limit of the symbolic Monte-Carlo method for the transport equation , 2004 .

[8]  J. A. Fleck,et al.  A random walk procedure for improving the computational efficiency of the implicit Monte Carlo method for nonlinear radiation transport , 1984 .

[9]  T. NKaoua,et al.  Solution of the Nonlinear Radiative Transfer Equations by a Fully Implicit Matrix Monte Carlo Method Coupled with the Rosseland Diffusion Equation via Domain Decomposition , 1991, SIAM J. Sci. Comput..

[10]  J. Keller,et al.  Asymptotic solution of neutron transport problems for small mean free paths , 1974 .

[11]  Edward W. Larsen,et al.  Asymptotic equilibrium diffusion analysis of time-dependent Monte Carlo methods for grey radiative transfer , 2004 .

[12]  R. LeVeque.  Finite Volume Methods for Hyperbolic Problems: Characteristics and Riemann Problems for Linear Hyperbolic Equations , 2002 .

[13]  J. Densmore.  ASYMPTOTIC ANALYSIS OF SPATIAL DISCRETIZATIONS IN IMPLICIT MONTE CARLO , 2009 .

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

[15]  Jim E. Morel,et al.  A Linear-Discontinuous Spatial Differencing Scheme forSnRadiative Transfer Calculations , 1996 .

[16]  Marvin L. Adams,et al.  Discontinuous Finite Element Transport Solutions in Thick Diffusive Problems , 2001 .

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

[18]  G. C. Pomraning,et al.  Linear Transport Theory , 1967 .

[19]  Eugene D. Brooks,et al.  Comparison of implicit and symbolic implicit Monte Carlo line transport with frequency weight vector extension , 2003 .

[20]  D. Mihalas,et al.  Foundations of Radiation Hydrodynamics , 1985 .

[21]  Edward W. Larsen,et al.  The Asymptotic Diffusion Limit of Discretized Transport Problems , 1992 .