Let be a finite element space of piecewise polynomials of degree subject to a quasi-uniform tetrahedral partition of a convex polyhedron , where denotes the mesh size of the tetrahedral partition. Let be the subspace of consisting of functions with zero boundary values.
A function is called discrete harmonic if it satisfies the following equation:
In this article, we establish the following result, called weak maximum principle of finite element methods.
A discrete harmonic function satisfies the following estimate:
satisfies the following estimate:
where the constant is independent of the mesh size .
As an application of the weak maximum principle, we show that the Ritz projection defined by
is stable in for finite elements of degree , i.e.
In the finite element literature, the maximum principle has attracted a lot of attention; see [7, 8, 24, 29, 30], to mention a few. However, the sufficient conditions for discrete maximum principle put serious restrictions on the geometry of the mesh. For piecewise linear elements in two-dimensions, the angles of the triangles must be less than , or the sum of opposite angles of the triangles that share an edge must be less than (for example, see [30, §5]). For quadratic elements in two dimensions, discrete maximum principe holds only for equilateral triangles . The situation in three dimensions is more complicated [4, 17, 18, 31], essentially it is hard to guarantee the discrete maximum principe even for piecewise linear elements. In this respect, there stands out the work of Schatz , who proved that a weak maximum principle in the sense of (1.2) holds for a wide class of finite elements on general quasi-uniform triangulation of any two dimensional polygonal domain . By utilizing the weak maximum principle, Schatz also established the stability of the Ritz projection in and norms. Such stability results have a wide range of applications, for example to pointwise error estimates of finite element methods for parabolic problems [21, 16, 20], Stokes systems , nonlinear problems [11, 10, 22], obstacle problems , optimal control problems [1, 2], to name a few. As far as we know,  is the only paper that establishes weak maximum principle and stability estimate (without the logarithmic factor) for the Ritz projection on nonsmooth domains.
In three dimensions the situation is less satisfactory. The stability of the Ritz projection in and norms are available on smooth domains [26, 28] and convex polyhedral domains [13, 19]. However, on convex polyhedral domains in , the -stability constant depends logarithmically on the mesh size , and it is not obvious how the logarithmic factor can be removed there. There are no results on the weak maximum principles in three dimensions even on smooth domains or convex polyhedra. The objective of this paper is to close this gap for convex polyhedral domains. In order to obtain the result, we have to modify the argument in  by considering error analysis in norm for some . In the case of convex polyhedral domains, the -norm based argument used in  would yield an additional logarithmic factor. Unfortunately, the current analysis does not allow us to extend the results to nonconvex polyhedral domains or graded meshes. These would be the subject of future research.
The paper is organized as follows. In section 2 we state some preliminary results that we use later in our arguments. In section 3, we reduce the proof of the weak discrete maximum principle to a specific error estimate. Section 4 is devoted to the proof of this estimate, which constitutes the main technical part of the paper. Finally, section 5, gives an application of the weak discrete maximum principle to showing the stability of the Ritz projection in norm uniformly in for higher order elements.
In the rest of this article, we denote by a generic positive constant, which may be different at different occurrences but will be independent of the mesh size .
2. Preliminary results
In this section, we present several well-known results that are used in our analysis. First result concerns global regularity of the weak solution to the problem
Let be a convex polyhedron. Then there exists a constant depending on such that for any and , the solution of (2.3) is in and
The next result addresses the problem (2.3) when the source function is supported in some part of . The following lemma traces the dependence of the stability constant on the diameter of the support.
For any bounded Lipschitz domain , there exist positive constants and (depending on ) such that for and such that with , the solution of (2.3) satisfies
For any there holds
If then we let be the weak solution of
The solution defined above satisfies
and, according to [15, Theorem B], there exists a positive constant such that
By using these properties, we have
The duality pairing estimate above implies the desired result. ∎
The next lemma concerns basic properties of harmonic functions on convex domains. The result is essentially the same as in [27, Lemma 8.3].
Let and be two subdomains satisfying , with
where is a positive constant. If and is harmonic on , i.e.
then the following estimates hold:
Finally, we need the best approximation property of the Ritz projection in norm. In , the best approximation property of the Ritz projection in norm was established on convex polyhedral domains. Together with the standard best approximation property in norm we obtain
for any . Extension of the above result to follows by duality (cf. [5, §8.5]). These can be summarized as below.
On a convex polyhedron , the following estimate holds for any fixed :
3. Basic estimates
In [26, Corollary 5.1], the following interior error estimate was established
for , where for , for and , with and . Choosing , and in the above estimate, we obtain that there exists a constant independent of such that
Let be a point satisfying
Otherwise, we have . In this case, the inverse inequality of finite element functions implies
Hence, either or , the following estimate holds:
To estimate the term on the right hand side of the inequality above, we use the following duality property:
which implies the existence of a function with the following properties:
For this function , we define to be the solution of the PDE problem (in weak form):
and let be the finite element solution of
Then is the Ritz projection of , satisfying
Let be the solution of the PDE problem (in weak form)
Then the continuous maximum principle of the PDE problem implies
Notice, that is the Ritz projection of , i.e.
Therefore, we have
where we have used (3.8) and the Hölder inequality in deriving the last inequality.
4. Estimate of
Let and for We define a sequence of subdomains
For each we denote to be a subdomain slightly bigger than , defined by
Let , with denoting the greatest integer not exceeding . Then
By using these subdomains defined above, we have
where the Hölder inequality and (4.13) were used in deriving the last inequality.
where we have used and in deriving the last inequality. Substituting the last inequality into (4) yields
where we have used and the following inequality in deriving the last inequality:
The inequality above follows from Lemma 2.3, the Hölder inequality and Sobolev embedding, i.e.
This proves that (4) holds for .
where the last inequality is due to the following Hölder inequality:
From (4) we see that
It remains to estimate . To this end, we let be a smooth cut-off function satisfying
By using (4
) and the interpolation inequality (for)
where can be an arbitrary positive number. By choosing with , we obtain
where we have used (4.19) in deriving the last inequality. Note that
Combining the last two estimates, we obtain
If for sufficiently large constant , then the last term can be absorbed by the left side. Hence, we have