Spectral discretization of time-dependent vorticity–velocity–pressure formulation of the Navier–Stokes equations

In this work, we propose a nonstationary Navier–Stokes problem equipped with an unusual boundary condition. The time discretization of such a problem is based on the backward Euler’s scheme. However, the variational formulation deduced from the nonstationary Navier–Stokes equations is discretized using the spectral method. We prove that the time semidiscrete problem and the full spectral discrete one admit at most one solution.


Introduction
We consider to be an open, bounded, and simply-connected domain of R d (d = 2, 3), and ∂ as its Lipschitz-continuous connected boundary. Let [0, T] be an interval in R where T is a positive real number. We denote x = (x, y) or x = (x, y, z) according to the dimension, and let n be the unit outward normal vector to on the boundary ∂ . We are interested in this paper in the following time-dependent Navier-Stokes system: where u and P are the unknowns denoting the velocity and pressure of the fluid, f represents the density of the body forces, and ν is the viscosity that we suppose to be a positive constant. The operator τ defines the boundary value of curlu in dimension d = 2 and the boundary tangential components of curlu in dimension d = 3. This problem with nonstandard boundary conditions was studied for the first time in the pioneering paper [1] where the domain is assumed to be convex in both dimension d = 2 and d = 3. The basic idea in [2,3] consists in introducing the vorticity ω = curlu as a new unknown (see also [4][5][6]) and the fact that the convection term can be written as: Consequently, problem (1) is fully equivalent to the following system: The dynamical pressure p is defined by the formula We assume that the initial velocity and vorticity satisfy the following conditions: div u 0 = 0 in and ω(x, 0) = ω 0 = curlu 0 in .
This type of problem has been handled in the stationary case by the finite element method [3,7] and by the spectral method [8,9]. However, the nonstationary case has been considered for a posteriori analysis of the finite discretization [10,11] and also for a spectral discretization of the linear Stokes problem [12]. Relying on the equivalent variational formulation of problem (2), we propose a semidiscrete scheme based on the backward Euler method. In dimension d = 2, we show that the time semidiscrete problem admits a solution with no restriction on the regularity of the domain, while in dimension d = 3, we assume that the domain has either a C 1,1 boundary or is a polyhedron with no reentrant corners. Note, however, that this result of existence is established only for a sufficiently large viscosity in dimension d = 3. For such a problem, we also study the uniqueness of the solution. We prove the existence of a discrete solution of the full spectral discrete problem. The paper is organized as follows: • In Sect. 2, we present the variational formulation and the analysis of the model. • Sect. 3 is devoted to the time semidiscrete problem and to the proof of existence and uniqueness results.
• In Sect. 4 we prove the well posedness of the full spectral discrete problem.

The variational formulation
To write the variational formulation of problem (2), we suggest the following Sobolev spaces: each of which is a Banach space equipped with the following norm and seminorm: We denote by (·, ·) the L 2 ( ) scalar product. We note that L 2 0 ( ) is the space of functions in L 2 ( ) which have a null integral on , and D( ) is the space of indefinitely differentiable functions with a compact support in . We consider the domain H(div, ) of the div operator, We remind (see [13, Chap. I, Sect. 2]) that the normal trace operator is defined from H(div, ) into H -1 2 (∂ ) such that for any scalar function χ in H(div, ), where ·, · is the duality product between H -1 2 (∂ ) and H 1 2 (∂ ). This allows us to introduce the kernel of the normal trace operator in H(div, ), H 0 (div, ) = {ϕ ∈ H(div, ); ϕ · n = 0 on ∂ }.
In order to handle some problems depending on time, we consider the following spaces. Let Z be a separable Banach space. We introduce C m (0, T; Z), the set of time C m class functions with values in Z; C m (0, T; Z) is a Banach space provided with the norm where ∂ i t ϕ is the partial derivative of order i in time of the function ϕ. Consider also the spaces Finally, we define also L(Z) to be the Banach space of the linear and continuous functions from Z to R provided with the norm We assume that the data f belongs to the space L 2 (0, T; H 0 (div, ) ), where H 0 (div, ) is the dual space of H 0 (div, ) (see [14] for more details about the space H 0 (div, ) ). Consider the following variational formulation: Find (ω, u, p) ∈ L 2 (0, T; H 0 (curl, )) × L 2 (0, T; H 0 (div, )) × L 2 (0, T; L 2 0 ( )) such that where ≺ ·, · is the duality product between H 0 (div, ) and H 0 (div, ); a(·, ·; ·), b(·, ·), and c(·, ·; ·) are bilinear forms defined as follows: On the other hand, the trilinear form N(·, ·; ·) is defined by By the same arguments used in our work about the stationary Stokes problem (see [9, Proposition 2.1]), we deduce the following

