Finite Element Spaces of Double Two-Forms With Polynomial Coefficients

TL;DR

This paper introduces finite element spaces of symmetric tensor products of two-forms with polynomial coefficients, enabling advanced numerical methods in elasticity and general relativity through higher order matrix field representations.

Contribution

It develops new finite element spaces of double two-forms with polynomial coefficients, expanding applications in elasticity and numerical relativity, and provides explicit geometric basis functions.

Findings

Spaces enable higher order finite element methods for elasticity.

Application to Riemann curvature tensor representation.

Explicit basis functions with geometric decomposition.

Abstract

We develop finite element spaces of symmetric tensor products of two-forms with polynomial coefficients. In three dimensions, these give higher order finite element spaces of matrix fields with normal-normal continuity, which have applications to the TDNNS method for elasticity, for example. In general dimension, these spaces can be used to represent the Riemann curvature tensor in numerical relativity. In many ways, our methods parallel Li's work generalizing Regge calculus to higher order, as Regge elements can be thought of as symmetric tensor products of one-forms. However, whereas the constant coefficient Regge space has one shape function per edge, the constant coefficient space of double-forms in our paper has one shape function per triangle and two shape functions per tetrahedron, so we must address the fact that there are shape functions of two different types. Like Li, we obtain an explicit geometrically decomposed basis of shape functions.