Goal-Oriented Regional Angular Adaptive Algorithm for the SN Equations

AbstractAngular discretization errors inherent in the discrete ordinates method are a major problem, especially for localized source problems and problems with strongly absorbing media or large-volume void regions, where angular discretization errors would be totally unacceptable. This paper proposes a regional angular adaptive algorithm together with a goal-oriented error estimate to solve the SN equations. Standard angular adaptive refinement techniques are based on estimated local errors. We compare an interpolated angular flux value against a calculated value to generate local errors. The adaptive quadrature sets can be created by subdividing a spherical quadrilateral into four spherical subquadrilaterals that have positive weights and can be locally refined. Techniques for mapping angular fluxes from one quadrature set to another are developed to transfer angular fluxes on the interfaces of different spatial regions. To provide a better detector response, local errors are weighted by the importance o...