CMFD acceleration of spatial domain-decomposed neutron transport problems

A significant limitation to parallelizing the solution of neutron transport problems is the need for sweeps across the entirety of the problem domain. Angular domain decomposition is common practice, as the equations for each direction are independent aside from their shared scattering/fission source. Accordingly, spatial domain decomposition does not naturally arise in the transport equations and is therefore less frequent in practice. In this paper, we show that a neutron transport domain can be straightforwardly divided into independent, parallelizable sweep regions, globally linked with the standard CMFD method, with an additional update equation. We verify, theoretically (via Fourier analysis) and computationally, that the convergence properties of this method are stable and nominally as rapid as standard CMFD. (authors)