arxiv: v1 [math.ap] 1 Jan 2015

Transcription

1 EXISTENCE AND UNIQUENESS OF WEAK SOLUTIONS OF THE COMPRESSIBLE SPHERICALLY SYMMETRIC NAVIER-STOKES EQUATIONS XIANGDI HUANG arxiv:5.262v [math.ap] Jan 25 Abstract. One of the most influential fundamental tools in harmonic analysis is Riesz transform. It maps L p functions to L p functions for any p, which plays an important role in singular operators. As an application in fluid dynamics, the norm equivalence between u L p and divu L p+ curlu L p is well established for p,. However, since Riesz operators sent bounded functions only to BMO functions, there is no hope to bound u L in terms of divu L + curlu L. As pointed out by Hoff[SIAM J. Math. Anal. 3726, No. 6, ], this is the main obstacle to obtain uniqueness of weak solutions for isentropic compressible flows. Fortunately, based on new observations, see Lemma 2.2, we derive an exact estimate for u L 2+/N divu L for any N-dimensional radially symmetric vector functions u. As a direct application, we give an affirmative answer to the open problem of uniqueness of some weak solutions to the compressible spherically symmetric flows in a bounded ball.. Introduction and main results We are concerned with the isentropic system of compressible Navier-Stokes equations which reads as ρ t + divρu=,. ρu t + divρu U+ P=µ U+ µ+λ divu, where t, x Ω R N N= 2, 3,ρ=ρt, x and U= Ut, x are the density and fluid velocity respectively, and P=Pρ is the pressure given by a state equation.2 Pρ=aρ γ with the adiabatic constant γ > and a positive constant a. The shear viscosity µ and the bulk oneλare constants satisfying the physical hypothesis.3 µ>, µ+ N 2 λ. The domainωis a bounded ball with a radius R, namely,.4 Ω= B R ={x R N ; x R< }. We study an initial boundary value problem for. with the initial condition.5 ρ, U, x=ρ, U x, x Ω, and the boundary condition.6 Ut, x=, t, x Ω, 2 Mathematics Subject Classification. 35Q3, 76N

2 2 XIANGDI HUANG and we are looking for the smooth spherically symmetric solution ρ, U of the problem.,.5,.6 which enjoys the form.7 ρt, x=ρt, x, Ut, x=ut, x x x. Then, for the initial data to be consistent with the form.7, we assume the initial data ρ, U also takes the form.8 ρ =ρ x, U = u x x x. In this paper, we further assume the initial density is uniformly positive, that is,.9 ρ =ρ x ρ>, x Ω for a positive constantρ. Then it is noted that as long as the classical solution of.,.5,.6 exists the density ρ is positive, that is, the vacuum never occurs. It is also noted that since the assumption.7 implies. Ut, x+ut, x=, x Ω, we necessarily have Ut, =also U =. There are many results about the existence of local and global strong solutions in time of the isentropic system of compressible Navier-Stokes equations when the initial density is uniformly positive refer to [, 3, 5, 6, 23 25, 28, 29] and their generalization [2 22, 27] to the full system including the conservation law of energy. On the other hand, for the initial density allowing vacuum, the local well-posedness of strong solutions of the isentropic system was established by Kim [7]. For strong solutions with spatial symmetries, the authors in [8] proved the global existence of radially symmetric strong solutions of the isentropic system in an annular domain, even allowing vacuum initially. However, it still remains open whether there exist global strong solutions which are spherically symmetric in a ball. The main difficulties lie on the lack of estimates of the density and velocity near the center. In the case vacuum appears, it is worth noting that Xin [3] established a blow-up result which shows that if the initial density has a compact support, then any smooth solution to the Cauchy problem of the full system of compressible Navier-Stokes equations without heat conduction blows up in a finite time. The same blowup phenomenon occurs also for the isentropic system. Indeed, Zhang-Fang [3],Theorem.8 showed that if ρ, U C [, T]; H k k>3 is a spherically symmetric solution to the Cauchy problem with the compact supported initial density, then the upper limit of T must be finite. On the other hand, it s unclear whether the strong classical solutions lose their regularity in a finite time when the initial density is uniformly away from vacuum. On the other hand, there are amount of literatures investigating the global existence of weak solutions to the compressible Navier-Stokes equations, such as finite energy solutions proposed and developed by Lions [9], Hoff [7] and Feireisl [4], etc. One remarkable result is due to Jiang-Zhang [4], where they prove a global existence for three-dimensional compressible spherically symmetric flows with γ >. However, whether their weak solution is unique remains a long-standing open problem. Meanwhile, Desjardins [3] built a more regular weak solution for three-dimensional torus. Inspired by his work and a new observation in Lemma 2.2, we will give some positive answer to the spherically symmetric flow in this paper.

3 SPHERICALLY SYMMETRIC FLOWS 3 In the spherical coordinates, the original system. under the assumption.7 takes the form ρ t + ρu r + N ρu r =,. ρu t + ρu 2 + Pρ + N ρu2=κ u r + N u r r r r whereκ=2µ+λ. Now, we consider the following Lagrangian transformation:.2 t=t, y= Then, it follows from. that ρt, s s N ds..3 y t = ρur N, r t = u, r y = ρr N, and the system. can be further reduced to ρ t +ρ 2 r N u y =,.4 r N u t + P y =κ ρr N u y where t, y [, M ] and M is defined by.5 M = R ρ r r N dr= according to the conservation of mass. Note that R.6 rt, =, rt, M =R. We denote by E the initial energy.7 E = R y ρt, r r N dr, u 2 ρ 2 + aργ r N dr, γ and define a cuboid Q t,y, for t and y [, M ], as.8 Q t,y = [, t] [y, M ]. From now on, we denote for p [, and a radially symmetric function f=f x.9 f L p Ω= f p dx p Ω where =ωn R f p r N p dr.2 ω N = N B, B stands for the volume of N-dimensional unit ball. We first prove a local existence of weak solutions to the compressible Navier-Stokes equations in Theorem.. The initial data are supposed to satisfy.8-.9 and ρ L Ω.2 U H Ω N.

4 4 XIANGDI HUANG Theorem.. AssumeΩ= B R is a bounded ball inr N for N= 2, 3 andγ>, then there exists T,+ ] and a weak solution ρ, U to the Navier-Stokes equations. in, T such that for all T< T, ρ L, T Ω, ρ U=ρU t + U U, U L 2, T Ω N,.22 divu L 2, T; L Ω, U L, T; L 2 Ω N L 2, T; L Ω Remark.. The key idea to establish local existence of weak solution with regularity.22 is to derive uniform upper bound of the density, which is analogous to Desjardins [3]. However, there are two obstacles in bounded domain. First of all, due to the lack of commutator estimates, whether there is local weak solution with higher regularity.22 for system. remains unknown for initial boundary value problem. We rewrite it in spherically coordinate and derive some new estimates to play a critical role as commutator estimates frequently used by [3, 9] for Cauchy problem and torus. On the other hand, general global finite energy weak solution withγ> was proved by Jiang-Zhang [4]. However, their weak solution is much less regular than.22 and leave a challenging problem on uniqueness of such weak solutions. The main value in Theorem. is to weak the assumption fromγ>3 in [3] toγ> for three-dimensional spherically symmetric flow, and further more, give an affirmative answer on uniqueness of weak solutions stemming from regularity class.22, which is our main issue in Theorem.2. The technical part lies on the combination of Caffarelli-Kohn-Nirenberg [2] inequalities with weights and pointwise estimates for radially symmetric functions, see Lemmas Remark.2. From many early works on the blowup criterion [ 2, 26] of strong solutions to the compressible Navier-Stokes equations, the uniform bound of the density induces the regularity ofρ U, U, divu and curlu, as indicated by the first three lines in Theorem.. Besides, one of the most important observations is Lemma 2.2, which gives desired bound for U L in terms of divu L. This is a key ingredient in proving the uniqueness of weak solutions illustrated in the following theorem. Theorem.2. Let ρ i, U i be two weak solutions of obtained by Theorem.. Then.23 ρ =ρ 2, U = U 2 a.e on, T Ω. Remark.3. There are many weak-strong uniqueness results [5, 6, 8] concerning compressible Navier-Stokes equations. However, weak-weak uniqueness incorporates more difficulties and need special attention. Up to now, the most far reaching result appears in [8]. As pointed out by Hoff [8], the main obstacle to prevent us from establishing uniqueness is whether we can prove U L L instead of U L BMO. Such a fact was first verified by Hoff [9], where the Lipschitz regularity of the velocity U was established with piecewise C α density. This is a only result concerning uniqueness of weak solutions. With the help of Lemma 2.2 and Theorem., we also give an affirmative answer to the spherically symmetric case. 2. Proofof Theorem. First we recall the following famous Caffarelli-Kohn-Nirenberg [2] inequalities with weights.

