• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    On the Global Well-Posedness of 3-D Boussinesq System with Variable Viscosity*

    2019-09-26 02:54:18HammadiABIDIPingZHANG

    Hammadi ABIDI Ping ZHANG

    (Dedicated to Professor Andrew J.Majda for the 70th birthday)

    Abstract In this paper,the authors first consider the global well-posedness of 3-D Boussinesq system,which has variable kinematic viscosity yet without thermal conductivity and buoyancy force,provided that the viscosity coefficient is sufficiently close to some positive constant in L∞and the initial velocity is small enough in(R3).With some thermal conductivity in the temperature equation and with linear buoyancy force θe3 on the velocity equation in the Boussinesq system,the authors also prove the global well-posedness of such system with initial temperature and initial velocity being sufficiently small in L1(R3)and(R3)respectively.

    Keywords Boussinesq systems,Littlewood-Paley theory,Variable viscosity,Maximal regularity of heat equation

    1 Introduction

    The purpose of this paper is to investigate the global well-posedness to the following threedimensional Boussinesq system with variable kinematic viscosity

    Here θ,u = (u1,u2,u3)stand for the temperature and velocity of the fluid respectively,anddenotes the deformation tensor,Π is a scalar pressure function,and the kinematic viscous coefficient μ(θ)is a smooth,positive and non-decreasing function on [0,∞).The thermal conductivity coefficient ν ≥0,and e3= (0,0,1),ε ≥0,εθe3denotes buoyancy force.Furthermore,in all that follows,we shall always denote |D|sto be the Fourier multiplier with symbol |ξ|sfor s ≥0.

    The Boussinesq system arises from a zeroth order approximation to the coupling between Navier-Stokes equations and the thermodynamic equations.It can be used as a model to describe many geophysical phenomena (see [28]).In the Boussinesq approximation of a large class of flow problems,thermodynamic coefficients such as kinematic viscosity,specific heat and thermal conductivity may be assumed to be constants,leading to a coupled system of parabolic equations with linear second order operators.

    However,there are some fluids such as lubricants or some plasma flow for which this is not an accurate assumption (see [30]),and a quasilinear parabolic system as follows has to be considered:

    One may check[17]and the references therein for more details about(1.2).Furthermore,under some technical assumptions,the global existence of weak solutions to (1.2)and in the case of constant viscosity,the uniqueness of such weak solutions in two space dimension was proved in[17].

    Recently the System(1.2)has attracted a lot of attentions in the field of mathematical fluid dynamics.In particular,in two space dimension,with F(θ)= θe2for e2= (0,1)in (1.2),Wang and Zhang [32] proved the global existence of smooth solutions to (1.2).In this case,even with φ(θ)=0 and μ(θ)=μ>0 in (1.2),Chae [11] and Hou,Li[23] independently proved the global existence of smooth solutions to (1.2),Hmidi and Keraani [20] proved the global existence of weak solutions to (1.2)with θ0,u0belonging to L2(R2)and the uniqueness of such solutions was proved for θ0,u0belonging to Hs(R2)for any s >0,the first author of this paper and Hmidi [3] established the global well-posedness of this system with initial data satisfyingWhen N ≥3,eN= (0,··· ,1),and F(θ)= θeN,φ(θ)= 0 and μ(θ)=μ>0 in (1.2),which corresponds to ν =ε=0 and μ(θ)=μ>0 in (1.1),Danchin and Paicu [14] proved the global well-posedness of this system withandfor p ∈[N,∞)provided that

    for some sufficiently small constant c.

    We should also mention that there are many studies on the so-called Boussinesq system with critical dissipation in two space dimension,which reads

    When ν = 0 and μ >0,the above system is called Boussinesq-Navier-Stokes system with critical dissipation,Hmidi,Keraani and Rousset [21] proved the global well-posedness of such system.When ν >0 andμ=0,the System(1.3)is called Boussinesq-Euler system with critical dissipation,Hmidi,Keraani and Rousset [22] proved the global well-posedness of this system.Very recently even the logarithmically critical Boussinesq system was investigated by Hmidi in[19].There are also studies to the global well-posedness of the anisotropic Boussinesq system(with partial thermal conductivity and partial kinematic viscosity)in two space dimension (see[10,15] for instance).

    On the other hand,Abidi[2]proved the global well-posedness of (1.2)in two space dimension under the assumptions that:φ(θ)= 0,F(θ)= 0,and the initial data satisfiesmoreover for some sufficiently small ε,there holds

    Furthermore,the authors of[6]established the global well-posedness of a 2-D Boussinesq system,which has variable kinematic viscosity and with thermal conductivity of|D|θ,with general initial data provided that the viscosity coefficient is sufficiently close to some positive constant in L∞norm.

    Motivated by [2,6] and the recent results of the authors of [4-5] concerning the global well-posedness of inhomogeneous Navier-Stokes system with variable density(see also[16,24]),we are going to investigate the global well-posedness to the following Boussinesq system with variable viscosity,which corresponds ν =ε=0 in (1.1):

    In all that follows,we always make the convention that,for any α >0,α+means any constant greater than α,and

    Theorem 1.1Letandbe a solenoidal vector filed for someThen there exists a sufficiently small constant ε0,and some small enough constant ε,which depends on,such that if

    (1.4)has a unique global solution (θ,u,?Π)with

    Remark 1.1Let us give the following remarks concerning this theorem:

    (1)We point out that the exact dependence of ε onwill be given by(5.1)and(5.49).Furthermore,compared with the 2-D result in [2],here we do not require any smallness condition on the initial temperature in Theorem 1.1.

    (2)The assumption that u0∈(R3)is to make sure that the solution decay to 0 as time t goes to ∞,which will be essential to obtain the a priori estimate of the velocity u in the space L1(R+,Lip(R3)).In fact,the exact decay rate of u(t)will be given by (5.31)and (5.37).The assumption thatis due to the variable viscosity,one may check (3)of Remark 1.2 for details.

    (3)Compared with the results in[4-5]and[24]for the inhomogeneous Navier-Stokes system with variable viscosity,here ε0is a uniform small positive constant,which does not depend on θ0.While in [4-5]and [24],the smallness condition for μ(ρ0)-1 is in some sense formulated as

    for some δ >0.In general,under the assumption that

    for some ε0sufficiently small,Desjardins[16]only proved the global existence of strong solutions for 2-D inhomogeneous Navier-Stokes system.Yet the uniqueness and regularities of such strong solutions are still open.Moreover,here we do not require the initial velocity u0∈H1as was assumed in [4-5,16,24].

    However,with linear buoyancy force θe3on the right-hand side of the velocity equation in(1.4),we still do not know how to prove such global well-posedness result as Theorem 1.1 for the corresponding system.Yet with some dissipation on the temperature equation,more precisely,for the system below,

    we have the following global well-posedness result.

    Theorem 1.2Let sand αandbe a solenoidal vector field.Then there exist sufficiently small constants η and η0,which depend on,such that if

    (1.9)has a unique global solution (θ,u)with

    Remark 1.2Let us mention the following facts about Theorem 1.2.

    (1)We remark that the exact smallness conditions for η and η0will be given by (7.1),(7.2),(7.12)and(7.20).This result in some sense extends the global well-posedness result in[14]with constant viscosity to the case of variable kinematic viscosity.As a matter of fact,the method of the proof to Theorem 1.2 is also motivated a lot from that in [14,27],namely,we need first to control the L∞(R+;L3,∞(R3))norm for the velocity field u before we deal with the evolution of the Besov norm for u.

    (2)We also remark that the global well-posedness of the System (1.9)is easier for larger s.However,we just choose s ≤1 for simplicity.And the reason why we choose s >3-is due to the facts that when we manipulate the globalfor the velocity in Subsection 7.2,we need θ(t1)belongs towhich is obtained by using the smooth effect of |D|sin the θ equation of (1.9),and psgiven by (5.39)satisfies s >,which is crucial for us to work the paraproduct estimate in (7.28).We do not claim the assumption thatin Theorem 1.2 is optimal in any sense.

    (3)We emphasize that the main difficulty in the proof of Theorems 1.1 and 1.2 is due to the variable viscosity.In this case,when we apply Littlewood-Paley theory and smoothing effect of heat semigroup to prove the global L1in time of the space Lipschize estimate for the velocity,we need some positive space derivative estimate of θ.Yet to propagate the positive space derivative estimate for θ,we require the the global L1in time of the space Lipschize estimate for the convection velocity u.Especially for the transport equation of (1.4),one has

    which makes impossible to close the a priori estimates.

    Let us complete this section with the notations we are going to use in this context.

    NotationsLet A,B be two operators,we denote [A;B] = AB -BA,the commutator between A and B.For ab,we mean that there is a uniform constant C,which may be different on different lines,such that a ≤Cb.We shall denote by (a | b)(or (a | b)L2)the L2(R3)inner product of a and b,and denote by (dj)j∈Z(resp.(cj)j∈Z)a generic element of?1(Z)(resp.?2(Z))so that ‖(dj)j∈Z‖?1(Z)=1 (resp.‖(cj)j∈Z‖?2(Z)=1).

    For X a Banach space and I an interval of R,we denote by C(I; X)the set of continuous functions on I with values in X,and by Cb(I; X)the subset of bounded functions of C(I; X).For q ∈[1,+∞],the notation Lq(I; X)stands for the set of measurable functions on I with values in X,such that t‖f(t)‖Xbelongs to Lq(I).Finally for any vector field v = (v1,v2,v3),we denoteand the Leray projection operatordiv.

    2 Strategies to the Proof of Theorems 1.1 and 1.2

    As the existence part of both Theorems 1.1 and 1.2 basically follows from the a priori estimates for smooth enough solutions of (1.4)and (1.9).We shall only outline the main steps in the derivation of the estimates.

    2.1 Strategy to the proof of Theorem 1.1

    By applying maximal regularity estimates for heat semi-group,we prove that under the assumption of (1.6),for smooth enough solution (θ,u)of (1.4)on [0,T*[,there holds

    for any p ∈[6,8].If one assumes moreover the smallness condition (5.1),we get,by using Littlewood-Paley analysis that T*≥1 and

    With (2.2),we can prove the propagation of regularities of (θ0,u0)for (θ,u)on [0,1],namely the estimate (5.14).This in particular ensures some t0∈]0,1[ so that there holds

    Due to (2.3),we can prove the following Desjardins type (see [16])energy estimates for t ∈]t0,T*[,

    for C given by (5.31).Based on (2.5),we deduce that for p ∈[6,8] and satisfying (5.39),

    By virtue of (2.4)and (2.6),we can prove the following key estimate by applying the smoothing effect of heat semi-group and Littlewood-Paley theory

    for any t ∈]0,T*[,where ‖u(t0)‖H1is determined by (2.3)andby (5.31).

    This is basically the contents of Section 5.With (2.7),we shall complete the proof to the existence part of Theorems 1.1 in Subsection 6.1 by constructing appropriate approximate solutions and passing to the limit.And finally the uniqueness part of Theorems 1.1 will be proved in Subsection 6.2 by applying Osgood lemma.

    2.2 Strategy to the proof of Theorem 1.2

    As in the proof of Theorem 1.1,under the smallness conditions (7.1)and (7.2),we shall prove that the corresponding local smooth solution (θ,u)of (1.9)on [0,T*[ satisfies T*≥1,and

    Thanks to (2.8),we get,by applying the smoothing effect of etΔand et|D|s,that there exists some t1∈],1[ so that

    Starting from t1,we write the velocity equation of (1.9)as

    from which,Lemma 4.2 and under the additional smallness assumption (7.20),we infer the following.

    Lemma 2.1Under the assumptions of Proposition 7.3 and for η4given by (7.20),we have

    Lemma 2.2Under the assumptions of Proposition 7.3,for 0 <s <1,one has for any t ∈]t1,T*[,

    where psis given by (7.21).

    In view of (2.10)and the above lemmas,we obtain the a priori estimate forfor any t <T*.The detail will be presented in Proposition 7.3.

    Finally the existence part of Theorem 1.2 will be proved by constructing appropriate approximate solutions and passing to the limit in Subsection 8.1.Whereas the uniqueness part of Theorems 1.2 will be proved in Subsection 8.2 through Osgood lemma type argument.

    3 Littlewood-Paley Analysis and Lorentz Spaces

    The proofs of Theorems 1.1 and 1.2 require Littlewood-Paley decomposition.Let us briefly explain how it may be built in the case x ∈RN(see e.g.[7]).Let φ be a smooth function supported in the annulusand χ(ξ)be a smooth function supporte d in the ballsuch that

    Then for u ∈S′h(RN)(see [7,Definition 1.26]),which means u ∈S′(RN)and

    we set

    we have the formal decomposition

    Moreover,the Littlewood-Paley decomposition satisfies the property of almost orthogonality:

    We recall now the definition of homogeneous Besov spaces and Bernstein type inequalities from [7].Similar definitions in the inhomogeneous context can be found in [7].

    Definition 3.1(see [7,Definition 2.15])Let (p,r)∈[1,+∞]2,s ∈R and u ∈S′h(RN),we set

    Lemma 3.1Let B be a ball and C an annulus of RN.A constant C exists so that for any positive real number δ,any non-negative integer k,any smooth homogeneous function σ of degree m,and any couple of real numbers (a,b)with b ≥a ≥1,there hold

    Lemma 3.2(see [7,Lemma 2.4])Let C be an annulus.A positive constant C exists so that for any p ∈]1,∞[ and any couple (t,λ)of positive real numbers,we have

    We also recall Bony’s decomposition from [9]:

    where

    In order to obtain a better description of the regularizing effect of the transport-diffusion equation,we need to use Chemin-Lerner type spacesfrom [7].

    Definition 3.2Let (r,λ,p)∈[1,+∞]3and T ∈]0,+∞].We defineas the completion of C([0,T]; S(RN))by the norm

    with the usual change if r =∞.For short,we just denote this space by

    To prove Theorem 1.2,we also need to use Lorentz space Lp,q(R3).For the convenience of the readers,we recall some basic facts on Lp,q(RN)from [18,25].

    Definition 3.3(see [18,Definition 1.4.6])For a measurable function f on RN,we define its non-increasing rearrangement by

    where μ denotes the usual Lebesgue measure.For (p,q)∈[1,+∞]2,the Lorentz space Lp,q(RN)is the set of functions f such that ‖f‖Lp,q<∞,with

    We remark that Lorentz spaces can also be defined by real interpolation from Lebesgue spaces (see for instance [25,Definition 2.3]):

    where 1 ≤p0<p <p1≤∞,β satisfiesand 1 ≤q ≤∞.

    Lemma 3.3(see [25,pages 18-20])Let 1 <p <∞and 1 ≤q ≤∞,we have

    · For 1 ≤p ≤∞and 1 ≤q1≤q2≤∞,we have

    Lemma 3.4(see [14,Lemma 3.9])For 1 <p <q ≤∞,one has

    As an application of the above basic facts on Littlewood-Paley theory and Lorentz spaces,we prove the following estimates.

    Lemma 3.5Let p ≥,s ∈]-1,∞[,α ∈[0,1[ and a,b ∈S(R3).We have

    ProofThe proof of the first inequality of (3.6)basically follows from that of (26)in [14].For completeness,we present its proof here.Using Bony’s decomposition (3.5),we write

    It follows from Lemma 3.1 that

    for some (cj,r)j∈Z∈?r(Z)so that ‖(cj,r)j∈Z‖?r(Z)=1.Similar estimate holds for Tba.

    Whereas applying Lemma 3.1 and Lemma 3.3 yields

    Hence by virtue of Lemma 3.4,we obtain the first inequality of (3.6).

    Along the same line,by using Bony decomposition (3.5),we write

    It follows from the classical commutator’s estimate (see [7])that

    While it is easy to observe that

    4 Some Technical Lemmas

    In this section,we shall collect some technical lemmas which will be used throughout this paper.The first one is concerning the definition of Besov spaces with negative indices through heat semi-group.

    Proposition 4.1(see [7,Theorem 2.34])Let s be a negative real number and (p,r)∈[1,∞]2.A constant C exists such that

    The other key ingredient used in this paper is the maximal Lp(Lq)regularity for the heat kernel.

    Lemma 4.1(see [25,Lemma 7.3])The operator A defined by

    is bounded from Lp(]0,T[;Lq(R3))to Lp(]0,T[;Lq(R3))for every T ∈(0,∞] and 1 <p,q <∞.Moreover,there holds

    Lemma 4.2Let 3 <p <∞.The operator B defined byis bounded fromfor every T ∈(0,∞],and there holds

    If moreover,p ∈]3,6],one has

    ProofNote that

    Applying Young’s inequality in the space variables yields

    where 1[0,t](t′)denotes the characteristic function on [0,t],from which and Hardy-Littlewood-Sobolev inequality,we conclude the proof of (4.1).

    It remains to prove the limiting case,i.e.,(4.2).Indeed it follows by a similar derivation of(4.4)that

    Note that for 3 <p ≤6,we have

    as a result,it comes out (4.2).This completes the proof of Lemma 4.2.

    Lemma 4.3Let s ∈]0,1] and u be a smooth solenoidal vector filed on [0,T].Let θ be a smooth enough solution of

    Then one has for any t ≤T,

    (1)for all p ∈[1,∞],‖θ(t)‖Lp≤‖θ0‖Lp;

    (3)for all p ∈]1,∞[ and r ∈[1,∞],

    (4)for any α ∈]0,1[,

    ProofPart (1)follows directly from [12].While by taking the L2inner product of (4.5)with θ and using div u=0,we get

    integrating the above inequality over [0,t] yields part (2)of the lemma.

    To deal with part(3),we first apply the dyadic operatorto(4.5)and then use a standard commutator’s process to write

    Taking L2inner product of (4.6)withand using Hlder inequality,we have

    Recalling from [19,33] the following generalized Bernstein inequality that

    for some p independent constant c.We thus obtain

    from which,we infer

    and hence

    However recalling from [22,Lemma 4.3] that

    Resuming the above inequality into (4.8)and using part (1)of the lemma gives rise to part (3)of the lemma.

    Finally we deduce from Lemma 3.5 and (4.6)that for any α ∈]0,1[,

    Applying Gronwall’s inequality gives rise to

    Whereas it follows from (4.6)and the proof of part (1)in [12] that

    which ensures

    Applying Gronwall’s lemma gives

    Substituting the above inequality into (4.9)leads to part (4)of the lemma,and we complete the proof of Lemma 4.3.

    Let us complete this section by recalling the following proposition from [1].

    Proposition 4.2(see [1,Proposition 3.3],see also [7,Theorem 3.37])Let p ∈]2,∞[ andLet u,v be two solenoidal vector field which satisfy u ∈C([0,T];and

    Then there exists a constant C such that

    Moreover,there holds

    5 The a Priori Estimates Related to the System (1.4)

    In this section,we shall establish the a priori estimates which will be used to prove the global existence part of Theorem 1.1.

    5.1 The short time estimates for smooth enough solutions of (1.4)

    Proposition 5.1Let (θ,u)be a smooth enough solution of (1.4)on [0,T*[.Then under the assumption of (1.6),for all t ∈[0,T*[,we have (2.1)for any p ∈[6,8].If moreover we assume that

    for some ε1sufficiently small,one has (2.2).

    ProofLetdiv be the Leray projection operator to the divergence free vector space.In order to prove (2.1),we first apply the operator P to the u equation of (1.4)to get

    or equivalently,

    Note that 3 <p <∞,we infer from Proposition 4.1 that

    from which and Lemma 4.2,we deduce from (5.3)that for any p ∈[6,8],

    Along the same line,since p ≥6,we deduce from Proposition 4.1 and Lemma 3.1 that

    from which,(5.3)and Lemma 4.1,we infer for any p ∈[6,8] that

    By summing up (5.4)and (5.5),for ε0being sufficiently small in (1.6),we write

    Then for ε sufficiently small in (1.6),we infer

    While we deduce from (4.2)for p=6 and (5.3)that

    Note that a similar proof of (4.1)also yields

    so that we deduce from (5.3)that

    On the other hand,it follows from Lemma 4.1 that

    which together with the fact (see [13])that

    implies

    Then we deduce from (5.3)that

    Thanks to the fact that

    and ‖1-μ(θ0)‖L∞?1,we get,by summing up (5.8)and (5.9)that

    On the other hand,it follows from Proposition 4.1 that

    Yet by virtue of the real interpolation with the pairs(see [8,Theorem 6.4.5,p.160]),we get

    As a result,it comes out

    and

    Then for ε sufficiently small in (1.6),we deduce from (5.10)that

    Inserting the above inequality into (5.7)yields

    This together with (5.6)concludes the proof of (2.1).

    On the other hand,by virtue of Lemma 3.2,we deduce from (5.3)that

    Yet applying Bony’s decomposition (3.5)and standard paraproduct estimates (see [7])leads to

    Along the same line,we write

    Applying Lemma 3.1 gives

    The same estimate holds for Td(μ(θ)-1).

    Resuming the above two inequalities into (5.11)and summing up the resulting inequality over j in Z,we arrive at

    Without loss of generality,we may assume that T*>1.We define

    We claim that under the assumptions of Proposition 5.1,t*=1.Otherwise,it follows from [2,Lemma 3.1] that

    and (2.1)ensures that

    Then taking t=t*in(5.12),for ε0,ε being sufficiently small in(1.6),we deduce from(2.1)and(5.12)that

    which together with (5.1)ensures

    In particular,if we take ε,ε1so small thatwe obtain

    which contradicts with(5.13)if t*<1,and this in turn shows that t*=1 and there holds(2.2).Hence we complete the proof of Proposition 5.1.

    Proposition 5.2Under the assumptions of Proposition 5.1,one has

    ProofWe first apply the dyadic operatorto(5.2)and then taking the L2inner product of the resulting equation withu that

    from which and Lemma 3.1,we infer

    Whereas using Bony’s decomposition (3.5)and div u=0,one has

    Applying commutator’s estimate (see [7,29])gives

    The same estimate holds for div.We thus obtain

    where we used (2.2)in the last step.

    Along the same line,we get,by using paraproduct estimates in [7] that

    However applying Lemma 3.1 gives

    And it follows from [7,Theorem 3.14] and (1.6),(2.2)that

    Therefore by virtue of (2.1)for p=12,we obtain

    Resuming the above estimate and (5.16)into (5.15),we write

    Taking ε0,ε small enough in (1.6)and ε1small enough in (5.1),we obtain the second line of(5.14).

    Finally similar to (5.17),we can prove

    This concludes the proof of the proposition.

    Corollary 5.1Under the assumptions of Proposition 5.1,we can find some t0∈]0,1[ such thatu(t0)∈H1(R3),moreover,there holds (2.3).

    ProofNote that since δ ∈]0,12[,we have

    from which and (5.14),we deduce that there exists some t0∈]0,1[ so that

    This together with (5.14)concludes the proof of (2.3).

    5.2 The propagation of H1 regularity for u

    As a convention in the remaining of this section,we shall always denote t0to be the positive time determined by Corollary 5.1.

    Proposition 5.3Let (θ,u)be a smooth enough solution of (1.4)on [0,T*[.Then under the assumptions of Proposition 5.1,we have (2.4)for t ∈]t0,T*[.

    ProofThe proof of this proposition is motivated by that of [16,Theorem 1] for 2-D inhomogeneous Navier-Stokes system and that of [5,Proposition 2.1] for 3-D inhomogeneous Navier-Stokes system.In fact,since div u=0,we get,by taking L2inner product of the velocity equation of (1.4)with u,that

    integrating the above inequality over [t0,t] and using (1.5)and div u = 0,we obtain the first line of (2.4).

    Whereas by taking the L2inner product of (5.2)with ?tu,we write

    Motivated by the derivation of (29)in [16],we get,by using integration by parts,that

    Using the θ equation of (1.4)and then integration by parts,we get

    Notice that

    Hence due to div u=0,we obtain

    which together with the velocity equation of (1.4)implies that

    Integrating the above inequality over [t0,t] and using again div u = 0 and ‖u‖L6≤C‖?u‖L2,we infer

    To deal with the pressure function Π,we get,by taking space divergence to the velocity equation of (1.4),that

    and

    from which,we deduce

    We thus deduce from (5.20)that

    On the other hand,it is easy to observe that

    from which and

    we infer

    Taking ε0sufficiently small in (1.6),we obtain for p ∈[3,6],

    Substituting the above inequality for p=4 into (5.22),we obtain

    Applying Gronwall’s lemma and using (2.1)for p=6,we infer

    This concludes the proof of (2.4).

    An immediate consequence of Proposition 5.3 is the following corollary concerning the estimate offor 3 ≤p ≤8.

    Corollary 5.2Under the assumptions of Proposition 5.1,for any p ∈[3,8] and any t ∈]t0,T*[,we have

    ProofFor p ∈[3,6],by virtue of (5.25),(2.1)and (2.4),we have for any t ∈]t0,T*[,

    which together with (1.6)yields the first inequality of (5.28).For p ∈[6,8],applying H?lder’s inequality gives rise to

    from which,(2.1)and the first inequality of(5.28),we conclude the proof of the second inequality of (5.28).

    Corollary 5.3LetThen under the assumption of Proposition 5.1,one has

    ProofThe proof of this proposition basically follows from that of Proposition 5.3.We first get,by multiplying (5.19)by 〈t〉 and then by using a similar derivation of (5.26),that

    By (2.1)and Gronwall’s lemma,we thus obtain

    which together with the second inequality of (1.6)leads to (5.30).

    5.3 Large time decay estimate for u

    Lemma 5.1Under the assumptions of Proposition 5.1,we have

    ProofWe first deduce from (5.18)that

    Applying Schonbek’s strategy in [31],we split the frequency space R3into two time-dependent domains:R3= S(t)∪S(t)c,wherefor somewhich will be chosen hereafter.Then we deduce from (5.32)that

    To deal with the term on the right hand side of (5.33),we write the velocity equation of (1.4)as

    Taking Fourier transform with respect to x variables gives rise to

    so that

    We now estimate term by term above.It is easy to observe that

    and

    and

    Resuming the above estimates into (5.34)and using (2.4)yields

    Inserting the resulting inequality into (5.33),we get forgiven by (5.31).

    from which,we infer

    Hence a similar derivation of (5.35)leads to

    where in the last step,we used once again that.This ensures (5.31)forand we completes the proof of Lemma 5.1.

    Proposition 5.4Under the assumptions of Proposition 5.1,we have

    for any t ∈[t0,T*[ andgiven by (5.31).

    ProofMultiplying (5.32)by〈t〉2δ-and then integrating the resulting inequality over[t0,t],we get,by applying (5.31),that

    On the other hand,by multiplying (5.19)by,we deduce,by a similar derivation of(5.26),that

    Applying Gronwall’s lemma and using (2.1)gives rise to

    This together with(2.3)and(5.38)implies(5.37),and we completes the proof of the proposition.

    We now present the key estimate of this section.

    Proposition 5.5Under the assumptions of Theorem 1.1,for anysatisfying

    and C given by (5.31),we have (2.6).

    ProofWe first deduce from (5.25)that

    from which,(2.1)and H?lder’s inequality,we get

    On the other hand,since p satisfies (5.39),we get,by applying H?lder’s inequality,that

    This together with (2.1)concludes the proof of the proposition.

    5.4 The L1(R+; )estimate for the velocity field

    The goal of the this section is to present the a priori L1(R+;(R3))estimate for the velocity field,which is the most important ingredient used in the proof of Theorem 1.1.

    Lemma 5.2Let ∈∈]0,[,let (θ,u,?Π)be a smooth enough solution of (1.4)on [0,T*[.Then under the assumptions of Theorem 1.1,we have

    for any t ∈]t0,T*[ andgiven by (5.31).

    ProofApplyingto (5.2)and using a standard commutator’s process,we write

    Throughout this paper,we always denoteand abbreviate d(u)as d.

    Taking L2inner product of (5.41)withand using Lemma 3.1,we obtain

    Using integration by parts and [14,Lemma A.5],one has

    for some uniform positive constant c.

    Whereas it follows from Lemma 3.1 that

    Thus,by taking ε0sufficiently small in (1.6),we deduce from (5.42)that

    This gives rise to

    Then by virtue of Definition 3.2,we infer

    In what follows,we shall handle term by term the right-hand side of (5.43).Firstly applying Bony’s decomposition (3.5)yields

    Applying [29,Lemma 1] gives

    And by applying Lemma 3.1,one has

    The same estimate holds foru.Hence,thanks to the interpolation inequality

    and Corollary 5.2 implies

    The same process along with (2.6)for p=6 ensures

    Substituting (5.44)and (5.45)into (5.43)leads to (5.40),and we complete the proof of Lemma 5.2.

    With Lemma 5.2,we can prove the a priori L1(R+;(R3))estimate for u.

    Proposition 5.6Under the assumptions of Theorem 1.1,one has (2.7)for any t ∈[0,T*[.

    ProofBy virtue of (2.2),we only need to prove the estimate offor any t <T*.As a matter of fact,for any integer N and p ∈]6,8] satisfying (5.39),we deduce from Lemma 3.1 and Lemma 5.2 that

    However as

    we get,by taking

    in (5.46),that

    Therefore,in view of (2.3)and (2.6),wheneveris so small that

    we infer from (5.48)that

    which together with (2.2)leads to (2.7).This completes the proof of the proposition.

    6 The Proof of Theorem 1.1

    6.1 Existence part of Theorem 1.1

    The proof to the existence part of Theorem 1.1 basically follows the following strategy.We begin by solving an appropriate approximate problem,and then we provide uniform estimates for such approximate solutions,and finally we prove the convergence of such approximate solution sequence to a solution of (1.4)through a standard compactness argument.

    Lemma 6.1Let (θ,u,?Π)be a smooth enough solution of (1.4)on [0,T*[.Then under the assumption that

    for ε0sufficiently small,one has for any t <T*,

    ProofWe first notice that u solves

    from which,we get,by using Proposition 4.2,that

    whereas it follows from product laws in Besov spaces that

    which together with (6.1)ensures that

    On the other hand,we get,by first applying the operatorto the θ equation of (1.4)and then taking the L2inner product of the resulting equation withthat

    As div u=0,by applying (2)of (3.6)for α=we write

    which together with

    and Definition 3.2 implies that

    Summing up (6.3)and (6.4),we obtain

    Yet due to (2.1)for p=6,we have

    Then applying Gronwall’s lemma to (6.5)and using the fact that

    we conclude the proof of (6.2).

    Now let us recall the following lemma from [1].

    Lemma 6.2(see [1,Lemma 4.2])Let s ∈R and (p,r)∈[1,∞[2.Let G ∈(R3).Then there exists Gn∈H∞(R3),such that for all ε >0,there exists some n0∈N such that

    If moreover G = (G1,G2,G3)and div G = 0,we can choose Gn= (Gn1,Gn2,Gn3)so that div Gn=0.

    Proof(The Existence Part of Theorem 1.1)By virtue of Lemma 6.2,we can find θn0,un0∈H∞(R3)for n ∈N so that

    Then according to [1,Theorem 1.1],we deduce that the System (1.4)with the initial data(θn0,un0)admits a unique local in time solution (θn,un,?Πn)on [0,T*n[ verifying

    Moreover,whenever ε0,ε are small enough in (6.6),we deduce from Proposition 5.6 that there exists a positive constant C0,which depends onso that

    from which and (5.47),Lemma 6.1,we infer that (θn,un,?Πn)is uniformly bounded infor any fixed t <T*n.This implies that T*n=∞.To prove that there is a subsequence of {(θn,un,?Πn)}n∈N,which converges to a solution(θ,u,?Π)of (1.4),which satisfies(1.7),we only need to use a standard compactness argument of Lions-Aubin’s lemma.Since this argument is rather classical,we shall not present the details here.One may check similar argument from page 582 to page 583 of [1] for details.

    6.2 Uniqueness part of Theorem 1.1

    Let (θi,ui,?Πi),for i=1,2 be two solutions of (1.4),which satisfy (1.7).We denote

    Then due to (1.4),the system for (δθ,δu,δ?Π)reads

    where δF is determined by

    We first deduce from the transport equation of (6.8)that

    Yet notice that

    and

    we infer

    with

    On the other hand,it follows from (6.8)and Proposition 4.2 that

    While we deduce from (6.9)and standard product laws in Besov spaces (see [7])that

    Observing that

    Resuming the above inequality into (6.12),we obtain for t ≤t1,

    As for α ≥0 and x ∈(0,1],there holds

    therefore,by virtue of (6.11),we eventually find for t ≤t1,

    Due to

    applying Osgood lemma to (6.14)leads to

    This proves the uniqueness part of Theorem 1.1 for t ≤t1.The global in time uniqueness of solutions to (1.4),which satisfies (1.7),can be obtained by a boot-strap argument.This completes the proof of Theorem 1.1.

    7 The a Priori Estimates Related to the System (1.9)

    In this section,we shall present the a priori estimates,which we need to prove Theorem 1.2,for the System (1.9).As in Section 5,we shall first present the short time a priori estimates.

    7.1 The short time estimates for smooth enough solutions of (1.9)

    Proposition 7.1Let (θ,u)be a smooth enough solution of (1.9)on [0,T*[.We assume that ‖are so small that

    and that

    for some sufficiently small constants η1,η2.Then one has (2.8).

    ProofWe first rewrite the velocity equation of (1.9)as

    Yet it follows from [19] that

    And it is easy to observe from the proof of[2,Lemma 3.1]that this lemma still holds for smooth enough solutions of (4.5),so that we have

    Whereas according to Definition 3.1,for any fixed integer N,one has

    where we used the fact thatin the last step.Taking N in the above inequality so that

    we obtain

    Along the same line,we have

    Hence by virtue of (7.1),we deduce from (7.5)that

    To estimate ‖u‖L∞1(L3,∞),we need the following lemma,which we admit for the time being.

    Lemma 7.1Under the assumptions of Proposition 7.1,one has

    Thanks to (7.2)and (7.8),we infer from (7.7)that

    which leads to (2.8),and this completes the proof of Proposition 7.1.

    Proposition 7.1 is proved provided that we present the proof of Lemma 7.1,which we give as follows.

    Proof of Lemma 7.1In fact,it follows from (7.3)and [26,Lemma 23] that

    Applying Lemma 4.2 for p=6 yields

    Whereas due to (7.3),we get,by a similar proof of (5.8)and (5.9),that

    and

    so that thanks to (7.1)and (7.6),one has

    Since for p ≥q,there holds (see [18])

    we have

    As a result,it comes out

    In particular,if

    we infer

    Finally it is easy to observe that

    Inserting the above two inequalities into (7.9),we arrive at

    Note from (7.2)that

    from which and (7.11),we conclude the proof of (7.8).

    Proposition 7.2Let (θ,u)be a smooth enough solution of (1.9)on [0,T*[.Then under the assumptions of Proposition 7.1 and

    for some sufficiently small η3,we have

    ProofLet us denoteThen we get,by a similar derivation of (6.3),that

    from which,(7.2)and (2.8),we infer

    However it follows from part (4)of Lemma 4.3 and Proposition 7.1 that

    So that we obtain

    Yet by using interpolation and Young’s inequalities,one has for any σ >0,there exists some positive constant Cσ>0 such that

    and

    Resuming the above two inequalities into (7.14),we conclude the proof of (7.13)if there holds(7.12)and the assumptions of Proposition 7.1.

    Corollary 7.1Let (θ,u)be a smooth enough solution of (1.9)on [0,T*[.Let η2and η3be given by (7.2)and (7.12)respectively.Then under the assumptions of Proposition 7.2,we can find some t1∈],1[ so that there holds (2.9).

    ProofWe first deduce from (7.13)that

    which ensures that there exists someso that

    Next we claim that

    Indeed we first infer from (7.4)that

    Note that by using (1)of (3.6)and (7.8),we have

    and

    Hence by virtue of (7.1),we obtain

    from which,and the para-product estimates (see e.g.[7]),we infer

    However,it follows from (2.8)that

    While part (4)of Lemma 4.3 and (7.13)ensures that

    Resuming the above estimates into (7.18)and using (7.13)leads to (7.17).

    On the other hand,it follows from(7.19)that there exists some t1∈],1[so that there holds the second inequality of (2.9).And (7.17)ensures the first inequality of (2.9).This completes the proof of Corollary 7.1.

    7.2 The L1(R+; )estimate for the velocity field

    As a convention in the remaining of this subsection,we always denote t1to be the positive time determined by Corollary 7.1.

    Proposition 7.3Let 3-<s ≤1,and (θ,u)be a smooth enough solution of the System(1.9)on [0,T*[.Then under the assumptions of Proposition 7.2 and

    for some small enough constant η4,and where

    we have

    which is bounded according to (2.9).

    Proof of Proposition 7.3We admit Lemmas 2.1 and 2.2 for the time being and continue the proof of Proposition 7.3.

    As a matter of fact,due to(2.8),we only need to prove(7.22)and(7.23)forIndeed by virtue of (2.10)and Lemma 3.2,we deduce that

    Choosing N in the above inequality so that

    we obtain

    Whereas applying (1)of (3.6)and para-product estimates (see [7] for instance)leads to

    and

    where in the last step,we used (1.5)so that

    Resuming the above estimates into (7.24)and using (7.1)and (2.11)gives rise to

    It remain to handle the last term in (7.26).In the case when s = 1,p1given by (7.21)equals 6,and then we get,by using para-product estimate,that

    Yet it follows from part (3)of Lemma 4.3 and (7.31)below that

    Whereas similar to (7.24),we infer from (2.10)that

    Note that para-product estimates (see [7])along with (7.27)and (7.31)below ensures that

    Hence by virtue of (7.1)and (7.31),we infer

    We thus obtain

    Resuming the above estimate into (7.26)leads to (7.22).

    Yet it follows from part (3)of Lemma 4.3 that

    And according to Definition 3.1 and Lemma 3.3,for any integer L,we write

    Taking L so that

    in the above inequality,we obtain

    from which,we infer

    Resuming the estimates (7.29),(7.30)into (7.28)and using Lemmas 2.1 and 2.2,we arrive at

    Substituting the above estimate into (7.26)and using (2.8)gives rise to (7.23).This completes the proof of Proposition 7.3.

    We now turn to the proof of Lemmas 2.1 and 2.2.

    Proof of Lemma 2.1Note that for s ∈]0,1],psgiven by (7.21)belongs to ]3,6].Then we get,by a similar derivation of (5.4)and (5.5)that

    and

    By virtue of Part (1)and Part (2)of Lemma 4.3,we have

    and

    Therefore,whenever there hold (7.1)and (7.20),we have

    On the other hand,since psgiven by (5.39)is greater than,we get,by applying Lemma 4.2,that

    Then by virtute of (7.9),(7.10)and (7.31),we conclude

    We thus infer from (7.20)that

    which leads to (2.11).

    Proof of Lemma 2.2Thanks to (7.3),we get,by applying Proposition 4.1 and Lemma 4.2,that

    Note that for q ∈]3,∞[,we have

    Applying the above inequality for q =gives

    Resuming the above estimate into (7.32)and using (2.11)and the fact that

    we obtain the first inequality of (2.12).

    Similar to the proof of (7.32),one has

    By using interpolation inequality that

    and (7.31),we obtain

    While applying Lemma 4.2 gives rise to

    Resuming the above two estimates into(7.33)and using the assumption(7.1),we get the second inequality of (2.12).This completes the proof of Lemma 2.2.

    8 Proof of Theorem 1.2

    8.1 Existence part of Theorem 1.2

    As in Subsection 6.1,we first present the following a priori estimate.

    Lemma 8.1Let (θ,u)be a sufficiently smooth solution of (1.9)on [0,T*[.We assume that there holds (7.1).Then for any α ∈]0,1[,one has

    ProofWe first write the u equation of (1.9)as

    from which and Proposition 4.2,we deduce

    Yet by using Bony’s decomposition,we have

    and it is easy to observe that

    Hence by virtue of (7.1)and (7.25),we obtain

    which together with part (4)of Lemma 4.3 implies that

    However,it follows from interpolation and Young’s inequalities that for any σ >0,there exists some Cσ>0 so that

    We thus infer from (8.2)that

    Applying Gronwall’s lemma to the above inequality leads to (8.1).

    We now turn to the existence part of Theorem 1.2.

    The existence part of Theorem 1.2We basically follow the same line as that in Subsection 6.1.Firstly by virtue of Lemma 6.2,we can findfor n ∈N so that

    Then according to [1,Theorem 1.1],we deduce that the System (1.9)with the initial data(θn0,un0)admits a unique local in time solution (θn,un,?Πn)on [0,T*n[ verifying

    and

    Moreover,whenever η0,η are small enough in(1.10),we deduce from Proposition 7.3 that there exists a positive constant C1,which depends onso that

    from which and Lemma 8.1,we infer that (θn,un,?Πn)is uniformly bounded infor any fixed t <T*n.This implies that

    To prove that there is a subsequence of {(θn,un,?Πn)}n∈N,which converges to a solution(θ,u,?Π)of (1.9),which satisfies (1.11),we need to use a standard compactness argument of Lions-Aubin’s lemma,which we shall not present the details here.One may check similar argument from page 582 to page 583 of [1] for details.

    8.2 Uniqueness part of Theorem 1.2

    Let (θi,ui,?Πi),for i=1,2 be two solutions of (1.9)which satisfy (1.11).We denote

    Then thanks to (1.9),the system for (δθ,δu,δ?Π)reads

    where δG is determined by

    for δF given by (6.9).

    Then similar to (6.10),we first deduce from the transport diffusion equation of (8.5)and part (1)of Lemma 4.3 that

    and

    While a similar derivation of (6.13)yields for some small enough positive time t2and for t ≤t2,

    Resuming the Estimate(8.6)into(8.7)ensures that for some small enough positive time t2and t ≤t2,

    With (8.8),we can follow the same line as that in Subsection 6.2 to complete the uniqueness part of Theorem 1.2.

    AcknowledgementPart of this work was done when we were visiting Morningside Center of Mathematics,CAS,in the summer of 2013.We appreciate the hospitality and the financial support from the Center.

    97超碰精品成人国产| 91精品一卡2卡3卡4卡| 亚洲成人精品中文字幕电影| 免费黄频网站在线观看国产| 色综合色国产| 久久国产乱子免费精品| 蜜臀久久99精品久久宅男| 亚洲av日韩在线播放| 国产极品天堂在线| 草草在线视频免费看| 亚洲欧美日韩无卡精品| 亚洲欧美日韩无卡精品| 丰满乱子伦码专区| 国产人妻一区二区三区在| 国产精品一区www在线观看| 内地一区二区视频在线| 久久久久国产网址| 哪个播放器可以免费观看大片| 免费播放大片免费观看视频在线观看| 精品国产三级普通话版| 啦啦啦中文免费视频观看日本| 亚洲av免费高清在线观看| 91精品国产九色| 日韩av免费高清视频| 国产免费一级a男人的天堂| 伊人久久精品亚洲午夜| 国产成人a区在线观看| 久热久热在线精品观看| 免费观看的影片在线观看| 天天一区二区日本电影三级| 18禁裸乳无遮挡动漫免费视频 | 亚洲av成人精品一二三区| 丝袜脚勾引网站| 国内精品美女久久久久久| 麻豆成人午夜福利视频| 成人国产av品久久久| 国产一区二区亚洲精品在线观看| 国产精品熟女久久久久浪| 一级av片app| 亚洲国产精品国产精品| 成人美女网站在线观看视频| 插逼视频在线观看| 亚洲综合色惰| 王馨瑶露胸无遮挡在线观看| 午夜精品国产一区二区电影 | 色吧在线观看| 插阴视频在线观看视频| 男女那种视频在线观看| 免费观看在线日韩| 国产欧美另类精品又又久久亚洲欧美| 大香蕉97超碰在线| 国产日韩欧美在线精品| 你懂的网址亚洲精品在线观看| 在线观看一区二区三区| 女人十人毛片免费观看3o分钟| 成人鲁丝片一二三区免费| 80岁老熟妇乱子伦牲交| 亚洲av国产av综合av卡| 青春草亚洲视频在线观看| 国产精品一区www在线观看| 亚洲av成人精品一区久久| 国产黄片美女视频| 我的女老师完整版在线观看| 亚洲婷婷狠狠爱综合网| 老师上课跳d突然被开到最大视频| 日韩国内少妇激情av| 亚洲精品乱码久久久久久按摩| 777米奇影视久久| 能在线免费看毛片的网站| 久久久亚洲精品成人影院| 男女下面进入的视频免费午夜| 综合色丁香网| 欧美 日韩 精品 国产| 尾随美女入室| 欧美bdsm另类| 18禁在线无遮挡免费观看视频| 国产精品久久久久久精品古装| 日韩一区二区视频免费看| 国产午夜福利久久久久久| 日韩不卡一区二区三区视频在线| 国模一区二区三区四区视频| 国产精品麻豆人妻色哟哟久久| 黄色日韩在线| 国产欧美另类精品又又久久亚洲欧美| 成人无遮挡网站| 狠狠精品人妻久久久久久综合| 麻豆成人av视频| av在线观看视频网站免费| 丰满乱子伦码专区| 91精品国产九色| 一个人看视频在线观看www免费| 国产精品国产三级专区第一集| 男女边吃奶边做爰视频| 国产精品一二三区在线看| 欧美高清成人免费视频www| 日韩人妻高清精品专区| 日韩强制内射视频| 欧美 日韩 精品 国产| 国产成人精品福利久久| 久久久久久久久大av| 能在线免费看毛片的网站| 蜜桃久久精品国产亚洲av| 91狼人影院| 国产亚洲5aaaaa淫片| 激情 狠狠 欧美| 亚洲精品成人av观看孕妇| 人妻 亚洲 视频| 亚洲精品aⅴ在线观看| 黑人高潮一二区| 国语对白做爰xxxⅹ性视频网站| 日日啪夜夜爽| 国产免费福利视频在线观看| 亚洲最大成人av| 51国产日韩欧美| 成年女人在线观看亚洲视频 | 最近的中文字幕免费完整| 80岁老熟妇乱子伦牲交| 国产一区二区三区av在线| av在线播放精品| 久久久久久久大尺度免费视频| 日韩,欧美,国产一区二区三区| 黄片wwwwww| 午夜福利视频精品| 成人综合一区亚洲| 亚洲欧美精品专区久久| 色播亚洲综合网| 国产91av在线免费观看| 亚洲第一区二区三区不卡| 日韩三级伦理在线观看| 秋霞在线观看毛片| 国产白丝娇喘喷水9色精品| 婷婷色综合大香蕉| 一本久久精品| 人人妻人人看人人澡| 午夜福利在线在线| 午夜福利视频精品| 欧美潮喷喷水| av国产免费在线观看| 啦啦啦中文免费视频观看日本| 九九爱精品视频在线观看| av福利片在线观看| 久久久精品94久久精品| 中文字幕亚洲精品专区| 大话2 男鬼变身卡| 国产精品精品国产色婷婷| 一级av片app| 一个人看的www免费观看视频| 久久久久精品久久久久真实原创| av女优亚洲男人天堂| 亚洲av中文字字幕乱码综合| 成人国产麻豆网| av国产久精品久网站免费入址| 一级爰片在线观看| 国产精品福利在线免费观看| 亚洲精品国产av成人精品| 国产探花在线观看一区二区| 精华霜和精华液先用哪个| 亚洲欧美日韩卡通动漫| 男人添女人高潮全过程视频| 精华霜和精华液先用哪个| 欧美xxxx黑人xx丫x性爽| 真实男女啪啪啪动态图| av播播在线观看一区| 欧美激情在线99| 新久久久久国产一级毛片| 身体一侧抽搐| 观看美女的网站| 高清毛片免费看| 欧美最新免费一区二区三区| a级毛片免费高清观看在线播放| 99热全是精品| 美女脱内裤让男人舔精品视频| 久久久国产一区二区| 韩国av在线不卡| 在线精品无人区一区二区三 | 卡戴珊不雅视频在线播放| 丝袜喷水一区| 3wmmmm亚洲av在线观看| 色婷婷久久久亚洲欧美| 91精品一卡2卡3卡4卡| 日韩制服骚丝袜av| 国产精品一二三区在线看| 夫妻性生交免费视频一级片| 一级二级三级毛片免费看| 国产乱人偷精品视频| 亚洲av欧美aⅴ国产| 日本av手机在线免费观看| 涩涩av久久男人的天堂| 男女边摸边吃奶| 69av精品久久久久久| 午夜精品一区二区三区免费看| 亚洲真实伦在线观看| 国产免费一区二区三区四区乱码| 欧美日本视频| 亚洲av福利一区| 又爽又黄无遮挡网站| 国产精品久久久久久久电影| 69人妻影院| tube8黄色片| 亚洲精品久久久久久婷婷小说| 欧美成人a在线观看| 欧美潮喷喷水| 亚洲伊人久久精品综合| 国产视频内射| 全区人妻精品视频| 久久6这里有精品| 大又大粗又爽又黄少妇毛片口| 精品人妻视频免费看| 久久久久九九精品影院| 亚洲最大成人中文| 国产精品人妻久久久影院| 久久久精品欧美日韩精品| 欧美高清性xxxxhd video| 91久久精品电影网| 欧美日韩国产mv在线观看视频 | 国产亚洲5aaaaa淫片| 成人无遮挡网站| 91精品一卡2卡3卡4卡| 七月丁香在线播放| 国产一级毛片在线| 欧美精品一区二区大全| 丰满少妇做爰视频| 黄色配什么色好看| av在线天堂中文字幕| 尤物成人国产欧美一区二区三区| av国产免费在线观看| 亚洲国产精品成人久久小说| 成年版毛片免费区| 美女xxoo啪啪120秒动态图| 高清在线视频一区二区三区| 免费高清在线观看视频在线观看| 亚洲av免费高清在线观看| 亚洲国产日韩一区二区| 成人特级av手机在线观看| 直男gayav资源| 亚洲精品视频女| 天天躁日日操中文字幕| 欧美激情国产日韩精品一区| 丰满乱子伦码专区| 国产亚洲av片在线观看秒播厂| 国产亚洲5aaaaa淫片| 2021少妇久久久久久久久久久| 日韩成人伦理影院| 精品久久国产蜜桃| 午夜精品一区二区三区免费看| 1000部很黄的大片| 久久国内精品自在自线图片| 国产一区二区三区综合在线观看 | 黄色视频在线播放观看不卡| 欧美三级亚洲精品| 亚洲精品影视一区二区三区av| 日韩在线高清观看一区二区三区| 一边亲一边摸免费视频| 日本一本二区三区精品| 国产伦精品一区二区三区视频9| 久久人人爽人人片av| 免费观看av网站的网址| 六月丁香七月| 在线精品无人区一区二区三 | 亚洲婷婷狠狠爱综合网| 国产毛片a区久久久久| 欧美国产精品一级二级三级 | 亚洲国产成人一精品久久久| 亚洲欧美日韩无卡精品| 建设人人有责人人尽责人人享有的 | 男人爽女人下面视频在线观看| 岛国毛片在线播放| 视频中文字幕在线观看| kizo精华| 欧美少妇被猛烈插入视频| 老女人水多毛片| 久久精品国产自在天天线| 99热全是精品| 男人舔奶头视频| 大香蕉97超碰在线| 欧美日韩一区二区视频在线观看视频在线 | 国产精品嫩草影院av在线观看| 国产精品不卡视频一区二区| 男女国产视频网站| 免费人成在线观看视频色| 九九久久精品国产亚洲av麻豆| 亚洲四区av| 看免费成人av毛片| 久久久久久久精品精品| 国产毛片在线视频| 亚洲精品国产成人久久av| 国产日韩欧美在线精品| 国产欧美另类精品又又久久亚洲欧美| 欧美成人午夜免费资源| 亚洲精品成人久久久久久| 久久久久国产精品人妻一区二区| 永久免费av网站大全| 一区二区三区免费毛片| 免费av不卡在线播放| 色视频在线一区二区三区| 欧美97在线视频| 日日撸夜夜添| 中文在线观看免费www的网站| 最新中文字幕久久久久| 香蕉精品网在线| 国产永久视频网站| 亚洲国产精品成人久久小说| 两个人的视频大全免费| 亚洲av免费高清在线观看| 在线a可以看的网站| 亚洲精品久久午夜乱码| 精品少妇黑人巨大在线播放| av在线天堂中文字幕| 国产在线男女| 综合色av麻豆| 在线播放无遮挡| 国产精品三级大全| 日本一本二区三区精品| 街头女战士在线观看网站| 一区二区三区四区激情视频| 毛片女人毛片| 亚洲精品国产av成人精品| 波多野结衣巨乳人妻| 免费av毛片视频| 午夜爱爱视频在线播放| 亚洲性久久影院| 亚洲综合色惰| 亚洲欧美日韩另类电影网站 | 一级黄片播放器| 搞女人的毛片| 91午夜精品亚洲一区二区三区| 噜噜噜噜噜久久久久久91| 国产精品久久久久久久久免| 最近中文字幕2019免费版| 精品久久久久久电影网| 免费大片18禁| videos熟女内射| 黄色一级大片看看| 精品国产三级普通话版| 午夜爱爱视频在线播放| 九草在线视频观看| xxx大片免费视频| 美女视频免费永久观看网站| 欧美日韩在线观看h| 国产精品国产三级国产专区5o| 国产爱豆传媒在线观看| 97精品久久久久久久久久精品| 三级国产精品欧美在线观看| 久久久久精品性色| 中文天堂在线官网| 狂野欧美白嫩少妇大欣赏| 男的添女的下面高潮视频| 中文乱码字字幕精品一区二区三区| 夜夜看夜夜爽夜夜摸| 97热精品久久久久久| 五月开心婷婷网| 校园人妻丝袜中文字幕| 少妇高潮的动态图| 亚洲婷婷狠狠爱综合网| 高清午夜精品一区二区三区| 免费观看的影片在线观看| 一级毛片黄色毛片免费观看视频| 精品久久久精品久久久| 欧美日韩一区二区视频在线观看视频在线 | 欧美 日韩 精品 国产| 亚洲欧美日韩无卡精品| 欧美日韩亚洲高清精品| 在线观看一区二区三区| 中文字幕免费在线视频6| av国产免费在线观看| 国产精品久久久久久久久免| 日韩不卡一区二区三区视频在线| 精品人妻一区二区三区麻豆| 18禁裸乳无遮挡动漫免费视频 | 嫩草影院精品99| 国产成人午夜福利电影在线观看| 亚洲欧美日韩无卡精品| 美女视频免费永久观看网站| 日韩精品有码人妻一区| 日日啪夜夜爽| 国产一区亚洲一区在线观看| 男女无遮挡免费网站观看| 又粗又硬又长又爽又黄的视频| 日韩视频在线欧美| 欧美国产精品一级二级三级 | 国产高清三级在线| 久久久国产一区二区| 联通29元200g的流量卡| 欧美另类一区| 中文资源天堂在线| 国产黄片视频在线免费观看| 大香蕉久久网| 一区二区三区乱码不卡18| 美女脱内裤让男人舔精品视频| 日韩免费高清中文字幕av| 肉色欧美久久久久久久蜜桃 | 国产亚洲5aaaaa淫片| 青春草国产在线视频| 自拍偷自拍亚洲精品老妇| 99热这里只有是精品在线观看| 国产精品成人在线| 香蕉精品网在线| 精品国产乱码久久久久久小说| 亚洲久久久久久中文字幕| 草草在线视频免费看| 97在线人人人人妻| 黑人高潮一二区| 精品一区二区三卡| 国产精品爽爽va在线观看网站| 久久亚洲国产成人精品v| 免费大片18禁| 三级经典国产精品| 男插女下体视频免费在线播放| 欧美三级亚洲精品| 国产成人aa在线观看| 成人毛片a级毛片在线播放| 久久97久久精品| 青春草视频在线免费观看| 街头女战士在线观看网站| 日本欧美国产在线视频| 日本免费在线观看一区| 在线免费十八禁| 国产乱人视频| 干丝袜人妻中文字幕| 日本色播在线视频| 尤物成人国产欧美一区二区三区| 一级片'在线观看视频| 人妻 亚洲 视频| 欧美性猛交╳xxx乱大交人| 国产探花极品一区二区| 欧美日韩国产mv在线观看视频 | 又粗又硬又长又爽又黄的视频| 人妻 亚洲 视频| 欧美性感艳星| 一个人看的www免费观看视频| 王馨瑶露胸无遮挡在线观看| 三级国产精品欧美在线观看| 日韩一区二区三区影片| 日本午夜av视频| 亚洲av在线观看美女高潮| 国产精品国产三级国产av玫瑰| 夜夜爽夜夜爽视频| 91久久精品国产一区二区三区| 欧美激情国产日韩精品一区| 精品久久久久久久末码| 国产女主播在线喷水免费视频网站| 日韩 亚洲 欧美在线| 亚洲欧美精品专区久久| 直男gayav资源| 国产精品人妻久久久影院| 亚洲精品国产色婷婷电影| 亚洲欧美成人精品一区二区| 少妇裸体淫交视频免费看高清| 女人久久www免费人成看片| 一区二区三区乱码不卡18| 永久免费av网站大全| 国产精品国产三级国产专区5o| 日韩中字成人| 丝袜喷水一区| 噜噜噜噜噜久久久久久91| 亚洲欧美日韩东京热| 成人无遮挡网站| 纵有疾风起免费观看全集完整版| 久久精品久久久久久久性| 国产一区二区在线观看日韩| av在线老鸭窝| 我要看日韩黄色一级片| 国产精品一区二区三区四区免费观看| 亚洲va在线va天堂va国产| 三级国产精品欧美在线观看| 性色av一级| 草草在线视频免费看| 午夜精品一区二区三区免费看| 久久ye,这里只有精品| 免费大片黄手机在线观看| 精品久久久精品久久久| 日韩av在线免费看完整版不卡| 亚洲精品456在线播放app| 美女国产视频在线观看| 联通29元200g的流量卡| 啦啦啦在线观看免费高清www| 国产高清三级在线| 一本色道久久久久久精品综合| 久热这里只有精品99| 激情 狠狠 欧美| 亚洲国产日韩一区二区| 大香蕉97超碰在线| 午夜福利在线在线| 精品99又大又爽又粗少妇毛片| 一本色道久久久久久精品综合| 人人妻人人看人人澡| 日本爱情动作片www.在线观看| 人妻制服诱惑在线中文字幕| 99久久精品国产国产毛片| 国产高潮美女av| 日本黄色片子视频| 亚洲av.av天堂| 另类亚洲欧美激情| 毛片一级片免费看久久久久| 国产午夜福利久久久久久| 黄色视频在线播放观看不卡| 国产精品偷伦视频观看了| 亚洲久久久久久中文字幕| 国产精品人妻久久久影院| 国产精品人妻久久久久久| 少妇猛男粗大的猛烈进出视频 | 建设人人有责人人尽责人人享有的 | 亚洲精品视频女| 久久久久久久国产电影| 大片电影免费在线观看免费| 日韩成人伦理影院| 亚洲婷婷狠狠爱综合网| 天天躁夜夜躁狠狠久久av| 蜜桃亚洲精品一区二区三区| 爱豆传媒免费全集在线观看| 欧美高清成人免费视频www| 亚洲最大成人av| 91狼人影院| 精品久久久久久久末码| 亚洲av中文av极速乱| 日本三级黄在线观看| 日韩av免费高清视频| 亚洲最大成人中文| 亚洲国产最新在线播放| 汤姆久久久久久久影院中文字幕| 伊人久久国产一区二区| 久久久精品欧美日韩精品| 夫妻性生交免费视频一级片| 久久人人爽人人片av| 精品久久久久久电影网| 亚洲精华国产精华液的使用体验| 六月丁香七月| 婷婷色麻豆天堂久久| 99久久人妻综合| 精品一区在线观看国产| 色网站视频免费| 国产精品国产三级国产专区5o| 亚洲综合色惰| 91狼人影院| 午夜福利视频1000在线观看| 亚洲欧美日韩无卡精品| 国产日韩欧美在线精品| 91在线精品国自产拍蜜月| 国内精品美女久久久久久| 免费看av在线观看网站| 777米奇影视久久| 亚洲国产日韩一区二区| 国产 精品1| 搡女人真爽免费视频火全软件| 精品久久久久久久人妻蜜臀av| 嫩草影院入口| 婷婷色麻豆天堂久久| 涩涩av久久男人的天堂| 亚洲av一区综合| 欧美xxⅹ黑人| 在线亚洲精品国产二区图片欧美 | 看免费成人av毛片| 人妻少妇偷人精品九色| 精品国产三级普通话版| 国产精品人妻久久久久久| 日韩欧美精品免费久久| 日本与韩国留学比较| 最近最新中文字幕大全电影3| 噜噜噜噜噜久久久久久91| 熟女电影av网| 亚洲成人一二三区av| 又粗又硬又长又爽又黄的视频| 国产毛片a区久久久久| 亚洲人与动物交配视频| 国产精品久久久久久精品电影| 欧美高清成人免费视频www| 日韩免费高清中文字幕av| 欧美日韩在线观看h| 精品一区二区三卡| 99精国产麻豆久久婷婷| 国产精品一二三区在线看| 亚洲精品色激情综合| 国产欧美亚洲国产| 国产美女午夜福利| 久久久久久久久大av| 久久久久久久久久成人| 亚洲国产色片| 九九久久精品国产亚洲av麻豆| 欧美日本视频| 国产精品嫩草影院av在线观看| 99热国产这里只有精品6| 午夜免费观看性视频| 丰满少妇做爰视频| 国产真实伦视频高清在线观看| 日韩电影二区| 好男人视频免费观看在线| av播播在线观看一区| 联通29元200g的流量卡| 王馨瑶露胸无遮挡在线观看| 中文欧美无线码| 成年女人在线观看亚洲视频 | 波野结衣二区三区在线| 老司机影院成人| 日日啪夜夜撸| 国精品久久久久久国模美| 又黄又爽又刺激的免费视频.| 欧美另类一区| 国产毛片a区久久久久| 毛片女人毛片| 婷婷色麻豆天堂久久| 99热这里只有精品一区| 自拍欧美九色日韩亚洲蝌蚪91 | 我要看日韩黄色一级片| 伦精品一区二区三区| 天天躁日日操中文字幕| 成年女人看的毛片在线观看| 97精品久久久久久久久久精品| 国产精品成人在线| 熟妇人妻不卡中文字幕| 色吧在线观看| 久久久欧美国产精品| 国产免费视频播放在线视频| 制服丝袜香蕉在线| 亚洲欧美日韩另类电影网站 | 在线看a的网站|