arxiv: v1 [math.ap] 2 Jan 2013

Transcription

1 GLOBAL UNIQUE SOLVABILITY OF INHOMOGENEOUS NAVIER-STOKES EQUATIONS WITH BOUNDED DENSITY MARIUS PAICU, PING ZHANG, AND ZHIFEI ZHANG arxiv:131.16v1 [math.ap] Jan 13 Abstract. In this paper, we prove the global existence an uniqueness of solution to -imensional for =,3) incompressible inhomogeneous Navier-Stokes equations with initial ensity being boune from above an below by some positive constants, an with initial velocity u H s R ) for s > in -D, or u H 1 R 3 ) satisfying u L u L being sufficiently small in 3-D. This in particular improves the most recent well-poseness result in [1], which requires the initial velocity u H R ) for the local well-poseness result, an a smallness conition on the fluctuation of the initial ensity for the global well-poseness result. Keywors: Inhomogeneous Navier-Stokes equations, well-poseness, Lagrangian coorinates. AMS Subject Classification ): 35Q3, 76D5 1. Introuction In this paper, we consier the global existence an uniqueness of the solution to the following -imensional for =, 3) incompressible inhomogeneous Navier-Stokes equations with initial ensity in R ) an having a positive lower boun: 1.1) t ρ+u ρ =, t,x) R + R, ρ t u+u u) u+ p =, ivu =, ρ,u) t= = ρ,u ), where ρ,u stan for the ensity an velocity of the flui respectively, p is a scalar pressure function, an the viscosity coefficient is suppose to be 1. Such a system escribes a flui which is obtaine by mixing two miscible fluis that are incompressible an that have ifferent ensities. It may also escribe a flui containing a melte substance. One may check [17] for the etaile erivation of this system. Given ρ R ), an u satisfying ivu =, ρ u L R ), Lions [17] see also [5, 19] an the references therein for an overview of results on weak solutions of 1.1)) prove that 1.1) has a global weak solution so that 1 ρt)ut) L + Moreover, for any α an β, the Lebesgue measure uτ) L τ 1 ρ u L. µ { x R ; α ρt,x) β } is inepenent of t. In imension two an uner the aitional assumption that ρ is boune below by a positive constant an u L R ), smoother weak solutions may be built. Their existence stems from a quasi-conservation law involving the norm of u,t);l R )) an of t u, p an u L,T);L R )) for any T <. For both types of weak solutions however, the problem of uniqueness has not been solve. Date: Dec., 1. 1

2 M. PAICU, P. ZHANG, AND Z. ZHANG Layženskaja an Solonnikov [16] first aresse the question of unique solvability of 1.1). More precisely, they consiere the system 1.1) in a boune omain Ω with homogeneous Dirichlet bounary conition for u. Uner the assumption that u W p,p Ω) p > ) is ivergence free an vanishes on Ω an that ρ C 1 Ω) is boune away from zero, then they [16] prove Global well-poseness in imension = ; Local well-poseness in imension = 3. If in aition u is small in W p,p Ω), then global well-poseness hols true. More recently, Danchin [8] establishe the well-poseness of the system 1.1) in the whole space R for small perturbations of some constant ensity. Abii, Gui an Zhang [3] investigate the large time ecay an global stability to any global smooth solutions of 1.1). Another important feature of 1.1) is the scaling invariant property: if ρ, u) is a solution of 1.1) associate to the initial ata ρ,u ), then ρλ t,λx),λuλ t,λx)) is also a solution of 1.1) associate to the initial ata ρ λx),λu λx)). A functional space for the ata ρ,u ) or for the solution ρ,u) is sai to be at the scaling of the equation if its norm is invariant uner the above transformation. In this framework, it has been state in [1, 7] that for the initial ata ρ,u ) satisfying ρ 1) Ḃ p p,1 R ), u Ḃ an that for a small enough constant c we have for any p [1,) p 1 ρ + u p 1 Ḃ p,1 Ḃ p,1 R ) with ivu = p 1 p,1 c, p existence of global solution ρ,u, p) with ρ 1 C b [, ); Ḃ p Ḃ 1 p 1 p,1 R )), an t u, u, p L 1 R + ;Ḃ p,1 R )); uniqueness in the above space if in aition p. p,1 R )), u C b [, ); These results have been somewhat extene in [] so that u belongs to a larger Besov space. Paicu an Zhang [18] further extene the well-poseness result in [] so that even if one component of the initial velocity is large, 1.1) still has a unique global solution. The smallness assumption for the initial ensity in [1, 7] has also been remove in [4], an the restriction of p [1,] for uniqueness result in [1, 7] has been remove recently in [9]. A byprouct in [4] implies the global existence of solutions to 1.1) in 3-D with initial ensity in R 3 ) an having a positive lower boun, an initial velocity being sufficiently small in H R 3 ). The authors [9] prove the global wellposeness of 1.1) provie that ρ 1 MḂ p 1 p,1 ) + u Ḃ p forsomesufficientlysmallconstantc,wheremḃ 1 p,1 R p ))enotesthemultiplierspaceofḃ 1 p,1 R ). This space in particular inclues initial ensities having small jumps across a C 1 interface. Again in the scaling invariant framework, the authors [13] prove the global existence of weak solutions to 1.1) provie that the initial ata satisfy the nonlinear smallness conition: ρ u h ) Ḃ 1+ p exp Cr u ) r c µ p,r Ḃ 1+ p p,r for some positive constants c,c r an 1 < p <, 1 < r <, where u h = u1 ) an u = u h,u ). With a little bit more regularity assumption on the initial velocity, they [13] also prove the uniqueness of such solutions. p 1 p,1 c,,,u 1

