On a general implementation of hh- and pp-adaptive curl-conforming finite elements

Edge (or N\'ed\'elec) finite elements are theoretically sound and widely used by the computational electromagnetics community. However, its implementation, specially for high order methods, is not trivial, since it involves many technicalities that are not properly described in the literature. To fill this gap, we provide a comprehensive description of a general implementation of edge elements of first kind within the scientific software project FEMPAR. We cover into detail how to implement arbitrary order (i.e., $p$-adaptive) elements on hexahedral and tetrahedral meshes. First, we set the three classical ingredients of the finite element definition by Ciarlet, both in the reference and the physical space: cell topologies, polynomial spaces and moments. With these ingredients, shape functions are automatically implemented by defining a judiciously chosen polynomial pre-basis that spans the local finite element space combined with a change of basis to automatically obtain a canonical basis with respect to the moments at hand. Next, we discuss global finite element spaces putting emphasis on the construction of global shape functions through oriented meshes, appropriate geometrical mappings, and equivalence classes of moments, in order to preserve the inter-element continuity of tangential components of the magnetic field. Finally, we extend the proposed methodology to generate global curl-conforming spaces on non-conforming hierarchically refined (i.e., $h$-adaptive) meshes with arbitrary order finite elements. Numerical results include experimental convergence rates to test the proposed implementation

Curl-conforming hierarchical vector bases for triangles and tetrahedra

A new family of hierarchical vector bases is proposed for triangles and tetrahedra. These functions span the curl-conforming reduced-gradient spaces of Nédélec. The bases are constructed from orthogonal scalar polynomials to enhance their linear independence, which is a simpler process than an orthogonalization applied to the final vector functions. Specific functions are tabulated to order 6.5. Preliminary results confirm that the new bases produce reasonably well-conditioned matrice

High-order finite elements on pyramids: approximation spaces, unisolvency and exactness

We present a family of high-order finite element approximation spaces on a pyramid, and associated unisolvent degrees of freedom. These spaces consist of rational basis functions. We establish conforming, exactness and polynomial approximation properties.Comment: 37 pages, 3 figures. This work was originally in one paper, then split into two; it has now been recombined into one paper, with substantial changes from both of its previous form