Moment Analysis of Angular APproximation Methods for Time-Dependent Radiation Transport

Abstract

We extend moment analysis, a technique developed for investigating the accuracy of discrete-ordinates spatial discretization schemes, to time-dependent radiation transport and apply it to several angular approximation methods. Specifically, we examine the diffusion approximation, the P 1/3 approximation, and three time-dependent generalizations of the simplified PNapproximation: the SP 2 , SP 3 , and SSP 3 approximations. We show that all of the these methods preserve the correct flux-weighted average of x but not the correct flux-weighted average of (x-xa)2, where x is the spatial variable and xais an arbitrary point. We also demonstrate that, for general cross sections and large elapsed time, the error in the flux-weighted average of (x-xa)2 is smallest in magnitude for the SP 2 and approximations. In addition, we present a simple improvement to the SP 2 approximation that allows this method to produce the correct flux-weighted average of (x-xa)2. We present numerical results that test this analysis. From these results, we find that the angular approximation methods with the most accurate solutions also have the most accurate values for the flux-weighted average of (x-xa)2. In particular, the SP 2 and SP 3 approximations are two of the most accurate methods at large elapsed times, while the improved SP 2 approximation is one of the most accurate methods at all times. We also observe, however, that an accurate value for the flux-weighted average of (x-xa)2 is not always accompanied by an accurate solution. Consequently, we conclude that an accurate flux-weighted average of (x-xa)2 is a necessary rather than sufficient condition for an overall accurate angular approximation method.