3 GLOBAL UNIQUE SOLVABILITY OF INHOMOGENEOUS NS EQUATIONS 3 In general when ρ R ) with a positive lower boun an u H R ), Danchin an Mucha [1] prove that the system 1.1) has a unique local solution. Furthermore, with the initial ensity fluctuation being sufficiently small, for any initial velocity u B 1 4, R ) L R ) in two space imensions, an u B q q,p R ) with 1 < p <, < q < an p 1 q, they also prove the global well-poseness of 1.1). On the other han, Hoff [14, 15] prove the global existence of small energy solutions to the isentropic compressible Navier-Stokes system. The main iea in [14, 15] is that with appropriate time weight see Remark 1. for etails), one can close the energy estimate for space erivatives of the velocity fiel even if the initial velocity only belongs to L R ). Motivate by [14, 15], we shall investigate the global well-poseness of 1.1) with less regular initial velocity than that in [1] an without the small fluctuation assumption on the initial ensity. We emphasize that the Lagrangian iea introuce in [9, 1] will also be essential for the proof of the uniqueness result here. Our main results in this paper can be liste as follows. Theorem 1.1. Let s >. Given the initial ata ρ,u ) satisfying 1.) < c ρ x) C < +, u H s R ), the system 1.1) has a unique global solution ρ,u) such that 1.3) c ρt,x) C for t,x) [,+ ) R, A t) C u L, A 1 t) C u H s exp{ C u 4 L }, A t) C 1+ u 8 L ) u H s + u 4 H s)exp{ C u 4 L }, for any t [,+ ). Here C is a constant epening on c,c, an A t),a 1 t), an A t) are efine by A T) ef = 1 T sup ρ ut,x) x+ u xt, t [,T] R R A 1 T) ef = 1 T sup σt) 1 s ut,x) t [,T] R x+ σt) 1 s ρ u t + u + p ) xt, R A T) ef = 1 sup + t [,T] T with σt) ef = min1,t). σt) s R ρ ut t,x) + ut,x) + p ) x R σt) s u t xt, Theorem 1.. Given the initial ata ρ,u ) satisfying 1.4) < c ρ x) C < +, u H 1 R 3 ), there exists a constant ε > epening only on C such that if 1.5) u L u L ε, the system 1.1) has a unique global solution ρ, u) which satisfies c ρt,x) C for t,x) [,+ ) R 3, 1.6) B t) ρ 1 u L, B 1 t) u L, B t) C 1+ u 4 L ) u L exp { C u L + u L ) },

4 4 M. PAICU, P. ZHANG, AND Z. ZHANG for any t [,+ ). Here C is a constant epening on C, an B t),b 1 t), an B t) are efine by T B T) ef = sup ρ ut,x) x+ u xt, t [,T] R 3 R 3 T B 1 T) ef = sup ut,x) x+ ρ ut + u + p ) xt, t [,T] R 3 R 3 B T) ef = sup σt) ρ ut t,x) + ut,x) + p ) T x+ σt) u t xt. t [,T] R 3 R 3 Remark 1.1. We shoul point out that we o not nee the lower boun assumption for the initial ensity in the existence part of Theorem 1.. Inee, the constant C in 1.6) is inepenent of c in 1.4). We can also prove the local existence an uniqueness solution to the system 1.1) even if the initial velocity oes not satisfy the smallness conition 1.5). One may check Remark.1 for etails. Remark 1.. The powers to the weight σt) in the energy functionals, A i t),b i t) for i = 1,, are motivate by the following observation: let s be a negative real number an p,r) [1, ], a constant C exists such that 1.7) C 1 f Ḃs t e s t f p,r L p L r R + ; t t ) C f Ḃ, p,r s see Theorem.34 of [6] for instance). In particular, if u H s R ) for s,1), u H s 1 R ) B, s 1 R ) an u B, s R ). Then accoring to 1.7), t 1 s e t u L C u H s an t s e t u L C u H s. This in some sense explains the weights in 1.3). Similarly, we can explain the weights in 1.6). Remark 1.3. We shoul point out that we can not irectly apply Theorem 1 of [1] concerning the uniqueness of solutions to 1.1) to conclue the uniqueness part of Theorem 1.1 an Theorem 1.. Yet the Lagrangian iea in [9, 1] can be successfully applie to prove the uniqueness part of both Theorems 1.1 an 1.. An the uniqueness result of Germain [11] can not be applie here either. The uniqueness result of [11] requires the ensity function satisfying ρ [,T];L R )), but here our ensity function only belongs to [,T] R ). Moreover, the velocity fiel in Theorems 1.1 an 1. oes not satisfy the time growth conition in [11], especially in Theorem Global solutions to 1.1) with large boune ensity The purpose of this section is to present the proof to the existence part of both Theorem 1.1 an Theorem Existence of the solution in -D. Proof to the existence part of Theorem 1.1. Let j ǫ be the stanar Frierich s mollifier. We efine An we choose ǫ so small that ρ ǫ = j ǫ ρ, u ǫ = j ǫ u. c ρǫ x) C, x R. With the initial ata ρ ǫ,uǫ ), the system 1.1) in -D has a unique global smooth solution ρǫ,u ǫ ). In what follows, we shall only present uniform energy estimates 1.3) for the approximate solutions ρ ǫ,u ǫ ). Then the existence part of Theorem 1.1 essentially follows from 1.3) for ρ ǫ,u ǫ ) an a stanar compactness argument. The uniqueness part of Theorem 1.1 will be prove in Section 3.

5 GLOBAL UNIQUE SOLVABILITY OF INHOMOGENEOUS NS EQUATIONS 5 To simplify the notations, we will omit the superscript ǫ in what follows. First of all, applying the basic L energy estimate to 1.1) gives.1) A t) = 1 ρ1 u L C u L for t R +. While it follows from the transport equation of 1.1) an 1.) that.) c ρt,x) C for t,x) R + R. To erive the estimate for A 1 t), we get, by taking the L inner prouct of the momentum equation of 1.1) with u t, that ρ u t x+ 1 u x = ρu u) u t x, R t R R from which, we infer à 1 t) ef = στ)ρ u t xτ + 1 R σt) ut,x) x R.3) t στ)ρu u) u t xτ + u xτ. R R In this subsection, we shall frequently use the following version of Gagliaro-Nirenberg inequality:.4) a L p R ) C a p By virtue of.) an.4), we obtain.5) C C R στ)ρu u) u t xτ L R ) a 1 p L R ) στ) uτ) L 4 uτ) L 4 ρ 1 ut τ) L τ for p <. στ) uτ) 1 L uτ) L uτ) 1 L ρ 1 ut τ) L τ C δ στ) uτ) L uτ) L uτ) L τ +δ C δ A t) στ) uτ) 4 L τ +δ στ) ρ 1 ut τ) L τ στ) uτ) L τ +δã1t), for any δ >, where C δ is a positive constant so that C δ as δ. Whereas thanks to 1.1), we write.6) u+ p = ρu t +u u), which along with the classical estimate on the Stokes system ensures that u L + p L ) ρu t L + ρu u L so that ρu t L + u 1 L u L u 1 L ),.7) u L + p L C ρ 1 u t L + u L u L ). Substituting.7) into.5) gives rise to.8) στ)ρu u) u t xτ Cδ A t) R στ) uτ) 4 L τ +CδÃ1t).

6 6 M. PAICU, P. ZHANG, AND Z. ZHANG Summing up.1),.3) an.7)-.8) an taking δ sufficiently small, we obtain à 1 t)+ στ) u xτ C u R L uτ) L Ã1τ)τ + u L ). Applying Gronwall s inequality gives.9) à 1 t)+ στ) u xτ C u R L exp { C u 4 } L. To obtain the estimate of A 1 t), we nee to use an interpolation argument. For this, we consier the linear momentum equation ρ t v +u v) v + p =, v,x) = v x). Then it follows from the same line to the proof of.9) an.7) that ρ v t + v + p )xτ + vt,x) x C v R R H exp { C u 1 4 } L, στ)ρ v t + v + p )xτ +σt) vt,x) x C v R R L exp { C u 4 } L. We efine the linear operator Tv = v. The above inequalities tell us that Tv L Cexp { C u 4 L } v H 1, Tv L Cσt) 1 exp { C u 4 L } v L, from which an Riesz-Thorin interpolation theorem [1], we infer v L = Tv L Cσt) 1+s exp { C u 4 L } v H s. Further, we efine a family of operators T z v = σt) z ρ 1 t vt,x) for Rez [,1]. Then we have T iy v L,t;L ) Cexp { C u 4 L } v H 1, T 1+iy v L,t;L ) Cexp { C u 4 L } v L, for any y R. Apply Stein interpolation theorem [1] to get στ) 1 s ρv 1 t L L,t) = T 1 sv L,t;L ) Cexp { C u 4 } L v H s. The other terms can be treate in a similar way. Therefore we have.1) A 1 t) C u H s exp{ C u 4 L }. Finally, we manipulate the H energy estimate for u. We first get, by taking the time erivative to the momentum equation of 1.1), that ρu tt +u u t ) u t + p t = ρ t u t +u u) ρu t u. Taking the L inner prouct of the above equation with u t, an then using integration by parts, we write 1 ρ u t x+ u t x t R R = ρ t u t x ρ t u u) u t x ρu t u) u t x, R R R

7 from which, we infer GLOBAL UNIQUE SOLVABILITY OF INHOMOGENEOUS NS EQUATIONS 7.11) Ã t) ef = 1 σt) s R ρ u t x+ = στ) s R ρ t u t xτ στ) s R u t xτ στ) s R ρu t u) u t xτ + s) ef =A+B +C+D. στ) s R ρ t u u) u t xτ στ) 1 s R ρ u t xτ It is obvious to check that.1) D s)a 1 t), an.13) C στ) s ρ u L u t L 4 τ στ) s u L u t L u t L τ στ) s u L ρ 1 ut L τ + 1 4Ãt). While noticing that ρ t = u ρ, we get by using integration by parts an.4) that.14) A = στ) s R ρu t u u t )xτ στ) s u L 4 ρ u t L u t L 4 τ στ) s u 1 L u 1 L u t 1 L u t 3 L τ στ) s u L u L u t L τ u L στ) s u L u t L τ + 1 4Ãt). στ) s u t L τ Along the same line, we write B as B = + + στ) s R ρu u) u u t )xτ ef =B 1 +B +B 3. στ) s R ρu u) u) u t xτ στ) s R ρu u) : )u u t xτ

8 8 M. PAICU, P. ZHANG, AND Z. ZHANG By virtue of Höler inequality an.4), one has B 1 which along with.7) implies στ) s ρ u L 8 u L 4 u t L τ στ) s u 1 L u L u 1 L u t L τ στ) s u L u 4 L u L τ Ãt), { B 1 στ) s u L ρ 1 u t L τ + The same argument gives rise to an B B 3 } στ) s u L u 6 L τ Ãt) στ) s u L ρ 1 ut L τ +C u 4 L A 1 t) Ãt). στ) s ρ 1 u L 4 u L 8 ρ 1 ut L τ στ) s u 1 L u L u 3 L ρ 1 ut L τ στ) s u L u L ρ 1 ut L τ + στ) s u L u L ρ 1 ut L τ +C u L στ) s u 3 L τ στ) s ρ 1 ut 3 L + u 3 L u 6 L ) τ στ) s u L ρ 1 ut L τ +A 1 t) + u 5 L A 1 t) ) Ãt), στ) s ρ 1 u u L ρ 1 u t L τ στ) s u L u L ρ 1 u t L τ στ) s u L ρ 1 u t 3 L τ + u L A 1 t) + u 8 L A 1 t) + Summing up the above estimates, we conclue that στ) s u 3 L u 4 L ρ 1 u t L τ ) στ) s u L ρ 1 u t L τ Ãt)..15) B ) 1+ u L στ) s u L ρu 1 t L τ +C 1+ u 8 L ) A1 t) Ãt).

9 GLOBAL UNIQUE SOLVABILITY OF INHOMOGENEOUS NS EQUATIONS 9 Combining.11) with.1)-.15), we obtain { 1+ u à t) ) L στ) s u L ρ 1 ut L τ + 1+ u 8 ) } L A 1 t)+a 1 t). Applying Gronwall s inequality an.1) leas to à t) C 1+ u 8 L ) u H s + u 4 H s)exp{ C u 4 L }, which together with.7) ensures that A t) à t)+ u L A 1t) ) 1+ u 8 L ) u H s + u 4 H s)exp{ C u 4 L }, This together with.1) an.1) completes the proof of 1.3)... Existence of the solution in 3-D. Proof to the existence part of Theorem 1.. By mollifying the initial ensity ρ, we euce from [4] that 1.1) has a unique global solution ρ ε,u ε ) provie that ε is small enough in 1.5). Then the existence part of Theorem 1. follows from the uniform estimate 1.6) for ρ ε,u ε ) an a stanar compactness argument. For simplicity, we only present the a priori estimates 1.6) for smooth enough solutions ρ,u) of 1.1). The uniqueness of such solution will be prove in Section 3. As a convention in the rest of this section, we shall always enote by C a constant epening on C in 1.4), which may be ifferent from line to line. First of all, it is easy to check from 1.1) an 1.4) that.16) c ρt,x) C for t,x) [,+ ) R 3, T ρ ut,x) x+ R 3 u xt = ρ u x. R 3 R 3 While we get by taking the L inner prouct of the momentum equation to 1.1) an u t that.17) ρ u t xτ + ut,x) x u R 3 R 3 L ρu u) u t xτ. R 3 In what follows, we nee to use Gagliaro-Nirenberg inequality 3.18) a L p R 3 ) C a p 1 3 L R 3 ) a By virtue of.18), one has R 3 ρu u) u t xτ 3 p L R 3 ) for p 6. ρ 1 u L 6 u L 3 ρ 1 u t L τ u 3 L u 1 L ρ 1 u t L τ u 3 L u L τ ρ 1 u t L τ. Whereas it follows from the momentum equation of 1.1) an classical estimates on the Stokes system that u L + p L ) ρ 1 u t L + u u L.19) ρ 1 u t L + u 3 ) 1 L + u L,

10 1 M. PAICU, P. ZHANG, AND Z. ZHANG so that R 3 ρu u) u t xτ C from which,.17) an.19), we infer u 6 L τ + 1.) B 1 t) u L +C u L B 1 t). ρ 1 ut L τ, Hence, as long as we choose ε small enough in 1.6), we obtain the estimate for B 1 t) in 1.6). We now turn to the estimate of B t). Inee along the same line to the proof of.11), we have 1 t σt) ρ u t x+ στ) u t xτ R 3 R 3 t = στ) ρ t u t xτ στ) ρ t u u) u t xτ.1) R 3 R 3 t στ) ρu t u) u t xτ + ρ u t xτ R 3 R 3 ef = E +F +G+H. It is obvious to observe that.) H B 1 t) u L. Whereas using ρ t = u ρ an integrating by parts, we get, by applying.18), that E = στ) ρu u t ) u t xτ R 3.3) Notice that στ) u ρ 1 u t L ρ 1 u t L τ στ) u 1 L u 1 L u t L ρ 1 u t L τ u L + u L ) στ) ρ 1 u t L τ G στ) ρ 1 u L 3 u t L 6 ρ 1 ut L τ στ) u 1 L u 1 L u t L ρ 1 ut L τ. στ) u t L τ..3) hols also for G. To eal with F in.1), we get, by using ρ t = u ρ once again an integrating by parts, to write F = στ) ρu u) u u t )xτ R 3 t + στ) ρu u) u) u t xτ R 3 t + στ) ρu u) : )u u t xτ ef = F 1 +F +F 3. R 3

11 GLOBAL UNIQUE SOLVABILITY OF INHOMOGENEOUS NS EQUATIONS 11 By virtue of 1.6) for B 1 t) an.18), one has F 1 Along the same line, we have F στ) ρ u L 6 u L 6 u t L τ στ) u L u L u t L τ στ) u 4 L u L τ B t) C u 6 L B t). The same estimate hols for F 3, as στ) ρ 1 u L 6 u L 6 ρ 1 u t L τ F 3 Therefore, we conclue that.4) F στ) u L u L ρ 1 u t L τ u L + u L ) στ) ρ 1 u t L + u L ) τ. στ) ρ 1 u u L ρ 1 u t L τ στ) u L u L ρ 1 u t L τ. u L + u L ) στ) ρ 1 u t L + u L ) τ +C u 6 L B t). On the other han, notice from.19) that σt) ut) L + p L ) Cσt) ρ 1 ut t) L + ut) 6 L ) C B t)+ u 6 L ), which together with.1)-.4) ensures that { B t) C u L + u ) 1 L στ) ρ ut L + u ) L τ + u 6 L + u L }, applying Gronwall s inequality gives rise to the estimate of B t) in 1.6). This completes the proof of 1.6). Remark.1. Along the same line to the erivation of.), we also get 1 t ut) L + ρ u t x = ρu u) u t x R 3 R 3 u 6 L + 1 ρ u t x, R 3 which gives.5) t ut) L + ρ u t x C u 6 R 3 L. Hence if the initial velocity u oes not satisfy 1.5), we euce from.5) that there exists a positive time T so that u T L ) + u t L T L ) C u L.

