High Continuity Basis’s Impact on Continuous Global L2 (CGL2) Recovery

Abstract

In the recovery-based estimates method, we employ a projection technique to recover a post-processed quantity (usually the stresses or the gradient computed from the FE-approximation). The error is estimated by taking the difference between the recovered quantity and the FE-solution. An easy procedure to implement is the continuous global L2 (CGL2) recovery initially used for a posteriori error estimation by Zienkiewicz and Zhu [1]. Kumar, Kvamsdal and Johannessen [2] developed CGL2 and Superconvergent Patch Recovery (SPR) error estimation methods applicable for adaptive refinement using LR B-splines [3] and observed very good results for both the CGL2 and the SPR recovery technique. However, Cai and Zhang reported in [4] a case of malfunction for the CGL2-recovery applied to second order triangular and tetrahedral Lagrange finite element. Here we will start out by presenting a motivational example that illustrates the benefits of using high regularity splines in the CGL2 based gradient recovery procedure compared to using the classical Lagrange FEM basis functions. We will then show the performance on some benchmark problems comparing the use of splines