Among a few known techniques the isoparametric version of the finite element method for meshes consisting of curved triangles or tetrahedra is the one most widely employed to solve PDEs with essential conditions prescribed on curved boundaries. It allows to recover optimal approximation properties that hold for elements of order greater than one in the energy norm for polytopic domains. However, besides a geometric complexity, this technique requires the manipulation of rational functions and the use of numerical integration. We consider a simple alternative to deal with Dirichlet boundary conditions that bypasses these drawbacks, without eroding qualitative approximation properties. In the present work we first recall the main principle this technique is based upon, by taking as a model the solution of the Poisson equation with quadratic Lagrange finite elements. Then we show that it extends very naturally to viscous incompressible flow problems. Although the technique applies to any higher order velocity–pressure pairing, as an illustration a thorough study thereof is conducted in the framework of the Stokes system solved by the classical Taylor–Hood method.
