Polytopal templates for semi-continuous vectorial finite elements of arbitrary order on triangulations and tetrahedralizations

The Hilbert spaces $H\left(\mathrm{curl}\right)$ and $H\left(\mathrm{div}\right)$ are employed in various variational problems formulated in the context of the de Rham complex in order to guarantee well-posedness. Seeing as the well-posedness follows automatically from the continuous setting to the discrete setting in the presence of commuting interpolants as per Fortin’s criterion, the construction of conforming subspaces becomes a crucial step in the formulation of stable numerical schemes. This work aims to introduce a novel, simple method of directly constructing semi-continuous vectorial base functions on the reference element via template vectors associated with the geometric polytopes of the element and an underlying $H^1$-conforming polynomial subspace. The base functions are then mapped from the reference element to the element in the physical domain via consistent Piola transformations. The method is defined in such a way, that the underlying $H^1$-conforming subspace can be chosen independently, thus allowing for constructions of arbitrary polynomial order. We prove a linearly independent construction of Nédélec elements of the first and second type, Brezzi–Douglas–Marini elements, and Raviart–Thomas elements on triangulations and tetrahedralizations. The application of the method is demonstrated with two examples in the relaxed micromorphic model.