In this paper, a unified gas kinetic particle (UGKP) method is developed for radiative transfer in both absorbing and anisotropic scattering media. This numerical method is constructed based on our theoretical work on the model reduction for an anisotropic scattering system. The macroscopic solver of this method directly solves the macroscopic anisotropic diffusion equations, eliminating the need to solve higher-order moment equations. The reconstruction of macroscopic scattering source in the microscopic solver, based on the multiscale equivalent phase function we proposed in this work, has also been simplified as one single scattering process, significantly reducing the computational costs. The proposed method has also the property of asymptotic preserving. In the optically thick regime, the proposed method solves the diffusion limit equations for an anisotropic system. In the optically thin regime, the kinetic processes of photon transport are simulated. The consistency and efficiency of the proposed method have been validated by numerical tests in a wide range of flow regimes. The novel equivalent scattering source reconstruction can be used for various transport processes, and the proposed numerical scheme is widely applicable in high-energy density engineering applications.
Keywords: anisotropic scattering; asymptotic preserving; radiative transfer system; scattering phase function.