Maximum Norm Error Estimates For Neumann Boundary Value Problems On Graded Meshes

This paper deals with a priori pointwise error estimates for the finite element solution of boundary value problems with Neumann boundary conditions in polygonal domains. Due to the corners of the domain, the convergence rate of the numerical solutions can be lower than in the case of smooth domains. As a remedy, the use of local mesh refinement near the corners is considered. In order to prove quasi-optimal a priori error estimates, regularity results in weighted Sobolev spaces are exploited. This is the first work on the Neumann boundary value problem where both the regularity of the data is exactly specified and the sharp convergence order h(2)vertical bar 1n h vertical bar in the case of piecewise linear finite element approximations is obtained. As an extension we show the same rate for the approximate solution of a semilinear boundary value problem. The proof relies in this case on the supercloseness between the Ritz projection to the continuous solution and the finite element solution.