12 1 M. PAICU, P. ZHANG, AND Z. ZHANG With the above estimate, we can obtain the estimate B 1 t) an B t) for t T as we i before. This implies the local existence of solutions to 1.1) in 3-D with large ata. 3. Uniqueness of the solution 3.1. More regularity of the solutions. Before we present the proof to the uniqueness part of both Theorem 1.1 an Theorem 1., we nee the following regularity results for the solutions of 1.1) obtaine in Section. Lemma 3.1. Let ρ,u, p) be the global solution of 1.1) obtaine in Theorem 1.. Then for any T R +, one has 3.1) T T T σt) ut) L 6 + pt) L 6 ) t C, ut) t Cmax T 1 4,T 1 ), σt) 1 ut) t C, for some constant C epening only on c,c in 1.4) an u H 1. Proof. We first get by taking iv to.6) that which along with 1.6) an.6) ensures that p = ) 1 iv { ρ t u+u u) }, 3.) ut) L q + pt) L q C t u+u u)t) L q for any q 1, ). However, it follows from Sobolev imbeing theorem that u u)t) L 6 ut) ut) L C ut) 1 L ut) 3 L, from which, 1.6) an 3.), we infer that T σt) ut) L 6 + p L 6 ) t { T C σt) u t L t+ sup ut) L σt) ut) L ) T } ut) L t C. t [,T] This proves the first part of 3.1). Then we euce from it, Gagliaro-Nirenberg inequality an 1.6) that T ut) t T u 1 L u 1 L 6 t u 1 T L T L ) T )1 σt) 1 t )1 σt) ut) T )1 4 L t σt) 6 1 t Cmax T 1 1 4,T ). Along the same line, we can also prove the estimate for T σt)1 ut) t.

13 GLOBAL UNIQUE SOLVABILITY OF INHOMOGENEOUS NS EQUATIONS 13 The -D version of the above lemma is more complicate, which we present as follows. Lemma 3.. Let ρ,u, p) be the global solution of 1.1) obtaine in Theorem 1.1. Then for any T R + an α [,1), one has 3.3) T T T σt) 1+α s t u, u)t) L 1 α ut) t Cmax T s 1+α),T 1 ), σt) 1 s 1+α ut) t C, where the constant C epens on α, c,c in 1.) an u H s. + pt) ) t C, L1 α Proof. Fist of all, for any β,1), we euce from Gagliaro-Nirenberg inequality that a L 1 β C a 1 β L a β L for any a H 1 R ), which along with 1.3) ensures that for any α,β,γ,1) satisfying β +γ = 1+α 1 s σt) 1+α) ρu u)t) L 1 α L T L T 3.4) C σt) 1 s 1+α) ut) ut) L 1 β L1 γ ut) 1 β C sup t [,T] Cc,C, u H s). σt) 1 s L β ut) β ) ut) 1 γ σt) 1 s L L γ ut) γ L L T Along the same line, we obtain the same estimate for σt) 1+α s t u L it follows from 3.) that 1 α ut) L 1 α + pt) L 1 α u t +u u)t) L 1 α, from which an 3.4), we obtain the first inequality of 3.3). Whereas we get, by using Gagliaro-Nirenberg inequality once again, that T ut) t T u α 1+α ut) ut) L α 1+α L T L ) 1 1+α L1 α t 1+α s σt) ut) L 1 1+α 1 α L T T L. On the other han, T )1 σt) 1+α s 1+α t which together with 1.3) ensures the secon inequality of 3.3). Along the same line, we can prove the last inequality in 3.3). 3.. Lagrangian formulation. As in [9, 1], we shall apply Lagrangian approach to prove the the uniqueness of the solutions. We remark that even with 3.1) an 3.3), the solution of 1.1) obtaine in Theorem 1.1 an Theorem 1. oes not satisfy the assumptions require by Theorem 1 of [1] concerning the uniqueness of solutions to 1.1). Fortunately, the iea use to prove Theorem 1 of [1] can be successfully applie here. Let ρ,u,p) be the solution of 1.1) obtaine in Theorem 1.1 an Theorem 1.. Then thanks to 3.1) an 3.3), we can efine the trajectory Xt,y) of ut,x) by t Xt,y) = ut,xt,y)), X,y) = y,,

14 14 M. PAICU, P. ZHANG, AND Z. ZHANG which leas to the following relation between the Eulerian coorinates x an the Lagrangian coorinates y: 3.5) Xt,y) = y + uτ, Xτ, y))τ. Moreover, we euce from 3.1) an 3.3) that we can take T small enough such that 3.6) T ut) t 1. Then for t T, Xt,y) is invertible with respect to y variables, an we enote by Yt, ) its inverse mapping. Let vt,y) ef = ut,x) = ut,xt,y)). One has 3.7) t vt,y) = t ut,x)+ut,x) ut,x), xi u j t,x) = yk v j t,y) xi y k for x = Xt,y), y = Yt,x). Let At,y) ef = Xt,y)) 1 = x Yt,x). So we have 3.8) x ut,x) = A t t,x) y vt,y) an ivut,x) = ivat,y)vt,y)). By the chain rule, we also have 3.9) iv y A ) = A t y. Here an in what follows, we always enote A t the transpose matrix of A. As in [1], we enote 3.1) u ef = A t y, iv u ef = iva ) an u ef = iv u u, ηt,y) ef = ρt,xt,y)), vt,y) ef = ut,xt,y)) an Πt,y) ef = pt,xt,y)). Notice that for any t >, the solution of 1.1) obtaine in Theorem 1.1 an Theorem 1. satisfies the smoothness assumption of Proposition in [1], so that η,v, Π) efine by 3.1) solves t η =, η 3.11) t v u v + u Π =, iv u v =, η,v) t= = ρ,u ), which is the Lagrangian formulation of 1.1). Now we transform the regularity information of the solution in the Eulerian coorinates into those in the Lagrangian coorinates. Lemma 3.3. Let ρ,u, p) be the global solution of 1.1) obtaine in Theorem 1. an η,v,π) be given by 3.1). Then for any t T etermine by 3.6), one has 3.1) τ 1 t v, v)τ) L 3 + Πτ) L 3 ) τ C, A t L 3 ) + τ 1 vτ) τ C, vτ) τ Ct 1 4, for some constant C epening only on c,c in 1.4) an u H 1.

15 GLOBAL UNIQUE SOLVABILITY OF INHOMOGENEOUS NS EQUATIONS 15 Proof. We first euce from 3.5), 3.6) an 3.1) that 3.13) y Xt, ) exp { x uτ) τ } e 1, which together with 3.1) an 3.1) implies that vt) t Ct 1 4 an τ 1 vτ) L τ C. Furthermore, thanks to et Xt,y) ) y = 1, an 1.6), 3.1), one has τ 1 4 Π L t L 3 ) p 1 L t L ) τ 1 p 1 L t L6 ) C. On the other han, it follows from the proof of 3.13) that yxt, ) L p exp { C x uτ) τ } uτ,xτ, )) L p τ uτ, ) L p τ for any p [1, ]. In particular, if we take p = 3 in the above inequality an use 3.1) to get from which an 3.1), we infer y Xt, ) L 3 Ct1 4 u 1 L t L ) τ 1 u 1 L t L6 ) Ct1 4, τ 4 1 v L t L 3 ) τ 1 4 x uτ,xτ, )) L t L3 ) yx t ) + τ 4 1 x u L t ) y X ) t L3 ) ) 1+ τ 1 4 x u L t L 3 ) 1+ xu 1 L t L ) τ 1 x u 1 ) L t L6 ) C. On the other han, thanks to 3.6), we have for t T 3.14) At,y) = x Yt,x) = I+ y Xt,y) I) ) 1 = 1) l y ut,xt,y))t ) l, for x = Xt,y), so that A t L 3 ) xuτ,xτ, )) L 1 t L 3 ) yx t ) Ct 1 4. Finally, it follows from 1.6), 3.1) an 3.7) that l= τ 1 4 τ v L t L 3 ) τ 4 1 τ u L t L 3 ) + τ 4u u 1 L t L 3 ) { τ u 1 L t L ) τ 1 τ u 1 L t L ) + u 1 t L ) τ 1 u 1 t L ) u 1 L t L ) u 1 L t L ) } C. This completes the proof of the lemma.

16 16 M. PAICU, P. ZHANG, AND Z. ZHANG Lemma 3.4. Let ρ,u, p) be the global solution of 1.1) obtaine in Theorem 1.1 an η,v,π) be given by 3.1). Then for any t T etermine by 3.6) an α < s, one has 3.15) τ 1+α s τ v, v)τ) L 1 α A t L1 α) Cts α, τ 1 s 1+α vτ) τ C, + Πτ) ) τ C, L1 α vτ) τ Ct for some constant C epening only on c,c in 1.) an u H s. Proof. The proof is similar to Lemma 3.3. We omit the etails. s 1+α), 3.3. The proof of the uniqueness. We first recall the following lemma from [1]. Lemma 3.5. Let η R ) be a time inepenent positive function, an be boune away from zero. Let R satisfy R t L,T) R ) an ivr L,T) R ). Then the following system has a unique solution v, Π) such that η t v v + Π = f t,x),t) R, ivv = ivr, v t= = v, v T L ) + v t, v, Π) L T L ) C v L + f,r t ) L T L ) + ivr L T L )), where C epens on infη an supη, but inepenent of T. Proof to the uniqueness parts of Theorems 1.1 an 1.. Let ρ i,u i, p i ),i = 1,, be two solutions of 1.1) obtaine in Theorem 1.1 an Theorem 1., an η i,v i, Π),i = 1,, be etermine by ef 3.1). We enote A i = Au i ) for i = 1,, an δv ef = v v 1,δΠ ef = Π Π 1 an δa = A A 1, then we euce from 3.11) that 3.16) with 3.17) We enote ρ t δv δv + δπ = δf 1 +δf, ivδv = ivδg, δv t= =, δf 1 ef = [ u1 )Π 1 u )Π ] = I A t ) δπ δa t Π 1, ef δf = [ ] u1 )v 1 u )v =iv [ I A A t ) δv +A 1A t 1 A A t ) v 1], δg ef = I A 1 )v 1 +I A )v =I A )δv δav 1. δet) ef = δv t L ) + t δv, δv, δπ) L t L ). Then we infer from Lemma 3.5 an 3.16) that δet) C δf 1 L t L ) + δf L t L ) + ivδg L t L ) + tδg L t L )).

17 We will show that 3.18) GLOBAL UNIQUE SOLVABILITY OF INHOMOGENEOUS NS EQUATIONS 17 δf 1 L t L ) + δf L t L ) + ivδg L t L ) + tδg L t L ) εt)δet), where the function εt) tens to zero as t goes to zero. With 3.18) being grante, we infer that δet) εt)δet), which ensures the uniqueness of solutions obtaine in Theorem 1.1 an Theorem 1. on a sufficiently small time interval [,T 1 ]. The uniqueness on the whole time [, ) can be obtaine by a bootstrap argument. Now let us turn to the proof 3.18). Inee thanks to 3.6), we can take the time T to be small enough so that T v i τ) τ 1, i = 1,. As a convention in the sequel, we shall always assume that t T. Thanks to 3.14), we write δat) = δvτ ) C j ) t 3.19) 1 Ck 1 j with C i t) ef = v i τ. k 1 j<k The proof of 3.18) will split into the following two cases. Proof of 3.18) in 3-D case. We first euce from 3.14) an 3.1) that I A t ) δπ L t L ) C y v t ) t δπ L t L ) 4δEt). Ct1 While it follows from 3.19) that δa t L 6 ) C which along with 3.1) implies δv t t L 6 ) Ct1 δv L t L ), δa t Π 1 L t L ) τ 1 4δAτ) t L 6 ) τ 1 4 Π 1 τ) L t L 3 ) Ct1 4δEt), so that thanks to 3.17), we obtain 3.) δf 1 L t L ) Ct1 4δEt). Next we hanle ivδg. We first get by applying 3.19) that 3.1) δat) C Hence, we have by 3.1) that δv τ +C δa v 1 ) L t L ) { τ C δv τ v L 1 t L ) + + C{ τ τ 1 4 Ct 1 4δEt). δv τ δv τ δv τ v L 1 t L ) v 1 + v ) τ } v L 1 t L ) δv τ t L ) τ 1 4 v1 L t ) δv τ t L 6 ) τ 1 4 v 1, v ) L t L 3 ) v 1 + v ) τ. 1 ) } 1+ τ 4 v1l t )

