Infinite series solution of the time-dependent radiative transfer equation in anisotropically scattering media

We solve the radiative transfer equation (RTE) in anisotropically scattering media as an infinite series. Each series term represents a distinct number of scattering events, with analytical solutions derived for zero and single scattering. Higher-order corrections are addressed through numerical calculations or approximations.

The RTE solution corresponds to Monte Carlo sampling of photon trajectories with fixed start and end points. Validated against traditional Monte Carlo simulations, featuring random end points, our solution demonstrates enhanced efficiency for both anisotropic and isotropic scattering functions, significantly reducing computational time and resources. The advantage of our method over Monte Carlo simulations varies with the position of interest and the asymmetry of light scattering, but it is typically orders of magnitude faster while achieving the same level of accuracy. The exploitation of hidden symmetries further accelerates our numerical calculations, enhancing the method’s overall efficiency.

In addition, we extend our analysis to the first and second moments of the photon’s flux, elucidating the transition between transport and diffusive regimes.