Distribution Iteration

The discrete ordinates method is widely used to solve the Boltzmann transport equation for neutral particle transport for many engineering applications. Source iteration is used to solve the discrete ordinates system of equations, but can be slow to converge in highly scattering problems. Synthetic acceleration techniques have been developed to address this shortcoming; however, recent research has shown synthetic acceleration to lose effectiveness or diverge for certain problems. LTC Wager introduced an alternative to source iteration and demonstrated it in slab geometry. Here the method is further developed, enhancing efficiency in various ways, and demonstrated in XY-geometry as well as slab geometry. It is shown to be efficient even for those problems for which diffusion-synthetic and transport-synthetic accelerations fail or are ineffective. The method has significant advantages for massively-parallel implementations.