Computing flux intensity factors by a boundary method for elliptic equations with singularities

A simple method for computing the flux intensity factors associated with the asymptotic solution of elliptic equations having a large convergence radius in the vicinity of singular points is presented. The Poisson and Laplace equations over domains containing boundary singularities due to abrupt change of the boundary geometry or boundary conditions are considered. The method is based on approximating the solution by the leading terms of the local symptotic expansion, weakly enforcing boundary conditions by minimization of a norm on the domain boundary in a least-squares sense. The method is applied to the Motz problem, resulting in extremely accurate estimates for the flux intensity factors. It is shown that the method converges exponentially with the number of singular functions and requires a low computational cost. Numerical results to a number of problems concerned with the Poisson equation over an L-shaped domain, and over domains containing multiple singular points, demonstrate accurate estimates for the flux intensity factors.