Bernstein-Bezier Bases for Tetrahedral Finite Elements
This work provides a unified, computationally efficient basis for high-order finite element methods on tetrahedral meshes, benefiting computational mechanics and electromagnetics.
The paper introduces new Bernstein-Bezier basis functions for H(curl)-, H(div)-, and L2-conforming finite elements on tetrahedra, enabling efficient matrix assembly without quadrature for constant coefficients. Numerical examples demonstrate their effectiveness.
We present a new set of basis functions for H(curl)-conforming, H(div)-conforming, and L2 -conforming finite elements of arbitrary order on a tetrahedron. The basis functions are expressed in terms of Bernstein polynomials and augment the natural H1 -conforming Bernstein basis. The basis functions respect the differential operators, namely, the gradients of the high-order H1 -conforming Bernstein-Bezier basis functions form part of the H(curl)-conforming basis, and the curl of the high-order, non-gradients H(curl)-conforming basis functions form part of the H(div)-conforming basis, and the divergence of the high-order, non-curl H(div)-conforming basis functions form part of the L2-conforming basis. Procedures are given for the efficient computation of the mass and stiffness matrices with these basis functions without using quadrature rules for (piece-wise) constant coefficients on affine tetrahedra. Numerical results are presented to illustrate the use of the basis to approximate representative problems.