A new discretization for mth-Laplace equationswith arbitrary polynomial degreesThis work was supported by the Berlin Mathematical School.

A new discretization for th-Laplace equations
with arbitrary polynomial degreesthanks: This work was supported by the Berlin Mathematical School.

M. Schedensack Institut für Numerische Simulation, Universität Bonn, Wegelerstraße 6, D-53115 Bonn, Germany

This paper introduces new mixed formulations and discretizations for th-Laplace equations of the form for arbitrary based on novel Helmholtz-type decompositions for tensor-valued functions. The new discretizations allow for ansatz spaces of arbitrary polynomial degree and the lowest-order choice coincides with the non-conforming FEMs of Crouzeix and Raviart for and of Morley for . Since the derivatives are directly approximated, the lowest-order discretizations consist of piecewise affine and piecewise constant functions for any Moreover, a uniform implementation for arbitrary is possible. Besides the a priori and a posteriori analysis, this paper proves optimal convergence rates for adaptive algorithms for the new discretizations.

Keywords th-Laplace equation, polyharmonic equation, non-conforming FEM, mixed FEM, adaptive FEM, optimality

AMS subject classification 31A30, 35J30, 65N30, 65N12, 74K20

1 Introduction

This paper considers th-Laplace equations of the form


for arbitrary Standard conforming FEMs require ansatz spaces in . To circumvent those high regularity requirements and resulting complicated finite elements, non-standard methods are of high interest [Mor68, EGH02, Bre12, GN11]. The novel Helmholtz decomposition of this paper decomposes any (tensor-valued) function in an th derivative and a symmetric part of a Curl. Given a tensor-valued function which satisfies in the weak sense, the projection of to the space of th derivatives then coincides with the th derivative of the exact solution of (1.1) (see Theorem 5.1 below). This results in novel mixed formulations and discretizations for (1.1). This approach generalises the discretizations of [Sch15, Sch16] from to .

The direct approximation of instead of enables low order discretizations; only first derivatives appear in the symmetric part of the Curl and so the lowest order approach only requires piecewise affine functions for any . In contrast to that, even interior penalty methods require piecewise quadratic [Bre12] resp. piecewise cubic [GN11] functions for resp. . Mnemonic diagrams in Figure 1 illustrate lowest-order standard conforming FEMs from [Žen70] and the lowest-order novel FEMs proposed in this work for . Since the proposed new FEMs differ only in the number of components in the ansatz spaces, an implementation of one single program, which runs for arbitrary order, is possible. In particular, the system matrices are obtained by integration of standard FEM basis functions.

Figure 1: Lowest order standard conforming [Cia78, Žen70] and novel FEMs for the problem for .

For and the lowest polynomial degree in the ansatz spaces, discrete Helmholtz decompositions of [AF89, CGH14] prove that the discrete solutions are piecewise gradients (resp. Hessians) of Crouzeix-Raviart [CR73] (resp. Morley [Mor68]) finite element functions and therefore the new discretizations can be regarded as a generalization of those non-conforming FEMs to higher polynomial degrees and higher-order problems. The generalization of [WX13] of the non-conforming Crouzeix-Raviart and Morley FEMs to is restricted to a space dimension .

In the context of the novel (mixed) formulations, the discretizations appear to be conforming. The new generalization to higher polynomial degrees proposed in this paper appears to be natural in the sense that the inherent properties of the lowest order discretization carry over to higher polynomial ansatz spaces, namely an inf-sup condition, the conformity of the method, and a crucial projection property (also known as integral mean property of the non-conforming interpolation operator).

Besides the a priori and a posteriori error analysis, this paper proves optimal convergence rates for an adaptive algorithm, which are also observed in the numerical experiments from Section 7.

The remaining parts of this paper are organised as follows. Section 2 introduces some notation while some preliminary results are proved in Section 3. The proposed discretization of (1.1) in Section 5 is based on a novel Helmholtz decomposition for higher derivatives which is stated and proved in Section 4. Section 6 introduces an adaptive algorithm and proves optimal convergence rates. Section 7 concludes the paper with numerical experiments on fourth- and sixth-order problems.

Throughout this paper, let be a bounded, polygonal, simply connected Lipschitz domain. Standard notation on Lebesgue and Sobolev spaces and their norms is employed with scalar product . Given a Hilbert space , let resp.  denote the space of functions with values in whose components are in resp. . The space of infinitely differentiable functions reads and the subspace of functions with compact support in is denoted with . The piecewise action of differential operators is denoted with a subscript . The formula represents an inequality for some mesh-size independent, positive generic constant ; abbreviates . By convention, all generic constants do neither depend on the mesh-size nor on the level of a triangulation but may depend on the fixed coarse triangulation and its interior angles.

2 Notation

This section introduces notation related to higher-order tensors and tensor-valued functions and triangulations.

Define the set of -tensors over by

and let denote the symmetric group, i.e., the set of all permutations of . Define the set of symmetric tensors by

The symmetric part of a tensor is defined by

for all , where denotes the number of elements in a set . For , the set coincides with the set of symmetric matrices, while for , the tensors consist of the four different components , , , and . Given -tensors and a vector , define the scalar product and the dot product by

