Abstract: In this study, a highly-accurate, conforming finite element method is developed and justified for the solution of the Dirichlet problem of the generalized Helmholtz equation on domains with re-entrant corners. The k − t h order Lagrange elements are used for the discretization of the variational form of the problem on exponentially compressed polar meshes employed in the neighbourhood of the corners whose interior angle is α π , α ≠ 1 ∕ 2 , and on the triangular and curved mesh formed in the remainder of the polygon. The exponentially compressed polar meshes are constructed such that they are transformed to square meshes using the Log-Polar transformation, simplifying the realization of the method significantly. For the error bound between the exact and the approximate solution obtained by the proposed method, an accuracy of O ( h k ) , h mesh size and k ≥ 1 an integer, is obtained in the H 1 -norm. Numerical experiments are conducted to support the theoretical analysis made. The proposed method can be applied for dealing with the corner singularities of general nonlinear parabolic partial differential equations with semi-implicit time discretization.