DIFFUSION SYNTHETIC ACCELERATION PRECONDITIONING FOR DISCONTINUOUS GALERKIN DISCRETIZATIONS OF SN TRANSPORT ON HIGH-ORDER CURVED MESHES.

This paper derives and analyzes new diffusion synthetic acceleration (DSA) pre-conditioners for the S_N transport equation when discretized with a high-order (HO) discontinuous Galerkin (DG) discretization. DSA preconditioners address the need to accelerate the S_N transport equation when the mean free path ε of particles is small and the condition number of the S_N transport equation scales like O(ε^-2). By expanding the S_N transport operator in " and employing a rigorous singular matrix perturbation analysis, we derive a DSA matrix that reduces to the symmetric interior penalty (SIP) DG discretization of the standard continuum diffusion equation when the mesh is first-order and the total opacity is constant. We prove that preconditioning the HO DG S_N transport equation with the SIP DSA matrix results in an O(ε) perturbation of the identity, and fixed-point iteration therefore converges rapidly for optically thick problems. However, the SIP DSA matrix is conditioned like O(ε^-1), making it difficult to invert for small ε. We further derive a new two-part, additive DSA preconditioner based on a continuous Galerkin discretization of diffusion-reaction, which has a condition number independent of ε, and prove that this DSA variant has the same theoretical efficiency as the SIP DSA preconditioner in the optically thick limit. The analysis is extended to the case of HO (curved) meshes, where so-called mesh cycles can result from elements both being upwind of each other (for a given discrete photon direction). In particular, we prove that performing two additional transport sweeps, with fixed scalar flux, in between DSA steps yields the same theoretical conditioning of fixed-point iterations as in the cycle-free case. Theoretical results are validated by numerical experiments on a HO, highly curved two- and three-dimensional meshes that are generated from an arbitrary Lagrangian-Eulerian hydrodynamics code, where the additional inner sweeps between DSA steps offer up to a 4x reduction in total number of sweeps required for convergence. [ABSTRACT FROM AUTHOR]