Finite Element Differential Forms Research Articles

Basis for high order divergence-free finite element spaces

A method classically used in the lower polynomial degree for the construction of a finite element basis of the space of divergence-free functions is here extended to any polynomial degree for a bounded domain without topological restrictions. The method uses graphs associated with two differential operators: the gradient and the divergence, and selects the basis using a spanning tree of the first graph. It can be applied for the two main families of degrees of freedom, weights and moments, used to express finite element differential forms.

Symmetry and Invariant Bases in Finite Element Exterior Calculus

We study symmetries of bases and spanning sets in finite element exterior calculus, using representation theory. We want to know which vector-valued finite element spaces have bases invariant under permutation of vertex indices. The permutations of vertex indices correspond to the symmetry group of the simplex. That symmetry group is represented on simplicial finite element spaces by the pullback action. We determine a natural notion of invariance and sufficient conditions on the dimension and polynomial degree for the existence of invariant bases. We conjecture that these conditions are necessary too. We utilize Djoković and Malzan’s classification of monomial irreducible representations of the symmetric group and show new symmetries of the geometric decomposition and canonical isomorphisms of the finite element spaces. Explicit invariant bases with complex coefficients are constructed in dimensions two and three for different spaces of finite element differential forms.

High order approximation of Hodge Laplace problems with local coderivatives on cubical meshes

In mixed finite element approximations of Hodge Laplace problems associated with the de Rham complex, the exterior derivative operators are computed exactly, so the spatial locality is preserved. However, the numerical approximations of the associated coderivatives are nonlocal and can be regarded as an undesired effect of standard mixed methods. For numerical methods with local coderivatives, a perturbation of low order mixed methods in the sense of variational crimes has been developed for simplicial and cubical meshes. In this paper we extend the low order method to all high orders on cubical meshes using a new family of finite element differential forms on cubical meshes. The key theoretical contribution is a generalization of the linear degree, in the construction of the serendipity family of differential forms, and this generalization is essential in the unisolvency proof of the new family of finite element differential forms.

Local finite element approximation of Sobolev differential forms

We address fundamental aspects in the approximation theory of vector-valued finite element methods, using finite element exterior calculus as a unifying framework. We generalize the Clément interpolant and the Scott-Zhang interpolant to finite element differential forms, and we derive a broken Bramble-Hilbert lemma. Our interpolants require only minimal smoothness assumptions and respect partial boundary conditions. This permits us to state local error estimates in terms of the mesh size. Our theoretical results apply to curl-conforming and divergence-conforming finite element methods over simplicial triangulations.

Charge-conserving hybrid methods for the Yang–Mills equations

The Yang-Mills equations generalize Maxwell's equations to nonabelian gauge groups, and a quantity analogous to charge is locally conserved by the nonlinear time evolution. Christiansen and Winther observed that, in the nonabelian case, the Galerkin method with Lie algebra-valued finite element differential forms appears to conserve charge globally but not locally, not even in a weak sense. We introduce a new hybridization of this method, give an alternative expression for the numerical charge in terms of the hybrid variables, and show that a local, per-element charge conservation law automatically holds.

Generalized Gaffney inequality and discrete compactness for discrete differential forms

We prove generalized Gaffney inequalities and the discrete compactness for finite element differential forms on s-regular domains, including general Lipschitz domains. In computational electromagnetism, special cases of these results have been established for edge elements with weakly imposed divergence-free conditions and used in the analysis of nonlinear and eigenvalue problems. In this paper, we generalize these results to discrete differential forms, not necessarily with strongly or weakly imposed constraints. The analysis relies on a new Hodge mapping and its approximation property. As an application, we show $$L^{p}$$ estimates for several finite element approximations of the scalar and vector Laplacian problems.

Trimmed serendipity finite element differential forms

We introduce the family of trimmed serendipity finite element differential form spaces, defined on cubical meshes in any number of dimensions, for any polynomial degree, and for any form order. The relation between the trimmed serendipity family and the (non-trimmed) serendipity family developed by Arnold and Awanou [Math. Comp. 83(288) 2014] is analogous to the relation between the trimmed and (non-trimmed) polynomial finite element differential form families on simplicial meshes from finite element exterior calculus. We provide degrees of freedom in the general setting and prove that they are unisolvent for the trimmed serendipity spaces. The sequence of trimmed serendipity spaces with a fixed polynomial order r provides an explicit example of a system described by Christiansen and Gillette [ESAIM:M2AN 50(3) 2016], namely, a minimal compatible finite element system on squares or cubes containing order r-1 polynomial differential forms.

Smoothed projections over weakly Lipschitz domains

We develop finite element exterior calculus over weakly Lipschitz domains. Specifically, we construct commuting projections from $L^p$ de Rham complexes over weakly Lipschitz domains onto finite element de Rham complexes. The projections satisfy uniform bounds for finite element spaces with bounded polynomial degree over shape-regular families of triangulations. Thus we extend the theory of finite element differential forms to polyhedral domains that are weakly Lipschitz but not strongly Lipschitz. As new mathematical tools, we use the collar theorem in the Lipschitz category, and we show that the degrees of freedom in finite element exterior calculus are flat chains in the sense of geometric measure theory.

Finite element exterior calculus for parabolic problems

