Research article Special Issues

Biogenic synthesis, characterization and effects of Mn-CuO composite nanocatalysts on Methylene blue photodegradation and Human erythrocytes

  • Each year more than 150, 000 tons of dyes are released in effluents by industries. These chemicals entities non-biodegradable and toxic can be removed from effluent by metallic nanomaterials. The aqueous extract of Manotes expansa leaves is used as reducing and stabilizing agent in the biogenic synthesis of Mn-CuO nanocomposites. The nanoparticles obtained were characterized using UV-visible spectroscopy, X-ray Diffraction (XRD), X-ray Fluorescence, Dynamic Light Scattering (DSL), and Scanning Electron Microscopy (SEM). The hemotoxicity of biosynthesized nanomaterials was assessed by evaluating their hemolytic activity using erythrocytes as a model system. The photocatalytic activity of Mn-CuO was carried out by photocatalytic degradation of Methylene Blue dye as a model. The results obtained by UV-vis spectroscopy showed a Plasmonic Surface Resonance band at 408 nm. XRD and X-ray fluorescence made it possible to identify the presence of particles of formula Mn0.53Cu0.21O having crystallized in a Hexagonal system (a = 3.1080 Å and c = 5.2020 Å). Spherical morphology and average height 49.34 ± 6.71 nm were determined by SEM and DSL, respectively. The hemolytic activity of biosynthesized nanomaterials revealed that they are not hemotoxic in vitro (% hemolysis 3.2%) and 98.3% of Methylene Blue dye was removed after 120 min under irradiation with solar light in the presence of Mn-CuO nanocomposites.

    Citation: Carlos N. Kabengele, Giresse N. Kasiama, Etienne M. Ngoyi, Clement L. Inkoto, Juvenal M. Bete, Philippe B. Babady, Damien S. T. Tshibangu, Dorothée D. Tshilanda, Hercule M. Kalele, Pius T. Mpiana, Koto-Te-Nyiwa Ngbolua. Biogenic synthesis, characterization and effects of Mn-CuO composite nanocatalysts on Methylene blue photodegradation and Human erythrocytes[J]. AIMS Materials Science, 2023, 10(2): 356-369. doi: 10.3934/matersci.2023019

    Related Papers:

    [1] Hongying Jiao, Shuhai Zhu, Jinguo Zhang . Existence of infinitely many solutions for critical sub-elliptic systems via genus theory. Communications in Analysis and Mechanics, 2024, 16(2): 237-261. doi: 10.3934/cam.2024011
    [2] Jizheng Huang, Shuangshuang Ying . Hardy-Sobolev spaces of higher order associated to Hermite operator. Communications in Analysis and Mechanics, 2024, 16(4): 858-871. doi: 10.3934/cam.2024037
    [3] Xiulan Wu, Yaxin Zhao, Xiaoxin Yang . On a singular parabolic p-Laplacian equation with logarithmic nonlinearity. Communications in Analysis and Mechanics, 2024, 16(3): 528-553. doi: 10.3934/cam.2024025
    [4] Erlend Grong, Irina Markina . Harmonic maps into sub-Riemannian Lie groups. Communications in Analysis and Mechanics, 2023, 15(3): 515-532. doi: 10.3934/cam.2023025
    [5] Leandro Tavares . Solutions for a class of problems driven by an anisotropic (p,q)-Laplacian type operator. Communications in Analysis and Mechanics, 2023, 15(3): 533-550. doi: 10.3934/cam.2023026
    [6] Velimir Jurdjevic . Time optimal problems on Lie groups and applications to quantum control. Communications in Analysis and Mechanics, 2024, 16(2): 345-387. doi: 10.3934/cam.2024017
    [7] Zhiyong Wang, Kai Zhao, Pengtao Li, Yu Liu . Boundedness of square functions related with fractional Schrödinger semigroups on stratified Lie groups. Communications in Analysis and Mechanics, 2023, 15(3): 410-435. doi: 10.3934/cam.2023020
    [8] Ming Liu, Binhua Feng . Grand weighted variable Herz-Morrey spaces estimate for some operators. Communications in Analysis and Mechanics, 2025, 17(1): 290-316. doi: 10.3934/cam.2025012
    [9] Xiaotian Hao, Lingzhong Zeng . Eigenvalues of the bi-Xin-Laplacian on complete Riemannian manifolds. Communications in Analysis and Mechanics, 2023, 15(2): 162-176. doi: 10.3934/cam.2023009
    [10] Eleonora Amoroso, Angela Sciammetta, Patrick Winkert . Anisotropic (p,q)-Laplacian problems with superlinear nonlinearities. Communications in Analysis and Mechanics, 2024, 16(1): 1-23. doi: 10.3934/cam.2024001
  • Each year more than 150, 000 tons of dyes are released in effluents by industries. These chemicals entities non-biodegradable and toxic can be removed from effluent by metallic nanomaterials. The aqueous extract of Manotes expansa leaves is used as reducing and stabilizing agent in the biogenic synthesis of Mn-CuO nanocomposites. The nanoparticles obtained were characterized using UV-visible spectroscopy, X-ray Diffraction (XRD), X-ray Fluorescence, Dynamic Light Scattering (DSL), and Scanning Electron Microscopy (SEM). The hemotoxicity of biosynthesized nanomaterials was assessed by evaluating their hemolytic activity using erythrocytes as a model system. The photocatalytic activity of Mn-CuO was carried out by photocatalytic degradation of Methylene Blue dye as a model. The results obtained by UV-vis spectroscopy showed a Plasmonic Surface Resonance band at 408 nm. XRD and X-ray fluorescence made it possible to identify the presence of particles of formula Mn0.53Cu0.21O having crystallized in a Hexagonal system (a = 3.1080 Å and c = 5.2020 Å). Spherical morphology and average height 49.34 ± 6.71 nm were determined by SEM and DSL, respectively. The hemolytic activity of biosynthesized nanomaterials revealed that they are not hemotoxic in vitro (% hemolysis 3.2%) and 98.3% of Methylene Blue dye was removed after 120 min under irradiation with solar light in the presence of Mn-CuO nanocomposites.



    The fundamental work of Hansen and Spies [4] modeled a two-layer beam with a structural damping due to the interfacial slip through the following system

    {ρφtt+G(ψφx)x=0,Iρ(3wψ)ttD(3wψ)xxG(ψφx)=0,IρwttDwxx+3G(ψφx)+4γw+4βwt=0, (1.1)

    where φ=φ(x,t) is the transverse displacement, ψ=ψ(x,t) is the rotation angle, w=w(x,t) is proportional to the amount of slip along the interface, 3wψ denotes the effective rotation angle. The physical quantities ρ,Iρ,G,D,β and γ are respectively: the density, mass moment of inertia, shear stiffness, flexural rigidity, adhesive damping and adhesive stiffness. Equation (1.1)3 describes the dynamics of the slip. For β=0, system (1.1) describes the coupled laminated beams without structural damping at the interface. In the recent result [1], Apalara considered the thermoelastic-laminated beam system without structural damping, namely

    {ρφtt+G(ψφx)x=0,Iρ(3sψ)ttD(3sψ)xxG(ψφx)=0,IρsttDsxx+3G(ψφx)+4γs+δθx=0,ρ3θtλθxx+δstx=0, (1.2)

    where (x,t)(0,1)×(0,+), θ=θ(x,t) is the difference temperature. The positive quantities γ,β,k,λ are adhesive stiffness, adhesive damping, heat capacity and the diffusivity respectively. The author proved that (1.2) is exponential stable provided

    Gρ=DIρ. (1.3)

    When β>0, the adhesion at the interface supplies a restoring force proportion to the interfacial slip. But this is not enough to stabilize system (1.1), see for instance [2]. To achieve exponential or general stabilization of system (1.1), many authors in literature have used additional damping. In this direction, Gang et al. [9] studied the following memory-type laminated beam system

    {ρφtt+G(ψφx)x=0,Iρ(3wψ)ttD(3wψ)xx+t0g(ts)(3wψ)xx(x,s)dsG(ψφx)=03Iρwtt3Dwxx+3G(ψφx)+4γw+4βwt=0 (1.4)

    and established a general decay result for more regular solutions and GρDIρ. Mustafa [15] also considered the structural damped laminated beam system (1.4) and established a general decay result provided Gρ=DIρ. Feng et al. [8] investigated the following laminated beam system

    {ρwtt+Gφx+g1(wt)+f1(w,ξ,s)=h1,IρξttGφDξxx+g2(ξt)+f2(w,ξ,s)=h2,Iρstt+GφDsxx+g3(st)+f2(w,ξ,s)=h3 (1.5)

    and established the well-posedness, smooth global attractor of finite fractal dimension as well as existence of generalized exponential attractors. See also, recent results by Enyi et al. [20]. We refer the reader to [5,6,7,11,13,14,17,18] and the references cited therein for more related results.

    In this present paper, we consider a thermoelastic laminated beam problem with a viscoelastic damping

    {ρwtt+G(ψwx)x=0,Iρ(3sψ)ttD(3sψ)xx+t0g(tτ)(3sψ)xx(x,τ)dτG(ψwx)=03Iρstt3Dsxx+3G(ψwx)+4γs+δθx=0,kθtλθxx+δsxt=0 (1.6)

    under initial conditions

    {w(x,0)=w0(x), ψ(x,0)=ψ0(x), s(x,0)=s0(x), θ(x,0)=θ0(x),  x[0,1],wt(x,0)=w1(x), ψt(x,0)=ψ1(x), st(x,0)=s1(x),                          x[0,1] (1.7)

    and boundary conditions

    {w(0,t)=ψx(0,t)=sx(0,t)=θ(0,t)=0,t[0,+),wx(1,t)=ψ(1,t)=s(1,t)=θx(1,t)=0,t[0,+). (1.8)

    In the system (1.6), the integral represents the viscoelastic damping, and g is the relaxation function satisfying some suitable assumptions specified in the next section. According to the Boltzmann Principle, the viscoelastic damping (see [21] for details) is represented by a memory term in the form of convolution. It acts as a damper to reduce the internal/external forces like the beam's weight, heavy loads, wind, etc., that cause undesirable vibrations.

    In most of the above works, the authors have established their decay result by including the structural damping along with other dampings. So, the natural question that comes to mind.

    Is it possible to obtain general/optimal decay result (decay rates that agrees with that of g) to the thermoelastic laminated beam system (1.6)–(1.8), in the absence of the structural damping.

    The novelty of this article is to answer this question in a consenting way, by using the ideas developed in [10] to establish general and optimal decay results for Problem 1.6. Moreover, we establish a weaker decay result in the case of a non-equal wave of speed propagation. To the best of our knowledge, there is no stability result for the latter in the literature.

    The rest of work is organized as follows: In Section 2, we recall some preliminaries and assumptions on the memory term. In Section 3, we state and prove the main stability result for the case equal-speed and in the case of non-equal-speed of propagation. We also give some examples to illustrate our findings. Finally, in Section 4, we give the proofs of the lemmas used our main results.

    In this section, we recall some useful materials and conditions. Through out this paper, C is a positive constant that may change through lines, .,. and .2 denote respectively the inner product and the norm in L2(0,1). We assume the relaxation function g obeys the assumptions:

    (G1). g:[0,+)(0,+) is a non-inecreasing C1 function such that

    g(0)>0,D0g(τ)dτ=l0>0. (2.1)

    (G2). There exist a C1 function H:[0,+)(0,+) which is linear or is strictly convex C2 function on (0,ϵ0), ϵ0g(0), with H(0)=H(0)=0 and a positive nonincreasing differentiable function ξ:[0,+)(0,+), such that

    g(t)ξ(t)H(g(t)),t0, (2.2)

    Remark 2.1. As in [10], we note here that, if H is a strictly increasing convex C2 function on (0,r], with H(0)=H(0)=0, then H has an extension ˉH, which is strictly increasing and strictly convex C2-function on (0,+). For example, ˉH can be defined by

    ˉH(s)=H(r)2s2+(H(r)H(r)r)s+H(r)H(r)r+H(r)2r2, s>r. (2.3)

    Let

    H1(0,1)={uH1(0,1)/u(0)=0},   ˉH1(0,1)={uH1(0,1)/u(1)=0},
    H2(0,1)={uH2(0,1)/uxH1(0,1)},   ˉH2(0,1)={uH2(0,1)/uxˉH1(0,1)}.

    The existence and regularity result of problem (1.6) is the following

    Theorem 2.1. Let (w0,3s0ψ0,s0,θ0)H1(0,1)×ˉH1(0,1)×ˉH1(0,1)×H1(0,1) and (w1,3s1ψ1,s1)L2(0,1)×L2(0,1)×L2(0,1) be given. Suppose (G1) and (G2) hold. Then problem (1.6) has a unique global weak solution (w,3sψ,s,θ) which satisfies

    wC(R+,H1(0,1))C1(R+,L2(0,1)), (3sψ)C(R+,ˉH1(0,1))C1(R+,L2(0,1)),
    sC(R+,ˉH1(0,1))C1(R+,L2(0,1)), θC(R+,L2(0,1))L2(R+,H1(0,1)).

    Furthermore, if (w0,(3s0ψ0),s0,θ0)H2(0,1)×ˉH2(0,1)×ˉH2(0,1)×H2(0,1)H1(0,1) and (w1,(3s1ψ1),s1)H1(0,1)×ˉH1(0,1)×ˉH1(0,1), then the solution of (1.6) satisfies

    wC(R+,H2(0,1))C1(R+,H1(0,1))C2(R+,L2(0,1)),
    (3sψ)C(R+,ˉH2(0,1))C1(R+,ˉH1(0,1))C2(R+,L2(0,1)),
    sC(R+,ˉH2(0,1))C1(R+,ˉH1(0,1))C2(R+,L2(0,1)),
    θC(R+,H2(0,1)H2(0,1))C1(R+,H1(0,1)).

    The proof of Theorem 2.1 can be established using the Galerkin approximation method as in [16]. Throughout this paper, we denote by the binary operator, defined by

    (gν)(t)=t0g(tτ)ν(t)ν(τ)22dτ,t0.

    We also define h(t) and Cα as follow

    h(t)=αg(t)g(t)   and   Cα=+0g2(τ)αg(τ)g(τ)dτ.

    The following lemmas will be applied repeatedly throughout this paper

    Lemma 2.1. For any function fL2loc([0,+),L2(0,1)), we have

    10(t0g(ts)(f(t)f(s))ds)2dx(1l0)(gf)(t), (2.4)
    10(x0f(y,t)dy)2dxf(t)22. (2.5)

    Lemma 2.2. Let vH1(0,1)  or  ˉH1(0,1), we have

    10(t0g(ts)(v(t)v(τ))dτ)2dxCp(1l0)(gv)(t), (2.6)

    where Cp>0 is the poincaré constant.

    Lemma 2.3. Let (w,3sψ,s,θ) be the solution of (1.6). Then, for any 0<α<1 we have

    10(t0g(tτ)((3sψ)x(τ)(3sψ)x(t))dτ)2dxCα(h(3sψ)x)(t). (2.7)

    Proof. Using Cauchy-Schwarz inequality, we have

    10(t0g(tτ)((3sψ)x(τ)(3sψ)x(t))dτ)2dx=10(t0g(tτ)h(tτ)h(tτ)((3sψ)x(τ)(3sψ)x(t))dτ)2dx(+0g2(τ)h(τ)ds)10t0h(tτ)((3sψ)x(τ)(3sψ)x(t))2dτdx=Cα(h(3sψ)x)(t). (2.8)

    Lemma 2.4. [12] Let F be a convex function on the close interval [a,b], f,j:Ω[a,b] be integrable functions on Ω, such that j(x)0 and Ωj(x)dx=α1>0. Then, we have the following Jensen inequality

    F(1α1Ωf(y)j(y)dy)1α1ΩF(f(y))j(y)dy. (2.9)

    In particular if F(y)=y1p,  y0,  p>1, then

    (1α1Ωf(y)j(y)dy)1p1α1Ω(f(y))1pj(y)dy. (2.10)

    Lemma 2.5. The energy functional E(t) of the system (1.6)-(1.8) defined by

    E(t)=12[ρwt22+3Iρst22+Iρ3stψt22+3Dsx22+Gψwx22]+12[(Dt0g(τ)dτ)3sxψx22+(g(3sxψx))(t)+4γs22+kθ22], (2.11)

    satisfies

    E(t)=12(g(3sxψx))(t)12g(t)3sxψx22λθx2212(g(3sxψx))(t)0,   t0. (2.12)

    Proof. Multiplying (1.6)1, (1.6)2, (1.6)3 and (1.6)4, respectively, by wt, (3stψt), st and θ, integrating over (0,1), and using integration by parts and the boundary conditions (1.7), we arrive at

    12ddt(ρwt22+Gψwx22)=G(ψwx),ψt, (2.13)
    12ddt[Iρ3stψt22+(Dt0g(τ)dτ)3sxψx22+(g(3sxψx))(t)]=G(ψwx),(3sψ)t+12(g(3sxψx))(t)12g(t)3sxψx22, (2.14)
    12ddt[3Iρst22+3Dsx22+4γs22]=3G(ψwx),stδθx,st, (2.15)

    and

    12ddt(kθ22)=λθx22+δθx,st. (2.16)

    Adding the equations (2.13)–(2.16), taking into account (G1) and (G2), we obtain (2.12) for regular solutions. The result remains valid for weak solutions by a density argument. This implies the energy functional is non-increasing and

    E(t)E(0),  t0.

    This section is subdivided into two. In the first subsection, we prove the stability result for equal-wave-speed of propagation, whereas in the second subsection, we focus on the stability result for non-equal-wave-speed of propagation.

    Our aim, in this subsection, is to prove an explicit, general and optimal decay rate of solutions for system (1.6)–(1.8). To achieve this, we define a Lyapunov functional

    L(t)=NE(t)+6j=1NjIj(t), (3.1)

    where N, Nj, j=1,2,3,4,5,6 are positive constants to be specified later and

    I1(t)=Iρ10(3sψ)tt0g(tτ)((3sψ)(t)(3sψ)(τ))dτdx,t0,
    I2(t)=3Iρ10sstdx+3ρ10wtx0s(y)dydx,I3(t)=3kIρ10θx0st(y)dydx,t0,
    I4(t)=ρ10wtwdx,I5(t)=Iρ10(3sψ)(3sψ)tdx,t0,
    I6(t)=3IρG10(ψwx)stdx3ρD10wtsxdx,I7(t)=10t0J(tτ)(3sxψx)2(τ)dτdx,t0,

    where

    J(t)=+tg(τ)dτ.

    The following lemma is very important in the proof of our stability result.

    Lemma 3.1. Suppose Gρ=DIρ. Under suitable choice of t0,N, Nj, j=1,2,3,4,5,6, the Lyapunov functional L satisfies, along the solution of (1.6)(1.8), the estimate

    L(t)β(wt22+st22+3stψt22+sx22+3wxψx22+ψwx22)β(s22+θx22)+12(g(3sxψx))(t), tt0 (3.2)

    and the equivalence relation

    α1E(t)L(t)α2E(t) (3.3)

    holds for some β>0, α1, α2>0.

    Proof. By virtue of assumption (3.1) and using h(t)=αg(t)g(t), it follows from Lemmas 2.5, 4.1-4.6 (see the Appendix for detailed derivations) that, for all tt0>0,

    L(t)[N4ρN2δ4]wt22[N3δIρ2N2C(1+1ϵ2)N6C(1+1ϵ1)]st223N2γs22[N1Iρg0N5IρN6ϵ1]3stψt22[3DN2N3ϵ3N4CN6C]sx22[N6G2N1ϵ2N3ϵ3N4Cϵ4N5C]ψwx22[N5l04N1ϵ1N4ϵ4]3sxψx22[λNN2CN3C(1+1ϵ3)N6C]θx22+Nα2(g(3sxψx))(t)[N2CCα(N5+N1(1+1ϵ1+1ϵ2))](h(3sxψx))(t). (3.4)

    Now, we choose

    N4=N5=1,  ϵ4=l08 (3.5)

    and select N1 large enough such that

    μ1:=N1Iρg0Iρ>0. (3.6)

    Next, we choose N6 large so that

    μ2:=N6G2C>0. (3.7)

    Also, we select N2 large enough so that

    μ3:=3DN2CN6C>0. (3.8)

    After fixing N1,N2,N6, and letting ϵ3=μ12N3, we then select ϵ1,ϵ2, and δ4 very small such that

    ρN2δ4>0,  μ1N6ϵ1>0,   μ4:=μ22N1ϵ2>0 (3.9)

    and select N3 large enough so that

    N3δIρ2N2C(1+1ϵ2)N6C(1+1ϵ1)>0. (3.10)

    Now, we note that αg2(s)h(s)=αg2(s)αg(s)g(s)<g(s); thus the dominated convergence theorem gives

    αCα=+0αg2(s)αg(s)g(s)ds0  as  α0. (3.11)

    Therefore, we can choose some 0<α0<1 such that for all 0<αα0,

    αCα<14C(1+N1(1+1ϵ1+1ϵ2)). (3.12)

    Finally, we select N so large enough and take α=1N So that

    λNN2CN3C(1+1ϵ3)N6C>0,N2CCα(1+N1(1+1ϵ1+1ϵ2))>0. (3.13)

    Combination of (3.6) - (3.13) yields the estimate (3.2). The equivalent relation (3.3) can be obtain easily by using Young's, Cauchy-Schwarz, and Poincaré's inequalities.

    Now, we state and prove our stability result for this subsection.

    Theorem 3.1. Assume Gρ=DIρ and (G1) and (G2) hold. Then, there exist positive constants a1 and a2 such that the energy solution (2.11) satisfies

    E(t)a2H11(a1tt0ξ(τ)dτ),   where  H1(t)=rt1τH(τ)dτ (3.14)

    and H1 is a strictly decreasing and strictly convex function on (0,r], with limt0H1(t)=+.

    Proof. Using the fact that g and ξ are positive, non-increasing and continuous, and H is positive and continuous, we have that for all t[0,t0]

    0<g(t0)g(t)g(0),  0<ξ(t0)ξ(t)ξ(0).

    Thus for some constants a,b>0, we obtain

    aξ(t)H(g(t))b.

    Therefore, for any t[0,t0], we get

    g(t)ξ(t)H(g(t))ag(0)g(0)ag(0)g(t) (3.15)

    and

    ξ(t)g(t)g(0)ag(t). (3.16)

    From (2.12) and (3.15), it follows that

    t00g(τ)(3sxψx)(t)(3sxψx)(tτ)22dτg(0)at00g(τ)(3sxψx)(t)(3sxψx)(tτ)22dτCE(t),  tt0. (3.17)

    From (3.2) and (3.17), we have

    L(t)βE(t)+12(g(3sxψx))(t)=βE(t)+12t00g(τ)(3sxψx)(t)(3sxψx)(tτ)22dτ+12tt0g(τ)(3sxψx)(t)(3sxψx)(tτ)22dτβE(t)CE(t)+12tt0g(τ)(3sxψx)(t)(3sxψx)(tτ)22dτ.

    Thus, we get

    L1(t)βE(t)+12tt0g(τ)(3sxψx)(t)(3sxψx)(tτ)22dτ, tt0, (3.18)

    where L1=L+CEE by virtue of (3.3). To finish our proof, we distinct two cases:

    Case 1: H(t) is linear. In this case, we multiply (3.18) by ξ(t), keeping in mind (2.12) and (G2), to get

    ξ(t)L1(t)βξ(t)E(t)+12ξ(t)tt0g(τ)(3sxψx)(t)(3sxψx)(tτ)22dτβξ(t)E(t)+12tt0ξ(τ)g(τ)(3sxψx)(t)(3sxψx)(tτ)22dτβξ(t)E(t)12tt0g(τ)(3sxψx)(t)(3sxψx)(tτ)22dτβξ(t)E(t)CE(t),  tt0. (3.19)

    Therefore

    (ξL1+CE)(t)βξ(t)E(t),    tt0. (3.20)

    Since ξ is non-increasing and L1E, we have

    L2=ξL1+CEE. (3.21)

    Thus, from (3.20), we get for some positive constant α

    L2(t)βξ(t)E(t)αξ(t)L2(t),   tt0. (3.22)

    Integrating (3.22) over (t0,t) and recalling (3.21), we obtain

    E(t)a1ea2tt0ξ(s)ds=a1H11(a2tt0ξ(s)ds).

    Case 2: H(t) is nonlinear. In this case, we consider the functional L(t)=L(t)+I7(t). From (3.2) and Lemma 4.7 (see the Appendix), we obtain

    L(t)dE(t),  tt0, (3.23)

    where d>0 is a positive constant. Therefore,

    dtt0E(s)dsL(t0)L(t)L(t0).

    Hence, we get

    +0E(s)ds<. (3.24)

    Using (3.24), we define p(t) by

    p(t):=ηtt0(3sxψx)(t)(3sxψx)(tτ)22dτ,

    where 0<η<1 so that

    p(t)<1,tt0. (3.25)

    Moreover, we can assume p(t)>0 for all tt0; otherwise using (3.18), we obtain an exponential decay rate. We also define q(t) by

    q(t)=tt0g(τ)(3sxψx)(t)(3sxψx)(tτ)22dτ.

    Then q(t)CE(t), tt0. Now, we have that H is strictly convex on (0,r] (where r=g(t0)) and H(0)=0. Thus,

    H(στ)σH(τ),  0σ1 and τ(0,r]. (3.26)

    Using (3.26), condition (G2), (3.25), and Jensen's inequality, we get

    q(t)=1ηp(t)tt0p(t)(g(τ))η(3sxψx)(t)(3sxψx)(tτ)22dτ1ηp(t)tt0p(t)ξ(τ)H(g(τ))η(3sxψx)(t)(3sxψx)(tτ)22dτξ(t)ηp(t)tt0H(p(t)g(τ))η(3sxψx)(t)(3sxψx)(tτ)22dτξ(t)ηH(ηtt0g(τ)η(3sxψx)(t)(3sxψx)(tτ)22dτ)=ξ(t)ηˉH(ηtt0g(τ)η(3sxψx)(t)(3sxψx)(tτ)22dτ), (3.27)

    where ˉH is the convex extention of H on (0,+) (see remark 2.1). From (3.27), we have

    tt0g(τ)η(3sxψx)(t)(3sxψx)(tτ)22dτ1ηˉH1(ηq(t)ξ(t)).

    Therefore, (3.18) yields

    L1(t)βE(t)+CˉH1(ηq(t)ξ(t)),   tt0. (3.28)

    For r0<r, we define L3(t) by

    L3(t):=ˉH(r0E(t)E(0))L1(t)+E(t)E(t)

    since L1E. From (3.28) and using the fact that

    E(t)0, ˉH(t)>0, ˉH(t)>0,

    we obtain for all tt0

    L3(t)=r0E(t)E(0)ˉH(r0E(t)E(0))L1(t)+ˉH(r0E(t)E(0))L1(t)+E(t)βE(t)ˉH(r0E(t)E(0))+CˉH(r0E(t)E(0))ˉH1(ηq(t)ξ(t))+E(t). (3.29)

    Let us consider the convex conjugate of ˉH denoted by ˉH in the sense of Young (see [3] page 61-64). Thus,

    ˉH(τ)=τ(ˉH)1(τ)ˉH[(ˉH)(τ)] (3.30)

    and ˉH satisfies the generalized Young inequality

    ABˉH(A)+ˉH(B). (3.31)

    Let A=ˉH(r0E(t)E(0)) and B=ˉH1(μz(t)ξ(t)), It follows from (2.12) and (3.29)-(3.31) that

    L3(t)βE(t)ˉH(r0E(t)E(0))+CˉH(ˉH(r0E(t)E(0)))+Cηq(t)ξ(t)+E(t)βE(t)ˉH(r0E(t)E(0))+Cr0E(t)E(0)ˉH(r0E(t)E(0))+Cηq(t)ξ(t)+E(t). (3.32)

    Next, we multiply (3.32) by ξ(t) and recall that r0E(t)E(0)<r and

    ˉH(r0E(t)E(0))=H(r0E(t)E(0)),

    we arrive at

    ξ(t)L3(t)βξ(t)E(t)H(r0E(t)E(0))+Cr0E(t)E(0)ξ(t)H(r0E(t)E(0))+Cηq(t)+ξ(t)E(t)βξ(t)E(t)H(r0E(t)E(0))+Cr0E(t)E(0)ξ(t)H(r0E(t)E(0))CE(t). (3.33)

    Let L4(t)=ξ(t)L3(t)+CE(t). Since L3E, it follows that

    b0L4(t)E(t)b1L4(t), (3.34)

    for some b0,b1>0. Thus (3.33) gives

    L4(t)(βE(0)Cr0)ξ(t)E(t)E(0)ξ(t)H(r0E(t)E(0)), tt0.

    We select r0<r small enough so that βE(0)Cr0>0, we get

    L4(t)mξ(t)E(t)E(0)ξ(t)H(r0E(t)E(0))=mξ(t)H2(E(t)E(0)),  tt0, (3.35)

    for some constant m>0 and H2(τ)=τH(r0τ). We note here that

    H2(τ)=H(r0τ)+r0tH(r0τ),

    thus the strict convexity of H on (0,r], yields H2(τ)>0,H2(τ)>0 on (0,r]. Let

    F(t)=b0L4(t)E(0).

    From (3.34) and (3.35), we obtain

    F(t)E(t) (3.36)

    and

    F(t)=a0L4(t)(t)E(0)m1ξ(t)H2(F(t)), tt0. (3.37)

    Integrating (3.37) over (t0,t), we arrive at

    m1tt0ξ(τ)dτtt0F(τ)H2(F(τ))dτ=1r0r0F(t0)r0F(t)1τH(τ)dτ. (3.38)

    This implies

    F(t)1r0H11(¯m1tt0ξ(τ)dτ),   where  H1(t)=rt1τH(τ)dτ. (3.39)

    Using the properties of H, we see easily that H1 is strictly decreasing function on (0,r] and

    limt0H1(t)=+.

    Hence, (3.14) follows from (3.36) and (3.39). This completes the proof.

    Remark 3.1. The stability result in (3.1) is general and optimal in the sense that it agrees with the decay rate of g, see [10], Remark 2.3.

    Corollary 3.2. Suppose Gρ=DIρ, and (G1), and (G2) hold. Let the function H in (G2) be defined by

    H(τ)=τp,  1p<2, (3.40)

    then the solution energy (2.11) satisfies

    E(t)a2exp(a1t0ξ(τ)dτ),  for p=1,E(t)C(1+tt0ξ(τ)dτ)1p1,  for 1<p<2 (3.41)

    for some positive constants a2,a1 and C.

    In this subsection, we establish another stability result in the case non-equal speeds of wave propagation. To achieve this, we consider a stronger solution of (1.6). Let (w,3sψ,s,θ) be the strong solution of problem (1.6)–(1.8), then differentiation of 1.6 with respect to t gives

    {ρwttt+G(ψwx)xt=0,Iρ(3sψ)tttD(3sψ)xxt+t0g(τ)(3sψ)xxt(x,tτ)dτ+g(t)(3s0ψ0)xxG(ψwx)t=03Iρsttt3Dsxxt+3G(ψwx)t+4γst+δθxt=0,kθttλθxxt+δsxtt=0, (3.42)

    where (x,t)(0,1)×(0,+) and (3sψ)xx(x,0)=(3s0ψ0)xx. The modified energy functional associated to (3.42) is defined by

    E1(t)=12[ρwtt22+3Iρstt22+Iρ3sttψtt22+3Dsxt22+Gψtwxt22]+12[4γst22+kθt22+(Dt0g(τ)dτ)3sxtψxt22+(g(3sxtψxt))(t)]. (3.43)

    Lemma 3.2. Let (w,3sψ,s,θ) be the strong solution of problem (1.6)-(1.8). Then, the energy functional (3.43) satisfies, for all t0

    E1(t)=12(g(3sxtψxt))(t)12g(t)3sxtψxt22g(t)(3sttψtt),(3s0ψ0)xxλθxt22 (3.44)

    and

    E1(t)C(E1(0)+(3s0ψ0)xx22). (3.45)

    Proof. The proof of (3.44) follows the same steps as in the proof of Lemma 2.5. From (3.44), it is obvious that

    E1(t)g(t)(3sttψtt),(3s0ψ0)xx.

    So, using Cauchy-Schwarz inequality, we obtain

    E1(t)Iρg(t)23sttψtt22+g(t)2Iρ(3s0ψ0)xx22g(t)E1(t)+g(t)2Iρ(3s0ψ0)xx22. (3.46)

    This implies

    ddt(E1(t)et0g(τ)dτ)et0g(τ)dτg(t)2Iρ(3s0ψ0)xx22g(t)2Iρ(3s0ψ0)xx22 (3.47)

    Integrating (3.47) over (0,t) yields

    E1(t)e+0g(τ)dτE1(t)et0g(τ)dτE1(0)+12Iρ(t0g(τ)dτ)(3s0ψ0)xx22E1(0)+12Iρ(+0g(τ)dτ)(3s0ψ0)xx22. (3.48)

    Hence, (3.45) follows.

    Remark 3.2. Using Young's inequality, we observe from (3.44) and (3.45) that

    λθxt22=E1(t)+12(g(3sxtψxt))(t)12g(t)3sxtψxt22g(t)(3sttψtt),(3s0ψ0)xxE1(t)g(t)(3sttψtt),(3s0ψ0)xxE1(t)+g(t)(3sttψtt22+(3s0ψ0)xx22)E1(t)+g(t)(2IρE1(t)+(3s0ψ0)xx22)C(E1(t)+c1g(t)) (3.49)

    for some fixed positive constant c1. Similarly, we obtain

    0(g(3sxtψxt))(t)C(E1(t)+c1g(t)). (3.50)

    As in the case of equal-wave-speed of propagation, we define a Lyapunov functional

    ˜L(t)=˜NE(t)+6j=1~NjIj(t)+~N6I8(t), (3.51)

    where ˜N, ~Nj, j=1,2,3,4,5,6, are positive constants to be specified later and

    I8(t)=3λδ(IρGρD)10θxwxdx.

    Lemma 3.3. Suppose GρDIρ. Then, under suitable choice of ˜N, ~Nj, j=1,2,3,4,5,6, the Lyapunov functional ˜L satisfies, along the solution of (1.6), the estimate

    ˜L(t)˜βE(t)+12(g(3sxψx))(t)+C(E1(t)+c1g(t)), tt0, (3.52)

    for some positive constants ˜β and c1.

    Proof. Following the proof of Lemma 3.1, we end up with (3.52).

    Lemma 3.4. Suppose assumptions (G1) and (G2) hold and the function H in (G2) is linear. Let (w,3sψ,s,θ) be the strong solution of problem (1.6)-(1.8). Then,

    ξ(t)(g(3sxtψxt))(t)C(E1(t)+c1g(t)),  t0, (3.53)

    where c1 is a fixed positive constant.

    Proof. Using (3.50) and the fact that ξ is decreasing, we have

    ξ(t)(g(3sxtψxt))(t)=ξ(t)t0g(tτ)((3sxtψxt)(t)(3sxtψxt)(τ)22)dτt0ξ(tτ)g(tτ)((3sxtψxt)(t)(3sxtψxt)(τ)22)dτt0g(tτ)((3sxtψxt)(t)(3sxtψxt)(τ)22)dτ=(g(3sxtψxt))(t)C(E1(t)+c1g(t)). (3.54)

    Our stability result of this subsection is

    Theorem 3.3. Assume (G1) and (G2) hold and GρDIρ. Then, there exist positive constants a1,a2 and t2>t0 such that the energy solution (2.11) satisfies

    E(t)a2(tt0)H12(a1(tt0)tt2ξ(τ)dτ),t>t2,   where  H2(τ)=τH(τ). (3.55)

    Proof. Case 1: H is linear. Multiplying (3.52) by ξ(t) and using (G1), we get

    ξ(t)˜L(t)˜βξ(t)E(t)+12ξ(t)(g(3sxψx))(t)+Cξ(t)(E1(t)+c1g(t))˜βξ(t)E(t)CE(t)Cξ(0)E1(t)+ξ(0)c1g(t),  tt0

    Using the fact that ξ non-increasing, we obtain

    (ξ˜L+CE+E1)(t)˜βξ(t)E(t)+c2g(t),  tt0.

    for some fixed positive constant c2. This implies

    ˜βξ(t)E(t)(ξ˜L+CE+E1)(t)+c2g(t),  tt0. (3.56)

    Integrating (3.56) over (t0,t), using the fact that E is non-increasing and the inequality (3.45), we arrive at

    ˜βE(t)tt0ξ(τ)dτ˜βtt0ξ(τ)E(τ)dτ(ξ˜L+CE+E1)(t)+(ξ˜L+CE+E1)(t0)+c2tt0g(τ)dτ(ξ˜L+CE+E1)(0)+C(3s0ψ0)xx22+c20g(τ)dτ=(ξ˜L+CE+E1)(0)+C(3s0ψ0)xx22+c2(Dl0). (3.57)

    Thus, we have

    E(t)Ctt0ξ(τ)dτ,  tt0. (3.58)

    Case II: H is nonlinear. First, we observe from (3.52) that

    ˜L(t)˜βE(t)+12(g(3sxψx))(t)+C(E1(t)+c1g(t))˜βE(t)+C((g(3sxψx))(t)+(g(3sxtψxt))(t))+C(E1(t)+c1g(t)),  tt0. (3.59)

    From (2.12), (3.16) and (3.50), we have for any tt0

    t00g(τ)(3sxψx)(t)(3sxψx)(tτ)22dτ+t00g(τ)(3sxtψxt)(t)(3sxtψxt)(tτ)22dτ1ξ(t0)t00ξ(τ)g(τ)(3sxψx)(t)(3sxψx)(tτ)22dτ+1ξ(t0)t00ξ(τ)g(τ)(3sxtψxt)(t)(3sxtψxt)(tτ)22dτg(0)aξ(t0)t00g(τ)(3sxψx)(t)(3sxψx)(tτ)22dτg(0)aξ(t0)t00g(τ)(3sxtψxt)(t)(3sxtψxt)(tτ)22dτC(E(t)+E1(t))+c2g(t), (3.60)

    where c2 is a fixed positive constant. Substituting (3.60) into (3.59), we obtain for any tt0

    ˜L(t)˜βE(t)C(E(t)+E1(t))+c3g(t)+Ctt0g(τ)(3sxψx)(t)(3sxψx)(tτ)22dτ+Ctt0g(τ)(3sxtψxt)(t)(3sxtψxt)(tτ)22dτ, (3.61)

    where c3 is a fixed positive constant. Now, we define the functional Φ by

    Φ(t)=σtt0tt0(3sxψx)(t)(3sxψx)(tτ)22dτ+σtt0tt0(3sxtψxt)(t)(3sxtψxt)(tτ)22dτ,  t>t0. (3.62)

    Using (2.11), (2.12), (3.43) and (3.45), we easily get

    1tt0tt0(3sxψx)(t)(3sxψx)(tτ)22dτ+1tt0tt0(3sxtψxt)(t)(3sxtψxt)(tτ)22dτ2tt0tt0((3sxψx)(t)22+(3sxψx)(tτ)22)dτ+2tt0tt0((3sxtψxt)(t)22+(3sxtψxt)(tτ)22)dτ4l0(tt0)tt0(E(t)+E(tτ)+E1(t)+E1(tτ))dτ8l0(tt0)tt0(E(0)+C(E1(0)+(3s0ψ0)xx22))dτ8l0(E(0)+C(E1(0)+(3s0ψ0)xx22))<,  t>t0. (3.63)

    This last inequality allows us to choose 0<σ<1 such that

    Φ(t)<1,   t>t0. (3.64)

    Hence forth, we assume Φ(t)>0, otherwise, we get immediately from (3.61)

    E(t)Ctt0,   t>t0.

    Next, we define the functional μ by

    μ(t)=tt0g(τ)(3sxψx)(t)(3sxψx)(tτ)22dτtt0g(τ)(3sxtψxt)(t)(3sxtψxt)(tτ)22dτ (3.65)

    and observe that

    μ(t)C(E(t)+E1(t))+c4g(t),     t>t0, (3.66)

    where c4 is a fixed positive constant. The fact that H is strictly convex and H(0)=0 implies

    H(ντ)νH(τ),  0ν1 and τ(0,r]. (3.67)

    Using assumption (G1), (3.67), Jensen’s inequality and (3.64), we get for any t>t0

    μ(t)=1Φ(t)tt0Φ(t)g(τ)(3sxψx)(t)(3sxψx)(tτ)22dτ1Φ(t)tt0Φ(t)g(τ)(3sxtψxt)(t)(3sxtψxt)(tτ)22dτ1Φ(t)tt0Φ(t)ξ(τ)H(g(τ))(3sxψx)(t)(3sxψx)(tτ)22dτ+1Φ(t)tt0Φ(t)ξ(τ)H(g(τ))(3sxtψxt)(t)(3sxtψxt)(tτ)22dτξ(t)Φ(t)tt0H(Φ(t)g(τ))(3sxψx)(t)(3sxψx)(tτ)22dτ+ξ(t)Φ(t)tt0H(Φ(t)g(τ))(3sxtψxt)(t)(3sxtψxt)(tτ)22dτξ(t)(tt0)σH(σtt0tt0g(τ)(Ω1(tτ)+Ω2(tτ))dτ)=ξ(t)(tt0)σˉH(σtt0tt0(Ω1(tτ)+Ω2(tτ))dτ), (3.68)

    where

    Ω1(tτ)=(3sxψx)(t)(3sxψx)(tτ)22,Ω2(tτ)=(3sxtψxt)(t)(3sxtψxt)(tτ)22

    and ˉH is the C2 strictly increasing and convex extension of H on (0,+). This implies

    tt0g(τ)(3sxψx)(t)(3sxψx)(tτ)22dτ+tt0g(τ)(3sxtψxt)(t)(3sxtψxt)(tτ)22dτ(tt0)σˉH1(σμ(t)ξ(t)(tt0)),  t>t0. (3.69)

    Thus, the inequality (3.61) becomes

    ˜L(t)˜βE(t)C(E(t)+E1(t))+c3g(t)+C(tt0)σˉH1(σμ(t)ξ(t)(tt0)),  t>t0. (3.70)

    Let ˜L1(t):=˜L(t)+C(E(t)+E1(t)). Then (3.70) becomes

    ˜L1(t)˜βE(t)+C(tt0)σˉH1(σμ(t)ξ(t)(tt0))+c3g(t),  t>t0. (3.71)

    For 0<r1<r, we define the functional ˜L2 by

    ˜L2(t):=ˉH(r1tt0.E(t)E(0))˜L1(t),,  t>t0. (3.72)

    From (3.71) and the fact that

    E(t)0, ˉH(t)>0, ˉH(t)>0,

    we obtain, for all t>t0,

    ˜L2(t)=(r1(tt0)2.E(t)E(0)+r1tt0.E(t)E(0))ˉH(r1tt0.E(t)E(0))˜L1(t)+ˉH(r1tt0.E(t)E(0))˜L1(t)˜βE(t)ˉH(r1tt0.E(t)E(0))+c3g(t)ˉH(r1tt0.E(t)E(0))+C(tt0)σˉH(r1tt0.E(t)E(0))ˉH1(σμ(t)ξ(t)(tt0)). (3.73)

    Let ˉH be the convex conjugate of ˉH as in (3.30) and let

    A=ˉH(r1tt0.E(t)E(0))  and  B=ˉH1(σμ(t)ξ(t)(tt0)).

    Then, (3.30), (3.31) and (3.73) yield, for all t>t0,

    ˜L2(t)˜βE(t)ˉH(r1tt0.E(t)E(0))+c3g(t)ˉH(r1tt0.E(t)E(0))+C(tt0)σˉH(ˉH(r1tt0.E(t)E(0)))+C(tt0)σ.σμ(t)ξ(t)(tt0)˜βE(t)ˉH(r1tt0.E(t)E(0))+c3g(t)ˉH(r1tt0.E(t)E(0))+Cr1E(t)E(0)ˉH(r1tt0.E(t)E(0))+Cμ(t)ξ(t)(˜βE(0)Cr1)E(t)E(0)ˉH(r1tt0.E(t)E(0))+Cμ(t)ξ(t)+c3g(t)ˉH(r1tt0.E(t)E(0)) (3.74)

    By selecting r1 small enough so that (˜βE(0)Cr1)>0, we arrive at

    ˜L2(t)˜β2E(t)E(0)ˉH(r1tt0.E(t)E(0))+Cμ(t)ξ(t)+c3g(t)ˉH(r1tt0.E(t)E(0)),  t>t0, (3.75)

    for some positive constant ˜β2.

    Now, multiplying (3.75) by ξ(t) and recalling that r1E(t)E(0)<r, we arrive at

    ξ(t)˜L2(t)˜β2ξ(t)E(t)E(0)ˉH(r1tt0.E(t)E(0))+Cμ(t)+c3g(t)ξ(t)H(r1tt0.E(t)E(0))˜β2ξ(t)E(t)E(0)H(r1tt0.E(t)E(0))C(E(t)+E1(t))+c4g(t)+c3g(t)H(r1tt0.E(t)E(0)),  t>t0. (3.76)

    Since r1tt00 as t, there exists t2>t0 such that r1tt0<r1, whenever t>t2. Using this fact and observing that H strictly increasing, and E and ξ are non-decreasing, we get

    H(r1tt0.E(t)E(0))H(r1),  t>t2. (3.77)

    Using (3.77), it follows from (3.76) that

    ˜L3(t)˜β2ξ(t)E(t)E(0)H(r1tt0.E(t)E(0))+c5g(t),  t>t2, (3.78)

    where ˜L3=(ξ˜L2+CE+CE1) and c5>0 is a constant. Using the non-increasing property of ξ, we have

    ˜β2ξ(t)E(t)E(0)H(r1tt0.E(t)E(0))˜L3(t)+c5g(t),  t>t2. (3.79)

    Using the fact that E is non-increasing and H>0 we conclude that the map

    tE(t)H(r1tt0.E(t)E(0))

    is non-increasing. Therefore, integrating (3.79) over (t2,t) yields

    ˜β2E(t)E(0)H(r1tt0.E(t)E(0))tt2ξ(τ)dτ˜β2tt2ξ(τ)E(τ)E(0)H(r1τt0.E(τ)E(0))dτ˜L3(t)+˜L3(t2)+c5tt2g(τ)dτ˜L3(t2)+c50g(τ)dτ=˜L3(t2)+c5(bl0),  t>t2. (3.80)

    Next, we multiply both sides of (3.80) by 1tt0, for t>t2, we get

    ˜β2(tt0).E(t)E(0)H(r1tt0.E(t)E(0))tt2ξ(τ)dτ˜L3(t2)+c5(bl0)tt0,   t>t2. (3.81)

    Since H is strictly increasing, then H2(τ)=τH(τ) is a strictly increasing function. It follows from (3.81) that

    E(t)a2(tt0)H12(a1(tt0)tt2ξ(τ)dτ),  t>t2.

    for some positive constants a1 and a2. This completes the proof.

    (1). Let g(t)=aebt, t0,  a, b>0 are constants and a is chosen such that (G1) holds. Then

    g(t)=abebt=bH(g(t))  with  H(t)=t.

    Therefore, from (3.14), the energy function (2.11) satisfies

    E(t)a2eαt,  t0, where  α=ba1. (3.82)

    Also, for H2(τ)=τ, it follows from (3.55) that, there exists t2>0 such that the energy function (2.11) satisfies

    E(t)Ctt2,   t>t2, (3.83)

    for some positive constant C.

    (2).Let g(t)=ae(1+t)b, t0,  a>0, 0<b<1 are constants and a is chosen such that (G1) holds. Then,

    g(t)=ab(1+t)b1e(1+t)b=ξ(t)H(g(t)),

    where ξ(t)=b(1+t)b1 and H(t)=t. Thus, we get from (3.14) that

    E(t)a2ea1(1+t)b,  t0. (3.84)

    Likewise, for H2(t)=t, then estimate (3.55) implies there exists t2>0 such that the energy function (2.11) satisfies

    E(t)C(1+t)b,   t>t2, (3.85)

    for some positive constant C.

    (3).Let g(t)=a(1+t)b, t0,  a>0, b>1 are constants and a is chosen in such a way that (G1) holds. We have

    g(t)=ab(1+t)b+1=ξ(a(1+t)b)b+1b=ξgq(t)=ξH(g(t)),

    where

    H(t)=tq,  q=b+1b  satisfying  1<q<2  and  ξ=ba1b>0.

    Hence, we deduce from (3.41) that

    E(t)C(1+t)b,  t0. (3.86)

    Furthermore, for H2(t)=qtq, estimate (3.55) implies there exists t2>0 such that the energy function (2.11) satisfies

    E(t)C(1+t)(b1)/(b+1),  t>t2, (3.87)

    for some positive constant C.

    In this section, we prove the functionals Li,i=18, used in the proof of our stability results.

    Lemma 4.1. The functional I1(t) satisfies, along the solution of (1.6)(1.8), for all tt0>0 and for any ϵ1,ϵ2>0, the estimate

    I1(t)Iρg023stψt22+ϵ13sxψx22+ϵ2ψwx22+CCα(1+1ϵ1+1ϵ2)(h(3sxψx))(t), (4.1)

    where g0=t00g(τ)dτt0g(τ)dτ.

    Proof. Differentiating I1(t), using (1.6)2 and integrating by part, we have

    I1(t)=Iρ10(3stψt)t0g(tτ)((3sψ)(t)(3sψ)(τ))dτdx+D(t)10(3sxψx)t0g(tτ)((3sxψx)(t)(3sxψx)(τ))dτdx+10(t0g(tτ)((3sxψx)(t)(3sxψx)(τ))dτ)2dxIρ(t0g(τ)dτ)10(3stψt)2dxG10(ψwx)t0g(tτ)((3sψ)(t)(3sψ)(τ))dτdx, (4.2)

    where D(t)=(Dt0g(τ)dτ). Now, we estimate the terms on the right hand-side of (4.2). Exploiting Young's and Poincaré's inequalities, Lemmas 2.1-2.6 and performing similar computations as in (2.8), we have for any ϵ1>0,

    D(t)10(3sxψx)t0g(tτ)((3sxψx)(t)(3sxψx)(τ))dτdxϵ1(3sxψx22+CCαϵ1(h(3sxψx))(t) (4.3)

    and

    10(t0g(tτ)((3sxψx)(t)(3sxψx)(τ))dτ)2dxCα(h(3sxψx))(t). (4.4)

    Also, for δ1>0, we have

    Iρ10(3stψt)t0g(tτ)((3sψ)(t)(3sψ)(τ))dτdx=Iρ10(3stψt)t0h(tτ)((3sψ)(t)(3sψ)(τ))dτdxIρα10(3stψt)t0g(tτ)((3sψ)(t)(3sψ)(τ))dτdxδ13stψt22+I2ρ2δ110(t0h(tτ)((3sψ)(t)(3sψ)(τ))dτ)2dx+α2I2ρ2δ110(t0g(tτ)((3sψ)(t)(3sψ)(τ))dτ)2dxδ13stψt22+I2ρ2δ1(t0h(τ)dτ)(h(3sψ))(t)+α2I2ρCα2δ1(h(3sψ))(t)δ13stψt22+C(Cα+1)δ1(h(3sψ)x)(t). (4.5)

    For the last term, we have

    G10(ψwx)t0g(tτ)((3sψ)(t)(3sψ)(τ))dτdxϵ2ψwx22+G2Cα4ϵ2(h(3sψ)x)(t). (4.6)

    Combination of (4.2)-(4.6) lead to

    I1(t)(Iρt0g(τ)dτδ1)3wtψt22+ϵ13sxψx22+ϵ2ψwx22+CCα(1+1δ1+1ϵ1+1ϵ2)(h(3sxψx))(t). (4.7)

    Since g(0)>0 and g is continuous. Thus for any tt0>0, we get

    t0g(τ)dτt00g(τ)dτ=g0>0. (4.8)

    We select δ1=Iρg02 to get (4.1).

    Lemma 4.2. The functional I2(t) satisfies, along the solution of (1.6)(1.8) and for any δ4>0, the estimate

    I2(t)3Dsx223γs22+δ4wt22+C(1+1δ4)st22+Cθx22,  t0. (4.9)

    Proof. Differentiation of I2(t), using (1.6)1 and (1.6)3 and integration by part, leads to

    I2(t)=3Iρst223Dsx224γs22δt0sθxdx+3ρ10wtx0st(y)dydx.

    Applying Cauchy-Schwarz and Young's inequalities and (2.5), we get for any δ4>0,

    I2(t)3Iρst223Dsx224γs22+γs22+δ24γθx22+δ4wt22+9ρ24δ410(x0st(y)dy)2dx3Dsx223γs22+δ4wt22+C(1+1δ4)st22+Cθx22.

    This completes the proof.

    Lemma 4.3. The functional I3(t) satisfies, along the solution of (1.6)(1.8) and for any ϵ3>0, the estimate

    I3(t)δIρ2st22+ϵ3sx22+ϵ3ψwx22+C(1+1ϵ3)θx22,  t0. (4.10)

    Proof. Differentiation of I3, using (1.6)3, (1.6)4 and integration by parts, yields

    I3(t)=3λIρ10θxstdx3Iρδst223kD10θsxdx+kδθ22+3kG10θx0(ψwy)(y)dydx+4γk10θt0s(y)dydx.

    Using Cauchy-Schwarz, Young's and Poincaré's inequalities together with Lemmas 2.1-2.6, we have

    I3(t)δ2st22+Cδ2θx223Iρδst22+ϵ32sx22+C(1+1ϵ3)θ22+ϵ310(x0(ψwy)(y)dy)2dx+ϵ3210(x0s(y)dy)2dxδ2st22+Cδ2θx223Iρδst22+ϵ3sx22+ϵ3ψwx22+C(1+1ϵ3)θx22.

    We choose δ2=5Iρδ2 to get (4.10).

    Lemma 4.4. The functional I4(t) satisfies, along the solution of (1.6)(1.8) and for any ϵ4>0, the estimate

    I4(t)ρwt22+ϵ43sxψx22+Csx22+Cϵ4ψwx22,  t0. (4.11)

    Proof. Using (1.6)1 and integration by parts, we have

    I4(t)=ρwt22G10(ψwx)wxdx.

    We note that wx=(ψwx)(3sψ)+3s to arrive at

    I4(t)=ρwt22+Gψwx22+G10(ψwx)(3sψ)dx3G10(ψwx)sdx.

    It follows from Young's and Poincaré's inequalities that

    I4(t)ρwt22+Gψwx22+ϵ43sψ22+Cϵ4ψwx22+3G2ψwx22+3G2s22ρwt22+Gψwx22+ϵ43sxψx22+Csx22+C(1+1ϵ4)ψwx22.

    This completes the proof.

    Lemma 4.5. The functional I5(t) satisfies, along the solution of (1.6)(1.8) and for any 0<α<1, the estimate

    I5(t)l043sxψx22+Iρ3stψt22+Cψwx22+CCα(h(3sxψx))(t). (4.12)

    Proof. Differentiating I5, using (1.6)2, we arrive at

    I5(t)=Iρ3stψt22(Dt0g(τ)dτ)3sxψx22+G10(3sψ)(ψwx)dx+10(3sxψx)t0g(tτ)((3sxψx)(x,τ)(3sxψx)(x,t))dτdx.

    Applying Lemmas 2.1-2.6, Cauchy-Schwarz, Young's and Poincaré's inequalities, we obtain any δ3>0

    I5(t)Iρ3stψt22l03sxψx22+δ33sxψx22+G24δ3ψwx22+l023sxψx22+12l0Cα(h(3sxψx))(t). (4.13)

    We select δ3=l04 and obtain the desired result.

    Lemma 4.6. The functional I6(t) satisfies, along the solution of (1.6)(1.8) and for any for any ϵ1, the estimate

    I6(t)G2ψwx22+ϵ13stψt22+C(1+1ϵ1)st22+Csx22+Cθx22+3(IρGρD)10wtsxtdx,  t0. (4.14)

    Proof. Differentiating I6(t), using (1.6)1 and (1.6)3 and integration by parts, we obtain

    I6(t)=3G2ψwx224γG10(ψwx)sdxδG10(ψwx)θxdx3IρGt0(3stψt)stdx+9IρGst22+3(IρGρD)10wtsxtdx. (4.15)

    Young's and Poincaré's inequalities give

    4γG10(ψwx)sdxG2ψwx22+4γ2Cpsx22,δG10(ψwx)θxdxG2ψwx22+δ24θx22,3IρGt0(3stψt)stdxϵ13stψt22+(3IρG)2ϵ1st22. (4.16)

    Substituting (4.16) into (4.15), we obtain (4.14). This completes the proof.

    Lemma 4.7. The functional I7(t) satisfies, along the solution of (1.6)(1.8), the estimate

    I7(t)3(Dl0)3sxψx2212(g(3sxψx))(t), t0. (4.17)

    Proof. Differentiate I7(t) and use the fact that J(t)=g(t) to get

    I7(t)=10t0J(tτ)(3sxψx)2(τ)dτdx+J(0)3sxψx22=(g(3sxψx))(t)+J(t)3sxψx22 210(3sxψx)t0g(tτ)((3sxψx)(τ)(3sxψx)(t))dx. (4.18)

    Using Cauchy-Schwarz and (G1), we have

    210(3sxψx)t0g(tτ)((3sxψx)(τ)(3sxψx)(t))2(Dl0)3sxψx22+t0g(τ)dτ2(Dl0)(g(3sxψx))(t)2(Dl0)3sxψx22+12(g(3sxψx))(t) (4.19)

    Thus, we get

    I7(t)2(Dl0)3sxψx2212(g(3sxψx))(t)+J(t)3sxψx22. (4.20)

    Since J is decreasing (J(t)=g(t)0), so J(t)J(0)=Dl0. Hence, we arrive at

    I7(t)3(Dl0)3sxψx2212(g(3sxψx))(t).

    The next lemma is used only in the proof of the stability result for nonequal-wave-speed of propagation.

    Lemma 4.8. Let (w,3sψ,s,θ) be the strong solution of problem (1.6). Then, for any positive numbers σ1,σ2,σ3, the functional I8(t) satisfies

    I8(t)3(IρGρD)10wtsxtdx+σ1wt22+σ2ψwx22+σ33sxψx22+Csx22+C(1+1σ1+1σ2+1σ3)θxt22,  tt0. (4.21)

    Proof. Differentiation of I8, using integration by part and the boundary condition give

    I8(t)=3λδ(IρGρD)10θxwxtdx+3λδ(IρGρD)10θxtwxdx=3λδ(IρGρD)[10θxxwtdx]+3λδ(IρGρD)10θxtwxdx. (4.22)

    We note that wx=(ψwx)(3sψ)+3s and from (1.6)4, λθxx=kθt+δsxt. So, (4.22) becomes

    I8(t)=3δ(IρGρD)k10θtwtdx3(IρGρD)10sxtwtdx+9λδ(IρGρD)10θxtsdx3λδ(IρGρD)10θxt(ψwx)dx3λδ(IρGρD)10θxt(3sψ)dx (4.23)

    Using Young's and Poincaré's inequalities, we have for any positive numbers σ1,σ2,σ3,

    3δ(IρGρD)10θtwtdxσ1wt22+Cσ1θxt22,3λδ(IρGρD)10θxt(ψwx)dxσ2ψwx22+Cσ2θxt22,3λδ(IρGρD)10θxt(3sψ)dxσ33sxψx22+Cσ3θxt22,9λδ(IρGρD)10θxtsdxCsx22+Cθxt22. (4.24)

    Substituting (4.24) into (4.23), we obtain (4.21).

    In this paper, we have established a general and optimal stability estimates for a thermoelastic Laminated system, where the heat conduction is given by Fourier's Law and memory as the only source of damping. Our results are established under weaker conditions on the memory and physical parameters. From our results, we saw that the decay rate is faster provided the wave speeds of the first two equations of the system are equal (see (1.3)). A similar result was established recently in [19] when the heat conduction is given by Maxwell-Cattaneo's Law. An interesting case is when the kernel memory term is couple with the first or third equations in system (1.6). Our expectation is that the stability in both cases will depend on the speed of wave propagation.

    The authors appreciate the continuous support of University of Hafr Al Batin, KFUPM and University of Sharjah. The first and second authors are supported by University of Hafr Al Batin under project #G-106-2020 . The third author is sponsored by KFUPM under project #S B181018.

    The authors declare no conflict of interest



    [1] Pinheiro LRS, Gradissimo DG, Xavier LP, et al. (2022) Degradation of azo dyes: Bacterial potential for bioremediation. Sustainability 14: 1510. https://doi.org/10.3390/su14031510. doi: 10.3390/su14031510
    [2] Rauf MA, Ashraf SS (2009) Review: Fundamental principles and application of heterogeneous photocatalytic degradation of dyes in solution. Chem Eng J 151: 10–18. https://doi.org/10.1016/j.cej.2009.02.026 doi: 10.1016/j.cej.2009.02.026
    [3] Zhuang Y, Zhu Q, Li G, et al. (2022) Photocatalytic degradation of organic dyes using covalent triazine-based framework. Mater Res Bulletin 146: 111619. https://doi/org/10/1016/j/materresbull.2021.111619. doi: 10.1016/j/materresbull.2021.111619
    [4] Sibhatu AS, Weldegebrieal KG, Sgaradevan S (2022) Photocatalytic activity of CuO nanoparticles for organic and inorganic pollutants removal in wastewater remediation. Chemosphere 300: 134623. https://doi.org/10.1016/j.chemosphere.2022.134623 doi: 10.1016/j.chemosphere.2022.134623
    [5] Fouda A, Salam S, Wassel AR, et al. (2020) Optimization of green biosynthesized visible light active CuO/ZnO nano-photocatalysts for the degradation of organic methylene blue dye. Hélion 6: e04896. https://doi.org/10.1016/jhelion.2020.e04896.
    [6] Lacombe S, Tran-thi T, Guillard C, et al. (2007) La photocalyse pour l'elimination des polluants. Actualités chimique 308: 79–93.
    [7] Liu X, C Chen, Zhao Z, et al. (2013) A review on the synthesis of manganese oxide nanomaterials and their applications on lithium-ion batteries. J Nanomater 2013: 736375. http://dx.doi.org/10.1155/2013/736375 doi: 10.1155/2013/736375
    [8] Naika HR, Lingaraju K, Manjunath K, et al. (2015) Green synthesis of CuO nanoparticles using Gloriosa superba L. extract and their antibacterial activity. J Taibah Univ Sci 9: 7–12. https://doi.org/10.1016/j.jtusci.2014.04.006. doi: 10.1016/j.jtusci.2014.04.006
    [9] Ahmad MM, Kotb HM, Mushta S, et al. (2022) Green synthesis of Mn + Cu bimetallic nanoparticles using vinca rosea extract and their antioxidant, antibacterial, and catalytic activities. Crystals 12: 72. https://doi.org/10.3390/cryst12010072. doi: 10.3390/cryst12010072
    [10] Basavegowda N, Baek K (2021) Multimetallic nanoparticles as alternative antimicrobial agents: Challenges and perspectives. Molecules 26: 912. https://doi.org/10.3390/molecules26040912. doi: 10.3390/molecules26040912
    [11] Iqbal M, Thebo AA, Shah AH, et al. (2016) Influence of Mn-doping on the photocatalytic and solar cell efficiency of CuO nanowires. Inorg Chem Commun 76: 71–76. https://doi.org/10.1016/j.inoche.2016.11.023. doi: 10.1016/j.inoche.2016.11.023
    [12] Pramothkumar A, Senthilkumar N, Mercy Gnana Malar KC, et al. (2019) A comparative analysis on the dye degradation efciency of pure, Co, Ni and Mn‑doped CuO nanoparticles. J Mater Sci-Mater El 30: 19043–19059. https://doi.org/10.1007/s10854-019-02262-4. doi: 10.1007/s10854-019-02262-4
    [13] Vindhya PS, Kavitha VT (2022) Leaf extract-mediated synthesis of Mn-doped CuO nanoparticles for antimicrobial, antioxidant and photocatalytic applications. Chem Pap. https://doi.org/10.1007/s11696-022-02631-0.
    [14] Kabengele CN, Kasiama GN, Ngoyi EM, et al. (2022) Secondary metabolites and mineral elements of Manotes expansa and Aframomum alboviolaceum leaves collected in the democratic republic of Congo. ARRB 37: 57–63.
    [15] Rizwana H, Alwhibi MS, Al-Judaie RA, et al. (2022) Sunlight-mediated green synthesis of silver nanoparticles using the berries of Ribes rubrum (Red Currants): characterization and evaluation of their antifungal and antibacterial activities. Molecules 27: 2186. https://doi.org/10.3390/molecules27072186. doi: 10.3390/molecules27072186
    [16] Chen LQ, Li Fang, Ling J, et al. (2015) Nanotoxicity of silver nanoparticles to red blood cells: size dependent adsorption, uptake, and hemolytic activity. Chem Res Toxicol 28: 501–509. https://doi.org/10.1021/tx500479 doi: 10.1021/tx500479
    [17] Pandey S, Singh S (2020) Eco-friendly nanocomposite and properties of manganese nanoparticles using UV-vis and IR fourier spectrum. IJISRT 5: 770–773.
    [18] Shah M, Fawcett D, Sharma S, et al. (2015) Review green synthesis of metallic nanoparticles via biological entities. Materials 8: 7278–7308. https://doi.org/10.3390/ma8115377. doi: 10.3390/ma8115377
    [19] El-seedi, El-Shabasy RM, Khalifa SAM, et al. (2019) Metal nanoparticles fabricated by green chemistry using natural extracts: biosynthesis, mechanisms, and applications. RSC Adv 24539–24559. https://doi.org/10.1039/C9RA02225B
    [20] Makarov VV, Love AJ, Sinitsyna OV, et al. (2014) Green nanotechnologies: Synthesis of metal nanoparticles using plants. Acta Naturae 6: 35–44. https://doi.org/10.32607/20758251-2014-6-1-35-44 doi: 10.32607/20758251-2014-6-1-35-44
    [21] Desai R, Mankad V, Gupta SG, et al. (2012) Size distribution of silver nanoparticles: UV-visible spectroscopic assessment. Nanosci Nanotechnol Let 4: 30–34. https://doi.org/10.1166/nnl.2012.1278 doi: 10.1166/nnl.2012.1278
    [22] Yeshchenko OA, Bondarchuk IS, Gurin VS, et al. (2013) Temperature dependence of the surface plasmon resonance in gold nanoparticles. Surf Sci 608: 275–281, http://dx.doi.org/10.1016/j.susc.2012.10.019. doi: 10.1016/j.susc.2012.10.019
    [23] Seifipour R, Nozari M, Pishkar L (2020) Green synthesis of silver nanoparticles using Tragopogon Collinus leaf extract and study of their antibacterial effects. JIOPM 30: 2926–2936. https://doi.org/10.1007/s10904-020-01441-9 doi: 10.1007/s10904-020-01441-9
    [24] Vidhu VK, Aromal SA, Philip D (2011) Green synthesis of silver nanoparticles using Macrotyloma uniform. Spectrochim Acta A 83: 392–397. https://doi.org/10.1016/j.saa.2011.08.051 doi: 10.1016/j.saa.2011.08.051
    [25] Berta L, Coman NA, Rusu A, et al. (2021) A review on plant-Mediated synthesis of Bimetallic nanoparticles, characterization and their biological applications. Materials 14: 7677. https://doi.org/10.3390/ma14247677 doi: 10.3390/ma14247677
    [26] Pinto VV, Ferreira MJ, Silva R, et al. (2010) Long time effect on the stability of silver nanoparticles in aqueous medium: effect of synthesis and storage conditions. Colloid Surface A 364: 19–25. https://doi.org/10.1016/j.colsurfa.2010.04.015 doi: 10.1016/j.colsurfa.2010.04.015
    [27] Azeez F, Al-Hetlani E, Arafa M (2018) The effect of surface charge on photocatalytic degradation of Methylene Blue dye using chargeable titania nanoparticles. Sci Rep 2018: 7104. https://doi.org/10.1038/s4158-018-15673-5. doi: 10.1038/s4158-018-15673-5
    [28] Taylor MG, Augustin N, Gounaris CE, et al. (2015) Catalyst design based on morphology and environment dependent adsorption on metal nanoparticles. ACS Catal 20155: 6296–6301. https://doi.org/10.1021/acscatal.5b01696 doi: 10.1021/acscatal.5b01696
    [29] Chanu LA, Singh WJ, Singh KJ, et al. (2019) Effect of operational parameters on the photocatalytic degradation of Methylene blue dye solution using manganese doped ZnO nanoparticles. Results Phys 12: 1230–1237. https://doi.org/10.1016/j.rinp.2018.12.089 doi: 10.1016/j.rinp.2018.12.089
    [30] Dobrovolskaia MA, Clogston JD, Neun BW (2008) Method for analysis of nanoparticle hemolytic properties in vitro. Nano Lett 8: 2180–2187. https://doi.org/10.1021/nl0805615 doi: 10.1021/nl0805615
    [31] Gabor F (2011) Characterization of nanoparticles intended for drug delivery. Sci Pharm 79: 701–702.
    [32] Gul A, Shaheen A, Ahmad I, et al. (2021) Green synthesis, characterization, enzyme inhibition, antimicrobial potential, and cytotoxic activity of plant mediated silver nanoparticle using Ricinus communis leaf and root extracts. Biomolecules 11: 206. https://doi.org/10.3390/biom11020206. doi: 10.3390/biom11020206
  • This article has been cited by:

    1. Yuxia Guo, Yichen Hu, Infinitely many solutions for Hamiltonian system with critical growth, 2024, 13, 2191-950X, 10.1515/anona-2023-0134
    2. Xingyue He, Chenghua Gao, Jingjing Wang, k-convex solutions for multiparameter Dirichlet systems with k-Hessian operator and Lane-Emden type nonlinearities, 2024, 13, 2191-950X, 10.1515/anona-2023-0136
    3. Hongying Jiao, Shuhai Zhu, Jinguo Zhang, Existence of infinitely many solutions for critical sub-elliptic systems via genus theory, 2024, 16, 2836-3310, 237, 10.3934/cam.2024011
    4. Jinli Yang, Jiajing Miao, Algebraic Schouten solitons of Lorentzian Lie groups with Yano connections, 2023, 15, 2836-3310, 763, 10.3934/cam.2023037
  • Reader Comments
  • © 2023 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(1964) PDF downloads(101) Cited by(3)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog