An hpr-mesh refinement algorithm for BEM

This paper describes a mesh refinement scheme for boundary element method in which the number of elements, the size of elements, element end locations and the element polynomial order are determined to meet the user specified accuracy. The use of grading function in conjunction with the L1 norm makes the mesh refinement scheme applicable to a variety of boundary element formulation and applications. The algorithm is stable for smooth, discontinuous, as well as singular density functions. Numerical results for mathematical test functions as well as for elastostatic problems demonstrate the viability and versatility of the algorithm for BEM applications.