and only if it is a solution of problem (4) in the sense of distributions.
To prove that problem (4) has a solution, we need to define which is the kernel of the bilinear form b(·, ·) and is the space of divergence-free functions in H 0 (div, ), since H 0 (div, ) is included in L 2 0 ( ). We also consider the kernel of the bilinear form c(·, ·; ·). We note that V and U are Hilbert spaces, since the bilinear forms b(·, ·) and c(·, ·; ·) are continuous. So we observe that if (ω, u, p) is a solution of problem (4), then (ω, u) is a solution of the following reduced problem: The main difficulty consists in proving the existence of a solution of problem (6). Moreover, for the three-dimensional case, in order to give a sense to the nonlinear term N(ω, u; v), we need to assume that the spaces H 0 (div, ) ∩ H(curl, ) and H(div, ) ∩ H 0 (curl, ) are compactly embedded in H 1 ( ) see ( [15], Theorem 2.17). Nonetheless, the proof is much simpler in dimension d = 2.

Assumption 1
In dimension d = 3, we assume that the domain has a C 1,1 boundary or is convex.

I, Corollary 2.4]):
There exists a constant β * > 0 such that When Assumption 1 is satisfied, the arguments for the proof of the existence of a solution of problems (6) and (4) [8]. We prefer to omit this proof since it is beyond the aim of this paper.

The time semidiscrete problem
In this section we study the discretization in time of problem (4) using the implicit Euler's method. We consider a partition of the interval . . , h n ) stands for the step of the partition where h n = t nt n-1 , and |h| = max 1≤n≤M h n .
For any data f ∈ L 2 (0, T; (H 0 (div, )) ), u 0 ∈ H 0 (div, ), and ω 0 ∈ H 0 (curl, ) satisfying condition (3), we consider the following problem: and for all 1 ≤ n ≤ M, where f n = f(·, t n ), Then, if (ω n , u n , p n ) is a solution of problem (8)-(9), (ω n , u n ) belongs to U and is a solution of the following reduced problem: The main difficulty now consists in showing that problem (10) admits a solution. We recall here the two properties (positivity and inf-sup condition) related to the bilinear form A(·, ·; ·) proved in [12, Lemma 1]: There exists a constant β > 0 such that A ω n , u n ; v > 0 (11) and Let now study the properties of non linear term N(·, ·; ·).
• In dimension 3, if Assumption 1 holds, using the Sobolev embedding theorem, the space U is compactly included in L r ( ) 3 × L s ( ) 3 for r = 4 and s = 4.
2) Since u ∈ V and equality (15) holds, we easily deduce the antisymmetry property (14). Equality (15) is proved by using Green formula: Proposition 2 Assume that the data f belongs to L 2 (0, t; H 0 (div, ) ) and that the initial vorticity-velocity (ω 0 , u 0 ) belongs to H 0 (curl, ) × H 0 (div, ) and satisfies (3). Knowing u n-1 at each time step n, problem (10) has a solution (ω n , u n ) in U for dimension d = 2. This solution satisfies for n ≥ 1: where c is a positive constant independent of n.
Proof We consider the function defined on U into its dual space by We derive from the continuity of the trilinear form N in (13) that the function is continuous on the space U. Having the following inequality proved in [12]: and the fact that N(ω, u; curlω) = 0 in dimension 2, we obtain where c 0 is the smallest constant such that Then if we consider the sphere with radius the duality product (ω n , u n ), (θ , v) is nonnegative on this sphere. The space U is a separable Hilbert space since it is included in L 2 ( ) × L 2 ( ) 2 . Then there exists an increasing sequence (U k ) k of finite-dimensional subspaces U k of U such that k (U k ) is dense in U. Using Brouwer's fixed point theorem [13, Chap. IV, Corollary 1.1], there exists (ω n k , u n k ) ∈ U k such that Since the sequence (ω n k , u n k ) is bounded by r, there exists a subsequence, still denoted by (ω n k , u n k ), that converges weakly to (ω n , u n ). Now, using the fact that the two operators ∀(θ , v) ∈ U: (ω, u) → ω × v and (ω, u) → θ × u are compact on U, we conclude in passing to the limit in equation (18) that Using the density of k (U k ) in U yields ∀(θ , v) ∈ U, ω n , u n , (θ , v) = 0.
Then we conclude that (ω n , u n ) is a solution of problem (10) for any w = θ + curlv ∈ V. Moreover, since the solution (ω n , u n ) is bounded by r, we obtain (16) by the same proof as for [12, Corollary 1].
In dimension d = 3, the existence of the solution of problem (10) is obtained for a large enough viscosity ν with respect to the norm of the operator L.

