Finite element exterior calculus: from Hodge theory to numerical stability

Abstract

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.

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

“…The setting of finite element exterior calculus (see [1,7,8]) provides a unified framework for the study of this problem. In this paper we discuss the construction of finite element subspaces of the domain H Λ k of the exterior derivative acting on differential k-forms, 0 ≤ k ≤ n, in any number n of dimensions.…”

Symmetric Matrix Fields in the Finite Element Method

“…The construction of discrete versions of the elasticity complex has been made further systematic in [20,45,46] where a procedure to construct discrete versions of the elasticity complex from discrete versions of the de Rham complex, an exact sequence connecting spaces of differential forms, was made explicit. This process has been further improved in the finite element exterior calculus [20,48]. Finally, we note that despite their relative complexity, conforming mixed finite elements with symmetric stress field are useful in certain situations [49].…”

Effect of ocular shape and vascular geometry on retinal hemodynamics: a computational model

“…For computational purposes exterior calculus has been discretized as finite element exterior calculus [3] and Discrete Exterior Calculus [13,25]. These discretizations are useful either when a mixed method (involving both velocities and pressures) is to be implemented or when the domain is not flat, as is the case in this paper.…”