Simulating radiative transfer with filtered spherical harmonics

This Letter presents a novel application of filters to the spherical harmonics (PN)(PN) expansion for radiative transfer problems in the high-energy-density regime. The filter, which is based on non-oscillatory spherical splines, preserves both the equilibrium diffusion limit and formal convergence properties of the unfiltered expansion. While the method requires further mathematical justification and computational studies, preliminary results demonstrate that solutions to the filtered PNPN equations are (1) more robust and less oscillatory than standard PNPN solutions and (2) more accurate than discrete ordinates solutions of comparable order. The filtered P7P7 solution demonstrates comparable accuracy to an implicit Monte Carlo solution for a benchmark hohlraum problem. Given the benefits of this method we believe it will enable more routine use of high-fidelity radiation-hydrodynamics calculations in the simulation of physical systems.