for all . The following definition summarizes some differential operators. Recall that, for a Hilbert space , the space (resp. ) denotes the space of (resp. ) functions with components in .

Definition 1 (differential operators).

Let and and define by and . Define the th derivative of , the derivative , the divergence , the Curl, , and the curl, by

for .

For , these definitions coincide with the row-wise application of , , , and . The scalar product of tensor-valued functions is defined by . Given such that there exists with

define the th order divergence of . The space is defined by

Define furthermore for

Remark 2.1.

Note that the existence of the th weak divergence does not imply the existence of any -th divergence for , e.g., for .

A shape-regular triangulation of a bounded, polygonal, open Lipschitz domain is a set of closed triangles such that and any two distinct triangles are either disjoint or share exactly one common edge or one vertex. Let denote the edges of a triangle and the set of edges in . Any edge is associated with a fixed orientation of the unit normal on (and denotes the unit tangent on ). On the boundary, is the outer unit normal of , while for interior edges , the orientation is fixed through the choice of the triangles and with and is the outer normal of on . In this situation, denotes the jump across . For an edge on the boundary, the jump across reads . For and , let

denote the set of piecewise polynomials and . Given a subspace , let denote the projection onto and let abbreviate . Given a triangle , let denote the square root of the area of and let denote the piecewise constant mesh-size with for all . For a set of triangles , let abbreviate

3 Results for tensor-valued functions

The main result of this section is Theorem 3.2, which proves that defines a norm on the space defined in (3.5) below and can, thus, be viewed as a generalized Korn inequality. The following theorem is used in the proof of Theorem 3.2. Recall the definition of the Curl and the symmetric part of a tensor from Section 2.

Theorem 3.1.

Any satisfies


The proof is subdivided in three steps.

Step 1. Let and with for all and for all , i.e.,

The combination of the definitions of and reads


Let be the multi-index with the same number of ones and the number of twos reduced by one and the multi-index with the same number of twos and the number of ones reduced by one, i.e.,

The symmetry of implies that if and if . Since the number of permutations such that is and the number of permutations such that is and since and , this implies that (3.1) equals


Step 2. This step applies [Neč67, Chap. 3, Thm. 7.6] and [Neč67, Chap. 3, Thm. 7.8] to operators defined below. Step 3 then proves a relation between these operators and the operator .

Define for , , and a multi-index

Furthermore, define for

with the multi-index notation . Then the matrix reads

If , the columns of this matrix are linear independent. Define the operators , , by

Then, the combination of [Neč67, Chap. 3, Thm. 7.6] with [Neč67, Chap. 3, Thm. 7.8] proves


Step 3. This step proves a relation between and for a proper choice of .

Define by setting for each the function with and (with for and for ). The symmetry of proves


With the notation from Step 1 it holds that and the definition of from Step 2 and (3.2) reveal

This leads to

This, (3.4), and an application of (3.3) implies the assertion. ∎

Define, for , the spaces


A computation reveals for , that the spaces and read


and for the space reads


The following theorem generalizes [CGH14, Lemma 3.3] from to higher-order tensors and states that defines a norm on . Note that .

Theorem 3.2.

Any satisfies


Assume for contradiction that the statement does not hold. Then there exists a sequence with

Since , Poincaré’s inequality implies that all components of are bounded in . Since is reflexive and compactly embedded in , there exists a subsequence (not relabelled) with a limit , in . This and Theorem 3.1 imply

The Poincaré inequality and the completeness of imply the existence of with in and thus . It holds that and, therefore, defines a bounded functional on . Hence,


Let . Since , the Cauchy inequality reveals

This and (3.8) lead to and therefore . This contradicts and, hence, implies the assertion. ∎

Remark 3.3 (dependency on the domain).

The proof by contradiction from Theorem 3.2 does not provide information about the dependency on the domain. A scaling argument reveals that it does not depend on the size of the domain, but it may depend on its shape.

4 Helmholtz decomposition for higher orders

This section proves a Helmholtz decomposition of tensors into th derivatives and the symmetric part of a Curl in Theorem 4.4. This is a generalization of the Helmholtz decomposition of [BNS07] for fourth-order problems (). The proof is based on Theorem 4.1 below, which characterizes th-divergence-free smooth functions as symmetric parts of Curls.

Theorem 4.1.

Let and with . Then there exists with


The proof is based on mathematical induction.

The base case is a classical result [Rud76]. Assume as induction hypothesis that the statement holds for , i.e., for all with there exists with .

The inductive step is split in five steps. Suppose that with .

Step 1. Then and . Let and . Recall the definition of the divergence from Definition 1. The symmetry of implies

Hence, . The induction hypothesis guarantees the existence of with .

Step 2. This step defines some with .

The definitions of and from Section 2 for tensors combine to


Define by


The definition of implies

Since if and only if , this equals

and, hence, the combination of the foregoing two displayed formulae with (4.1) leads to . The combination with Step 1 proves .

Step 3. Since , the base case (applied “row-wise” to ) guarantees the existence of with .

Step 4. This step shows .

Let be fixed and let and be the number of ones and twos. Then


are the numbers of ones and twos in . Define the index set