NEW NUMERICAL INTEGRATION SCHEMES FOR APPLICATIONS OF GALERKIN BEM TO 2-D PROBLEMS

We consider hypersingular integral formulation of some elasticity and potential boundary value problems on 2-D domains. In particular, we consider all integrals whose evaluation is required when the equations are solved by a Galerkin BEM based on piecewise polynomial approximants of arbitrary local degrees. In order to compute these integrals, we use very e cient formulas recently proposed, which require the user to de ne a mesh, not necessarily uniform, on the boundary and specify the local degrees of the approximant. These rules are quite suitable for the construction of h{p version of the BEM. Implementation of h−, p− and h{p methods are applied to some classical problems and several numerical results are presente