Numerical approximation of a singularly perturbed contact problem Original Research Article

As a simplified model for contact problems, we study a mixed Neumann-Robin boundary value problem for the Laplace operator in a smooth domain in R2. The Robin condition contains a small parameter ε inducing boundary layers of corner type at the transition points as proved in [4]. We present an integral equation for the numerical solution of this problem together with estimates of the error. We investigate the improvement obtained if we add to the discrete space some special functions.