5 SPHERICALLY SYMMETRIC FLOWS 5 Lemma 2. Caffarelli-Kohn-Nirenberg. There exists a positive constant such that the following inequality holds for all u C Rn 2. x γ u L r C x α Du a x β u a L p L q if and only if the following relations hold: 2.2 r +γ n = a p +α + a n q +β n 2.3 α σ if a>, and 2.4 α σ if a> and p +α = n r +γ n. satisfying p, q, r>, a 2.5 p +α n, q +β n, r +γ n >, γ=aσ+ aβ Furthermore, on any compact set in which 2.5 and α σ hold, the constant C is bounded. The following lemma is essential in proving Theorems. and.2. Lemma 2.2. AssumeΩis either a bounded ball B R with radius R or the whole space R N, then for any N-dimensional radially symmetric vector functions Ux=u x x x for x Ω, we have the following estimates We have 2.6 N divu L Ω U L Ω 2+ N divu L Ω. 2 For p [,, we have u r N p r N p ω p 2.7 where u r N divu L p Ω +N N p r N p ω p N 2.8 ω N,r = N B. Proof. Set r = x, obviously, divu L p Ω 2.9 divu= u r + N u, curlu=. r Denote by 2. u r + N u r = F. It follows from 2. that 2. r N u r = r N F,

6 6 XIANGDI HUANG which gives 2.2 r N u= Consequently, s N Fds F L s N ds= rn N F L. 2.3 u r N F L. And 2.4 u r L F L + N u r L 2N N F L. On the other hand, as r xi = x i, we have r 2.5 It immediately implies xi u x j r =u r x i x j r 2 + u x j r x i = x ix j r u r+ δ 2 i j x ix j r u 2 r. 2.6 U L Ω u r L Ω+ 2 u r L Ω 2+ N F L Ω. For p,, r,, one obtains 2.7 r N u= s N Fds p p F p s N ds s N ds Hence, rn N p ω p N F L p Ω 2.8 u N p r N p ω p N F L p Ω Therefore, 2.7 follows from 2.8 and 2.9. This finishes the proof of Lemma. We only prove the case when N= 3 since the case N= 2 is even simpler. Throughout of this section, we assume that ρ, U with the form.7 is the solution to the initial boundary value problem.,.5,.6 in [, T] Ω, and we denote by C generic positive constants only depending on the initial data and time T. We give a sketch of proof of Theorem., since the main idea can be borrowed from Desjardins [3]. Denote 2.9 Φt=+ ρ L Ω+ P 2 L 2 Ω + U 2 L 2 Ω. Our main procedure is to derive the following estimate for t<, 2.2 Φ C+ C exp C exp χφλsds. for some positive increasing smooth funtion χx and integrable function λs. The existence of T < follows immediately from 2.2.

7 SPHERICALLY SYMMETRIC FLOWS 7 Say denoting byζt the right hand side of 2.2, we conclude that d 2.2 dt ζt Cπζtλt, for some smooth functionπ. Thus, we have 2.22 ζt du πu so that there exists T such that for all T< T <, 2.23 Φ C T. λsds, We first have the following basic energy estimate. Since the proof is standard, we omit it. Lemma 2.3. It holds for any t T, 2.24 ρ U aργ dx+κ γ Ω or equivalently, 2.25 Denoted by R ρ u2 R 2 + aργ γ u 2 ρ 2 + aργ γ Ω r N dr+κ U 2 dxdτ R r N dr=e G=κdivU P Ω ρ U 2 2 2r u2 u + r N drdτ r 2 as effective flux. The momentum equations. 2 can be rewritten as 2.27 ρ U= G. + aργ dx, γ The main difficulty lies in the bound of the density. We will work it in spherical coordinate as system.4 as follows. Lemma 2.4. There exists a smooth positive increasing function χx and an integrable functionλt such that 2.28 ρt, y C exp Cβt+Cexp Cβt+C χφλsds, t, y [, ] [, M ]. where 2.29 βt= P 2 L 2 Ω + U 2 L 2 Ω Remark 2.. From the work of [3], refer to 3, we can prove 2.3 βt C exp C χφλsds. Therefore, the left work is concentrated on Lemma 2.4.

8 8 XIANGDI HUANG Proof. Step. In view of.4, it holds 2.3 κlogρ ty =κ ρ t ρ y= κρr N u y y = r N u t p y = r N u t p y N u2 r N. Thus, integrating 2.3 over, t, y, we deduce that ρt, y y κ log ρt, =κ logρ y ρ + r N u, z r N ut, z dz ps, ps, y ds y N u2 s, z r N which is equivalent to ρt, y y ρt, =ρ y κ ρ exp r N u, z r N ut, z dz exp κ ps, ps, y ds 2.33 y exp κ N u2 s, z dzds. r N = ρ y ρ Π3 i= Ψ i. dzds, We first deal withψ andψ 3. Step 2. It follows from Lemma 2.2 with N= 3, energy equality 2.25 andγ> that 2.34 y r N u dy= where we used the following fact. First recall that ρ u dr C U L 2 Ω C U L 2 Ω ρ 6γ s 2 ds 6γ Cr 3γ 6γ γ P U L 6 Ω L 2 Ω C+ C P 2 + L 6 Ω C U 2 L 2 Ω C+ Cβt+C P 2 dt, 2.35 P 6 t + divp 6 u+6γ Pdivu=, one immediately has 2.36 P 6 P L 6 Ω 6 + C L 6 Ω Ω ρs 2 ds s 3γ+2 6γ 6γ ds P 6 divu dxds C+ C P 2 ds. Hence, 2.37 Ψ,Ψ C exp Cβt+C P 2 ds.

9 SPHERICALLY SYMMETRIC FLOWS 9 Step 3. Similarly, recall Lemma 2.2 and G=κdivU P=κF P 2.38 Hence, r 2 u= s 2 Fds F 6 s 4 ds s 5 8 ds Cr 3 6 F 6 s 4 ds u 2 /r Cr 2 3 F 6 s 4 ds 3 Cr 2 3 G 6 s 4 ds 3 + Cr ρ 6γ. Also, we will use the following CKN inequality easily from Lemma s 2 3 G L 6,r C s G 2 3 L 2,r sg 3 L 2,r or equivalently, 2.4 G 6 s 4 ds C G 4 L 2 B,r G 2 L 2 B,r Indeed, in view of 2., take n=,γ= 2 3, r=6,α=β=>σ= 2, p=q=2, a= 2 3. Now we are ready to give estimates forψ y r N u 2 r dydτ= ρ u 2 dsdτ, r, R s C C C = C C C C r 3 ρ F 6 s 4 ds 3 ρs 2 3 dsdτ G 6 s 4 ds 3 + r 3 ρ 2γ+ dτ 3 ρ G 2 G 3 4 L 2 B,r H B,r + ρ 2γ+ 3 ρ G 2 G 3 4 L 2 B,r H B,r + ρ 2γ+ dτ dτ + ρ U 2 + ρ 3 L 2 B,r U 2 3 G 4 L 2 L 2 B,r + ρ 2γ+ dτ + ρ G 2 L 2 B,r + U 2 L 2 B,r ρ ρ 6γ+ dτ + ρ 2 U 2 + L 2 B,r C ρ 5 + U 2 + ρ 6γ+ L 2 dτ The estimates of ρ 2 U L 2 then follows similarly as 3 in Desjardins [3], one concludes that 2.43 Ψ 3 Ψ 3 C exp C χφλsds.

10 XIANGDI HUANG Step 4. We can rewrite 2.33 as 2.44 ρt, y=ptut, y exp κ where 2.45 Pt= and ps, y ds ρt, t κ ρ exp ps, ds 2.46 Ut, y=ρ yψ Ψ 3 On the other hand, it follows from 2.44 that d γ t 2.47 dt exp ps, y ds = aγ γ t κ κ ρt, yγ exp ps, y ds κ = aγ κ PtUt, yγ, which implies 2.48 exp t ps, yds = + aγ κ κ PsUs, y γ ds /γ. Next, we are in a position to estimatept. First, observe that 2.49 M dy ρt, y = In view of 2.44 and 2.48, we have 2.5 ρt, y = Then,Pt can be estimated as 2.5 R N N Pt= = C M M M + aγ κ + aγ κ Pt ρt, y dy Ut, y dy R r N dr= RN N. PtUt, y PsUs, yγ ds. /γ PsUs, yγ ds /γ Ut, y + C aγ M /γ sup κ /γ Ut, y sup U t, y Ps ds γ dy Q T, Q T, CM sup U t, y Q T, /γ + CM sup Ut, y sup U t, y Ps ds γ. Q T, Q T, dy

11 SPHERICALLY SYMMETRIC FLOWS Using 2.49 and takingγ-th power on both sides of 2.5, we have γ γ M γ Pt γ C sup U t, y E Q T, 2.52 γ + C sup U t, y γ sup Ut, y Ps γ ds. Q T, Q T, Therefore, by Gronwall s inequality, we deduce from 2.52 that E γ Pt C sup U t, y M Q 2.53 T, exp sup CT U t, y γ sup Ut, y γ. Q T, Q T, Finally, recalling 2.5, we have 2.54 ρt, y Pt sup Ut, y, Q T, Plugging 2.37,2.43 and 2.46 into , thus finishes the proof of Lemma 2.4. Therefore, Theorem. can be done easily from the work of [3] and Lemma Proofof Theorem.2 The remaining part is dedicated to the proof of Theorem.2. Proof. Set v i = ρ i, Ui = u i x r. We prove Theorem.2 for the system v t = r N u y, 3. r N r N u y u t + p y =κ v where p=av γ. The two weak solutions enjoy the following regularity <ρ ρ i ρ, 3.2 u i r, ui r L2 L. This is equivalent to say <v v i v, 3.3 r N u i y, ui r L2 L. Indeed, 3.4 r N u i y= r N u i rr y = v i u i r L 2 L. y

12 2 XIANGDI HUANG Denote by 3.5 Λ=v v 2,,Θ=u u 2. We have 3.6 Λ 2v, Λ, y=θ, y=. and 3.7 Λ t = r N Θ y, 3.8 r N Θ t + p p 2 y =κ rn u y v =κ rn Θ y v rn u 2 y v 2 y Λ v v 2 rn u 2 y y. Multiplying 3.7 byλand integrating over, y, we have d 3.9 Λ 2 dy ǫ r N Θ y 2 dy+cǫ Λ 2 dy. 2 dt Similarly, multiplying 3.8 by r N Θ and integrating over, y, we also obtain d Θ 2 dy+ r N Θ 2 dt v 3. y 2 dy ǫ r N Θ y 2 dy+cǫ Λ 2 dy+cv, v r N u 2 y 2 L Λ 2 dy. where we have used Young s inequality and 3. p p 2 = v γ v 2 γ Cv, v Λ. Recalling the fact 3.3 and.3 that 3.2 r N u 2 y = r N u 2 y + N vu2 r L2 L. Collecting and choosing ǫ small enough, Gronwall s inequality immediately imply 3.3 Λ=Θ=. Acknowledgment. This research is supported by President Fund of Academy of Mathematics Systems Science, CAS, No.24-cjrwlzx-hxd and National Natural Science Foundation of China, No and References [] H. Beira da Veiga, Long time behavior for one-dimensional motion of a general barotropic viscous fluid. Arch. Rational Mech. Anal. 8989, 4-6. [2] Caffarelli L, Kohn R, Nirenberg L, First order interpolation inequalities with weights. Compositio. Math , [3] B. Desjardins, Regularity of weak solutions of the compressible isentropic navier-stokes equations. Comm. Partial Diff Eqs , No.5, [4] E. Feireisl, A. Novotny, H. Petzeltová, On the existence of globally defined weak solutions to the Navier-Stokes equations. J. Math. Fluid Mech. 3 2, [5] E. Feireisl, Weak-strong uniqueness property for the full Navier-Stokes-Fourier system. Arch. Rational Mech. Anal ,

13 SPHERICALLY SYMMETRIC FLOWS 3 [6] P. Germain, Weak-strong uniqueness for the isentropic compressible Navier-Stokes system.j. Math. Fluid Mech. 32, No., [7] D. Hoff, Discontinuous solutions of the Navier-Stokes equations for multidimensional flows of heat-conducting fluids. Arch. Rational Mech. Anal , [8] D. Hoff, Uniqueness of weak solutions of the Navier-Stokes equations of multidimensional, compressible flow.siam J. Math. Anal. 3726, No. 6, [9] D. Hoff, Dynamics of Singularity Surfaces for Compressible, Viscous Flows in Two Space Dimensions.Comm. Pure. Appl Math. VOL.LV22, [] Huang, X. D., Li, J., Xin, Z. P.: Blowup criterion for viscous barotropic flows with vacuum states. Comm. Math. Phys. 32, [] Huang, X. D., Li, J., Xin, Z. P.: Serrin type criterion for the three-dimensional viscous compressible flows. SIAM J. Math. Anal. 432, [2] Huang, X. D., Xin, Z. P.: A Blow-Up Criterion for Classical Solutions to the Compressible Navier- Stokes Equations, Sci. in China. 5332, [3] N. Itaya, On the Cauchy problem for the system of fundamental equations describing the movement of compressible viscous fluids,kodai Math. Sem. Rep , 6-2. [4] Jiang. S, Zhang. P, On Spherically Symmetric Solutions of the Compressible Isentropic NavierC- Stokes Equations,Comm. Math. Phys , [5] Vaigant, V. A.; Kazhikhov. A. V. On existence of global solutions to the two-dimensional Navier- Stokes equations for a compressible viscous fluid. Sib. Math. J , no.6, [6] A.V. Kazhikhov, Stabilization of solutions of an initial-boundary-value problem for the equations of motion of a barotropic viscous fluid. Differ. Equ. 5979, [7] Y. Cho, H. Kim, On classical solutions of the compressible Navier-Stokes equations with nonnegative initial densities. Manuscript Math. 2 26, [8] H.J, Choe, H. Kim, Global existence of the radially symmetric solutions of the Navier-Stokes equations for the isentropic compressible fluids. Math. Meth. Appl. Sci. 2825, -8. [9] P. L. Lions, Mathematical topics in fluid mechanics. Vol. 2. Compressible models, Oxford University Press, New York, 998. [2] A. Matsumura, T. Nishida, The initial value problem for the equations of motion of compressible and heat-conductive fluids, Proc. Japan Acad. Ser. A Math. Sci , [2] A. Matsumura, T. Nishida, The initial value problem for the equations of motion of viscous and heat-conductive gases, J. Math. Kyoto Univ. 2 98,67-4. [22] A. Matsumura, T. Nishida, The initial boundary value problems for the equations of motion of compressible and heat-conductive fluids, Comm. Math. Phys , [23] J. Nash, Le probleme de Cauchy pour les equations differentielles dn fluide general, Bull. Soc. Math. France 9 962, [24] R. Salvi, I. Straskraba, Global existence for viscous compressible fluids and their behavior as t. J. Fac. Sci. Univ. Tokyo Sect. IA, Math , 7-5. [25] V.A. Solonnikov, Solvability of the initial boundary value problem for the equation of a viscous compressible fluid, J. Sov. Math. 4 98, [26] Sun, Y. Z., Wang, C., Zhang, Z. F.: A Beale-Kato-Majda Blow-up criterion for the 3-D compressible Navier-Stokes equations. J. Math. Pures Appl., 952, [27] A. Tani, On the first initial-boundary value problem of compressible viscous fluid motion, Publ. Res. Inst. Math. Sci. Kyoto Univ. 3 97, [28] A. Valli, An existence theorem for compressible viscous fluids, Ann. Mat. Pura Appl. IV 3 982, 97-23; [29] A. Valli, Periodic and stationary solutions for compressible Navier-Stokes equations via a stability method, Ann. Scuola Norm. Sup. Pisa Cl. Sci. 983, [3] Z. P. Xin, Blowup of smooth solutions to the compressible Navier-Stokes equation with compact density. Comm. Pure Appl. Math , [3] Zhang, T; Fang D. Compressible flows with a density-dependent viscosity coefficient. SIAM J. Math. Analysis, 4 29, no.6, Xiangdi Huang NCMIS, AMSS, Chinese Academy of Sciences, Beijing

14 4 XIANGDI HUANG 9, People s Republic of China xdhuang@amss.ac.cn