Energy-Dependent, Self-Adaptive Mesh h(p)-Refinement of an Interior-Penalty Scheme for a Discontinuous Galerkin Isogeometric Analysis Spatial Discretization of the Multi-Group Neutron Diffusion Equation with Dual-Weighted Residual Error Measures

S. G. WilsonM. D. EatonJ. Kópházia Nuclear Engineering Group,Department of Mechanical Engineering,City and Guilds Building,Imperial College London,London,UKb Rolls-Royce plc,Derby,UKc Institute of Nuclear Techniques,Budapest University of Technology and Economics,Budapest,Hungary
DOI: https://doi.org/10.1080/23324309.2024.2334277
IF: 0.348
2024-04-15
Journal of Computational and Theoretical Transport
Abstract:Energy-dependent self-adaptive mesh refinement algorithms are developed for a symmetric interior-penalty scheme for a discontinuous Galerkin spatial discretization of the multi-group neutron diffusion equation using NURBS-based isogeometric analysis (IGA). The spatially self-adaptive algorithms employ both mesh ( h ) and polynomial degree ( p ) refinement. The discretized system becomes increasingly ill-conditioned for increasingly large penalty parameters; and there is no gain in accuracy for over penalization. Therefore, optimized penalty parameters are rigorously calculated, for general element types, from a coercivity analysis of the bilinear form. Local mesh refinement allows for a better allocation of computational resources; and thus, more accuracy per degree of freedom. Two a posteriori interpolation-based error measures are proposed. The first heuristically minimizes local contributions to the discretization error, which becomes competitive for global quantities of interest (QoIs). However, for localized QoIs, over energy-dependent meshes, certain multi-group components may become under-resolved. The second employs duality arguments to minimize important error contributions, which consistently and reliably reduces the error in the QoI.
mathematics, applied,physics, mathematical
What problem does this paper attempt to address?