Research article Special Issues

Fractional infinite time-delay evolution equations with non-instantaneous impulsive

  • This dissertation is regarded to investigate the system of infinite time-delay and non-instantaneous impulsive to fractional evolution equations containing an infinitesimal generator operator. It turns out that its mild solution is existed and is unique. Our model is built using a fractional Caputo approach of order lies between 1 and 2. To get the mild solution, the families associated with cosine and sine which are linear strongly continuous bounded operators, are provided. It is common to use Krasnoselskii's theorem and the Banach contraction mapping principle to prove the existence and uniqueness of the mild solution. To confirm that our results are applicable, an illustrative example is introduced.

    Citation: Ahmed Salem, Kholoud N. Alharbi. Fractional infinite time-delay evolution equations with non-instantaneous impulsive[J]. AIMS Mathematics, 2023, 8(6): 12943-12963. doi: 10.3934/math.2023652

    Related Papers:

    [1] A.G. Ibrahim, A.A. Elmandouh . Existence and stability of solutions of ψ-Hilfer fractional functional differential inclusions with non-instantaneous impulses. AIMS Mathematics, 2021, 6(10): 10802-10832. doi: 10.3934/math.2021628
    [2] Muneerah Al Nuwairan, Ahmed Gamal Ibrahim . Nonlocal impulsive differential equations and inclusions involving Atangana-Baleanu fractional derivative in infinite dimensional spaces. AIMS Mathematics, 2023, 8(5): 11752-11780. doi: 10.3934/math.2023595
    [3] Thabet Abdeljawad, Pshtiwan Othman Mohammed, Hari Mohan Srivastava, Eman Al-Sarairah, Artion Kashuri, Kamsing Nonlaopon . Some novel existence and uniqueness results for the Hilfer fractional integro-differential equations with non-instantaneous impulsive multi-point boundary conditions and their application. AIMS Mathematics, 2023, 8(2): 3469-3483. doi: 10.3934/math.2023177
    [4] Mohamed Adel, M. Elsaid Ramadan, Hijaz Ahmad, Thongchai Botmart . Sobolev-type nonlinear Hilfer fractional stochastic differential equations with noninstantaneous impulsive. AIMS Mathematics, 2022, 7(11): 20105-20125. doi: 10.3934/math.20221100
    [5] Muneerah AL Nuwairan, Ahmed Gamal Ibrahim . The weighted generalized Atangana-Baleanu fractional derivative in banach spaces- definition and applications. AIMS Mathematics, 2024, 9(12): 36293-36335. doi: 10.3934/math.20241722
    [6] Zhilin Li, Guoping Chen, Weiwei Long, Xinyuan Pan . Variational approach to p-Laplacian fractional differential equations with instantaneous and non-instantaneous impulses. AIMS Mathematics, 2022, 7(9): 16986-17000. doi: 10.3934/math.2022933
    [7] Velusamy Kavitha, Mani Mallika Arjunan, Dumitru Baleanu . Non-instantaneous impulsive fractional-order delay differential systems with Mittag-Leffler kernel. AIMS Mathematics, 2022, 7(5): 9353-9372. doi: 10.3934/math.2022519
    [8] Anil Kumar, Muslim Malik, Mohammad Sajid, Dumitru Baleanu . Existence of local and global solutions to fractional order fuzzy delay differential equation with non-instantaneous impulses. AIMS Mathematics, 2022, 7(2): 2348-2369. doi: 10.3934/math.2022133
    [9] Noorah Mshary, Hamdy M. Ahmed, Ahmed S. Ghanem . Existence and controllability of nonlinear evolution equation involving Hilfer fractional derivative with noise and impulsive effect via Rosenblatt process and Poisson jumps. AIMS Mathematics, 2024, 9(4): 9746-9769. doi: 10.3934/math.2024477
    [10] Ye Li, Biao Qu . Mild solutions for fractional non-instantaneous impulses integro-differential equations with nonlocal conditions. AIMS Mathematics, 2024, 9(5): 12057-12071. doi: 10.3934/math.2024589
  • This dissertation is regarded to investigate the system of infinite time-delay and non-instantaneous impulsive to fractional evolution equations containing an infinitesimal generator operator. It turns out that its mild solution is existed and is unique. Our model is built using a fractional Caputo approach of order lies between 1 and 2. To get the mild solution, the families associated with cosine and sine which are linear strongly continuous bounded operators, are provided. It is common to use Krasnoselskii's theorem and the Banach contraction mapping principle to prove the existence and uniqueness of the mild solution. To confirm that our results are applicable, an illustrative example is introduced.



    Fractional differential equations are becoming a considerably more important and popular topic. In order to specify the so-called fractional differential equations, the conventional integer order derivative is generalized to arbitrary order. Fractional differential equations have been extensively employed to explain a variety of physical processes because of the effective memory function of the fractional derivative, such as seepage flow in porous media, fluid dynamic, and traffic models. Also, there are several applications of fractional differential equations in control theory, polymer rheology, aerodynamics, physics, chemistry, biology, and other exciting conceptual advancements (see [1,2,3,4] and their references).

    In the real world, there are numerous processes and phenomena that are influenced by transient external factors as they evolve. When compared to the whole duration of the occurrences and processes being researched, their duration is tiny. As a result, it is reasonable to believe that these exterior impacts are "instantaneous", or take the shape of impulses. Differential equations including the impulse effect, or impulsive differential equations, seem to be a plausible explanation of the known evolution processes of various real-world issues. The impulsive differential equations have been studied as the subject of numerous excellent monographs [5,6].

    Differential equations are used as representations for many processes in applied sciences research. There are a variety of classical mechanics that experience abrupt changes in their states, such as biological systems (heartbeats and blood flow), mechanical systems with impact, radiophysics, pharmacokinetics, population dynamics, mathematical economics, ecology, industrial robotics, biotechnology processes, etc [7,8]. The systems of differential equations with impulses are suitable mathematical models for such phenomena. Impulsive differential equations essentially have three parts: An impulse equation simulates an impulsive leap that is described by a jump function at the moment of impulse occurs, a continuous-time differential equation determines the state of the system between impulses, and jump criteria identifies a set of jump occurrences [9,10,11].

    Furthermore, there are various models that have been developed in many fields like biology, economics, and materials science where the rate of change at time t depends not only on the system's current state but also on its history over a period of time [tτ,t] [12,13,14,15]. These models have evolutionary equations with delay, which describe them mathematically. Equations with infinite delay are produced by the more generic if we take τ=.

    Physics provides a compelling justification for studying the nonlocal partial differential equation. Fractional derivatives in space and time are used in abstract partial differential equations such as fractional diffusion equations. They may be used to simulate anomalous diffusion, in which a particle plume spreads differently than the traditional diffusion equation would suggest. By replacing the second-order space derivative in the classic diffusion equation by an infinitesimal generator operator of strongly continuous C0 semi-group or cosine functions, the time fractional evolution equation is derived [16,17,18].

    In [19], Kumar and Pandey attempted to examine the results of the existence of a solution to a class of FDEs (the fractional calculus due to Atangana-Baleanu) of the sort

    {ABCDρv(r)=Fv(r)+h(r,v(r)),rmj=0(sj,rj+1],v(r)=δj(r,v(r)),rmj=1(rj,sj],v(0)=v0f(v),

    for all vD(F) (the domain of F), where ABCDρ is the ABC fractional derivative of order ρ(0,1), F:D(F)XX is a generator of ρ-resolvent operator {Sρ(r)}r0 on a Banach space (X,), J=[0,d], 0=r0<r1<r2<<rm<rm+1 and sj(rj,rj+1) for all j=1,2,,m;mN. The functions h:mk=0(sj,rj+1]×XX and f:XX are given continuous functions and δj:(rj,sj]×XX are non-instantaneous impulsive functions for each j=1,2,,m;mN and voX.

    Recently, Kavitha et al. [20] examined the existence of solutions for a class of non-instantaneous impulses and infinite delay of fractional differential equations within the Mittag-Leffler kernel of the kind

    {ABCDvp(ξ)=Bp(ξ)+F(ξ,pξ),ξml=0(sl,ξl+1],p(ξ)=κl(ξ,pξ),ξml=1(ξl,sl],p(ξ)=ϕ(ξ),ξ(,0],

    where the fractional order v(0,1), B:D(B)EE is infinitesimal generator of an ρ-resolvent operator {Sρ(ξ)}ξ0 on a Banach space (E,, K=[0,b], 0=ξ0=s0<ξ1s1<ξ2<<ξq<ξq+1=b and sl(ξl,ξl+1) for all l=1,2,,q;qN. The function F:ql=0(sl,ξl+1]×AE satisfies Caratheodory conditions and the functions κl:(ξl,sl]×AE are non-instantaneous impulsive functions for each l=1,2,,q. They considered that pξ:(,0]E such that pξ(x)=p(ξ+x) for all x0 and ϕA where A is an abstract phase space.

    In light of the foregoing, in this publication, we examine the existence results for a class of fractional-order non-instantaneous impulses functional evolution equations.

    Consider the fractional semilinear evolution of the form

    {u(t)=ϕ(t),t(,0],cDαtu(t)=Au(t)+h(t,ut),tmk=0(sk,tk+1],u(t)=μk(t,u(t)),tmk=1(tk,sk],u(t)=ξk(t,u(t)),tmk=1(tk,sk],u(0)=u0, (1.1)

    where cDαt is the fractional derivative due to Caputo of order 1<α<2 and J=[0,a] is operational interval. Here, h:J×PhX is a given function satisfying some assumptions that will be determined later, where Ph is an abstract phase space and X is a Banach space. The functions μk,ξkC((tk,sk]×X;X) for all k=1,2,m;mN reflect the impulsive circumstances and 0=t0=s0<t1s1t2<<tmsmtm+1=a are pre-fixed numbers. The history function ut:(,0]X is an element of Ph and defined by ut(θ)=u(t+θ),θ(,0].

    The closed operator A is an infinitesimal generator of a uniformly bounded family of strongly continuous cosine operators {R(t)}tR, which is defined on a Banach space X. The Banach space of continuous and bounded functions from (,a] into X provided with the topology of uniform convergence is denoted by C=Ca((,a],X) with the norm

    uC=supt(,a]|u(t)|.

    As {R(t)}tR is a cosine family on X, then there exists ϖ1 [21] such that

    R(t)ϖ. (1.2)

    The rest of the text is structured as follows. In section 2, we give some fundamental concepts and lemmas related to our study. In section 3, we formulate the mild solution of (1.1) by considering that operator A is an infinitesimal generator of strongly continuous cosine functions {R(t)}tR. By using the fixed-point theorem, Section 4 presents our study outcomes. An instance is provided in Section 5 to be an application.

    In this section, we present some concepts and definitions related to the components of the research paper such as fractional calculus, cosine and sine families operators, and abstract phase space. Also, some lemmas that give helpful results to prove the main results of this contribution, are provided.

    The definitions of R-L integral and Caputo derivative and the important related lemma are introduced as follow.

    Definition 2.1. (Caputo derivative [22]) Let 1<q;N and x:[a,b]R (<a<b<) be nth continuously differentiable function. Then, the left derivative of fractional order q due to Caputo is presented as

    cDqax(s)=1Γ(q)sa(st)q1x()(t)dt,s[a,b].

    Definition 2.2. (Riemann-Liouville fractional integral [22]) The left R-L fractional integral of the integrable function x over the interval [a,b] is derived as

    Iqax(s)=1Γ(q)sa(st)q1x(t)dt,q>0,s[a,b].

    Lemma 2.1. [23] Let N,1<q and x(s) be nth continuously differentiable function over the interval [a,b]. Then,

    IqacDqax(s)=x(s)+a0+a1(sa)++a1(sa)1,s[a,b].

    Definition 2.3. (Atangana-Baleanu fractional derivative in Caputo sense [19]) ABC-derivative for the order α[0,1] and x(s)H1(a,b),a<b is given by

    ABCDαx(s)=M(α)(1α)saEα[α(sr)α1α]x(r)dr,

    where Eα() and H1 are the Mittag-Leffler function and the non-typical Banach space defined, respectively, as

    Eα(z)=i=0ziΓ(αi+1),(α)>0,zC,H1(Ω)={η(s)|η(s),Dρη(s)L2(Ω),ρ1}.

    Let A be an infinitesimal generator of a uniformly bounded family of strongly continuous cosine operators {R(t)}tR which is defined on a Banach space X. We collect some basic properties of a cosine family and its relations with the operator A and the associated sine family.

    Definition 2.4. [24] Consider {R(t)}tR is a one parameter family of bounded linear operators mapping the Banach space XX. It is referred to a strongly continuous cosine family if and only if

    (ⅰ) R(0)=I;

    (ⅱ) R(s+t)+R(st)=2R(s)R(t) for all s,tR;

    (ⅲ) The function tR(t)x is a continuous on R for any xX.

    The sine family {T(t)}tR is correlated to the strongly continuous cosine family {R(t)}tR, it is characterized by

    T(t)x=t0R(s)xds,xX,tR.

    Lemma 2.2. [24] Consider A is an infinitesimal generator of a strongly continuous cosine family {R(t)}tR on Banach space X such that R(t)Meξ|t|,tR. Then, for λ>ξ and (ξ2,)ρ(A) (the resolvent set of A), we have

    λR(λ2;A)x=0eλtR(t)xdt,R(λ2;A)x=0eλtT(t)xdt,xX,

    where the operator R(λ;A)=(λIA)1 is the resolvent of the operator A and λρ(A).

    In this case, the operator A is defined by

    Ax=limt0d2dt2R(t)x,xD(A),

    where D(A)={xX:R(t)xC2(R,X)} is the domain of the operator A. Clearly the infinitesimal generator A is densely defined operator in X and closed.

    In the sequel to present our results, we need the following:

    Definition 2.5. Suppose that τ>0, the Mainardi's Wright-type function is defined as

    Mϱ(τ)=n=0(τ)nn!Γ(1ϱ(n+1)),ϱ(0,1),τC,

    and achieves

    Mϱ(τ)0,0θξMϱ(θ)dθ=Γ(1+ξ)Γ(1+ϱξ),ξ>1.

    The abstract phase space Ph is demonstrated by convenient way [25,26]. Let h=C((,0],[0,)) with 0h(t)dt<. Then, for any c>0, we can define the set

    P={A:[c,0]X,A is bounded and measurable},

    and establish the space P with the norm

    AP=sups[c,0]|A(s)|,for allAP.

    Let us define the space

    Ph={A:(,0]Xsuch that for any c>0,A|[c,0]Pand0h(t)supts0A(s)dt<}.

    If Ph is furnished with the norm

    APh=0h(t)supts0A(s)dt,APh,

    then (Ph,Ph) is a Banach space. Next, we introduce the available space

    ¯Ph={v:(,a]Xsuch thatv|[0,a]C((tk,tk+1],X),v|(,0]=ϕPh},

    which has the norm

    x¯Ph=sups[0,a]v(s)+ϕPh,x¯Ph.

    Definition 2.6. [27] If v:(,a]X,a>0, such that ϕPh. The situations listed below are accurate for all τ[0,a],

    1) vτPh;

    2) Two functions, ζ1(τ),ζ2(τ)>0, are such that ζ1(τ):[0,)[0,) is a continuous function and ζ2(τ):[0,)[0,) is a locally bounded function which are independent to v() whereas

    vτPhζ1(τ)sup0<s<τv(s)+ζ2(τ)ϕPh;

    3) v(τ)HvτPh, where H>0 is a constant.

    Before introducing the mild solution of evolution Eq (1.1), we have to establish the following Lemmas.

    Lemma 3.1. Let Iαs be the left R-L integral of order α and f(t) is integrable function defined for ts0. Then,

    seλtIαsf(t)dt=λαseλtf(t)dt.

    Proof. From the Definition 2.2 and the rule of converting double integral to single integral, we get

    seλtIαsf(t)dt=seλtts(tτ)α1f(τ)dτdt=1Γ(α)sf(τ)dττeλt(tτ)α1dt=1Γ(α)sf(τ)eλτdτ0tα1eλtdt=λαseλtf(τ)dτ.

    The proof is over.

    Lemma 3.2. Let 1<α2 and h:JX be an integrable function. Then, the mild solution to our problem (1.1) possess the form

    u(t)={ϕ(t),t(,0],Rq(t)ϕ(0)+t0Rq(s)uods+t0(ts)q1Tq(t,s)h(s)ds,t[0,t1],μk(t,u(t)),tmk=1(tk,sk],Rq(tsk)μk(sk,u(sk))+tskRq(ysk)ξk(sk,u(sk))dy+tsk(ty)q1Tq(ty)h(y)dy,tmk=1(sk,tk+1],

    where 1/2<q=α21,

    Rq(t)=0Mq(θ)R(tqθ)dθ,Tq(t,s)=q0θMq(θ)T((ts)qθ)dθ,

    and Mq is a probability density function defined by Definition 2.5.

    Proof. Using Lemma 2.1 with operating by Iαsr on both sides to the fractional differential equation in (1.1), we arrive at

    u(t)=Iαsk[Au(t)+h(t)]+c1,k(tsk)+c0,k, (3.1)

    where c1,k,c0,kR,k=0,1,,m are constants to be determined.

    ● For t[0,t1]: By taking ρ1 to the results given in Lemma 5 in [28], we have

    u(t)=Rq(t)ϕ(0)+t0Rq(s)uods+t0(ts)q1Tq(t,s)h(s)ds.

    ● For t(t1,s1]: We obtain

    u(t)=μ1(t,u(t)) andu(t)=ξ1(t,u(t)).

    ● For t(s1,t2]: The problem (1.1) becomes

    cDαs1u(t)=Au(t)+h(t),u(s1)=μ1(s1,u(s1)),u(s1)=ξ1(s1,u(s1)).

    In this interval, Eq (3.1) becomes

    u(t)=Iαs1[Au(t)+h(t)]+c1,1(ts1)+c0,1.

    Considering the past impulsive conditions, we get

    c0,1=μ1(s1,u(s1))andc1,1=ξ1(s1,u(s1)),

    which imply that

    u(t)=Iαs1[Au(t)+h(t)]+ξ1(s1,u(s1))(ts1)+μ1(s1,u(s1)).

    Multiplying both sides by eλt followed by integrating from s1 to , we achieve

    U(λ)=λα{AU(λ)+H(λ)}+λ1eλs1μ1(s1,u(s1))+λ2eλs1ξ1(s1,u(s1)),

    where

    U(λ)=s1u(t)eλtdtandH(λ)=s1h(t)eλtdt.

    Given that (λαIA)1 exists, then λαρ(A). We obtain

    U(λ)=(λαIA)1{λα1eλs1μ1(s1,u(s1))+λα2eλs1ξ1(s1,u(s1))+H(λ)}=λq1eλs10eλqtR(t)μ1(s1,u(s1))dt+λq2eλs10eλqtR(t)ξ1(s1,u(s1))dt+0eλqtT(t)H(λ)dt.

    Let Ψq(θ)=qθq+1Mq(θq) be defined for θ(0,) and q(12,1). Then,

    0epθΨq(θ)dθ=epq,

    which can be used to calculate the first term with replacing t by sq as

    λq1eλs10eλqtR(t)μ1(s1,u(s1))dt=q0(λs)q1e(λs)qR(sq)eλs1(μ1s1,u(s1))ds=1λ0dds(e(λs)q)R(sq)eλs1(μ1s1,u(s1))ds=00θΨq(θ)eλsθR(sq)eλs1(μ1s1,u(s1))dθds=0eλ(x+s1){0Ψq(θ)R((xθ)q)μ1(s1,u(s1))dθ}dx=0eλ(x+s1){0Mq(θ)R(xqθ)μ1(s1,u(s1))dθ}dx=0eλ(x+s1)Rq(x)μ1(s1,u(s1))dx=s1eλtRq(ts1)μ1(s1,u(s1))dt.

    By using Lemma 3.1 with α=1, we get

    λq20eλqtRq(t)eλs1ξ1(s1,u(s1))dt=s1eλt{ts1Rq(ys1)ξ1(s1,u(s1))dy}dt.

    Finally, we can write

    0eλqtT(t)H(λ)dt=q0e(λs)qT(sq)sq1H(λ)ds=q00eλsθΨq(θ)T(sq)sq1H(λ)dθds=q0s10θqeλxΨq(θ)T((xθ)q)xq1eλyh(y)dθdydx=q0s10eλ(x+y)θMq(θ)T(xqθ)xq1h(y)dθdydx=qs1y0eλtθMq(θ)T((ty)qθ)(ty)q1h(y)dθdtdy=s1yeλt(ty)q1Tq(ty)h(y)dtdy=s1eλt{ts1(ty)q1Tq(ty)h(y)dy}dt.

    In conclusion, we can write

    s1eλtu(t)dt=s1eλt{Rq(ts1)μ1(s1,u(s1))+ts1Rq(ys1)ξ1(s1,u(s1))dy+ts1(ty)q1Tq(ty)h(y)dy}dt.

    Therefore, by taking the inverse Laplace transform, we have

    u(t)=Rq(ts1)μ1(s1,u(s1))+ts1Rq(ys1)ξ1(s1,u(s1))dy+ts1(ty)q1Tq(ty)h(y)dy.

    ● For t(sk,tk+1],k=2,3,,m: In a similar manner, we can write

    u(t)=Rq(tsk)μk(sk,u(sk))+tskRq(ysk)ξk(sk,u(sk))dy+tsk(ty)q1Tq(ty)h(y)dy.

    Consequently, we get the solution from the earlier (1.1). Direct calculations show that the opposite results are true. The proof is completed.

    Remark 3.1. [28] From linearity of R(t) and T(t) for all t0, it is clearly to deduce that Rq(t) and Tq(t,s) are also linear operators where 0<s<t. Therefore, the proofs of all next Lemmas are same when taking ρ approaches 1.

    Lemma 3.3. [28] The following estimates for Rq(t) and Tq(t,s) are verified for any fixed t0 and 0<s<t

    |Rq(t)x|ϖ|x|and|Tq(t,s)x|ϖaqΓ(2q)|x|.

    Lemma 3.4. [28] The operators Rq(t) and Tq(s,t) are strongly continuous for every 0<s<t and t>0.

    Lemma 3.5. [28] Assume that R(t) and T(t,s) are compact for every 0<s<t. Then, the operators Rq(t) and Tq(s,t) are compact for every 0<s<t.

    Define the operator N:¯Ph¯Ph as follows

    N(u)(t)={ϕ(t),t(,0],Rq(t)ϕ(0)+t0Rq(y)uody+t0(ty)q1Tq(t,y)h(y,uy)dy,t[0,t1],μk(t,u(t)),tmk=1(tk,sk],Rq(tsk)μk(sk,u(sk))+tskRq(ysk)ξk(sk,u(sk))dy,+tsk(ty)q1Tq(ty)h(y,uy)dy,tmk=1(sk,tk+1].

    Let ϰ():(,a]X be the function denoted by

    ϰ(t)={ϕ(t),t(,0],0,t(0,a].

    Plainly, ϰ(0)=ϕ(0). For each zC([0,a],X) with z(0)=0, we indicate by ϑ to the function defined as

    ϑ(t)={0,t(,0],z(t),t[0,a].

    If u() satisfies that u(t)=N(u)(t) for all t(,a], we can decompose that u(t)=ϑ(t)+ϰ(t), t(,a], it denotes ut=ϑt+ϰt for every t(,a] and the function z() satisfies

    z(t)={Rq(t)ϕ(0)+t0Rq(y)uody+t0(ty)q1Tq(t,y)h(y,ϑy+ϰy)dy,t[0,t1],μk(t,ϑ+x),tmk=1(tk,sk],Rq(tsk)μk(sk,ϑ+ϰ)+tskRq(ysk)ξk(sk,ϑ+ϰ)dy,+tsk(ty)q1Tq(ty)h(ϑy+ϰy)dy,tmk=1(sk,tk+1].

    Set the space Υ={zC([0,a],X),z(0)=0} equipped the norm

    zΥ=supt[0,a]z(t).

    Therefore, (Υ,Υ) is a Banach space. Assume that the operator G:ΥΥ is formulated as follows:

    G(z)(t)={Rq(t)ϕ(0)+t0Rq(y)uody+t0(ty)q1Tq(t,y)h(y,ϑy+ϰy)dy,t[0,t1],μk(t,ϑ+x),tmk=1(tk,sk],Rq(tsk)μk(sk,ϑ+ϰ)+tskRq(ysk)ξk(sk,ϑ+ϰ)dy,+tsk(ty)q1Tq(ty)h(ϑy+ϰy)dy,tmk=1(sk,tk+1].

    The operator N seems to have a fixed point is equivalent to G has a fixed point. Thus, we proceed to prove that G has a fixed point.

    Now, we make the following assumptions:

    (E1) The function h:[0,a]×PhX is a continuous and μk,ξk:[tk,sk]×XX are continuous functions for all k=1,2,,m;mN.

    (E2) There is a constant Ω>0 satisfying

    h(t,ut)h(t,vt)ΩutvtPh.

    (E3) There exist δk,δk>0;k=1,2,,m;mN such that

    μk(t,u)δkandξk(t,u)δk.

    (E4) There are positive constants Dk,Dk,k=1,2,,m;mN such that

    μk(t,u1)μk(t,u2)Dku1u2,
    ξk(t,u1)ξk(t,u2)Dku1u2.

    (E5) There exists a continuous function g(t):[0,a][0,) such that, for any (t,ut)[0,a]×Ph, it satisfies

    h(t,ut)g(t)utPh.

    The brief constants that will be utilized later to streamline handling, are listed as follow

    E(q)=tq+11qΓ(2q),Ek(q)=a(ask)qqΓ(2q),B=ϖΩζ1,¯B=ϖgζ1E(q),Bk=ϖHζ1[Dk+Dk(ask)],¯Bk=ϖgζ1Ek(q),O=ϖ(ϕ(0)+uot1+E(q)Ωζ2ϕPh+E(q)c),Ok=ϖ(δk+δk(ask)+Ek(q)Ωζ2ϕPh+Ek(q)c),Q=ϖ(ϕ(0)+uot1+E(q)ζ2gϕPh),Qk=ϖ(δk+δk(ask)+Ek(q)ζ2gϕPh)

    where k=1,2,,m;mN.

    Lemma 4.1. Assume that the requirement (E2) is met by c=maxt[0,a]|h(t,0)|. Ponder about the expressions ζ1=supt[0,a]ζ1(t) and ζ2=supt[0,a]ζ2(t) where ζ1() and ζ2() are established in Definition 2.6. Then,

    h(t,ϑt+ϰt)Ω(ζ1zΥ+ζ2ϕPh)+c.

    Proof. Regarding Definition 2.6 and the presumption (E2). Then,

    h(t,ϑt+ϰt)=h(t,ϑt+ϰt)h(t,0)+h(t,0)h(t,ϑt+ϰt)h(t,0)+h(t,0)Ωϑt+ϰtΥ+cΩ(ζ1(t)supt[0,a]ϑ(t)+ζ2(t)ϕPh)+cΩ(ζ1(t)zΥ+ζ2(t)ϕPh)+cΩ(ζ1zΥ+ζ2ϕPh)+c.

    This ends the proof.

    Lemma 4.2. Suppose that the statement (E5) is satisfied with g=supt[0,a]g(t). Let ζ1=supt[0,a]ζ1(t) and ζ2=supt[0,a]ζ2(t) where ζ1() and ζ2() are outlined in Definition 6. Then,

    h(t,ϑt+ϰt)(t),

    where

    =supt[0,a](t)=supt[0,a]{g(t)(ζ1(t)zΥ+ζ2(t)ϕPh)}=g(ζ1zΥ+ζ2ϕPh).

    Proof. By the same way in Lemma 4.1, we can easily reach the desired result.

    Theorem 4.1. Consider the assertions (E1)(E4) hold and

    Λ=maxk{BE(q),DkHζ1,Bk+BEk(q)}.

    Then, the fractional evolution equation with non-instantaneous impulsive (1.1) has a unique mild solution on (,a] if Λ<1.

    Proof. To show that the operator G maps bounded subset of Υ into bounded subset in Υ, we set

    Υr={zΥ:zΥr},

    where

    rmaxk{O1BE(q),δk,Ok1BEk(q)}.

    Then, for any zΥr and in spite of (E2) and (E3) and Lemma 4.1. Correspondingly, three situations are taken into consideration.

    Case Ⅰ. Whenever t[0,t1], we have

    G(z)(t)Υϖ(ϕ(0)+uot1)+ϖt1Γ(2q)t0(ts)q1h(y,ϑy+ϰy)dyϖ[ϕ(0)+uot1+tqt1qΓ(2q){Ω(ζ1zΥ+ζ2ϕPh)+c}]ϖ[ϕ(0)+uot1+E(q){Ω(ζ1zΥ+ζ2ϕPh)+c}]O+E(q)BzΥO+E(q)Brr.

    Case Ⅱ. Whenever t(tk,sk],k=1,,m;mN, we have

    G(z)(t)Υ=μk(t,ϑ+ϰ)δk.

    Case Ⅲ. Whenever t(sk,sk+1],k=1,,m;mN, we have

    G(z)(t)Υϖ[δk+δk(tsk)+aΓ(2q)tsk(ty)q1h(y,ϑy+ϰy)dy]ϖ[δk+δk(ask)+a(ask)qqΓ(2q){Ω(ζ1zΥ+ζ2ϕPh)+c}]ϖ[δk+δk(ask)+Ek(q){Ω(ζ1zΥ+ζ2ϕPh)+c}]Ok+Ek(q)Brr.

    For the aforementioned, we acquire G(z)(t)Υr. Thus, the operator G maps bounded subset into bounded subset in Υ.

    Now, we prove that the operator G is a contraction mapping. Certainly, consider z,zΥ. Then, there still are the subsequent situations.

    Case Ⅰ. For any t[0,t1], we have

    G(z)(t)G(z)(t)ϖt1Γ(2q)t0(ts)q1h(y,ϑy+ϰy)h(y,ϑy+ϰy)dyΩϖt1Γ(2q)t0(ts)q1ϑyϑyPhdyΩϖt1Γ(2q)ζ1zzΥt0(ts)q1dyΩϖtq+11qΓ(2q)ζ1zzΥ=BE(q)ϑϑΥ.

    Case Ⅱ. For any t(tk,sk],k=1,,m;mN, we have

    G(z)(t)G(z)(t)=μk(t,ϑ+ϰ)μk(t,ϑ+ϰ)DkϑϑΥDkHztztPhDkHζ1zzΥ=DkHζ1ϑϑΥ.

    Case Ⅲ. For any t(sk,sk+1],k=1,,m;mN, we have

    G(z)(t)G(z)(t)Rq(tsk)μk(sk,ϑ+ϰ)μk(sk,ϑ+ϰ)+tskRq(ysk)ξk(sk,ϑ+ϰ)ξk(sk,ϑ+ϰ)dy+tsk(ty)q1Tq(ty)h(ϑy+ϰy)h(ϑy+ϰy)dyϖDkϑϑΥ+ϖDktskϑ(y)ϑ(y)Υdy+ϖaΓ(2q)Ωtsk(ty)q1ϑyϑyPhdyϖDkHzyzyPh+ϖDkHtskzyzyPhdy+ϖaΓ(2q)Ωtsk(ty)q1ϑyϑyPhdyϖHζ1[Dk+DkEk(q)]zzΥ+ϖaΓ(2q)Ωζ1tsk(ty)q1zzΥdyϖHζ1[Dk+Dk(ask)+a(ask)qΩ qΓ(2q)H]zzΥ=(Bk+BEk(q))zzΥ=(Bk+BEk(q))ϑϑΥ.

    For the aforementioned, we may write

    G(z)(t)G(z)(t)ΥΛϑϑΥ.

    Amid the existing circumstances, Λ<1 shows that the operator G is a contraction. This suggests that the problem (1.1) has a unique solution on (,a] relying on the Banach contraction mapping principle.

    Remark 4.1. In viewing our problem, it is very difficult to obtain the exact solution and so it is useful to investigate some properties of the solutions, especially the uniqueness. The previous theorem shaw that the mild solution of the problem (1.1) is unique under the assumptions (E1)(E4) and Λ<1. This enables us to apply our results to real-life problem or phenomena as in the last section.

    Assume that the operator G is divided as a sum of the two operators Gi:ΥΥ,i=1,2 as

    G=G1(z)+G2(z) (4.1)

    where,

    G1(z)(t)={Rq(t)ϕ(0)+t0Rq(y)uody+t0(ty)q1Tq(t,y)h(y,ϑy+ϰy)dy,t[0,t1],0,tmk=1(tk,sk],tsk(ty)q1Tq(ty)h(ϑy+ϰy)dy,tmk=1(sk,tk+1]

    and

    G2(z)(t)={0,t[0,t1],μk(t,ϑ+x),tmk=1(tk,sk],Rq(tsk)μk(sk,ϑ+ϰ)+tskRq(ysk)ξk(sk,ϑ+ϰ)dy,tmk=1(sk,tk+1].

    Theorem 4.2. Suppose the hypotheses (E1)and(E3)(E5) are correct. Then the fractional evolution equation with non-instantaneous impulsive (1.1) has at least one mild solution on (,a] if Δ<1 where Δ is given by

    Δ=maxk{¯B,¯Bk}.

    Proof. Let the operators G1 and G2 be defined as (4.1). Setting g=supt[0,a]|g(t)|. Let us define the closed ball Υρ={zΥ:zΥρ} with radius

    ρmaxk{Q1¯B,δk,Qk1¯Bk}.

    Then, for u,vΥρ, we claim that G1(z)(u)+G2(z)(v)ρ which concludes that G1(u)+G1(v)Υρ. To verify our claiming, we show that G maps bounded sets of Υ into bounded sets in Υ, for any ρ0. Then for any zΥρ and in light of (E3),(E5) and Lemma 4.2, we have three cases

    Case Ⅰ. For any t[0,t1], we have

    G(z)(t)Υϖ(ϕ(0)+uot1)+ϖt1Γ(2q)t0(ts)q1h(y,ϑy+ϰy)dyϖ[ϕ(0)+uot1+t1Γ(2q)t0(ts)q1{g(y)(ζ1(y)zΥ+ζ2(y)ϕPh)}dy]ϖ[ϕ(0)+uot1+E(q){g(ζ1zΥ+ζ2ϕPh)}]Q+E(q)¯BzΥQ+¯Bρρ.

    Case Ⅱ. For any t(tk,sk],k=1,,m;mN, we have

    G(z)(t)Υ=μk(t,ϑ+ϰ)δkρ.

    Case Ⅲ. For any t(sk,sk+1],k=1,,m;mN, we have

    G(z)(t)Υϖ[δk+δk(tsk)+aΓ(2q)tsk(ty)q1h(y,ϑy+ϰy)dy]ϖ[δk+δk(ask)+aΓ(2q)tsk(ty)q1{g(y)(ζ1(y)zΥ+ζ2(y)ϕPh)}dy]ϖ[δk+δk(ask)+a(ask)qqΓ(2q){g(ζ1zΥ+ζ2ϕPh)}]ϖ[δk+δk(ask)+Ek(q){g(ζ1zΥ+ζ2ϕPh)}]Qk+¯Bkρρ.

    By virtue of the above, we obtain G(z)(t)Υρ. Thus the operator G maps bounded sets into bounded sets in Υ.

    The following step is to confirm that the operator G2 maps bounded sets into equicontinuous sets in Υ. In light of the situation (E1), G2 is continuous. The following scenarios are therefore possible.

    Case Ⅰ. For each tkγ1<γ2sk and zΥρ, we have

    G2(z)(γ2)G2(z)(γ1)μk(γ2,ϑ+ϰ)μk(γ1,ϑ+ϰ).

    Due to the continuity of μ(t,u(t)). It is clear that the above inequality approaches zero when letting γ2γ1.

    Case Ⅱ. For any skγ1<γ2tk+1,k=1,,m;mN and zΥρ, we have

    G2(z)(γ2)G2(z)(γ1)δkRq(γ2sk)Rq(γ1sk)+δkγ2γ1Rq(ysk)dyδkRq(γ2sk)Rq(γ1sk)+δkϖ(γ2γ1).

    Due to compactness of operator Rq(y) and Tq(t,y) (see Lemma 3.5), we infer that G2(z)(γ1)G2(z)(γ2)0 as γ2γ1. Thus, G2 is a relatively compact on Υ. By Arezela Ascoli Theorem the operator G2 is completely continuous on Υρ. The only thing left to do is provide evidence that G1 is a contraction mapping. Thus, two cases are thought about.

    Case Ⅰ. For any t[0,t1],k=1,,m;mN and zΥρ, we have

    G1(z)(t)G1(z)(t)ϖt1Γ(2q)t0(ty)q1h(y,ϑy+ϰy)h(y,ϑy+ϰy)dyϖt1Γ(2q)t0(ty)q1g(y)ϑyϑyPhdygϖt1Γ(2q)ζ1zzΥt0(ts)q1dygϖtq+11qΓ(2q)ζ1zzΥ=¯BϑϑΥ.

    Case Ⅱ. For any t(sk,tk+1],k=1,,m;mN, we have

    G1(z)(t)G1(z)(t)tsk(ty)q1Tq(ty)h(ϑy+ϰy)h(ϑy+ϰy)dyϖaΓ(2q)tsk(ty)q1g(y)ϑyϑyPhdyaΓ(2q)ϖgζ1zzΥtsk(ty)q1dya(ask)q qΓ(2q)ϖgζ1zzΥ=¯BkzzΥ.

    As a sense, the fractional evolution equation with non-instantaneous impulsive (1.1) has at least one mild solution on Υ, according to the Krasnoselskii Theorem. The evidence is now complete.

    Remark 4.2. Also, it is useful to investigate the existence of the solution instead of the its uniqueness. The theorem above shaw that the mild solution of the problem (1.1) exists under the assumptions (E1)(E3) and (E5) with Δ<1.

    Presume the following fractional wave equation with impulsive effect and infinite delay

     {u(t,x)=13sint,t(,0],x[0,π],cDαtu(t,x)=2x2u(t,x)+h(t,ut),t(0,25](45,1],x[0,π],u(t,x)=17t32+15sinu(t),t(25,45]x[0,π],u(t)=314t12+18cosu(t),t(25,45]x[0,π],u(t,0)=u(t,π)=0,t[0,1],u(0,x)=32ex3,x[0,π].

    Consider that

    J=[0,1],0=t0=s0<t1=25<s1=45<t2<1=a,α=32q=34,u0=32ex3,A=2x2,x[0,π],H=116.Whileζ1(t)=35t32ζ1=35{supt(0,1]t32}ζ1(1)=35.If we takeϖ=1Tq(t,s)1Γ(34),0<s<t1.

    Case Ⅰ. Banach fixed point theorem.

    In order to explain Theorem 4.1, we obtain:

    h(t,ut)=t38t+1+19ut. (5.1)

    Clearly, h:[0,1]×PhR is continuous and satisfying, for ut,vtPh, that

    h(t,ut)h(t,vt)19utvtPh,

    it suggests that Ω=19. For all t(25,45] and u,vR, we get

    μ(t,u)μ(t,v)15sinusinv15uv,ξ(t,u)ξ(t,v)18cosucosv18uv.

    As you can see, the Theorem 4.1 condition (E4) is satisfied with

    Dk=15andDk=18.

    In summary, we have

    Λ=maxk{BE(q),DkHζ1,Bk+BEk(q)}={0.0202,0.0075,0.0384}=0.0384<1.

    Thus all assumptions of this theorem are verified. Therefore, the problem (1.1) has a unique mild solution on (,1].

    Case Ⅱ. Krasnoselskii's theorem.

    To realize Theorem 4.2, take h(t,ut) as given in (5.1). Therefore, g(t)=t38t+1 is increasing function which admits the hypothesis (E5) with

    gg(1)=182.

    These calculate that

    Δ=maxk{¯B,¯Bk}=maxk{0.0161,0.0239}=0.0239<1.

    Since every requirements of Theorem 4.2 are met, it follows that there exists at least one mild solution of (1.1) on (,1].

    We analyzed a set of impulsive fractional evolution equations with infinite delay in the current work. Current functional analysis methodologies serve as the foundation for our conclusions. By using the unbounded operator A as the generator of the strongly continuous cosine family, we were able to suggest a mild solution for the suggested problem. In the instance of problem (1.1), we had two successful outcomes: While the second argument focuses on whether there are solutions for the given problem, the first argument concentrated on the existence and uniqueness of the solution.

    The first result, which is built on a Banach fixed point theorem, provides criteria for ensuring that the problem at hand has no prior solutions by requiring the usage of h(t,ut) to satisfy the classic Lipschitz condition.

    The second argument was based on a Krasnoselskii's theorem, which allows h(t,ut) to behave as h(t,ut)g(t)utPh. The instruments used by fixed point theory in the scenario with simple assumptions. Finally, a numerical example that examines a function that satisfies all the prerequisites was provided to illustrate our conclusion.

    In the next paper, we will study the controllability of mild solution to fractional evolution equations with an infinite time-delay and nonlocal condition by applying Krasnoselskii's theorem in the compactness case and the Sadvskii and Kuratowski measure of noncompactness.

    The Deanship of Scientific Research (DSR) at King Abdulaziz University (KAU), Jeddah, Saudi Arabia has funded this project, under grant no. (KEP-PhD: 34-130-1443).

    The authors declare that they have no conflicts of interest.



    [1] D. Baleanu, A. Jajarmi, H. Mohammadi, S. Rezapour, A new study on the mathematical modelling of human liver with Caputo-Fabrizio fractional derivative, Chaos Soliton. Fract., 134 (2020), 109705. https://doi.org/10.1016/j.chaos.2020.109705 doi: 10.1016/j.chaos.2020.109705
    [2] B. Ghanabri, S. Djilali, Mathematical analysis of a fractional-order predator-prey model with prey social behavior and infection developed in predator population, Chaos Soliton. Fract., 138 (2020), 109960. https://doi.org/10.1016/j.chaos.2020.109960 doi: 10.1016/j.chaos.2020.109960
    [3] A. Salem, Existence results of solutions for ant-periodic fractional Langevin equation, J. Appl. Anal. Comput., 10 (2020), 2557–2574. https://doi.org/10.11948/20190419 doi: 10.11948/20190419
    [4] A. Salem, N. Mshary, On the existence and uniqueness of solution to fractional-order Langevin equation, Adv. Math. Phys., 2020 (2020), 8890575. https://doi.org/10.1155/2020/8890575 doi: 10.1155/2020/8890575
    [5] T. V. An, N. D. Phu, N. V. Hoa, A survey on non-instantaneous impulsive fuzzy differential equations involving the generalized Caputo fractional derivative in the short memory case, Fuzzy Set. Syst., 443 (2022), 160–197. https://doi.org/10.1016/j.fss.2021.10.008 doi: 10.1016/j.fss.2021.10.008
    [6] A. Salem, R. Babusail, Finite-time stability in nonhomogeneous delay differential equations of fractional Hilfer type, Mathematics, 10 (2022), 1520. https://doi.org/10.3390/math10091520 doi: 10.3390/math10091520
    [7] A. Salem, A. Al-dosari, Positive solvability for conjugate fractional differential inclusion of (k,nk) type without continuity and compactness, Axioms, 10 (2021), 170. https://doi.org/10.3390/axioms10030170 doi: 10.3390/axioms10030170
    [8] A. Salem, L. Almaghamsi, F. Alzahrani, An infinite system of fractional order with p-Laplacian operator in a tempered sequence space via measure of noncompactness technique, Fractal Fract., 5 (2021), 182. https://doi.org/10.3390/fractalfract5040182 doi: 10.3390/fractalfract5040182
    [9] R. Agarwal, R. Almeida, S. Hristova, D. O'Regan, Non-instantaneous impulsive fractional differential equations with state dependent delay and practical stability, Acta Math. Sci., 41 (2021), 1699–1718. https://doi.org/10.1007/s10473-021-0518-1 doi: 10.1007/s10473-021-0518-1
    [10] A. Salem, S. Abdullah, Non-instantaneous impulsive BVPs involving generalized Liouville-Caputo derivative, Mathematics, 10 (2022), 291. https://doi.org/10.3390/math10030291 doi: 10.3390/math10030291
    [11] S. Asawasamrit, Y. Thadang, S. K. Ntouyas, J. Tariboon, Non-instantaneous impulsive boundary value problems containing Caputo fractional derivative of a function with respect to another function and Riemann-Stieltjes fractional integral boundary conditions, Axioms, 10 (2021), 130. https://doi.org/10.3390/axioms10030130 doi: 10.3390/axioms10030130
    [12] A. Salem, A. Al-Dosari, Hybrid differential inclusion involving two multi-valued operators with nonlocal multi-valued integral condition, Fractal Fract., 6 (2022), 109. https://doi.org/10.3390/fractalfract6020109 doi: 10.3390/fractalfract6020109
    [13] S. Zhaoa, M. Song, Stochastic impulsive fractional differential evolution equations with infinite delay, Filomat, 31 (2017), 4261–4274. https://doi.org/10.2298/FIL1713261Z doi: 10.2298/FIL1713261Z
    [14] A. Salem, B. Alghamdi, Multi-strip and multi-point boundary conditions for fractional Langevin equation, Fractal Fract., 4 (2020), 18. https://doi.org/10.3390/fractalfract4020018 doi: 10.3390/fractalfract4020018
    [15] A. Salem, N. Mshary, Coupled fixed point theorem for the generalized Langevin equation with four-point and Strip conditions, Adv. Math. Phys., 2022 (2022), 1724221. https://doi.org/10.1155/2022/1724221 doi: 10.1155/2022/1724221
    [16] R. Saadati, E. Pourhadi, B. Samet, On the PC-mild solutions of abstract fractional evolution equations with non-instantaneous impulses via the measure of noncompactness, Bound. Value Probl., 2019 (2019), 19. https://doi.org/10.1186/s13661-019-1137-9 doi: 10.1186/s13661-019-1137-9
    [17] A. El. Mfadel, S. Melliani, A. Kassidi, M. Elomari, Existence of mild solutions for nonlocal Ψ-Caputo-type fractional evolution equations with nondense domain, Nonautonomous Dyn. Syst., 9 (2022), 272–289. https://doi.org/10.1515/msds-2022-0157 doi: 10.1515/msds-2022-0157
    [18] A. Salem, L. Almaghamsi, Existence solution for coupled system of Langevin fractional differential equations of caputo type with Riemann-Stieltjes integral boundary conditions, Symmetry, 13 (2021), 2123. https://doi.org/10.3390/sym13112123 doi: 10.3390/sym13112123
    [19] A. Kumar, D. N. Pandey, Existence of mild solution of Atangana-Baleanu fractional differential equations with non-instantaneous impulses and with non-local conditions, Chaos Soliton. Fract., 132 (2020), 109551. https://doi.org/10.1016/j.chaos.2019.109551 doi: 10.1016/j.chaos.2019.109551
    [20] V. Kavitha, M. M. Arjunan, D. Baleanu, Non-instantaneous impulsive fractional-order delay differential systems with Mittag-Leffler kernel, AIMS Math., 7 (2022) 9353–9372. https://doi.org/10.3934/math.2022519 doi: 10.3934/math.2022519
    [21] M. Sova, Cosine operator functions, Warszawa: Instytut Matematyczny Polskiej Akademi Nauk, 1966.
    [22] I. Podlubny, Fractional differential equations, 1999.
    [23] A. A. Kilbas, H. M. Srivastava, J. J. Trujillo, Theory and applications of fractional differential equations, Elsevier, 2006. https://doi.org/10.1016/s0304-0208(06)x8001-5
    [24] C. C. Travis, G. F. Webb, Cosine families and abstractnonlinear second order differential equations, Acta Math. Acad. Sci. Hung., 32 (1978), 75–96. https://doi.org/10.1007/BF01902205 doi: 10.1007/BF01902205
    [25] X. Zhang, X. Huanga, Z. Liu, The existence and uniqueness of mild solutions for impulsive fractional equations with nonlocal conditions and infinite delay, Nonlinear Anal. Hybri., 4 (2010), 775–781. https://doi.org/10.1016/j.nahs.2010.05.007 doi: 10.1016/j.nahs.2010.05.007
    [26] M. Benchohra, J. Henderson, S. K. Ntouyas, A. Ouahaba, Existence results for fractional order functional differential equations with infinite delay, J. Math. Anal. Appl., 338 (2008), 1340–1350. https://doi.org/10.1016/j.jmaa.2007.06.021 doi: 10.1016/j.jmaa.2007.06.021
    [27] J. Hale, J. Kato, Phase space for retarded equations with infinite delay, Funkcial. Ekvac., 1978.
    [28] A. Salem, K. N. Alharbi, H. M. Alshehri, Fractional evolution equations with infinite time delay in abstract phase space, Mathematics, 10 (2022), 1332. https://doi.org/10.3390/math10081332 doi: 10.3390/math10081332
  • This article has been cited by:

    1. Kholoud N. Alharbi, Hyers–Ulam and Hyers–Ulam–Rassias Stability for a Class of Fractional Evolution Differential Equations with Neutral Time Delay, 2025, 17, 2073-8994, 83, 10.3390/sym17010083
    2. Ahmed Salem, Rania Al-Maalwi, Fractional Evolution Equation with Nonlocal Multi-Point Condition: Application to Fractional Ginzburg–Landau Equation, 2025, 14, 2075-1680, 205, 10.3390/axioms14030205
  • 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(1532) PDF downloads(93) Cited by(2)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog