Nuclear Science and Engineering / Volume 164 / Number 3 / March 2010 / Pages 248-263
A Fourier analysis is conducted for the discrete ordinates, or S
The results for the unaccelerated algorithm show that convergence of IPBJ can degrade, leading in particular to stagnation of the generalized minimum residual method with restart parameter m, GMRES(m), in problems containing optically thin subdomains. The IPBJ algorithm preconditioned with TTSA can be effective, provided the parameter is properly tuned for a given scattering ratio c, but is potentially unstable. Compared to TTSA, MTSA is less sensitive to the choice of , more effective for the same computational effort, measured in terms of the effective scattering ratio c′, and it is unconditionally stable.