18 18 M. PAICU, P. ZHANG, AND Z. ZHANG Along the same line, one has I A ) δv ) L t L ) { t C v τ t L 3 ) δv L t L6 ) + Ct 1 4δEt). Thus, it follows from 3.9) an 3.17) that 3.) ivδg L t L ) 4δEt). Ct1 To eal with δf, we write v τ t ) δv L t L ) } 3.3) A A t A 1A t 1 ) v 1 = δai A t )+I A 1)δA t) v 1 +δa t +δa) v 1. It is easy to check that iv [ δai A t ) ) ] v 1 L t L ) τ 1 4 δa t L ) I A t ) τ 1 4 v 1 L ) + τ 1 4 δa t L 6 ) A t L 3 ) τ 1 4 v1 L ) + τ 1 4δA t L 6 ) I A t ) τ 1 4 v 1 L L 3 ), which along with 3.1) an 3.1) implies that iv [ δai A t ) ) ] v 1 L t L ) { Ct4 1 1 δv L t L ) v L 1 t ) τ 4 v 1 L ) + τ 1 4 ) v 1 L t L 3 ) } + A t L 3 ) τ 1 4 v 1 L ) CtδEt). 1 The same estimate hols for iv [ I A 1 )δa t) v 1 ] an iv [ δa t +δa) v 1 ]. Hence, To hanle iv [ I A A t ) δv], we write iv [ A A t A 1A t 1 ) v 1] L t L ) Ct1 δet). 3.4) I A A t ) δv = { I A )I A t ) I A t ) I A ) } δv. Then we get, by using 3.1), that iv [ I A A t ) δv ] L t L ) { } C A t L 3 ) I A t ) δv L t L 6 ) + I A ) t ) δv L t L ) Ct 1 4δEt). Summing up, we obtain 3.5) δf L t L ) Ct1 4δEt).

19 GLOBAL UNIQUE SOLVABILITY OF INHOMOGENEOUS NS EQUATIONS 19 Finally, let us turn to the estimate of t δg L t L ). Thanks to 3.19), we have t [δav 1 ] L t L ) { δv)v 1 L t L ) + δv τ v 1, v ) v L 1 t L ) + δv τ t v 1 } L t L ) { δv t L ) v 1 L t ) + τ 1 4 δv τ t L 6 ) τ 4 v 1 1, v ) L t ) v 1 t L 3 ) + τ 1 4 t v 1 L t L )) } 3, which together with 1.6) an 3.1) ensures that t [δav 1 ] L t L ) { u 1 1 L t L ) u 1 1 L t L ) δv t L ) +t 1 4 δv L t L ) τ 1 4 v1, v ) L t ) u 1 1 t L ) u 1 1 t L ) + τ ) } 1 4 t v 1 L t L 3 ) Along the same line, one has t 1 4 δv t L ) +Ct 1 4 δv L t L ) Ct1 4 δet). t [I A )δv] L t L ) { v δv L t L ) + Hence, we arrive at 3.6) + v τ t δv L t L ) v τ v δv L t L ) } { δv t L 6 ) v L t L 3 ) + v L 1 t ) t 1 4 δet). v L t L 3 ) δv ) } t L6 ) + t δv L t L ) t δg L t L ) Ct1 4 δet). Then 3.18) follows from 3.), 3.), 3.5) an 3.6). Proof of 3.18) in -D case. In what follows, we shall always take < α < s. By virtue of 3.14) an 3.15), we have I A t ) δπ L t L ) C an by Gagliaro-Nirenberg inequality, one has 3.7) δa t L α) which along with 3.15) implies that y v t ) t δπ L t L ) Ct s 1+α) δet), δv t t L α) δv α L δv 1 α L t Ct 1+α δet), δa t 1+α s Π 1 L t L ) τ 1+α s δaτ) t Lα) τ Π 1 τ) L t L1 α) δet). Cts As a consequence, we obtain 3.8) δf 1 L t L ) Ct s 1+α) δet).

20 M. PAICU, P. ZHANG, AND Z. ZHANG While ue to 3.1), we get, by using 3.7) an 3.15), that { τ δa v 1 ) L t L ) C δv τ v L 1 t L ) + an we also have + C δv τ { τ 1+α s 1+α) + τ 1+α s 1+ τ 1+α s 1+α) v 1 L t ) Ct s 1+α) δet), I A ) δv ) L t L ) { t C v τ t L 1 α) δv L t L α) + Ct s 1+α) δet). Here we use the fact that 3.9) Hence, we obtain δv τ v L 1 t L ) } v 1 + v ) τ v 1 L t L ) δv τ 1+α s τ 1+α) v t L ) 1 L t ) δv τ 1+α s t L α) τ v 1, v ) L t L 1 α) ) } v τ t ) δv L t L ) δv L t L α) Ctα δv α t L ) δv 1 α L t L ) Ctα δet). 3.3) ivδg L t L ) Ct s 1+α) δet). Now we hanle δf L t L ). Inee it follows from 3.15) an 3.7) that iv [ δai A t )) v 1 ] L t L ) τ 1+α s 1+α) δa t L ) I A t ) τ 1+α) v 1+α s 1 L ) + τ 1+α s 1+α)δA t Lα) A 1+α s t L 1 α) τ 1+α) v 1 L ) + τ 1+α s δa t L α) I A t ) τ 1+α s v 1 L L 1 α) Ct s 1+α) δet). The same estimate hols for the remaining terms in3.3). Hence, we get iv [ A A t A 1 A t 1) v 1 ] L t L ) Ct s 1+α) δet). Whereas thanks to 3.4), we get, by using 3.15), that iv I A A t ) δv ) L t L ) C A t L1 α) I A t ) δv L t Lα) + I A ) t ) ) δv L t L ) Ct So we obtain s 1+α) δet). 3.31) δf L t L ) Ct s 1+α) δet). }

