Static two-grid mixed finite-element approximations to the Navier-Stokes equations
A two-grid scheme based on mixed finite-element approximations to the incompressible Navier-Stokes equations is introduced and analyzed. In the first level the standard mixed finite-element approximation over a coarse mesh is computed. In the second level the approximation is postprocessed by solving a discrete Oseen-type problem on a finer mesh. The two-level method is optimal in the sense that, when a suitable value of the coarse mesh diameter is chosen, it has the rate of convergence of the standard mixed finite-element method over the fine mesh. Alternatively, it can be seen as a postprocessed method in which the rate of convergence is increased by one unit with respect to the coarse mesh. The analysis takes into account the loss of regularity at initial time of the solution of the Navier-Stokes equations in absence of nonlocal compatibility conditions. Some numerical experiments are shown.
We consider the incompressible Navier–Stokes equations
in a bounded domain () with a smooth boundary subject to homogeneous Dirichlet boundary conditions on . In (1), is the velocity field, the pressure, the diffusion coefficient and a given force field.
In this paper we study the following two-grid mixed finite-element method for the spatial discretization of the above equations. First, for the solution of the fully nonlinear Navier-Stokes equations (1) corresponding to a given initial condition
the mixed finite-element approximation over a coarse mesh of diameter is computed. Then, for any time , the postprocessed approximation is obtained as the mixed finite-element approximation over a finer mesh () to the following steady Oseen-type problem:
In this paper we prove that, in terms of the fine mesh diameter , this two-grid technique is of optimal order in the sense that, for appropriate choices of the coarse mesh diameter , the method has the same rate of convergence of standard mixed finite element approximations in the fine mesh. On the other hand, for a suitable value of the discretization parameter , the rate of convergence of the postprocessed approximation in terms of increases by one unit the rate of convergence of the coarse standard approximation. The improvement in precision is achieved in both the norm for the velocity and the norm for the pressure in the case of linear, quadratic and cubic elements. For other than linear elements the rate of convergence in the norm of the velocity is also increased by one unit. We remark that time evolution is performed only at the coarse mesh whereas at the fine grid the time appears only as a parameter (see equation (1)), thus the name of static two-grid method.
Two-grid or two-level methods are a well established technique for nonlinear steady problems, see . In ,  several two-level methods are considered to approximate the steady Navier-Stokes equations. They require solving a nonlinear system over a coarse mesh and, depending on the algorithm chosen, one Stokes problem, one linear Oseen problem or one Newton step over the fine mesh. The corresponding algorithms obtain the optimal rate of convergence in the fine mesh for appropriate choices of the coarse mesh diameter .
In the case of nonlinear evolutionary equations, two-grid techniques have been proposed and studied in , , , . In these methods, as opposed to the method studied in the present paper, time evolution is also performed over the fine mesh. The advantage of the method studied in the present paper is that since the time integration is only carried out on the coarse mesh, computations on the fine grid can be done at selected target time levels where an improved approximation is desired, with the corresponding reduction of computing time, specially if these target time levels are sufficiently spaced in time. For this reason, although some of the two-grid methods that incorporate the evolution in time of the fine mesh approximation are more accurate, the method we present can still be more efficient in terms of computational effort for a given error level.
Two-grid techniques that integrate in time only on the coarse level have previously been developed in ,  (see also ) for spectral methods, and later extended to mixed finite-element formulations in , , . In all these works the two grid method is referred to as postprocessed Galerkin method, and, instead of (1), the approximation is found as an approximation to the following Stokes problem
This two-grid method will be termed standard postprocessed method, to differentiate it to that studied in the present paper, which will be termed new postprocessed method. Both, the standard and the new postprocessed methods, have the same rate of convergence. However, as already noted in  for nonlinear convection-diffusion problems, the new postprocessing technique produces more accurate approximations than the standard postprocessed method, for moderate to small values of the diffusion parameter . This will also be the case in the numerical experiments in the present paper for moderate values of the Reynolds number.
In the present paper we take into account the loss of regularity suffered by the solutions of the Navier-Stokes equations at the initial time in the absence of nonlocal compatibility conditions. Thus, for the analysis, we do not assume the solution to have more than second-order spatial derivatives bounded in up to initial time , since demanding further regularity requires the data to satisfy nonlocal compatibility conditions unlikely to be fulfilled in practical situations , . Due to the loss of regularity at , the best error bound that we can obtain is . For this reason we do not analyze higher than cubic finite elements. The same limit in the rate of convergence was found in  for standard mixed finite-element approximations and in ,  for two-grid schemes.
In practice, any method to numerically solve evolutionary equations needs of some time discretization procedure. For brevity reasons, we have preferred to present the method in a semidiscrete manner without reference to any particular time discretization. However, we emphasize that being static, the method we present can be applied exactly in the same form, with any time discretization. The analysis of fully discrete procedures can be developed along the same lines that appear in , .
The rest of the paper is as follows. In Section 2 we introduce some preliminaries and notation. In Section 3 we carry out the error analysis of the new method. Finally, some numerical experiments are shown in the last section.
2 Preliminaries and notations
We will assume that is a bounded domain in , of class , for . When dealing with linear elements ( below) may also be a convex polygonal or polyhedral domain. We consider the Hilbert spaces
endowed with the inner product of and , respectively. For integer and , we consider the standard spaces, , of functions with derivatives up to order in , and . We will denote by the norm in , and will represent the norm of its dual space. We consider also the quotient spaces with norm .
We recall the following Sobolev’s imbeddings : For , there exists a constant such that
For , (3) holds with .
The following inf-sup condition is satisfied (see ), there exists a constant such that
Let be the projection onto . We denote by the Stokes operator on :
We shall assume that is a strong solution up to time , so that
for some constants and . We shall also assume that there exists a constant such that
Finally, we shall assume that for some
so that, according to Theorems 2.4 and 2.5 in , there exist positive constants and such that the following bounds hold:
where and for some . Observe that for , we can take and . For simplicity, we will take these values of and .
Let , be a family of partitions of suitable domains , where is the maximum diameter of the elements , and are the mappings of the reference simplex onto .
Let , we consider the finite-element spaces
where denotes the space of polynomials of degree at most on .
and the so-called mini-element  when , where , and . Here, is spanned by the bubble functions , , defined by , if and 0 elsewhere, where denote the barycentric coordinates of . For these elements a uniform inf-sup condition is satisfied (see ), that is, there exists a constant independent of the mesh grid size such that
The approximate velocity belongs to the discrete divergence-free space
which is not a subspace of .
Let be the solution of a Stokes problem with right-hand side , we will denote by the so-called Stokes projection (see ) defined as the velocity component of the solution of the following problem: find such that
The following bound holds for :
where the constant depends on the constant in the inf-sup condition (9).
where for any .
3 The new postprocessed method
The postprocessing technique we propose is a two-level or two-grid method. In the first level, we choose a coarse mesh of size and compute the mixed finite-element approximation to defined by (14)-(15). In the second level, the discrete velocity and pressure are postprocessed by solving the following linear Oseen problem: find , , satisfying for all and
Equations (18)-(19) can also be solved over a higher order mixed finite-element space over the same grid. For simplicity in the exposition we will only consider the case in which we refine the mesh at the postprocessing step.
Let us observe that projecting equation (18) over the discretely-free space , and avoiding for simplicity the dependence on in the notation, we get that satisfies
We proceed to show that is coercive which implies that there exists a unique function satisfying (20). Let us also observe that once a unique is found, using the inf-sup condition (9) one easily obtains the existence and uniqueness of the pair satisfying (18)-(19).
Let be the bilinear form defined in (21). Then, there exists a constant such that for the following bound holds:
where if the dimension is , and if .
To prove the coercivity of we follow [25, p. 2042]. Let us first observe that for any
Let be the orthogonal projection of over , so that applying standard finite-element theory  and interpolation theory on Hilbert spaces (see e. g. [31, § II.2] we have , for . Taking into account that the velocity satisfies then
Following [25, p. 2042] we get
We introduce now a linearized problem that will be used in the proof of Theorem 1 where we state the rate of convergence of the new method. Let be the velocity in the solution of (1)-(2). We will denote by the solution of the following linearized problem
in the domain subject to homogeneous Dirichlet boundary conditions. Let us observe that since the divergence of is zero the bilinear form:
associated to this problem is continuous and coercive. Since the solution of (25) satisfies
by the Lax-Milgram theorem there exists a unique solution . Due to (4) there exists also a unique pressure .
We will assume in the sequel that both problem (25) and its dual problem satisfy the regularity assumption
In the following lemma we state the rate of convergence of the mixed finite-element approximation to the solution of (25) defined as follows: find such that
Let be the solution of (25) and let be its mixed finite-element approximation. Then, the following bounds hold for
Since applying (12) we conclude is bounded by the righ-hand side of (29). The bound (30) for the pressure is readily obtained by means of the auxiliary value . Subtracting (27) from (29) and applying the inf-sup condition (9) one easily gets
We are left with the task of proving the bound for the norm of the error in the velocity. We will argue by duality. Let us observe that
Let us fix and let us denote by the solution of the linearized dual problem
As stated before we assume that this problem satisfies the regularity assumption (26), so that
We will denote by the mixed finite-element approximations to . Reasoning exactly as before and applying (35) we obtain
Integrating by parts we reach
Then, to conclude, it only remains to bound which by (31) is equal to . Let us decompose
Then, integrating by parts in the last term
so that the proof is finished.
We now state some results that will be use to get the rate of convergence of the new postprocessed method. The proof of the following lemma can be found in [15, Lemma 4] for the case and in [10, Lemma 5.1] for .
The proof of the following lemma can be found in [10, p. 226].
We end this section with a theorem that states the rate of convergence of the new postprocessed method.
Let us consider the linearized problem (25) with right hand side . Then, the solution of (25) is the solution of (1)–(2). Let us denote by its mixed finite-element approximation, that is the solution of (27)–(28). This approximation satisfy the error bounds (29) and (30) for . Let us decompose and . To bound the first terms in these two decompositions we will apply (29) and (30). In the rest of the proof we deal with the other two terms.