Reliable anisotropic-adaptive discontinuous Galerkin method for simplified PN approximations of radiative transfer

In this paper we present the analysis for the error estimator for radiative transfer problems presented in Giani and Seaid (2016) where we showed the capabilities of the error estimator to accurately drive the adaptivity to resolve steep boundary layers, which are among the difficulties that most numerical methods fail to resolve accurately. We prove reliability for the error estimator in terms of a global upper bound of the error measured in the natural norm. We present a series of numerical experiments to test the efficiency of this approach within a fully automated hp-adaptive refinement algorithm.