A New Proof of the Existence of Suitable Weak Solutions and Other Remarks for the Navier-Stokes Equations

Show more

1. Introduction

The main objective of this paper is to provide a new proof of the existence of suitable weak solutions to the Navier-Stokes equations. Specifically, we show that the semi-discrete and the completely discrete semi-implicit Euler schemes lead to families of approximate solutions that converge to a weak solution that is suitable in the sense of Caffarelli, Kohn and Nirenberg [2] .

We will be concerned with the 3D Navier-Stokes equations completed with initial and Dirichlet boundary conditions in bounded domains $\Omega \times \left(0,T\right)$ (as usual, Ω is the spatial domain, a regular, bounded and connected open set in ${\mathbb{R}}^{3}$ “filled” by the fluid particles; $\left(0,T\right)$ is the time observation interval).

The key concept of suitable weak solution was introduced in [2] . In few words, this is a weak solution satisfying a local energy inequality. Due to its definition, it is expected that suitable solutions are regular (and unique). However, up to now, this is unknown. The best we can prove is that the set of singular points of a suitable solution is small, in the sense that has Hausdorff dimension ≤1. As shown below, in order to improve this result, we would need an estimate that we do not have at hand at present. A similar analysis can be performed in the context of the Boussinesq system; see [4] .

The estimate in [2] of the Hausdorff dimension relies on some technical results asserting that adequate criteria, applied to suitable solutions in a given space-time region, imply the regularity of the points in a subregion. During the last years, several authors have tried to improve or weaken these criteria and some achievements have been obtained:

・ In Seregin [5] , a family of sufficient conditions that contains the Caffarelli-Kohn-Nirenberg condition as a particular case is introduced. Its formulation is given in terms of functionals invariant with respect to scale transformations.

・ In Vasseur [6] , an interesting criterion appears: if we normalize the solution and the sum of the associated kinetic and viscous energies and the L^{p} norm of the pressure is small enough, we get regularity. The proof of this assertion is inspired by a method of De Giorgi designed to prove the regularity of elliptic equations, see [7] .

・ In Wolf [8] , the author provides a notion of local pressure. It permits to estimate the integrals involving the pressure in terms of the velocity and, again, deduce regularity. The method is interesting and can be adapted to get partial regularity results for other systems, such as the equations of quasi-Newtonian or Boussinesq (heat conducting) fluids.

・ Finally, in Choe, Wolf and Yang [9] , an improved version of the Caffarelli-Kohn-Nirenberg criterion is furnished, using ideas from [5] .

Note that, in a recent paper, Buckmaster and Vicol [10] have proved that, for a very weak class of distributional solutions in spatially periodic domains, non-uniqueness occurs.

The techniques employed in this paper can be applied to many other approximation schemes that lead to energy inequalities, as those in [11] [12] [13] [14] . More precisely, we first use the well known energy estimates, together with appropriate interpolation results and recall that the approximate solutions converge to a weak solution $\left(u,p\right)$ . Then, we analyze the role of the pressure p; this reduces in fact to a detailed study of the behavior of the time derivate of the velocity field. This way, we are able to take (lower) limits in the local energy identities satisfied by the approximate solutions and deduce that $\left(u,p\right)$ is suitable.

Our results can be compared to other previous proofs of existence: the one in the Appendix in [2] (based on the construction of a family of time delayed linear approximations), the main result in Da Veiga [15] (relying on regularization with vanishing fourth-order terms), the main result in Guermond [3] (where Faedo-Galerkin techniques are employed) and, also, the results by Berselli and Spirito [16] [17] , where the Voigt approximations and the artificial compressibility method are shown to converge.

We think that our results can be useful at least from two points of view. First, a new (relatively simple) constructive argument is used to prove the existence of suitable solutions. Then, there is some practical interest: by inspection of the behavior of the computed approximations in a prescribed region, we may try to deduce if the related points are regular. In other words, checking whether or not the Caffarelli-Kohn-Nirenberg criteria are satisfied on the computed numerical solutions can serve to identify or discard singular points. Based on this idea, we will present in a forthcoming paper several numerical experiments for which interesting conclusions can be obtained.

The plan of the paper is the following:

・ In Section 2, we review the main results in the papers [2] and [18] . In particular, we explain why suitable solutions are relevant in the context of the regularity problem.

・ In Section 3, we recall the Euler approximation schemes and we establish the convergence to a suitable solution of the Navier-Stokes equations.

・ Finally, Section 4 is devoted to some additional comments and open questions.

2. Background: The Basic Results by Caffarelli, Kohn and Nirenberg

In the sequel, we denote by
$\left|\text{\hspace{0.05em}}\cdot \text{\hspace{0.05em}}\right|$ and
$\left(\cdot ,\cdot \right)$ the usual L^{2} norm and scalar product, respectively. The symbol C will be used to denote a generic positive constant.

2.1. The Main Properties of Suitable Solutions

In this section, we will recall the main contributions of Caffarelli, Kohn and Nirenberg, see [2] . In this reference, the best results known to date in relation to the regularity of the Navier-Stokes equations are established.

Let $\Omega \subset {\mathbb{R}}^{3}$ be a nonempty, regular, bounded and connected open set and assume that $T>0$ . Let us set $Q:=\Omega \times \left(0,T\right)$ and $\Sigma :=\partial \Omega \times \left(0,T\right)$ . We will consider local and global solutions to the Navier-Stokes equations in three dimensions

$\{\begin{array}{l}{u}_{t}+u\cdot \nabla u-\Delta u+\nabla p=f\left(x,t\right),\\ \nabla \cdot u=0,\end{array}$ (1)

where $f=\left({f}^{1},{f}^{2},{f}^{3}\right)$ verifies $f\in {L}^{q}{\left(\Omega \times \left(0,T\right)\right)}^{3}$ with $q\ge 2$ and $\nabla \cdot f=0$ .

At a local level, we will consider solutions in sets of the form $G\times \left(a,b\right)$ , where $G\subset \Omega $ is open and connected and $\left(a,b\right)\subset \left(0,T\right)$ :

Definition 2.1 Let the open set $D:=G\times \left(a,b\right)$ be given. It will be said that the couple $\left(u,p\right)$ is a weak solution to the Navier-Stokes equations, Equation (1) in D if the following holds:

・ $u\in {L}^{2}\left(0,T;{H}^{1}{\left(G\right)}^{3}\right)\cap {L}^{\infty}\left(0,T;{L}^{2}{\left(G\right)}^{3}\right)$ and $p\in {L}^{5/3}\left(G\times \left(a,b\right)\right)$ .

・ u and p satisfy the Navier-Stokes equations in Equation (1) in the distributional sense in D.

On the other hand, for the definition of a global solution, it will be convenient to use the spaces H and V, with

$H:=\left\{v\in {L}^{2}{\left(\Omega \right)}^{3}:\nabla \cdot v=0\text{\hspace{0.17em}}\text{in}\text{\hspace{0.17em}}\Omega ,v\cdot n=0\text{\hspace{0.17em}}\text{on}\text{\hspace{0.17em}}\partial \Omega \right\},$

$V:=\left\{v\in {H}_{0}^{1}{\left(\Omega \right)}^{3}:\nabla \cdot v=0\text{\hspace{0.17em}}\text{in}\text{\hspace{0.17em}}\Omega \right\}.$

Let us assume that ${u}_{0}\in H$ and let us consider the initial-boundary value problem

$\{\begin{array}{l}{u}_{t}+u\cdot \nabla u-\Delta u+\nabla p=f\left(x,t\right),\left(x,t\right)\in Q,\\ \nabla \cdot u=0,\left(x,t\right)\in Q,\\ u\left(x,t\right)=0,\left(x,t\right)\in \Sigma ,\\ u\left(x,0\right)={u}_{0}\left(x\right),x\in \Omega .\end{array}$ (2)

Definition 2.2 It will be said that the couple $\left(u,p\right)$ is a (global) weak solution to Equation (2) if the following holds:

・ $u\in {L}^{2}\left(0,T;V\right)\cap {L}^{\infty}\left(0,T;H\right)$ and $p\in {L}^{5/3}\left(Q\right)$ .

・ u and p satisfy the Navier-Stokes equations in Equation (2) in the distributional sense in Q

・ $u\left(x,0\right)={u}_{0}\left(x\right)$ a.e. in Ω.

It is known that any couple $\left(u,p\right)$ satisfying the previous first and second points also verifies

$u\in {L}^{10/3}{\left(Q\right)}^{3}\cap {C}_{w}^{0}\left(\left[0,T\right];H\right),\text{\hspace{0.17em}}{u}_{t}\in {L}^{4/3}\left(0,T;{V}^{\prime}\right).$

In particular, u can be viewed as a well-defined H-valued function and the third assertion in Definition 2.2 has a sense as an equality in H.

In order to understand the role and relevance of the terms in the estimates that follow, it is convenient to associate a dimension to each variable in Equation (2). Note that, if the pair $\left(u,p\right)$ is a weak solution to the Navier-Stokes equations in $D=G\times \left(a,b\right)$ , then for each $\lambda >0$ the functions

${u}_{\lambda}\left(x,t\right):=\lambda u\left(\lambda x,{\lambda}^{2}t\right)$ and ${p}_{\lambda}\left(x,t\right):={\lambda}^{2}p\left(\lambda x,{\lambda}^{2}t\right)$

solve a similar problem in ${D}_{\lambda}:={\lambda}^{-1}G\times \left({\lambda}^{-2}a,{\lambda}^{-2}b\right)$ , with force ${f}_{\lambda}:={\lambda}^{3}f\left(\lambda x,{\lambda}^{2}t\right)$ . Thus, for any integer k, we say that a variable or a linear differential operator is of dimension k if it is non-dimensionalized when it is multiplied by ${\lambda}^{-k}$ , where λ is a characteristic length. We can affirm that

・ ${x}_{i}$ has dimension 1 and t is of dimension 2,

・ ${u}^{i}$ has dimension −1 and p has dimension −2,

・ f has dimension −3,

・ ${\partial}_{i}$ has dimension −1 and ${\partial}_{t}$ has dimension 2,

so that all the terms of the motion equation in Equation (2) have dimension −3.

The analysis of the existence of a weak solution to Equation (2) can be found for instance in [19] and [20] . Now, we will speak of the regularity problem, that is, the possible regularity properties of the weak solution. To this purpose, let us consider the following definition:

Definition 2.3 Let $\left(u,p\right)$ be a weak solution to Equation (1) in $D=G\times \left(a,b\right)$ and let $\left({x}_{0},{t}_{0}\right)\in D$ be given. It will be said that $\left({x}_{0},{t}_{0}\right)$ is a singular point if u is not ${L}^{\infty}$ in any neighborhood of $\left({x}_{0},{t}_{0}\right)$ , that is, there are no r and C such that $\left|u\left(x,t\right)\right|\le C$ for $\left(x,t\right)$ a.e. in $B\left(\left({x}_{0},{t}_{0}\right);r\right)$ . The remaining points, those where u is locally bounded, will be called regular points.

According to a result by Serrin [21] , it is known that, if $\left(u,p\right)$ is a weak solution to Equation (2) and $\left({x}_{0},{t}_{0}\right)\in D$ is a regular point for u, then u coincides a.e. with a ${C}^{\infty}$ function in a neighborhood of $\left({x}_{0},{t}_{0}\right)$ . This gives an idea of how interesting can be to get a description of the set S of singular points.

In fact, Serrin proved that, in order to have u of class ${C}^{\infty}$ near $\left({x}_{0},{t}_{0}\right)$ , one just needs an estimate of the kind ${L}^{r}$ in time and ${L}^{s}$ in space, with sufficiently large r and s. Note that, in [22] , it is shown that a weaker condition is sufficient for ${C}^{\infty}$ regularity. Note also that, in accordance with the results in [23] , if one component of the velocity field is essentially bounded in a region, there is no singular point in a subregion.

The first papers devoted to describe S are due to Scheffer [1] [24] [25] . There, some estimates of the size of the set were given in terms of appropriate Hausdorff measures. Actually, the main result in [1] is the following:

Theorem 2.4 Assume that $f\equiv 0$ . There exists a weak solution to Equation (2) whose associated singular set S satisfies:

${\mathcal{H}}^{5/3}\left(S\right)<+\infty $ and ${\mathcal{H}}^{1}\left(S\cap \left(\Omega \times \left\{t\right\}\right)\right)<+\infty $ uniformly in t.

Here, ${\mathcal{H}}^{k}$ denotes the usual Hausdorff k-dimensional measure in ${\mathbb{R}}^{4}$ .

This result was improved by Caffarelli, Kohn and Nirenberg in [2] in several directions. There, the authors used a particular class of weak solution, denoted suitable weak solution or simply suitable solution, according to the following definition:

Definition 2.5 Let $D=G\times \left(a,b\right)$ be a cylinder in ${\mathbb{R}}^{3}\times \mathbb{R}$ . It is said that $\left(u,p\right)$ is a suitable weak solution to the Navier-Stokes equations in D if it satisfies points 1 and 2 of Definition 2.1 and, furthermore, the following generalized energy inequality property: for any $\varphi \in {C}_{0}^{\infty}\left(D\right)$ with $\varphi \ge 0$ ,

$2{\displaystyle {\iint}_{D}{\left|\nabla u\right|}^{2}\varphi}\le {\displaystyle {\iint}_{D}\left({\left|u\right|}^{2}\left({\varphi}_{t}+\Delta \varphi \right)+\left({\left|u\right|}^{2}+2p\right)u\cdot \nabla \varphi +2\left(u\cdot f\right)\varphi \right)}\text{\hspace{0.05em}}.$

Then, the authors of [2] introduced the so called “parabolic” Hausdorff measure ${\mathcal{P}}^{1}$ , as follows:

・ First, for any small $\delta >0$ and any $X\subset {\mathbb{R}}^{4}$ , they set

${\mathcal{P}}_{\delta}^{1}\left(X\right):=\mathrm{inf}\left\{\underset{i\ge 1}{{\displaystyle \sum}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{r}_{i}:X\subset {\displaystyle \underset{i\ge 1}{\cup}{Q}_{{r}_{i}}},{r}_{i}<\delta \right\}.$

Here, each ${Q}_{{r}_{i}}$ is a parabolic cylinder, that is, a set of the form

${Q}_{r}:=\left\{\left(x,t\right)\in {\mathbb{R}}^{4}:\left|x-\stackrel{\xaf}{x}\right|\le r,\left|t-\stackrel{\xaf}{t}\right|\ge {r}^{2}\right\}$

for some $\left(\stackrel{\xaf}{x},\stackrel{\xaf}{t}\right)\in Q$ .

・ Then, for any $X\subset {\mathbb{R}}^{4}$ , they set

${\mathcal{P}}^{1}\left(X\right):=\underset{\delta \to {0}^{+}}{\mathrm{lim}}{\mathcal{P}}_{\delta}^{1}(X)$

With the help of ${\mathcal{P}}^{1}$ , a local partial regularity result can be established for any suitable solution:

Theorem 2.6 Let $\left(u,p\right)$ be a suitable solution to Equation (1) in D. Then the associated singular set satisfies ${\mathcal{P}}^{1}\left(S\right)=0$ .

This result improves Theorem 2.4 in several aspects: first, it has local character; then, it allows a rather general force term f; finally, it gives a better estimate of the Hausdorff dimension of S, since one has

${\mathcal{H}}^{1}\left(X\right)\le C{\mathcal{P}}^{1}\left(X\right),\text{\hspace{0.17em}}\forall X\subset {\mathbb{R}}^{4}$

for some $C>0$ .

In the sequel, for any $\left(x,t\right)$ and $r>0$ , we will denote by ${Q}_{r}\left(x,t\right)$ the following parabolic cylinder:

${Q}_{r}\left(x,t\right):=\left\{\left(y,\tau \right):\left|y-x\right|<r,t-{r}^{2}<\tau <t\right\}.$

For the proof of Theorem 2.6, we need two results. The first one is the following:

Proposition 2.7 Suppose that $\left(u,p\right)$ is a suitable weak solution to Equation (1) in ${Q}_{1}:={Q}_{1}\left(0,0\right)$ and $f\in {L}^{q}{\left({Q}_{1}\right)}^{3}$ with $q>5/2$ . There exist ${\u03f5}_{1}>0$ , ${C}_{1}>0$ and ${\u03f5}_{2}={\u03f5}_{2}\left(q\right)>0$ such that, if

${\iint}_{{Q}_{1}}\left({\left|u\right|}^{3}+\left|u\right|\left|p\right|\right)}+{\displaystyle {\int}_{-1}^{0}{\left({\displaystyle {\int}_{\left|x\right|<1}\left|p\right|\text{d}x}\right)}^{5/4}\text{d}t}\le {\u03f5}_{1$ (3a)

and

${\iint}_{{Q}_{1}}{\left|f\right|}^{q}}\le {\u03f5}_{2$ (3b)

then

$\left|u\left(x,t\right)\right|\le {C}_{1}\text{\hspace{0.17em}}\text{\hspace{0.17em}}a.e.\text{\hspace{0.17em}}\text{in}\text{\hspace{0.17em}}{Q}_{1/2}:={Q}_{1/2}\left(0,0\right).$ (3c)

In particular, $\left(0,0\right)$ is a regular point.

Proposition 2.7 shows that the sizes of the data have an influence on the regularity of suitable solutions. Now, if we introduce

$M\left(r\right):=\frac{1}{{r}^{2}}{\displaystyle {\iint}_{{Q}_{r}}\left({\left|u\right|}^{3}+\left|u\right|\left|p\right|\right)}+{r}^{-13/4}{\displaystyle {\int}_{t-{r}^{2}}^{t}{\left({\displaystyle {\int}_{\left|y-x\right|<r}\left|p\right|\text{d}y}\right)}^{5/4}\text{d}\tau}$ (4)

and

${F}_{q}\left(r\right):={r}^{3q-5}{\displaystyle {\iint}_{{Q}_{r}}{\left|f\right|}^{q}},$ (5)

taking into account the dimensions of these quantities, we can easily deduce the following:

Corollary 2.8 Suppose that $\left(u,p\right)$ is a suitable solution to Equation (1) in the cylinder ${Q}_{r}\left(x,t\right)$ and $f\in {L}^{q}{\left({Q}_{r}\left(x,t\right)\right)}^{3}$ , with $q>5/2$ . Then, if $M\left(r\right)\le {\u03f5}_{1}$ and ${F}_{q}\left(r\right)\le {\u03f5}_{2}$ , one has

$\left|u\right|\le {C}_{1}{r}^{-1}\text{\hspace{0.17em}}\text{\hspace{0.17em}}a.e.\text{\hspace{0.17em}}\text{\hspace{0.05em}}\text{in}\text{\hspace{0.17em}}{Q}_{r/2}\left(x,t\right)$

and, consequently, every point in ${Q}_{r/2}\left(x,t\right)$ is regular.

Let us now set

${Q}_{r}^{*}\left(x,t\right):=\left\{\left(y,\tau \right):\left|y-x\right|<r,t-\frac{7}{8}{r}^{2}<\tau <t+\frac{1}{8}{r}^{2}\right\}.$

The second fundamental result used in the proof of Theorem 2.6 is the following:

Proposition 2.9 Let $\left(u,p\right)$ be a suitable solution to Equation (1) in a neighborhood of $\left(x,t\right)$ . There exists ${\u03f5}_{3}>0$ such that, if

$\underset{r\to 0}{\mathrm{lim}\mathrm{sup}}\frac{1}{r}{\displaystyle {\iint}_{{Q}_{r}^{*}\left(x,t\right)}{\left|\nabla u\right|}^{2}}\le {\u03f5}_{3},$

then $\left(x,t\right)$ is a regular point.

For the proofs of Propositions 2.7 and 2.9, Caffarelli, Kohn and Nirenberg used the generalized energy inequality with well chosen test functions ϕ; a simpler proof is given in [18] . Then, Theorem 2.6 is deduced from these results by contradiction using a covering lemma and the usual energy estimates.

2.2. Sketch of the Proofs of Theorems 2.4 and 2.6

Theorem 2.6 is a consequence of Proposition 2.9. The argument is explained below.

Consider first the proof of the fact that S has Hausdorff dimension less than or equal to 5/3, that is, Theorem 2.4. Using Corollary 2.8 and a covering lemma, we can easily see that, for each $\delta >0$ , S can be covered by a family of parabolic cylinders $\left\{{Q}_{{r}_{i}}^{*}\left({x}_{i},{t}_{i}\right)\right\}$ such that ${r}_{i}<\delta $ , the ${Q}_{{r}_{i}/5}^{*}\left({x}_{i},{t}_{i}\right)$ are mutually disjoint and

$\frac{1}{{r}^{2}}{\displaystyle {\iint}_{{Q}_{{r}_{i}/5}^{*}\left({x}_{i},{t}_{i}\right)}\left({\left|u\right|}^{3}+\left|u\right|\left|p\right|\right)}+{r}^{-13/4}{\displaystyle {\int}_{{t}_{i}-7{r}_{i}^{2}/8}^{{t}_{i}+{r}_{i}^{2}/8}{\left({\displaystyle {\int}_{\left|x-{x}_{i}\right|<{r}_{i}}\left|p\right|\text{d}y}\right)}^{5/4}\text{d}\tau}>C{\u03f5}_{1}$ (6)

for all i. Using Hölder’s inequality, we deduce that

${r}_{i}^{-5/3}{\displaystyle {\iint}_{{Q}_{{r}_{i}/5}^{*}\left({x}_{i},{t}_{i}\right)}\left({\left|u\right|}^{10/3}+{\left|p\right|}^{5/3}\right)}\ge C(\u03f51)$

and, therefore,

${{\displaystyle \sum}}^{\text{}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{r}_{i}^{5/3}\le {\displaystyle {\iint}_{\cup {Q}_{{r}_{i}/5}^{*}\left({x}_{i},{t}_{i}\right)}\left({\left|u\right|}^{10/3}+{\left|p\right|}^{5/3}\right)}\le C.$

Taking $\delta \to 0$ , we find that ${\mathcal{P}}^{5/3}\left(S\right)=0$ , whence we see in particular that the Hausdorff dimension of S is at most 5/3.

To show that ${\mathcal{P}}^{1}\left(S\right)=0$ by a similar method, instead of the integral of $\left({\left|u\right|}^{10/3}+{\left|p\right|}^{5/3}\right)$ , we need a global quantity of dimension 1. This is furnished by Proposition 2.9. Indeed, this result allows to replace Equation (6) by

$\frac{1}{{r}_{i}}{\displaystyle {\iint}_{{Q}_{{r}_{i}/5}^{*}\left({x}_{i},{t}_{i}\right)}{\left|\nabla u\right|}^{2}}>C\left({\u03f5}_{3}\right)$ (7)

and, this way, we are led to the estimate

${{\displaystyle \sum}}^{\text{}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{r}_{i}\le C{\displaystyle {\iint}_{\cup {Q}_{{r}_{i}/5}^{*}\left({x}_{i},{t}_{i}\right)}{\left|\nabla u\right|}^{2}},$

whence we conclude that ${\mathcal{P}}^{1}\left(S\right)=0$ .

It is natural to ask if we can get a better estimate of the dimension of S. In other words, can we find $k<1$ such that ${\mathcal{P}}^{k}\left(S\right)=0$ ? Unfortunately, this question has not been answered up to now. Actually, the answer does not seem simple and is related to the possibility of demonstrating an additional estimate of the (suitable) weak solutions of order less than 1.

It is important to note that the assumption $f\in {L}^{q}{\left(Q\right)}^{3}$ with $q>5/2$ is mainly needed to prove Proposition 2.7. On the other hand, note that, in Theorem 2.6, Caffarelli, Kohn and Nirenberg chose to estimate the measure ${\mathcal{P}}^{1}$ of the set S, instead of the standard measure ${\mathcal{H}}^{1}$ . Both definitions are special cases of a construction made by Carathéodory that is detailed in [26] .

The argument used by Caffarelli, Kohn and Nirenberg is valid for any suitable solution. In the Appendix of [2] , they prove the existence of such a solution. Thus, the following holds:

Theorem 2.10 Suppose that ${u}_{0}\in V$ , $f\in {L}^{q}{\left(Q\right)}^{3}$ with $q>5/2$ and $\nabla \cdot f=0$ in Q. Then, there exists at least one suitable weak solution $\left(u,p\right)$ to the Navier-Stokes equations in Q satisfying

$u\left(t\right)\to {u}_{0}$ weakly in H as $t\to 0$ .

In addition, one has:

$\begin{array}{l}{\displaystyle {\int}_{\Omega \times \left\{t\right\}}{\left|u\right|}^{2}\varphi}+2{\displaystyle {\int}_{0}^{t}{\displaystyle {\int}_{\Omega}{\left|\nabla u\right|}^{2}\varphi}}\\ \le {\displaystyle {\int}_{\Omega}{\left|{u}_{0}\right|}^{2}\varphi \left(x,0\right)}+{\displaystyle {\int}_{0}^{t}{\displaystyle {\int}_{\Omega}\left({\left|u\right|}^{2}\left({\varphi}_{t}+\Delta \varphi \right)+\left({\left|u\right|}^{2}+2p\right)u\cdot \nabla \varphi +2\left(u\cdot f\right)\varphi \right)}}\end{array}$ (8)

for all functions $\varphi \in \mathcal{D}\left(\Omega \times \left[0,T\right]\right)$ with $\varphi \ge 0$ and $\varphi =0$ near $\partial \Omega \times \left(0,T\right)$ .

2.3. On the Existence of Suitable Weak Solutions

The existence of a suitable weak solution to Equation (2) is established in [2] by introducing a family of linear approximated problems and checking that the generalized energy inequalities are satisfied in the limit. A second proof is given in [3] , using Faedo-Galerkin approximations. In both cases, the main difficult point is passing to the limit in the term $pu\cdot \nabla \varphi $ in the right-hand side of the inequality. This requires nontrivial estimates on the pressure. In particular, Guermond [3] is able to achieve by reproducing for the discrete pressure some a priori estimates similar to the estimates of Sohr and Von Wahl in [27] .

3. Some Convergence Results

3.1. The Convergence of the Semi-Approximate Problems

In this section, we will give a new proof of Theorem 2.10. To do this, we will apply the semi-implicit Euler scheme to produce a family of approximations to the Navier-Stokes problem Equation (2). We will see that, at least for a subsequence, we have convergence to a suitable weak solution.

The scheme is the following. We take N large enough (the number of time steps) and we define the time step size $\tau :=T/N$ , the instants ${t}^{m}:=m\tau $ and the approximations

${f}^{m}:=\frac{1}{\tau}{\displaystyle {\int}_{{t}^{m-1}}^{{t}^{m}}f\left(x,t\right)\text{d}t}$ , ${u}^{m}\approx u\left(\cdot ,{t}^{m}\right)$ and ${p}^{m}\approx p\left(\cdot ,{t}^{m}\right)$ , with ${u}^{0}={u}_{0}$ and

$\{\begin{array}{l}\frac{{u}^{m+1}-{u}^{m}}{\tau}+\left({u}^{m}\cdot \nabla \right){u}^{m+1}-\Delta {u}^{m+1}+\nabla {p}^{m+1}={f}^{m+1},x\in \Omega \\ \nabla \cdot {u}^{m+1}=0,x\in \Omega ;{\displaystyle {\int}_{\Omega}{p}^{m+1}\text{d}x}=0,\\ {u}^{m+1}=0,x\in \partial \Omega ,\end{array}$ (9)

for $m=0,1,\cdots ,N-1$ .

First of all, let us check that the ${u}^{m}$ are well defined:

Lemma 3.1 The Euler scheme in Equation (9) is well defined. In other words, for every $m\ge 0$ , there exists a unique solution $\left({u}^{m+1},{p}^{m+1}\right)$ to Equation (9).

The proof is immediate by induction. We only need to note that for each

Find
$u\in V$ such that
$1/\tau \left(u,v\right)+\left(\left(w\cdot \nabla \right)u,v\right)+\left(\nabla u,\nabla v\right)=\left(g,v\right)\text{\hspace{0.17em}}\forall v\in V$ , Now, let us see that the
${u}^{m}$ are uniformly bounded in the L^{2} norm. We have:

$\left(\frac{{u}^{m+1}-{u}^{m}}{\tau},{u}^{m+1}\right)+\left(\left({u}^{m}\cdot \nabla \right){u}^{m+1},{u}^{m+1}\right)+\left(\nabla {u}^{m+1},\nabla {u}^{m+1}\right)=\left({f}^{m+1},{u}^{m+1}\right),$ (10)

which can be rewritten in the form

$\frac{1}{2}\left({\left|{u}^{m+1}\right|}^{2}-{\left|{u}^{m}\right|}^{2}\right)+\frac{1}{2}{\left|{u}^{m+1}-{u}^{m}\right|}^{2}+\tau {\left|\nabla {u}^{m+1}\right|}^{2}=\tau \left({f}^{m+1},{u}^{m+1}\right).$ (11)

Using the Cauchy-Schwarz and Young inequalities, we easily get that

$\frac{1}{2}\left({\left|{u}^{m+1}\right|}^{2}-{\left|{u}^{m}\right|}^{2}\right)+\tau {\left|\nabla {u}^{m+1}\right|}^{2}\le C{\left|{f}^{m+1}\right|}^{2}+\frac{\tau}{2}{\left|\nabla {u}^{m+1}\right|}^{2}.$ (12)

Hence,

${\left|{u}^{n+1}\right|}^{2}+\tau {\displaystyle {\sum}_{m=0}^{n}{\left|\nabla {u}^{m+1}\right|}^{2}}\le C{\left|{f}^{m+1}\right|}^{2}+{\left|{u}^{0}\right|}^{2}\le C$ (13)

for all n and, certainly, ${u}^{m}$ is uniformly bounded in H.

Using this Euler scheme, we can construct the approximate solutions of the Navier-Stokes system. More precisely, let us introduce the functions ${u}_{N}$ and ${u}_{N}^{*}$ as follows:

・ ${u}_{N}:\left[0,T\right]\mapsto V$ is the unique continuous piecewise linear function satisfying

${u}_{N}\left({t}^{m}\right)={u}^{m}$ for $m=0,1,\cdots ,N$ .

・ ${u}_{N}^{*}:\left[0,T\right]\mapsto V$ is the piecewise constant function characterized by

${u}_{N}^{*}\left(t\right)={u}^{m+1}$ in $\left({t}^{m},{t}^{m+1}\right]$ for $m=0,1,\cdots ,N-1$ .

In a similar way, we can introduce the approximate pressures ${p}_{N}^{*}$ and forces ${f}_{N}^{*}$ (again piecewise constant). The following holds:

Lemma 3.2 For any N and almost every $t\in \left(0,T\right)$ , one has

$\{\begin{array}{l}{u}_{N,t}+\left({u}_{N}^{*}(t-\tau )\cdot \nabla \right){u}_{N}^{*}-\Delta {u}_{N}^{*}+\nabla {p}_{N}^{*}={f}_{N}^{*},\\ \nabla \cdot {u}_{N}^{*}=0.\end{array}$ (14)

We can now present the main result of this section. It is related to the convergence of ${u}_{N}$ and ${u}_{N}^{*}$ towards a suitable weak solution to the Navier-Stokes equation:

Theorem 3.3 After eventual extraction of a subsequence, the functions ${u}_{N}^{*}$ converge weakly in ${L}^{2}\left(0,T;V\right)$ , weakly-* in ${L}^{\infty}\left(0,T;H\right)$ , strongly in ${L}^{2}{\left(Q\right)}^{3}$ and a.e. in Q towards a suitable weak solution to Equation (2) as $N\to +\infty $ .

For the proof of Theorem 3.3, it will be convenient to recall the following well known lemma (for instance, see the proof in [19] ):

Lemma 3.4 Let u be a function satisfying $u\in {L}^{2}\left(0,T;V\right)$ and ${u}_{t}\in {L}^{2}\left(0,T;{V}^{\prime}\right)$ . Then, u is a.e. equal to a continuous function from $\left[0,T\right]$ into H. In addition, the function $t\mapsto {\left|u\left(t\right)\right|}^{2}$ is absolutely continuous and

$\frac{\text{d}}{\text{d}t}{\left|u\left(t\right)\right|}^{2}=2\langle {u}_{t}\left(t\right),u\left(t\right)\rangle $ a.e. in $\left(0,T\right)$ ,

where $\langle \cdot ,\cdot \rangle $ denotes the duality product in ${V}^{\prime}\times V$ .

Proof of Theorem 3.3:

Let us first try to find the spaces where the functions ${u}_{N}^{*}$ , ${u}_{N}$ and ${p}_{N}^{*}$ are uniformly bounded. This is classical and very well known, but we will give the details for completeness.

Consider the Equation (11). Let us fix N and n with $0\le n\le N-1$ and let us carry out summation in m, from 0 to n. The following is obtained:

$\begin{array}{l}\frac{1}{2}{\left|{u}^{n+1}\right|}^{2}+\frac{1}{2}\underset{m=0}{\overset{n}{{\displaystyle \sum}}}{\left|{u}^{m+1}-{u}^{m}\right|}^{2}+\tau \underset{m=0}{\overset{n}{{\displaystyle \sum}}}{\left|\nabla {u}^{m+1}\right|}^{2}\\ =\tau \underset{m=0}{\overset{n}{{\displaystyle \sum}}}\left({f}^{m+1},{u}^{m+1}\right)+\frac{1}{2}{\left|{u}^{0}\right|}^{2}.\end{array}$

Obviously, this can also be written in the form

$\begin{array}{l}\frac{1}{2}{\left|{u}_{N}^{*}\left(t\right)\right|}^{2}+\frac{1}{2}\underset{m=0}{\overset{n}{{\displaystyle \sum}}}{\left|{u}^{m+1}-{u}^{m}\right|}^{2}+\underset{m=0}{\overset{n}{{\displaystyle \sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{\displaystyle {\int}_{{t}^{m}}^{{t}^{m+1}}{\left|\nabla {u}_{N}^{*}\left(t\right)\right|}^{2}\text{d}t}\\ =\underset{m=0}{\overset{n}{{\displaystyle \sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{\displaystyle {\int}_{{t}^{m}}^{{t}^{m+1}}\left({f}_{N}^{*}\left(t\right),{u}_{N}^{*}\left(t\right)\right)\text{d}t}+\frac{1}{2}{\left|{u}^{0}\right|}^{2}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\forall t\in \left({t}^{m},{t}^{m+1}\right].\end{array}$

Therefore,

$\frac{1}{2}{\left|{u}_{N}^{*}(t)\right|}^{2}+{\displaystyle {\int}_{0}^{{t}^{n+1}}{\left|\nabla {u}_{N}^{*}\left(t\right)\right|}^{2}\text{d}t}\le {\displaystyle {\int}_{0}^{{t}^{n+1}}\left({f}_{N}^{*}\left(t\right),{u}_{N}^{*}\left(t\right)\right)\text{d}t}+\frac{1}{2}{\left|{u}^{0}\right|}^{2}$ (15)

and, from the Cauchy-Schwarz and Young inequalities, we easily see that

${\left|{u}_{N}^{*}\left(t\right)\right|}^{2}+{\displaystyle {\int}_{0}^{T}{\left|\nabla {u}_{N}^{*}\left(t\right)\right|}^{2}\text{d}t}\le {\left|{u}^{0}\right|}^{2}+C{\Vert {f}_{N}^{*}\Vert}_{{L}^{2}\left(Q\right)}^{2}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\forall t\in \left[0,T\right].$

This means that

${u}_{N}^{*}$ is uniformly bounded in ${L}^{2}\left(0,T;V\right)$ and ${L}^{\infty}\left(0,T;H\right)$ . (16)

On the other hand, it can also be deduced from Equation (9) that

${\int}_{0}^{T}{\left|{u}_{N}\left(t\right)-{u}_{N}^{*}\left(t\right)\right|}^{2}\text{d}t}\le \tau \underset{m=0}{\overset{N-1}{{\displaystyle \sum}}}{\left|{u}^{m+1}-{u}^{m}\right|}^{2}\le C\tau ,$

whence ${\Vert {u}_{N}^{*}-{u}_{N}\Vert}_{{L}^{2}\left(Q\right)}^{2}\le C\tau $ .

To estimate ${u}_{N}$ , we use its definition and the fact that, for any $t\in \left({t}^{m},{t}^{m+1}\right)$ ,

$\left|{u}_{N}\left(t\right)\right|\le \left|{u}^{m}\right|+\left|{u}^{m+1}\right|$ and $\left|\nabla {u}_{N}\left(t\right)\right|\le \left|\nabla {u}^{m}\right|+\left|\nabla {u}^{m+1}\right|$ .

Accordingly, we also have that

${u}_{N}$ is uniformly bounded in ${L}^{2}\left(0,T;V\right)$ and ${L}^{\infty}\left(0,T;H\right)$ . (17)

Now, from classical interpolation results, we deduce that ${u}_{N}^{*}$ and ${u}_{N}$ is uniformly bounded in

${L}^{r}\left(0,T;{L}^{\frac{6r}{3r-4}}{\left(\Omega \right)}^{3}\right)$ $\forall r\in \left[2,+\infty \right]$ . (18)

It is well known that the estimates of Equation (17) and Equation (18) allow us to prove that the ${u}_{N}$ belong to and are uniformly bounded in the Sobolev spaces of fractional order ${H}^{\gamma}\left(0,T;H\right)$ for $0<\gamma <1/4$ ; see for example [19] . Therefore, as a consequence of Aubin-Lions’ Theorem, the ${u}_{N}$ belong to a compact set of ${L}^{2}\left(Q\right)$ .

As a consequence, at least for a subsequence (again indexed by N), we must have:

$\{\begin{array}{l}{u}_{N}\to u\text{\hspace{0.17em}}\text{weakly}\text{\hspace{0.17em}}\text{in}\text{\hspace{0.17em}}{L}^{2}\left(0,T;V\right)\text{\hspace{0.17em}}\text{and}\text{\hspace{0.17em}}\text{weakly-}\ast \text{\hspace{0.17em}}\text{in}\text{\hspace{0.17em}}{L}^{\infty}\left(0,T;H\right),\\ {u}_{N}\to u\text{\hspace{0.17em}}\text{strongly}\text{\hspace{0.17em}}\text{in}\text{\hspace{0.17em}}{L}^{2}{\left(Q\right)}^{3}\text{\hspace{0.17em}}\text{and}\text{\hspace{0.17em}}\text{a}\text{.e}\text{.}\text{\hspace{0.17em}}\text{in}\text{\hspace{0.17em}}Q.\end{array}$ (19)

This is enough to pass to the limit in Equation (14) and deduce that u is a weak solution of Equation (2). Note that it can also be assumed that

${u}_{N}\to u$ strongly in ${L}^{r}\left(0,T;{L}^{q}{\left(\Omega \right)}^{3}\right)$ for all $2<r<+\infty $ , $1\le q<6r/\left(3r-4\right)$ . (20)

To show that u is suitable, we have to give new estimates. To this purpose, we will use some regularity results that, as those in [3] , play the role of the Sohr and Wahl’s estimates in [27] .

For $0<s<1$ , the space ${H}^{s}\left(\Omega \right):={\left[{H}^{1}\left(\Omega \right),{L}^{2}\left(\Omega \right)\right]}_{s}$ can be defined by the method of real interpolation between ${H}^{1}\left(\Omega \right)$ and ${L}^{2}\left(\Omega \right)$ , i.e. the so-called K-method of Lions and Peetre [28] ; see also [29] and [30] . We will denote by ${H}_{0}^{s}$ the closure of $\mathcal{D}\left(\Omega \right)$ in ${H}^{s}\left(\Omega \right)$ . For any $s<0$ , the space ${H}^{-s}\left(\Omega \right)$ and the corresponding norm ${\Vert \text{\hspace{0.05em}}\cdot \text{\hspace{0.05em}}\Vert}_{{H}^{-s}}$ are defined by duality and, in particular,

${\Vert v\Vert}_{{H}^{-s}}:=\underset{w\in \mathcal{D}\left(\Omega \right)\backslash \left\{0\right\}}{\mathrm{sup}}\frac{\left(v,w\right)}{{w}_{{H}^{s}}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\forall v\in {L}^{2}\left(\Omega \right).$

We will look for a uniform estimate of ${u}_{N,t}$ in a space of the form ${L}^{a}\left(0,T;{H}^{-\sigma}{\left(\Omega \right)}^{3}\right)$ . This way, by applying De-Rham’s Lemma (see [31] ), we will get a bound of ${p}_{N}^{*}$ in ${L}^{a}\left(0,T;{H}^{1-\sigma}{\left(\Omega \right)}^{3}\right)$ and we will be able to take limits in the generalized energy inequality.

Note that, for all m, one has ${u}^{m}={w}^{m}+{z}^{m}$ , where the ${w}^{m}$ and the ${z}^{m}$ are respectively given by

$\{\begin{array}{l}\frac{1}{\tau}\left({w}^{m+1}-{w}^{m}\right)+A{w}^{m+1}=0,\\ {w}^{0}={u}^{0}\end{array}$ (21)

and

$\{\begin{array}{l}\frac{1}{\tau}\left({z}^{m+1}-{z}^{m}\right)+A{z}^{m+1}={F}^{m+1},\\ {z}^{0}=0,\end{array}$ (22)

where ${F}^{m+1}={f}^{m+1}-\left({u}^{m}\cdot \nabla \right){u}^{m+1}$ and A is the Stokes operator. Recall that $A:D\left(A\right)\subset H\mapsto H$ , with

$D\left(A\right)={H}^{2}{\left(\Omega \right)}^{3}\cap V,Av=P\left(-\Delta v\right)\text{\hspace{0.17em}}\text{\hspace{0.17em}}\forall v\in D(A)$

(here, $P:{L}^{2}{\left(\Omega \right)}^{3}\mapsto H$ is the orthogonal projector). Also, recall that there exists an orthogonal basis of V formed by eigenfunctions ${\xi}_{j}$ ,

$A{\xi}_{j}={\lambda}_{j}{\xi}_{j},{\xi}_{j}\in V,\left|{\xi}_{j}\right|=1,{\lambda}_{j}\nearrow \infty $

and

$D\left({A}^{r}\right)=\left\{v\in H:\underset{j\ge 1}{{\displaystyle \sum}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{\lambda}_{j}^{2r}{\left|\left(v,{\xi}_{j}\right)\right|}^{2}<+\infty \right\}$

for all $r\ge 0$ .

In the sequel, we will consider the functions ${w}_{N},{w}_{N}^{*},{z}_{N}$ and ${z}_{N}^{*}$ , whose definitions can be obtained from the ${z}^{m}$ and the ${w}^{m}$ in a way similar to ${u}_{N}$ and ${u}_{N}^{*}$ .

First, note that

${w}^{m}={\left(Id.+\tau A\right)}^{-m}{u}_{0}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\forall m=0,1,\cdots ,N,$

whence

$\begin{array}{c}{\Vert {w}_{N,t}\Vert}_{{L}^{2}\left(Q\right)}^{2}=\tau \underset{m=0}{\overset{N-1}{{\displaystyle \sum}}}{\left|A{\left(Id.+\tau A\right)}^{-\left(m+1\right)}{u}_{0}\right|}^{2}\\ =\tau \underset{m=0}{\overset{N-1}{{\displaystyle \sum}}}\text{\hspace{0.17em}}\underset{j\ge 1}{{\displaystyle \sum}}\frac{{\lambda}_{j}^{2}}{{\left(1+\tau {\lambda}_{j}\right)}^{2\left(m+1\right)}}{\left|\left({u}_{0},{\xi}_{j}\right)\right|}^{2}\\ =\tau \underset{j\ge 1}{{\displaystyle \sum}}\left(\underset{m=0}{\overset{N-1}{{\displaystyle \sum}}}\frac{1}{{\left(1+\tau {\lambda}_{j}\right)}^{2\left(m+1\right)}}\right){\lambda}_{j}^{2}{\left|\left({u}_{0},{\xi}_{j}\right)\right|}^{2}\\ \le {\displaystyle {\sum}_{j\ge 1}\frac{{\lambda}_{j}^{2}}{{\left(1+\tau {\lambda}_{j}\right)}^{2\left(m+1\right)}}{\left|\left({u}_{0},{\xi}_{j}\right)\right|}^{2}}={\left|{u}_{0}\right|}^{2}.\end{array}$ (23)

Therefore,

${w}_{N,t}$ and $A{w}_{N}^{*}$ are uniformly bounded in ${L}^{2}{\left(Q\right)}^{3}$ . (24)

Let us now see what can be said of ${z}_{N,t}$ and $A{z}_{N}^{*}$ . For all $m\ge 1$ , we have

${z}^{m}=\underset{l=1}{\overset{m}{{\displaystyle \sum}}}{\left(Id.+\tau A\right)}^{-\left(m+1-l\right)}{F}^{l}$

Let $s,\sigma \in \left(0,1\right)$ be such that with $\sigma >s$ . Then

${\Vert A{z}^{m+1}\Vert}_{{H}^{-\sigma}}\le \tau \underset{l=1}{\overset{m+1}{{\displaystyle \sum}}}{\Vert A{\left(Id.+\tau A\right)}^{-\left(m+2-l\right)}\Vert}_{\mathcal{L}\left({H}^{-s};{H}^{-\sigma}\right)}{\Vert {F}^{l}\Vert}_{{H}^{-s}}=\tau \underset{l=1}{\overset{m}{{\displaystyle \sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{a}_{m-l}{b}_{l},$

where the ${a}_{n}$ and the ${b}_{l}$ are given by

${a}_{n}={\Vert A{\left(Id.+\tau A\right)}^{-\left(m+2-l\right)}\Vert}_{\mathcal{L}\left({H}^{-s};{H}^{-\sigma}\right)},{b}_{l}={\Vert {F}^{l}\Vert}_{{H}^{-s}}.$

We will apply the following result, that must be viewed as a discrete version of the well known Young inequality for convolution products:

Lemma 3.5 Let us assume that $k\ge 1$ , $a\in {l}^{p}$ and $b\in {l}^{q}$ . Then, if $r\in \left[1,+\infty \right]$ and

$\frac{1}{p}+\frac{1}{q}=1+\frac{1}{r},$

one has

${\left({\displaystyle {\sum}_{n=1}^{k}{\left|{\displaystyle {\sum}_{l=1}^{n}{a}_{n-l}{b}_{l}}\right|}^{r}}\right)}^{1/r}\le {\left({\displaystyle {\sum}_{n=1}^{k}{\left|{a}_{n}\right|}^{p}}\right)}^{\frac{1}{p}}{\left({\displaystyle {\sum}_{n=1}^{k}{\left|{b}_{n}\right|}^{q}}\right)}^{1/q}$ (25)

for all $k\ge 1$ .

The proof of this result can be found in [32] . Using Lemma 2.5 with $r=a$ , $p=1$ and $q=a$ , we find that

${\left(\tau \underset{n=1}{\overset{N}{{\displaystyle \sum}}}{\Vert A{z}^{n+1}\Vert}_{{H}^{-\sigma}}^{a}\right)}^{1/a}\le {\left({\tau}^{1+a}\underset{n=1}{\overset{N}{{\displaystyle \sum}}}{\left|\underset{l=1}{\overset{n}{{\displaystyle \sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{a}_{n-l}{b}_{l}\right|}^{a}\right)}^{\frac{1}{a}}\le \left(\underset{n=1}{\overset{N}{{\displaystyle \sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{a}_{n}\right){\left(\underset{n=1}{\overset{N}{{\displaystyle \sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{b}_{n}^{a}\right)}^{1/a}.$

From the estimates in Equation (16) already obtained for
${u}_{N}^{*}$ , it is immediate that, for any
$a\in \left[1,2\right]$ ,
${F}_{N}^{*}$ is uniformly bounded in
${L}^{a}\left(0,T;{L}^{3a/\left(4a-2\right)}{\left(\Omega \right)}^{3}\right)$ and, consequently, also in
${L}^{a}\left(0,T;{H}^{-\left(5a-4\right)/\left(2a\right)}{\left(\Omega \right)}^{3}\right)$ . Thus, choosing

${\Vert {F}_{N}^{*}\Vert}_{{L}^{a}\left(0,T;{H}^{-s}\right)}={\left(\tau \underset{m=1}{\overset{N}{{\displaystyle \sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{b}_{m}^{a}\right)}^{1/a}\le C\left(a\right).$

On the other hand, for any smooth z, one has

${\Vert A{\left(Id.+\tau A\right)}^{-\left(n+2\right)}z\Vert}_{{H}^{-\sigma}}^{2}=\underset{j\ge 1}{{\displaystyle \sum}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{\lambda}_{j}^{-\sigma}\frac{{\lambda}_{j}^{2}}{{\left(1+\tau {\lambda}_{j}\right)}^{2\left(n+2\right)}}{\left|z,{\xi}_{j}\right|}^{2}\le \left[\underset{j}{\mathrm{sup}}\frac{{\lambda}_{j}^{2\left(1-\u03f5\right)}}{{\left(1+\tau {\lambda}_{j}\right)}^{2\left(n+2\right)}}\right]{\Vert z\Vert}_{{H}^{-s}}^{2},$

where $\u03f5=\left(\sigma -s\right)/2$ . Therefore, recalling the definition of the ${a}_{n}$ , we deduce that

${a}_{n}\le \frac{C\left(\u03f5\right)}{{\left(n\tau \right)}^{1-\u03f5}},\text{\hspace{0.17em}}\text{\hspace{0.17em}}\tau \underset{n=1}{\overset{m}{{\displaystyle \sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{a}_{n}\le C\left(\u03f5\right){\displaystyle {\int}_{0}^{T}\frac{\text{d}s}{{s}^{1-\u03f5}}}\le C(\u03f5)$

and, finally,

${\Vert A{z}_{N}^{*}\Vert}_{{L}^{a}\left(0,T;{H}^{-\sigma}\right)}\le C\left({\displaystyle {\sum}_{m=1}^{N}\tau {a}_{m}}\right){\left(\tau {\displaystyle {\sum}_{m=1}^{N}{\Vert {F}^{m}\Vert}_{{H}^{-s}}^{a}}\right)}^{\frac{1}{a}}\le {\Vert C{F}_{N}^{*}\Vert}_{{L}^{a}\left(0,T;{H}^{-s}\right)}$ (26)

Note that this estimate is valid for all $a\in \left[1,2\right]$ , with $s=\left(5a-4\right)/\left(2a\right)$ and $\sigma \in \left(0,1\right)$ , $\sigma >s$ .

Obviously, the same estimate is valid for ${\Vert {z}_{N,t}\Vert}_{{L}^{a}\left(0,T;{H}^{-\sigma}\right)}$ . This proves that

$\begin{array}{l}{z}_{N,t}\text{\hspace{0.17em}}\text{and}\text{\hspace{0.17em}}A{z}_{N}^{*}\text{\hspace{0.17em}}\text{areuniformlyboundedin}\text{\hspace{0.17em}}{L}^{a}\left(0,T;{H}^{-\sigma}{\left(\Omega \right)}^{3}\right)\\ \forall a\in \left[0,1\right],\forall \sigma >s=\left(5a-4\right)/\left(2a\right).\end{array}$ (27)

In view of Equation (24) and Equation (27) and recalling that ${u}_{N}^{*}={w}_{N}^{*}+{z}_{N}^{*}$ , it follows that $A{u}_{N}^{*}$ and ${u}_{N,t}^{*}$ are also uniformly bounded in ${L}^{a}\left(0,T;{H}^{-\sigma}{\left(\Omega \right)}^{3}\right)$ .

Now, from De-Rham’s Lemma [31] , we see that ${p}_{N}^{*}$ is uniformly bounded in ${L}^{a}\left(0,T;{H}^{1-\sigma}\left(\Omega \right)\right)$ , which is continuously embedded in ${L}^{a}\left(0,T;{L}^{6/\left(1+2\sigma \right)}\left(\Omega \right)\right)$ . In particular, for $a=5/3$ , we have $s=13/10$ and we can take σ as close as desired to s, which gives $6/\left(1+2\sigma \right)$ as close as desired to 5/3.

As a consequence of these estimates, after extracting a new sequence (if this is needed), we see that

${p}_{N}^{*}$ converges weakly in ${L}^{5/3}\left(0,T;{L}^{\beta}\left(\Omega \right)\right)$ $\forall \beta <5/3$ . (28)

Let us check that the local energy inequality holds for u and p.

If we multiply the Equation (14) by the function ${u}_{N}^{*}\varphi $ , where $\varphi \in {C}_{0}^{\infty}\left(\Omega \times \left[0,T\right]\right)$ is non-negative and we integrate in space, we have:

$\begin{array}{l}{\displaystyle {\int}_{\Omega}{u}_{N,t}\cdot {u}_{N}^{*}\varphi}+{\displaystyle {\int}_{\Omega}\left({u}_{N}^{*}\left(t-\tau \right)\cdot \nabla \right){\left|{u}_{N}^{*}\right|}^{2}\varphi}+{\displaystyle {\int}_{\Omega}\left(-\Delta {u}_{N}^{*}\right)\cdot {u}_{N}^{*}\varphi}+{\displaystyle {\int}_{\Omega}\nabla {p}_{N}^{*}\cdot {u}_{N}^{*}\varphi}\\ ={\displaystyle {\int}_{\Omega}{f}_{N}^{*}\cdot {u}_{N}^{*}\varphi}.\end{array}$ (29)

If $t\in \left[0,T\right]$ , there exists n such that $t\in \left({t}^{n},{t}^{n+1}\right]$ and then, using Lemma 3.4, one has:

$\begin{array}{c}{\displaystyle {\int}_{\Omega}{u}_{N,t}\cdot {u}_{N}^{*}\varphi}={\displaystyle {\int}_{\Omega}{u}_{N,t}\cdot {u}_{N}\varphi}+{\displaystyle {\int}_{\Omega}{u}_{N,t}\cdot \left({u}_{N}^{*}-{u}_{N}\right)\varphi}\\ =\frac{1}{2}\frac{\text{d}}{\text{d}t}{\displaystyle {\int}_{\Omega}{\left|{u}_{N}\right|}^{2}\varphi}-\frac{1}{2}{\displaystyle {\int}_{\Omega}{\left|{u}_{N}\right|}^{2}{\varphi}_{t}}+{\displaystyle {\int}_{\Omega}{u}_{N,t}\cdot \left({u}_{N}^{\text{*}}-{u}_{N}\right)\varphi}\text{\hspace{0.05em}}\text{.}\end{array}$

Note moreover that ${\int}_{\Omega}{u}_{N,t}\cdot \left({u}_{N}^{*}-{u}_{N}\right)\varphi}\ge 0$ , because ${u}_{N}^{*}-{u}_{N}$ is by definition equal to $\left({t}^{n+1}-t\right){u}_{N,t}$ .

• Also,

${\int}_{\Omega}\left({u}_{N}^{*}\left(t-\tau \right)\cdot \nabla \right){\left|{u}_{N}^{*}\right|}^{2}\varphi}=-\frac{1}{2}{\displaystyle {\int}_{\Omega}{u}_{N}^{*}\left(t-\tau \right){\left|{u}_{N}^{*}\right|}^{2}\cdot \nabla \varphi},$

• On the other hand,

${\int}_{\Omega}\left(-\Delta {u}_{N}^{*}\right)\cdot {u}_{N}^{*}\varphi}={\displaystyle {\int}_{\Omega}{\left|\nabla {u}_{N}^{*}\right|}^{2}\varphi}-\frac{1}{2}{\displaystyle {\int}_{\Omega}{\left|{u}_{N}^{*}\right|}^{2}\Delta \varphi}.$

• Finally,

${\int}_{\Omega}\nabla {p}_{N}^{*}\cdot {u}_{N}^{*}\varphi}=-{\displaystyle {\int}_{\Omega}{p}_{N}^{*}\cdot {u}_{N}^{*}\nabla \varphi}.$

Consequently, we see from Equation (29) that

$\begin{array}{l}\frac{1}{2}\frac{\text{d}}{\text{d}t}{\displaystyle {\int}_{\Omega}{\left|{u}_{N}\right|}^{2}\varphi}+{\displaystyle {\int}_{\Omega}{\left|\nabla {u}_{N}^{*}\right|}^{2}\varphi}\le \frac{1}{2}{\displaystyle {\int}_{\Omega}{\left|{u}_{0}\right|}^{2}\varphi}+{\displaystyle {\iint}_{\Omega \times \left(0,t\right)}{\left|{u}_{N}\right|}^{2}{\varphi}_{t}}\\ \text{\hspace{0.05em}}+{\displaystyle {\iint}_{\Omega \times \left(0,t\right)}\left[\left(\frac{1}{2}{u}_{N}^{*}\left(t-\tau \right){\left|{u}_{N}^{*}\right|}^{2}+{p}_{N}^{*}{u}_{N}^{*}\right)\cdot \nabla \varphi +\frac{1}{2}{\left|{u}_{N}^{*}\right|}^{2}\Delta \varphi +{f}_{N}^{*}\cdot {u}_{N}^{*}\varphi \right]}\end{array}$

If we integrate in time, we find that

$\begin{array}{l}{\displaystyle {\int}_{\Omega}{\left|{u}_{N}\right|}^{2}\varphi}+2{\displaystyle {\iint}_{\Omega \times \left(0,t\right)}{\left|\nabla {u}_{N}^{*}\right|}^{2}\varphi}\le {\displaystyle {\int}_{\Omega}{\left|{u}_{0}\right|}^{2}\varphi}+{\displaystyle {\iint}_{\Omega \times \left(0,t\right)}{\left|{u}_{N}\right|}^{2}{\varphi}_{t}}\\ \text{\hspace{0.05em}}+{\displaystyle {\iint}_{\Omega \times \left(0,t\right)}\left[\left({u}_{N}^{*}\left(t-\tau \right){\left|{u}_{N}^{*}\right|}^{2}+2{p}_{N}^{*}{u}_{N}^{*}\right)\cdot \nabla \varphi +{\left|{u}_{N}^{*}\right|}^{2}\Delta \varphi +2{f}_{N}^{*}\cdot {u}_{N}^{*}\varphi \right]}\end{array}$

Thanks to the energy estimates Equation (16), we can take the lower limit in the left-hand side. On the other hand, thanks to Equation (19), Equation (20) and Equation (28), we can take limits in all the terms in the right; for example, since ${u}_{N}^{*}$ converges strongly in ${L}^{5/2}\left(0,T;{L}^{19/5}{\left(\Omega \right)}^{3}\right)$ and ${p}_{N}^{*}$ converges weakly in ${L}^{5/3}\left(0,T;{L}^{14/9}\left(\Omega \right)\right)$ , we see that ${p}_{N}^{*}{u}_{N}^{*}\cdot \nabla \varphi $ converges weakly in ${L}^{1}\left(Q\right)$ towards $pu\cdot \nabla \varphi $ .

The final result is that

$\begin{array}{l}2{\displaystyle {\iint}_{\Omega \times \left(0,t\right)}{\left|\nabla u\right|}^{2}\varphi}\\ \le \underset{\Omega}{{\displaystyle \int}}{\left|{u}_{0}\right|}^{2}\varphi +{\displaystyle {\iint}_{\Omega \times \left(0,t\right)}\left({\left|u\right|}^{2}\left({\varphi}_{t}+\Delta \varphi \right)+\left({\left|u\right|}^{2}+2p\right)u\cdot \nabla \varphi +2\left(u\cdot f\right)\varphi \right)},\end{array}$ (30)

as desired.

3.2. The Convergence of the Fully Discretized Problems

In this section, we will argue as in [3] and we will check that the approximate solutions obtained via the semi-implicit Euler discrete scheme, used together with an appropriate approximation in space, converge to a suitable solution to Equation (2).

As before, let us introduce N, $\tau :=T/N$ and the ${t}^{m}:=m\tau $ . We will also consider two families of finite dimensional spaces ${\left\{{X}_{h}\right\}}_{h>0}$ and ${\left\{{P}_{h}\right\}}_{h>0}$ with the ${X}_{h}\subset {H}_{0}^{1}{\left(\Omega \right)}^{3}$ and the ${P}_{h}\subset {L}^{2}\left(\Omega \right)$ such that

$\{\begin{array}{l}{\mathrm{inf}}_{{v}_{h}\in {X}_{h}}{\Vert v-{v}_{h}\Vert}_{{H}^{1}}\to 0\text{\hspace{0.17em}}\text{\hspace{0.17em}}\forall v\in {H}_{0}^{1}{\left(\Omega \right)}^{3},\\ {\mathrm{inf}}_{{q}_{h}\in {P}_{h}}{\Vert q-{q}_{h}\Vert}_{{L}^{2}}\to 0\text{\hspace{0.17em}}\text{\hspace{0.17em}}\forall q\in {L}^{2}\left(\Omega \right),\end{array}$ (31)

and the $\left({X}_{h},{P}_{h}\right)$ are uniformly compatible, in the sense that there exists a constant $\mu >0$ independent of h such that the following inf-sup conditions are satisfied:

${\mathrm{inf}}_{{q}_{h}\in {P}_{h}\backslash \left\{0\right\}}{\mathrm{sup}}_{{v}_{h}\in {X}_{h}\backslash \left\{0\right\}}\frac{\left(\nabla {q}_{h},{v}_{h}\right)}{{\Vert {v}_{h}\Vert}_{{H}^{-1-s}}{\Vert {q}_{h}\Vert}_{{H}^{s}}}\ge \mu .$ (32)

Now, we consider the approximations

${u}_{h}^{m}=u\left(\cdot ,{t}^{m}\right)\in {X}_{h},{p}_{h}^{m}=p\left(\cdot ,{t}^{m}\right)\in {P}_{h}$ ,

with ${u}_{h}^{0}={u}_{0h}$ (the orthogonal projection of ${u}_{0}$ on ${X}_{h}$ ) and

$\{\begin{array}{l}\left(\frac{{u}_{h}^{m+1}-{u}_{h}^{m}}{\tau},{v}_{h}\right)+\left(\left({u}_{h}^{m}\cdot \nabla \right){u}_{h}^{m+1},{v}_{h}\right)+\left(\nabla {u}_{h}^{m+1},\nabla {v}_{h}\right)+\left(\nabla {p}_{h}^{m+1},{v}_{h}\right)\\ =\left({f}_{h}^{m+1},{v}_{h}\right),\forall {v}_{h}\in {X}_{h},\\ \left({q}_{h},\nabla \cdot {u}_{h}^{m+1}\right)=0,\\ \left({u}_{h}^{m+1},{p}_{h}^{m+1}\right)\in \left({X}_{h},{P}_{h}\right),\end{array}$ (33)

for $m=0,1,\cdots ,N-1$ . The following result, which is a consequence of Equation (32), gives coherence to our scheme:

As before, the ${u}_{h}^{m}$ and ${p}_{h}^{m}$ serve to construct approximate solutions to the Navier-Stokes system. Thus, we define the functions ${u}_{N,h},{u}_{N,h}^{*},{p}_{N,h}^{*}$ etc. similarly to ${u}_{N},{u}_{N}^{*},{p}_{N}^{*}$ etc. The main result of this section is the following:

Theorem 3.7 After eventual extraction of a subsequence, the functions ${u}_{N,h}^{*}$ converge weakly in ${L}^{2}\left(0,T;V\right)$ , weakly-* in ${L}^{\infty}\left(0,T;H\right)$ , strongly in ${L}^{2}{\left(Q\right)}^{3}$ and a.e. in Q towards a suitable weak solution to Equation (2) as $N\to +\infty $ and $h\to 0$ .

Sketch of the proof:

Arguing as in the proof of the Theorem 3.3, it can be seen that the ${u}_{N,h}$ and the ${u}_{N,h}^{*}$ are uniformly bounded in ${L}^{2}\left(0,T;V\right)$ and ${L}^{\infty}\left(0,T;H\right)$ and, furthermore, ${\Vert {u}_{N,h}-{u}_{N,h}^{*}\Vert}_{{L}^{2}\left(Q\right)}^{2}\le C\tau $ . As in [19] , we can also prove that the ${u}_{N,h}$ are uniformly bounded in ${H}^{\gamma}\left(0,T;H\right)$ for any $\gamma \in \left(0,1/4\right)$ . Consequently, at least for a subsequence (still indexed with N and h), one has:

$\{\begin{array}{l}{u}_{N,h}\to u\text{\hspace{0.17em}}\text{weakly}\text{\hspace{0.17em}}\text{in}\text{\hspace{0.17em}}{L}^{2}\left(0,T;V\right)\text{\hspace{0.17em}}\text{andweakly-}\ast \text{\hspace{0.17em}}\text{in}\text{\hspace{0.17em}}{L}^{\infty}\left(0,T;H\right),\\ {u}_{N,h}\to u\text{\hspace{0.17em}}\text{strongly}\text{\hspace{0.17em}}\text{in}\text{\hspace{0.17em}}{L}^{2}{\left(Q\right)}^{3}\text{\hspace{0.17em}}\text{and}\text{\hspace{0.17em}}\text{a}\text{.e}\text{.}\text{\hspace{0.17em}}\text{in}\text{\hspace{0.17em}}Q.\end{array}$ (34)

Also,

${u}_{N,h}\to u$ strongly in ${L}^{r}\left(0,T;{L}^{q}{\left(\Omega \right)}^{3}\right)$ for all $2<r<+\infty $ , $1\le q<6r/\left(3r-4\right)$ . (35)

As before, this is enough to pass to the limit and deduce that u is a weak solution of Equation (2).

In order to prove that u is suitable, we can argue as in Guermond [3] . Here, we need the spaces ${\stackrel{\u02dc}{H}}_{0}^{s}\left(\Omega \right):={\left[{L}^{2}\left(\Omega \right),{H}_{0}^{1}\left(\Omega \right)\right]}_{s}$ for $s\in \left(0,1\right)$ and their dual spaces ${\stackrel{\u02dc}{H}}_{0}^{-s}\left(\Omega \right)$ .

The following estimates are established in [3] :

・ For any $\alpha \in \left[1/4,1/2\right]$ and any $\delta <\stackrel{\xaf}{\delta}\left(\alpha \right)=2\left(1+\alpha \right)/5$ , one has

${\Vert {u}_{N,h,t}\Vert}_{{H}^{\delta -1}\left(0,T;{\stackrel{\u02dc}{H}}^{-\alpha}\right)}+{\Vert {u}_{N,h}^{*}\Vert}_{{H}^{\delta}\left(0,T;{\stackrel{\u02dc}{H}}^{-\alpha}\right)}\le C\left(\alpha \right),$

・ For any and $s\in \left[1/2,7/10\right]$ any, one also has

${\Vert {u}_{N,h,t}\Vert}_{{H}^{-r}\left(0,T;{\stackrel{\u02dc}{H}}^{-s}\right)}+{\Vert {p}_{N,h}^{*}\Vert}_{{H}^{-r}\left(0,T;{\stackrel{\u02dc}{H}}^{1-s}\right)}\le C\left(s\right),$

As a consequence, it can be assumed that the ${p}_{N,h}^{*}$ converge weakly (for instance) in ${H}^{-r}\left(0,T;{H}^{3/8}\left(\Omega \right)\right)$ for all $r>7/16$ and the ${u}_{N,h}^{*}$ converge strongly in ${H}^{\delta}\left(0,T;{\stackrel{\u02dc}{H}}^{-\alpha}{\left(\Omega \right)}^{3}\right)$ for all $\alpha <3/8$ and $\delta <11/20$ . This is sufficient to ensure that ${p}_{N,h}^{*}{u}_{N,h}^{*}\cdot \nabla \varphi $ converges weakly in ${L}^{1}{\left(Q\right)}^{3}$ towards $pu\cdot \nabla \varphi $ .

Hence, arguing as in the final part of the proof of Theorem 3.3, it is not difficult to check that the limit $\left(u,p\right)$ of the $\left({u}_{N,h}^{*},{p}_{N,h}^{*}\right)$ is a suitable weak solution to Equation (2).

This ends the proof.

4. Some Additional Comments and Questions

4.1. The Same Results Hold for the Boussinesq System

The Boussinesq system is the following:

$\{\begin{array}{l}{u}_{t}-\Delta u+\left(u\cdot \nabla \right)u+\nabla p=f+\theta k,\left(x,t\right)\in Q,\\ \nabla \cdot u=0,\left(x,t\right)\in Q,\\ {\theta}_{t}+u\cdot \nabla \theta -\Delta \theta =g,\left(x,t\right)\in Q,\\ u\left(x,t\right)=0,\theta \left(x,t\right)=0,\left(x,t\right)\in \Sigma ,\\ u\left(x,0\right)={u}_{0}\left(x\right),\theta \left(x,0\right)={\theta}_{0}\left(x\right),x\in \Omega .\end{array}$ (36)

We assume here that

${u}_{0}\in V,{\theta}_{0}\in {H}_{0}^{1}\left(\Omega \right),f\in {L}^{2}\left(0,T;{H}^{-1}{\left(\Omega \right)}^{3}\right),k\in {\mathbb{R}}^{3}$ and $g\in {L}^{2}\left(0,T;{L}^{2}\left(\Omega \right)\right)$ .(37)

As in [4] , we can speak of weak solutions to Equation (36) and, also, of suitable weak solutions to the previous equations in any set of the form $D=G\times \left(a,b\right)$ , with $G\subset {\mathbb{R}}^{3}$ a connected open set.

The results in Section 3 can be extended to this framework. Thus, we can for instance consider the semi-implicit Euler scheme

$\{\begin{array}{l}\frac{{u}^{m+1}-{u}^{m}}{\tau}+\left({u}^{m}\cdot \nabla \right){u}^{m+1}-\Delta {u}^{m+1}+\nabla {p}^{m+1}={f}^{m+1}+{\theta}^{m+1}k,\\ \nabla \cdot {u}^{m+1}=0,\\ \frac{{\theta}^{m+1}-{\theta}^{m}}{\tau}+{u}^{m}\cdot \nabla {\theta}^{m+1}-\Delta {\theta}^{m+1}={g}^{m+1}\end{array}$ (38)

and prove that, at least for a subsequence, the corresponding $\left({u}_{N},{p}_{N},{\theta}_{N}\right)$ and $\left({u}_{N}^{*},{p}_{N}^{*},{\theta}_{N}^{*}\right)$ converge, in an appropiate sense, to a suitable weak solution $\left(u,p,\theta \right)$ .

4.2. Possible Extensions to Other Systems

It would be interesting to prove similar results to Theorems 3.3 and 3.7 for the solutions to the variable-density Navier-Stokes equations

$\{\begin{array}{l}{\rho}_{t}+\nabla \cdot \left(\rho u\right)=0,\left(x,t\right)\in Q,\\ \rho \left({u}_{t}+\left(u\cdot \nabla \right)u\right)-\Delta u+\nabla p=\rho f,\left(x,t\right)\in Q\\ \nabla \cdot u=0,\left(x,t\right)\in Q,\\ u\left(x,t\right)=0,\left(x,t\right)\in \Sigma ,\\ u\left(x,0\right)={u}_{0}\left(x\right),\rho \left(x,0\right)={\rho}_{0}\left(x\right),x\in \Omega .\end{array}$ (39)

However, this is not clear at present. Note that the “reasonable” definition of a suitable weak solution should involve the following property: for any $\varphi \in D\left(Q\right)$ with $\varphi \ge 0$ ,

$2{\displaystyle {\iint}_{D}}{\left|\nabla u\right|}^{2}\varphi \le {\displaystyle {\iint}_{D}\left({\left|u\right|}^{2}\left(\rho {\varphi}_{t}+\Delta \varphi \right)+\left(\rho {\left|u\right|}^{2}+2p\right)\left(u\cdot \nabla \varphi \right)+2\rho \left(u\cdot f\right)\varphi \right)}.$

But, unfortunately, the apparent lack of regularity of p makes it difficult to prove this.

4.3. Extensions to Other Approximation Schemes for the Navier-Stokes Equations

As we already said, Theorems 3.3 and 3.7 can be adapted to many other time approximation schemes. Among then, let us simply recall the following:

・ Crank-Nicholson scheme:

$\frac{{u}^{m+1}-{u}^{m}}{\tau}+\left({u}^{m}\cdot \nabla \right){u}^{m+1}-\Delta \left(\frac{{u}^{m+1}+{u}^{m}}{2}\right)+\nabla {p}^{m+1}={f}^{m+1},\nabla \cdot {u}^{m+1}=0.$

・ Gear scheme:

$\frac{3{u}^{m+1}-4{u}^{m}+{u}^{m-1}}{2\tau}+\left({u}^{m}\cdot \nabla \right){u}^{m+1}-\Delta {u}^{m+1}+\nabla {p}^{m+1}={f}^{m+1},\nabla \cdot {u}^{m+1}=0.$

θ-scheme: For α and β such that $0<\alpha ,\beta <1$ and $\alpha +\beta =1$ , we compute $\left({u}^{n+\theta},{p}^{n+\theta}\right)$ , then ${u}^{n+1-\theta}$ and finally $\left({u}^{n+1},{p}^{n+1}\right)$ as follows:

$\frac{{u}^{n+\theta}-{u}^{n}}{\theta \Delta t}-\alpha \nu \Delta {u}^{n+\theta}+\nabla {p}^{n+\theta}={f}^{n+\theta}+\beta \nu \Delta {u}^{n}-\left({u}^{n}\cdot \nabla \right){u}^{n},\nabla \cdot {u}^{n+\theta}=0,$

$\frac{{u}^{n+1-\theta}-{u}^{n+\theta}}{\left(1-2\theta \right)\Delta t}-\beta \nu \Delta {u}^{n+1-\theta}+\left({u}^{n+1-\theta}\cdot \nabla \right){u}^{n+1-\theta}={f}^{n+\theta}+\alpha \nu \Delta {u}^{n+\theta}-\nabla {p}^{n+\theta}.$

$\frac{{u}^{n+1}-{u}^{n+1-\theta}}{\theta \Delta t}-\alpha \nu \Delta {u}^{n+1}+\nabla {p}^{n+1}={f}^{n+1}+\beta \nu \Delta {u}^{n+1-\theta}-\left({u}^{n+1-\theta}\cdot \nabla \right){u}^{n+1-\theta},\nabla \cdot {u}^{n+1}=0$ .

As we mentioned in Section 1, it would be interesting to establish an analog of Propositions 2.7 and 2.9 for a family of approximated solutions. This should help to detect or discard the occurrence of singular points just observing the results of appropriate numerical experiments.

Acknowledgements

E.F.-C. was partially supported by MINECO (Spain), Grant MTM2016-76990-P. I.M.-G. was partially supported by CEI (Junta de Andalucía, Spain), Grant FQM-131, VI PPIT (Universidad de Sevilla, Spain).

References

[1] Scheffer, V. (1977) Hausdorff the Navier-Stokes Equations. Communications in Mathematical Physics, 55, 97-112.

https://doi.org/10.1007/BF01626512

[2] Caffarelli, L., Kohn, R. and Nirenberg, L. (1982) Partial Regularity of Suitable Weak Solutions of the Navier-Stokes Equations. Communications on Pure and Applied Mathematics, 35, 771-831.

https://doi.org/10.1002/cpa.3160350604

[3] Guermond, J.-L. (2007) Faedo-Galerkin Weak Solutions of the Navier-Stokes Equations with Dirichlet Boundary Conditions Are Suitable. Journal de Mathématiques Pures et Appliquées, 88, 87-106.

https://doi.org/10.1016/j.matpur.2007.04.009

[4] Boling, G. and Guangwei, Y. (1996) On the Suitable Weak Solutions for the Cauchy-Problem of the Boussinesq Equations. Nonlinear Analysis Theory, Methods and Applications, 26, 1367-1385.

[5] Seregin, G.A. (2007) On the Local Regularity of Suitable Weak Solutions of the Navier-Stokes Equations. Uspekhi Matematicheskikh Nauk, 62, 149-168.

[6] Vasseur, A.F. (2007) A New Proof of Partial Regularity of Solutions to Navier-Stokes Equations. Nonlinear Differential Equations and Applications NoDEA, 14, 753-785.

[7] De Giorgi, E. (1957) Sulla differenziabilità e l’analiticità delle estremali degli integrali multipli regolari. Mem. Accad. Sci. Torino. Cl. Sci. Fis. Mat. Nat., Vol. 3, 25-43.

[8] Wolf, J. (2015) On the Local Regularity of Suitable Weak Solutions to the Generalized Navier-Stokes Equations. Annali dell’Universita di Ferrara, 61, 149-171.

https://doi.org/10.1007/s11565-014-0203-6

[9] Choe, H.J., Wolf, J. and Yang, M. (2018) A New Local Regularity Criterion for Suitable Weak Solutions of the Navier-Stokes Equations in Terms of the Velocity Gradient. Mathematische Annalen, 370, 629-647.

[10] Buckmaster, T. and Vicol, V. (2017) Nonuniqueness of Weak Solutions to the Navier-Stokes Equations.

[11] Dean, E.J. and Glowinski, R. (1993) On Some Finite Element Methods for the Numerical Simulation for Incompressible Viscous Flow. Cambridge University Press, New York, 109-150.

[12] Douglas, J. and Rachford, H.H. (1956) On the Solution of the Heat Conduction Problem in 2 and 3 Space Variables. Transactions of the American Mathematical Society, 82, 421-439.

https://doi.org/10.1090/S0002-9947-1956-0084194-4

[13] Glowinski, R. (2013) Numerical Methods for Fluids (Part 3). North-Holland, Vol. IX.

[14] Peaceman, D.H. and Rachford, H.H. (1955) The Numerical Solution of Parabolic and Elliptic Differential Equations. Journal of the Society for Industrial and Applied Mathematics, 3, 28-41.

https://doi.org/10.1137/0103003

[15] Da Veiga, H.B. (1985) On the Suitable Weak Solution to the Navier-Stokes Equations in the Whole Space. Journal de Mathématiques Pures et Appliquées, 64, 77-86.

[16] Berselli, L.C. and Spirito, S. (2017) Suitable Weak Solutions to the 3D Navier-Stokes Equations Are Constructed with the Voigt Approximation. Journal of Differential Equations, 262, 3285-3316.

[17] Berselli, L.C. and Spirito, S. (2018) On the Construction of Suitable Weak Solutions to the 3D Navier-Stokes Equations in a Bounded Domain by an Artificial Compressibility Method. Communications in Contemporary Mathematics, 20, Article ID: 1650064.

[18] Lin, F.H. (1998) A New Proof of the Caffarelli-Kohn-Nirenberg Theorem. Communications on Pure and Applied Mathematics, 51, 241-257.

https://doi.org/10.1002/(SICI)1097-0312(199803)51:3<241::AID-CPA2>3.0.CO;2-A

[19] Temam, R. (1977) Navier-Stokes Equations. Theory and Numerical Analysis. Studies in Mathematics and Applications, North-Holland Publishing Co., Amsterdam, 2.

[20] Lions, P.-L. (1996) Mathematical Topics in Fluid Mechanics, 1: Incompresible Models. Clarendon Press, Oxford.

[21] Serrin, J. (1962) On the Interior Regularity of Weak Solutions of the Navier-Stokes Equations. Archive for Rational Mechanics and Analysis, 9, 187-195.

https://doi.org/10.1007/BF00253344

[22] Wolf, J. (2010) A New Criterion for Partial Regularity of Suitable Weak Solutions to the Navier-Stokes Equations. In: Rannacher, A.S.R., Ed., Advances in Mathematical Fluid Mechanics, Springer, Berlin, 613-630.

[23] Neustupa, J. and Penel, P. (1999) Regularity of a Suitable Weak Solution to the Navier-Stokes Equations as a Consequence of Regularity of One Velocity Component. In: Applied Nonlinear Analysis, Kluwer/Plenum, New York, 391-402.

[24] Scheffer, V. (1976) Partial Regularity of Solutions to the Navier-Stokes Equations. Pacific Journal of Mathematics, 66, 535-552.

https://doi.org/10.2140/pjm.1976.66.535

[25] Scheffer, V. (1980) The Navier-Stokes Equations on a Bounded Domain. Communications in Mathematical Physics, 73, 1-42.

https://doi.org/10.1007/BF01942692

[26] Federer, H. (1969) Geometric Measure Theory. Springer-Verlag, New York.

[27] Sohr, H. and Wahl, W. (1986) On the Regularity of the Pressure of Weak Solutions of Navier-Stokes Equations. Archiv der Mathematik, 46, 428-439.

https://doi.org/10.1007/BF01210782

[28] Lions, J.-L. and Peetre, J. (1964) Sur une classe d’espaces d’interpolation. Publications mathématiques de l’IHéS, 19, 5-68.

https://doi.org/10.1007/BF02684796

[29] Lions, J.-L. and Magenes, E. (1968) Problèmes aux limites non homogènes et applications. Vol. 1, Dunod, Paris.

[30] Adams, R.A. and Fournier, J.J.F. (2003) Sobolev Spaces. Journal of Pure and Applied Mathematics, 140, 305.

[31] Simon, J. (1990) Non-Homogeneous Viscous Incompressible Fluids: Existence of Velocity, Density and Pressure. SIAM Journal on Mathematical Analysis, 21, 1093-1117.

https://doi.org/10.1137/0521061

[32] Galperin, Y. (2014) Young’s Convolution Inequalities for Weighted Mixed (Quasi-) Norm Spaces. Journal of Inequalities and Special Functions, 5, 1-12.