New robust nonconforming finite elements of higher order

We show that existing quadrilateral nonconforming finite elements of higher order exhibit a reduction in the order of approximation if the sequence of meshes is still shape-regular but consists no longer of asymptotically affine equivalent mesh cells. We study second order nonconforming finite elements as members of a new family of higher order approaches which prevent this order reduction. We present a new approach based on the enrichment of the original polynomial space on the reference element by means of nonconforming cell bubble functions which can be removed at the end by static condensation. Optimal estimates of the approximation and consistency error are shown in the case of a Poisson problem which imply an optimal order of the discretization error. Moreover, we discuss the known nonparametric approach to prevent the order reduction in the case of higher order elements, where the basis functions are defined as polynomials on the original mesh cell. Regarding the efficient treatment of the resulting linear discrete systems, we analyze numerically the convergence of the corresponding geometrical multigrid solvers which are based on the canonical full order grid transfer operators. Based on several benchmark configurations, for scalar Poisson problems as well as for the incompressible Navier–Stokes equations (representing the desired application field of these nonconforming finite elements), we demonstrate the high numerical accuracy, flexibility and efficiency of the discussed new approaches which have been successfully implemented in the FeatFlow software (www.featflow.de). The presented results show that the proposed FEM-multigrid combinations (together with discontinuous pressure approximations) appear to be very advantageous candidates for efficient simulation tools, particularly for incompressible flow problems.