A posteriori goal-oriented bounds for the Poisson problem using potential and equilibrated flux reconstructions: Application to the hybridizable discontinuous Galerkin method
We present a general framework to compute upper and lower bounds for linear-functional outputs of the exact solutions of the Poisson equation based on reconstructions of the field variable and flux for both the primal and adjoint problems. The method is devised from a generalization of the complemen...
| Authors: | , , , |
|---|---|
| Format: | article |
| Publication Date: | 2021 |
| Country: | España |
| Institution: | Universitat Politècnica de Catalunya (UPC) |
| Repository: | UPCommons. Portal del coneixement obert de la UPC |
| Language: | English |
| OAI Identifier: | oai:upcommons.upc.edu:2117/359138 |
| Online Access: | https://hdl.handle.net/2117/359138 https://dx.doi.org/10.1016/j.cma.2021.114088 |
| Access Level: | Open access |
| Keyword: | Applied mathematics Exact/guaranteed/strict bounds for quantities of interest Output bounds Goal-oriented error estimation Adaptivity Potential and equilibrated flux reconstructions Hybridizable discontinuous Galerkin method (HDG) Matemàtica aplicada Àrees temàtiques de la UPC::Matemàtiques i estadística |
| Summary: | We present a general framework to compute upper and lower bounds for linear-functional outputs of the exact solutions of the Poisson equation based on reconstructions of the field variable and flux for both the primal and adjoint problems. The method is devised from a generalization of the complementary energy principle and the duality theory. Using duality theory, the computation of bounds is reduced to finding independent potential and equilibrated flux reconstructions. A generalization of this result is also introduced allowing to derive alternative guaranteed bounds from nearly-arbitrary flux reconstructions (only zero-order equilibration is required). This approach is applicable to any numerical method used to compute the solution. In this work, the proposed approach is applied to derive bounds for the hybridizable discontinuous Galerkin (HDG) method. An attractive feature of the proposed approach is that superconvergence on the bound gap is achieved, yielding accurate bounds even for very coarse meshes. Numerical experiments are presented to illustrate the performance and convergence of the bounds for the HDG method in both uniform and adaptive mesh refinements. |
|---|