An h-narrow band finite-element method for elliptic equations on implicit surfaces

In this article we define a finite-element method for elliptic partial differential equations (PDEs) on curves or surfaces, which are given implicitly by some level set function. The method is specially designed for complicated surfaces. The key idea is to solve the PDE on a narrow band around the surface. The width of the band is proportional to the grid size. We use finite-element spaces that are unfitted to the narrow band, so that elements are cut off. The implementation nevertheless is easy. We prove error estimates of optimal order for a Poisson equation on a surface and provide numerical tests and examples.