A method for solving the Neumann problem for the Poisson equation on the exterior of a poligon

This paper is concerned with a numerical method for solving the problem Δu=f in Ωc (=intR2\Ω), (du/dn)|Γ=g, where ΩR2 is a polygon and Γ is the boundary of Ω. The method is based on coupling finite and boundary element techniques. To compensate for the loss of smoothness of the solution u near the corners of the polygon Ω we refine the triangulation without changing the number of triangles. We apply the affine triangular Lagrangean element of degree kN and the Lagrangean boundary element of degree k−1 to obtain the optimal order of convergence via the Galerkin projection.