Ragnar Winther

Douglas N. Arnold, Richard S. Falk, Ragnar Winther

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 Nedelec 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.

Jeonghun J. Lee, Kent-Andre Mardal, Ragnar Winther

Biot's consolidation model in poroelasticity has a number of applications in science, medicine, and engineering. The model depends on various parameters, and in practical applications these parameters ranges over several orders of magnitude. A current challenge is to design discretization techniques and solution algorithms that are well behaved with respect to these variations. The purpose of this paper is to study finite element discretizations of this model and construct block diagonal preconditioners for the discrete Biot systems. The approach taken here is to consider the stability of the problem in non-standard or weighted Hilbert spaces and employ the operator preconditioning approach. We derive preconditioners that are robust with respect to both the variations of the parameters and the mesh refinement. The parameters of interest are small time-step sizes, large bulk and shear moduli, and small hydraulic conductivity.

Douglas N. Arnold, Gerard Awanou, Ragnar Winther

This paper presents a nonconforming finite element approximation of the space of symmetric tensors with square integrable divergence, on tetrahedral meshes. Used for stress approximation together with the full space of piecewise linear vector fields for displacement, this gives a stable mixed finite element method which is shown to be linearly convergent for both the stress and displacement, and which is significantly simpler than any stable conforming mixed finite element method. The method may be viewed as the three-dimensional analogue of a previously developed element in two dimensions. As in that case, a variant of the method is proposed as well, in which the displacement approximation is reduced to piecewise rigid motions and the stress space is reduced accordingly, but the linear convergence is retained.

Trygve Baerland, Jeonghun J. Lee, Kent-Andre Mardal et al.

We discuss the construction of robust preconditioners for finite element approximations of Biot's consolidation model in poroelasticity. More precisely, we study finite element methods based on generalizations of the Hellinger-Reissner principle of linear elasticity, where the stress tensor is one of the unknowns. The Biot model has a number of applications in science, medicine, and engineering. A challenge in many of these applications is that the model parameters range over several orders of magnitude. Therefore, discretization procedures which are well behaved with respect to such variations are needed. The focus of the present paper will be on the construction of preconditioners, such that the preconditioned discrete systems are well-conditioned with respect to variations of the model parameters as well as refinements of the discretization. As a byproduct, we also obtain preconditioners for linear elasticity that are robust in the incompressible limit.

Snorre Harald Christiansen, Ragnar Winther

Eigenvalue problems for semidefinite operators with infinite dimensional kernels appear for instance in electromagnetics. Variational discretizations with edge elements have long been analyzed in terms of a discrete compactness property. As an alternative, we show here how the abstract theory can be developed in terms of a geometric property called the vanishing gap condition. This condition is shown to be equivalent to eigenvalue convergence and intermediate between two different discrete variants of Friedrichs estimates. Next we turn to a more practical means of checking these properties. We introduce a notion of compatible operator and show how the previous conditions are equivalent to the existence of such operators with various convergence properties. In particular the vanishing gap condition is shown to be equivalent to the existence of compatible operators satisfying an Aubin-Nitsche estimate. Finally we give examples demonstrating that the implications not shown to be equivalences, indeed are not.

Kent-Andre Mardal, Joachim Schöberl, Ragnar Winther

A uniform inf-sup condition related to a parameter dependent Stokes problem is established. Such conditions are intimately connected to the construction of uniform preconditioners for the problem, i.e., preconditioners which behave uniformly well with respect to variations in the model parameter as well as the discretization parameter. For the present model, similar results have been derived before, but only by utilizing extra regularity ensured by convexity of the domain. The purpose of this paper is to remove this artificial assumption. As a byproduct of our analysis, in the two dimensional case we also construct a new projection operator for the Taylor-Hood element which is uniformly bounded in $L^2$ and commutes with the divergence operator. This construction is based on a tight connection between a subspace of the Taylor-Hood velocity space and the lowest order Nedelec edge element.

Jeonghun J. Lee, Ragnar Winther

The standard mixed finite element approximations of Hodge Laplace problems associated with the de Rham complex are based on proper discrete subcomplexes. As a consequence, the exterior derivatives, which are local operators, are computed exactly. However, the approximations of the associated coderivatives are nonlocal. In fact, this nonlocal property is an inherent consequence of the mixed formulation of these methods, and can be argued to be an undesired effect of these schemes. As a consequence, it has been argued, at least in special settings, that more local methods may have improved properties. In the present paper, we construct such methods by relying on a careful balance between the choice of finite element spaces, degrees of freedom, and numerical integration rules. Furthermore, we establish key convergence estimates based on a standard approach of variational crimes.

Richard S. Falk, Ragnar Winther

We construct projections from the space of differential k-forms which belong to L2 and whose exterior derivative also belongs to L2, to finite dimensional subspaces of piecewise polynomial differential forms defined on a simplicial mesh. These projections have the properties that they commute with the exterior derivative and are bounded independent of the mesh size. Unlike some other recent work in this direction, the projections are also locally defined in the sense that they are defined by local operators on overlapping macroelements, in the spirit of the Clement interpolant.

Douglas N. Arnold, Richard S. Falk, Ragnar Winther

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 the continuous problem. After a brief introduction to finite element methods, the discretization methods we consider, we develop an abstract Hilbert space framework for analyzing stability and convergence. In this framework, the differential complex is represented by a complex of Hilbert spaces and stability is obtained by transferring Hodge theoretic structures from the continuous level to the discrete. We show stable discretization is achieved if the finite element spaces satisfy two hypotheses: they form a subcomplex and there exists a bounded cochain projection from the full complex to the subcomplex. Next, 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 to the equations of elasticity. Background material is included to make the presentation self-contained for a variety of readers.

Douglas N. Arnold, Gerard Awanou, Ragnar Winther

We construct finite element subspaces of the space of symmetric tensors with square-integrable divergence on a three-dimensional domain. These spaces can be used to approximate the stress field in the classical Hellinger--Reissner mixed formulation of the elasticty equations, when standard discontinous finite element spaces are used to approximate the displacement field. These finite element spaces are defined with respect to an arbitrary simplicial triangulation of the domain, and there is one for each positive value of the polynomial degree used for the displacements. For each degree, these provide a stable finite element discretization. The construction of the spaces is closely tied to discretizations of the elasticity complex, and can be viewed as the three-dimensional analogue of the triangular element family for plane elasticity previously proposed by Arnold and Winther.

Douglas N. Arnold, Richard S. Falk, Ragnar Winther

In this paper, we construct new finite element methods for the approximation of the equations of linear elasticity in three space dimensions that produce direct approximations to both stresses and displacements. The methods are based on a modified form of the Hellinger--Reissner variational principle that only weakly imposes the symmetry condition on the stresses. Although this approach has been previously used by a number of authors, a key new ingredient here is a constructive derivation of the elasticity complex starting from the de Rham complex. By mimicking this construction in the discrete case, we derive new mixed finite elements for elasticity in a systematic manner from known discretizations of the de Rham complex. These elements appear to be simpler than the ones previously derived. For example, we construct stable discretizations which use only piecewise linear elements to approximate the stress field and piecewise constant functions to approximate the displacement field.