Superconvergent patch recovery and a posteriori error estimation technique in adaptive isogeometric analysis

Abstract

Highlights • A posteriori error estimation methodology for adaptive isogeometric analysis using LR B-splines is developed. • Superconvergent patch recovery of gradient field on adaptive meshes is developed. • Algorithm for computation of true superconvergent points on non-uniform adaptive meshes is provided. • Numerical tests verify that the developed error estimator are highly efficient and asymptotically exact. Abstract In this article, we address adaptive methods for isogeometric analysis based on local refinement guided by recovery based a posteriori error estimates. Isogeometric analysis was introduced a decade ago and an impressive progress has been made related to many aspects of numerical methods and advanced applications. However, related to adaptive mesh refinement guided by a posteriori error estimators, rather few attempts are pursued besides the use of classical residual based error estimators. In this article, we explore a feature common for Isogeometric analysis (IGA), namely the use of structured tensorial meshes that facilitates superconvergence behavior of the gradient in the Galerkin discretization. By utilizing the concept of structured mesh refinement using LR B-splines, our aim is to facilitate superconvergence behavior for locally refined meshes as well. Superconvergence behavior matches well with the use of recovery based a posteriori estimator in the Superconvergent Patch Recovery (SPR) procedure. However, to our knowledge so far, the SPR procedure has not been exploited in the IGA community. We start out by addressing the existence of derivative superconvergent points in the computed finite element solution based on B-splines and LR B-splines for an elliptic model problem (1D and 2D Poisson). Then, we present some recovery procedures for improving the derivatives (or gradient) of the isogeometric finite element solution where the SPR procedure will be the main focus. In particular, we show that our SPR procedure for the improvement of derivatives fulfills the desired consistency criteria. At the end, we develop a posteriori error estimator where the improved gradient obtained from the proposed recovery procedures is used. Numerical results are presented to illustrate the efficiency of using SPR procedure for the improvement of derivatives (or gradient) of computed solution in isogeometric analysis. Then the proposed a posteriori error estimator based adaptive refinement methodology is tested to solve smooth and non-smooth elliptic benchmark problems. The focus is put on whether optimal convergence rates are obtained in the computed solution or not, as well as the effectivity index of the proposed error estimators.