Goal-oriented spatial adaptivity for the S N equations on unstructured triangular meshes

• Published in: Annals of nuclear energy 2011-06, Vol.38 (6), p.1373-1381
• Main Author: Lathouwers, D.
• Format: Article
• Language: English
• Subjects: A posteriori error analysis; Adaptive mesh refinement; Adjoints; Algorithms; Economics; Errors; Finite element; Forward problem; Heterogeneity; Mathematical analysis; Mathematical models; Neutron transport
• ISSN: 0306-4549; 1873-2100
• DOI: 10.1016/j.anucene.2011.01.038

► We investigate goal-oriented adaptivity for the S N equations. ► The method is based on the use an adjoint weighted error measure. ► The algorithm is applied to unstructured triangular meshes. ► The resulting algorithm is efficient in terms of the elements required. In this paper we investigate a goal-oriented adaptive algorithm for the S N equations on unstructured meshes. The method is based on a dual-weighted residual approach where an appropriate adjoint problem is formulated and solved in order to obtain the importance of residual errors in the forward problem on the specific goal of interest. The forward residuals and the adjoint function are combined to obtain both economical finite element meshes tailored to the solution of the target functional as well as providing reliable error estimates. The suitability of the approach is demonstrated by performing a series of numerical experiments in two shielding geometries exhibiting strong heterogeneity and flux gradients. Tests are performed with varying target functionals. The effectiveness of the method is shown by comparing the number of elements required for a specified accuracy with that necessary using uniform refinement to give an indication what savings can be achieved by adaptive refinement.