Loading [MathJax]/jax/element/mml/optable/MathOperators.js
Research article Special Issues

Uniform bound of the highest-order energy for three dimensional inhomogeneous incompressible elastodynamics

  • We are concerned with the time growth of the highest-order energy of three-dimensional inhomogeneous incompressible isotropic elastodynamics. Utilizing Klainerman's generalized energy method, refined weighted estimates, and the Keel-Smith-Sogge estimate [J. Anal. Math., 87: 265-279, 2002], it is justified that the highest-order generalized energy is uniformly bounded for all time.

    Citation: Xiufang Cui, Xianpeng Hu. Uniform bound of the highest-order energy for three dimensional inhomogeneous incompressible elastodynamics[J]. Communications in Analysis and Mechanics, 2025, 17(2): 429-461. doi: 10.3934/cam.2025018

    Related Papers:

    [1] Senhao Duan, Yue MA, Weidong Zhang . Conformal-type energy estimates on hyperboloids and the wave-Klein-Gordon model of self-gravitating massive fields. Communications in Analysis and Mechanics, 2023, 15(2): 111-131. doi: 10.3934/cam.2023007
    [2] Yuxuan Chen . Global dynamical behavior of solutions for finite degenerate fourth-order parabolic equations with mean curvature nonlinearity. Communications in Analysis and Mechanics, 2023, 15(4): 658-694. doi: 10.3934/cam.2023033
    [3] Shuai Xi . A Kato-type criterion for the inviscid limit of the nonhomogeneous NS equations with no-slip boundary condition. Communications in Analysis and Mechanics, 2024, 16(4): 896-909. doi: 10.3934/cam.2024039
    [4] Shuyue Ma, Jiawei Sun, Huimin Yu . Global existence and stability of temporal periodic solution to non-isentropic compressible Euler equations with a source term. Communications in Analysis and Mechanics, 2023, 15(2): 245-266. doi: 10.3934/cam.2023013
    [5] Farrukh Dekhkonov, Wenke Li, Weipeng Wu . On the time-optimal control problem for a fourth order parabolic equation in the three-dimensional space. Communications in Analysis and Mechanics, 2025, 17(2): 413-428. doi: 10.3934/cam.2025017
    [6] Jiangbo Han, Runzhang Xu, Chao Yang . Continuous dependence on initial data and high energy blowup time estimate for porous elastic system. Communications in Analysis and Mechanics, 2023, 15(2): 214-244. doi: 10.3934/cam.2023012
    [7] Yang Liu . Global attractors for a nonlinear plate equation modeling the oscillations of suspension bridges. Communications in Analysis and Mechanics, 2023, 15(3): 436-456. doi: 10.3934/cam.2023021
    [8] Chunyou Sun, Junyan Tan . Attractors for a Navier–Stokes–Allen–Cahn system with unmatched densities. Communications in Analysis and Mechanics, 2025, 17(1): 237-262. doi: 10.3934/cam.2025010
    [9] Henryk Żołądek . Normal forms, invariant manifolds and Lyapunov theorems. Communications in Analysis and Mechanics, 2023, 15(2): 300-341. doi: 10.3934/cam.2023016
    [10] Haifa El Jarroudi, Mustapha El Jarroudi . Asymptotic behavior of a viscous incompressible fluid flow in a fractal network of branching tubes. Communications in Analysis and Mechanics, 2024, 16(3): 655-699. doi: 10.3934/cam.2024030
  • We are concerned with the time growth of the highest-order energy of three-dimensional inhomogeneous incompressible isotropic elastodynamics. Utilizing Klainerman's generalized energy method, refined weighted estimates, and the Keel-Smith-Sogge estimate [J. Anal. Math., 87: 265-279, 2002], it is justified that the highest-order generalized energy is uniformly bounded for all time.



    The motion dynamics of incompressible isotropic elastodynamics is characterized as a wave system in Lagrangian coordinates, which inherently satisfies the null condition. Based on this structure, a series of studies have been conducted to establish the global well-posedness of classical solutions to this system; see[1,2]. However, these studies reveal a certain time growth for the highest-order generalized energy. In this paper, we investigate the time growth of the Sobolev norm for classical solutions to three-dimensional inhomogeneous incompressible isotropic elastodynamics with small initial perturbation and establish the uniform bound for the highest-order energy.

    Before presenting the main result of this paper, we briefly review related known results. For three-dimensional elastic waves, John [3] proved that the genuine nonlinearity condition leads to singularity formation even for arbitrarily small spherically symmetric displacement. We also refer readers to [4] regarding large displacement singularity. The existence of almost global solutions was established in [5,6] for three-dimensional quasilinear wave equations with sufficiently small initial data. Significant contributions toward global existence were independently made by Sideris [7,8] and Agemi [9] under the assumption that nonlinearity satisfies the null condition in three dimensions. In terms of three-dimensional incompressible elastodynamics, the only waves presented in the isotropic systems are shear waves, which are linearly degenerate. The global existence of a solution was demonstrated by Sideris and Thomases in [1,2] through two different methods. It is more challenging to obtain the global existence for the two-dimensional incompressible elastodynamics due to the weaker dispersive decay. In [10], the authors proved almost global existence for a two-dimensional incompressible system in Eulerian coordinates. By introducing the concept of strong null condition and observing that the incompressible elastodynamics automatically satisfies such strong null structure in Lagrangian coordinates, Lei [11] successfully proved the global well-posedness for two-dimensional incompressible elastodynamics by the method of Klainerman and Alinhac's ghost weight method [12]. We also see [13] for a different approach using the spacetime resonance method. All the aforementioned works considered the homogeneous fluids. In [14], the authors established the global well-posedness for the three-dimensional inhomogeneous incompressible elastodynamics in Lagrangian coordinates. It is noteworthy that the upper bound of the highest-order generalized energy in those studies depends on time. Utilizing the Klainerman's generalized energy method, an analysis of the inherent structure of the system and the ghost weight method, [15,16] established the uniform bound for the highest-order generalized energy estimates for two-dimensional and three-dimensional incompressible elastodynamics, respectively. Based on the above foundational works, it is natural to verify the uniform bound for the highest-order generalized energy for three-dimensional inhomogeneous incompressible isotropic elastodynamics. To establish the time growth of the Sobolev norm of classical solution, two novel methods are presented in this paper. First, based on the Sobolev embedding inequality and the structure of the system, the refined decay rates were derived for the solution in the domain away from the light cone. Second, we apply the KSS-type estimate to overcome the difficulties posed by insufficient time decay resulting from density perturbation.

    This paper is organized as follows. In Section 2, we introduce the system of three-dimensional inhomogeneous isotropic elastodynamics and define the notations utilized throughout this paper. Besides, the main result along with several useful lemmas are presented in this section. The energy estimates are discussed in Section 3.

    We first formulate the inhomogeneous isotropic elastodynamics and denote some notations that are used frequently in this paper.

    For any given smooth flow map X(t,x), we call it incompressible if

    Ωdx=ΩtdX,Ωt={X(t,x)|xΩ}

    for any smooth bounded connected domain Ω, which yields that

    det(X)=1.

    Denote

    X(t,x)=x+v(t,x).

    Simple calculation shows that the incompressible condition is equivalent to

    v+12[(trv)2tr(v)2]+det(v)=0. (2.1)

    Without loss of generality, we assume that the density of fluid is a small perturbation around the constant state 1, that is, ρ(x)=1+η(x). For the inhomogeneous isotropic material, the motion of the elastic fluid in the Lagrangian coordinate is determined by

    L(X;T,Ω)=T0Ω(12ρ(x)|tX|2W(X)+p(t,x)[det(X)1])dxdt. (2.2)

    Here W(X)C is the strain energy function. p(t,x) is a Lagrangian multiplier that is used to force the flow maps to be incompressible. To simplify the presentation, we only study the typical Hookean case for which the strain energy functional is given by

    W(X)=12|X|2.

    By calculating the variation of (2.2), we obtain the equation

    ρ2tvΔv=(X)Tp. (2.3)

    Now, we introduce the following derivative vector fields

    t=0,=(1,2,3)and=(0,1,2,3).

    The scaling operator is denoted by

    S=tt+rr.

    Here, the radial derivative is defined by r=xr, r=|x|. The angular momentum operators are denoted by

    Ω=x.

    In the application, we usually use the modified rotational operators and scaling operator; that is, for any vector v and scalar p and ρ, we set

    ˜Sp=Sp,˜Sρ=Sρ,˜Sv=(S1)v,˜Ωip=Ωip,˜Ωiρ=Ωiρ,˜Ωiv=Ωiv+Uiv,i=1,2,3,

    where

    U1=e2e3e3e2,U2=e3e1e1e3,U3=e1e2e2e1.

    Let

    Γ=(Γ1,,Γ8)=(,˜Ω,˜S)

    and for any multi-index α=(α1,α2,,α8)N8, we denote

    Γα=Γα11Γα88.

    We apply the derivatives Γα to the equations (2.1) and (2.3), and then the three-dimensional inhomogeneous isotropic elastodynamics can be written as

    ρ2tΓαvΔΓαv=(X)TΓαpβ+γ=α|β||α|CβαΓγη2tΓβvβ+γ=α|β||α|Cβα(X)T(Γγv)T(2tΓβvΔΓβv+ι1+ι2=βCι1βΓι2η2tΓι1v)=:Nα (2.4)

    with the incompressible condition

    Γαv+β+γ=αi,j=1,2,3,i<jCβα(iΓβvijΓγvjiΓβvjjΓγvi)+β+γ+ι=αCβαCγαβ|1Γβv12Γβv13Γβv11Γβv22Γβv23Γβv21Γβv32Γβv33Γβv3|=0. (2.5)

    Here the binomial coefficient Cβα is given by

    Cβα=α!β!(αβ)!.

    We denote the Klainerman's generalized energy by

    Eκ(t)=|α|κ1R3(|tΓαv|2+|Γαv|2)dx.

    We also define the weighted energy norm

    Xκ(t)=|α|κ2R3tr2|2Γαv|2dx.

    To describe the space of initial data, we introduce the time-independent analogue of Γ as

    Λ=(,rr1,˜Ω)

    and the space of initial data is defined by

    HκΛ={(f,g):|α|κ1(ΛαfL2+ΛαfL2+ΛαgL2)<}.

    As the first step to investigate our problem, we introduce the following lemma, which helps us to solve the additional terms resulting from density perturbation. Let

    LEκ(T)=|α|κ1T0Rnr1+2μr2μ(|Γαv|+|Γαv|r)2dxdt

    with μ(0,12) and μ>μ. Without loss of generality, we choose μ=14 and μ=12 in this paper. In the case of μ=μ, we see LEκ(T) is the KSS norm, and we denote it by KSSκ(T).

    Lemma 1. Let f0=[r/(1+r)]2μ, fk=r/(r+2k) with k1, μ(0,1/2), and v be the solution to the equation 2tvc2Δv+hababv=N in [0,T]×Rn with hab=hba, 0a,bn|hab|min(1,c2)/2 for any integer n3. Then there exists a positive constant C0 that depends only on the dimension n such that

    sup0tTRn|v|2(t)dx+LE1(T)+(ln(2+T))1KSS1(T)C0Rn|v|2(0)dx+C0T0Rn[(|h|+|h|r12μr2μ)|v|(|v|+|v|r)]dxdt+C0|T0RntvNdxdt|+C0supk0|T0Rnfk(rv+n12rv)Ndxdt|,

    where |h|=na,b=0|hab| and |h|=na,b,c=0|chab|.

    This lemma can be found in [17]. See also [18,19] and references therein.

    Based on the previous statement, we are ready to show the main result of this paper.

    Theorem 1. Let W(X)=12|X|2 be an isotropic Hookean strain energy function and (v0,tv0)HκΛ with κ12. Let C0>0 be given constant in Lemma 1. Suppose v0 satisfies the structural constraint condition (2.1) and

    Eκ(0)=|α|κ1(tΛαv02L2+Λαv02L2)ε.

    If

    rΛαηL2δfor|α|κ,

    then there exist two sufficiently small constants, ε0,δ0 and constant C1, that depend only on κ and C0 such that if εε0 and δδ0, the system (2.3) has a unique global classical solution that satisfies

    Eκ(t)C1ε

    uniformly for all t[0,+).

    In this part, we establish several lemmas that are crucial for the energy estimates. Throughout this paper, we denote =(1+||2)12. The notation fg stands for fCg for some generic constant C>0, which may vary from line to line. In the process of deriving the energy estimates, we usually separate the whole integration domain R3 into two parts:

    R={xR3:rt/8},Rc={xR3:r>t/8}.

    We first recall the Sobolev-type inequalities, which were justified in [8].

    Lemma 2. For any vC0(R3)3, r=|x| and ˜r=|y|, we have

    r12v(x)L|α|1˜ΩαvL2,rv(x)L|α|1r˜Ωαv12L2(|y|r)˜Ωαv12L2(|y|r),rtrv(x)L|α|1t˜r˜r˜ΩαvL2(|y|r)+|α|2t˜r˜ΩαvL2(|y|r).

    The following lemma concerns the dispersive decay of solutions in the domain away from the light cone.

    Lemma 3. For any vH2(R3), there holds

    tvL(rt/8)vL2(rt/4)+trvL2(rt/4)+tr2vL2(rt/4).

    Proof. The proof can be found in Lemma 4.3 in [16].

    Let s=8rt. We introduce a radial cutoff function ξ(s)C0 that satisfies

    ξ(s)={1,s1,0,s2.

    It is easy to observe from Lemmas 2 and 3 that

    tvLtvL(R)+r(1ξ(s))vLvL2(rt/4)+trvL2(rt/4)+tr2vL2(rt/4)+ι1,|ι2|1ι1r˜Ωι2[(1ξ(s))v]L2E123(t)+X123(t). (2.6)

    To control the weighted energy norm by the generalized energy norm, the pressure must be estimated via the system (2.4)-(2.5). A similar proof can be found in [14]. For a self-contained presentation, we include its proof below.

    Lemma 4. For any integer κ6 and multi-index α satisfying |α|κ2, if E[κ/2]+3(t) and rΛαηL2 are small. Then there holds

    tΓαpL2+tρ2tΓαvΔΓαvL2δX12κ(t)+E12κ(t)X12[κ/2]+3(t)+E12[κ/2]+3(t)X12κ(t).

    Proof. It observes from (2.4) that

    Γαp=(X)T(ρ2tΓαvΔΓαv)β+γ=α|β||α|Cβα(X)TΓγη2tΓβvβ+γ=α|β||α|Cβα(Γγv)T(2tΓβvΔΓβv+ι1+ι2=βCι1βΓι2η2tΓι1v). (2.7)

    We apply Δ1 to the above equality and take the L2 norm. By the L2 boundness of the Riesz operator, one has

    ΓαpL2Δ1(ρ2tΓαvΔΓαv)L2+(v)T(ρ2tΓαvΔΓαv)L2+β+γ=α|β||α|(X)TΓγη2tΓβvL2+β+γ=α|β||α|(Γγv)T(2tΓβvΔΓβv+ι1+ι2=βCι1βΓι2η2tΓι1v)L2. (2.8)

    Special attention is paid to the first term on the right-hand side of (2.8). We apply the derivative operator (2tΔ) to the incompressible condition (2.5). Without loss of generality, we assume at least one order derivative operator of (2tΔ) works on the first component of the velocity field in the second line of (2.5). Then we obtain

    (2tΔ)(Γαv)=β+γ=αi,j=1,2,3,i<jCβα(2tΔ)(iΓβvijΓγvjiΓβvjjΓγvi)β+γ+ι=αCβαCγαβ(2tΔ)|1Γβv12Γβv13Γβv11Γβv22Γβv23Γβv21Γβv32Γβv33Γβv3|=β+γ=αm+n=1,i=0,1,2,3Cβα(m+1iΓβv12niΓγv2+m+1iΓβv13niΓγv32miΓβv1n+1iΓγv23miΓβv1n+1iΓγv3m+1iΓβv11niΓγv2+m+1iΓβv23niΓγv3+1miΓβv1n+1iΓγv23miΓβv2n+1iΓγv3m+1iΓβv22niΓγv3m+1iΓβv11niΓγv3+2miΓβv2n+1iΓγv3+1miΓβv1n+1iΓγv3)+β+γ+ι=αm2+m2+m3=1,i=0,1,2,3CβαCγαβ[m1+1iΓβv1(2m2iΓγv23m3iΓιv33m2iΓγv22m3iΓιv33m2iΓγv21m3iΓιv31m2iΓγv23m3iΓιv31m2iΓγv22m3iΓιv32m2iΓγv21m3iΓιv3)].

    Based on the above equality, we handle the first term on the right-hand side of (2.8) as follows

    Δ1(ρ2tΓαvΔΓαv)L2Δ1[(ρ2tΔ)(Γαv)]L2+i=1,2,3Δ1(η2tΓαvi)L2(1+vL)vLρ2tΓαvΔΓαvL2+β+γ=α,|β||α|i,j=1,2,3,ij2ΓβviΓγvjL2+β+γ+ι=α,|β||α|i,j,k=1,2,3,ijk2ΓβviΓγvjΓιvkL2+β+γ=α,|β||α|i,j=1,2,3,ijΔ1(η2ΓβviΓγvj)L2+β+γ+ι=α,|β||α|i,j,k=1,2,3,ijkΔ1(η2ΓβviΓγvjΓιvk)L2+i=1,2,3Δ1(η2tΓαvi)L2+i,j=1,2,3ijΔ1(η2tΓαvivj)L2+i,j,k=1,2,3ijkΔ1(η2tΓαvivjvk)L2. (2.9)

    By (2.8), (2.9), the Sobolev embedding inequality, and the smallness of ηL3, we have

    ΓαpL2(1+vL)vLρ2tΓαvΔΓαvL2+β+γ=α,|β||α|i,j=1,2,3,ij2ΓβviΓγvjL2+β+γ+ι=α,|β||α|1,j,k=1,2,3,ijk2ΓβviΓγvjΓιvkL2+β+γ=α|β||α|Γγη2tΓβvL2+β+γ=α|β||α|(Γγv)T(2tΓβvΔΓβv+ι1+ι2=βCι1βΓι2η2tΓι1v)L2+i=1,2,3η2tΓαviL65+i,j=1,2,3ijη2tΓαvivjL65+i,j,k=1,2,3ijη2tΓαvivjvkL65.

    By (2.4), the above inequality and the smallness of E[κ/2]+3(t), one obtains

    ΓαpL2+ρ2tΓαvΔΓαvL2β+γ=α,|β||α|i,j=1,2,3,ij2ΓβviΓγvjL2+β+γ+ι=α,|β||α|i,j,k=1,2,3,ijk2ΓβviΓγvjΓιvkL2+β+γ=α|β||α|Γγη2tΓβvL2+β+γ=α|β||α|(Γγv)T(2tΓβvΔΓβv+ι1+ι2=βCι1βΓι2η2tΓι1v)L2+i=1,2,3η2tΓαviL65+i,j=1,2,3ijη2tΓαvivjL65+i,j,k=1,2,3ijkη2tΓαvivjvkL65.

    We deduce from Lemma 2 that

    tΓαpL2+tρ2tΓαvΔΓαvL2β+γ=α|γ|<|β|<|α|tr2ΓβvL2rΓγvL+β+γ=α|β||γ|rtr2ΓβvLΓγvL2+β+γ+ι=α|γ|,|ι|<|β|<|α|tr2ΓβvL2rΓγvLΓιvL+β+γ+ι=α|β|,|ι||γ|rtr2ΓβvLΓγvL2ΓιvL+β+γ+ι=α|β|,|γ||ι|rtr2ΓβvLΓγvLΓιvL2+β+γ=α|γ|<|β|<|α|rΓγηLtr2tΓβvL2+β+γ=α|β||γ|ΓγηL2rtr2tΓβvL+rηL3(1+vL)2tr2tΓαvL2+β+γ=α|γ|<|β|<|α|rΓγvLtr(2tΓβvΔΓβv+ι1+ι2=βCι1βΓι2η2tΓι1v)L2+β+γ=α|β||γ|ΓγvL2rtr(2tΓβvΔΓβv+ι1+ι2=βCι1βΓι2η2tΓι1v)LδX12κ(t)+E12κ(t)X12[κ/2]+3(t)+E12[κ/2]+3(t)X12κ(t).

    It completes the proof.

    As an application of the above result, we establish the estimate of weighted energy.

    Lemma 5. Let vHκΓ(R3) be the solution to the system (2.3) with the constraint condition (2.1). For any integer κ6, if Eκ(t) and rΛαηL2 are small, then we have

    Xκ(t)Eκ(t).

    Proof. For any multi-index α satisfying |α|κ2, we apply Γα to Lemma 3.3 in [11] to get

    Xκ(t)Eκ(t)+t2tΓαvΔΓαv2L2.

    It follows from the above inequality and Lemma 4 that

    Xκ(t)Eκ(t)+δ2Xκ(t)+E[κ/2]+3(t)Xκ(t)+Eκ(t)X[κ/2]+3(t).

    By the smallness of δ and Eκ(t), we arrive at the lemma.

    In preparation for the energy estimates, more detailed analysis of pressure is needed. In what follows, we always assume that Eκ(t) and rΛαηL2 are small.

    Lemma 6. For any integer κ8 and multi-index α satisfying |α|κ1, we have

    ΓαpL2E12κ(t) (2.10)

    and

    tΓαpL2(Rc)+tρ2tΓαvΔΓαvL2(Rc)E12κ(t). (2.11)

    Proof. Following the calculations in Lemma 4, we arrive at (2.9). Special attention is paid to the last three terms on the right-hand side of (2.9). From (2.4), one has

    2tΓαv=ρ1ΔΓαvρ1(X)TΓαpβ+γ=α|β||α|Cβαρ1Γγη2tΓβvβ+γ=α|β||α|Cβαρ1(X)T(Γγv)T(2tΓβvΔΓβv+ι1+ι2=βCι1βΓι2η2tΓι1v).

    We use the above equality and the Sobolev embedding inequality to solve the last three terms on the right-hand side of (2.9) by

    Δ1[ρ1η(1+v)2ΔΓαv]L2+Δ1[ρ1η(1+v)2(X)TΓαp]L2+β+γ=α|β||α|Δ1[ρ1η(1+v)2(X)T(Γγv)T(2tΓβvΔΓβv+ι1+ι2=βCι1βΓι22tΓι1v)]L2+β+γ=α|β||α|Δ1[ρ1η(1+v)2Γγη2tΓβv]L2η(1+v)2ΓαvL2+i+j+k=1i(ρ1)jηk(1+v)2ΓαvL65+η(1+v)2(X)TΓαpL65+β+γ=α|β||α|η(1+v)2(X)T(Γγv)T(2tΓβvΔΓβv+ι1+ι2=βCι1βΓι2η2tΓι1v)L65+β+γ=α|β||α|η(1+v)2Γγη2tΓβvL65.

    Utilizing (2.4), (2.8), the above estimate, and the Sobolev embedding inequality, we arrive at

    ΓαpL2+ρ2tΓαvΔΓαvL2β+γ=α|γ|<|β|<|α|2ΓβvL2ΓγvL+β+γ=α|β||γ|2ΓβvLΓγvL2+β+γ+ι=α|γ|,|ι|<|β|<|α|2ΓβvL2ΓγvLΓιvL+β+γ+ι=α|β|,|ι||γ|2ΓβvLΓγvL2ΓιvL+β+γ+ι=α|β|,|γ||ι|2ΓβvLΓγvLΓιvL2+β+γ=α|γ|<|β|<|α|ΓγηL2tΓβvL2+β+γ=α|β||γ|ΓγηL22tΓβvL+1|ι|3ιηL2ΓαvL2+β+γ=α|γ|<|β|<|α|ΓγvL2tΓβvΔΓβv+ι1+ι2=βCι1βΓι2η2tΓι1vL2+β+γ=α|β||γ|ΓγvL22tΓβvΔΓβv+ι1+ι2=βCι1βΓι2η2tΓι1vLδE12κ(t)+E12κ(t)E12[κ/2]+4(t). (2.12)

    The smallness of δ and Eκ(t) leads to (2.10).

    To verify (2.11), we use (2.12), Lemma 2, and the smallness of Eκ(t) to get

    tΓαpL2(Rc)+tρ2tΓαvΔΓαvL2(Rc)β+γ=α|γ|<|β|<|α|2ΓβvL2rΓγvL+β+γ=α|β||γ|r2ΓβvLΓγvL2+β+γ+ι=α|β|,|ι||γ|2ΓβvLΓγvL2rΓιvL+β+γ+ι=α|β|,|γ||ι|2ΓβvLrΓγvLΓιvL2+β+γ+ι=α|γ|,|ι|<|β|<|α|2ΓβvL2rΓγvLΓιvL+β+γ=α|γ|<|β|<|α|rΓγηL2tΓβvL2+β+γ=α|β||γ|ΓγηL2r2tΓβvL+1|ι|3rιηL2ΓαvL2+β+γ=α|γ|<|β|<|α|rΓγvL2tΓβvΔΓβv+ι1+ι2=βCι1βΓι2η2tΓι1vL2+β+γ=α|β||γ|ΓγvL2r(2tΓβvΔΓβv+ι1+ι2=βCι1βΓι2η2tΓι1v)LE12κ(t)E12[κ/2]+4(t)+δE12κ(t).

    The smallness of δ and Eκ(t) implies (2.11).

    In the subsequent part, we present the improved decay properties for the third-order spatial derivatives of unknown variables in the domain away from the light cone.

    Lemma 7. For any integer κ10 and multi-index α satisfying |α|[κ/2], it holds that

    t23ΓαvL2(R)E12[κ/2]+5(t).

    Proof. We apply the derivative operator to the equation (2.4) to get

    ΔΓαv=η2tΓαv+ρ2tΓαv+(X)TΓαp+(X)TΓαp+β+γ=α|β||α|CβαΓγη2tΓβv+β+γ=α|β||α|CβαΓγη2tΓβv+β+γ=α|β||α|Cβα[(X)T(Γγv)T(2tΓβvΔΓβv+ι1+ι2=βCι1βΓι2η2tΓι1v)].

    By multiplying the above equality by t2ξ(s) and taking the L2 inner product, one has

    R3|t2ξ(s)ΔΓαv|2dx7R3|t2ξ(s)η2tΓαv|2dx+7R3|t2ξ(s)ρ2tΓαv|2dx+7R3|t2ξ(s)(X)TΓαp|2dx+7R3|t2ξ(s)(X)T2Γαp|2dx+β+γ=α|β||α|7R3|Cβαt2ξ(s)Γγη2tΓβv|2dx+β+γ=α|β||α|7R3|Cβαt2ξ(s)Γγη2tΓβv|2dx+β+γ=α|β||α|7R3|Cβαt2ξ(s)[(X)T(Γγv)T(2tΓβvΔΓβv+ι1+ι2=βCι1βΓι2η2tΓι1v)]|2dx. (2.13)

    Since for any fH2(R3), the following Sobolev embedding inequality holds true

    fLf12L22f12L2. (2.14)

    By Lemma 2, the first term on the right-hand side of (2.13) is estimated by

    7R3|t2ξ(s)η2tΓαv|2dxR3|trξ(s)ηrtΓαv|2dx+R3|tξ(s)ηt˜SΓαv|2dx(ξ(s)rtrrtΓαv2L+tξ(s)t˜SΓαv2L)η2L2(rtrrtΓαv2L+(tξ(s)t˜SΓαv)L22(tξ(s)t˜SΓαv)L2)η2L2[rtrrtΓαv2L+(ξ(s)t˜SΓαvL2+ξ(s)trt˜SΓαvL2)(t1ξ(s)t˜SΓαvL2+ξ(s)t˜SΓαvL2+ξ(s)trt2˜SΓαvL2)]η2L2δ2E[κ/2]+4(t). (2.15)

    By a similar argument, we handle the second term on the right-hand side of (2.13) as follows

    7R3|t2ξ(s)ρ2tΓαv|2dx7R3|trξ(s)ρrtΓαv|2dx+CR3(|tξ(s)ρt˜SΓαv|2+|tξ(s)ρtΓαv|2)dx7R3|r2ξ(s)ρ2rΓαv|2dx+CR3(|tξ(s)ρt˜SΓαv|2+|tξ(s)ρtΓαv|2+|ξ(s)ρrr˜SΓαv|2+|ξ(s)ρrrΓαv|2)dx7R3|r2ξ(s)ρ3Γαv|2dx+C(trΓαv2L2(R)+tr˜SΓαv2L2(R))716R3|t2ξ(s)3Γαv|2dx+CE[κ/2]+3(t). (2.16)

    In view of (2.7), the smallness of and , one has

    We substitute the above inequality into the third term on the right-hand side of (2.13) to get

    We come back to the fourth term on the right-hand side of (2.13). We apply the divergence operator to the equality (2.7) and multiply on both sides of the resulting equality. By taking the inner product, one has

    (2.17)

    For the first two terms on the right-hand side of (2.17), we deduce from Lemma 3 that

    For the third and fourth terms on the right-hand side of (2.17), we have

    (2.18)

    By the definition of and the Sobolev embedding inequality, we solve the first term on the right-hand side of (2.18) by

    In terms of (2.5) and the definition of , the last two terms on the right-hand side of (2.18) are estimated by

    Along the same line, the last term on the right-hand side of (2.17) is dealt with by

    Combining all the estimates, Lemma 4, and the fact

    the fourth term on the right-hand side of (2.13) is estimated by

    We employ the similar method as (2.15) and (2.16) to estimate the fifth and sixth terms on the right-hand side of (2.13) by

    For the last term on the right-hand side of (2.13), we have

    Collecting the above estimates and the fact

    we complete the proof.

    By Lemma 7, we obtain the following estimate.

    Lemma 8. For any integer and multi-index satisfying , there holds

    (2.19)
    (2.20)

    Proof. By the definition of and Lemma 7, we have

    (2.21)

    The inequalities (2.14), (2.21), and Lemma 7 yield that

    To consider (2.20), the definition of , combined with inequalities (2.14) and (2.21) and Lemmas 3 and 7, implies that

    which implies the desired.

    Before concluding this section, we formulate the following two lemmas, which are utilized in the process of deriving energy estimates.

    Lemma 9. For any integer and multi-index satisfying , we have

    Proof. We separate two cases to consider this lemma. For the case , we use the Sobolev embedding inequality (2.14) to get

    For the integer satisfying , one has

    (2.22)

    We estimate the first term on the right-hand side of (2.22) by

    By (2.4), we solve the second term on the right-hand side of (2.22) as follows

    (2.23)

    For the first term on the right-hand side of (2.23), by Lemmas 3 and 4, we have

    Adopting the same method as was used in (2.17), we estimate the second term on the right-hand side of (2.23) by

    For the third term on the right-hand side of (2.23), by the definition of and (2.6), we obtain

    The last term on the right-hand side of (2.23) is estimated by

    We verify the case . By Lemmas 2 and 6, we have

    Collecting all of the estimates together, we justify the lemma.

    Lemma 10. For any integer and multi-index satisfying , one has

    Proof. We deduce from the definition of that

    By Lemma 2 and (2.6), one has

    This completes the proof.

    This section is devoted to the energy estimates. For any integer and multi-index satisfying , we apply Lemma 1 to the system (2.4)-(2.5) to get

    (3.1)

    For the second term on the right-hand side of (3.1), by Lemma 2, one has

    By utilizing (2.4), we formulate the third term on the right-hand side of (3.1) as follows

    (3.2)

    We handle term by term on the right-hand side of (3.2). For the first term on the right-hand side of (3.2), since is composed of elements of the form , , and the constant , where , it follows that

    (3.3)

    For the first term on the right-hand side of (3.3), by (2.6) and Lemma 6, we have

    To handle the second and third terms on the right-hand side of (3.3), we apply (2.11) and (2.19) to show

    Along the same line, the last term on the right-hand side of (3.3) is estimated by

    For the second term on the right-hand side of (3.2), by the definition of and (2.6), one has

    It is left to estimate the last term on the right-hand side of (3.2). Two cases are considered. By Lemmas 2, 4, and (2.6), we solve the case by

    Utilizing Lemmas 9 and 10, we solve the case by

    By summing up the above estimates, we deduce that

    We continue to handle the last term on the right-hand side of (3.1). The identity (2.4) yields that

    (3.4)

    We use Lemmas 9 and 10 to handle the first two terms on the right-hand side of (3.4) by

    We rewrite the last term on the right-hand side of (3.4) as follows

    (3.5)

    For the first two terms on the right-hand side of (3.5), the Lemma 6 and (2.19) imply that

    We use (2.7) to formulate the last term on the right-hand side of (3.5) by

    (3.6)

    By (2.19) and Lemma 6, we estimate the first term on the right-hand side of (3.6) by

    (3.7)

    The second and third terms on the right-hand side of (3.6) can be solved using the same method employed by the first two terms of (3.4).

    For the last term on the right-hand side of (3.6), we observe that

    (3.8)

    In view of (2.5), we write the first term on the right-hand side of (3.8) by

    (3.9)

    The first two terms on the right-hand side of (3.9) are dealt with by the same method as (3.7). By (2.6), the third and fourth terms on the right-hand side of (3.9) are estimated by

    To consider the fifth term on the right-hand side of (3.9), we separate two cases to consider it. For the case , by (2.20), we have

    To consider the case , we use the integration by parts to get

    (3.10)

    For the first term on the right-hand side of (3.10), we have

    We use (2.5) to estimate the second term on the right-hand side of (3.10) by

    For the last term on the right-hand side of (3.9), it follows from the Sobolev embedding inequality, (2.6), and (2.19), that

    To solve the second term on the right-hand side of (3.8), we rewrite it as follows

    (3.11)

    In view of (2.5), we formulate the first term on the right-hand side of (3.11) by

    (3.12)

    Since the second term has analogous estimates to the first term of (3.12), it suffices to concentrate on the first term. We observe that

    (3.13)

    Here, we restrict our analysis to the first term on the right-hand side of (3.13), as the remaining two terms have similar estimates.

    (3.14)

    For the first term on the right-hand side of (3.14), we deduce from (2.4) that

    (3.15)

    We formulate the first term on the right-hand side of (3.15) as follows

    (3.16)

    The first three terms on the right-hand side of (3.16) can be estimated by

    We utilize (2.19) to solve the last two terms on the right-hand side of (3.16) by

    Along the same line, the second term on the right-hand side of (3.15) can be handled by

    The same estimates hold for the last two terms on the right-hand side of (3.14).

    Applying (2.4), we formulate the second term on the right-hand side of (3.11) as follows

    (3.17)

    The first term on the right-hand side of (3.17) is solved by

    The calculations in (2.12) imply that

    (3.18)

    Substituting (3.18) into the second term on the right-hand side of (3.17), we have

    (3.19)

    For the third term on the right-hand side of (3.17), we have

    (3.20)

    We use (2.4) to formulate the first term on the right-hand side of (3.20) as follows

    (3.21)

    For the first term on the right-hand side of (3.21), we have

    The second term on the right-hand side of (3.21) can be solved as (3.19). We employ the analogous method utilized for the first two terms on the right-hand side of (3.4) to solve the last two terms on the right-hand side of (3.21).

    For the second term on the right-hand side of (3.20), we have

    For the last term on the right-hand side of (3.17), by Lemmas 2 and 6, we have

    Combining all the estimates, we conclude that

    where is some positive constant. By the smallness of , , and the standard continuity method, we arrive at the main result.

    All authors contributed equally.

    The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.

    The work of the first author was partially supported by the National Natural Science Foundation of China under Grants 12401278. The work of the second author was partially supported by the RFS grant and GRF grants from the Research Grants Council (Project Nos. PolyU 11302021, 11310822, and 11302523). The authors would like to thank the research center for nonlinear analysis at PolyU for the opportunity of discussions and encouragement.

    The authors declare there is no conflict of interest.



    [1] T. C. Sideris, B. Thomases, Global existence for three-dimensional incompressible isotropic elastodynamics via the incompressible limit, Comm. Pure Appl. Math., 58 (2005), 750–788. https://doi.org/10.1002/cpa.20049 doi: 10.1002/cpa.20049
    [2] T. C. Sideris, B. Thomases, Global existence for three-dimensional incompressible isotropic elastodynamics, Comm. Pure Appl. Math., 60 (2007), 1707–1730. https://doi.org/10.1002/cpa.20196 doi: 10.1002/cpa.20196
    [3] F. John, Blow-up for quasilinear wave equations in three spaces dimensions, Comm. Pure Appl. Math., 34 (1981), 29–51. https://doi.org/10.1002/cpa.3160340103 doi: 10.1002/cpa.3160340103
    [4] A. S. Tahvildar-Zadeh, Relativistic and nonrelativistic elastodynamics with small shear strains, Ann. Inst. H. Poincaré Phys. Théor., 69 (1998), 275–307.
    [5] F. John, Almost global existence of elastic waves of finite amplitude arising from small initial disturbance, Comm. Pure Appl. Math., 41 (1988), 615–666. https://doi.org/10.1002/cpa.3160410507 doi: 10.1002/cpa.3160410507
    [6] S. Klainerman, T. C. Sideris, On almost global existence for nonrelativistic wave equations in 3D, Comm. Pure Appl. Math., 49 (1996), 307–321. https://doi.org/10.1002/(SICI)1097-0312(199603)49:3%3C307::AID-CPA4%3E3.0.CO;2-H doi: 10.1002/(SICI)1097-0312(199603)49:3%3C307::AID-CPA4%3E3.0.CO;2-H
    [7] T. C. Sideris, The null condition and global existence of nonlinear elastic waves, Invent. Math., 123 (1996), 323–342. https://doi.org/10.1007/s002220050030 doi: 10.1007/s002220050030
    [8] T. C. Sideris, Nonresonance and global existence of prestressed nonlinear elastic waves, Ann. of Math., 151 (2000), 849–874. https://doi.org/10.2307/121050 doi: 10.2307/121050
    [9] R. Agemi, Global existence of nonlinear elastic waves, Invent. Math., 142 (2000), 225–250. https://doi.org/10.1007/s002220000084 doi: 10.1007/s002220000084
    [10] Z. Lei, T. C. Sideris, Y. Zhou, Almost global existence for 2-D incompressible isotropic elastodynamics, Trans. Amer. Math. Soc., 367 (2015), 8175–8197. https://doi.org/10.1090/tran/6294 doi: 10.1090/tran/6294
    [11] Z. Lei, Global well-posedness of incompressible elastodynamics in two dimensions, Comm. Pure Appl. Math., 69 (2016), 2072–2106. https://doi.org/10.1002/cpa.21633 doi: 10.1002/cpa.21633
    [12] S. Alinhac, The null condition for quasilinear wave equations in two space dimensions Ⅰ, Invent. Math., 145 (2001), 597–618. https://doi.org/10.1007/s002220100165 doi: 10.1007/s002220100165
    [13] X. Wang, Global existence for the 2D incompressible isotropic elastodynamics for small initial data, Ann. Henri Poincaré., 18 (2017), 1213–1267. https://doi.org/10.1007/s00023-016-0538-x doi: 10.1007/s00023-016-0538-x
    [14] X. Cui, S. Yin, Global existence of inhomogeneous incompressible elastodynamics in three dimensions, SIAM J. Math. Anal., 50 (2018), 4721–4751. https://doi.org/10.1137/17M1143113 doi: 10.1137/17M1143113
    [15] Y. Cai, Uniform bound of the higest-order energy of the 2D incompressible elastodynamics, SIAM J. Math. Anal., 55 (2023), 5893–5918. https://doi.org/10.1137/22M1519894 doi: 10.1137/22M1519894
    [16] Z. Lei, F. Wang, Uniform bound of the highest energy for the three dimensional incompressible elastodynamics, Arch. Ration. Mech. Anal., 216 (2015), 593–622. https://doi.org/10.1007/s00205-014-0815-0 doi: 10.1007/s00205-014-0815-0
    [17] C. Wang, X. Yu, Global existence of null- form wave equations on small asymptotically Euclidean manifolds, J. Funct. Anal., 266 (2014), 5676–5708. https://doi.org/10.1016/j.jfa.2014.02.028 doi: 10.1016/j.jfa.2014.02.028
    [18] H. Hidano, C. Wang, K. Yokoyama, On almost global existence and local well posedness for some 3-D quas-ilinear wave equations, Adv. Differential Equations, 17 (2012), 267–306. https://doi.org/10.57262/ade/1355703087 doi: 10.57262/ade/1355703087
    [19] J. Metcalfe, C. Sogge, Long-time existence of quasilinear wave equations exterior to star-sharped obstacles via energy methods, SIAM J. Math. Anal., 38 (2006), 188–209. https://doi.org/10.1137/050627149 doi: 10.1137/050627149
  • Reader Comments
  • © 2025 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(230) PDF downloads(28) Cited by(0)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog