Diffusion Synthetic Acceleration for polytopic discretisations of Boltzmann transport
For computational transport, this provides a robust acceleration scheme for polytopic DG discretizations, though the improvement is incremental over existing DSA methods.
The paper studies diffusion synthetic acceleration (DSA) for Boltzmann transport discretized with polytopic DG, comparing SIP and MIP variants. MIP-based DSA remains robust across parameters, with convergence factors below 0.6 in challenging regimes.
We present a computational study of diffusion synthetic acceleration (DSA) for the monoenergetic, isotropically scattering $S_N$ transport equations, discretised in space by a polytopic discontinuous Galerkin method. Using a discrete ordinates angular discretisation, we construct the DSA correction with an interior-penalty diffusion operator and compare a classical symmetric interior penalty (SIP) formulation with a modified interior penalty (MIP) variant, together with homogeneous Dirichlet and Marshak (Robin) diffusion boundary conditions imposed weakly in the DG framework. We quantify the observed convergence behaviour of the resulting source iteration across variations in optical thickness, scattering ratio, angular quadrature, mesh refinement, polynomial degree and mesh anisotropy on families of bounded Voronoi meshes. The results show that MIP-based DSA remains robust across the parameter ranges tested, whereas SIP-based DSA can lose robustness in the intermediate regime. In challenging optically thick, highly scattering settings, the observed convergence factors for the MIP-based schemes are typically below $0.6$.