The Quasi-Diffusion method for solving transport problems in planar and spherical geometries

Abstract

The Quasi-Diffusion (QD) method is a nonlinear algorithm for solving linear transport problems. Because the QD method utilizes both a transport sweep and a diffusion calculation within each iteration, it is operationally more complex than the Source Iteration (SI) method, which utilizes only a transport sweep within each iteration. However, the QD method often converges rapidly and with high accuracy, especially for optically thick regions with scattering ratios close to unity; these are the regions for which acceleration is most needed. A difficulty with the QD method is that because it is nonlinear, every scalar flux iterate must be positive at each point in the system. Also, the formulation of diffusion boundary conditions to optimize accuracy and speed of convergence is not obvious. In this paper, both of these issues is considered. Specifically, a new formulation of the QD method in spherical geometry is proposed to guarantee positivity of the analytic solution, and new diffusion boundary conditions are proposed in planar and spherical geometry that lead to more accurate and efficient solution algorithms. Also, ways to accurately and positively discretize the transport and diffusion equations, and extensive numerical results are presented. 14 refs., 1 fig., 9 tabs.