This paper dealt with a laminated beam system along with structural damping, past history, distributed delay, and in the presence of both temperatures and micro-temperatures effects. The damping terms left the system dissipative. Employing the semigroup approach, we established the existence and uniqueness of the solution. Additionally, with the help of convenient assumptions on the kernel, we demonstrated a general decay result for the solution of the considered system, with no constraints regarding the speeds of wave propagation. The main aim was to address how specific behaviors of the system were related to memory and delays. We aimed to investigate the joint impact of an infinite memory, distributed delay and micro-temperature effects on the system. We found a new relationship between the decay rate of solution and the growth of g at infinity. The objective was to find studies that use no- trivial results and their applications to relevant problems from mathematical physics.
Citation: Abdelkader Moumen, Fares Yazid, Fatima Siham Djeradi, Moheddine Imsatfia, Tayeb Mahrouz, Keltoum Bouhali. The influence of damping on the asymptotic behavior of solution for laminated beam[J]. AIMS Mathematics, 2024, 9(8): 22602-22626. doi: 10.3934/math.20241101
[1] | Tijani A. Apalara, Aminat O. Ige, Cyril D. Enyi, Mcsylvester E. Omaba . Uniform stability result of laminated beams with thermoelasticity of type Ⅲ. AIMS Mathematics, 2023, 8(1): 1090-1101. doi: 10.3934/math.2023054 |
[2] | Hicham Saber, Fares Yazid, Fatima Siham Djeradi, Mohamed Bouye, Khaled Zennir . Decay for thermoelastic laminated beam with nonlinear delay and nonlinear structural damping. AIMS Mathematics, 2024, 9(3): 6916-6932. doi: 10.3934/math.2024337 |
[3] | Fatima Siham Djeradi, Fares Yazid, Svetlin G. Georgiev, Zayd Hajjej, Khaled Zennir . On the time decay for a thermoelastic laminated beam with microtemperature effects, nonlinear weight, and nonlinear time-varying delay. AIMS Mathematics, 2023, 8(11): 26096-26114. doi: 10.3934/math.20231330 |
[4] | Cyril Dennis Enyi, Soh Edwin Mukiawa . Dynamics of a thermoelastic-laminated beam problem. AIMS Mathematics, 2020, 5(5): 5261-5286. doi: 10.3934/math.2020338 |
[5] | Adel M. Al-Mahdi, Maher Noor, Mohammed M. Al-Gharabli, Baowei Feng, Abdelaziz Soufyane . Stability analysis for a Rao-Nakra sandwich beam equation with time-varying weights and frictional dampings. AIMS Mathematics, 2024, 9(5): 12570-12587. doi: 10.3934/math.2024615 |
[6] | Soh Edwin Mukiawa . Well-posedness and stabilization of a type three layer beam system with Gurtin-Pipkin's thermal law. AIMS Mathematics, 2023, 8(12): 28188-28209. doi: 10.3934/math.20231443 |
[7] | Houssem Eddine Khochemane, Ali Rezaiguia, Hasan Nihal Zaidi . Exponential stability and numerical simulation of a Bresse-Timoshenko system subject to a neutral delay. AIMS Mathematics, 2023, 8(9): 20361-20379. doi: 10.3934/math.20231038 |
[8] | Hasan Almutairi, Soh Edwin Mukiawa . On the uniform stability of a thermoelastic Timoshenko system with infinite memory. AIMS Mathematics, 2024, 9(6): 16260-16279. doi: 10.3934/math.2024787 |
[9] | Soh E. Mukiawa, Tijani A. Apalara, Salim A. Messaoudi . Stability rate of a thermoelastic laminated beam: Case of equal-wave speed and nonequal-wave speed of propagation. AIMS Mathematics, 2021, 6(1): 333-361. doi: 10.3934/math.2021021 |
[10] | Said Mesloub, Hassan Altayeb Gadain, Lotfi Kasmi . On the well posedness of a mathematical model for a singular nonlinear fractional pseudo-hyperbolic system with nonlocal boundary conditions and frictional damping terms. AIMS Mathematics, 2024, 9(2): 2964-2992. doi: 10.3934/math.2024146 |
This paper dealt with a laminated beam system along with structural damping, past history, distributed delay, and in the presence of both temperatures and micro-temperatures effects. The damping terms left the system dissipative. Employing the semigroup approach, we established the existence and uniqueness of the solution. Additionally, with the help of convenient assumptions on the kernel, we demonstrated a general decay result for the solution of the considered system, with no constraints regarding the speeds of wave propagation. The main aim was to address how specific behaviors of the system were related to memory and delays. We aimed to investigate the joint impact of an infinite memory, distributed delay and micro-temperature effects on the system. We found a new relationship between the decay rate of solution and the growth of g at infinity. The objective was to find studies that use no- trivial results and their applications to relevant problems from mathematical physics.
Mathematical modeling is indispensable in engineering, natural science, and applied mathematics to capture the effects of both memory and delay ingrained in the studied actualities. To this end, the inclusion of both of them is often simplified for presentation purposes, as a specific description of basic operations can be intricate for mathematical manipulation. A key question to address is how certain behaviors are related to memory and delays. In this study, we investigate the joint impact of an infinite memory, distributed delay, and micro-temperature effects on the system (1.1).
In the current work, we study the following thermoelastic laminated beam, together with structural damping, infinite memory, distributed delay, and micro-temperatures effects:
{ϱϖtt+G(ϕ−ϖx)x+γθx=0,Iϱ(3ψ−ϕ)tt−D(3ψ−ϕ)xx−G(ϕ−ϖx)−mθ+drx+∫∞0g(s)(3ψ−ϕ)xx(x,t−s)ds=0,3Iϱψtt−3Dψxx+3G(ϕ−ϖx)+4δψ+4βψt+4∫ς2ς1|μ2(ς)|ψt(x,t−ς)dς=0,cθt−k0θxx+m(3ψ−ϕ)t+γϖtx+k1rx=0,αrt−k2rxx+k3r+k1θx+d(3ψ−ϕ)tx=0, | (1.1) |
where
(x,ς,t)∈(0,1)×(ς1,ς2)×R+, |
and the initial and boundary conditions are given by
{ϖ(x,0)=ϖ0,ψ(x,0)=ψ0,ϕ(x,0)=ϕ0,θ(x,0)=θ0,r(x,0)=r0,x∈(0,1),ϖt(x,0)=ϖ1,ψt(x,0)=ψ1,ϕt(x,0)=ϕ1,x∈(0,1),ϖx(0,t)=ϕ(0,t)=ψ(0,t)=θ(0,t)=r(0,t)=0,t>0,ϖ(1,t)=ϕx(1,t)=ψx(1,t)=θx(1,t)=r(1,t)=0,ψt(x,−t)=f0(x,t)t>0. | (1.2) |
Here, ϖ denotes the transverse displacement, ϕ represents the rotation angle, ψ is relative to the amount of slip occurring along the interface, θ is the temperature difference and r is the micro-temperature vector. The coefficients δ,β,ϱ,Iϱ,G, and D, are positive constants representing the adhesive stiffness, the adhesive damping parameter, the density, the shear stiffness, the flexible rigidity and the mass moment of inertia, respectively. We denote by the positive constants c,k0,k1,k2,k3,d,γ,α,m, the physical parameters describing the coupling between the various constituents of the materials.
Herein, ς1,ς2 are positive numbers such that 0<ς1<ς2, and μ2 is an L∞ function satisfying the following assumption:
● The function μ2:[ς1,ς2]→R is bounded and it fulfills
β−∫ς2ς1|μ2(ς)|dς>0. |
To motivate our work, let us recall some earlier related results. For the problems with the Timoshenko system with/without thermal law, one can see the works [4,6,7,16,23,24] and for problems related to thermoelasticity, we mention for instance [8,10,11,13,17,18].
We start with the laminated beam model, which has become quite popular, and both scientists and engineers are interested in it. This model is a pertinent study topic, because of the wide industry applicability of such materials. Hansen and Spies in [12] were the first to introduce the following beam with two layers by developing this mathematical model:
{ρ1ϖtt+G(ϕ−ϖx)x=0,ρ2(3ψ−ϕ)tt−G(ϕ−ϖx)−D(3ψ−ϕ)xx=0,ρ3ψtt+G(ϕ−ϖx)+43γψ+43βψt−Dψxx=0. | (1.3) |
The laminated beam equations have produced some results so far, most of which are focused on the system's stability and existence. Provided that the assumption of equal wave speeds holds, it was demonstrated that system (1.3) is exponentially stable, when linear damping terms are incorporated in two of the three equations. However, if they are included in the three equations, then the system decays exponentially with no restriction on the speeds of wave propagations, see, for instance [1,22].
Lately, a renewed focus on investigating the asymptotic behavior of the solutions of several thermoelastic laminated beams has grown. For more details about this topic the reader may consult [2,9,20].
The thermoelastic laminated beam problem together with nonlinear weights and time-varying delay was the study topic of Nonato et al. in [20], where the authors considered two cases (with and without the structural damping) and proved an exponential decay result for both of them. Distributed delay is one of the main damping factors in our model. It is used to model systems in which there is a delay of uncertain duration. The physical interpretation of this term differs from the delayed differential equation, as it can take several values. For example, in incoming signals, distributed delay shortens the setup and lengthens the hold time. Even moderate distributed delay likely makes setup time negative on those inputs that are directly connected to the register.
The infinite memory is a critical aspect in addressing problems, and it has been explored in various contexts such as the work of Liu and Zhao [14], in which they considered a thermoelastic laminated beam model with past history. The authors managed to establish both exponential and polynomial stabilities, depending on the kernel function for the system involving structural damping and with no constraint on the wave speeds. Moreover, concerning the system in the absence of structural damping, they were able to establish both exponential and polynomial stabilities, in case of equal wave speeds and lack of exponential stability in the opposite case.
The time delays problems are one of the most significant and active research areas recently. Numerous studies have demonstrated that delay can lead to instability unless certain conditions are incorporated, and it also can lead to distinct solutions that differ from those found in prior studies. Therefore, the issue of stability for systems that involve delay is highly crucial. To learn more about this term, we refer the reader to the following papers [3,5,21].
In [19], Nicaise and Pignotti made a study on the following wave equation, together with linear frictional damping and internal distributed delay:
utt−Δu+μ1ut+a(x)∫τ2τ1μ2(s)ut(t−s)ds, in Ω×(0,∞), |
and assuming that
‖a‖∞∫τ2τ1μ2(s)ds<μ1, |
the authors managed to prove that the solution is exponentially stable.
Problem (1.1) is considered as a delayed system, and it is also called hereditary systems, posteffect systems, and deviating argument. Distributed delay is a physical phenomenon which is found in a multitude of applications: Many real systems whose temporal evolution is not defined from a simple vector of state (expressed in the present tense) but depends irreducibly on the history of the system. This situation is encountered in the cases-numerous-where a transport of matter, energy or information generates a "dead time" in the reaction: in information and communication technologies (high-speed communication networks, control of networked systems, quality of service in Moving Picture Experts Group (MPEG) video transmissions, tele-operated systems, parallel computing, realtime computing in robotics), in population dynamics and epidemiology (gestation or incubation time), and in mechanics (viscoelasticity). Even if the process does not intrinsically contain a post-effect, its control chain can introduce distributed delays (for example, if the sensors require a significant acquisition/transmission time). For these reasons, it seems reasonable to consider distributed delay as a universal characteristic of the interaction between man and nature (hence, of sciences for engineers). The aim of our study then concerns the interaction between the different damping terms which intervene in the qualitative properties of the energy associated to the system. Before this analysis, we must ensure the existence of unique solution and then we can pass to see the asymptotic behavior of the solution with respect to damping terms. We used classical semigroup theory to find nontrivial results regarding the well-posedness of solutions. Then, under minimal restrictions on the kernel, we found qualitative properties of the solution by contracting an appropriate Lyapnov functional. The main goal is to present fundamental and new techniques for modern models applying science and technology that can stimulate research interest for exploration of mathematical applications in real life sciences.
The rest of the current paper is structured this way: In Section 2, we provide some resources required for our research, then highlight our major results. In Section 3, we establish the well-posedness of the system. In Section 4, we introduce some fundamental lemmas required in the proof later. In Section 5, we demonstrate our general decay result.
In this section, we provide some materials required in the proof later, then state our major results.
● (A1) Let g:R+→R+ be a C1 function which satisfies
g(0)>0,D−g0=ˉl>0, where g0:=∫∞0g(s)ds. | (2.1) |
● (A2) There exists a strictly increasing convex function G:R+→R+ of class C1(R+)∩C2(]0,+∞[) which satisfies
{G(0)=G′(0)=0,limt→+∞G′(t)=+∞, |
such that
sups∈R+∫∞0g(s)G−1(−g′(s))ds+∫∞0g(s)G−1(−g′(s))ds<+∞. |
Now, we present the following useful inequalities.
Lemma 2.1. The following inequalities are valid,
∫10[∫∞0g(s)((3ψ−ψ)(t)−(3ψ−ϕ)(t−s))ds]2dx≤c1(g⋄(3ψ−ϕ)x)(t), | (2.2) |
∫10[∫∞0g′(s)((3ψ−ϕ)x(t)−(3ψ−ϕ)x(t−s))ds]2dx≤−g(0)(g′⋄(3ψ−ϕ)x)(t), | (2.3) |
∫10[∫∞0g(s)((3ψ−ϕ)x(t)−(3ψ−ϕ)x(t−s))ds]2dx≤g0(g⋄(3ψ−ϕ)x)(t), | (2.4) |
∫10[∫∞0g′(s)((3ψ−ϕ)(t)−(3ψ−ϕ)(t−s))ds]2dx≤−c2(g′⋄(3ψ−ϕ)x)(t), | (2.5) |
where c1,c2>0, and
(g⋄v)(t)=∫10∫∞0g(s)(v(x,t)−v(x,t−s))2dsdx. |
Let us start by introducing (see [19])
{ηt(x,s)=(3ψ−ϕ)(x,t)−(3ψ−ϕ)(x,t−s),S(x,p,ς,t)=ψt(x,t−ςp), | (2.6) |
where
(x,p,ς,s,t)∈((0,1))2×(ς1,ς2)×R+×R+. |
Then, the variables ηt and S surely satisfy
{ηtt+ηts=(3ψ−ϕ)t,ςSt(x,p,ς,t)+Sp(x,p,ς,t)=0,S(x,0,ς,t)=ψt(x,t). | (2.7) |
Hence, system (1.1) can be rewritten as
{ϱϖtt+G(ϕ−ϖx)x+γθx=0,Iϱ(3ψ−ϕ)tt−D(3ψ−ϕ)xx−G(ϕ−ϖx)−mθ+drx+∫∞0g(s)(3ψ−ϕ)xx(t−s)ds=0,3Iϱψtt−3Dsxx+3G(ϕ−ϖx)+4δψ+4βψt+4∫ς2ς1|μ2(ς)|S(x,1,ς,t)dς=0,cθt−k0θxx+m(3ψ−ϕ)t+γϖtx+k1rx=0,αrt−k2rxx+k3r+k1θx+d(3ψ−ϕ)tx=0,ηtt+ηts=(3ψ−ϕ)t,ςSt+Sp=0. | (2.8) |
Certainly, system (2.8) is depending on the initial and boundary conditions below
{ϖ(x,0)=ϖ0,ψ(x,0)=ψ0,ϕ(x,0)=ϕ0,θ(x,0)=θ0,r(x,0)=r0,x∈(0,1),ϖt(x,0)=ϖ1,ψt(x,0)=ψ1,ϕt(x,0)=ϕ1,x∈(0,1),ϖx(0,t)=ϕ(0,t)=ψ(0,t)=θ(0,t)=r(0,t)=0,t>0,ϖ(1,t)=ϕx(1,t)=ψx(1,t)=θx(1,t)=r(1,t)=0,ψt(x,−t)=f0(x,t)t>0,ηt(0,s)=ηtx(1,s)=0,ηt(x,0)=0,η0(x,s)=η0(x,s),t,s>0,S(x,p,ς,0)=f0(x,pς),x,p∈(0,1),ς∈(ς1,ς2),t,s>0. | (2.9) |
Now, let
{ζ=3ψ−ϕ,ζ(0,t)=ζx(1,t)=0,ζ(x,0)=ζ0,ζt(x,0)=ζ1,(x,t)∈(0,1)×R+. |
Then, system (2.8) is equivalent to
{ϱϖtt+G(3ψ−ζ−ϖx)x+γθx=0,Iϱζtt−Dζxx−G(3ψ−ζ−ϖx)−mθ+drx+∫∞0g(s)ζxx(t−s)ds=0,3Iϱψtt−3Dψxx+3G(3ψ−ζ−ϖx)+4δψ+4βψt+4∫ς2ς1|μ2(ς)|S(x,1,ς,t)dς=0,cθt−k0θxx+mζt+γϖtx+k1rx=0,αrt−k2rxx+k3r+k1θx+dζtx=0,ηtt+ηts=ζt,ςSt+Sp=0. | (2.10) |
Taking advantage of (2.6), we can rewrite the second equation of (2.10) as
Iϱζtt−ˉlζxx−G(3ψ−ζ−ϖx)−mθ+drx−∫∞0g(s)ηtxx(x,s)ds=0. |
At this step, let us introduce the vector function U=(ϖ,u,ζ,ν,ψ,y,θ,r,ηt,S)T, with
u=ϖt,ν=ζt,y=ψt, |
then, system (2.10) becomes
{ddtU(t)=AU(t),t>0,U(0)=U0=(ϖ0,ϖ1,ζ0,ζ1,ψ0,ψ1,θ0,r0,η0,f0)T, | (2.11) |
here, A:D(A)⊂H:→H stands for a linear operator indicated by
AU=(u−1ϱ(G(3ψ−ζ−ϖx)x+γθx)ν1Iϱ(ˉlζxx+G(3ψ−ζ−ϖx)+mθ−drx+∫∞0g(s)ηtxx(x,s)ds)y1Iϱ(Dψxx−G(3ψ−ζ−ϖx)−43δψ−43βy−43∫ς2ς1|μ2(ς)|S(x,1,ς,t)dς)1c(k0θxx−mν−γux−k1rx)1α(k2rxx−k3r−k1θx−dνx)ν−ηts−1ςSp). |
Now, we shall consider the ensuing energy space
H=˜J1∗(0,1)×L2(0,1)×J1∗(0,1)×L2(0,1)×J1∗(0,1)×L2(0,1)×L2(0,1)×L2(0,1)×Lg×L2((0,1)×(0,1)×(ς1,ς2)), |
where
J1∗(0,1)={φ∈H1(0,1):φ(0)=0},˜J1∗(0,1)={φ∈H1(0,1):φ(1)=0},J2∗(0,1)=H2(0,1)∩J1∗(0,1),˜J2∗(0,1)=H2(0,1)∩˜J1∗(0,1), |
and
Lg={φ:R+→J1∗(0,1),∫10∫∞0g(s)φ2xdsdx<∞}. |
For the space Lg, we take the following inner product
⟨φ1,φ2⟩Lg=∫10∫∞0g(s)φ1xφ2xdsdx. |
Furthermore, we consider the following domain
Lg(R+,J1∗(0,1))={ηt∈Lg,ηts∈Lg,ηt(x,0)=0}. |
Then, we introduce
⟨U,ˉU⟩H=ϱ∫10uˉudx+Iϱ∫10νˉνdx+3Iϱ∫10yˉydx+c∫10θˉθdx+α∫10rˉrdx+ˉl∫10ζxˉζxdx+G∫10(3ψ−ζ−ϖx)(3ˉψ−ˉζ−ˉϖx)dx+4δ∫10ψˉψdx+3D∫10ψxˉψxdx+∫10∫10g(s)ηtx(x,t)ˉηtx(x,s)dsdx+4∫10∫10∫ς2ς1ς|μ2(ς)|SˉSdςdpdx. | (2.12) |
We deduce that H together with (2.12) is a Hilbert space, once we do that, we define D(A) by
D(A)={U∈H:ϖ∈˜J2∗(0,1);ζ,ψ∈J2∗(0,1);u∈˜J1∗(0,1);ν,y∈J1∗(0,1),θ∈J1∗(0,1),θt∈L2(0,1);r∈H2(0,1)∩H10(0,1),ηt∈Lg(R+,J1∗(0,1));S,Sp∈L2((0,1)×(0,1)×(ς1,ς2)),S(x,0,ς,t)=yϖx(0,t)=ζx(1,t)=ψx(1,t)=θx(1,t)=ηtx(1,s)=0.}. |
Obviously, D(A) is dense in H.
Now, we are ready to state our results.
Theorem 2.1. Let U0∈D(A), then problem (2.9)-(2.10) admits a unique solution
U∈C(R+,D(A))∩C1(R+,H). |
In addition, if U0∈H, then
U∈C(R+,H). |
We give the energy of the solution of problem (2.8)-(2.9) by
E(t)=12∫10{ϱϖ2t+G(ϕ−ϖx)2+Iϱ(3ψt−ϕt)2+ˉl(3ψx−ϕx)2+3Iϱψ2t+3Dψ2x+4δψ2+cθ2+αr2}dx+12(g⋄(3ψ−ϕ)x)(t)+2∫10∫10∫ς2ς1ς|μ2(ς)|S2(x,p,ς,t)dςdpdx. | (2.13) |
Then, we have the following stability result.
Theorem 2.2. Let (ϖ,ϕ,ψ,θ,r,ηt,S) be the solution of (2.8)-(2.9), suppose that (T), (A1) and (A2) hold. Then, for any initial data U0∈D(A) satisfying, for some p0≥0,
∫10η20x(x,s)dx≤p0,for alls>0, | (2.14) |
there exist positive constants α1,α2, and α3, such that
E(t)≤α1G−1∗(α2t+α3), | (2.15) |
where
G−1∗(t)=∫∞tdsG0(s),G0(t)=tG′(ϵ0t),for allϵ0≥0. |
In this part, we utilize the semigroup approach to prove our well-posedness result.
Proof of Theorem 2.1. Let's us establish the dissipativity of A. By (2.12) and for any U∈D(A), we have
⟨AU,U⟩H=−4β∫10y2dx−k3∫10r2dx−k2∫10r2xdx−k0∫10θ2xdx−4∫10∫ς2ς1|μ2(ς)|yS(x,1,ς,t)dςdx−4∫10∫10∫ς2ς1|μ2(ς)|SpSdςdpdx+12(g′⋄ζx)(t)≤0. |
One can notice that
−4∫10∫10∫ς2ς1|μ2(ς)|SpSdςdpdx=−2∫10∫ς2ς1∫10|μ2(ς)|∂pS2dpdςdx=−2∫10∫ς2ς1|μ2(ς)|S2(x,1,ς,t)dςdx+2∫10∫ς2ς1|μ2(ς)|S2(x,0,ς,t)dςdx. | (3.1) |
Applying Youg's inequality, we obtain
−4∫10∫ς2ς1|μ2(ς)|yS(x,1,ς,t)dςdx≤2(∫ς2ς1|μ2(ς)|dς)∫10y2dx+2∫10∫ς2ς1|μ2(ς)|S2(x,1,ς,t)dςdx, |
therefore, by (T) and given S(x,0)=y, we end up with
⟨AU,U⟩H=−4(β−∫ς2ς1|μ2(ς)|dς)∫10y2dx−k3∫10r2dx−k2∫10r2xdx−k0∫10θ2xdx+12(g′⋄ζx)(t)≤0. |
Thereby, A is dissipative.
Thereafter, we establish the surjectivity of (I−A), that is, we show that
∀H=(h1,h2,h3,h4,h5,h6,h7,h8,h9,h10)T∈H,∃U∈D(A):(I−A)U=H. | (3.2) |
Now, we have
{ϖ−u=h1,ϱu+G(3ψ−ζ−ϖx)x+γθx=ϱh2,ζ−ν=h3,Iϱν−ˉlζxx−G(3ψ−ζ−ϖx)−mθ+drx−∫∞0g(s)ηtxx(x,s)ds=Iϱh4,ψ−y=h5,3Iϱy−3Dψxx+3G(3ψ−ζ−ϖx)+4δψ+4βy+4∫ς2ς1|μ2(ς)|S(x,1,ς,t)dς=3Iϱh6,cθ−k0θxx+mν+γux+k1rx=ch7,(α+k3)r−k2rxx+k1θx+dνx=αh8,ηt−ν+ηts=h9,ςS+Sp=ςh10. | (3.3) |
Solving (3.3)10 and using S(x,0,ς,t)=y(x,t), we find
S(x,p,ς,t)=y(x,t)e−ςp+ςe−ςp∫p0eςσh10(x,σ,ς,t)dσ. |
Hence,
S(x,1,ς,t)=y(x,t)e−ς+ςe−ς∫10eςσh10(x,σ,ς,t)dσ. | (3.4) |
Now, we solve Eq (3.3)9, and we find
ηt=e−s∫s0eσ(ν+h9(σ))dσ. | (3.5) |
Inserting (3.5), (3.4), and
{u=ϖ−h1,ν=ζ−h3,y=ψ−h5, |
into (3.3)2, (3.3)4, (3.3)6, (3.3)7 and (3.3)8, we get
{ϱϖ+G(3ψ−ζ−ϖx)x+γθx=ϱ(h1+h2),Iϱζ−(ˉl+∫∞0(1−e−s)g(s)ds)ζxx−G(3ψ−ζ−ϖx)−mθ+drx=Iϱ(h3+h4)+˜h,μ1ψ−3Dψxx+3G(3ψ−ζ−ϖx)=˜μ1h5+3Iϱh6−4∫ς2ς1|μ2(ς)|ςe−ς∫10eςσh10dσdς,cθ−k0θxx+mζ+γϖx+k1rx=γh1x+ch7+mh3,(α+k3)r−k2rxx+k1θx+dζx=αh8+dh3x, | (3.6) |
where
˜h=∫∞0g(s)∫s0eσ−s(h9−h3)xxdσds, |
μ1=3Iϱ+4δ+4β+4∫ς2ς1e−ς|μ2(ς)|dς, |
and
˜μ1=3Iϱ+4β+4∫ς2ς1e−ς|μ2(ς)|dς. |
We take the following variational formulation to solve (3.6):
Q((ϖ,ζ,ψ,θ,r),(ˉϖ,ˉζ,ˉψ,ˉθ,ˉr))=L(ˉϖ,ˉζ,ˉψ,ˉθ)),∀(ˉϖ,ˉζ,ˉψ,ˉθ,ˉr))∈X, | (3.7) |
where,
X=˜J1∗(0,1)×J1∗(0,1)×J1∗(0,1)×L2(0,1)×H10(0,1), |
is a Hilbert space endowed with
‖(ϖ,ζ,ψ,θ,r)‖2X=‖3ψ−ζ−ϖx‖22+‖ϖ‖22+‖ζx‖22+‖ψx‖22+‖θx‖22+‖r‖22+‖rx‖22. |
As a part of this step, we provide definitions for both the bilinear form Q:X×X→R and the linear form L:X→R as follows:
Q((ϖ,ζ,ψ,θ,r),(ˉϖ,ˉζ,ˉψ,ˉθ,ˉr))=ϱ∫10ϖˉϖdx+Iϱ∫10ζˉζdx+μ1∫10ψˉψdx+c∫10θˉθdx+(α+k3)∫10rˉrdx+k2∫10rxˉrxdx+γ∫10(θxˉϖ+ϖxˉθ)dx+k0∫10θxˉθxdx+G∫10(3ψ−ζ−ϖx)(3ˉψ−ˉζ−ˉϖx)dx+(ˉl+∫∞0(1−e−s)g(s)ds)∫10ζxˉζxdx+3D∫10ψxˉψxdx+d∫10(rxˉζ+ζxˉr)dx+k1∫10(rxˉθ+ˉrθx)dx+m∫10(ζˉθ−ˉζθ)dx, | (3.8) |
and
L(ˉϖ,ˉζ,ˉψ,ˉθ,ˉr)=ϱ∫10ˉϖ(h1+h2)dx+Iϱ∫10ˉζ(h3+h4)dx+∫10ˉζ˜hdx+∫10ˉθ(γh1x+mh3+ch7)dx+∫10ˉr(αh8+dh3x)dx+∫10ˉψ[˜μ1h5+3Iϱh6−4∫ς2ς1|μ2(ς)|ςe−ς∫10eςσh10dσdς]dx. |
We can easily prove the continuity of Q and L. Moreover, from (3.8) together with integration by parts, we arrive at
Q((ϖ,ζ,ψ,θ,r),(ϖ,ζ,ψ,θ,r))=ϱ∫10ϖ2dx+Iϱ∫10ζ2dx+μ1∫10ψ2dx+c∫10θ2dx+(α+k3)∫10r2dx+k2∫10r2xdx+k0∫10θ2xdx+G∫10(3ψ−ζ−ϖx)2dx+3D∫10ψ2xdx+(ˉl+∫∞0(1−e−s)g(s)ds)∫10ζ2xdx≥M‖(ϖ,ζ,ψ,θ,r)‖2X,M>0. |
From this, we conclude the coercivity of Q. It follows from the Lax-Milgram lemma that (3.6) admits a unique solution satisfying
ϖ∈˜J1∗(0,1), |
ζ,ψ∈J1∗(0,1), |
θ∈L2(0,1), |
and
r∈H10(0,1). |
If we substitute ϖ,ζ, and ψ into (3.3)1, (3.3)3 and (3.3)5, we find
u∈˜J1∗(0,1), |
and
ν,y∈J1∗(0,1). |
In addition, taking (ˉζ,ˉψ,ˉθ,ˉr)≡(0,0,0,0)∈(J1∗(0,1))2×L2(0,1)×H10(0,1), (3.7) becomes
G∫10ˉϖϖxxdx=∫10ˉϖ(ϱϖ+3Gψx−Gζx+γθx−ϱ(h1+h2))dx, | (3.9) |
for all ˉϖ∈˜J1∗(0,1), which indicates that
Gϖxx=ϱϖ+3Gψx−Gζx+γθx−ϱ(h1+h2)∈L2(0,1). | (3.10) |
The standard elliptic regularity implies that
ϖ∈˜J2∗(0,1). |
We note that (3.9) remains valid for ˉφ∈C1([0,1])⊂˜J1∗(0,1), that is ˉφ(1)=0. Then, we obtain
G∫10ˉφxϖxdx=∫10ˉφ(−ϱϖ−3Gψx+Gζx−γθx+ϱ(h1+h2))dx. |
Integrating by parts, it follows that
ϖx(0)ˉφ(0)=0, for all ˉφ∈C1([0,1]). |
Hence
ϖx(0)=0. |
Likewise, we show that
(ζ,ψ)∈(J2∗(0,1))2,θ∈J1∗(0,1),r∈H2(0,1)∩H10(0,1), and ζx(1)=ψx(1)=θx(1)=0. |
The standard elliptic regularity guarantees the existence of a unique U∈D(A) which fulfills (3.2). Thereby, A is surjective.
As a consequence, we infer that A is a maximal dissipative operator. Then, the well-posedness result follows using Lumer-Philips theorem [15].
The main purpose of this section is to establish the essential practical lemmas required to prove our stability results. To attain this goal, we apply a specific approach known as the multiplier technique, which enables us to prove the stability results of problem (2.8). Nevertheless, this method necessitates creating an appropriate Lyapunov functional equivalent to the energy and we will clarify on this in the next section. To simplify matters, we will employ χ∗>0 to represent a generic constant.
Lemma 4.1. Let (ϖ,ϕ,ψ,θ,r,ηt,S) be the solution of (2.8) and (2.9), then, the energy functional satisfies
ddtE(t)≤−m0∫10ψ2tdx−k0∫10θ2xdx−k2∫10r2xdx−k3∫10r2dx+12(g′⋄(3ψx−ϕx))(t)≤0,wherem0>0. | (4.1) |
Proof. As a start, we multiply (2.8)1, (2.8)2, (2.8)3, (2.8)4 and (2.8)5 by ϖt,(3ψt−ϕt),ψt,θ and r respectively, then, we integrate over (0,1) and use integration by parts together with boundary conditions (2.9) and (2.6) to find
12ddt∫10{ϱϖ2t+G(ϕ−ϖx)2+Iϱ(3ψt−ϕt)2+ˉl(3ψx−ϕx)2+3Iϱψ2t+3Dψ2x+4δψ2+cθ2+αr2}dx+4β∫10ψ2tdx+k0∫10θ2xdx+k2∫10r2xdx+k3∫10r2dx−∫10(3ψ−ϕ)t∫∞0g(s)ηtxx(x,s)dsdx+4∫10∫ς2ς1ψt|μ2(ς)|S(x,1,ς,t)dςdx=0. | (4.2) |
It follows from the sixth equation in (2.8) and the integration by parts that
∫10(3ψ−ϕ)t∫∞0g(s)ηtxx(x,s)dsdx=∫∞0g(s)(∫10ηttηtxx(x,s)dx)ds+∫∞0g(s)(∫10ηtsηtxx(x,s)dx)ds=−12ddt(g⋄(3ψx−ϕx))(t)+12(g′⋄(3ψx−ϕx))(t). | (4.3) |
Applying Young's inequality, we find
∫10∫ς2ς1ψt|μ2(ς)|S(x,1,ς,t)dςdx≤12∫10∫ς2ς1|μ2(ς)|S2(x,1,ς,t)dςdx+12(∫ς2ς1|μ2(ς)|dς)∫10ψ2tdx. | (4.4) |
Next, we multiply (2.8)7 by S|μ2(ς)| and integrate the result over (0,1)×(0,1)×(ς1,ς2). We get
12ddt∫10∫10∫ς2ς1ς|μ2(ς)|S2(x,p,ς,t)dςdpdx=−∫10∫10∫ς2ς1|μ2(ς)|SpS(x,p,ς,t)dςdpdx=−12∫10∫10∫ς2ς1|μ2(ς)|∂pS2(x,p,ς,t)dςdpdx=12(∫ς2ς1|μ2(ς)|dς)∫10ψ2tdx−12∫10∫ς2ς1|μ2(ς)|S2(x,1,ς,t)dςdx, | (4.5) |
which, together with (4.2)–(4.4) and (T) gives us
ddtE(t)≤−4(β−∫ς2ς1|μ2(ς)|dς)∫10ψ2tdx−k0∫10θ2xdx−k2∫10r2xdx−k3∫10r2dx+12(g′⋄(3ψx−ϕx))(t)≤0. |
We have then reached the desired result.
Lemma 4.2. Consider the functional
I1(t):=−Iϱ∫10(3ψt−ϕt)∫∞0g(s)((3ψ−ϕ)(t)−(3ψ−ϕ)(t−s))dsdx, | (4.6) |
then, it satisfies
I′1(t)≤−Iϱg02∫10(3ψt−ϕt)2dx+ϵ1∫10(3ψx−ϕx)2dx+ϵ1∫10(ϕ−ϖx)2dx+ϵ1∫10θ2xdx+χ∗∫10r2dx+χ∗(1+1ϵ1)(g⋄(3ψx−ϕx))(t)−χ∗(g′⋄(3ψx−ϕx))(t),∀ϵ1>0. | (4.7) |
Proof. First, we notice that
∂t(∫∞0g(s)((3ψ−ϕ)(t)−(3ψ−ϕ)(t−s))ds)=∂t(∫t−∞g(t−s)((3ψ−ϕ)(t)−(3ψ−ϕ)(s))ds)=∫t−∞g′(t−s)((3ψ−ϕ)(t)−(3ψ−ϕ)(s))ds+∫t−∞g(t−s)(3ψ−ϕ)t(t)ds=∫∞0g′(s)((3ψ−ϕ)(t)−(3ψ−ϕ)(t−s))ds+g0(3ψ−ϕ)t(t). | (4.8) |
Next, we proceed by differentiating I1(t) and using both (2.8)2 and relation (4.8), then, integrating by parts, we get
F′1(t)=−Iϱ∫10(3ψtt−ϕtt)∫∞0g(s)((3ψ−ϕ)(t)−(3ψ−ϕ)(t−s))dsdx−Iϱ∫10(3ψt−ϕt)∂∂t(∫∞0g(s)((3ψ−ϕ)(t)−(3ψ−ϕ)(t−s))dsdx)=D∫10(3ψx−ϕx)∫∞0g(s)((3ψ−ϕ)x(t)−(3ψ−ϕ)x(t−s))dsdx−G∫10(ϕ−ϖx)∫∞0g(s)((3ψ−ϕ)(t)−(3ψ−ϕ)(t−s))dsdx−m∫10θ∫∞0g(s)((3ψ−ϕ)(t)−(3ψ−ϕ)(t−s))dsdx−Iϱg0∫10(3ψt−ϕt)2dx−Iϱ∫10(3ψt−ϕt)∫∞0g′(s)((3ψ−ϕ)(t)−(3ψ−ϕ)(t−s))dsdx−d∫10r∫∞0g(s)((3ψ−ϕ)x(t)−(3ψ−ϕ)x(t−s))dsdx−∫10(∫∞0g(s)(3ψ−ϕ)x(x,t−s)ds)×(∫∞0g(s)((3ψ−ϕ)x(t)−(3ψ−ϕ)x(t−s))ds)dx. | (4.9) |
The last term in (4.9) can be rewritten as
−∫10(∫∞0g(s)(3ψ−ϕ)x(x,t−s)ds)(∫∞0g(s)((3ψ−ϕ)x(t)−(3ψ−ϕ)x(t−s))ds)dx=∫10(∫∞0g(s)((3ψ−ϕ)x(t)−(3ψ−ϕ)x(t−s))ds)2dx−g0∫10(3ψ−ϕ)x(∫∞0g(s)((3ψ−ϕ)x(t)−(3ψ−ϕ)x(t−s))ds)dx. | (4.10) |
Now, replacing (4.10) into (4.9), leads to
F′1(t)=ˉl∫10(3ψx−ϕx)∫∞0g(s)((3ψ−ϕ)x(t)−(3ψ−ϕ)x(t−s))dsdx−G∫10(ϕ−ϖx)∫∞0g(s)((3ψ−ϕ)(t)−(3ψ−ϕ)(t−s))dsdx−m∫10θ∫∞0g(s)((3ψ−ϕ)(t)−(3ψ−ϕ)(t−s))dsdx−Iϱg0∫10(3ψt−ϕt)2dx−Iϱ∫10(3ψt−ϕt)∫∞0g′(s)((3ψ−ϕ)(t)−(3ψ−ϕ)(t−s))dsdx−d∫10r∫∞0g(s)((3ψ−ϕ)x(t)−(3ψ−ϕ)x(t−s))dsdx+∫10(∫∞0g(s)((3ψ−ϕ)x(t)−(3ψ−ϕ)x(t−s))ds)2dx. |
Finally, applying Young's inequality and making use of Lemma 2.1, we obtain (4.7).
Lemma 4.3. Consider the functional
I2(t):=−cϱ∫10ϖt(∫1xθ(y)dy)dx, |
then, it satisfies
I′2(t)≤−γϱ2∫10ϖ2tdx+χ∗∫10(3ψt−ϕt)2dx+ϵ2∫10(ϕ−ϖx)2dx+χ∗∫10r2dx+χ∗(1+1ϵ2)∫10θ2xdx,∀ϵ2>0. | (4.11) |
Proof. Simple calculations, using (2.8)1, (2.8)4 and integration by parts, we get
I′2(t)=−cϱ∫10ϖtt(∫1xθ(y)dy)dx−cϱ∫10ϖt(∫1xθt(y)dy)dx=cG∫10(ϕ−ϖx)θdx+k0ϱ∫10θxϖtdx+γc∫10θ2dx−γϱ∫10ϖ2tdx−k1ϱ∫10rϖtdx+mϱ∫10ϖt∫1x(3ψt−ϕt)(y)dydx. |
Now, thanks to Young, Poincaré's and Cauchy–Schwarz inequalities, we get, for any ϵ2>0,
I′2(t)≤−γϱ2∫10ϖ2tdx+χ∗∫10(3ψt−ϕt)2dx+ϵ2∫10(ϕ−ϖx)2dx+χ∗∫10r2dx+χ∗(1+1ϵ2)∫10θ2xdx. |
The proof is then completed.
Lemma 4.4. Consider the functional
I3(t):=ϱ∫10ϖtϖdx+ϱ∫10ϕ(∫x0ϖt(y)dy)dx, | (4.12) |
then, it satisfies
I′3(t)≤−G2∫10(ϕ−ϖx)2dx+ϱ∫10(3ψt−ϕt)2dx+3ϱ2∫10ϖ2tdx+χ∗∫10θ2xdx+9ϱ∫10ψ2tdx. | (4.13) |
Proof. We differentiate I3, using (2.8)1 together with integration by parts, to get
I′3(t)=ϱ∫10ϖ2tdx+ϱ∫10ϖttϖdx+ϱ∫10ϕt(∫x0ϖt(y)dy)dx+ϱ∫10ϕ(∫x0ϖtt(y)dy)dx=ϱ∫10ϖ2tdx−G∫10(ϕ−ϖx)xϖdx−γ∫10ϖθxdx+ϱ∫10ϕt(∫x0ϖt(y)dy)dx−G∫10(ϕ−ϖx)ϕdx−γ∫10θϕdx=ϱ∫10ϖ2tdx−G∫10(ϕ−ϖx)2dx−γ∫10(ϕ−ϖx)θdx+ϱ∫10ϕt(∫x0ϖt(y)dy)dx. |
Notice that
∫10ϕ2tdx≤2∫10(3ψt−ϕt)2dx+18∫10ψ2tdx. |
By Young, Poincaré's, and Cauchy-Schwarz inequalities, we easily prove (4.13).
Lemma 4.5. Consider the functional
I4(t):=Iϱ∫10(3ψ−ϕ)t(3ψ−ϕ)dx, | (4.14) |
then, it satisfies
I′4(t)≤−ˉl2∫10(3ψx−ϕx)2dx+Iϱ∫10(3ψt−ϕt)2dx+χ∗∫10(r2+θ2x)dx+χ∗∫10(ϕ−ϖx)2dx+χ∗(g⋄(3ψx−ϕx))(t). | (4.15) |
Proof. We proceed by differentiating the functional I4 and using Eq (2.8)2 together with integration by parts, which leads to
I′4(t)=Iϱ∫10(3ψ−ϕ)tt(3ψ−ϕ)dx+Iϱ∫10(3ψt−ϕt)2dx=Iϱ∫10(3ψt−ϕt)2dx−ˉl∫10(3ψx−ϕx)2dx+G∫10(3ψ−ϕ)(ϕ−ϖx)dx+m∫10(3ψ−ϕ)θdx+d∫10(3ψ−ϕ)xrdx−∫10(3ψ−ϕ)x∫∞0g(s)((3ψ−ϕ)x(t)−(3ψ−ϕ)x(t−s))dsdx. | (4.16) |
By virtue of Young's inequality and (2.4), we have
I′4(t)≤−ˉl2∫10(3ψx−ϕx)2dx+Iϱ∫10(3ψt−ϕt)2dx+χ∗∫10(r2+θ2x)dx+χ∗∫10(ϕ−ϖx)2dx+C1∫10[∫∞0g(s)((3ψ−ϕ)x(t)−(3ψ−ϕ)x(t−s))ds]2dx≤−ˉl2∫10(3ψx−ϕx)2dx+Iϱ∫10(3ψt−ϕt)2dx+χ∗∫10(r2+θ2x)dx+χ∗∫10(ϕ−ϖx)2dx+χ∗(g⋄(3ψx−ϕx))(t). |
This completes the proof of (4.15).
Lemma 4.6. Consider the functional
I5(t):=3Iϱ∫10ψtψdx+2β∫10ψ2dx, | (4.17) |
then, it satisfies the estimate
I′5(t)≤−2δ∫10ψ2dx−3D∫10ψ2xdx+3Iϱ∫10ψ2tdx+χ∗∫10(ϕ−ϖx)2dx+χ∗∫10∫ς2ς1|μ2(ς)|S2(x,1,p,t)dςdx. | (4.18) |
Proof. Simple computations using Eq (2.8)3 and integration by parts, yield
I′5(t)=3Iϱ∫10ψ2tdx−3D∫10ψ2xdx−4δ∫10ψ2dx−3G∫10(ϕ−ϖx)ψdx−4∫10∫ς2ς1ψ|μ2(ς)|S(x,1,p,t)dςdx. | (4.19) |
Employing Young's inequality, we conclude (4.18).
Lemma 4.7. Consider the functional
I6(t):=∫10∫10∫ς2ς1ςe−ςp|μ2(ς)|S2(x,p,ς,t)dςdpdx, | (4.20) |
then, it satisfies
I′6(t)≤β∫10ψ2tdx−m1∫10∫ς2ς1|μ2(ς)|S2(x,1,ς,t)dςdx−m1∫10∫10∫ς2ς1ς|μ2(ς)|S2(x,p,ς,t)dςdpdx, | (4.21) |
where m1 is a positive constant.
Proof. Taking the derivative of I6 and using (2.8)7 and S(x,0,t)=ψt, we have
I′6(t)=−2∫10∫10∫ς2ς1e−ςp|μ2(ς)|SpS(x,p,ς,t)dςdpdx=−∫10∫10∫ς2ς1ςe−ςp|μ2(ς)|S2(x,p,ς,t)dςdpdx−∫10∫ς2ς1|μ2(ς)|{e−ςS2(x,1,ς,t)−ψ2t(x,t)}dςdx. |
From e−ς≤e−ςp≤1, where 0<p<1, we arrive at
I′6(t)≤−∫10∫10∫ς2ς1ςe−ς|μ2(ς)|S2(x,p,ς,t)dςdpdx+(∫ς2ς1|μ2(ς)|dς)∫10ψ2t(x,t)dx−∫10∫ς2ς1e−ς|μ2(ς)|S2(x,1,ς,t)dςdx. |
Since −e−ς is an increasing function, then
−e−ς≤−e−ς2, forall ς∈[ς1,ς2]. |
Hence, if we denote m1=e−ς2 and use (T), we easily prove (4.21).
Let us now prove our stability result by using the lemmas in Section 4.
Proof of Theorem 2.2. We proceed by introducing a Lyapunov functional
L(t)=NE(t)+6∑j=1NjIj(t), | (5.1) |
where constants N,Nj>0,j=1,⋯,6, will be chosen later.
From (5.1), we write
|L(t)−NE(t)|≤IϱN1∫10|(3ψ−ϕ)t∫∞0g(s)((3ψ−ϕ)(t)−(3ψ−ϕ)(t−s))ds|dx+cϱN2∫10|ϖt∫1xθ(y)dy|dx+ϱN3∫10|ϖtϖ|dx+ϱN3∫10|ϕ∫x0ϖt(y)dy|dx+IϱN4∫10|(3ψ−ϕ)t(3ψ−ϕ)|dx+3IϱN5∫10|ψtψ|dx+2βN5∫10ψ2dx+N6∫10∫10∫ς2ς1ςe−ςp|μ2(ς)|S2(x,p,ς,t)dςdpdx. |
Thanks to Young, Cauchy-Schwarz and Poincaré's inequalities, we get
|L(t)−NE(t)|≤ϑ1E(t), where ϑ1>0, |
i.e.,
(N−ϑ1)E(t)≤L(t)≤(N+ϑ1)E(t). | (5.2) |
Now, differentiating the Lyapunov functional L(t), using (4.1), (4.7), (4.11), (4.13), (4.15), (4.18), and (4.21), and fixing
N4=N5=1,ϵ1=ˉl4N1,ϵ2=GN34N2. |
We find
ddtL(t)≤−(γϱ2N2−3ϱ2N3)∫10ϖ2tdx−(Iϱg02N1−χ∗N2−ϱN3−Iϱ)∫10(3ψt−ϕt)2dx−(m0N−9ϱN3−βN6−3Iϱ)∫10ψ2tdx−ˉl4∫10(3ψx−ϕx)2dx−(G4N3−(ˉl4+2χ∗))∫10(ϕ−ϖx)2dx−2δ∫10ψ2dx−(k0N−χ∗(1+N2N3)N2−χ∗N3−χ∗−ˉl4)∫10θ2xdx−3D∫10ψ2xdx−k2N∫10r2xdx−(k3N−χ∗N1−χ∗N2−χ∗)∫10r2dx−m1N6∫10∫10∫ς2ς1ς|μ2(ς)|S2(x,p,ς,t)dςdpdx−(m1N6−χ∗)∫10∫ς2ς1|μ2(ς)|S2(x,1,ς,t)dςdx+(N2−χ∗N1)(g′⋄(3ψx−ϕx))(t)+(χ∗(1+4N1ˉl)N1+χ∗)(g⋄(3ψx−ϕx))(t). | (5.3) |
Next, we choose our coefficients in (5.3), in a way that, they all except the last two become negative. We start by selecting N6 big enough so that
m1N6−χ∗>0, |
then, we take N3 fairly wide, such that
G4N3−(ˉl4+2χ∗)>0, |
after that, we choose N2 large enough, so that
γϱ2N2−3ϱ2N3>0, |
now, we select N1 sufficiently large such that
Iϱg02N1−χ∗N2−ϱN3−Iϱ>0. |
We can now select N large enough so that we have (5.2) and
{12N−χ∗N1>0,m0N−9ϱN3−βN6−3Iϱ>0,k3N−χ∗N1−χ∗N2−χ∗>0,k0N−χ∗(1+N2N3)N2−χ∗N3−χ∗−ˉl4>0. |
Hence, relation (5.3) becomes
ddtL(t)≤−ϑ2∫10{ϖ2t+(ϕ−ϖx)2+(3ψt−ϕt)2+ψ2+(3ψx−ϕx)2+ψ2x+ψ2t+θ2+r2}dx−ϑ2∫10∫10∫ς2ς1ςe−ςp|μ2(ς)|S2(x,p,ς,t)dςdpdx+ϑ3(g⋄(3ψx−ϕx))(t),ϑ2,ϑ3>0. | (5.4) |
Now, exploiting (2.13) and Poincaré's inequality, we obtain
E(t)≤ϑ4∫10{ϖ2t+(ϕ−ϖx)2+(3ψt−ϕt)2+ψ2+(3ψx−ϕx)2+ψ2x+ψ2t+θ2+r2}dx+ϑ4∫10∫10∫ς2ς1ςe−ςp|μ2(ς)|S2(x,p,ς,t)dςdpdx+ϑ4(g⋄(3ψx−ϕx))(t), where ϑ4>0, |
from which
−∫10{ϖ2t+(ϕ−ϖx)2+(3ψt−ϕt)2+ψ2+(3ψx−ϕx)2+ψ2x+ψ2t+θ2+r2}dx−∫10∫10∫ς2ς1ςe−ςp|μ2(ς)|S2(x,p,ς,t)dςdpdx−(g⋄(3ψx−ϕx))(t)≤−ϑ5E(t), | (5.5) |
where ϑ5>0. Thereby, if we combine (5.5) and (5.4), we have
ddtL(t)≤−ϑ6E(t)+ϑ7(g⋄(3ψx−ϕx))(t), where ϑ6,ϑ7>0. | (5.6) |
Next, we multiply (5.6), by
G′(ϵ0E(t)E(0)), |
we find
G′(ϵ0E(t)E(0))ddtL(t)≤−ϑ6G′(ϵ0E(t)E(0))E(t)+ϑ7G′(ϵ0E(t)E(0))(g⋄(3ψx−ϕx))(t). | (5.7) |
Now, we estimate the last term in (5.7) and use both (A2) and (2.14), we find
ϑ7G′(ϵ0E(t)E(0))(g⋄(3ψx−ϕx))(t)≤−ϑ′7E′(t)+ϑ′7ϵ0G0(E(t)E(0)),ϑ′7>0. | (5.8) |
We insert (5.8) in (5.7) and set ϵ0=ϑ6E(0)2ϑ′7, we get
G′(ϵ0E(t)E(0))ddtL(t)+ϑ′7E′(t)≤−ΓG0(E(t)E(0)),Γ>0. | (5.9) |
We consider now the functional
{L}_1(t): = {G}^\prime\left( \frac{\epsilon_0 {E}(t)}{ {E}(0)} \right) {L}(t)+\vartheta_7' {E}(t). |
It is clear that
{L}_1(t) \sim {E}(t), |
moreover, noticing that {E}^\prime(t)\leq0, \; {G}^{\prime \prime}(t) > 0, we obtain
\begin{equation} \frac{d}{dt} {L}_1(t)\leq -\Gamma {G}_0\left( \frac{ {E}(t)}{ {E}(0)} \right). \end{equation} | (5.10) |
Next, we present the functional
{L}_2(t): = b_1\frac{ {L}_1(t)}{ {E}(0)}\sim {E}(t), \quad \text{ such that } |
\begin{cases} {L}_2(t)\leq 1, \\ \frac{d}{dt} {L}_2(t) \leq -\alpha_2 {G}_0( {L}_2(t)), \end{cases} |
where, \alpha_2 is a positive constant, therefore,
{G}^\prime_*( {L}_2(t))\geq \alpha_2. |
We integrate over (0, t) to find
{L}_2(t)\leq {G}^{-1}_*(\alpha_2 t+\alpha_3), |
from which, we deduce that
{E}(t)\leq \alpha_1 {G}^{-1}_*(\alpha_2 t+\alpha_3), |
where, \alpha_1 and \alpha_3 are positive constants. The proof is then completed.
The article is about the laminated beam system along with structural damping, past history, distributed delay, and in the presence of both temperatures and micro-temperature effects introduced in (1.1). By the semigroup approach, we established the existence and uniqueness of the solution which can be considered as the first main result. In addition, as a second novelty, a general decay result for the solution unusually with no constraints regarding the speeds of wave propagation is found. This last new result is considered, as far as we know, the first similar result in the literature for such a system, where we succeed to improve the earlier works known for the case of finite history, to the case of infinite history. The relaxation function becomes intended to satisfy a broader class of relaxation functions.
We mention here that the distributed delay in our system makes a good interaction between the past history and the other damping terms of system (1.1). This type of damping gives more information and qualitative properties on the solution and also its impact on stability is very important as it is shown in the requirement of Theorem 2.2. Of course, the other terms (both temperatures and micro-temperature effects) act as balances in the stability of the system.
Fares Yazid and Fatima Siham Djerad: Writing—original draft preparation; Abdelkader Moumen and Moheddine Imsatfia, Tayeb Mahrouz: Writing—review and editing; Keltoum Bouhali: Supervision. All authors have read and approved the final version of the manuscript for publication.
The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.
The authors extend their appreciation to the Deanship of Research and Graduate Studies at King Khalid University for funding this work through small group research project under grant number RGP1/21/45.
The authors declare that there is no conflict of interest.
[1] |
M. S. Alves, R. N. Monteiro, Exponential stability of laminated Timoshenko beams with boundary/internal controls, J. Math. Anal. Appl., 482 (2020), 123516. https://doi.org/10.1016/j.jmaa.2019.123516 doi: 10.1016/j.jmaa.2019.123516
![]() |
[2] |
T. A. Apalara, On the stability of a thermoelastic laminated beam, Acta Math. Sci., 39 (2019), 1517–1524. https://doi.org/10.1007/s10473-019-0604-9 doi: 10.1007/s10473-019-0604-9
![]() |
[3] | L. Bouzettouta, S. Zitouni, K. Zennir, A. Guesmia, Stability of Bresse system with internal distributed delay, J. Math. Comput. Sci., 17 (2017), 92–118. |
[4] | A. Choucha, D. Ouchenane, K. Zennir, B. Feng, Global well-posedness and exponential stability results of a class of Bresse-Timoshenko-type systems with distributed delay term, Math. Meth. Appl. Sci., 2020 (2020). https://doi.org/10.1002/mma.6437. |
[5] |
N. Doudi, S. Boulaaras, Global existence combined with general decay of solutions for coupled Kirchhoff system with a distributed delay term, RACSAM Rev. R. Acad. A, 114 (2020), 1–31. https://doi.org/10.1007/s13398-020-00938-9 doi: 10.1007/s13398-020-00938-9
![]() |
[6] |
H. Dridi, B. Feng, K. Zennir, Stability of Timoshenko system coupled with thermal law of Gurtin-Pipkin affecting on shear force, Appl. Anal., 101 (2022), 1–15. https://doi.org/10.1080/00036811.2021.1883591 doi: 10.1080/00036811.2021.1883591
![]() |
[7] |
H. Dridi, K. Zennir, Well-posedness and energy decay for some thermoelastic systems of Timoshenko type with Kelvin-Voigt damping, SeMA J., 78 (2021), 385–400. https://doi.org/10.1007/s40324-021-00239-0 doi: 10.1007/s40324-021-00239-0
![]() |
[8] |
A. S. El-Karamany, M. A. Ezzat, On the phase-lag Green-Naghdi thermoelasticity theories, Appl. Math. Model., 40 (2016), 5643–5659. https://doi.org/10.1016/j.apm.2016.01.010 doi: 10.1016/j.apm.2016.01.010
![]() |
[9] |
D. Fayssal, Well posedness and stability result for a thermoelastic laminated beam with structural damping, Ric. Mat., 2022 (2022), 1–25. https://doi.org/10.1007/s11587-022-00708-2 doi: 10.1007/s11587-022-00708-2
![]() |
[10] | E. I. Grigolyuk, Nonlinear behavior of shallow rods, Dokl. Akad. Nauk, 348 (1996), 759–763. |
[11] | E. I. Grigolyuk, E. A. Lopanitsyn, Large axisymmetric flexures of thin short shells of revolution under small deformations, Dokl. Akad. Nauk, 346 (1996), 753–756. |
[12] |
S. W. Hansen, R. D. Spies, Structural damping in laminated beams due to interfacial slip, J. Sound Vib., 204 (1997), 183–202. https://doi.org/10.1006/jsvi.1996.0913 doi: 10.1006/jsvi.1996.0913
![]() |
[13] |
H. E. Khochemane, General stability result for a porous thermoelastic system with infinite history and microtemperatures effects, Math. Meth. Appl. Sci., 45 (2022), 1538–1557. https://doi.org/10.1002/mma.7872 doi: 10.1002/mma.7872
![]() |
[14] |
W. Liu, W. Zhao, Stabilization of a thermoelastic laminated beam with past history, Appl. Math. Optim., 80 (2019), 103–133. https://doi.org/10.1007/s00245-017-9460-y doi: 10.1007/s00245-017-9460-y
![]() |
[15] | Z. Liu, S. Zheng, Semi-groups associated with dissipative systems, CRC Press, 1999. |
[16] |
A. Moumen, D. Ouchenane, A. Choucha, K. Zennir, S. A. Zubair, Exponential stability of Timoshenko system in thermoelasticity of second sound with a memory and distributed delay term, Open Math., 19 (2022), 1636–1647. https://doi.org/10.1515/math-2021-0117 doi: 10.1515/math-2021-0117
![]() |
[17] | V. F. Nesterenko, Dynamics of heterogeneous materials, Springer, New York, 2001. |
[18] |
V. F. Nesterenko, Propagation of nonlinear compression pulses ingranular media, J. Appl. Mech. Tech. Phys., 24 (1984), 733–743. https://doi.org/10.1007/BF00905892 doi: 10.1007/BF00905892
![]() |
[19] |
S. Nicaise, C. Pignotti, Stabilization of the wave equation with boundary or internal distributed delay, Differ. Integral. Equ., 21 (2008), 935–958. https://doi.org/10.57262/die/1356038593 doi: 10.57262/die/1356038593
![]() |
[20] | C. Nonato, C. Raposo, B. Feng, Exponential stability for a thermoelastic laminated beam with nonlinear weights and time-varying delay, Asymptotic Anal., 126 (2022), 157–185. |
[21] |
E. Pişkin, J. Ferreira, H. Yuksekkaya, M. Shahrouzi, Existence and asymptotic behavior for a logarithmic viscoelastic plate equation with distributed delay, Int. J. Nonlinear. Anal., 13 (2022), 763–788. https://doi.org/10.22075/IJNAA.2022.24639.2797 doi: 10.22075/IJNAA.2022.24639.2797
![]() |
[22] |
C. A. Raposo, O. V. Villagran, J. E. Muñoz Rivera, M. S. Alves, Hybrid laminated Timoshenko beam, J. Math. Phys., 58 (2017), 101512. https://doi.org/10.1063/1.4998945 doi: 10.1063/1.4998945
![]() |
[23] |
S. P. Timoshenko, On the correction for shear of the differential equation for transverse vibrations of prismatic bars, Lond. Edinb. Dublin. Philos. Mag. J. Sci., 41 (1921), 744–746. https://doi.org/10.1080/14786442108636264 doi: 10.1080/14786442108636264
![]() |
[24] |
X. Tian, O. Zhang, Stability of a Timoshenko system with local Kelvin-Voigt damping, Z. Angew. Math. Phys., 68 (2017), 20. https://doi.org/10.1007/s00033-016-0765-5 doi: 10.1007/s00033-016-0765-5
![]() |