Higher order Bernstein-Bézier and Nédélec finite elements for the relaxed micromorphic model

The relaxed micromorphic model is a generalized continuum model that is well-posed in the space $X = [H^1]^3 \times [H(\textrm{curl})]^3$. Consequently, finite element formulations of the model rely on $H^1$-conforming subspaces and N\'ed\'elec elements for discrete solutions of the corresponding variational problem. This work applies the recently introduced polytopal template methodology for the construction of N\'ed\'elec elements. This is done in conjunction with Bernstein-B\'ezier polynomials and dual numbers in order to compute hp-FEM solutions of the model. Bernstein-B\'ezier polynomials allow for optimal complexity in the assembly procedure due to their natural factorization into univariate Bernstein base functions. In this work, this characteristic is further augmented by the use of dual numbers in order to compute their values and their derivatives simultaneously. The application of the polytopal template methodology for the construction of the N\'ed\'elec base functions allows them to directly inherit the optimal complexity of the underlying Bernstein-B\'ezier basis. We introduce the Bernstein-B\'ezier basis along with its factorization to univariate Bernstein base functions, the principle of automatic differentiation via dual numbers and a detailed construction of N\'ed\'elec elements based on Bernstein-B\'ezier polynomials with the polytopal template methodology. This is complemented with a corresponding technique to embed Dirichlet boundary conditions, with emphasis on the consistent coupling condition. The performance of the elements is shown in examples of the relaxed micromorphic model.