Efficient Assembly of H(div) and H(curl) Conforming Finite Elements

Abstract

In this paper, we discuss how to efficiently evaluate and assemble general finite element variational forms on H(div) and H(curl). The proposed strategy relies on a decomposition of the element tensor into a precomputable reference tensor and a mesh-dependent geometry tensor. Two key points must then be considered: the appropriate mapping of basis functions from a reference element, and the orientation of geometrical entities. To address these issues, we extend here a previously presented representation theorem for affinely mapped elements to Piola-mapped elements. We also discuss a simple numbering strategy that removes the need to contend with directions of facet normals and tangents. The result is an automated, efficient, and easy-to-use implementation that allows a user to specify finite element variational forms on H(div) and H(curl) in close to mathematical notation.

References in corpus (7)

• A compiler for variational forms
• Optimisations for quadrature representations of finite element tensors through automated code generation
• Automating the Finite Element Method
• Automated code generation for discontinuous Galerkin methods
• Optimizing the Evaluation of Finite Element Matrices
• Topological Optimization of the Evaluation of Finite Element Matrices
• Benchmarking Domain-Specific Compiler Optimizations for Variational Forms

Cited by in corpus (23)

• Optimisations for quadrature representations of finite element tensors through automated code generation
• A Traveling-Wave Solution for Bacterial Chemotaxis with Growth
• Energy- and enstrophy-conserving schemes for the shallow-water equations, based on mimetic finite elements
• Compatible finite element spaces for geophysical fluid dynamics
• A FEniCS-Based Programming Framework for Modeling Turbulent Flow by the Reynolds-Averaged Navier-Stokes Equations
• A parallel edge orientation algorithm for quadrilateral meshes
• Variational inequality approach to enforce the non-negative constraint for advection-diffusion equations
• On a general implementation of $h$- and $p$-adaptive curl-conforming finite elements
• A Mixed Mimetic Spectral Element Model of the Rotating Shallow Water Equations on the Cubed Sphere
• Compatible-Strain Mixed Finite Element Methods for 3D Compressible and Incompressible Nonlinear Elasticity
• On orienting edges of unstructured two- and three-dimensional meshes
• Petrov-Galerkin flux upwinding for mixed mimetic spectral elements, and its application to geophysical flow problems
• Toward free-surface flow simulations with correct energy evolution: an isogeometric level-set approach with monolithic time-integration
• A mixed variational principle in nonlinear elasticity using Cartan's moving frame and implementation with finite element exterior calculus
• Bringing Trimmed Serendipity Methods to Computational Practice in Firedrake
• Curl-curl conforming elements on tetrahedra
• Design and analysis of an exactly divergence-free hybridized discontinuous Galerkin method for incompressible flows on meshes with quadrilateral cells
• Efficient N-to-M Checkpointing Algorithm for Finite Element Simulations
• Normal-normal continuous symmetric stresses in mixed finite element elasticity
• An evolving space framework for Oseen equations on a moving domain
• A cut finite element method for the Darcy problem
• An hp-hierarchical framework for the finite element exterior calculus
• FIAT: improving performance and accuracy for high-order finite elements