Orthogonality constrained gradient reconstruction for superconvergent linear functionals

The post-processing of the solution of variational problems discretized with Galerkin finite element methods is particularly useful for the computation of quantities of interest. Such quantities are generally expressed as linear functionals of the solution and the error of their approximation is bounded by the error of the solution itself. Several a posteriori recovery procedures have been developed over the years to improve the accuracy of post-processed results. Nonetheless such recovery methods usually deteriorate the convergence properties of linear functionals of the solution and, as a consequence, of the quantities of interest as well. The paper develops an enhanced gradient recovery scheme able to both preserve the good qualities of the recovered gradient and increase the accuracy and the convergence rates of linear functionals of the solution