On discontinuous Galerkin and discrete ordinates approximations for neutron transport equation and the critical eigenvalue

摘要

The objective of this paper is to give a mathematical framework fora fully discrete numerical approach for the study of the neutron transport equation in a cylindrical domain (container model). More specifically, we consider the discontinuous Galerkin (DG) finite element method for spatial approximation of the mono-energetic, critical neutron transport equation in an infinite cylindrical domain wiRM3ith a pnolygon~al convex cros-section Q. The velocity discretization relies on a special quadrature rule developed to give optimal estimates in discrete ordinate parameters compatible with the quasi-uniform spatial mesh. We use interpolation spaces and derive optimal error estimates, up to maximal available regularity, for the fully discrete scalar flux. Finally we employ a duality argument and prove superconvergence estimates for the critical eigenvalue.