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  .
We will be concerned with the 3D Navier-Stokes equations completed with initial and Dirichlet boundary conditions in bounded domains (as usual, Ω is the spatial domain, a regular, bounded and connected open set in “filled” by the fluid particles; is the time observation interval).
The key concept of suitable weak solution was introduced in  . 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  .
The estimate in  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  , 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  , an interesting criterion appears: if we normalize the solution and the sum of the associated kinetic and viscous energies and the Lp 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  .
・ In Wolf  , 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  , an improved version of the Caffarelli-Kohn-Nirenberg criterion is furnished, using ideas from  .
Note that, in a recent paper, Buckmaster and Vicol  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     . 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 . 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 is suitable.
Our results can be compared to other previous proofs of existence: the one in the Appendix in  (based on the construction of a family of time delayed linear approximations), the main result in Da Veiga  (relying on regularization with vanishing fourth-order terms), the main result in Guermond  (where Faedo-Galerkin techniques are employed) and, also, the results by Berselli and Spirito   , 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  and  . 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 and the usual L2 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  . In this reference, the best results known to date in relation to the regularity of the Navier-Stokes equations are established.
Let be a nonempty, regular, bounded and connected open set and assume that . Let us set and . We will consider local and global solutions to the Navier-Stokes equations in three dimensions
where verifies with and .
At a local level, we will consider solutions in sets of the form , where is open and connected and :
Definition 2.1 Let the open set be given. It will be said that the couple is a weak solution to the Navier-Stokes equations, Equation (1) in D if the following holds:
・ and .
・ 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
Let us assume that and let us consider the initial-boundary value problem
Definition 2.2 It will be said that the couple is a (global) weak solution to Equation (2) if the following holds:
・ and .
・ u and p satisfy the Navier-Stokes equations in Equation (2) in the distributional sense in Q
・ a.e. in Ω.
It is known that any couple satisfying the previous first and second points also verifies
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 is a weak solution to the Navier-Stokes equations in , then for each the functions
solve a similar problem in , with force . 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 , where λ is a characteristic length. We can affirm that
・ has dimension 1 and t is of dimension 2,
・ has dimension −1 and p has dimension −2,
・ f has dimension −3,
・ has dimension −1 and 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  and  . 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 be a weak solution to Equation (1) in and let be given. It will be said that is a singular point if u is not in any neighborhood of , that is, there are no r and C such that for a.e. in . The remaining points, those where u is locally bounded, will be called regular points.
According to a result by Serrin  , it is known that, if is a weak solution to Equation (2) and is a regular point for u, then u coincides a.e. with a function in a neighborhood of . 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 near , one just needs an estimate of the kind in time and in space, with sufficiently large r and s. Note that, in  , it is shown that a weaker condition is sufficient for regularity. Note also that, in accordance with the results in  , 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    . There, some estimates of the size of the set were given in terms of appropriate Hausdorff measures. Actually, the main result in  is the following:
Theorem 2.4 Assume that . There exists a weak solution to Equation (2) whose associated singular set S satisfies:
and uniformly in t.
Here, denotes the usual Hausdorff k-dimensional measure in .
This result was improved by Caffarelli, Kohn and Nirenberg in  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 be a cylinder in . It is said that 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 with ,
Then, the authors of  introduced the so called “parabolic” Hausdorff measure , as follows:
・ First, for any small and any , they set
Here, each is a parabolic cylinder, that is, a set of the form
for some .
・ Then, for any , they set
With the help of , a local partial regularity result can be established for any suitable solution:
Theorem 2.6 Let be a suitable solution to Equation (1) in D. Then the associated singular set satisfies .
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
for some .
In the sequel, for any and , we will denote by the following parabolic cylinder:
For the proof of Theorem 2.6, we need two results. The first one is the following:
Proposition 2.7 Suppose that is a suitable weak solution to Equation (1) in and with . There exist , and such that, if
In particular, 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
taking into account the dimensions of these quantities, we can easily deduce the following:
Corollary 2.8 Suppose that is a suitable solution to Equation (1) in the cylinder and , with . Then, if and , one has
and, consequently, every point in is regular.
Let us now set
The second fundamental result used in the proof of Theorem 2.6 is the following:
Proposition 2.9 Let be a suitable solution to Equation (1) in a neighborhood of . There exists such that, if
then 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  . 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 , S can be covered by a family of parabolic cylinders such that , the are mutually disjoint and
for all i. Using Hölder’s inequality, we deduce that
Taking , we find that , whence we see in particular that the Hausdorff dimension of S is at most 5/3.
To show that by a similar method, instead of the integral of , we need a global quantity of dimension 1. This is furnished by Proposition 2.9. Indeed, this result allows to replace Equation (6) by
and, this way, we are led to the estimate
whence we conclude that .
It is natural to ask if we can get a better estimate of the dimension of S. In other words, can we find such that ? 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 with 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 of the set S, instead of the standard measure . Both definitions are special cases of a construction made by Carathéodory that is detailed in  .
The argument used by Caffarelli, Kohn and Nirenberg is valid for any suitable solution. In the Appendix of  , they prove the existence of such a solution. Thus, the following holds:
Theorem 2.10 Suppose that , with and in Q. Then, there exists at least one suitable weak solution to the Navier-Stokes equations in Q satisfying
weakly in H as .
In addition, one has:
for all functions with and near .
2.3. On the Existence of Suitable Weak Solutions
The existence of a suitable weak solution to Equation (2) is established in  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  , using Faedo-Galerkin approximations. In both cases, the main difficult point is passing to the limit in the term in the right-hand side of the inequality. This requires nontrivial estimates on the pressure. In particular, Guermond  is able to achieve by reproducing for the discrete pressure some a priori estimates similar to the estimates of Sohr and Von Wahl in  .
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 , the instants and the approximations
, and , with and
First of all, let us check that the are well defined:
Lemma 3.1 The Euler scheme in Equation (9) is well defined. In other words, for every , there exists a unique solution to Equation (9).
The proof is immediate by induction. We only need to note that for each
Find such that , Now, let us see that the are uniformly bounded in the L2 norm. We have:
which can be rewritten in the form
Using the Cauchy-Schwarz and Young inequalities, we easily get that
for all n and, certainly, 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 and as follows:
・ is the unique continuous piecewise linear function satisfying
・ is the piecewise constant function characterized by
in for .
In a similar way, we can introduce the approximate pressures and forces (again piecewise constant). The following holds:
Lemma 3.2 For any N and almost every , one has
We can now present the main result of this section. It is related to the convergence of and towards a suitable weak solution to the Navier-Stokes equation:
Theorem 3.3 After eventual extraction of a subsequence, the functions converge weakly in , weakly-* in , strongly in and a.e. in Q towards a suitable weak solution to Equation (2) as .
For the proof of Theorem 3.3, it will be convenient to recall the following well known lemma (for instance, see the proof in  ):
Lemma 3.4 Let u be a function satisfying and . Then, u is a.e. equal to a continuous function from into H. In addition, the function is absolutely continuous and
a.e. in ,
where denotes the duality product in .
Proof of Theorem 3.3:
Let us first try to find the spaces where the functions , and 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 and let us carry out summation in m, from 0 to n. The following is obtained:
Obviously, this can also be written in the form
and, from the Cauchy-Schwarz and Young inequalities, we easily see that
This means that
is uniformly bounded in and . (16)
On the other hand, it can also be deduced from Equation (9) that
To estimate , we use its definition and the fact that, for any ,
Accordingly, we also have that
is uniformly bounded in and . (17)
Now, from classical interpolation results, we deduce that and is uniformly bounded in
It is well known that the estimates of Equation (17) and Equation (18) allow us to prove that the belong to and are uniformly bounded in the Sobolev spaces of fractional order for ; see for example  . Therefore, as a consequence of Aubin-Lions’ Theorem, the belong to a compact set of .
As a consequence, at least for a subsequence (again indexed by N), we must have:
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
strongly in for all , . (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  , play the role of the Sohr and Wahl’s estimates in  .
For , the space can be defined by the method of real interpolation between and , i.e. the so-called K-method of Lions and Peetre  ; see also  and  . We will denote by the closure of in . For any , the space and the corresponding norm are defined by duality and, in particular,
We will look for a uniform estimate of in a space of the form . This way, by applying De-Rham’s Lemma (see  ), we will get a bound of in and we will be able to take limits in the generalized energy inequality.
Note that, for all m, one has , where the and the are respectively given by
where and A is the Stokes operator. Recall that , with
(here, is the orthogonal projector). Also, recall that there exists an orthogonal basis of V formed by eigenfunctions ,
for all .
In the sequel, we will consider the functions and , whose definitions can be obtained from the and the in a way similar to and .
First, note that
and are uniformly bounded in . (24)
Let us now see what can be said of and . For all , we have
Let be such that with . Then
where the and the are given by
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 , and . Then, if and
for all .
The proof of this result can be found in  . Using Lemma 2.5 with , and , we find that
From the estimates in Equation (16) already obtained for
, it is immediate that, for any
is uniformly bounded in
and, consequently, also in
. Thus, choosing
On the other hand, for any smooth z, one has
where . Therefore, recalling the definition of the , we deduce that
Note that this estimate is valid for all , with and , .
Obviously, the same estimate is valid for . This proves that
In view of Equation (24) and Equation (27) and recalling that , it follows that and are also uniformly bounded in .
Now, from De-Rham’s Lemma  , we see that is uniformly bounded in , which is continuously embedded in . In particular, for , we have and we can take σ as close as desired to s, which gives 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
converges weakly in . (28)
Let us check that the local energy inequality holds for u and p.
If we multiply the Equation (14) by the function , where is non-negative and we integrate in space, we have:
If , there exists n such that and then, using Lemma 3.4, one has:
Note moreover that , because is by definition equal to .
• On the other hand,
Consequently, we see from Equation (29) that
If we integrate in time, we find that
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 converges strongly in and converges weakly in , we see that converges weakly in towards .
The final result is that
3.2. The Convergence of the Fully Discretized Problems
In this section, we will argue as in  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, and the . We will also consider two families of finite dimensional spaces and with the and the such that
and the are uniformly compatible, in the sense that there exists a constant independent of h such that the following inf-sup conditions are satisfied:
Now, we consider the approximations
with (the orthogonal projection of on ) and
for . The following result, which is a consequence of Equation (32), gives coherence to our scheme:
As before, the and serve to construct approximate solutions to the Navier-Stokes system. Thus, we define the functions etc. similarly to etc. The main result of this section is the following:
Theorem 3.7 After eventual extraction of a subsequence, the functions converge weakly in , weakly-* in , strongly in and a.e. in Q towards a suitable weak solution to Equation (2) as and .
Sketch of the proof:
Arguing as in the proof of the Theorem 3.3, it can be seen that the and the are uniformly bounded in and and, furthermore, . As in  , we can also prove that the are uniformly bounded in for any . Consequently, at least for a subsequence (still indexed with N and h), one has:
strongly in for all , . (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  . Here, we need the spaces for and their dual spaces .
The following estimates are established in  :
・ For any and any , one has
・ For any and any, one also has
As a consequence, it can be assumed that the converge weakly (for instance) in for all and the converge strongly in for all and . This is sufficient to ensure that converges weakly in towards .
Hence, arguing as in the final part of the proof of Theorem 3.3, it is not difficult to check that the limit of the 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:
We assume here that
As in  , 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 , with 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
and prove that, at least for a subsequence, the corresponding and converge, in an appropiate sense, to a suitable weak solution .
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
However, this is not clear at present. Note that the “reasonable” definition of a suitable weak solution should involve the following property: for any with ,
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:
・ Gear scheme:
θ-scheme: For α and β such that and , we compute , then and finally as follows:
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.
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).
 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.
 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.
 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.
 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.
 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.
 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.
 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.
 Lin, F.H. (1998) A New Proof of the Caffarelli-Kohn-Nirenberg Theorem. Communications on Pure and Applied Mathematics, 51, 241-257.
 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.
 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.