Proposition 3
Assume that the data f belongs to L 2 (0, T; H 0 (div, ) ) and that the initial vorticity-velocity (ω 0 , u 0 ) belongs to H 0 (curl, ) × H 0 (div, ) and satisfies (3). If Assumption 1 holds and there exists a constant c * such that c * ν -2 L L(H 0 (div, )) < 1, then, knowing u n-1 at each time step n, problem (10) has a solution (ω n , u n ) ∈ U. This solution satisfies n j=1 where c is a positive constant independent of n.
Proof We consider an iterative sequence ((ω n k , u n k )) k such that (ω n 0 , u n 0 ) = (0, 0) and (ω n k , u n k ) is the solution in U of the following problem: From the properties of the bilinear form A(·, ·; ·) in (11)- (12) and the continuity of the trilinear form N(·, ·; ·) in (13), problem (21) has a unique solution; see [9]. Let where M is the continuity constant of N(·, ·; ·). By induction on k and by an appropriate choice of c * in (19), the sequence ((ω n k , u n k )) k is bounded by r: We have for any k ≥ 2, Thanks to (17), (13), and (23), together with (22), we obtain that The sequence ((ω n k , u n k )) k is a Cauchy sequence in U, so it converges to a pair (ω n , u n ). By passing to the limit in (21), it is readily checked that (ω n , u n ) is a solution of problem (10). Moreover, the solution (ω n , u n ) is bounded by r, then we obtain (20) by the same proof as for [12,Corollary 1].
By using (17) and taking v = u n + curlω n in the previous line, we have Using the antisymmetry property of the trilinear form N(·, ·; ·) (14) gives N(ω n 1 , u n ; u n ) = 0. Applying the continuity of N(·, ·; ·) in (13) together with (16) and (20) gives Thus, if condition (19) holds, we obtain that ω n and u n are equal to zero. Let (ω n , u n ) be a solution of problem (8)- (10). We consider for any v ∈ H 0 (div, ), a linear and continuous functional which vanishes on the space V. Then, according to the inf-sup condition (7), there exists a unique p n in L 2 0 ( ) such that ∀v ∈ H 0 (div, ), ϒ n (v) = -div v, p n and p n L 2 ( ) ≤

