An Integral-Equation-Based Variance Reduction Method for Accelerated Monte Carlo Simulations

Abstract

In this work, we introduce a novel variance reduction approach utilising the integral formulation of the radiative transfer equation to calculate the radiance in a planar symmetric slab geometry. Due to its integral nature, our method offers a fundamental advantage over well-established variance reduction methods such as the local estimate technique. As opposed to the local estimate procedure, photons add to the overall radiance not only at specific points of interaction but also throughout each consecutive path element; hence, our variance reduction approach can be thought of as an integral local estimate method. This facilitates a substantial enhancement in statistical efficiency, especially in scenarios where only a small number of scattering events or a high attenuation along the detection paths is to be anticipated. To evaluate the overall performance of the integral approach, we incorporated it into a self-developed GPU-accelerated Monte Carlo software, together with a conventional local estimate implementation adapted to slab geometry for a comprehensive comparison.