TY - JOUR

T1 - Analysis of hypersingular residual error estimates in boundary element methods for potential problems

AU - Menon, Govind

AU - Paulino, Glaucio H.

AU - Mukherjee, Subrata

N1 - Funding Information:
Prof. Paulino acknowledges the support from National Science Foundation under Grant No. CMS-9713008. The authors thank Prof. Lars Wahlbin, Department of Mathematics, and Prof. Steven Vavasis, Department of Computer Science, both from Cornell University, for elucidating discussions on the content of this paper. The idea behind the superconvergence argument of Section 4.1.2 is due to Prof. Wahlbin.

PY - 1999/5/27

Y1 - 1999/5/27

N2 - A novel iteration scheme, using boundary integral equations, is developed for error estimation in the boundary element method. The iteration scheme consists of using the boundary integral equation for solving the boundary value problem and iterating this solution with the hypersingular boundary integral equation to obtain a new solution. The hypersingular residual r is consistently defined as the difference in the derivative quantities on the boundary, i.e. r=∂φ(1)/∂n-∂μ(2)/∂n where φ is the potential and (∂μ/∂n)(i),i = 1,2, is the flux obtained by solution (i). Here, i = 1 refers to the boundary integral equation, and i = 2 refers to the hypersingular boundary integral equation. The hypersingular residual is interpreted in the sense of the iteration scheme defined above and it is shown to provide an error estimate for the boundary value problem. Error-hypersingular residual relations are developed for Dirichlet and Neumann problems, which are shown to be limiting cases of the more general relation for the mixed boundary value problem. These relations lead to global bounds on the error. Four numerical examples, involving Galerkin boundary elements, are given, and one of them involves a physical singularity on the boundary and preliminary adaptive calculations. These examples illustrate important features of the hypersingular residual error estimate proposed in this paper.

AB - A novel iteration scheme, using boundary integral equations, is developed for error estimation in the boundary element method. The iteration scheme consists of using the boundary integral equation for solving the boundary value problem and iterating this solution with the hypersingular boundary integral equation to obtain a new solution. The hypersingular residual r is consistently defined as the difference in the derivative quantities on the boundary, i.e. r=∂φ(1)/∂n-∂μ(2)/∂n where φ is the potential and (∂μ/∂n)(i),i = 1,2, is the flux obtained by solution (i). Here, i = 1 refers to the boundary integral equation, and i = 2 refers to the hypersingular boundary integral equation. The hypersingular residual is interpreted in the sense of the iteration scheme defined above and it is shown to provide an error estimate for the boundary value problem. Error-hypersingular residual relations are developed for Dirichlet and Neumann problems, which are shown to be limiting cases of the more general relation for the mixed boundary value problem. These relations lead to global bounds on the error. Four numerical examples, involving Galerkin boundary elements, are given, and one of them involves a physical singularity on the boundary and preliminary adaptive calculations. These examples illustrate important features of the hypersingular residual error estimate proposed in this paper.

UR - http://www.scopus.com/inward/record.url?scp=0032688898&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=0032688898&partnerID=8YFLogxK

U2 - 10.1016/S0045-7825(98)00297-7

DO - 10.1016/S0045-7825(98)00297-7

M3 - Article

AN - SCOPUS:0032688898

VL - 173

SP - 449

EP - 473

JO - Computer Methods in Applied Mechanics and Engineering

JF - Computer Methods in Applied Mechanics and Engineering

SN - 0374-2830

IS - 3-4

ER -