A hybrid transport-diffusion method for Monte Carlo radiative-transfer simulations

Discrete Diffusion Monte Carlo (DDMC) is a technique for increasing the efficiency of Monte Carlo particle-transport simulations in diffusive media. If standard Monte Carlo is used in such media, particle histories will consist of many small steps, resulting in a computationally expensive calculation. In DDMC, particles take discrete steps between spatial cells according to a discretized diffusion equation. Each discrete step replaces many small Monte Carlo steps, thus increasing the efficiency of the simulation. In addition, given that DDMC is based on a diffusion equation, it should produce accurate solutions if used judiciously. In practice, DDMC is combined with standard Monte Carlo to form a hybrid transport-diffusion method that can accurately simulate problems with both diffusive and non-diffusive regions. In this paper, we extend previously developed DDMC techniques in several ways that improve the accuracy and utility of DDMC for nonlinear, time-dependent, radiative-transfer calculations. The use of DDMC in these types of problems is advantageous since, due to the underlying linearizations, optically thick regions appear to be diffusive. First, we employ a diffusion equation that is discretized in space but is continuous in time. Not only is this methodology theoretically more accurate than temporally discretized DDMC techniques, but it also has the benefit that a particle's time is always known. Thus, there is no ambiguity regarding what time to assign a particle that leaves an optically thick region (where DDMC is used) and begins transporting by standard Monte Carlo in an optically thin region. Also, we treat the interface between optically thick and optically thin regions with an improved method, based on the asymptotic diffusion-limit boundary condition, that can produce accurate results regardless of the angular distribution of the incident Monte Carlo particles. Finally, we develop a technique for estimating radiation momentum deposition during the DDMC simulation, a quantity that is required to calculate correct fluid motion in coupled radiation-hydrodynamics problems. With a set of numerical examples, we demonstrate that our improved DDMC method is accurate and can provide efficiency gains of several orders of magnitude over standard Monte Carlo.

[1]  Jeffery D. Densmore.  Interface methods for hybrid Monte Carlo-diffusion radiation-transport simulations , 2006 .

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

[3]  Jim E. Morel,et al.  Asymptotically correct angular distributions for Monte Carlo-diffusion interfaces , 2006 .

[4]  Marvin L. Adams,et al.  Finite Element Solutions of Second-Order Forms of the Transport Equation at the Interface Between Diffusive and Non-diffusive Regions , 2001 .

[5]  B. Mercier.  Application of accretive operators theory to the radiative transfer equations , 1987 .

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

[7]  M. Kozmanov,et al.  The maximum principle for a system of equations of energy and non-stationary radiation transfer , 1983 .

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

[9]  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 .

[10]  T. Evans,et al.  1-D Equilibrium Discrete Diffusion Monte Carlo , 2000 .

[11]  Jeffery D. Densmore,et al.  Discrete Diffusion Monte Carlo for grey Implicit Monte Carlo simulations. , 2005 .

[12]  G. C. Pomraning,et al.  Numerical transport and diffusion methods in radiative transfer , 1992 .

[13]  N. A. Gentile.  Implicit Monte Carlo diffusion---an acceleration method for Monte Carlo time-dependent radiative transfer simulations , 2000 .

[14]  Emissivity of discretized diffusion problems , 2006 .

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

[16]  E. Lewis,et al.  Computational Methods of Neutron Transport , 1993 .

[17]  R. Chaplin.  NUCLEAR REACTOR THEORY , 2022 .

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

[19]  J. Giorla,et al.  A random walk method for solving radiative transfer equations , 1987 .

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

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