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.