High-order spherical harmonics–nodal collocation scheme for the numerical solution of the time-dependent radiative transfer equation

Abstract

The validity of the spherical harmonics-nodal collocation approximation to the stationary Boltzmann equation has been established in multi-dimensional neutron transport problems. This is a high-order approximation method that allows a coarse spatial discretization without losing accuracy. We extend the method to solve the time-dependent radiative transfer equation for absorbing media with anisotropic scattering, also incorporating to the model reflecting boundary conditions, due to the refractive index m ismatching. The formalism is then applied to numerical test cases in one and two spatial dimensions that, using typical values in optical tomography for the physical parameters, check the accuracy and convergence of the method.