Processing math: 100%
Research article

Attractors for a quasilinear viscoelastic equation with nonlinear damping and memory

  • Received: 24 August 2020 Accepted: 11 October 2020 Published: 21 October 2020
  • MSC : 35L72, 35B41, 35B35

  • In this paper, the long time behavior of a quasilinear viscoelastic equation with nonlinear damping is considered. Under suitable assumptions, the existence of global attractors is established.

    Citation: Xiaoming Peng, Yadong Shang. Attractors for a quasilinear viscoelastic equation with nonlinear damping and memory[J]. AIMS Mathematics, 2021, 6(1): 543-563. doi: 10.3934/math.2021033

    Related Papers:

    [1] Jiangwei Zhang, Yongqin Xie . Asymptotic behavior for a class of viscoelastic equations with memory lacking instantaneous damping. AIMS Mathematics, 2021, 6(9): 9491-9509. doi: 10.3934/math.2021552
    [2] Peipei Wang, Yanting Wang, Fei Wang . Indirect stability of a 2D wave-plate coupling system with memory viscoelastic damping. AIMS Mathematics, 2024, 9(7): 19718-19736. doi: 10.3934/math.2024962
    [3] Abdelbaki Choucha, Salah Boulaaras, Asma Alharbi . Global existence and asymptotic behavior for a viscoelastic Kirchhoff equation with a logarithmic nonlinearity, distributed delay and Balakrishnan-Taylor damping terms. AIMS Mathematics, 2022, 7(3): 4517-4539. doi: 10.3934/math.2022252
    [4] Keltoum Bouhali, Sulima Ahmed Zubair, Wiem Abedelmonem Salah Ben Khalifa, Najla ELzein AbuKaswi Osman, Khaled Zennir . A new strict decay rate for systems of longitudinal $ m $-nonlinear viscoelastic wave equations. AIMS Mathematics, 2023, 8(1): 962-976. doi: 10.3934/math.2023046
    [5] Abdelbaki Choucha, Muajebah Hidan, Bahri Cherif, Sahar Ahmed Idris . Growth of solutions with $ L^{2(p+2)} $-norm for a coupled nonlinear viscoelastic Kirchhoff equation with degenerate damping terms. AIMS Mathematics, 2022, 7(1): 371-383. doi: 10.3934/math.2022025
    [6] Qian Li . General and optimal decay rates for a viscoelastic wave equation with strong damping. AIMS Mathematics, 2022, 7(10): 18282-18296. doi: 10.3934/math.20221006
    [7] Xiaoming Peng, Xiaoxiao Zheng, Yadong Shang . Lower bounds for the blow-up time to a nonlinear viscoelastic wave equation with strong damping. AIMS Mathematics, 2018, 3(4): 514-523. doi: 10.3934/Math.2018.4.514
    [8] Tingting Liu, Tasneem Mustafa Hussain Sharfi, Qiaozhen Ma . Time-dependent asymptotic behavior of the solution for evolution equation with linear memory. AIMS Mathematics, 2023, 8(7): 16208-16227. doi: 10.3934/math.2023829
    [9] Adel M. Al-Mahdi, Mohammad M. Al-Gharabli, Maher Nour, Mostafa Zahri . Stabilization of a viscoelastic wave equation with boundary damping and variable exponents: Theoretical and numerical study. AIMS Mathematics, 2022, 7(8): 15370-15401. doi: 10.3934/math.2022842
    [10] Qian Li, Yanyuan Xing . General and optimal decay rates for a system of wave equations with damping and a coupled source term. AIMS Mathematics, 2024, 9(10): 29404-29424. doi: 10.3934/math.20241425
  • In this paper, the long time behavior of a quasilinear viscoelastic equation with nonlinear damping is considered. Under suitable assumptions, the existence of global attractors is established.


    In this paper we investigate the long-time dynamics of solutions for the quasilinear viscoelastic equation with nonlinear damping and memory

    {|ut|ρuttΔuttαΔu+tμ(ts)Δu(s)ds+f(u)+g(ut)=h(x),u|Ω=0,u(x,0)=u0(x),ut(x,0)=u1(x), (1.1)

    where Ω is a bounded domain of RN(N1) with smooth boundary Ω, u0 is the prescribed past history of u.

    Problem (1.1) can be seen as an extension, accounting for memory effects in the material, of equations of the form

    f(ut)uttΔuΔutt=0. (1.2)

    This equation is interesting not only from the point of view of PDE general theory, but also due to its applications in Mechanics. In the case f(ut) is a constant, Eq (1.2) has been used to model extensional vibrations of thin rods (see Love [1,Chapter 20]). In the case f(ut) is not a constant, Eq (1.2) can model materials whose density depends on the velocity ut. We refer the reader to Fabrizio and Morro [2] for several other related models.

    When ρ=0 and Δutt is dropped in (1.1), the related problem has been extensively studied and several results concerning the global existence, decay of global solution and finite time blow up have been established. In this direction, we refer the readers to see Ref. [3,4,5,6,7,8,9,10,11,12] and the references therein.

    Now let us recall some results concerning quasilinear viscoelastic wave equations. In [13], Cavalcanti et al. studied the following equation with Dirichlet boundary conditions

    |ut|ρuttΔuttΔu+t0g(ts)Δu(s)dsγΔut=0. (1.3)

    A global existence result for γ0, as well as an exponential decay for γ>0, has been established. This last result has been extended by Messaoudi and Tatar [14] to the case γ=0.

    In [15], Messaoudi and Tatar studied the following equation

    |ut|ρuttΔuttΔu+t0g(ts)Δu(s)ds=b|u|p2u. (1.4)

    By introducing a new functional and using a potential well method, they obtained the global existence of solutions and the uniform decay of the energy if the initial data are in some stable set. In the case b=0 in (1.4), Messaoudi and Tatar [16] proved the exponential decay of global solutions to (1.4), without smallness of initial data, considering only the dissipation effect given by the memory. Liu [17] proved that for certain initial data in the stable set, the solution decays exponentially, and for certain initial data in the unstable set, the solution blows up in finite time.

    Replacing strong damping by weak damping in Eq (1.3), several authors have studied the energy decay rates of the related problems like

    |ut|ρuttΔuttΔu+t0g(ts)Δu(s)ds+h(ut)=0. (1.5)

    When h(ut)=ut, Han and Wang [18] investigated the global existence and exponential stability of the energy for solutions for Eq (1.5). When h(ut)=|ut|mut(m>0), the same authors [19] proved the general decay of energy for Eq (1.5). Later, Park and Park [20] established the general decay for Eq (1.5) with general nonlinear weak damping.

    Now, we list some important literature on the nonlinear evolution equation with hereditary memory and variable density. Araújo et al. [21] studied the following equation

    |ut|ρuttΔuttαΔu+tμ(ts)Δu(s)dsγΔut+f(u)=h(x), (1.6)

    and proved the global existence, uniqueness and exponential stability of solutions and existence of the global attractor. Subsequently, Qin et al. [22,23] proved the upper semicontinuity of pullback attractors and the existence of uniform attractors by assuming f(u)=0 and taking a frictional damping ut instead of strong damping Δut. However, their argument for uniqueness rely on the differentiability of the map σ(x)=|x|ρ at zero, which introduces the further restriction ρ>1.

    Lately, the authors [24] established an existence, uniqueness and continuous dependence result for the weak solutions to the semigroup generated for the system (1.6) in a three-dimensional space when ρ[0,4] and f has polynomial growth of (at most) critical order 5. Then, based on the [24], the same authors [25] established the existence of the global attractor of optimal regularity for Eq (1.6) when ρ[0,4). Recently, the authors [26] proved that the sole weak dissipation (γ=0) given by the memory term is enough to ensure existence and optimal regularity of the global attractor. Leuyacc and Parejas [27] proved the upper semicontinuity of global attractors when ρ0+ in (1.6). Li and Jia [28] proved the existence of a global solution by means of the Galerkin method, establish the exponential stability result and the polynomial stability result when the kernel μ(s) satisfies μ(s)k1μq(s), 1q<3/2.

    Motivated by the works above mentioned, our aim is to present the existence of global attractors for the problem (1.1).

    As in[29,30,31], we shall introduce a new variable ηt to the system which corresponds to the relative displacement history.

    Let us define

    ηt(x,s)=u(x,t)u(x,ts),(x,s)Ω×R+, t0.

    Note that

    ηtt(x,s)=ηts(x,s)+ut(x,s).

    Thus, the original memory term can be rewritten as

    tμ(ts)Δu(s)ds=0μ(s)Δu(ts)ds=0μ(s)dsΔu0μ(s)Δηt(s)ds,

    and Eq (1.1) becomes

    |ut|ρuttΔutt(α0μ(s)ds)Δu0μ(s)Δηt(s)ds+f(u)+g(ut)=h(x).

    Assuming for simplicity that α0μ(s)ds=1, we have the new system

    {|ut|ρuttΔuttΔu0μ(s)Δηt(s)ds+f(u)+g(ut)=h(x),ηtt(x,s)=ηts(x,s)+ut(x,s),u|Ω=0,ηt|Ω=0,u(x,0)=u0(x),ut(x,0)=u1(x),η0(x,s)=η0(x,s), (1.7)

    where

    {u0(x)=u0(x,0),xΩ,u1(x)=tu0(x,t)|t=0,xΩ,η0(x,s)=u0(x,0)u0(x,s),(x,s)Ω×R+.

    We begin with precise hypotheses on problem (1.7). Assume that

    0<ρ<4N2ifN3andρ>0ifN=1,2, (2.1)

    Concerning the source term f:RR, we assume that

    f(0)=0,|f(u)f(v)|c0(1+|u|p+|v|p)|uv|, u,vR, (2.2)

    where c0>0 and

    0<p<4N2ifN3andp>0ifN=1,2. (2.3)

    In addition, we assume that

    f(u)uF(u)0, uR, (2.4)

    where F(u)=u0f(s)ds.

    The damping function gC1(R) is a non-decreasing function with g(0)=0 and satisfies the polynomial condition

    g(s)0,|g(u)g(v)|c1(1+|u|q+|v|q)|uv|, u,vR, (2.5)

    where c1>0 and

    0<q4N2ifN3andq>0ifN=1,2. (2.6)

    With respect to the memory component, we assume that

    μC1(R+)L1(R+),μ(s)0,0μ(s)<, (2.7)

    and there exist k0,k1>0 such that

    0μ(s)ds=k0, (2.8)

    and

    μ(s)+k1μ(s)0, sR+. (2.9)

    As usual, p denotes the Lp-norms as well as (,) denotes either the L2-inner product. Let λ1>0 be the first eigenvalue of Δ in H10(Ω).

    In order to consider the relative displacement ηt as a new variable, one introduces the weighted L2space

    M=L2μ(R+;H10(Ω))={ξ:R+H10(Ω)|0μ(s)ξ(s)22ds<},

    which is a Hilbert space endowed with inner product and norm

    (ξ,ζ)M=0μ(s)(ξ,ζ)dsandξ2M=0μ(s)ξ22ds,

    respectively.

    Next let us introduce the phase space

    H=H10(Ω)×H10(Ω)×M,

    endowed with the norm

    zH=(u,v,η)2H=u22+v22+η2M.

    Then, the energy of problem (1.7) is given by

    E(t)=1ρ+2utρ+2ρ+2+12u22+12ut22+12ηt2M+Ω(F(u)hu)dx. (2.10)

    According to the arguments [24] with slightly modified, we can obtain the following well-posedness result.

    Theorem 2.1. Assume the assumptions (2.1)-(2.7) hold. If initial data z0=(u0,u1,η0)H and hL2(Ω), then the problem (1.7) admits a unique global solution

    z=(u,ut,ηt)C([0,T],H), (2.11)

    satisfying

    uL(R+;H10(Ω)),utL(R+;H10(Ω)),uttL(R+;H10(Ω)),ηtL(R+;M).

    Remark 1. The well-posedness of problem (1.7) given by Theorem 2.1 implies that the one-parameter family of operators S(t):HH defined by

    S(t)z0=(u(t),ut(t),ηt(t))=z,t0, (2.12)

    where z=(u(t),ut(t),ηt(t)) is the unique weak solution of the system (1.7), satisfies the semigroup properties

    S(0)=IandS(t+s)=S(t)S(s),t,s0,

    and defines a nonlinear C0-semigroup. Then problem (1.7) can be viewed as a nonlinear infinite dynamical system (H,S(t)).

    Now we give the following result concerning the global attractors.

    Theorem 2.2. Assume the assumptions (2.1)-(2.7) hold and hL2(Ω). Then the dynamical system (H,S(t)) generated by (1.7) has a compact global attractor AH.

    Before presenting our results we recall some fundamentals of the theory of infinite-dimensional dynamical systems which can be founded in the book by Chueshov and Lasiecka [32,33].

    Theorem 3.1. A dissipative dynamical system (X,S(t)) has a compact global attractor if and only if it is asymptotically smooth.

    The proof of asymptotic smoothness property can be very delicate. Here we use the following “compensated compactness” result [32,34].

    Theorem 3.2. Let (X,S(t)) be a dynamical system on a complete metric space X endowed with a metric d. Suppose that for any bounded positively invariant set BX and for any ε>0, there exists T=T(ε,B) such that

    S(T)xS(T)yXε+ΦT(x,y), x,yB,

    where ΦT:B×BR satisfies

    lim infnlim infmΦT(zn,zm)=0, (3.1)

    for any sequence {zn}nN in B. Then S(t) is asymptotic smooth in X.

    In the sequel we will apply the abstract results presented above to prove Theorem 2.2. Firstly, we show that the dynamical system (H,S(t)) is dissipative. The next step is to verify the asymptotic smoothness. Then the existence of a compact global attractor is guaranteed by Theorem 3.1. In what follows, the generic positive constants will be denoted as C, while Q() will stand for a generic increasing positive function.

    In this section, our aim is to show that the dynamical system (H,S(t)) is dissipative. To this aim, we first give some priori estimates used later.

    Proposition 3.1. For any initial data z0 with z0HR, we have the uniform estimate

    u22+ut22+ηt2M+utt22Q(R),t0.

    Proof. Multiplying (1.7) by (ut,ηt), we obtain

    E(t)=(g(ut),ut)(ηts,ηt)M. (3.2)

    Owing to (2.7), one can easily see that

    (ηts,ηt)M=12Ω(0μ(s)dds|ηt(s)|2ds)dx=12Ω(0μ(s)|ηt(s)|2ds)dx. (3.3)

    Combining (3.3) and (3.2), we have

    E(t)=(g(ut),ut)+120μ(s)ηt(s)22ds. (3.4)

    Since μ(s) is decreasing, g0 and g(0)=0, we get

    E(t)0,

    and consequently

    E(t)E(0).

    Applying Young inequality yields

    Ωhudx14u22+1λ1h22.

    It follows promptly from (2.10) that

    E(t)1ρ+2utρ+2ρ+2+14u22+12ut22+12ηt2M+ΩF(u)dx1λ1h22.

    Then making use of (2.4) we obtain

    1ρ+2utρ+2ρ+2+14u22+12ut22+12ηt2ME(t)+1λ1h22E(0)+1λ1h22.

    This means that

    utρ+2ρ+2+u22+ut22+ηt2MQ(R),t0. (3.5)

    A multiplication of (1.7) by utt gives

    Ω|ut|ρu2ttdx+utt22=Ωuuttdx(ηt,utt)MΩf(u)uttdxΩg(ut)uttdx+Ωhuttdx. (3.6)

    Next, we estimate each term individually. By Hölder inequality, Poincaré inequality and Young inequality, we have

    Ωuuttdx16utt22+32u22,
    (ηt,utt)M16utt22+3k02ηt2M,

    and

    Ωhuttdx16utt22+32λ1h22.

    Using Hölder inequality, Poincaré inequality, Young inequality, and (3.5), we have

    Ωf(u)uttdxC(up+2+up+1p+2)uttp+2C(u2+up+12)utt216utt22+Q(R),

    and

    Ωg(ut)uttdxC(utq+2+utq+1q+2)uttq+2C(ut2+utq+1q+2)utt216utt22+Q(R).

    Substituting all the above inequalities into (3.6) and taking (3.5) into account, we derive that

    Ω|ut|ρu2ttdx+16utt22Q(R).

    This means that

    utt22Q(R). (3.7)

    In light of (3.5) and (3.7), we obtain the Proposition 3.1.

    Now, we introduce the following two functionals

    Φ(t)=1ρ+1Ω|ut|ρutudx+Ωutudx, (3.8)

    and

    Ψ(t)=Ω(Δut1ρ+1|ut|ρut)0μ(s)ηt(s)dsdx. (3.9)

    Lemma 3.3. There exists a positive constant C1, such that

    Φ(t)E(t)14u22+C1ut22+2ρ+1utρ+2ρ+2C10μ(s)ηt(s)22ds. (3.10)

    Proof. A multiplication of the first equation of (1.7) by u gives

    Φ(t)=Ω(|ut|ρuttΔutt)udx+1ρ+1utρ+2ρ+2+ut22=u22+0μ(s)(ΩΔηt(s)u(t)dx)dsΩf(u)udxΩg(ut)udx+Ωhudx+1ρ+1utρ+2ρ+2+ut22.

    By Hölder inequality and Cauchy inequality, we obtain

    0μ(s)(ΩΔηt(s)u(t)dx)dsu(t)20μ(s)ηt2ds18u22+2k0ηt2M.

    Using Hölder inequality, Young inequality and Sobolev inequality, taking into account (2.5) and (2.6), we arrive at

    Ωg(ut)udxc1Ω(1+|ut|q)|ut||u|dxC(1+utqq+2)utq+2uq+2C(1+utq2)ut2u218u22+C(1+ut2q2)ut2218u22+Q(R)ut22.

    Combining the last two estimates, we end up with

    Φ(t)1ρ+1utρ+2ρ+234u22+Cut22+2k0ηt2MΩf(u)udx+Ωhudx.

    Besides, in light of (2.9), we get

    ηt2M1k10μ(s)ηt(s)22ds. (3.11)

    Using (2.10), (2.4) and (3.11) yields

    Φ(t)E(t)14u22+2ρ+1utρ+2ρ+2+C1ut22C10μ(s)ηt(s)22ds.

    This lemma is complete.

    Lemma 3.4. For any δ1>0, there exists C2>0 such that

    Ψ(t)δ1C2u(t)22(k0δ1C2)ut(t)22k0ρ+1ut(t)ρ+2ρ+2+k0h22C20μ(s)ηt22ds. (3.12)

    Proof. Taking the time derivative of Ψ(t), in light of the first equation of (1.7), we get

    Ψ(t)=Ω(|ut|ρutt+Δutt)(0μ(s)ηt(s)ds)dx+Ω(|ut|ρutρ+1+Δut)(0μ(s)ηtt(s)ds)dx=Ω(Δu0μ(s)Δηt(s)ds+g(ut)+f(u)h)(0μ(s)ηt(s)ds)dx+Ω(|ut|ρutρ+1+Δut)(0μ(s)ηtt(s)ds)dx.

    Next, we will estimate the right side of the above identity. Integrating by parts with respect to x and using Young inequality, we obtain

    ΩΔu(0μ(s)ηt(s)ds)dx=Ωu(0μ(s)ηt(s)ds)dxδ1u22+k04δ1ηt2M,

    and

    Ω(0μ(s)Δηt(s)ds)(0μ(s)ηt(s)ds)dx=ΩNj=1(0μ(s)ηt(s)xjds)2dxk0ΩNj=1(0μ(s)|ηt(s)xj|2ds)dx=k0ηt2M.

    Applying (2.5), Hölder inequality, Sobolev embedding inequality, Young inequality and Proposition 3.1, we obtain

    Ωg(ut)(0μ(s)ηt(s)ds)dx=0μ(s)(Ωg(ut)ηt(s)dx)dsC0μ(s)(1+ut(t)qq+2)ut(t)q+2ηt(s)q+2dsC0(1+ut(t)q2)ut(t)2μ(s)ηt2dsδ1Q(R)ut22+14δ1ηt2M.

    Analogously, but using (2.2) instead of (2.5), we have

    Ωf(u)(0μ(s)ηt(s)ds)dx=0μ(s)(Ωf(u)ηt(s)dx)dsC0μ(s)(1+u(t)pp+2)u(t)p+2ηt(s)p+2dsC0μ(s)(1+u(t)p2)u(t)2ηt2dsδ1Q(R)u22+14δ1ηt2M.

    Using Hölder inequality, Young's inequality and Sobolev embedding inequality, we get

    Ωh(0μ(s)ηt(s)ds)dx0μ(s)h2ηt(s)2dsk0h22+Cηt(s)2M.

    On the other hand, since

    0μ(s)ηtt(s)ds=0μ(s)ηts(s)ds+0μ(s)ut(t)ds=0μ(s)ηt(s)ds+k0ut(t),

    we find

    Ω(|ut|ρutρ+1+Δut)(0μ(s)ηtt(s)ds)dx=k0ut(t)22k0ρ+1ut(t)ρ+2ρ+2+0μ(s)(ΩΔut(t)ηt(s)dx)ds1ρ+10μ(s)(Ω|ut(t)|ρut(t)ηt(s))ds.

    Applying Hölder inequality, Young's inequality and Sobolev embedding inequality, we obtain

    0μ(s)(ΩΔut(t)ηt(s)dx)ds0μ(s)ut(t)2ηt(s)2dsδ1ut(t)22μ(0)4δ10μ(s)ηt(s)22ds,

    and

    1ρ+10μ(s)(Ω|ut(t)|ρut(t)ηt(s))ds1ρ+10μ(s)utρ+1ρ+2ηtρ+2dsδ1μ(0)ρ+1ut2(ρ+1)ρ+214δ1(ρ+1)0μ(s)ηt2ρ+2dsδ1Cut2(ρ+1)2C0μ(s)ηt22dsδ1Q(R)ut22C0μ(s)ηt22ds.

    Collecting all the above inequalities and (3.11), we end up with

    Ψ(t)δ1C2u(t)22(k0δ1C2)ut(t)22k0ρ+1ut(t)ρ+2ρ+2+C2h22C20μ(s)ηt22ds.

    This lemma is complete.

    Lemma 3.5 (Absorbing set). Under the hypotheses of Theorem 2.2, the dynamical system (H,S(t)) corresponding to problem (1.7) has a bounded absorbing set.

    Proof. Let us consider a perturbed functional

    L(t)=ME(t)+εΦ(t)+Ψ(t) (3.13)

    where Φ(t) and Ψ(t) are the same functional defined in (3.8) and (3.9).

    Since

    Ωhudx14u22+1λ1h22,

    we get

    E(t)1ρ+2utρ+2ρ+2+14u22+12ut22+12ηt2M+ΩF(u)dx1λ1h22.

    Then making use of (2.4) we obtain

    14(u22+ut22+ηt2M)E(t)+1λ1h22. (3.14)

    Now we claim that there exist three constants α1,α2,b>0 such that

    α1E(t)bh22L(t)α2E(t)+bh22. (3.15)

    To prove this, we first note that

    |Φ(t)|C0(ut22+u22)4C0(E(t)+1λ1h22),
    |Ψ(t)|C0(ut22+u22)4C0(E(t)+1λ1h22).

    Hence there exists a constant b>0 such that

    |εΦ(t)+Ψ(t)|b(E(t)+h22).

    Then taking M>b, we get (3.15) with α1=Mb and α2=Mb.

    Combining (3.4), (3.13), Lemma 3.3 and 3.4, we arrive at

    L(t)=ME(t)+εΦ(t)+Ψ(t)εE(t)(k0δ1C2εC1)ut22(ε4δ1C2)u22+k0h22+[M2(εC1+C2)]0μ(s)ηt22ds.

    Choosing ε,δ1>0 small enough and M>0 sufficiently large such that

    k0δ1C2εC1>0,ε4δ1C2>0,M2(εC1+C2)>0,

    then we have

    L(t)εE(t)+k0h22,

    which together with (3.15) implies that

    L(t)εα2L(t)+(bεα2+k0)h22.

    Integrating the above inequality over [0,t], we can derive

    L(t)L(0)eεα2t+(1eεα2t)(b+α2k0ε)h22.

    Using again (3.15) yields

    E(t)α2α1E(0)eεα2t+1α1(2b+α2k0ε)h22.

    Recalling (3.14), we obtain

    (u,ut,ηt)2H4α2α1E(0)eεα2t+R20,

    where

    R20=1α1(2b+α2k0ε)h22+4bh22.

    This shows that any closed ball ˉB(0,R) with R>R0 is a bounded absorbing set of (H,S(t)). The proof of Lemma 3.5 is now complete.

    As a straight consequence of Lemma 3.5, we have that the solutions of problem (1.7) are globally bounded provided initial data lying in bounded sets BH. Namely, let z=(u,ut,ηt) be a solution of (1.7) with initial data z0=(u0,u1,η0) in a bounded set B. Then one has

    zHCB,t0, (3.16)

    where CB is a constant depending on B. Lemma 3.5 also ensures the existence of bounded positively invariant sets.

    Lemma 3.6 (Stabilizability). Under the hypotheses of Theorem (2.2), given a bounded set BH, let z1=(u,ut,ηt) and z2=(v,vt,ξt) be two weak solutions of problem (1.7) such that z1(0)=(u0,u1,η0) and z2(0)=(v0,v1,ξ0) are in B. Then

    z1(t)z2(t)2HCBeγtz1(0)z2(0)2H+CBt0eγ(tτ)(w(τ)ρ+2+w(τ)p+2+wt(τ)ρ+2+wt(τ)p+2)dτ,

    where w=uv and γ,CB are positive constants depending on B.

    Proof. Let us write w=uv and ζt=ηtξt. Then (w,ζt) satisfies

    {Δw+Δwtt+0μ(s)Δζt(s)ds=|ut|ρutt|vt|ρvtt+f(u)f(v)+g(ut)g(vt),ζtt=ζts+wt,w(0)=u0v0,wt(0)=u1v1,ζ0=η0ξ0. (3.17)

    Now we consider the functional

    Ew(t)=12wt22+12w22+12ζt2M (3.18)

    and its perturbation

    G(t)=MEw(t)+εϕ(t)+ψ(t),ε>0, (3.19)

    where

    ϕ(t)=ΩΔwt(t)w(t)dx,
    ψ(t)=Ω[Δwt1ρ+1(|ut|ρut|vt|ρvt)](0μ(s)ζt(s)ds)dx.

    We divide the remaining of the proof into five steps. Hereafter, we use CB to denote several positive constants.

    Step 1. For M>0 sufficiently large, there exists ε0 such that

    M2Ew(t)G(t)3M2Ew(t),t0,ε(0,ε0]. (3.20)

    To prove this, we first observe that

    |ϕ(t)|12w22+12wt22Ew(t).

    Besides, using Hölder inequality, Sobolev inequality, Young inequality and (3.16), we can derive that

    |ψ(t)|0μ(s)wt2ζt(s)2ds+2ρ0μ(s)(utρρ+2+vtρρ+2)wtρ+2ζtρ+2dsk02wt22+12ζt2M+CB0μ(s)(utρ2+vtρ2)wt2ζt2dsk02wt22+12ζt2M+CBwt22+12ζt2MCBEw(t).

    Collecting the above two estimates and (3.19), we obtain

    (MCBε)Ew(t)G(t)(M+CB+ε)Ew(t).

    Now let us put ε0=M/2CB. Then for all εϵ0, the inequality (3.20) holds.

    Step 2. There exists a constant C3>0 such that

    Ew(t)120μ(s)ζt(s)22ds+C3(wtρ+2+wtp+2). (3.21)

    In fact, differentiating (3.18) with respect to t and using (3.17), we have

    Ew(t)=Ω(|vt|ρvtt|ut|ρutt)wtdxΩ(f(u)f(v))wtdxΩ(g(ut)g(vt))wtdx+120μ(s)ζt(s)22ds. (3.22)

    Using Hölder inequality, Sobolev inequality, estimate (3.7) and (3.16), we have

    Ω(|vt|ρvtt|ut|ρutt)wtdxΩ|vt|ρ|vtt||wt|dx+Ω|ut|ρ|utt||wt|dx(vtρρ+2vttρ+2+utρρ+2uttρ+2)wtρ+2CB(vtρ2vtt2+utρ2utt2)wtρ+2CBwtρ+2.

    Combining (2.2), Hölder inequality, Sobolev inequality, and estimate (3.16) yields

    Ω(f(u)f(v))wtdxCB(1+upp+2+vpp+2)wp+2wtp+2CB(1+up2+vp2)w2wtp+2CBwtp+2.

    Since g is a non-decreasing function, this yields

    Ω(g(ut)g(vt))wtdx0.

    Inserting last three estimates into (3.22), we arrive at

    Ew(t)120μ(s)ζt(s)22ds+CB(wtρ+2+wtp+2).

    Step 3. There exists C4>0 such that

    ϕ(t)Ew(t)14w22+C4(wρ+2+wp+2)+C4wt22C40μ(s)ζt(s)22ds. (3.23)

    Taking the derivative of ϕ(t), it follows from (3.17) that

    ϕ(t)=ΩΔwttwdx+wt22=Ω(|vt|ρvtt|ut|ρutt)wdxΩ(f(u)f(v))wdxΩ(g(ut)g(vt))wdx+0μ(s)ΩΔζt(s)wdxdsw22+wt22. (3.24)

    From Hölder inequality, Sobolev embedding, estimate (3.7) and (3.16), we have

    Ω(|vt|ρvtt|ut|ρutt)wdxΩ|vt|ρ|vtt||w|dx+Ω|ut|ρ|utt||w|dx(vtρρ+2vttρ+2+utρρ+2uttρ+2)wρ+2CB(vtρ2vtt2+utρ2utt2)wρ+2CBwρ+2.

    Applying (2.2), Hölder inequality, Sobolev embedding, and estimate (3.16), we get

    Ω(f(u)f(v))wdx(1+upp+2+vpp+2)wp+2wp+2CB(1+up2+vp2)w2wp+2CBwp+2.

    Similarly, in light of (2.5), Hölder inequality, Young inequality, Sobolev embedding, and estimate (3.16), we obtain

    Ω(g(ut)g(vt))wdxCB(1+utqq+2+vqq+2)wtq+2wq+2CBwt2w218w22+CBwt22.

    Using Young inequality gives

    0μ(s)ΩΔζt(s)wdxds18w22+2k0ζt2M.

    Combining these six last estimates with (3.18) we end up with

    ϕ(t)Ew(t)14w22+CBwt22+(12+2k0)ζtM+CB(wρ+2+wp+2),

    which together with (3.11) implies that inequality (3.23) holds for some C4>0.

    Step 4. For any δ2>0, there exists C5>0 such that

    ψ(t)k04wt22+2δ2w22C50μ(s)ζt22ds. (3.25)

    Taking derivative of ψ(t) and using (3.17), we derive that

    ψ(t)=ΩΔw0μ(s)ζt(s)dsdxΩ0μ(s)Δζt(s)ds0μ(s)ζt(s)dsdx+Ω(f(u)f(v))0μ(s)ζt(s)dsdx+Ω(g(ut)g(vt))0μ(s)ζt(s)dsdx+ΩΔwt0μ(s)ζtt(s)dsdx1ρ+1Ω(|ut|ρut|vt|ρvt)0μ(s)ζtt(s)dsdx=6i=1Ai. (3.26)

    Applying Hölder inequality, Young inequality, Sobolev embedding, estimate (3.7) and (3.16), we get

    A1δ2w22+k04δζt2M, (3.27)
    A2k0ζt2M, (3.28)
    A33p(1+upp+2+vpp+2)wp+20μ(s)ζt(s)p+2dsCB(1+up2+vp2)w20μ(s)ζt(s)2dsδ2w22+CB4δ2ζt2M, (3.29)
    A42q(1+uqq+2+vqq+2)wtq+20μ(s)ζt(s)q+2dsCB(1+uq2+vq2)wt20μ(s)ζt(s)2dsk04wt22+CBζt2M. (3.30)

    From (3.17), one can easily see that

    0μ(s)ζtt(s)ds=0μ(s)ζts(s)ds+0μ(s)wtds=k0wt+0μ(s)ζt(s)ds.

    Hence in light of Young inequality we obtain

    A5k0wt220μ(s)wt2ζt(s)2ds3k04wt22μ(0)k00μ(s)ζt(s)22ds. (3.31)

    By the monotonicity of function x|x|ρx (ρ>0), we get

    1ρ+1Ω(|ut|ρut|vt|ρvt)wtdx0.

    Using Hölder inequality, Young inequality and estimate (3.16), we obtain

    A62ρ0μ(s)(utρρ+2+vtρρ+2)wt(t)ρ+2ηt(s)ρ+2ds2ρ0μ(s)(utρ2+vtρ2)wt(t)2ηt(s)2dsμ(0)4wt(t)22CB0μ(s)ηt(s)22ds. (3.32)

    Inserting (3.27)-(3.32) into (3.26), (3.11) we end up with

    ψ(t)k04wt22+2δ2w22C50μ(s)ζt22ds.

    Step 5. Combining (3.21), (3.23), (3.25) with (3.19), we have

    G(t)=MEw(t)+εϕ(t)+ψ(t)εEw(t)(k02εC4)wt22(ε42δ2)w22+(M2εC4C5)0μ(s)ζt22ds+C(B,M,ε)(wρ+2+wp+2+wtρ+2+wtp+2). (3.33)

    Firstly we fix ε>0 such that εC4<k0/2. Then taking δ2>0 such that δ2<ε/8. For fixed ε and δ2, we choose M>0 so large that M>2(εC4+C5). Then (3.33) along with (3.20) give

    G(t)εEw(t)+CB(wρ+2+wp+2+wtρ+2+wtp+2)2ε3MG(t)+CB(wρ+2+wp+2+wtρ+2+wtp+2). (3.34)

    Integrating (3.34) over (0,t) with respect to t, we get

    G(t)G(0)e2ε3Mt+CBt0e2ε3M(tτ)(wρ+2+wp+2+wtρ+2+wtp+2)dτ,

    which together with (3.20) implies that

    Ew(t)3Ew(0)eγt+CBt0eγ(tτ)(wρ+2+wp+2+wtρ+2+wtp+2)dτ,

    where γ=2ε/3M is a positive constant. Notice that the functional Ew(t) is equivalent to the norm of H, the proof is complete.

    Lemma 3.7 (Asymptotic smoothness). Under the hypotheses of Theorem 2.2, the dynamical system corresponding to problem (1.7) is asymptotic smooth.

    Proof. Let B be a bounded subset of H positively invariant with respect to S(t). Let S(t)z1(0)=(u,ut,ηt) and S(t)z2(0)=(v,vt,ξt) be two solutions for problem (1.7) corresponding to initial data z1(0),z2(0)B. Given ε>0, we can choose T>0 so large that CBeγt<ε. We claim that there exists constant CBT>0 such that

    z1z2Hε+ΦT(z1(0),z2(0)), z1(0),z2(0)B, (3.35)

    with

    ΦT(z10,z20)=CBT(T0(u(τ)v(τ)2ρ+2+u(τ)v(τ)2p+2+ut(τ)vt(τ)2ρ+2+ut(τ)vt(τ)2p+2)dτ)12. (3.36)

    Indeed, from Lemma 3.6, we have

    z1(T)z2(T)HCBeγT+CB(T0e2γ(tτ)dτ)12(T0(u(τ)v(τ)2ρ+2+u(τ)v(τ)2p+2+ut(τ)vt(τ)2ρ+2+ut(τ)vt(τ)2p+2)dτ)12CBeγT+CBT(T0(u(τ)v(τ)2ρ+2+u(τ)v(τ)2p+2+ut(τ)vt(τ)2ρ+2+ut(τ)vt(τ)2p+2)dτ)12,

    and consequently (3.35) and (3.36) hold.

    We are left to prove that ΦT satisfies (3.1). Indeed, given a sequence of initial data zn=(un0,un1,ηn0)B, we write S(t)zn=(un(t),unt(t),ηn,t). Since B is invariant by S(t), t0, it follows that (un(t),unt(t),ηn,t) uniformly bounded in H. Namely,

    (un,unt,ηn,t) is bounded in C([0,T];H10(Ω)×H10(Ω)×M),T>0.

    Then by compact embedding H10(Ω)Lρ+2(Ω) and H10(Ω)Lp+2(Ω), there exists a subsequence (un,unt,ηn,t) such that

    un and untconverges strongly in C([0,T];Lρ+2(Ω));un and untconverges strongly in C([0,T];Lp+2(Ω)).

    Therefore,

    limnlimmT0(un(τ)um(τ)2ρ+2+unt(τ)umt(τ)2ρ+2+un(τ)um(τ)2p+2+unt(τ)umt(τ)2p+2)dτ=0,

    which implies (3.1) holds. Then asymptotic smoothness follows from Theorem 3.2.

    Proof of Theorem 2.2. We first note that Lemmas 3.5 and 3.7 imply that (H,S(t)) is a dissipative dynamical system which is asymptotically smooth. Then the existence of a compact global attractor A to problem (1.7) in the phase space H follows from Theorem 3.1.

    The authors are grateful to the referees for the constructive comments and kind suggestions.

    All authors declare no conflicts of interest in this paper.



    [1] A. E. H. Love, A treatise on the mathematical theory of elasticity, New York: Dover, 1944.
    [2] M. Fabrizio, A. Morro, Mathematical problems in linear viscoelasticity, Philadelphia: SIAM, 1992.
    [3] J. E. Muñoz Rivera, Asymptotic behaviour in linear viscoelasticity, Quart. Appl. Math., 52 (1994), 628-648.
    [4] J. E. Muñoz Rivera, E. C. Lapa, R. Barreto, Decay rates for viscoelastic plates with memory, J. Elasticity, 44 (1996), 61-87. doi: 10.1007/BF00042192
    [5] M. Aassila, M. M. Cavalcanti, J. A. Soriano, Asymptotic stability and energy decay rates for solutions of the wave equation with memory in a star-shaped domain, SIAM J. Control Optim., 38 (2000), 1581-1602. doi: 10.1137/S0363012998344981
    [6] M. M. Cavalcanti, V. N. D. Cavalcanti, T. F. Ma, J. A. Soriano, Global existence and asymptotic stability for viscoelastic problems, Differential Integral Equations, 15 (2002), 731-748.
    [7] M. M. Cavalcanti, H. P. Oquendo, Frictional versus viscoelastic damping in a semilinear wave equation, SIAM J. Control Optim., 42 (2003), 1310-1324. doi: 10.1137/S0363012902408010
    [8] A. Guesmia, S. A. Messaoudi, A general decay result for a viscoelastic equation in the presence of past and finite history memories, Nonlinear Anal. Real World Appl., 13 (2012), 476-485. doi: 10.1016/j.nonrwa.2011.08.004
    [9] S. A. Messaoudi, General decay of solutions of a viscoelastic equation, J. Math. Anal. Appl., 341 (2008), 1457-1467. doi: 10.1016/j.jmaa.2007.11.048
    [10] S. A. Messaoudi, Blow up and global existence in a nonlinear viscoelastic wave equation, Math. Nachr., 260 (2003), 58-66. doi: 10.1002/mana.200310104
    [11] S. A. Messaoudi, Blow-up of positive-initial-energy solutions of a nonlinear viscoelastic hyperbolic equation, J. Math. Anal. Appl., 320 (2006), 902-915. doi: 10.1016/j.jmaa.2005.07.022
    [12] J. Y. Park, J. R. Kang, Global attractor for hyperbolic equation with nonlinear damping and linear memory, Sci. China Math., 53 (2010), 1531-1539. doi: 10.1007/s11425-010-3110-z
    [13] M. M. Cavalcanti, V. N. D. Cavalcanti, J. Ferreira, Existence and uniform decay for a non-linear viscoelastic equation with strong damping, Math. Methods Appl. Sci., 24 (2001), 1043-1053. doi: 10.1002/mma.250
    [14] S. A. Messaoudi, N. Tatar, Exponential and polynomial decay for a quasilinear viscoelastic equation, Nonlinear Anal., 68 (2008), 785-793. doi: 10.1016/j.na.2006.11.036
    [15] S. A. Messaoudi, N. Tatar, Global existence and uniform stability of solutions for a quasilinear viscoelastic problem, Math. Methods Appl. Sci., 30 (2007), 665-680. doi: 10.1002/mma.804
    [16] S. A. Messaoudi, N. Tatar, Exponential decay for a quasilinear viscoelastic equation, Math. Nachr., 282 (2009), 1443-1450. doi: 10.1002/mana.200610800
    [17] W. J. Liu, General decay and blow-up of solution for a quasilinear viscoelastic problem with nonlinear source, Nonlinear Anal., 73 (2010), 1890-1904. doi: 10.1016/j.na.2010.05.023
    [18] X. S. Han, M. X. Wang, Global existence and uniform decay for a nonlinear viscoelastic equation with damping, Nonlinear Anal., 70 (2009), 3090-3098. doi: 10.1016/j.na.2008.04.011
    [19] X. S. Han, M. X. Wang, General decay of energy for a viscoelastic equation with nonlinear damping, Math. Methods Appl. Sci., 32 (2009), 346-358. doi: 10.1002/mma.1041
    [20] J. Y. Park, S. H. Park, General decay for quasiliear viscoelastic equations with nonlinear weak damping, J. Math. Phys., 50 (2009), 083505. doi: 10.1063/1.3187780
    [21] R. O. Araújo, T. F. Ma, Y. M. Qin, Long-time behavior of a quasilinear viscoelastic equation with past history, J. Differ. Equations, 254 (2013), 4066-4087. doi: 10.1016/j.jde.2013.02.010
    [22] Y. M. Qin, B. W. Feng, M. Zhang, Uniform attractors for a non-autonomous viscoelastic equation with a past history, Nonlinear Anal., 101 (2014), 1-15. doi: 10.1016/j.na.2014.01.006
    [23] Y. M. Qin, J. P. Zhang, L. L. Sun, Upper semicontinuity of pullback attractors for a nonautonomous viscoelastic equation, Appl. Math. Comput., 223 (2013), 362-376.
    [24] M. Conti, E. M. Marchini, V. Pata, A well posedness result for nonlinear viscoelastic equations with memory, Nonlinear Anal., 94 (2014), 206-216. doi: 10.1016/j.na.2013.08.015
    [25] M. Conti, E. M. Marchini, V. Pata, Global attractors for nonlinear viscoelastic equations with memory, Commun. Pure Appl. Anal., 15 (2016), 1893-1913. doi: 10.3934/cpaa.2016021
    [26] M. Conti, T. F. Ma, E. M. Marchini, P. N. Seminario Huertas, Asymptotics of viscoelastic materials with nonlinear density and memory effects, J. Differ. Equations, 264 (2018), 4235-4259.
    [27] Y. R. S. Leuyacc, J. L. C. Parejas, Upper semicontinuity of global attractors for a viscoelastic equations with nonlinear density and memory effects, Math. Methods Appl. Sci., 42 (2019), 871-882. doi: 10.1002/mma.5389
    [28] F. S. Li, Z. Q. Jia, Global existence and stability of a class of nonlinear evolution equations with hereditary memory and variable density, Bound. Value Probl., 2019 (2019), 37. doi: 10.1186/s13661-019-1152-x
    [29] C. M. Dafermos, Asymptotic stability in viscoelasticity, Arch. Ration. Mech. Anal., 37 (1970), 297-308. doi: 10.1007/BF00251609
    [30] C. Giorgi, J. E. Muñoz Rivera, V. Pata, Global attractors for a semilinear hyperbolic equation in viscoelasticity, J. Math. Anal. Appl., 260 (2001), 83-99. doi: 10.1006/jmaa.2001.7437
    [31] V. Pata, A. Zucchi, Attractors for a damped hyperbolic equation with linear memory, Adv. Math. Sci. Appl., 11 (2001), 505-529.
    [32] I. Chueshov, I. Lasiecka, Von Karman evolution equations, New York: Springer-Verlag, 2010.
    [33] I. Chueshov, Dynamics of quasi-stable dissipative systems, New York: Springer, 2015.
    [34] I. Chueshov, I. Lasiecka, Long-time behavior of second oreder evolution equations with nonlinear damping, Mem. Amer. Math. Soc., 912 (2008), 912.
  • This article has been cited by:

    1. Jiangwei Zhang, Zhiming Liu, Jianhua Huang, Upper semicontinuity of optimal attractors for viscoelastic equations lacking strong damping, 2022, 0003-6811, 1, 10.1080/00036811.2022.2088532
    2. 玉晶 任, Random Attractors of Viscoelastic Equation with Strong Damping, 2022, 12, 2160-7583, 1597, 10.12677/PM.2022.1210174
  • Reader Comments
  • © 2021 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(3749) PDF downloads(176) Cited by(2)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog