Premium
The superconvergent patch recovery and a posteriori error estimates. Part 1: The recovery technique
Author(s) -
Zienkiewicz O. C.,
Zhu J. Z.
Publication year - 1992
Publication title -
international journal for numerical methods in engineering
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.421
H-Index - 168
eISSN - 1097-0207
pISSN - 0029-5981
DOI - 10.1002/nme.1620330702
Subject(s) - superconvergence , quadratic equation , norm (philosophy) , a priori and a posteriori , mathematics , projection (relational algebra) , finite element method , rate of convergence , smoothing , convergence (economics) , simple (philosophy) , mathematical optimization , algorithm , computer science , geometry , statistics , engineering , law , computer network , philosophy , channel (broadcasting) , structural engineering , epistemology , political science , economics , economic growth
This is the first of two papers concerning superconvergent recovery techniques and a posteriori error estimation. In this paper, a general recovery technique is developed for determining the derivatives (stresses) of the finite element solutions at nodes. The implementation of the recovery technique is simple and cost effective. The technique has been tested for a group of widely used linear, quadratic and cubic elements for both one and two dimensional problems. Numerical experiments demonstrate that the recovered nodal values of the derivatives with linear and cubic elements are superconvergent. One order higher accuracy is achieved by the procedure with linear and cubic elements but two order higher accuracy is achieved for the derivatives with quadratic elements. In particular, an O ( h 4 ) convergence of the nodal values of the derivatives for a quadratic triangular element is reported for the first time. The performance of the proposed technique is compared with the widely used smoothing procedure of global L 2 projection and other methods. It is found that the derivatives recovered at interelement nodes, by using L 2 projection, are also superconvergent for linear elements but not for quadratic elements. Numerical experiments on the convergence of the recovered solutions in the energy norm are also presented. Higher rates of convergence are again observed. The results presented in this part of the paper indicate clearly that a new, powerful and economical process is now available which should supersede the currently used post‐processing procedures applied in most codes.