The full discrete problem
Hereinafter, we suppose that is a square [-1, 1] 2 in dimension d = 2 or a cube [-1, 1] 3 in dimension d = 3. To define the discrete spaces, we follow the same idea proposed by Nédélec in the case of the finite element method (see [20,Sect. 2]). If d = 2 (resp. d = 3), let P nm ( ) (resp. P nms ( )) be the space of the restriction on of the polynomials of degree n in the x direction and m in the y direction (resp. and s in the z direction). Let N ≥ 2 be an integer. We consider the spaces: which approximates the velocity in H 0 (div, ) and which approximates the vorticity in H 0 (curl, ). This space is rather different conforming to the dimension (see Remark 1). Lastly, for the approximation of pressure in L 2 0 ( ), we propose the space Having ξ 0 = -1 and ξ N = 1, we consider the N -1 nodes ξ j , 1 ≤ j ≤ N -1, roots of the polynomial L N , where L N is the Legendre polynomial of degree N and the N +1 weights ρ j , 0 ≤ j ≤ N , of the Gauss-Lobatto quadrature formula. Let P n (-1, 1) the space of restrictions to ] -1, 1[ of polynomials with degree ≤ n, then ∀φ ∈ P 2N-1 (-1, 1), We remind the important following property, see [21] for its proof: Based on formula (25), we define the following discrete scalar product on P N ( ): For a continuous functions ϕ and ψ on¯ , Finally, we consider I N to be the Lagrange interpolating operator at the nodes (ξ i , ξ j ), 0 ≤ i, j ≤ N , in dimension d = 2, and at the nodes (ξ i , ξ j , ξ k ), 0 ≤ i, j, k ≤ N , in dimension d = 3, with values in the space P N ( ). We assume that the data f is continuous on × [0, T]. We construct the discrete problem from problem (8)-(9) by means of the Galerkin method combined with numerical integration.
If u 0 N = I N (u 0 ), knowing u n-1 , we find (ω n N , u n N , where the bilinear forms A N (·, ·; ·), b N (·, ·), and c N (·, ·; ·) are defined by By formula (25) combined with Cauchy-Schwarz inequality, we prove that the bilinear forms A N (·, ·; ·), b N (·, ·), and c N (·, ·; ·) are respectively continuous on ( N is linear and continuous on D N . Moreover, as a consequence of the exactness property (24), the bilinear forms b(·, ·) and b N (·, ·) coincide on D N × M N . While the trilinear form N N (·, ·; ·) is defined as follows: To prove that problem (26) has a solution, we consider the kernel We remark that V N is the space of divergence-free polynomials in D N (if q N = div v N in (27)). Let also the kernel We note that the discrete kernel U N is not contained in U in the general case; see [9,Corollary 3.2]. We consider herein the following reduced discrete problem: If u 0 N = I N (u 0 ) and knowing u n-1 , find (ω n N , u n N ) ∈ U N such that for 1 ≤ n ≤ M, The existence of a solution of problem (28) is proved by the same arguments as for the continuous reduced problem (10), for instance.

Proposition 4
For any data f continuous on¯ × [0, T] and knowing u n-1 N at each time step n, problem (28) admits a solution (ω n N , u n N ) in U N . Moreover, this solution satisfies, for where c is a positive constant independent of N and n.
Proof We introduce the mapping N defined from U N into its dual space by We equip U N with the following norm: The mapping N is continuous since the space U N has a finite dimension. Next, noting by the same arguments as for Lemma 1 from the property of antisymmetry that N N (ω n N , u n N ; u n N ) = 0, we have Then using (25) with the definition of U N gives For any u n N in V N , it follows from [9, Lemmas 3.4 and 3.5] that there exists τ n N in C N such that u n N = curlτ n N and Using once more (25) and inserting τ n N in the definition of the space U N , we have Then we remark that N (ω n N , u n N ), (ω n N , u n N ) is nonnegative on the sphere of U N of radius r N . So by Brouwer's fixed point theorem (see [13, Chap. IV, Corollary 1.1]), we conclude that problem (28) has a solution (ω n N , u n N ) in U N . Moreover, the solution (ω n N , u n N ) is bounded by r N , so we obtain (29) by the same proof as that of [12,Proposition 5)].
The following inf-sup condition is proved in [9, Lemma 3.9]: There exists a positive constant γ , independent of N , such that the form b N (·, ·) satisfies the inf-sup condition: The full existence result is deduced from Proposition 4 and the above inf-sup condition. The proof follows the same arguments as for Theorem 1. Remark 2 Note that the previous existence result still holds when N N (·, ·; ·) is replaced by N(·, ·; ·) in problem (26). This means in practice that a more precise quadrature formula, exact on P 3N-1 ( ), is used to evaluate the integrals that appear in the treatment of the nonlinear term.

Conclusion
This work concerns the numerical analysis of the implicit Euler scheme in time and the spectral discretization in space of the nonstationary vorticity-velocity-pressure formulation of the Navier-Stokes equations. We prove using Brouwer's fixed point theorem that the new discrete formulation has at most one solution. The study of the error, the algorithm solution, and the numerical implementation of these results will be the subject of our forthcoming work.