Crouzeix-Raviart elements on simplicial meshes in $d$ dimensions
Nis-Erik Bohne, Patrick Ciarlet, Stefan Sauter
TL;DR
This work advances CR non-conforming finite elements to general polynomial order and dimension on simplicial meshes by constructing explicit CR bases using orthogonal polynomials on simplices and face bubbles, and by proving a direct sum decomposition with conforming spaces. It introduces bidual DOFs and a local interpolation operator, enabling stable, local, and consistent approximations that inherit the conforming companion’s approximation properties. In two dimensions with odd $k$, it further develops edge-based DOFs and a practical quasi-interpolation operator; the paper also establishes nonexistence results for split DOFs in higher dimensions ($d\ge3$) and $k\neq1$, delineating the limits of local DOF constructions. Overall, the framework yields robust, high-order, non-conforming CR spaces with explicit bases and biorthogonal DOFs suitable for efficient numerical implementation and analysis.
Abstract
In this paper we introduce Crouzeix-Raviart elements of general polynomial order $k$ and spatial dimension $d\geq2$ for simplicial finite element meshes. We give explicit representations of the non-conforming basis functions and prove that the conforming companion space, i.e., the conforming finite element space of polynomial order $k$ is contained in the Crouzeix-Raviart space. We prove a direct sum decomposition of the Crouzeix-Raviart space into (a subspace of) the conforming companion space and the span of the non-conforming basis functions. Degrees of freedom are introduced which are bidual to the basis functions and give rise to the definition of a local approximation/interpolation operator. In two dimensions or for $k=1$, these freedoms can be split into simplex and $\left( d-1\right) $ dimensional facet integrals in such a way that, in a basis representation of Crouzeix-Raviart functions, all coefficients which belong to basis functions related to lower-dimensional faces in the mesh are determined by these facet integrals. It will also be shown that such a set of degrees of freedom does \textbf{not} exist in higher space dimension and $k>1$.
