English

Diffusion Synthetic Acceleration for polytopic discretisations of Boltzmann transport

Numerical Analysis 2026-04-22 v1 Numerical Analysis Computational Physics

Abstract

We present a computational study of diffusion synthetic acceleration (DSA) for the monoenergetic, isotropically scattering SNS_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.60.6.

Keywords

Cite

@article{arxiv.2604.18771,
  title  = {Diffusion Synthetic Acceleration for polytopic discretisations of Boltzmann transport},
  author = {Ansar Calloo and Matthew Evans and François Madiot and Tristan Pryer},
  journal= {arXiv preprint arXiv:2604.18771},
  year   = {2026}
}

Comments

25 pages, 9 figures