On the density patch problem for the 2-D inhomogeneous Navier-Stokes equations
Abstract.
In this paper, we first construct a class of global strong solutions for the 2-D inhomogeneous Navier-Stokes equations under very general assumption that the initial density is only bounded and the initial velocity is in . With suitable assumptions on the initial density, which includes the case of density patch and vacuum bubbles, we prove that Lions’ s weak solution is the same as the strong solution with the same initial data. In particular, this gives a complete resolution of the density patch problem proposed by Lions [25]: for the density patch data with a smooth bounded domain , the regularity of is preserved by the time evolution of Lions’s weak solution.
1. Introduction
We consider the inhomogeneous incompressible Navier-Stokes equations in :
(1.1) |
where stand for the density and velocity of the fluid respectively, and is a scalar pressure function.
Ladyženskaja and Solonnikov [21] first addressed the question of unique solvability of (1.1) in a bounded domain . Under the assumption that is divergence free and vanishes on and is bounded away from zero, they proved the global well-posedness in dimension , and local well-posedness in dimension . We refer to [1, 2, 3, 4, 9, 11, 18, 27] and references therein for the well-posedness results in the critical functional framework. In all these works, the density has to be continuous, bounded and bounded away from zero.
In the case when with a positive lower bound and , Kazhikov [20] proved the global existence of weak solution to the system (1.1), and Simon [30] removed the lower bound assumption on . While with , Danchin and Mucha [12] proved that the system (1.1) has a unique local in time solution. Paicu, Zhang and the fourth author [28] proved the global well-posedness of (1.1) with bounded away from zero and for any , and small in dimension 3. Furthermore, Zhang [31] established the global existence of strong solution with bounded away from zero and small in the critical Besov space . See [14, 10] for recent important progress in this direction.
1.1. Lions’s density patch problem
In this paper, we only focus on the case of , since the global regularity is still a famous open problem even for the 3-D incompressible Navier-Stokes equations.
Given , and satisfying , Lions [25] proved that the system (1.1) has a global weak solution satisfying the energy inequality
(1.2) |
Moreover, for any , the Lebesgue measure is independent of . Furthermore, Lions proposed the following open question(see page 34 in [25]):
We would like to mention another interesting open question: suppose that for a smooth domain , i.e., a patch of a homogeneous incompressible fluid “surrounded” by the vacuum(or a bubble of vacuum embedded in the fluid). Then, Theorem 2.1 yields at least one global weak solution and (2.17) implies that, for all for some set such that . In this case, (2.1)-(2.2) can be reformulated as a somewhat complicated free boundary problem. It is also very natural to ask whether the regularity of is preserved by the time evolution.
The density patch type problem was first solved by Liao and Zhang [23, 24] for the initial density of the form
where is any pair of positive constants and is a bounded, simply connected domain with boundary, see [15, 16] for related works. When , the first breakthrough is due to Danchin and Mucha [13], where they constructed a class of strong solutions for and without extra compatibility condition in the case when the fluid domain is either a bounded domain or the Torus; in particular, for the density patch data, the regularity of is preserved by the time evolution of strong solution they constructed. Prange and Tan [29] extend some results in [13] to the case of ; moreover, for the density with a bubble of vacuum(i.e., ), they can prove the weak-strong uniqueness between Lions’s weak solution and strong solution they constructed, thus the regularity of is preserved by the time evolution of Lions’s weak solution in such case. However, for the density patch type data(i.e., ), [29] still need to impose the compatibility condition (1.3) on the initial data and a weighted bound for the weak solution. Thus, the following original version of Lions’ s density patch problem remains open:
For the density patch data with a smooth bounded domain , whether the regularity of is preserved by the time evolution of Lions’s weak solution?
1.2. Main results
The goal of this paper is twofolds. First of all, we prove the global existence of strong solution under very general assumption on the initial density.
Theorem 1.1.
Given the initial data satisfying , the system (1.1) has a global solution satisfying the weak energy inequality (1.2), , and the following properties:
-
•
;
-
•
, here ;
-
•
for all .
Remark 1.1.
-
(1)
Compared with the existence part in [29], here we only require and 111It’s natural to assume , otherwise the problem is trivial and the solution is , .. For the case of density patch data(FFV-2 condition in Theorem A, [29]), except for 222Indeed, they assume the initial velocity satisfies and ., [29] also imposed the following key compatibility condition
(1.3) - (2)
- (3)
-
(4)
As for any , there exists a unique continuous flow associated with globally defined on by
(1.4) Then .
To solve the density patch problem, we need to explore the further regularity of strong solution under one of the following two assumptions:
-
(H1)
has a compact support, thus .
-
(H2)
there exist such that for all .
Remark 1.2.
Indeed, (H1) corresponds to the density patch problem, and (H2) allows vacuum bubbles, especially countable equal size vacuum bubbles. While, the following condition in [29]
only allows finitely many vacuum bubbles with an equal non-zero measure.
Theorem 1.2.
Under the assumptions of Theorem 1.1, if satisfies (H1) or (H2), then for all , the solution satisfies .
As a corollary, we can prove that the regularity of is preserved by the time evolution of strong solution with the initial density .
Corollary 1.1.
Remark 1.3.
The second goal of this paper gives a complete resolution of Lions’s density patch problem via proving the weak-strong uniqueness and Corollary 1.1.
Theorem 1.3.
Remark 1.4.
At this stage, we only prove the uniqueness of weak solutions for satisfying (H1) or (H2). This gives a partial answer to Lions’s open question about the uniqueness of weak solutions in two dimensions (see page 31 in [25]).
1.3. Sketch of the proof
First of all, it is relatively standard to obtain the following energy estimates
and for ,
Here and in what follows, the constant depends only on and .
Based on an important observation due to an algebraic fact if for a matrix , we can further show that for (see Lemma 3.2)
However, due to the presence of vacuum of the density, it is highly nontrivial to obtain and from the above energy estimates, which are crucial to prove the existence of strong solution and weak-strong uniqueness.
To prove , we use the following propagation property of lower bound of the density in the sense that if , then for ,
Let’s emphasize that this property also holds for weak solutions(see Lemma 2.2). This lower bound ensures that there exists a such that . On the other hand, from the energy estimates, we can deduce that
Thus, we can conclude the following bound
This bound is completely new and may be of independent interest. See Lemma 3.3 for the detail.
To prove , the key ingredient is to establish the following pointwise inequality of , which is also completely new:
where we further have
See Lemma 3.5 for the detail.
To prove the regularity propagation of the density patch or vacuum bubbles, we need to improve the regularity of . Indeed, if satisfies (H1) or (H2), then for all , we have . For this, the key point is to use the following important inequality(see Lemma 2.4)
such that
For the weak-strong uniqueness, we follow the argument of Prange and Tan [29]. Let be the strong solution obtained in Theorem 1.1 and be the weak solution provided by Lions in [25] with same initial data . We denote . First of all, we have
where Using a duality argument, we can show that
where is the solution of the linear backward parabolic system (4.16). By Lemma 2.4, we can prove that
Compared with [29], the main innovation is to control the trouble term , which is based on the trilinear estimates in Lemma 4.1 and Lemma 4.2, in a way as
Thus, we can conclude
Then Osgood’s lemma ensures that on .
Notations.
-
(1)
is the material derivative of the velocity.
-
(2)
.
-
(3)
.
-
(4)
For matrices and , we denote .
-
(5)
is the Lebesgue measure of for .
-
(6)
We shall always denote to be an absolute constant which may vary from line to line. The dependence of the constant will be explicitly indicated if there are any exceptions.
2. Some important properties of weak solutions
In this section, we prove some important properties of weak solutions only using the density equation and the energy inequality (1.2). These properties are also very useful for strong solutions.
Lemma 2.1.
For any , we have
Proof.
As , there exist and such that
(2.1) |
Lemma 2.2.
Assume (2.1) and . Then we have
Proof.
Fix and . For any , we choose a bump function such that , and . By Lemma 2.1, we have
It follows from the support properties of that
hence
Letting completes the proof. ∎
Lemma 2.3.
If , then there holds
Proof.
Let be such that for and . Then for . For any , we define for , then , and , thus . By Lemma 2.1, we have
Therefore, letting and using for completes the proof.∎
Lemma 2.4.
Let and . If satisfies (H1), then there exists a constant depending only on , such that (for )
(2.2) |
If satisfies (H2), then there exists a constant depending only on , , , , , such that (for )
(2.3) |
3. Global existence and regularity of strong solution
This section is devoted to the global existence of strong solution and the regularity propagation of density patch and vacuum bubbles.
3.1. A priori energy estimates
Let be a smooth solution to the system (1.1) on satisfying .
First of all, (1.2) and Lemma 2.1–Lemma 2.4 hold for strong solutions. The following result has been proved in [26], Lemma 3.2.
Lemma 3.1.
There exists a positive constant depending only on and such that for , there holds
Motivated by Lemma 3.3 in [26], we derive the estimates for the material derivative of the velocity, which are crucial for the higher order estimates of the velocity. Here we have an important observation due to an algebraic fact if for a matrix . Then we can remove the assumption in [26].
Lemma 3.2.
There exists a positive constant depending only on , and such that for ,
Proof.
Applying the material derivative to the momentum equation of (1.1), we get 333Repeated indices represent summation.
Taking the inner product with to the above equation, we get by integration by parts and that
(3.1) |
where
The first two terms on the right hand side can be bounded by
(3.2) |
As for , using integration by parts, we get (as , )
We also have
and
Notice that . Then by symmetry, we obtain
(3.3) |
Moreover, it follows from that
(3.4) |
Here denotes the Hardy space.
3.2. bounds of and
In order to prove the existence of strong solution, we need to obtain a bound of .
Lemma 3.3.
Assume (2.1). Then there exists a constant depending only on , , , and such that
Proof.
By Lemma 3.19 and Lemma 3.2, we have
which ensure that
(3.9) |
Thus, for any , we have (as )
(3.10) |
where
which along with (3.9) and (3.10) ensures that
(3.11) |
Therefore, by (3.11), , and , we arrive at
This completes the proof of the lemma. ∎
For the estimate of , we need the following lemma.
Lemma 3.4.
Let and be smooth functions satisfying 444 , .
then there exists a constant independent of such that
Proof.
By the translation and dilation, it suffices to show that
(3.12) |
for some constant independent of and . We split the proof into three cases.
Case 1. . Let be such that and , and let be given by for all . Thanks to
we have
where .
Using the integration by parts and support properties of , we get
and
Combining above estimates, we obtain (3.12).
Case 2. . Notice that . For , we let
Firstly, we have
(3.13) | ||||
Since for , then
so we get
(3.14) |
and similarly
(3.15) |
Note that , then
(3.16) |
On the other hand, using integration by parts, we have
(3.17) |
where we have used the fact that
Combining (3.13), (3.14), (3.15), (3.16), (3.17) and letting , we obtain (3.12).
Case 3. . Notice that . The proof is similar to Case 2. ∎
Now we are in a position to give a priori estimate of .
Lemma 3.5.
Assume that and is a smooth solution of the following equation
(3.18) |
Then there is a constant depending only on , and such that
(3.19) | ||||
(3.20) |
Proof.
Let , then . By (3.18), we have , and thus . Note that
where , and . Hence, by the fact and , we obtain
As for , then by Lemma 3.4, we get the desired estimate (3.19).
Next we prove (3.20). We claim that there exists a constant such that
(3.21) |
for any . Assume (3.21) for the moment. Let . Then we infer that
Now it follows from (1.2), Lemma 3.19, Lemma 3.2 and Hölder’s inequality that
This proves (3.20).
Finally, it remains to prove (3.21). Fix and . We emphasize that in the proof of (3.21), the constant is independent of and . Let . By Poincaré inequality we get that
(3.22) |
As a consequence, for all , if we let be such that , then
(3.23) | ||||
On the other hand, as , we get by Hölder inequality and interpolation inequality in the Lorentz space that
and by (3.22), we have
hence
(3.24) | ||||
Thanks to (3.23) and (3.24), we have
(3.25) |
3.3. Proof of global existence
In this subsection, we prove Theorem 1.1. We follow the argument in [13]. The idea is to take advantage of classical result to construct smooth solutions corresponding to smoothed-out approximate data with no vacuum, then to pass to the limit. More precisely, we consider
such that
As , there exist and such that . Then for small enough, thus (2.1) holds for and replaced by . We also have for small enough, and . In light of the classical strong solution theory for the system (1.1), there exists a unique global smooth solution corresponding to data . Thus, the triple satisfies all the a priori estimates of previous subsections uniformly with respect to . In particular, by Lemma 3.19, Lemma 3.2 and Lemma 3.3, we have , and
By the definition, we have
for any . Thus, we obtain
Similarly, we have
Let then for small enough. By Lemma 2.2 we have for , . Then by Lemma A.2 we have for , . Thus for , which also holds for as the left hand side is increasing in .
By Lemma 3.4 in [13], we know that for all
(3.26) |
where depending only on and . By (1.2) and Lemma 3.19, we also have
By the interpolation with (3.26), we can deduce that
This implies that up to a subsequence, in . By Theorem 2.5 in [25], one can show that
Thus, we can pass to the limit for all the terms in the definition of weak solutions, and the weak energy inequality can be obtained by Fatou property. Moreover Lemma 3.19, Lemma 3.2, (3.20) and Lemma 3.3 also hold for the limit solution . Furthermore, the moment equation is fulfilled in the strong sense, i.e.,
for some pressure function with the estimate in Lemma 3.19.
3.4. Regularity propagation of density patch and vacuum bubbles
Proof of Theorem 1.2.
Now we prove Corollary 1.1.
Proof.
Assume that corresponds to the level set of some function with regularity and on . Then we have with , and on .
Next we check that satisfies (H1) or (H2). In fact, if , then , which implies that satisfies (H1) as is bounded. If , then , and thus there exist such that for all we have
which means that satisfies (H2).
4. Weak-strong uniqueness
Throughout this section, let be the strong solution obtained in Theorem 1.1 and be the weak solution provided by Lions in [25] with the same initial data . Furthermore, we assume that satisfies (H1) or (H2).
4.1. Trilinear estimates
Let be fixed such that , and we denote for .
Lemma 4.1.
There exists such that for any with , and any we have
(4.1) | ||||
(4.2) |
Proof.
Lemma 4.2.
For any , there exists such that for any and with , for any we have
(4.4) | ||||
(4.5) |
Proof.
Let be fixed such that and , and let for and . For , we denote
Then Poincaré’s inequality implies that
(4.6) |
for some absolute constant independent of and . We claim that
(4.7) | ||||
(4.8) | ||||
(4.9) |
for some constant depending only on .
Now we prove (4.4) by admitting (4.7), (4.8) and (4.9), whose proof will be given as soon as we finish the proof of (4.4). We write . First of all, by (4.3) we have
(4.10) |
Proof of (4.8).
By we have . Thus . Since , we have . If then and there exists such that for . Thus , which contradicts . Thus and .∎
4.2. Proof of weak-strong uniqueness
Proof of Theorem 1.3.
Here we follow the proof from Prange and Tan ([29], Proposition 4.1). We focus on the case when is small. We denote and . Then the system for () reads (here )
(4.13) |
We set so that . Testing the first equation of (4.13) against and using the fact , following the arguments in [29], for all we have
(4.14) |
where
(4.15) |
Next we use a duality argument to control the difference in . Let be the solution of the following linear backward parabolic system
(4.16) |
Then there holds
(4.17) | ||||
Testing the equation of in the system (4.13) by yields that
(4.18) |
Firstly, we have (by Lemma 2.4 and (LABEL:v_estimate))
where depends only on and . Hence, by adjusting to smaller if necessary, it follows from (4.18) that
(4.19) |
Using (4.14), by adjusting to smaller if necessary, we have for . Taking in Lemma 4.1 if satisfies (H2) (or in Lemma 4.2 if satisfies (H1)), and using Lemma 2.4 and (4.15), we obtain
(4.20) |
Thus, it follows from (4.14), (LABEL:v_estimate), (4.19) and (4.20) that
Let , then by Lemma 3.2 we have
hence , , and we also have
Furthermore, the function is increasing near and satisfies . Hence we can apply Osgood’s lemma ([5], Lemma 3.4) to conclude that on . By (4.19), (4.20) and (4.15), we have
Now, a direct estimate of (4.13) gives , and by for all (see Lemma 2.4), we know that a.e. in .
The uniqueness on the whole time interval can be obtained by a bootstrap argument. ∎
Appendix A Technical lemmas
The following lemma comes from Theorem B.1 in [25].
Lemma A.1.
For and , there exists a positive constant depending only on and such that for all ,
(A.1) |
The following technical lemma is similar to Proposition A.4 in [29], which will be used constantly.
Lemma A.2.
Assume that and . Then for all , there exists a constant depending only on , and such that for ,
Proof.
This completes the proof. ∎
Lemma A.3.
Let be as in Subsection 4.1. There exists a constant such that for all and all , there holds
(A.3) |
where . In particular, for all , all satisfying , we have
(A.4) |
Proof.
Lemma A.4.
There exists a contant such that for all , there holds
(A.5) |
Acknowledgments.
D. Wei is partially supported by the National Key R&D Program of China under the grant 2021YFA1001500. Z. Zhang is partially supported by NSF of China under Grant 12288101.
References
- [1] H. Abidi, Équation de Navier-Stokes avec densité et viscosité variables dans l’espace critique, Rev. Mat. Iberoam., 23 (2007), 537-586.
- [2] H. Abidi and G. Gui, Global well-posedness for the 2-D inhomogeneous incompressible Navier-Stokes system with large initial data in critical spaces, Arch. Ration. Mech. Anal., 242 (2021), 1533–1570.
- [3] H. Abidi, G. Gui and P. Zhang, On the wellposedness of -D inhomogeneous Navier-Stokes equations in the critical spaces, Arch. Ration. Mech. Anal., 204 (2012), 189-230.
- [4] H. Abidi and M. Paicu, Existence globale pour un fluide inhomogéne, Ann. Inst. Fourier (Grenoble), 57 (2007), 883-917.
- [5] H. Bahouri, J. Y. Chemin and R. Danchin, Fourier Analysis and Nonlinear Partial Differential Equations, Grundlehren der Mathematischen Wissenschaften, vol. 343 (2011), Springer-Verlag, Berlin, Heidelberg.
- [6] A. L. Bertozzi and P. Constantin, Global regularity for vortex patches, Comm. Math. Phys., 152 (1993), 19–28.
- [7] J. Y. Chemin, Persistance de structures géométriques dans les fluides incompressibles bidimensionnels, Ann. Sci. École Norm. Sup. (4), 26 (1993), 517–542.
- [8] H. Choe and H. Kim, Strong solutions of the Navier–Stokes equations for nonhomogeneous incompressible fluids, Commun. PDE, 28(2003), 1183–201
- [9] R. Danchin, Density-dependent incompressible viscous fluids in critical spaces, Proc. Roy. Soc. Edinburgh Sect. A, 133 (2003), 1311-1334.
- [10] R. Danchin, Global well-posedness for 2D inhomogeneous viscous flows with rough data via dynamic interpolation, arXiv:2404.02541.
- [11] R. Danchin and P. B. Mucha, A Lagrangian approach for the incompressible Navier-Stokes equations with variable density, Comm. Pure Appl. Math., 65 (2012), 1458-1480.
- [12] R. Danchin and P. B. Mucha, Incompressible flows with piecewise constant density, Arch. Ration. Mech. Anal., 207 (2013), 991-1023.
- [13] R. Danchin and P. B. Mucha, The incompressible Navier-Stokes equations in vacuum, Comm. Pure Appl. Math., 72 (2019), 1351-1385.
- [14] R. Danchin and S. Wang, Global unique solutions for the inhomogeneous Navier-Stokes equations with only bounded density in critical regularity spaces, Comm. Math. Phys., 399 (2023), 1647–1688.
- [15] R. Danchin and X. Zhang, On the persistence of Hölder regular patches of density for the inhomogeneous Navier-Stokes equations, J. Éc. polytech. Math., 4 (2017), 781–811
- [16] F. Gancedo and E. García-Juárez, Global regularity of 2D density patches for inhomogeneous Navier-Stokes, Arch. Ration. Mech. Anal. 229 (2018), 339–360.
- [17] X. Huang, J. Li and Z. Xin, Global well-posedness of classical solutions with large oscillations and vacuum to the three-dimensional isentropic compressible Navier-Stokes equations, Comm. Pure Appl. Math., 65(2012), 549–585.
- [18] J. Huang, M. Paicu and P. Zhang, Global wellposedness to incompressible inhomogeneous fluid system with bounded density and non-Lipschitz velocity, Arch. Ration. Mech. Anal., 209 (2013), 631–682.
- [19] X. Huang and Y. Wang, Global strong solution with vacuum to the two-dimensional density-dependent Navier–Stokes system, SIAM J. Math. Appl., 46(2014), 1771–88
- [20] A. Kazhikhov, Solvability of the initial-boundary value problem for the equations of the motion of an inhomogeneous viscous incompressible fluid, Dokl. Akad. Nauk SSSR, 216 (1974), 1008-1010.
- [21] O. A. Ladyženskaja and V.A. Solonnikov, The unique solvability of an initial-boundary value problem for viscous incompressible inhomogeneous fluids. (Russian) Boundary value problems of mathematical physics, and related questions of the theory of functions, 8, Zap. Naučn. Sem. Leningrad. Otdel. Mat. Inst. Steklov., 52 (1975), 52-109, 218-219.
- [22] X. Liao and S. Zodji, Global-in-time well-posedness of the compressible Navier-Stokes equations with striated density, arXiv:2405.11900.
- [23] X. Liao and P. Zhang, On the global regularity of the two-dimensional density patch for inhomogeneous incompressible viscous flow, Arch. Ration. Mech. Anal., 220 (2016), 937–981.
- [24] X. Liao and P. Zhang, Global regularity of 2D density patches for viscous inhomogeneous incompressible flow with general density: low regularity case, Comm. Pure Appl. Math., 72(2019), 835–884.
- [25] P. L. Lions, Mathematical topics in fluid mechanics. Vol. 1. Incompressible models, Oxford Lecture Series in Mathmatics and its Applications, 3. Oxford Science Publications. The Clarendon Press, Oxford University Press, New York, 1996.
- [26] B. Lü, X. Shi and X. Zhong, Global existence and large time asymptotic behavior of strong solutions to the Cauchy problem of 2D density-dependent Navier-Stokes equations with vacuum, Nonlinearity, 31 (2018), no.6, 2617-2632.
- [27] M. Paicu and P. Zhang, Global solutions to the -D incompressible inhomogeneous Navier-Stokes system, J. Funct. Anal., 262 (2012), 3556-3584.
- [28] M. Paicu, P. Zhang and Z. Zhang, Global unique solvability of inhomogeneous Navier-Stokes equations with bounded density, Comm. Partial Differential Equations, 38 (2013), 1208-1234.
- [29] C. Prange and J. Tan, Free boundary regularity of vacuum states for incompressible viscous flows in unbounded domains, arXiv:2310.09288v3.
- [30] J. Simon, Nonhomogeneous viscous incompressible fluids: existence of velocity, density, and pressure, SIAM J. Math. Anal., 21 (1990), 1093-1117.
- [31] P. Zhang, Global Fujita-Kato solution of the -D inhomogeneous incompressible Navier-Stokes system, Adv. Math., 363 (2020), 107007, 43 pp.