On a general implementation of $h$- and $p$-adaptive curl-conforming finite elements

Abstract

Edge (or Nédélec) 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.

References in corpus (8)

• The aggregated unfitted finite element method for elliptic problems
• Efficient Assembly of H(div) and H(curl) Conforming Finite Elements
• Orientation Embedded High Order Shape Functions for the Exact Sequence Elements of All Shapes
• Fast MATLAB assembly of FEM matrices in 2D and 3D: Edge elements
• H(div) and H(curl)-conforming VEM
• On orienting edges of unstructured two- and three-dimensional meshes
• An example of explicit implementation strategy and preconditioning for the high order edge finite elements applied to the time-harmonic Maxwell's equations
• A generic finite element framework on parallel tree-based adaptive meshes

Cited by in corpus (10)

• Finite element interpolated neural networks for solving forward and inverse problems
• Robust and scalable h-adaptive aggregated unfitted finite elements for interface elliptic problems
• Distributed-memory parallelization of the aggregated unfitted finite element method
• The aggregated unfitted finite element method on parallel tree-based adaptive meshes
• Lagrange and $H(\mathrm{curl},\mathcal{B})$ based Finite Element formulations for the relaxed micromorphic model
• Tailored meshing for parallel 3D electromagnetic modeling using high-order edge elements
• A tutorial-driven introduction to the parallel finite element library FEMPAR v1.0.0
• A generic finite element framework on parallel tree-based adaptive meshes
• Compatible finite element interpolated neural networks
• Scalable solvers for complex electromagnetics problems