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