Pollution error in the h-version of the finite element method and the local quality of the recovered derivatives

被引:49
作者
Babuska, I [1 ]
Strouboulis, T [1 ]
Gangaraj, SK [1 ]
Upadhyay, CS [1 ]
机构
[1] TEXAS A&M UNIV,DEPT AEROSP ENGN,COLLEGE STN,TX 77843
基金
美国国家科学基金会;
关键词
D O I
10.1016/S0045-7825(96)01013-4
中图分类号
T [工业技术];
学科分类号
08 ;
摘要
In this paper, we address the quality of the solution derivatives which are recovered from finite element solutions by local averaging schemes. As an example, we consider the Zienkiewicz-Zhu superconvergent patch-recovery scheme (the ZZ-SPR scheme), and we study its accuracy in the interior of the mesh for finite element approximations of solutions of Laplace's equation in polygonal domains. We will demonstrate the following: (1) In general, the accuracy of the derivatives recovered by the ZZ-SPR or any other local averaging scheme may not be higher than the accuracy of the derivatives computed directly from the finite element solution. (2) If the mesh is globally adaptive (i.e. it is nearly equilibrated in the energy-norm) then we can, practically always, gain in accuracy by employing the recovered derivatives instead of the derivatives computed directly from the finite element solution. (3) It is possible to guarantee that the recovered solution-derivatives have higher accuracy than the derivatives computed directly from the finite element solution, in any patch of elements of interest, by employing a mesh which is adaptive only with respect to the patch of interest (i.e. it is nearly equilibrated in a weighted energy-norm). (4) In practice, we are often interested in obtaining highly accurate derivatives (or heat-fluxes, stresses, etc.) only in a few critical regions which are identified by a preliminary analysis. A grid which is adaptive only with respect to the critical regions of interest may be much more economical for this purpose because it may achieve the desired accuracy by employing substantially fewer degrees of freedom than a globally adaptive grid which achieves comparable accuracy in the critical regions.
引用
收藏
页码:1 / 37
页数:37
相关论文
共 26 条
[1]   A UNIFIED APPROACH TO A POSTERIORI ERROR ESTIMATION USING ELEMENT RESIDUAL METHODS [J].
AINSWORTH, M ;
ODEN, JT .
NUMERISCHE MATHEMATIK, 1993, 65 (01) :23-50
[2]  
Aminpour M. A., 1992, ADAPTIVE MULTILEVEL, V157, P119
[3]  
AMINPOUR MA, 1992, P 33 AIAA STRUCT S 1, V1, P109
[4]   RELIABLE STRESS AND FRACTURE-MECHANICS ANALYSIS OF COMPLEX COMPONENTS USING A H-P VERSION OF FEM [J].
ANDERSSON, B ;
FALK, U ;
BABUSKA, I ;
VONPETERSDORFF, T .
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, 1995, 38 (13) :2135-2163
[5]  
[Anonymous], 1996, Numerical Methods for Partial Differential Equations
[6]   A-POSTERIORI ESTIMATION AND ADAPTIVE-CONTROL OF THE POLLUTION ERROR IN THE H-VERSION OF THE FINITE-ELEMENT METHOD [J].
BABUSKA, I ;
STROUBOULIS, T ;
UPADHYAY, CS ;
GANGARAJ, SK .
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, 1995, 38 (24) :4207-4235
[7]   ETA-PERCENT-SUPERCONVERGENCE OF FINITE-ELEMENT APPROXIMATIONS IN THE INTERIOR OF GENERAL MESHES OF TRIANGLES [J].
BABUSKA, I ;
STROUBOULIS, T ;
UPADHYAY, CS .
COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING, 1995, 122 (3-4) :273-305
[8]  
BABUSKA I, 1984, NUMER MATH, V44, P75, DOI 10.1007/BF01389757
[9]   A MODEL STUDY OF ELEMENT RESIDUAL ESTIMATORS FOR LINEAR ELLIPTIC PROBLEMS - THE QUALITY OF THE ESTIMATORS IN THE INTERIOR OF MESHES OF TRIANGLES AND QUADRILATERALS [J].
BABUSKA, I ;
STROUBOULIS, T ;
UPADHYAY, CS ;
GANGARAJ, SK .
COMPUTERS & STRUCTURES, 1995, 57 (06) :1009-1028
[10]  
Babuska I., 1980, Computational Methods in Nonlinear Mechanics. Proceedings of the TICOM Second International Conference, P67