1 Introduction
Equilibrated a posteriori error estimators have attracted much interest recently due to the guaranteed reliability bound with the reliability constant being one. This property implies that they are perfect for discretization error control on both coarse and fine meshes. Error control on coarse meshes is important but difficult for computationally challenging problems.
For the conforming finite element approximation, a mathematical foundation of equilibrated estimators is the Prager-Synge identity [31] that is valid in (see Section 3). Based on this identity, various equilibrated estimators have been studied recently by many researchers (see, e.g., [27, 22, 29, 20, 21, 6, 3, 33, 10, 12, 13, 34, 17, 14]). The key ingredient of the equilibrated estimators for the continuous elements is local recovery of an equilibrated (locally conservative) flux in the space through the numerical flux. By using a partition of unity, Ladevèze and Leguillon [27] initiated a local procedure to reduce the construction of an equilibrated flux to vertex patch based local calculations. For the continuous linear finite element approximation to the Poisson equation in two dimensions, an equilibrated flux in the lowest order Raviart-Thomas space was explicitly constructed in [10, 12]. This explicit approach does not lead to robust equilibrated estimator with respect to the coefficient jump without introducing a constraint minimization (see [17]). The constraint minimization on each vertex patch may be efficiently solved by first computing an equilibrated flux and then calculating a divergence free correction. For recent developments, see [14] and references therein.
The purpose of this paper is to develop and analyze equilibrated a posteriori error estimators for discontinuous elements including both nonconforming and discontinuous Galerkin elements. To do so, the first and the essential step is to extend the Prager-Synge identity to be valid for piecewise functions. This will be done by establishing a generalized Prager-Synge inequality (see Theorem 3.1) that contains an additional term measuring the distance between and piecewise . Moreover, by using a Helmholtz decomposition, we will be able to show that the inequality becomes an identity in two dimensions (see Lemma 3.4). A non-optimal inequality similar to ours was obtained earlier by Braess, Fraunholz, and Hoppe in [11] for the Poisson equation with pure Dirichlet boundary condition. Based on the generalized Prager-Synge inequality and an equivalent form (see Corollary 3.2), the construction of an equilibrated a posteriori error estimator for discontinuous finite element solutions is reduced to recover an equilibrated flux in and to recover either a potential function in
or a curl free vector-valued function in
.Recovery of equilibrated fluxes for discontinuous elements has been studied by many researchers. For discontinuous Garlerkin (DG) methods, equilibrated fluxes in Raviart-Thomas (RT) spaces were explicitly reconstructed in [2] for linear elements and in [23] for higher order elements. For nonconforming finite element methods, existing explicit equilibrated flux recoveries in RT spaces seem to be limited to the linear Crouzeix-Raviart (CR) and the quadratic Fortin-Soulie elements by Marini [28] (see [1] in the context of estimator) and Kim [26], respectively. For higher order nonconforming elements, a local reconstruction procedure was proposed by Ainsworth and Rankin in [4] through solving element-wise minimization problems. The recovered flux is not in the RT spaces. Nevertheless, the resulting estimator provides a guaranteed upper bound.
In this paper, we will introduce a fully explicit post-processing procedure for recovering an equilibrated flux in the RT space of index for the nonconforming elements of odd order of . Currently, we are not able to extend our recovery technique to even orders. This is because our recovery procedure heavily depends on the finite element formulation and the properties of the nonconforming finite element space; moreover, structure of the nonconforming finite element spaces of even and odd orders are fundamentally different.
Recovery of a potential function in for discontinuous elements was studied by some researchers (see, e.g., [4, 2, 11]). Local approaches for recovering equilibrated flux in [10, 12, 17, 13, 14] may be directly applied (at least in two dimensions) for computing an approximation to the gradient in the curl-free space. (As mentioned previously, this approach requires solutions of local constraint minimization problems over vertex patches.) The resulting a posteriori error estimator from either the potential or the gradient recoveries may be proved to be locally efficient. Nevertheless, to show independence of the efficiency constant on the jump, we have to assume that the distribution of the diffusion coefficient is quasi-monotone (see [30]).
In this paper, we will employ a simple averaging technique over edges to recover a gradient in . Due to the fact that the recovered gradient is not necessarily curl free, the reliability constant of the resulting estimator is no longer one. However, it turns out that the curl free constraint is not essential and, theoretically we are able to prove that the resulting estimator has the robust local reliability as well as the robust local efficiency without the quad-monotone assumption. This is compatible with our recent result in [15] on the residual error estimator for discontinuous elements.
This paper is organized as follows. The diffusion problem and the finite element mesh are introduced in Section 2. The generalized Prager-Synge inequality for piecewise functions are established in Section 3. Explicit recoveries of an equilibrated flux and a gradient and the resulting a posteriori error estimator for discontinuous elements are described in Section 4. Global reliability and local efficiency of the estimator are proved in Section 5. Finally, numerical results are presented in Section 6.
2 Model problem
Let be a bounded polygonal domain in , with Lipschitz boundary , where . For simplicity, assume that . Considering the diffusion problem:
(2.1) |
with boundary conditions
where and are the respective divergence and gradient operators; is the outward unit vector normal to the boundary; and are given scalar-valued functions; and the diffusion coefficient
is symmetric, positive definite, and piecewise constant full tensor with respect to the domain
. Here we assume that the subdomain, for , is open and polygonal.We use the standard notations and definitions for the Sobolev spaces. Let
Then the corresponding variational problem of (2.1) is to find such that
(2.2) |
where is the inner product on the domain . The subscript is omitted when .
2.1 Triangulation
Let be a finite element partition of that is regular, and denote by the diameter of the element . Furthermore, assume that the interfaces,
do not cut through any element . Denote the set of all edges of the triangulation by
where is the set of interior element edges, and and are the sets of boundary edges belonging to the respective and . For each , denote by the length of and by a unit vector normal to . Let and be the two elements sharing the common edge such that the unit outward normal of coincides with . When , is the unit outward normal to and denote by the element having the edge .
3 Generalized Prager-Synge inequality
For the conforming finite element approximation, the foundation of the equilibrated a posteriori error estimator is the Prager-Synge identity [31]. That is, let be the solution of (2.1), then
for all and for all , where is the so-called equilibrated flux space defined by
Here, denotes the space of all vector-valued functions whose divergence are in . The Prager-Synge identity immediately leads to
(3.1) |
Choosing to be the conforming finite element approximation, then (3.1) implies that
(3.2) |
is a reliable estimator with the reliability constant being one.
We now proceed to establish a generalization of (3.1) for piecewise functions with applications to nonconforming and discontinuous Galerkin finite element approximations. To this end, denote the broken space with respect to by
Define be the discrete gradient operator on such that for any
Theorem 3.1.
Let be the solution of (2.1). In both two and three dimensions, for all , we have
(3.3) |
Proof.
A suboptimal result for the Poisson equation () with pure Dirichlet boundary condition is proved in [11] by Braess, Fraunholz, and Hoppe:
Let be the space of all vector-valued functions whose curl are in , and denote its curl free subspace by
where denotes the tangent vector(s).
Corollary 3.2.
Let be the solution of (2.1). In both two and three dimensions, for all , we have
(3.5) |
In the remaining section, we prove that, in two dimensions, the inequality (3.3) in Theorem 3.1 is indeed an equality. For each , in two dimensions, assume that , then denote by the unit vector tangent to and by and the start and end points of , respectively, such that . Let
For a vector-valued function , define the curl operator by
For a scalar-valued function , define the formal adjoint operator of the curl by
For a fixed , there exist unique and for the following Helmholtz decomposition (see, e.g., [4]) such that
(3.6) |
and and satisfy
and
respectively. It is easy to see that and are orthogonal with respect to the inner product, which yields
(3.7) |
Lemma 3.3.
Let be a fixed function in and and be the corresponding Helmholtz decomposition of given in (3.6). We have
(3.8) |
Proof.
Lemma 3.4.
Let be the solution of (2.1). In two dimensions, for all , we have
(3.9) |
Remark 3.5.
For each , denote by and
the maximal and minimal eigenvalues of
, respectively. For each , let , , and if and if . To this end, letAssume that each local matrix
is similar to the identity matrix in the sense that its maximal and minimal eigenvalues are almost of the same size. More precisely, there exists a moderate size constant
such thatNevertheless, the ratio of global maximal and minimal eigenvalues, , is allowed to be very large.
For a function , denote its traces on by and and the jump of across the edge by
In the following lemma, we show the relationship between the nonconforming error and the residual based error of solution jump on edges. It is noted that the constant is robust with respect to the coefficient jump.
Lemma 3.6.
Let be a fixed function in . In two dimensions, there exists a constant that is independent of the jump of the coefficient such that
(3.10) |
Proof.
4 Error estimators and indicators
4.1 NC finite element approximation
For the convenience of readers, in this subsection we introduce the nonconforming finite element space and its properties.
Let and be the spaces of polynomials of degree less than or equal to on the element and , respectively. Define the nonconforming finite element space of order on the triangulation by
(4.1) |
and its subspace by
The spaces defined above are exactly the same as those defined in [19] for , [24] for , [18] for and , [4] for general odd order, and [32, 5] for general order. Then the nonconforming finite element approximation of order is to find such that
(4.2) |
Below we describe basis functions of and their properties. To this end, for each , let for and for . Denote by the set of all interior Lagrange points in with respect to the space and by the nodal basis function corresponding to , i.e.,
where is the Kronecker delta function. For each , let be the th order Gauss-Legendre polynomial on such that . Note that is an odd or even function when is odd or even. Hence, for odd and for even .
For odd
, the set of degrees of freedom of
(see Lemma 2.1 in [4]) can be given by(4.3) |
for all and
(4.4) |
for all . Define the basis function satisfying
(4.5) |
for and , and the basis function satisfying
(4.6) |
for and . Then the nonconforming finite element space is the space spanned by all these basis functions, i.e.,
Lemma 4.1.
For all , the basis functions have support on and vanish on the boundary of , i.e.,
Proof.
For each , denote by the set of all edges of . For each , denote by the union of all elements that share the common edge ; and define a sign function on the set (when is a boundary edge, let ) such that
Lemma 4.2.
For all , the basis functions have support on , and their restrictions on has the following representation:
(4.7) |
when is odd, and
(4.8) |
when is even.
Proof.
By (4.6), it is easy to see that support of is . Since , there exist constants such that
Using (4.6) and the orthogonality of , it is obvious that
and, hence,
(4.9) |
By (4.6), it is also easy to see that there exists constant for each such that
(4.10) |
Since is an odd function for all and is continuous in and , (4.10) implies that
(4.11) |
Combining the facts that for odd and that for even , (4.9), and (4.11), we have
which, together with (4.9), leads to the formulas of in (4.7) and (4.8). Finally, for each , in (4.10) can be directly computed based on the continuity of in and . This completes the proof of the lemma.
∎
Remark 4.3.
As a consequence of Lemma 4.2, the basis function is continuous on the edge , i.e., for all ; moreover, vanishes at end points of , i.e., , for odd .
Lemma 4.4.
Let be an edge of . Assume that . Then we have that
(4.12) |
Moreover, if for all , then on .
4.2 Equilibrated flux recovery
In this subsection, we introduce a fully explicit post-processing procedure for recovering an equilibrated flux. To this end, define by
where is the projection onto . For simplicity, assume that the Neumann data is a piecewise polynomial of degree less than or equal to , i.e., for all .
Denote the conforming Raviart-Thomas (RT) space of index with respect to by
where . Let
Comments
There are no comments yet.