The singularities that arise in elliptic boundary value problems are treated
locally by a singular function boundary integral method. This method extracts
the leading singular coefficients from a series expansion that describes the
local behavior of the singularity. The method is fitted into the framework of
the widely used boundary element method (BEM), forming a hybrid technique, with
the BEM computing the solution away from the singularity. Results of the hybrid
technique are reported for the Motz problem and compared with the results of
the standalone BEM and Galerkin/finite element method (GFEM). The comparison is
made in terms of the total flux (i.e. the capacitance in the case of
electrostatic problems) on the Dirichlet boundary adjacent to the singularity,
which is essentially the integral of the normal derivative of the solution. The
hybrid method manages to reduce the error in the computed capacitance by a
factor of 10, with respect to the BEM and GFEM