A High-Order Nonlinear Diffusion Acceleration for the SN Equations Discretized with the Discontinuous FEM I: Theory and Numerical Results
Recently, nonlinear diffusion acceleration (NDA) methods have received widespread attention in the radiation transport community because of their inherent efficiency (Willert 2014, Wang 2014, Park 2012). However, NDA is in general only conditionally stable \cite{Larsen2014}, while diffusion synthetic acceleration methods (DSA) are unconditionally stable if suitable discretization methods for the transport and diffusion systems are chosen \cite{Larsen1982}. The presented work introduces a new generalized NDA method for solving the transport equation discretized with $S_N$ (discrete ordinates) and discontinuous Finite Elements. In contrast to common coarse-mesh finite difference methods (CMFD), where the diffusion system is discretized using constant shape functions, the presented NDA method discretizes the diffusion system using the same (higher order) finite element basis functions as those for the transport system. This paper focuses on the solution of the steady-state transport equation with the low-order diffusion equation employing the same mesh and the same spatial shape functions, but the NDA methods currently under development for the Rattle$S_N$ake radiation transport solver\cite{Wang2013} based on the MOOSE\cite{Gaston2009} framework is more general: it allows for the diffusion system to be discretized using an arbitrary group structure and arbitrary mesh nested in their $S_N$ system counterparts, and arbitrary combinations of FEM basis functions.