In this paper, we consider the extension of the finite element exterior calculus from elliptic problems, in which the Hodge Laplacian is an appropriate model problem, to parabolic problems, for which we take the Hodge heat equation as our model problem. The numerical method we study is a Galerkin method based on a mixed variational formulation and using as subspaces the same spaces of finite element differential forms which are used for elliptic problems. We analyze both the semidiscrete and a fully-discrete numerical scheme.

Finite element differential forms on curvilinear cubic meshes and their approximation properties

We study the approximation properties of a wide class of finite element differential forms on curvilinear cubic meshes in n dimensions. Specifically, we consider meshes in which each element is the image of a cubical reference element under a diffeomorphism, and finite element spaces in which the shape functions and degrees of freedom are obtained from the reference element by pullback of differential forms. In the case where the diffeomorphisms from the reference element are all affine, i.e., mesh consists of parallelotopes, it is standard that the rate of convergence in L2 exceeds by one the degree of the largest full polynomial space contained in the reference space of shape functions. When the diffeomorphism is multilinear, the rate of convergence for the same space of reference shape function may degrade severely, the more so when the form degree is larger. The main result of the paper gives a sufficient condition on the reference shape functions to obtain a given rate of convergence.

Finite element differential forms on cubical meshes

We develop a family of finite element spaces of differential forms defined on cubical meshes in any number of dimensions. The family contains elements of all polynomial degrees and all form degrees. In two dimensions, these include the serendipity finite elements and the rectangular BDM elements. In three dimensions they include a recent generalization of the serendipity spaces, and new H ( c u r l ) H(\mathrm {curl}) and H ( d i v ) H(\mathrm {div}) finite element spaces. Spaces in the family can be combined to give finite element subcomplexes of the de Rham complex which satisfy the basic hypotheses of the finite element exterior calculus, and hence can be used for stable discretization of a variety of problems. The construction and properties of the spaces are established in a uniform manner using finite element exterior calculus.

Finite element exterior calculus: from Hodge theory to numerical stability

This article reports on the confluence of two streams of research, one emanating from the fields of numerical analysis and scientific computation, the other from topology and geometry. In it we consider the numerical discretization of partial differential equations that are related to differential complexes so that de Rham cohomology and Hodge theory are key tools for exploring the well-posedness of the continuous problem. The discretization methods we consider are finite element methods, in which a variational or weak formulation of the PDE problem is approximated by restricting the trial subspace to an appropriately constructed piecewise polynomial subspace. After a brief introduction to finite element methods, we develop an abstract Hilbert space framework for analyzing the stability and convergence of such discretizations. In this framework, the differential complex is represented by a complex of Hilbert spaces, and stability is obtained by transferring Hodge-theoretic structures that ensure well-posedness of the continuous problem from the continuous level to the discrete. We show stable discretization is achieved if the finite element spaces satisfy two hypotheses: they can be arranged into a subcomplex of this Hilbert complex, and there exists a bounded cochain projection from that complex to the subcomplex. In the next part of the paper, we consider the most canonical example of the abstract theory, in which the Hilbert complex is the de Rham complex of a domain in Euclidean space. We use the Koszul complex to construct two families of finite element differential forms, show that these can be arranged in subcomplexes of the de Rham complex in numerous ways, and for each construct a bounded cochain projection. The abstract theory therefore applies to give the stability and convergence of finite element approximations of the Hodge Laplacian. Other applications are considered as well, especially the elasticity complex and its application to the equations of elasticity. Background material is included to make the presentation self-contained for a variety of readers.

Geometric decompositions and local bases for spaces of finite element differential forms

We study the two primary families of spaces of finite element differential forms with respect to a simplicial mesh in any number of space dimensions. These spaces are generalizations of the classical finite element spaces for vector fields, frequently referred to as Raviart–Thomas, Brezzi–Douglas–Marini, and Nédélec spaces. In the present paper, we derive geometric decompositions of these spaces which lead directly to explicit local bases for them, generalizing the Bernstein basis for ordinary Lagrange finite elements. The approach applies to both families of finite element spaces, for arbitrary polynomial degree, arbitrary order of the differential forms, and an arbitrary simplicial triangulation in any number of space dimensions. A prominent role in the construction is played by the notion of a consistent family of extension operators, which expresses in an abstract framework a sufficient condition for deriving a geometric decomposition of a finite element space leading to a local basis.

Finite element differential forms

AbstractA differential form is a field which assigns to each point of a domain an alternating multilinear form on its tangent space. The exterior derivative operation, which maps differential forms to differential forms of the next higher order, unifies the basic first order differential operators of calculus, and is a building block for a great variety of differential equations. When discretizing such differential equations by finite element methods, stable discretization depends on the development of spaces of finite element differential forms. As revealed recently through the finite element exterior calculus, for each order of differential form, there are two natural families of finite element subspaces associated to a simplicial triangulation. In the case of forms of order zero, which are simply functions, these two families reduce to one, which is simply the well‐known family of Lagrange finite element subspaces of the first order Sobolev space. For forms of degree 1 and of degree n − 1 (where n is the space dimension), we obtain two natural families of finite element subspaces, unifying many of the known mixed finite element spaces developed over the last decades. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)