Abstract. In [8] they introduced a new finite element method for accurate numerical solutions of Poisson equations with corner singularities, which is useful for the problem with known stress intensity factor. They consider the Poisson equations with homogeneous Dirichlet boundary condition, compute the finite element solution using standard FEM and use the extraction formula to compute the stress intensity factor, then they pose a PDE with a regular solution by imposing the nonhomogeneous boundary condition using the computed stress intensity factor, which converges with optimal speed. From the solution they could get accurate solution just by adding the singular part. This approach works for the case when we have the reasonably accurate stress intensity factor. In this paper we consider Poisson equations defined on a domain with a concave corner with Neumann boundary conditions. First we compute the stress intensity factor using the extraction formular, then find the regular part of the solution and the solution.