21 GLOBAL UNIQUE SOLVABILITY OF INHOMOGENEOUS NS EQUATIONS 1 Finally we eal with t δg L t L ). As in the 3-D case, we have t [ δav1 ] L t L ) { δv)v 1 L t L ) + + δv τ t v 1 L t L ) By 3.15), the first term on the right han sie is boune by δv τ v 1, v ) v 1 L t L ) δv t L ) v 1 L t ) CδEt) v 1 1 t L ) v 1 1 L 1 t L ) Cts 4 δet), an the thir term is boune by τ 1+α s δv τ 1+α s t L α) τ t v 1 L t L1 α) Cts δet), an the secon term is boune by τ 1+α s 1+α) 1 s)α δv τ 1+α s 1 s)α t L α) τ 1+α) v 1, v ) L t ) τ v 1 t L 1 α) Ct s1+α+α ) 1+α) δet). On the other han, it follows from Gagliaro-Nirenberg inequality that τ 1 s)α δv t Lα) τ 1 s δv α t L ) δv 1 α t L ) }. t sα t δv α L t L ) δv 1 α t L ) Ctsα δet), τ 1 s)α v L t L1 α) v 1 α 1 s τ L t L ) v α L t L ) C, from which an 3.15), we infer that [ t I A )δv ] L t L ) { C v δv L t L ) + C v τ v δv L t L ) + v τ t δv L t L ) { τ 1 s)α 1 s)α δv t Lα) τ v L t L1 α) + v L 1 t ) τ 1 s)α 1 s)α δv t Lα) τ v L t L1 α) + ) } tδv L t L ) Ct sα δet). Therefore, we arrive at 3.3) t δg L t L ) Ctsα δet). Then 3.18) in the -D case follows from 3.8)-3.3). } Acknowlegments. The authors woul like to thank the anonymous referees for many profitable suggestions. Part of this work was one when M. Paicu was visiting Morningsie Center of the Chinese Acaemy of Sciences in the Spring of 1. We appreciate the hospitality of MCM an the financial support from the Chinese Acaemy of Sciences. P. Zhang is partially supporte by NSF of China uner Grant an 19317, the one hunre talents plan from Chinese Acaemy of Sciences uner Grant GJHZ89 an innovation grant from National Center for Mathematics an Interisciplinary Sciences. Z. Zhang is partially supporte by NSF of China uner Grant an

22 M. PAICU, P. ZHANG, AND Z. ZHANG References [1] H. Abii, Équation e Navier-Stokes avec ensité et viscosité variables ans l espace critique, Rev. Mat. Iberoam., 3) 7), [] H. Abii an M. Paicu, Existence globale pour un fluie inhomogéne, Ann. Inst. Fourier Grenoble), 57 7), [3] H. Abii, Guilong Gui an Ping Zhang, Stability to the global large solutions of the 3-D inhomogeneous Navier- Stokes equations, Comm. Pure. Appl. Math., 64 11), [4] H. Abii, G. Gui an P. Zhang, On the wellposeness of 3 D inhomogeneous Navier-Stokes equations in the critical spaces, Arch. Ration. Mech. Anal., 4 1), [5] S. N. Antontsev, A. V. Kazhikhov, an V. N. Monakhov, Bounary value problems in mechanics of nonhomogeneous fluis, Translate from the Russian. Stuies in Mathematics an its Applications,. North-Hollan Publishing Co., Amsteram, 199. [6] H. Bahouri, J. Y. Chemin an R. Danchin, Fourier analysis an nonlinear partial ifferential equations, Grunlehren er mathematischen Wissenschaften 343, Springer-Verlag Berlin Heielberg, 11. [7] R. Danchin, Density-epenent incompressible viscous fluis in critical spaces, Proc. Roy. Soc. Einburgh Sect. A, 133 3), [8] R. Danchin, Local an global well-poseness results for flows of inhomogeneous viscous fluis, Av. Differential Equations, 9, 4, [9] R. Danchin an P. B. Mucha, A Lagrangian approach for the incompressible Navier-Stokes equations with variable ensity, Comm. Pure. Appl. Math., 65 1), [1] R. Danchin an P. B. Mucha, Incompressible flows with piecewise constant ensity, arxiv: v1. [11] P. Germain, Strong solutions an weak-strong uniqueness for the nonhomogeneous Navier-Stokes system, J. Anal. Math., 15 8), [1] L. Grafakos, Classical an Moern Fourier analysis, GTM 5, Springer, New York, 9. [13] J. Huang, M. Paicu an P. Zhang, Global wellposeness to incompressible inhomogeneous flui system with boune ensity an non-lipschitz velocity, preprint 1. [14] D. Hoff, Global solutions of the Navier-Stokes equations for multiimensional compressible flow with iscontinuous initial ata, J. Differential Equations, ), [15] D. Hoff, Dynamics of singularity surfaces for compressible, viscous flows in two space imensions, Comm. Pure Appl. Math., 55 ), [16] O. A. Layženskaja an V. A. Solonnikov, The unique solvability of an initial-bounary value problem for viscous incompressible inhomogeneous fluis. Russian) Bounary value problems of mathematical physics, an relate questions of the theory of functions, 8, Zap. Naučn. Sem. Leningra. Otel. Mat. Inst. Steklov. LOMI), ), 5 19, [17] P. L. Lions, Mathematical topics in flui mechanics. Vol. 1. Incompressible moels, Oxfor Lecture Series in Mathematics an its Applications, 3. Oxfor Science Publications. The Clarenon Press, Oxfor University Press, New York, [18] M. Paicu an P. Zhang, Global solutions to the 3-D incompressible inhomogeneous Navier-Stokes system, J. Funct. Anal., 6 1), [19] J. Simon, Nonhomogeneous viscous incompressible fluis: existence of velocity, ensity, an pressure, SIAM J. Math. Anal., 1, 199, M. PAICU) Université Boreaux 1, Institut e Mathématiques e Boreaux, F-3345 Talence Ceex, France aress: marius.paicu@math.u-boreaux1.fr P. ZHANG) Acaemy of Mathematics & Systems Science an Hua Loo-Keng Key Laboratory of Mathematics, The Chinese Acaemy of Sciences, Beijing 119, CHINA aress: zp@amss.ac.cn Z. ZHANG) School of Mathematical Science, Peking University, Beijing 1871, P. R. CHINA aress: zfzhang@math.pku.eu.cn