Processing math: 100%
Research article

On the Cauchy problem of 3D nonhomogeneous micropolar fluids with density-dependent viscosity

  • Received: 17 June 2024 Revised: 18 July 2024 Accepted: 26 July 2024 Published: 01 August 2024
  • MSC : 35Q35, 76D03

  • In this paper, we considered the global well-posedness of strong solutions to the Cauchy problem of three-dimensional (3D) nonhomogeneous incompressible micropolar fluids with density-dependent viscosity and vacuum. Based on the energy method, some key a priori exponential decay-in-time rates of strong solutions are obtained. As a result, the existence and large-time asymptotic behavior of strong solutions in the whole space R3 are established, provided that the initial mass is sufficiently small. Note that this result is proven without any compatibility conditions.

    Citation: Mingyu Zhang. On the Cauchy problem of 3D nonhomogeneous micropolar fluids with density-dependent viscosity[J]. AIMS Mathematics, 2024, 9(9): 23313-23330. doi: 10.3934/math.20241133

    Related Papers:

    [1] Ru Bai, Tiantian Chen, Sen Liu . Global stability solution of the 2D incompressible anisotropic magneto-micropolar fluid equations. AIMS Mathematics, 2022, 7(12): 20627-20644. doi: 10.3934/math.20221131
    [2] Jianxia He, Ming Li . Existence of global solution to 3D density-dependent incompressible Navier-Stokes equations. AIMS Mathematics, 2024, 9(3): 7728-7750. doi: 10.3934/math.2024375
    [3] Abdulaziz H. Alharbi, Ahmed G. Salem . Analytical and numerical investigation of viscous fluid-filled spherical slip cavity in a spherical micropolar droplet. AIMS Mathematics, 2024, 9(6): 15097-15118. doi: 10.3934/math.2024732
    [4] Dayong Huang, Guoliang Hou . Blowup criterion for the Cauchy problem of 2D compressible viscous micropolar fluids with vacuum. AIMS Mathematics, 2024, 9(9): 25956-25965. doi: 10.3934/math.20241268
    [5] Mingyu Zhang . Regularity and uniqueness of 3D compressible magneto-micropolar fluids. AIMS Mathematics, 2024, 9(6): 14658-14680. doi: 10.3934/math.2024713
    [6] Mingyu Zhang . On the Cauchy problem of compressible Micropolar fluids subjected to Hall current. AIMS Mathematics, 2024, 9(12): 34147-34183. doi: 10.3934/math.20241627
    [7] Muhammad Naqeeb, Amjad Hussain, Ahmad Mohammed Alghamdi . Blow-up criteria for different fluid models in anisotropic Lorentz spaces. AIMS Mathematics, 2023, 8(2): 4700-4713. doi: 10.3934/math.2023232
    [8] Chenchen Lu, Lin Chen, Shaoyong Lai . Local well-posedness and blow-up criterion to a nonlinear shallow water wave equation. AIMS Mathematics, 2024, 9(1): 1199-1210. doi: 10.3934/math.2024059
    [9] Ahmed G. Salem, Turki D. Alharbi, Abdulaziz H. Alharbi, Anwar Ali Aldhafeeri . Impact of a spherical interface on a concentrical spherical droplet. AIMS Mathematics, 2024, 9(10): 28400-28420. doi: 10.3934/math.20241378
    [10] Ahmad Mohammad Alghamdi, Sadek Gala, Jae-Myoung Kim, Maria Alessandra Ragusa . The anisotropic integrability logarithmic regularity criterion to the 3D micropolar fluid equations. AIMS Mathematics, 2020, 5(1): 359-375. doi: 10.3934/math.2020024
  • In this paper, we considered the global well-posedness of strong solutions to the Cauchy problem of three-dimensional (3D) nonhomogeneous incompressible micropolar fluids with density-dependent viscosity and vacuum. Based on the energy method, some key a priori exponential decay-in-time rates of strong solutions are obtained. As a result, the existence and large-time asymptotic behavior of strong solutions in the whole space R3 are established, provided that the initial mass is sufficiently small. Note that this result is proven without any compatibility conditions.



    The nonhomogeneous incompressible micropolar fluids with density-dependent viscosity ([12,20]) in R3 read as follows:

    {ρt+div(ρu)=0,(ρu)t+div(ρuu)+P(ρ)=div((μ(ρ)+ζ)u)+2ζrotw,(ρw)t+div(ρuw)+4ζw=μΔw+(μ+λ)divw+2ζrotu,divu=0, (1.1)

    where ρ,u=(u1,u2,u3), w=(w1,w2,w3), and P denote the fluid density, velocity, micro-rotational velocity, and pressure, respectively. The viscosity coefficient μ(ρ) satisfies

    μC1[0,),μ(ρ)α>0, (1.2)

    for some positive constant α, while the constants μ and λ are the angular viscosities satisfying μ>0 and μ+λ0, and the constant ζ>0 denotes the dynamic micro-rotation viscosity.

    In this paper, we consider the Cauchy problem of (1.1)–(1.2) with the far-field behavior

    (ρ,u,w)(x,t)(0,0,0)as|x| (1.3)

    and the initial conditions

    (ρ,ρu,ρw)(x,0)=(ρ0,ρ0u0,ρ0w0)(x)withxR3. (1.4)

    A micropolar fluid system is the study of fluids that exhibit micro-rotational effects and micro-rotational inertia and can be viewed as non-Newtonian fluids. It can be used to describe many phenomena that appear in a large number of complex fluids, such as suspensions, animal blood, and liquid crystals. The micropolar fluid system reduces to the Navier-Stokes equations when there is no microstructure (ζ=0 and w=0) and has been discussed by many mathematicians (see [1,2,5,6,7,8,9,10,13,17,19], and references therein).

    When it comes to the case that ζ0 and w0, there have been substantial developments on the global regularity problem concerning nonhomogeneous micropolar fluids (1.1) with constant viscosity μ. When the initial density is strictly away from vacuum, Braz e Silva and his cooperators [4] investigated the global existence and uniqueness of solutions for 3D nonhomogeneous asymmetric fluids by using an approach and Lagrangian coordinates under suitable initial conditions. Qian-Chen-Zhang [15] studied the global existence of weak and strong solutions to 3D nonhomogeneous incompressible asymmetric fluid equations. For initial velocities sufficiently small in the critical Besov space, global Fujita-Kato type solutions with initial density in the bounded function space and that have a positive lower bound are obtained, and this result extends the classical one on the life-span by Leray. Subsequently, Qian-He-Zhang [16] investigated the global existence and uniqueness of the solutions for the 2D inhomogeneous incompressible asymmetric fluids, with the initial (angular) velocity being located in sub-critical Sobolev spaces Hs(R2)(0<s<1) and the initial density being bounded from above and below by some positive constants. In particular, the uniqueness of the solution in [16] is also obtained without any more regularity assumptions on the initial density. When the initial density contains a vacuum state, Braz e Silva and Santos [3] established the existence of global in-time weak solutions for the equations of asymmetric incompressible fluids with variable density. Zhang-Zhu [18] proved the global existence of strong solutions under the condition of the following compatibility:

    {(μ+ζ)Δu0+P(ρ0)2ζrotw0=ρ1/20g1,μΔw0(μ+λ)divw0+4ζw02ζrotu0=ρ1/20g2,

    for some (P(ρ0),g1,g2)H1×L2×L2.

    When we consider the case of μ=μ(ρ), Liu-Zhong [12] showed that the initial boundary value problem of 2D nonhomogeneous micropolar fluids with density-dependent viscosity has a global and unique strong solution under the assumption of the smallness of μ(ρ0)Lq. Qian-Qu [14] investigated the 3D inhomogeneous incompressible asymmetric fluids system and proved local well-posedness for initial velocity in the critical Besov space ˙B3/pp,1 for 1<p<6 and initial density ρ0 satisfying that ρ01 is in the critical Besov space and that ρ0 is bounded away from zero. Zhong [20] also considered the same model in 3D cases and established the global existence and uniqueness of strong solutions, provided that the initial energy is sufficiently small. It is worth noting that there is no need to impose some compatibility condition on the initial data.

    It should be noted that although the large initial velocity is allowed in [12], it excludes large oscillations of the initial density. A natural question arises: where can we establish the global strong solutions to the 3D Cauchy problem of (1.1)–(1.4) not only with large initial velocity but also allowing large oscillations of the initial density? In fact, this is the main aim of this paper.

    Before stating the main results, we set

    fdxR3fdx,(f,g)LpfLp+gLp.

    For 1r and β>0, we denote the standard homogeneous and inhomogeneous Sobolev spaces as follows:

    {Lr=Lr(R3),Wk,r=Wk,r(R3),Hk=Wk,2,Dk,r=Dk,r(R3)={vL1loc(R3)|kvLr(R3)},Dk=Dk,2,D1={vL6(R3)|vL2(R3)},C0,σ={fC0|divf=0},D10,σ=¯C0,σclosureinthenormofD1.

    The main result can be stated as follows:

    Theorem 1.1. For constant ˉρ and any given number q(3,6), assume that the initial data (ρ0,u0,w0) satisfies

    0ρ0ˉρ,ρ0L1H1,μ(ρ0)Lq,u0D10,σ,w0H10. (1.5)

    Then there exists a positive constant ε, depending only on ζ,μ,λ,ˉρ,q,α, βsup[0,ˉρ]μ(ρ), μ(ρ0)Lq, u0L2 and w0L2 such that if

    m0ρ0L1ε, (1.6)

    then the problem (1.1)–(1.4) possesses a unique global strong solution (ρ,u,w) in R3×(0,) satisfying that for any 0<t<T< and s(3,q),

    {0ρˉρ,ρL(0,;L1H1)C([0,);L1H1),(tu,tw)L(0,;H1)L2(0,;W1,s),tPL(0,;L2)L2(0,;Ls),(tu,tw)C([0,);H1),(ρu,ρw)C([0,);L2),(tρ1/2ut,tρ1/2wt)L(0,;L2),(tut,twt)L2(0,;L2). (1.7)

    Moreover,

    sup0tTμ(ρ)Lq2μ(ρ0)Lq,

    and for any t1,

    ut(,t)2H1+wt(,t)2H1+P(,t)2L2+ρ1/2ut(,t)2L2+ρ1/2wt(,t)2L2Ceσt,

    where σ=3σ1(π2)4/3ρ01L3/2 with σ1=min{α,μ}.

    Remark 1.1. It is worth noting that Theorem 1.1 holds for arbitrarily large initial velocity with a smallness only on the initial mass, which generalizes the result of [14], where they need the smallness assumption on u0˙B3/pp,1 with 1<p<6.

    The rest of this paper is organized as follows: In Section 2, we collect some elementary facts and inequalities that will be used later. Section 3 is devoted to a priori estimates, and Theorem 1.1 is also proved in Section 3.

    In this section, we list some auxiliary lemmas that will be used later. First of all, we start with the local existence of strong solutions that can be obtained from similar arguments as used in [13,20], and we omit the details.

    Lemma 2.1. Assume that (ρ0,u0,w0) satisfies (1.5). Then, there exists a small positive time T0 such that the problem (1.1)–(1.4) has a unique strong solution (ρ,u,w) on R3×(0,T0].

    Next, the following well-known Gagliardo-Nirenberg inequality will be used more frequently later (see [11]).

    Lemma 2.2. Let p[2,3s3s] for s[2,3), or p[2,] for s=3, and let q(1,),r(3,). There exists some generic constant C>0 may depend on s and r such that for fL2D1,s0 and gLqD1,r0, we have

    fpLpCfp3s(p2)/(5s6)L2f3s(p2)/(5s6)Ls, (2.1)

    and

    gLCgq(r3)/(3r+q(r3))Lqg3r/(3r+q(r3))Lr. (2.2)

    The following regularity results on the Stokes equations will be used for the derivations of higher-order a priori estimates (see [10]).

    Lemma 2.3. For constants q(3,6),α>0 and β>0, in addition to (1.2), assume that μ(ρ) satisfies

    μ(ρ)Lq,0<αμ(ρ)β<.

    Then, if GLr with r(2,q), there exists some positive constant C depending only on α,β,r and q such that the unique weak solution (u,P)D10,σ×L2 to the following problem

    {div((μ(ρ)+ζ)u)+P=G,xR3,divu=0,xR3,u(x)0,|x|,

    satisfies

    2uL2+PL2CGL2(1+μ(ρ)q/(q3)Lq), (2.3)

    and

    2uLr+PLrCGLr(1+μ(ρ)q(5r6)/2r(q3)Lq). (2.4)

    In this section, we will establish some necessary a priori bounds of local strong solutions (ρ,u,w) to the Cauchy problem (1.1)–(1.4), whose existence is guaranteed by Lemma 2.1. Thus, let T>0 be a fixed time and (ρ,u,w) be the smooth solution to (1.1)–(1.4) on R3×(0,T] with smooth initial data (ρ0,u0,w0) satisfying (1.5).

    We have the following key a priori estimates on (ρ,u,w).

    Proposition 3.1. There exists some positive constant ε0 depending only on q,ζ,ˉρ,α,β,μ, ρ0L3/2, u0L2,w0L2, and μ(ρ0)Lq such that if (ρ,u,w) is a smooth solution of (1.1)–(1.4) on R3×(0,T] satisfying

    supt[0,T]μ(ρ)Lq4μ(ρ0)Lq, (3.1)

    the following estimates hold

    supt[0,T]μ(ρ)Lq2μ(ρ0)Lq, (3.2)

    provided that

    m0ε0. (3.3)

    The proof of Proposition 3.1 consists of Lemmas 3.1–3.4 and is to be completed by the end of this section. Throughout this section, for simplicity, we denote by C or Ci(i=1,2,) the generic positive constants, which may depend on q,ζ,ˉρ,α,β,μ,λ,ρ0L3/2,u0L2 and w0L2, but are independent of time T>0 and m0.

    We begin with the following estimates:

    Lemma 3.1. Let (ρ,u,w) be a smooth solution of (1.1)–(1.4) on R3×(0,T] satisfying (1.5). Then one has

    supt[0,T]ρLpCρ0Lp,forp[1,], (3.4)
    supt[0,T](ρ1/2u2L2+ρ1/2w2L2)+T0(αu2L2+μw2L2)dt+T0ζrotu2w2L2dtCm2/30 (3.5)

    and

    supt[0,T]eσt(ρ1/2u+ρ1/2w2L2)+T0eσt(αu2L2+μw2L2)dt+T0eσtζrotu2w2L2dtCm2/30, (3.6)

    where σ=3σ1(π2)4/3ρ01L3/2 with σ1=min{α,μ}.

    Proof. First, Eq (3.4) can be shown by standard arguments ([8]).

    Next, in order to prove (3.5), we multiply (1.1)2 and (1.1)3 by u and w, respectively, and in integrating the resulting equations by parts over R3, we get after adding them together and using (1.1)4 that

    12ddt(ρ1/2u,ρ1/2w)2L2+μ(ρ)1/2u2L2+μw2L2+(μ+λ)divw2L2+ζrotu2w2L2=0. (3.7)

    Integrating (3.7) over [0,T] gives

    supt[0,T](ρ1/2u,ρ1/2w)2L2+T0(αu2L2+μw2L2+ζrotu2w2L2)dtCρ0L3/2(u02L6+w02L6), (3.8)

    which, together with (2.1) and (3.4), yields (3.5).

    Finally, we notice from (2.1), (3.4), and (3.8) and Hölder's inequality that

    ρ1/2u2L2+ρ1/2w2L2ρL3/2(u2L6+w2L6)13(2π)4/3ρ0L3/2(u2L2+w2L2), (3.9)

    where we have used the following fact:

    f2L613(2π)4/3f2L2,foranyfD1.

    Combining (3.7) with (3.9), one has

    12ddt(ρ1/2u2L2+ρ1/2w2L2)+σ(ρ1/2u2L2+ρ1/2w2L2)0,

    where σ=3σ1(π2)4/3ρ01L3/2 with σ1=min{α,μ}. By using Gronwall's inequality, one has

    ρ1/2u2L2+ρ1/2w2L2e2σt(ρ1/20u02L2+ρ1/20w02L2). (3.10)

    Multiplying (3.7) by eσt and using (3.10) show that

    ddt[eσt(ρ1/2u,ρ1/2w)2L2]+2eσt(αu2L2+μw2L2+ζrotu2w2L2)σeσt(ρ1/2u2L2+ρ1/2w2L2)σeσt(ρ1/20u02L2+ρ1/20w02L2). (3.11)

    Integrating the above inequality over [0,T] leads to (3.6).

    Remark 3.1. Evidently, we can infer from (1.2) and (3.4) that

    0<αμ(ρ)βmax0ρˉρμ(ρ)<. (3.12)

    Lemma 3.2. Let the condition of (3.1) be in force, then there exists some positive constant ε1, depending only on q,ζ,ˉρ,α,β,μ,λ,μ(ρ0)Lq,u0L2 and w0L2 such that if

    supt[0,T](αu2L2+μw2L2+(μ+λ)divw2L2+ζrotu2w2L2)+T0(ρ1/2ut2L2+ρ1/2wt2L2)dt4(βu02L2+μw02L2+(μ+λ)divw02L2+ζrotu02w02L2), (3.13)

    then

    supt[0,T](αu2L2+μw2L2+(μ+λ)divw2L2+ζrotu2w2L2)+T0(ρ1/2ut2L2+ρ1/2wt2L2)dt2(βu02L2+μw02L2+(μ+λ)divw02L2+ζrotu02w02L2), (3.14)

    provided

    m0ε1.

    Moreover, for i=1,2,3 and σ, as in Lemma 3.1, one has

    supt[0,T][ti(u2L2+w2L2)]+T0ti(ρ1/2ut2L2+ρ1/2wt2L2)dtCm2/30, (3.15)

    and

    supt[0,T][eσt(u2L2+w2L2)]+T0eσt(ρ1/2ut2L2+ρ1/2wt2L2)dtCm2/30. (3.16)

    Proof. The Eq (1.1)1 can be written as

    [μ(ρ)]t+uμ(ρ)=0. (3.17)

    Multiplying (1.1)2 by ut and integrating the resulting equation over R3, one can deduce from (3.17) that

    12ddt(μ1/2(ρ)u2L2+ζrotu2L2)+ρ1/2ut2L2=2ζ(rotuwdx)t2ζrotuwtdxρuuutdx12uμ(ρ)|u|2dx, (3.18)

    where we have used the fact that Δu+rot(rotu)=divu=0.

    Multiplying (1.1)3 by wt and integrating by parts over R3 show that

    12ddt(μw2L2+(μ+λ)divw2L2+4ζw2L2)+ρ1/2wt2L2=2ζrotuwtdxρuwwtdx, (3.19)

    Combining (3.18) with (3.19), yields

    ddt(μ1/2(ρ)u2L2+μw2L2+(μ+λ)divw2L2+ζrotu2w2L2)+2ρ1/2ut2L2+2ρ1/2wt2L2=2ρuuutdx2ρuwwtdxuμ(ρ)|u|2dx=3i=1Ii. (3.20)

    Now, we estimate Ii(i=1,2,3.) as follows: Hölder's inequality, together with (2.1), (2.2), and (3.4), gives

    I112ρ1/2ut2L2+2ρLu2L6u2L312ρ1/2ut2L2+C(ˉρ)u3L22uL2. (3.21)

    Similarly,

    I212ρ1/2wt2L2+2ρLu2Lw2L212ρ1/2wt2L2+C(ˉρ)w2L2uL22uL2, (3.22)

    and

    I3C(β)uL62uL2uL3C(β)u3/2L22u3/2L2. (3.23)

    Putting (3.21)–(3.23) into (3.20), one has

    ddt(μ1/2(ρ)u2L2+μw2L2+(μ+λ)divw2L2+ζrotu2w2L2)+32(ρ1/2ut2L2+ρ1/2wt2L2)C(u3L2+w2L2uL2)2uL2+Cu3/2L22u3/2L2. (3.24)

    We know from (1.1)2 that (u,P) satisfies the following system:

    {div[(μ(ρ)+ζ)u]+P=ρutρuu+2ζrotw,xR3,divu=0,xR3,u(x)0,|x|,

    Taking G=ρutρuu+2ζrotw and r=2 in (2.3), we deduce from (3.1), (3.4), and Cauchy-Schwarz's inequality that

    2uL2+PL2Cρutρuu+2ζrotwL2(1+μ(ρ)q/(q3)Lq)Cρ1/2Lρ1/2utL2+CρLuL6uL3+CwL2Cρ1/2utL2+Cu3/2L22u1/2L2+CwL2122uL2+Cρ1/2utL2+Cu3L2+CwL2,

    which gives

    2uL2+PL2Cρ1/2utL2+Cu3L2+CwL2. (3.25)

    Taking (3.25) into (3.24), one deduces from (3.13) that

    ddt(μ1/2(ρ)u2L2+μw2L2+(μ+λ)divw2L2+ζrotu2w2L2)+(ρ1/2ut2L2+ρ1/2wt2L2)C(u2L2+w2L2). (3.26)

    Integrating (3.26) over [0,T], then using (3.5) and (3.12), shows that

    supt[0,T](αu2L2+μw2L2+(μ+λ)divw2L2+ζrotu2w2L2)+T0(ρ1/2ut2L2+ρ1/2wt2L2)dtC1m2/30+(βu02L2+μw02L2+(μ+λ)divw02L2+ζrotu02w02L2), (3.27)

    where the positive constant C1 depends only on q,ζ,ˉρ,α,β,μ,λ,μ(ρ0)Lq,u0L2, and w0L2, but is independent of time T>0 and m0. Taking

    m0ε1(βu02L2+μw02L2+(μ+λ)divw02L2+ζrotu02w02L2C1)3/2.

    Thus, we obtain (3.14).

    Multiplying (3.26) by t gives that

    ddt[t(μ1/2(ρ)u2L2+μu2L2+(μ+λ)divw2L2+ζrotu2w2L2)]+t(ρ1/2ut2L2+ρ1/2wt2L2)C(u2L2+w2L2+rotu2w2L2)+Ct(u2L2+w2L2). (3.28)

    Integrating (3.28) over [0,T], we obtain from (3.5) and (3.6) that

    supt[0,T][t(u2L2+w2L2)]+T0t(ρ1/2ut2L2+ρ1/2wt2L2)dtCT0(u2L2+w2L2+rotu2w2L2)dt+CT0t(u2L2+w2L2)dtCm2/30+Csupt[0,T](teσt)T0eσt(u2L2+w2L2)dtCm2/30,

    which implies that (3.15) holds for i=1. For i=2,3, we can take a similar approach to obtain the results of (3.15).

    Next, multiplying (3.26) by eσt, one has

    ddt[eσt(μ1/2(ρ)u2L2+μu2L2+(μ+λ)divw2L2+ζrotu2w2L2)]+eσt(ρ1/2ut2L2+ρ1/2wt2L2)Ceσt(u2L2+w2L2+rotu2w2L2). (3.29)

    Integrating (3.29) over [0,T], one can deduce from (3.6) that

    supt[0,T][eσt(u2L2+w2L2)]+T0eσt(ρ1/2ut2L2+ρ1/2wt2L2)dtCT0eσt(u2L2+w2L2+rotu2w2L2)dtCm2/30.

    Thus, we obtain (3.16).

    Lemma 3.3. Under the condition of (3.1), then for i=1,2,3,

    supt[0,T][ti(ρ1/2ut2L2+ρ1/2wt2L2)]+T0ti(ut2L2+wt2L2)dtCm2/30. (3.30)

    Moreover, for σ as in Lemma 3.1 and δ(t)min{1,t},

    supt[δ(T),T][eσt(ρ1/2ut2L2+ρ1/2wt2L2)]+Tδ(T)eσt(ut2L2+wt2L2)dtC. (3.31)

    Proof. Operating t to (1.1)2 and (1.1)3, respectively, we infer from (1.1)1 that

    ρutt+ρuutdiv((μ(ρ)+ζ)u)t+Pt=ρutu+(uρ)(ut+uu)+2ζrotwt, (3.32)

    and

    ρwtt+ρuwtμΔwt(μ+λ)divwt+4ζwt=ρutw+(uρ)(wt+uw)+2ζrotut. (3.33)

    Multiplying (3.32) and (3.33) by ut and wt, respectively, and then integrating by parts on R3, we deduce from (3.17) that

    12ddt(ρ1/2ut2L2+ρ1/2wt2L2)+(μ(ρ)+ζ)1/2ut2L2+μwt2L2+4ζwt2L2+(μ+λ)divwt2L2=[ρutu+(uρ)(ut+uu)]utdx+[ρutw+(uρ)(wt+uw)]wtdx+2ζ(rotwtut+rotutwt)dx+uμ(ρ)uutdx=4i=1Ni. (3.34)

    Now, we estimate Ni(i=1,2,3,4) as follows: Thanks to divu=0, we deduce from (2.1), (2.2), (3.4), (3.14) and Cauchy-Schwarz's inequality that

    N1Cρ1/2Lρ1/2utL3utL2uL6+CρL2uL2utL6u2L6+CρLutL2uL6u2L6+CρLutL6uL6uL2uL6+ρ1/2ut2L4uL2Cρ3/4Lρ1/2ut1/2L2uL2ut3/2L2+CρLutL22uL2u2L2α6ut2L2+Cρ1/2ut2L2+C2u2L2, (3.35)

    and

    N2Cρ1/2Lρ1/2wtL3wtL2uL6+CρL2wL2wtL6u2L6+CρLwtL2wL6u2L6+CρLwtL6wL6uL2uL6+ρ1/2wtL3wL2utL6ρ1/2LCρ3/4Lρ1/2wt1/2L2uL2wt3/2L2+CρLwtL22wL2u2L2+Cρ3/4Lρ1/2wt1/2L2wL2wt1/2L2utL2α6ut2L2+μ2wt2L2+Cρ1/2wt2L2+C2w2L2. (3.36)

    Cauchy-Schwarz's inequality gives

    N34ζwt2L2+ζut2L2. (3.37)

    The inequalities (2.2) and (3.1), together with Cauchy-Schwarz's inequality, show that

    N4CuLμ(ρ)LquL2q/(q2)utL2Cu13q+12L22u3q+12L2utL2α6ut2L2+Cu2L2+C2u4L2. (3.38)

    Putting (3.35)–(3.38) into (3.34) and using (3.25) yields

    ddt(ρ1/2ut2L2+ρ1/2wt2L2)+ut2L2+wt2L2C(ρ1/2ut2L2+ρ1/2wt2L2)+C(2u2L2+2w2L2)+C2u4L2C(ρ1/2ut2L2+ρ1/2wt2L2)+C(u2L2+w2L2)+Cρ1/2ut2L2ρ1/2ut2L2+C2w2L2. (3.39)

    To deal with the last term on the right side of (3.39), we first multiply (1.1)3 by w and then integrate the resulting equation by parts on R3 to get that

    μw2L2+(μ+λ)divw2L2+4ζw2L2wL2ρwt+ρuw2ζrotuL2ζw2L2+C(ρwt2L2+ρuw2L2+u2L2),

    which implies

    wL2+wL2C(ρwtL2+ρuwL2+uL2). (3.40)

    On the other hand, one can deduce from (1.1)3 that

    2wL2C(ρwt+ρuw2ζrotu+4ζwL2)CwL2+C(ρwtL2+ρuwL2+uL2),

    which, together with (3.40), (2.1), and (3.4), yields

    wH2C(ρwtL2+ρuwL2+uL2)Cρ1/2wtL2ρ1/2L+CρLuL6wL3+CuL212wH2+Cρ1/2wtL2+Cu2L2wL2+CuL2,

    thus

    wH2Cρ1/2wtL2+Cu2L2wL2+CuL2. (3.41)

    Taking (3.41) into (3.39), we infer from (3.14) that

    ddt(ρ1/2ut2L2+ρ1/2wt2L2)+ut2L2+wt2L2C(ρ1/2ut2L2+ρ1/2wt2L2)+C(u2L2+w2L2)+Cρ1/2ut2L2ρ1/2ut2L2. (3.42)

    Multiplying (3.42) by t, one has

    ddt[t(ρ1/2ut2L2+ρ1/2wt2L2)]+t(ut2L2+wt2L2)Cρ1/2ut2L2[t(ρ1/2ut2L2+ρ1/2wt2L2)]+Ct(u2L2+w2L2+ρ1/2ut2L2+ρ1/2wt2L2)+C(ρ1/2ut2L2+ρ1/2wt2L2),

    which, together with (3.14), (3.15) and Gronwall's inequality, gives that (3.30) holds for i=1. Using the same methods, we can show that (3.30) for i=2,3.

    Next, in order to prove (3.31), we first multiply (3.42) by eσt, then

    ddt[eσt(ρ1/2ut2L2+ρ1/2wt2L2)]+eσt(ut2L2+wt2L2)Cρ1/2ut2L2[eσt(ρ1/2ut2L2+ρ1/2wt2L2)]+Ceσt(u2L2+w2L2+ρ1/2ut2L2+ρ1/2wt2L2),

    which, combining (3.6), (3.14), and (3.16) with Gronwall's inequality, gives that (3.31) holds.

    Lemma 3.4. Under the condition of (3.1), then

    T0uLdtCm1/2r0, (3.43)

    where 3<r<q with q(3,6).

    Proof. It follows from (1.1)2, (2.1), (2.2), (2.4), and (3.1) that for any r(3,min{6,q})

    2uLr+PLrCρutρuu+2ζrotwLr(1+μ(ρ)q(5r6)/2r(q3)Lq)Cρ1/2Lρ1/2ut(6r)/2rL2ρ1/2ut(3r6)/2rL6+CρL6r/(6r)uLuL6+Cw(6r)/2rL2w(3r6)/2rL6Cρ1/2ut(6r)/2rL2ut(3r6)/2rL2+Cm(6r)/6r0u1/2L22u3/2L2+Cw(6r)/2rL22w(3r6)/2rL2, (3.44)

    thus

    uLCu(2r6)/(5r6)L22u3r/(5r6)LrCuL2+C2uLrCuL2+Cρ1/2ut(6r)/2rL2ut(3r6)/2rL2+Cm(6r)/6r0u1/2L22u3/2L2+Cw(6r)/2rL22w(3r6)/2rL2, (3.45)

    The inequality (3.6), together with Hölder's inequality, shows that

    T0uL2dt=T0eσt/2uL2eσt/2dtC(T0eσtu2L2dt)1/2(T0eσtdt)1/2Cm1/30. (3.46)

    For any T(0,1], we can show from (3.30) and Hölder's inequality that

    T0ρ1/2ut(6r)/2rL2ut(3r6)/2rL2dtCsupt[0,T](tρ1/2ut2L2)(6r)/4r(T0tut2L2dt)(3r6)/4r×(T0t2r/(r+6)dt)(r+6)/4rCm1/30. (3.47)

    For T>1, one can deduce from (3.30), (3.31), and Hölder's inequality that

    T1ρ1/2ut(6r)/2rL2ut(3r6)/2rL2dtCsupt[0,T](eσtρ1/2ut2L2)(6r)/4r(T0t3ut2L2dt)(3r6)/4r×(T1t9(r2)/(r+6)dt)(r+6)/4rCm(3r6)/6r0. (3.48)

    Similarly,

    T0w(6r)/2rL22w(3r6)/2rL2dtCm(3r6)/6r0+Cm1/30. (3.49)

    It follows from (3.5), (3.13), and (3.25) that

    T0m(6r)/6r0u1/2L22u3/2L2dtCm(6r)/6r0(T0u2L2dt)1/4(T0(ρ1/2ut2L2+u6L2+w2L2)dt)3/4Cm1/r0+Cm1/2+1/r0. (3.50)

    Due to 3<r<q,q(3,6), one has

    12r<3r66r<13.

    Integrating (3.45) over t[0,T] and using (3.46)–(3.50), we can obtain (3.43).

    With Lemmas 3.1–3.4 at hand, we are in a position to prove Proposition 3.1.

    Proof of Proposition 3.1. It follows from (3.17) that

    [μ(ρ)]t+u2μ(ρ)+uμ(ρ)=0.

    Multiplying the above equation by q|μ(ρ)|q2μ(ρ) and integrating the resulting equations on R3, we can obtain

    ddtμ(ρ)qLq=qu2μ(ρ)|μ(ρ)|q2μ(ρ)dxquμ(ρ)|μ(ρ)|q2μ(ρ)dx.

    Due to divu=0, then

    qu2μ(ρ)|μ(ρ)|q2μ(ρ)dx=|μ(ρ)|qdivudx=0.

    Thus

    (μ(ρ)qLq)tquLμ(ρ)qLq,

    which implies that

    (μ(ρ)Lq)tuLμ(ρ)Lq,

    which, together with (3.43) and Gronwall's inequality, shows that

    supt[0,T]μ(ρ)Lqexp{T0uLdt}μ(ρ0)Lqexp{C2m1/2r0}μ(ρ0)Lq2μ(ρ0)Lq, (3.51)

    provided

    m0ε0min{ε1,(ln2C2)2r}.

    Thus, we complete the proof of Proposition 3.1.

    Lemma 3.5. Under the condition of (3.1), then

    supt[0,T](ρL2+ρtL3/2)C. (3.52)

    Proof. Similar to the method of (3.51), we can deduce from (1.1)1 that

    supt[0,T]ρL2C. (3.53)

    Hölder's inequality, together with (1.1)1, (3.14), and (3.53), yields

    ρtL3/2=uρL3/2ρL2uL6CρL2uL2C.

    Thus, we complete the proof of Lemma 3.5.

    Lemma 3.6. Under the condition of (3.1), then for 3<r<q with q(3,6), the following estimates hold:

    supt[0,T]t((u,w)2H1+P2L2)+T0t((u,w)2W1,r+P2Lr)dtC, (3.54)

    and

    supt[δ(T),T]eσt((u,w)2H1+P2L2)C, (3.55)

    for σ as in Lemma 3.1 and δ(t) as in Lemma 3.3.

    Proof. It follows from (3.14), (3.15), (3.25), (3.30), and (3.41) that

    u2H1+w2H1+P2L2C(ρ1/2ut2L2+ρ1/2wt2L2)+C(u2L2+w2L2),

    thus

    supt[0,T](tu2H1+w2H1+P2L2)Ct(ρ1/2ut2L2+ρ1/2wt2L2)+Ct(u2L2+w2L2)C. (3.56)

    Hence, we can use the same methods to obtain (3.55).

    By virtue of r(3,q) with q(3,6), one can obtain from (3.44) that

    u2W1,r+P2Lr+w2W1,rC(u2Lr+w2Lr)+C(2u2Lr+P2Lr)+C2w2LrC(2u2L2+2w2L2+u2L2+w2L2)+C(2u2Lr+P2Lr)+C2w2LrC(ρ1/2ut2L2+ρ1/2wt2L2+u2L2+w2L2+w2L2)+C(ρ1/2ut(6r)/rL2ut(3r6)/rL2+uL22u3L2+wL22w3L2)+C(ρ1/2wt(6r)/rL2wt(3r6)/rL2+uL22uL22w2L2)C(ρ1/2ut2L2+ρ1/2wt2L2)+C(ut2L2+wt2L2)+C(uL22u3L2+wL22w3L2)+C(u2L2+w2L2+w2L2),

    which, together with (3.5), (3.15), (3.25), (3.30), (3.41), and (3.56), yields

    T0t((u,w)2W1,r+P2Lr)dtC.

    Thus, we complete the proof of Lemma 3.6.

    Proof of Theorem 1.1. Similar to the standard arguments in [10,18,20], with all the a priori estimates established in Sections 3 at hand, we can immediately obtain our main results.

    The author declares she has not used Artificial Intelligence (AI) tools in the creation of this article.

    We sincerely thank the Associate Editor and the anonymous referees for their carefully reading and helpful suggestions that led to the improvement of the paper.

    M.Y. Zhang was partially supported by the Natural Science Foundation of Shandong Province of China (Grant No. ZR2021QA049) and the Science and Technology Project of Weifang (2022GX006).

    The author declares that she has no conflict of interest.



    [1] H. Abidi, G. L. Gui, P. Zhang, On the decay and stability to global solutions of the 3D inhomogeneous Navier-Stokes equations, Commun. Pure Appl. Math., 64 (2011), 832–881. https://doi.org/10.1002/cpa.20351 doi: 10.1002/cpa.20351
    [2] H. Abidi, P. Zhang, On the global well-posedness of 2-D inhomogeneous incompressible Navier–Stokes system with variable viscous coefficient, J. Differ. Equations, 259 (2015), 3755–3802. https://doi.org/10.1016/j.jde.2015.05.002 doi: 10.1016/j.jde.2015.05.002
    [3] P. Braz e Silva, E. G. Santos, Global weak solutions for variable density asymmetric incompressible fluids, J. Math. Anal. Appl., 387 (2012), 953–969. https://doi.org/10.1016/j.jmaa.2011.10.015 doi: 10.1016/j.jmaa.2011.10.015
    [4] P. Braz e Silva, F. W. Cruz, M. Loayza, M. A. Rojas-Medar, Global unique solvability of nonhomogeneous asymmetric fluids: a Lagrangian approach, J. Differ. Equations, 269 (2020), 1319–1348. https://doi.org/10.1016/j.jde.2020.01.001 doi: 10.1016/j.jde.2020.01.001
    [5] B. Desjardins, Regularity results for two-dimensional flows of multiphase viscous fluids, Arch. Rational Mech. Anal., 137 (1997), 135–158. https://doi.org/10.1007/s002050050025 doi: 10.1007/s002050050025
    [6] R. Danchin, Local and global well-posedness results for flows of inhomogeneous viscous fluids, Adv. Differential Equations, 9 (2004), 353–386. https://doi.org/10.57262/ade/1355867948 doi: 10.57262/ade/1355867948
    [7] R. Danchin, P. B. Mucha, A Lagrangian approach for the incompressible Navier-Stokes equations with variable density, Commun. Pure Appl. Math., 65 (2012), 1458–1480. https://doi.org/10.1002/cpa.21409 doi: 10.1002/cpa.21409
    [8] C. He, J. Li, B. Q. Lü, Global well-posedness and exponential stability of 3D Navier-Stokes equations with density-dependent viscosity and vacuum in unbounded domains, Arch. Rational Mech. Anal., 239 (2021), 1809–1835. https://doi.org/10.1007/s00205-020-01604-5 doi: 10.1007/s00205-020-01604-5
    [9] J. C. Huang, M. Paicu, P. Zhang, Global solutions to 2D inhomogeneous Navier-Stokes system with general velocity, J. Math. Pures Appl., 100 (2013), 806–831. https://doi.org/10.1016/j.matpur.2013.03.003 doi: 10.1016/j.matpur.2013.03.003
    [10] X. D. Huang, Y. Wang, Global strong solution of 3D inhomogeneous Navier-Stokes equations with density-dependent viscosity, J. Differ. Equations, 259 (2015), 1606–1627. https://doi.org/10.1016/j.jde.2015.03.008 doi: 10.1016/j.jde.2015.03.008
    [11] O. A. Ladyzhenskaya, V. A. Solonnikov, N. N. Ural'ceva, Linear and quasilinear equations of parabolic type, Providence: American Mathematical Society, 1968. https://doi.org/10.1090/mmono/023
    [12] L. X. Liu, X. Zhong, Global existence and expinential decay of strong solutions for 2D nonhomogeneous micropolar fluids with density-dependent viscosity, J. Math. Phys., 62 (2021), 061508. https://doi.org/10.1063/5.0055689 doi: 10.1063/5.0055689
    [13] B. Q. Lü, S. S. Song, On local strong solutions to the three-dimensional nonhomogeneous incompressible Navier-Stokes equations with density-dependent viscosity and vacuum, Nonlinear Anal. Real, 46 (2019), 58–81. https://doi.org/10.1016/j.nonrwa.2018.09.001 doi: 10.1016/j.nonrwa.2018.09.001
    [14] C. Y. Qian, Y. Qu, Global well-posedness for 3D incompressible inhomogeneous asymmetric fluids with density-dependent viscosity, J. Differ. Equations, 306 (2022), 333–402. https://doi.org/10.1016/j.jde.2021.10.045 doi: 10.1016/j.jde.2021.10.045
    [15] C. Y. Qian, H. Chen, T. Zhang, Global existence of weak solutions for 3D incompressible inhomogeneous asymmetric fluids, Math. Ann., 386 (2023), 1555–1593. https://doi.org/10.1007/s00208-022-02427-3 doi: 10.1007/s00208-022-02427-3
    [16] C. Y. Qian, B. B. He, T. Zhang, Global well-posedness for 2D inhomogeneous asymmetric fluids with large initial data, Sci. China Math., 67 (2024), 527–556. https://doi.org/10.1007/s11425-022-2099-1 doi: 10.1007/s11425-022-2099-1
    [17] J. W. Zhang, Global well-posedness for the incompressible Navier-Stokes equations with density-dependent viscosity coefficient, J. Differ. Equations, 259 (2015), 1722–1742. https://doi.org/10.1016/j.jde.2015.03.011 doi: 10.1016/j.jde.2015.03.011
    [18] P. X. Zhang, M. X. Zhu, Global regularity of 3D nonhomogeneous incompressible micropolar fluids, Acta Appl. Math., 161 (2019), 13–34. https://doi.org/10.1007/s10440-018-0202-1 doi: 10.1007/s10440-018-0202-1
    [19] L. Zhou, C. L. Tang, Global well-posedness to the 3D Cauchy problem of nonhomogeneous Navier-Stokes equations with density-dependent viscosity and large initial velocity, J. Math. Phys., 64 (2023), 111509. https://doi.org/10.1063/5.0144133 doi: 10.1063/5.0144133
    [20] X. Zhong, Global well-posedness for 3D nonhomogeneous micropolar fluids with density-dependent viscosity, Bull. Malays. Math. Sci. Soc., 46 (2023), 23. https://doi.org/10.1007/s40840-022-01399-6 doi: 10.1007/s40840-022-01399-6
  • Reader Comments
  • © 2024 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(951) PDF downloads(71) Cited by(0)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog