In order to investigate the asphalt mixture enhancement due to the granular anti-rutting additives, the multiple stress creep recovery (MSCR) test and wheel tracking test were performed to examine the effects of three anti-rutting additives on the high-temperature performance of binders and mixtures, respectively. Further, the interaction mechanism between anti-rutting additives and asphalt binder was revealed via a fluorescence microscopy (FM) test. The results indicate that the incorporation of anti-rutting additives causes a substantial increase in elasticity behavior for binders, along with a massive enhancement of dynamic stability for mixtures. Further, the enhancement of performance is not conclusively determined by the compatibility between the additive and asphalt. Instead, unmelted viscous-flow additives contribute to the enhancement of mixtures by interlocking, filling voids, cementing and wrapping the aggregate, which are essential mechanisms of mixture modification. This study contributes to selecting appropriate additives in engineering and enhancing anti-rutting additives based on their modification mechanism.
Citation: Gang Tian, Ying Gao, Conglin Chen, Qin Ye. Study on high-temperature performance evaluation of asphalt mixtures with anti-rutting additives[J]. Electronic Research Archive, 2024, 32(1): 160-173. doi: 10.3934/era.2024008
[1] | 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 |
[2] | Weerawat Sudsutad, Chatthai Thaiprayoon, Sotiris K. Ntouyas . Existence and stability results for $ \psi $-Hilfer fractional integro-differential equation with mixed nonlocal boundary conditions. AIMS Mathematics, 2021, 6(4): 4119-4141. doi: 10.3934/math.2021244 |
[3] | 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 |
[4] | 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 |
[5] | Xiaoming Wang, Rizwan Rizwan, Jung Rey Lee, Akbar Zada, Syed Omar Shah . Existence, uniqueness and Ulam's stabilities for a class of implicit impulsive Langevin equation with Hilfer fractional derivatives. AIMS Mathematics, 2021, 6(5): 4915-4929. doi: 10.3934/math.2021288 |
[6] | J. Vanterler da C. Sousa, E. Capelas de Oliveira, F. G. Rodrigues . Ulam-Hyers stabilities of fractional functional differential equations. AIMS Mathematics, 2020, 5(2): 1346-1358. doi: 10.3934/math.2020092 |
[7] | Murugesan Manigandan, R. Meganathan, R. Sathiya Shanthi, Mohamed Rhaima . Existence and analysis of Hilfer-Hadamard fractional differential equations in RLC circuit models. AIMS Mathematics, 2024, 9(10): 28741-28764. doi: 10.3934/math.20241394 |
[8] | Thabet Abdeljawad, Sabri T. M. Thabet, Imed Kedim, Miguel Vivas-Cortez . On a new structure of multi-term Hilfer fractional impulsive neutral Levin-Nohel integrodifferential system with variable time delay. AIMS Mathematics, 2024, 9(3): 7372-7395. doi: 10.3934/math.2024357 |
[9] | 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 |
[10] | Ugyen Samdrup Tshering, Ekkarath Thailert, Sotiris K. Ntouyas . Existence and stability results for a coupled system of Hilfer-Hadamard sequential fractional differential equations with multi-point fractional integral boundary conditions. AIMS Mathematics, 2024, 9(9): 25849-25878. doi: 10.3934/math.20241263 |
In order to investigate the asphalt mixture enhancement due to the granular anti-rutting additives, the multiple stress creep recovery (MSCR) test and wheel tracking test were performed to examine the effects of three anti-rutting additives on the high-temperature performance of binders and mixtures, respectively. Further, the interaction mechanism between anti-rutting additives and asphalt binder was revealed via a fluorescence microscopy (FM) test. The results indicate that the incorporation of anti-rutting additives causes a substantial increase in elasticity behavior for binders, along with a massive enhancement of dynamic stability for mixtures. Further, the enhancement of performance is not conclusively determined by the compatibility between the additive and asphalt. Instead, unmelted viscous-flow additives contribute to the enhancement of mixtures by interlocking, filling voids, cementing and wrapping the aggregate, which are essential mechanisms of mixture modification. This study contributes to selecting appropriate additives in engineering and enhancing anti-rutting additives based on their modification mechanism.
A non-instantaneous impulsive differential equation is due to Hernándaz et. al.[1], and is used to describe impulsive action, which stays active on a finite time interval. Hilfer [2] introduced a fractional derivative, which is a generalization for Riemann-Liouville fractional derivative and Caputo fractional derivative. Many works have been appeared studying various models involving fractional differential with instantaneous and non-instantaneous impulses and providing solutions to those models. For example, Saravanakumar et al. [3] analyzed the existence of mild solution of non instantaneous impulsive for Hilfer fractional stochastic differential equations driven by fractional Brownian motion, . Shu et al.[4] presented a right formula of mild solutions to a fractional semilinear evolution equation generated by a sectorial operator, and its order belongs to the intervals (0,1) and (1,2), Wang et al. [5] studied the global attracting solutions to non-instantaneous impulsive differential inclusions containing Hilfer fractional, and Ngo et al.[6] presented a formula of solution for a non-instantaneous impulsive differential equation containing ψ−Hilfer derivative with lower limit of the fractional derivative at zero. For more works on non-instantaneous impulsive differential equations and inclusions, we refer to [7,8,9,10,11,12,13].
Moreover, Ulam problem [14] has been attracted by many researchers. We highlight some recent works on the existence and Hyers-Ulam stability of solutions for fractional differential equations. Guo et al.[15] investigated the existence and Hyers-Ulam stability of mild solution for an impulsive Riemann-Liouville fractional neutral functional stochastic differential equation with infinite delay of order between one and two, Guo et al.[16] proved the existence and Hyers-Ulam stability of the almost periodic solution to fractional differential equations with impulse involving fractional Brownian, Wang et al. [11] presented the generalized Ulam-Hyers stability for a non-instantaneous impulsive differential inclusions containing the Caputo derivative and Vanterler et al. [17] studied, in finite dimensional Banach spaces, the stability of a Volterra integro-differential equation containing ψ− Hilfer derivative in the sense of Ulam-Hyers. More recently, Vanterler et al.[18] investigated, in finite dimensional Banach spaces, the δ−Ulam-Hyers-Rassias stability for a non-instantaneous impulsive fractional differential equation containing ψ−Hilfer derivative, Benchohra et al. [19] established, in finite dimensional spaces, the existence and stability of solutions for an implicit fractional differential equations with Riemann-Liouville fractional derivative, and Kumar et al. [20] studied the existence and stability of solution for a fractional differential equation with non-instantaneous integral impulses. Very recently, Ben Mahlouf et al.[21] given sufficient conditions to guarantee the existence and stability of solutions for generalized nonlinear fractional differential equations of order α∈(1,2), Elsayed et al.[22] established the existence and stability of boundary value problem for differential equation with Hilfer-Katugampola fractional derivative. For more papers on Ulam-Hyers stability of solutions, we refer to [23,24,25,26,27,28,29,30,31].
It is worth noting that, when the considered problem contains non-instantaneous impulses, there are two approaches in the literature, one by keeping the lower limit of the fractional derivative at zero [6,11,17,18], and the other by switching it at the impulsive points [5,10] Motivated by the above cited work, we prove two existence results of solutions, in infinite dimensional Banach spaces, for a non-instantaneous impulsive fractional differential inclusion involving ψ−Hilfer derivative with delay and we switch the lower limit of the fractional derivative at the impulsive points, and then we study the ψ−generalized Ulam-Hyers stability.
Let E be a real Banach space, J=[0,b], b>0, J∗=(0,b], r>0, 0<ϑ<1, 0≤ν≤1, μ=ϑ+ν−ϑν, Ψ:[−r,0]→E a continuous function except a finite number of discontinuity points s≠0 such that Ψ(0)=0, all values x(s+), and x(s−) are finite, ψ∈C1([0,b],R) be increasing, ψ´(ϱ)>0,;∈J, and Dϑ,ν,ψsi+ be the ψ−Hilfer derivative with lower limit at si of orderϑ and type ν. Moreover, 0=s0< ϱ1<s1<ϱ2<⋯<ϱn<sn<ϱn+1=b, and I1−μ,ψs+ix(s+i)=limϱ→s+i I1−μ,ψs+ix(ϱ) , F:J×E→2E−{ϕ} is a multifunction, and gi:[ϱi,si]×E⟶Ei=1,2,…,n. Finally, for any ϱ∈J, τ(ϱ):H→Θ, x(θ)=x(ϱ+θ),θ∈[−r,0]x∈H, where Θ and H will be introduced in the next section. In this paper, we establish existence results of solutions of the following ψ−Hilfer fractional non-instantaneous impulsive differential inclusions with delay:
{Dϑ,ν,ψs+ix(ϱ)∈F(ϱ,τ(ϱ)x)),a.e.ϱ∈∪i=ni=0(si,ϱi+1]x(ϱ)=Ψ(ϱ),ϱ∈[−r,0]I1−μ,ψ0+x(0+)=Ψ(0),x(ϱ+i)=gi(ϱi,x(ϱ−i)),i=1,......n,x(ϱ)=gi(ϱ,x(ϱ−i)),ϱ∈(ϱisi],i=1,......n,I1−μ,ψs+ix(s+i)=gi(si,x(ϱ−i)),i=1,......n. | (1) |
Then, we investigate the ψ-generalized Ulam-Hyers stability of Problem (1). To achieve our aim, we present a relation between the solutions of this problem and the corresponding fractional integral equation (Lemma5).
To make a comparison between the present paper objectives and other relevant recent papers, we refer to the following:
1-Abbas et al.[31] proved the existence of solutions and studied Ulam-Hyers-Rassias stability of problem (1) in the absence of both delay and impulses effect, E=R, ψ(ϱ)=lnϱ,ϱ∈[1,b], F is a single-valued function.
2-Benchohra [19] investigated the existence and Ulam stability for nonlinear implicit differential equations with Riemann-Liouville derivative, which is including in Hilfer derivative
3-Ngo et al.[6] presented a formula of solution for a non-instantaneous impulsive differential equation containing ψ−Hilfer derivative with lower limit of the fractional derivative at zero and in the absence of delay.
4-Vanterler el al.[17] established the existence and stability of solutions for Problem (1) in the absence of delay. and when E=R, the lower limit of the fractional derivative at zero F is a single-valued function
5-Wang et al. [10] considered Problem (1) in the absence of delay, when ψ(ϱ)=ϱ and without studying the stability of solutions.
6-Wang et al.[11] consider a non-instantaneous impulsive semilinear differential inclusions containing Caputo derivative and in the absence of delay.
To clarify the novelty and contribution of this study, we refer to, in this paper, we present a relation between a solution of Problem (1) and the corresponding fractional integral equation (Lemma5), provide two methods to demonstrate the existence of solutions for Problem (1), then, investigate the ψ−generalized Ulam-Hyers stability of solutions. Because our considered problem contains ψ−Hilfer fractional derivative, non-instantaneous impulses with the lower limit of the fractional derivative switches at the impulsive points, presence of delay, and the right hand side is a multi-valued function, therefore, this study generalize recent results, as it is shown above, such as [6,10,11,17,19,31]. In addition, there isn't work in the literature, on ψ−Hilfer fractional non-instantaneous impulsive differential inclusions, in infinite dimensional spaces, in the presence of delay, and the lower limit of the fractional derivative switches at the impulsive points. Moreover, the technique presented in this paper can be used to study the existence and Ulam-Hyers stability of solutions or mild solutions for the problems considered in [3,4,15,16,20,21,22] to the case when, there are impulses and delay on the system, the right hand side is a multi-valued function and involving ψ−Hilfer fractional derivative.
In section 2, we prove some properties for ψ−fractional integral and ψ−fractional derivative, then we present, in Lemma 3.3, a relation between any solution of problem (1) and the corresponding fractional integral equation. In section 3, we prove an existence result of Problem (1). In section 4, we give another existence result of (1), then we investigate the ψ−generalized Ulam-Hyers stability of solutions. In the last section, examples are given to clarify the possibility of applicability of our assumptions.
Let Pck(E) be the family of non-empty convex and compact subsets of E. Since the given problem containing Hilfer derivative we need to the the spaces:
C1−μ,ψ(J,E):={x∈C(J∗,E):(ψ(.)−ψ(a))1−μx(.)∈C([a,b],E)}, |
and
Cn1−μ,ψ(J,E):={x∈Cn−1(J,E):[1ψ′(ϱ)ddϱ]n x∈C1−μ,ψ(J,E),n∈N. |
Obviously C1−μ,ψ(J,E) and Cn1−μ,ψ(J,E) are Banach spaces with norms
||x||C1−μ,ψ(J,E):=supϱ∈J|| (ψ(ϱ) −ψ(0))1−μx(ϱ)||, |
and
||x||Cn1−μ,ψ(J,E):=k=n−1∑k=1||x(k)||C(J,E)+||[1ψ′(ϱ)ddϱ]n x||C1−μ,ψ(J,E). |
Because Problem (1) involving impulses effect we recall the Banach space:
PC1−μ,ψ(J,E):={x:J∗→E, (ψ(.)−ψ(sk))1−μx(.)∈C(Jk,E),k=0,1,..,n,limϱ→s+k(ψ(ϱ)−ψ(sk))1−μx(ϱ)exists,x∈C(ϱi,E), andlimϱ→ϱ+ix(ϱ)exist,i=1,2,..,n}, |
endowed with the norm
||x||PC1−μ,ψ(J,E):=max{supϱ∈¯Jkk=0,1,,,,n,(ψ(ϱ)−ψ(sk))1−μ||x(ϱ)||E, supϱ∈¯ϱii=1,...,n||x(ϱ)||E}, |
where Jk=(sk,ϱk+1], ¯Jk=[sk,ϱk+1] (k=0,1,..,n),ϱi=(ϱi,si] and ¯ϱi=[ϱi,si] (i=1,2,..,n),
Next, the function χPC1−μ(J,E):Pb(PC1−μ,ψ(J,E))→[0,∞),given by by
χPC1−μ,ψ(J,E)(D):=max{maxk=0,1,..,nχC(¯Jk,E)(D|¯Jk),maxi=1,..,nχC(¯ϱi,E)(D|¯ϱi)} |
is a measure of noncompactness on PC1−μ(J,E), where
D∣¯Jk:={h∗∈C(¯Jk,E):h∗(ϱ)=(ψ(ϱ)−ψ(sk))1−μh(ϱ),ϱ∈Jk,h∗(sk)=limϱ→s+kh∗(ϱ),h∈D}, |
and
D∣¯ϱi:={h∗∈C(¯ϱi,E):h∗(ϱ)=h(ϱ),ϱ∈ϱi,h∗(ϱi)=h(ϱ+i),h∈D}. |
In the sequel, Iq,ψa+ denotes to the ψ−Riemann-Liouville fractional integral operator of order q with the lower limit at a, Dϑ,ψa+f to the ψ−Riemann-Liouville fractional derivative operator of order ϑ with the lower limit at a and cDϑ,ψa+f to the ψ−Caputo fractional derivative of order ϑ with the lower limit at a for f ∈AC1,ψ([a,b],E), where
AC1,ψ([a,b],E):={x:J→E,[1ψ′(ϱ)ddϱ] x∈AC(J,E)}. |
If f ∈C1,ψ([a,b],E):={x∈C(J,E):[1ψ′(ϱ)ddϱ] x∈C(J,E)}, then
cDϑ,ψa+f(ϱ):=I1−ϑ,ψa+[1ψ′(ϱ)ddϱ] f(ϱ),ϱ∈[a,b]. |
If ψ(ϱ)=ϱ, we obtain the Caputo fractional derivative, and if ψ(ϱ)=lnϱ, we obtain the Caputo-Hadamard fractional derivative. The following remark and more information about ψ−fractional integral and derivative can be found in [32,33,34,35]
Remark 2.1. If q=1, then I1,ψa+f(ϱ)=∫ϱaψ′(s)f(s) ds, and hence
1ψ′(ϱ)ddϱ I1,ψa+f(ϱ)=1ψ′(ϱ)ddϱ∫ϱa ψ′(s)f(s) ds=f(ϱ),a.e.for ϱ∈[a,b]. |
In the following lemma we give an important for Iϑ,ψa+, which we need later.
Lemma 2.1. Let 0<η≤ϑ. Then Iϑ,ψa+is bounded from Cη,ψ([a,b],E) into C([a,b],E) and
Iϑ,ψa+f(a)=limϱ→a+Iϑ,a+f(ϱ)=0,∀f∈Cη,ψ([a,b],E). |
Proof. The assumption f∈Cη,ψ([a,b],E) leads to (ψ(ϱ)−ψ(a))ηf(ϱ) is continuous on [a,b], and hence there is M>0 such that ||(ψ(ϱ)−ψ(a))ηf(ϱ)||≤M,∀ϱ∈ [a,b]. As a consequence, ||f(ϱ)||≤(ψ(ϱ)−ψ(a))−ηM,∀ϱ∈ (a,b]. Then
||Iϑ,ψa+f(ϱ)||=||∫ϱa (ψ(ϱ)−ψ(s))ϑ−1ψ′(s)Γ(ϑ)f(s) ds||≤MΓ(ϑ) |Iϑ,ψa+(ψ(ϱ)−ψ(a))−η|. |
In virtu of Lemma2 in[35], Iϑ,ψa+(ψ(ϱ)−ψ(a))−η=Γ(1−η)Γ(1−η+ϑ)(ψ(ϱ)−ψ(a))ϑ−η, and hence
||Iϑ,ψa+f(ϱ)||≤MΓ(ϑ)Γ(1−η)Γ(1−η+ϑ)(ψ(ϱ)−ψ(a))ϑ−η |
Since 0<η≤ϑ, we get limϱ→a+Iϑ,ψa+f(ϱ)=0, which means that Iϑ,ψa+f(ϱ) is continuous. Moreover, Iϑ,ψa+is bounded from Cη,ψ([a,b],E) into C([a,b],E).
Let us recall the definition of ψ−Hilfer fractional derivative.
Definition 1. [2] Let f∈L1([a,b],E) be such that I(1−ν)(1−ϑ),ψa+f∈AC1,ψ([a,b],E). The ψ−Hilfer fractional derivative of order 0<ϑ<1and type 0≤ν≤1 and with lower limit at a for a function f:[a,b]→E is defined by
Dϑ,ν,ψa+f(ϱ)=Iν(1−ϑ),ψa+[1ψ′(ϱ)ddϱ] (I(1−ν)(1−ϑ),ψa+f)(ϱ),ϱ∈[a,b], |
Denote
Cμ1−μ,ψ(J,E):={x∈C1−μ,ψ(J,E),Dμ,ψa+x∈C1−μ,ψ(J,E)}, |
Cϑ,ν1−μ,,ψ(J,E):={x∈C1−μ(J,E),Dϑ,ν,ψa+x∈C1−μ,ψ(J,E)}, |
PCμ1−μ,ψ(J,E):={x∈PC1−μ(J,E),Dμ,ψs+kx|Jk∈C1−μ,ψ(Jk,E),k=0,1,..,n}, |
and
PCϑ,ν1−μ,ψ(J,E):={x∈PC1−μ,ψ(J,E),Dϑ,ν,ψs+kx|Jk∈C1−μ,ψ(Jk,E), k=0,1,..,n}. |
Notice that, the operator Dϑ,ν,ψa+ can be written as:
Let
Dϑ,ν,ψa+f(ϱ)=Iν(1−ϑ),ψa+[1ψ′(ϱ)ddϱ](I1−μa+f(ϱ)=Iν(1−ϑ),ψa+Dμ,ψa+f(ϱ), μ=ϑ+ν−ϑν. |
So, if f∈Cμ1−μ,ψ([a,b],E), then by Lemma 2.1, Dϑ,ν,ψa+f(ϱ) exists ϱ∈[a,b].
Remark 2.2. Since Dϑ,ν,ψa+x=(Iν(1−ϑ),ψa+Dμ,ψa+x)(ϱ), it follows from lemma1. that Cμ1−μ,ψ([a,b],E) ⊆Cϑ,ν1−μ,ψ([a,b],E). Similarly, PCμ1−μ,ψ([a,b],E) ⊆PCϑ,ν1−μ,ψ([a,b],E).
Now, since our considered problem contains a delay we need to present the following spaces:
1- The normed space
Θ={z:[−r,0]→E suchthatz hasafinitenumberofdiscontinuitypointss≠0,allvaluesz(s+),and z(s−) arefinite}, |
endowed with the norm:
||z||Θ=∫0−r||z(s)||ds. |
2-The metric space (the space of solutions)
H:={x:[−r,b]→E,x|[−r,0]=Ψ, x|J∗∈PC1−μ,ψ(J,E)}, |
where the metric function is given by:
dH(x,y):=supϱ∈J||x(ϱ)−y(ϱ)||. |
3- The Banach space
H:={x:[−r,b]→Esuchthat x(ϱ)=0,∀ϱ∈[−r,0],x|J∗∈PC1−μ,ψ(J,E)}, |
endowed with the norm:
||x||H:=||x|J∗||PC1−μ,ψ(J,E)+||x|[−r,0]||Θ=||x|J∗||PC1−μ,ψ(J,E). |
Remark 2.3. (i) If x∈H, then x(0−)=0 and x(0+)=limϱ→0+(ψ(ϱ)−ψ(0))1−μx(ϱ).
(ⅱ) If x∈H, then x(0−)=Ψ(0) and x(0+−)=limϱ→0+(ψ(ϱ)−ψ(0))1−μx(ϱ). So, if Ψ(0)=0, then x will be continuous at zero.
It is easily seen that the function:
χH(B):=χPC1−μ(J,E){x|J:x∈B}, |
define a measure of noncompactness on H, where B is a bounded subset of H.
We need to the following fixed point for multi-valued functions.
Lemma 2.2. ([36], Theorem 3.1) Let W be a closed convex subset of a Banach space X and ϱ:W→Pc(W). Suppose that ϱ is closed, ϱ(D) is relatively compact, whenever D is compact, and that, for some x0∈W, one has
B⊆W, B=conv({x0}∪ϱ(B)), ¯B=¯C withC⊆B countable⟹B isrelativelycompact. |
Then, there is a fixed point for ϱ.
In this section, we demonstrate the existence of solutions of Problem (1). For any x∈H let
S1F(.,τ(.)x)={z∈L1(J,E):z(ϱ)∈F(ϱ,τ(.)x),a.e. for ϱ∈J }. |
and
Iν(1−ϑ),ψ(PC1−μ,ψ(J,E))={f:J→E, thereish∈PC1−μ,ψ(J,E)suchthatf(ϱ)=Iν(1−ϑ),ψsi+h(ϱ),ϱ∈Jk, k=0,1,2,..,n}. |
and for any x∈H let
¯x(ϱ):={Ψ(ϱ),ϱ∈[−r,0],x(ϱ),ϱ∈(0,b]. |
In order to derive the relation between any solution for Problem (1) and the corresponding fractional integral equation, we need to the following essential Lemmas.
Lemma 3.1. Let 0<ϑ<1,η∈[0,1).If f∈Cη,ψ([a,b],E) and I1−ϑ,ψa+f∈C1η,ψ([a,b],E), then
Iϑ,ψa+Dϑ,ψa+f(ϱ)=f(ϱ)−(ψ(ϱ)−ψ(a))ϑ−1I1−ϑ,ψa+f(a)Γ(ϑ),ϱ∈[a,b]. |
Proof. Since f∈Cη,ψ([a,b],E), then I1−ϑ,ψa+f(ϱ) is defined for ϱ∈[a,b]. Moreover, the assumption I1−ϑ,ψa+f∈C1η,ψ([a,b],E) implies to Dϑ,ψa+f(ϱ)∈Cη,ψ([a,b],E), and hence Iϑ,ψa+Dϑ,ψa+f(ϱ) and Iϑ+1,ψa+Dϑ,ψa+f(ϱ) are well defined for ϱ∈[a,b]. Observe that
Iϑ,ψa+Dϑ,ψa+f(ϱ)=1ψ′(ϱ)ddϱ(∫ϱa ψ′(s)Iϑ,ψa+Dϑ,ψa+f(s)ds)=1ψ′(ϱ)ddϱI1,ψa+Iϑ,ψa+Dϑ,ψa+f(ϱ)=1ψ′(ϱ)ddϱIϑ+1,ψa+Dϑ,ψa+f(ϱ),ϱ∈[a,b]. |
Now,
Iϑ+1,ψa+Dϑ,ψa+f(ϱ)=1Γ(ϑ+1)∫ϱa (ψ(ϱ)−ψ(s))ϑψ′(s)(1ψ′(s)dds)(I1−ϑ,ψa+f(s))ds=1Γ(ϑ+1)∫ϱa (ψ(ϱ)−ψ(s))ϑdds(I1−ϑ,ψa+f(s))ds. |
By integration by parts, we get
Iϑ+1,ψa+Dϑ,ψa+f(ϱ)=−(ψ(ϱ)−ψ(a))ϑI1−ϑ,ψa+f(a)Γ(ϑ+1)+1Γ(ϑ)∫ϱa (ψ(ϱ)−ψ(s))ϑ−1ψ′(s)I1−ϑ,ψa+f(s)ds=−(ψ(ϱ)−ψ(a))ϑI1−ϑ,ψa+f(a)Γ(ϑ+1)+Iϑ,ψa+I1−ϑ,ψa+f(ϱ)=I1,ψa+f(ϱ)−(ψ(ϱ)−ψ(a))ϑI1−ϑ,ψa+f(a)Γ(ϑ+1). |
It follows from Remark 2.1, that
Iϑ,ψa+Dϑ,ψa+f(ϱ)=1ψ′(ϱ)ddϱI1,ψa+f(ϱ)−1ψ′(ϱ)ddϱ(ψ(ϱ)−ψ(a))ϑ(I1−ϑ,ψa+f(a))Γ(ϑ+1)]=f(ϱ)−(ψ(ϱ)−ψ(a))ϑ−1I1−ϑ,ψa+f(a)Γ(ϑ). |
Lemma 3.2. Let α∈(0,1), η>0 and f∈Cη,ψ([a,b],E). Then Dα,ψa+Iα,ψa+f(ϱ)=f(ϱ),a.e. for ϱ∈[a,b].
Proof. In view of Remark 2.1, we get for a.e.ϱ∈[a,b]
Dα,ψa+Iα,ψa+f(ϱ)=[1ψ′(ϱ)ddϱ] I1−α,ψa+Iα,ψa+f(ϱ)=[1ψ′(ϱ)ddϱ] I1,ψa+f(ϱ)=f(ϱ). |
Now, we give in the following lemma the relation between any solution for Problem (1) and the corresponding fractional integral equation.
Lemma 3.3. Let 0<ϑ<1, 0≤ν≤1, μ=ϑ+ν−ϑν and x∈PCν(1−ϑ),ψ1−μ,ψ([0,b],E), gi:(ϱi,si]→E,i=1,......n, is continuous. The following hold.
(1) The function y:(0,b]→E given by
y(ϱ)={(ψ(ϱ)−ψ(0))μ−1Γ(μ)Ψ(0)+Iϑ,ψ0+x(ϱ),ϱ∈(0,ϱ1]gi(ϱ,x(ϱ−i)),ϱ∈(ϱi,si],i=1,......n,(ψ(ϱ)−ψ(si))μ−1Γ(μ)gi(si,x(ϱ−i))+Iϑ,ψsi+x(ϱ),ϱ∈(siϱi+1],i=1,......n, | (2) |
belongs to PCμ1−μ,ψ([0,b],E), Dϑ,ν,ψsi+y(ϱ) exists for any ϱ∈(si ϱi+1], i=0,1,......n, and verifies the ψ−Hilfer fractional problem:
{Dϑ,ν,ψs+iy(ϱ)=x(ϱ),ϱ∈∪i=ni=0(si,ϱi+1]I1−μ0+y(0+)=Ψ(0),y(ϱ+i)=gi(ϱi,x(ϱ−i)),i=1,......n,y(ϱ)=gi(ϱ,x(ϱ−i)),ϱ∈(ϱisi],i=1,......n,I1−μs+iy(s+i)=gi(si,x(ϱ−i)),i=1,......n, | (3) |
(2) If y∈PCμ1−μ,ψ([0,b],E) is a solution of (3), then y satisfies (2).
Proof. Let ϱ∈(0,ϱ1]. In view of (2)
y(ϱ)=(ψ(ϱ)−ψ(0))μ−1Γ(μ)Ψ(0)+Iϑ,ψ0+x(ϱ). | (4) |
According to Lemma 3.1 in[35], Dμ,ψ0+(ψ(ϱ)−ψ(a))μ−1=0. Then, by applying Dμ,ψa+ to both side of (4), it yields
D0+y(ϱ)=(1ψ′(ϱ)ddϱ)I(1−μ),ψ0+Iϑ,ψ0+x(ϱ)=(1ψ′(ϱ)ddϱ)I1−(μ−ϑ),ψ0+x(ϱ)=Dμ−ϑ,ψ0+x(ϱ)=Dν(1−ϑ),ψ0+x(ϱ). | (5) |
Observe that the assumption x∈PCν(1−ϑ),ψ1−μ,ψ([0,b],E) implies to Dν(1−ϑ),ψ0+x|J0∈C1−μ,ψ([0,ϱ1],E). It follows from (5) that Dμ,ψa+y∈C1−μ,ψ([0,ϱ1],E) and hence y|J0∈Cμ1−μ,ψ([0,ϱ1],E). Consequently y|J0∈Cϑ,ν1−μ,ψ([0,ϱ1],E) (see, Remark 2), and this assures that Dϑ,ν,ψ0+y(ϱ) is well defined for ϱ∈J0. Now, since Dν(1−ϑ),ψ0+x|J0∈C1−μ,ψ([0,ϱ1],E), then (1ψ′(ϱ)ddϱ)I1−ν(1−ϑ),ψ0+x|J0∈C1−μ,ψ([0,ϱ1],E). So, I1−ν(1−ϑ),ψ0+x|J0∈C11−μ,ψ([0,ϱ1],E). As a result from Lemma (3) one obtains
Dϑ,ν,ψ0+y(ϱ)=Iν(1−ϑ)0+Dμ,ψ0+y(ϱ)=Iν(1−ϑ)0+Dν(1−ϑ),ψ0+x(ϱ)=x(ϱ)−(ψ(ϱ)−ψ(0))ν(1−ϑ)−1Γ(ν(1−ϑ))(I1−ν(1−ϑ),ψ0+x)(0). |
Since I1−ν(1−ϑ),ψ0+x∈C1−μ,ψ([0,b],E) and 1−μ<1−ν(1−ϑ), then, by Lemma(1) (I1−ν(1−ϑ),ψ0+x)(0)=0. So, Dϑ,ν,ψ0+y(ϱ)=x(ϱ), ϱ∈[0,ϱ1]. It remains to demonstrate that y satisfies I1−μ0+y(0+)=ya. To do this, apply I1−μ,ψa+ to both side of (4)
I1−μ,ψ0+y(ϱ)=Ψ(0)Γ(μ)I1−μ,ψ0+(ψ(ϱ)−ψ(0))μ−1(0)+I1−μ,ψ0+Iϑ,ψa+x(ϱ)=Ψ(0)+I1−μ+ϑ,ψ0+Iϑ,ψa+x(ϱ). |
Because x∈C1−μ,ψ([0,ϱ1],E) and 1−μ<1−μ+ϑ, then, by Lemma(1),I1−μ+ϑ,ψ0+Iϑ,ψ0+x(0)=0. So, I1−μ,ψ0+y(0)=Ψ(0).
Similarly, we can show that for ϱ∈(si ϱi+1],i=1,......n,we have Dϑ,ν,ψs+iy(ϱ)=x(ϱ) and I1−μs+iy(s+i)=gi(si,x(ϱ−i)), and hence y is a solution for(3).
(2) Let y∈PCμ1−μ,ψ(J,E) be a solution of (3). Let i=0. Then y∈C1−μ,ψ(J0,E) and (1ψ′(ϱ)ddϱ)I1−μ,ψ0+y|J0∈C1−μ,ψ(J0,E). Thus, I1−μ,ψsiy|J0∈C11−μ,ψ(J0,E). By applying Lemma (3) it yields
Iμ,ψ0+Dμ,ψ0+y(ϱ)=y(ϱ)−(ψ(ϱ)−ψ(0))μ−1Γ(μ)(I1−μ,ψ0+x)(0)=y(ϱ)−(ψ(ϱ)−ψ(0))μ−1Γ(μ)Ψ(0),ϱ∈(0,ϱ1). | (6) |
Next, applying Iϑ,ψ0+ to both side of the equation Dϑ,ν,ψ0+y(ϱ)=x(ϱ), we get from (6)
Iϑ,ψ0+x(ϱ)=Iϑ,ψ0+Dϑ,ν,ψ0+y(ϱ)=Iϑ,ψ0+Iμ−ϑ,ψ0+Dμ,ψ0+y(ϱ)=Iμ,ψ0+Dμ,ψ0+y(ϱ)=y(ϱ)−(ψ(ϱ)−ψ(0))μ−1Γ(μ)Ψ(0),ϱ∈(0,ϱ1]. |
So,
y(ϱ)=(ψ(ϱ)−ψ(0))μ−1Γ(μ)Ψ(0)+Iϑ,ψ0+x(ϱ),ϱ∈(0,ϱ1]. |
Similarly, we can show that for i=1,......n,
y(ϱ)=(ψ(ϱ)−ψ(si))μ−1Γ(μ)gi(si,x(ϱ−i))+Iϑ,ψ0+x(ϱ),ϱ∈(siϱi+1]. |
Now, based on Lemma 3.3, we can give the concept of solutions for problem (1).
Definition 2. A function ¯x∈H is called a mild solution of (1) if there is f∈S1F(,τ(ϱ)¯x) such that ¯x satisfies the integral equation
¯x(ϱ)={Ψ(ϱ),ϱ∈[−r,0](ψ(ϱ)−ψ(0))μ−1Γ(μ)Ψ(0)+Iϑ,ψ0+f(ϱ),ϱ∈(0,ϱ1]gi(ϱ,x(ϱ−i)),ϱ∈(ϱi,si],i=1,......n,(ψ(ϱ)−ψ(si))μ−1Γ(μ)gi(si,x(ϱ−i))+Iϑ,ψs+if(ϱ),ϱ∈(siϱi+1],i=1,......n. |
Remark 3.1. The function ¯x is not necessarily continuous at the points si,i=0,1,......n. But if Ψ(0)=0, then it will be continuous at s0, and if gi(si,x(ϱ−i))=0, it will be continuous at si,i=1,......n.
In the following, we present our first existence result of solutions for Problem (1).
Theorem 3.1. Let F:J×Θ→Pck(E) be a multifunction, Ψ∈Θ, and gi:[ϱi,si]×E→E (i=1,2,....n). We assume the following conditions
(F1) For every z∈H, S1F(.,τ(.)z) is not empty subset of Iν(1−ϑ),ψ(PC1−μ,ψ(J,E)) and for almost every ϱ∈J,z⟶F(ϱ,τ(ϱ)z) is upper semicontinuous.
(F2) There is a φ∈L1q(I,R+),(0<q<ϑ) such that for any z∈H
‖F(ϱ,τ(ϱ)¯z‖≤φ(ϱ)(1+||z||H),a.e,ϱ∈J |
(F3) There is a ς∈L1q(I,R+),(0<q<ϑ) with the property that for any bounded subset D⊆H, any k=0,1,2,..,n, and a.e., for ϱ∈Jk
χE(F(ϱ,{τ(ϱ)¯x:x∈D))≤(ψ(ϱ)−ψ(sk))1−μς(ϱ)supθ∈[−r,0]χE{(τ(ϱ))¯x(θ):x∈D}, |
and
2κ1Γ(ϑ) η||ς||L1q(I,R+)<1, | (7) |
where η=δq(ψ(b)−ψ(0))ϑ−q(ϑ−q1−q)1−q, δ=maxs∈Jψ′(s), κ1=(ψ(b)−ψ(0))1−μ and χ is the Hausdorff measure of noncompactness on E.
(H1) for every i=1,2,…,n, gi:[ϱi,si]×E→E is uniformly continuous on bounded sets and for any ϱ∈J, gi (ϱ,.) maps any bounded subset of E to a relatively compact subset and there exists a positive constant hi such that for any x∈E
‖gi(ϱ,x))‖≤hi(ψ(ϱi)−ψ(si−1))1−μ‖x‖, ϱ∈(ϱi,si],x∈E. |
Then Problem (1) has a mild solution provided that
κ1Γ(ϑ)||φ||L1q(J,R+)η+h∗+h∗Γ(μ)<1, | (8) |
where, h∗=∑i=ni=0hi.
Proof. We define a multioperator Φ:H→P(H) as follows: let x∈H, then due to (F1) there is f ∈S1F(.,τ(.)¯x), and hence we can define y∈Φ(x) if and only if
y(ϱ)={0,ϱ∈[−r,0](ψ(ϱ)−ψ(0))μ−1Γ(μ)Ψ(0)+Iϑ,ψ0+f(ϱ),ϱ∈(0,ϱ1]gi(ϱ,x(ϱ−i)),ϱ∈(ϱi,si],i=1,......n,(ψ(ϱ)−ψ(0))μ−1Γ(μ)gi(si,x(ϱ−i))+Iϑ,ψs+if(ϱ),ϱ∈(siϱi+1],i=1,......n. | (9) |
Let us clarify that a point x is a fixed point for Φ if and only if ¯x is a solution for (1). Let x be a fixed point to Φ. Then
x(ϱ)={0,ϱ∈[−r,0](ψ(ϱ)−ψ(0))μ−1Γ(μ)Ψ(0)+Iϑ,ψ0+f(ϱ),ϱ∈(0,ϱ1]gi(ϱ,x(ϱ−i)),ϱ∈(ϱi,si],i=1,......n,(ψ(ϱ)−ψ(0))μ−1Γ(μ)gi(si,x(ϱ−i))+Iϑ,ψs+if(ϱ),ϱ∈(siϱi+1],i=1,......n. | (10) |
where f ∈S1F(.,τ(.)¯x). Therefore,
¯x(ϱ)={Ψ(ϱ),ϱ∈[−r,0](ψ(ϱ)−ψ(0))μ−1Γ(μ)Ψ(0)+Iϑ,ψ0+f(ϱ)ds,ϱ∈(0,ϱ1]gi(ϱ,x(ϱ−i)),ϱ∈(ϱi,si],i=1,......n,(ψ(ϱ)−ψ(0))μ−1Γ(μ)gi(si,x(ϱ−i))+Iϑ,ψs+if(ϱ),ϱ∈(siϱi+1],i=1,......n. |
This means ¯x satisfies (4), and hence it is a solution for (1). Similarly, it is easy to see that if ¯x satisfies (4), then x is a fixed point for Φ. So we prove, by application Lemma 3.3, that Φ has a fixed point. Obviously the values of Φ are convex.
Step1. We demonstrate that there is a n∈N with Φ(Bn)⊆Bn,where Bn={x∈H:‖x‖H≤n}. Suppose that for any natural number n, there are xn,yn∈H with yn∈Φ(xn),‖xn‖H≤n and ‖yn‖H>n. Then, there is fn∈ S1F(.,τ(.)¯xn), n≥1 such that
yn(ϱ)={0,ϱ∈[−r,0](ψ(ϱ)−ψ(0))μ−1Γ(μ)Ψ(0)+Iϑ,ψ0+fn(ϱ)∈(0,ϱ1]gi(ϱ,xn(ϱ−i)),ϱ∈(ϱi,si],i=1,......n,(ψ(ϱ)−ψ(si))μ−1Γ(μ)gi(si,xn(ϱ−i))+Iϑ,ψs+ifn(ϱ),ϱ∈(siϱi+1],i=1,......n. |
In view of (F2), we get for almost ϱ∈J
||fn(ϱ)||≤φ(ϱ)(1+‖xn‖H)≤φ(ϱ)(1+n). | (11) |
Then, if for almost ϱ∈(0,ϱ1], we get from Holder's inequality and (12)
(ψ(ϱ)−ψ(0))1−μ||yn (ϱ)||≤Ψ(0)Γ(μ)+(1+n)(ψ(b)−ψ(0))1−μΓ(ϑ)∫ϱ0(ψ(ϱ)−ψ(s))ϑ−1ψ′(s)φ(ϱ)ds |
Now, from Holder's inequality
∫ϱ0(ψ(ϱ)−ψ(s))ϑ−1ψ′(s)φ(ϱ)ds≤||φ||L1q(J,R+)[∫ϱ0(ψ(ϱ)−ψ(s))ϑ−11−q(ψ′(s))11−qds]1−q≤||φ||L1q(J,R+)[∫ϱ0(ψ(ϱ)−ψ(s))ϑ−11−qψ′(s)(ψ′(s))q1−qds]1−q≤δq||φ||L1q(J,R+)[∫ϱ0(ψ(ϱ)−ψ(s))ϑ−11−qψ′(s)ds]1−q≤δq(ψ(b)−ψ(0))ϑ−q(ϑ−q1−q)1−q=η. | (12) |
Then,
(ψ(ϱ)−ψ(0))1−μ||yn (ϱ)||≤Ψ(0)Γ(μ)+(1+n)(ψ(b)−ψ(0))1−μΓ(ϑ)∫ϱ0(ψ(ϱ)−ψ(s))ϑ−1ψ′(s)φ(ϱ)ds≤Ψ(0)Γ(μ)+κ1Γ(ϑ)||φ||L1q(J,R+)(1+n)η. | (13) |
If ϱ∈(ϱi si],=1,2,…,n, then by (H1)
||yn(ϱ)||=supϱ∈(ϱisi]||gi(ϱ,xn(ϱ−i))||≤h∗(ψ(ϱi) −ψ(si−1))1−μ||xn(ϱ−i)||≤h∗||xn||H≤hn. | (14) |
Similar as in (13), we get for almost ϱ∈(si ϱi+1],i=1,2,…,n.
(ψ(ϱ)−ψ(si))1−μ||yn(ϱ)||≤||gi(si,xn(ϱ−i))||Γ(μ)+ κ1Γ(ϑ)||φ||L1q(J,R+)(1+n)η≤h∗nΓ(μ)+ κ1Γ(ϑ)||φ||L1q(J,R+)(1+n)η. |
It follows from this inequality, (13) and (14), that
n<||yn||H≤Ψ(0)Γ(μ)+h∗n+h∗nΓ(μ)+κ1Γ(ϑ)||φ||L1q(J,R+)(1+n)η. |
By dividing both side by n and taking the limit as n→∞, one obtains
1<κ1Γ(ϑ)||φ||L1q(J,R+)η+h∗+h∗Γ(μ), |
which contradicts with (8). Thus, there is a natural number n0 such that Φ(Bn0)⊆Φ(Bn0).
Step2. Φ is closed on Bn0.
Let (xn)n≥1, (yn)n≥1be two sequences in Bn0 with xn→x in Bn0,yn→y in H and yn∈Φ(xn);n≥1. Then, there is fn ∈S1F(.,τ(.)¯xn) such that
yn(ϱ)={0,ϱ∈[−r,0](ψ(ϱ)−ψ(0))μ−1Γ(μ)Ψ(0)+1Γ(ϑ)∫ϱ0(ψ(ϱ)−ψ(s))ϑ−1ψ′(s)fn(s)ds,ϱ∈(0,ϱ1]gi(ϱ,xn(ϱ−i)),ϱ∈(ϱi,si],i=1,......n,(ψ(ϱ)−ψ(si))μ−1Γ(μ)gi(si,xn(ϱ−i))+1Γ(ϑ)∫ϱsi(ψ(ϱ)−ψ(s))ϑ−1ψ′(s)fn(s),ϱ∈(siϱi+1],i=1,......n. | (15) |
According to (11), it follows that ||fn(ϱ)||≤φ(ϱ)(1+n),a.e.ϱ∈J, and hence (fn)n≥1 is bounded in L1q(J,E). Since L1q(J,E) is reflexive, then, by using Mazur's Lemma, there is sequence, (zj)j≥1, of convex combinations of (fn)n≥1 converging strongly to f in L1(J,E) as j→∞. Notice that, by (F2) again, for every ϱ∈J,s∈(0,ϱ] and every n≥1
||(ϱ−s)ϑ−1fn(s)||≤ |ϱ−s)|ϑ−1φ(ϱ)(1+n)∈L1((0,ϱ],R+). |
Let
˜yn(ϱ)={0,ϱ∈[−r,0](ψ(ϱ)−ψ(0))μ−1Γ(μ)Ψ(0)+1Γ(ϑ)∫ϱ0(ψ(ϱ)−ψ(s))ϑ−1ψ′(s)zn(s)ds,ϱ∈(0,ϱ1]gi(ϱ,xn(ϱ−i)),ϱ∈(ϱi,si],i=1,......n,(ψ(ϱ)−ψ(si))μ−1Γ(μ)gi(si,xn(ϱ−i))+1Γ(ϑ)∫ϱsi(ψ(ϱ)−ψ(s))ϑ−1ψ′(s)zn(s),ϱ∈(siϱi+1],i=1,......n. | (16) |
Clearly ˜yn(ϱ)→y(ϱ), ϱ∈J and zn(ϱ)→f(ϱ), for almost ϱ∈J. Also, τ(ϱ)¯xn→τ(ϱ)¯x;ϱ∈J, and hence, by the upper semicontinuity of F(ϱ,.);a.e. ϱ∈J, it follows that f(ϱ)∈F(ϱ,τ(ϱ)¯x),a.e.[[37]. Theorem1, Sec. 4, Ch.1 ]. Therefore, by, the uniform continuity of gi(si,.) on bounded subsets and by passing to the limit as n→∞ in (15) we obtain, from the Lebesgue dominated convergence theorem y∈Φ(x).
Step3. We show that M∣¯Jk (k=0,1,...,n) and M∣¯ϱi (i=1,2,....,n) are equicontinuous, where
M∣¯Jk={z:¯Jk→E, z(ϱ)=(ψ(ϱ)−ψ(sk))1−μy(ϱ),ϱ∈Jk,z(sk)=limϱ→sk(ψ(ϱ)−ψ(sk))1−μy(ϱ), y∈Φ(x),x∈Bn0}, |
and
M∣¯ϱi={y∗∈C(¯ϱi,E):y∗(ϱ)=y(ϱ),ϱ∈(ϱi,si],y∗(ϱi)=y(ϱ+i),y∈Φ(x),x∈Bn0}. |
Case 1.let z∈M∣¯J0.Then there are x∈Bn0and y∈Φ(x) such that for ϱ∈(0,ϱ1],
y(ϱ)=(ψ(ϱ)−ψ(0))1−μΨ(0)Γ(μ)+1Γ(ϑ)∫ϱ0(ψ(ϱ)−ψ(s))ϑ−1ψ′(s)f(s)ds, |
z(ϱ)=(ψ(ϱ)−ψ(0))1−μy(ϱ) and z(0)=limϱ→0+(ψ(ϱ)−ψ(0))1−μy(ϱ), where f ∈S1F(.,τ(.)¯x). It follows
limδ→0+z(0+δ)=limδ→0+z(δ)=limδ→0+(ψ(δ)−ψ(0))1−μy(δ)=z(0). |
Let ϱ1,ϱ2 be two points in (0,ϱ1] be such that ϱ1<ϱ2. Then,
‖z(ϱ2)−z(ϱ1)‖≤||(ψ(ϱ2)−ψ(0))1−μΓ(ϑ)∫ϱ20(ψ(ϱ2)−ψ(s))ϑ−1ψ′(s)f(s)ds−(ψ(ϱ1)−ψ(0))1−μΓ(ϑ)∫ϱ10(ψ(ϱ1)−ψ(s))ϑ−1ψ′(s)f(s)ds||≤(ψ(ϱ2)−ψ(0))1−μΓ(ϑ)||∫ϱ2ϱ1(ψ(ϱ2)−ψ(s))ϑ−1ψ′(s)f(s)ds||+||∫ϱ10(ψ(ϱ1)−ψ(s))ϑ−1ψ′(s)||f(s)||ds× |(ψ(ϱ2)−ψ(0))1−μΓ(ϑ)−(ψ(ϱ1)−ψ(0))1−μΓ(ϑ)|+(ψ(ϱ1)−ψ(0))1−μΓ(ϑ)×||∫ϱ10[(ψ(ϱ2)−ψ(s))ϑ−1ψ′(s)−(ψ(ϱ1)−ψ(s))ϑ−1ψ′(s)]||f(s)||ds||. |
By the absolute integral of the Lebesgue integral and Holder's inequality, it yields from (12)
limϱ2→ϱ1(ψ(ϱ2)−ψ(0))1−μΓ(ϑ)∫ϱ2ϱ1(ψ(ϱ2)−ψ(s))ϑ−1ψ′(s)||f(s)||ds≤κ1Γ(ϑ)(1+n0)limϱ2→ϱ1∫ϱ2ϱ1(ψ(ϱ2)−ψ(s))ϑ−1ψ′(s)φ(s)ds≤κ1Γ(ϑ)(1+n0)=0, |
independent of x.
Since ψ is continuous, we get by (12)
||∫ϱ10(ψ(ϱ1)−ψ(s))ϑ−1ψ′(s)Γ(ϑ)f(s)ds||× |(ψ(ϱ2)−ψ(0))1−μ−(ψ(ϱ1)−ψ(0))1−μ|≤(1+n0)∫ϱ10(ψ(ϱ1)−ψ(s))ϑ−1ψ′(s)Γ(ϑ)φ(s)ds×limϱ2→ϱ1|(ψ(ϱ2)−ψ(0))1−μ−(ψ(ϱ1)−ψ(0))1−μ|=0 |
Moreover,
limϱ2→ϱ1(ψ(ϱ1)−ψ(0))1−μΓ(ϑ)×||∫ϱ10[(ψ(ϱ2)−ψ(s))ϑ−1ψ′(s)−(ψ(ϱ1)−ψ(s))ϑ−1ψ′(s)]f(s)ds||≤κ1Γ(ϑ)(1+n0)×limϱ2→ϱ1∫ϱ10[(ψ(ϱ2)−ψ(s))ϑ−1ψ′(s)−(ψ(ϱ1)−ψ(s))ϑ−1ψ′(s)]φ(s)ds≤κ1Γ(ϑ)(1+n0)||φ||L1q(I,R+)×limϱ2→ϱ1[∫ϱ10[(ψ(ϱ2)−ψ(s))ϑ−1ψ′(s)−(ψ(ϱ1)−ψ(s))ϑ−1ψ′(s)]11−qds]1−q |
Put ¯ω=ϑ−11−q∈(−1,0), then for s<ϱ1<ϱ2, we have (ψ(ϱ1)−ψ(s))¯ω≥(ψ(ϱ2)−ψ(s))¯ω. By applying Lemma 3.1 in [43] and taking into account 1−q∈(0,1), we get
|[(ψ(ϱ1)−ψ(s))¯ω]1−q−[(ψ(ϱ2)−ψ(s))¯ω]1−q| ≤[(ψ(ϱ1)−ψ(s))¯ω−(ψ(ϱ2)−ψ(s))¯ω]1−q. |
Then,
|(ψ(ϱ1)−ψ(s))ϑ−1−(ψ(ϱ2)−ψ(s))ϑ−1|≤ [(ψ(ϱ1)−ψ(s))¯ω−(ψ(ϱ2)−ψ(s))¯ω]1−q. |
This leads to
|(ψ(ϱ1)−ψ(s))ϑ−1−(ψ(ϱ2)−ψ(s))ϑ−1|11−q≤(ψ(ϱ1)−ψ(s))¯ω−(ψ(ϱ2)−ψ(s))¯ω. |
Therefore,
limϱ2→ϱ1[∫ϱ10[(ψ(ϱ2)−ψ(s))ϑ−1ψ′(s)−(ψ(ϱ1)−ψ(s))ϑ−1ψ′(s)]11−qds]1−q≤limϱ2→ϱ1[∫ϱ10[(ψ(ϱ2)−ψ(s))ϑ−1−(ψ(ϱ1)−ψ(s))ϑ−1]11−q(ψ′(s))11−qds]1−q≤limϱ2→ϱ1[∫ϱ10[(ψ(ϱ1)−ψ(s))¯ω−(ψ(ϱ2)−ψ(s))¯ω](ψ′(s))11−qds]1−q≤δqlimϱ2→ϱ1[∫ϱ10[(ψ(ϱ1)−ψ(s))¯ω−(ψ(ϱ2)−ψ(s))¯ω]ψ′(s)ds]1−q=δqlimϱ2→ϱ1[(ψ(ϱ1)−ψ(s))¯ω+1¯ω+1−(ψ(ϱ2)−ψ(s))¯ω+1¯ω+1]1−q=0. |
Thus,
limϱ2→ϱ1(ψ(ϱ1)−ψ(0))1−μΓ(ϑ)×||∫ϱ10[(ψ(ϱ2)−ψ(s))ϑ−1ψ′(s)−(ψ(ϱ1)−ψ(s))ϑ−1ψ′(s)]f(s)ds||=0 |
independent of x.
Case 2. Let i∈{1,2,..,n} be fixed and set Ki={x(ϱ−i):x∈Bn0}.ObviouslyKi is bounded subset of E.
Let y∈M∣ϱi. Then, there is x∈Bn0 such that
y(ϱ)=gi(ϱ,x(ϱ−i)),ϱ∈(ϱi,si] |
It follows, from the uniform continuity of gi on the bounded set [ϱi,si]×Ki, that for ϱ1,ϱ2∈(ϱi,si]
limϱ2→ϱ1||y(ϱ2)−y(ϱ1)||=limϱ2→ϱ1||gi(ϱ2,x(ϱ−i))−gi(ϱ1,x(ϱ−i))||=0, |
independent of x. When ϱ=ϱi, let δ>0 be such that ϱi+δ∈(ϱi,si] and λ>0 with ϱi<λ<ϱi+δ≤si. Then, we have
limδ→0+‖y∗(ϱi+δ)−y∗(ϱi)‖=limδ→0+limλ→ϱ+i‖y(ϱi+δ)−y(λ)‖=0. |
Case 3. let k=1,…,n be fixed, z∈M∣¯Jk. Then, there are x∈Bn0and y∈Φ(x) such that for ϱ∈ϱ∈(sk,ϱk+1],
y(ϱ)=(ψ(ϱ)−ψ(sk))μ−1Γ(μ)gk(sk,xn(ϱ−k))+1Γ(ϑ)∫ϱsk(ψ(ϱ)−ψ(s))ϑ−1ψ′(s)zn(s), |
z(ϱ)=(ψ(ϱ)−ψ(sk))1−μy(ϱ) and z(sk)=limϱ→sk+(ψ(ϱ)−ψ(sk))1−μy(ϱ), where f ∈S1F(.,τ(.)¯x). It follows
limδ→0+z(sk+δ)=limδ→0+(ψ(sk+δ)−ψ(sk))1−μy(sk+δ)=limϱ→sk+(ψ(ϱ)−ψ(sk))1−μy(ϱ)=z(sk). |
Next, let ϱ1,ϱ2∈(sk,ϱk+1](ϱ1< ϱ2).By arguing as in case 1, one can show that
limϱ2→ϱ1‖z(ϱ+δ)−z(ϱ)‖=0, |
independent of x.
Step 4. Let K⊆Bn0, K=conv({0}∪Φ(K)), ¯Z=¯C with C⊆K countable. We have to show that K is relatively compact in H Let D={yn:n≥1}⊆Φ(K) with C⊆conv({x0}∪D), xn∈K with yn∈Φ(xn).This means that, there is fn∈S1F(.,τ(.)¯xn) such that the relation (15) holds. Observe that, from (F3) it holds for a.e.s∈J0
χE{fn(s):n≥1}≤(ψ(ϱ)−ψ(0))1−μς(ϱ)χ{F(s,τ(s)¯xn):n≥1}≤ζ(s)(ψ(ϱ)−ψ(0))1−μsupθ∈[−r,0] χ{¯xn(s+θ):n≥1}≤ζ(s)(ψ(ϱ)−ψ(0))1−μsupδ∈[−r,s] χ{¯xn(δ):n≥1}. |
Because supδ∈[−r,0] χ{¯xn(δ):n≥1}=supδ∈[−r,0] χ{Ψ(δ)}=0. Thus,
χE{fn(s):n≥1}≤ζ(s)(ψ(ϱ)−ψ(0))1−μ supδ∈[0,s]χ{xn(δ):n≥1}=ζ(s)supδ∈[0,s]χ{(ψ(ϱ)−ψ(0))1−μxn(δ):n≥1}≤ζ(s)χPC1−μ,ψ(J,E){xn:n≥1}≤ζ(s)χH(Bn−1). | (17) |
According to the definition χH(Z), one obtains
χH(Z)=χH(¯Z)=χH(¯C)=χH(C))≤χH(conv({x0}∪D))=χH(D)=max{maxk=0,1,..,nχC(¯Jk,E)(D|¯Jk),maxi=1,..,nχC(¯ϱi,E)(D|¯ϱi)}. | (18) |
In view of Step 3, D|¯Ji and D|¯ϱi are equicontinuous, consequently
χH(K)≤max{maxi=0,1,..,nmaxϱ∈¯Jkχ{y∗n(ϱ):n≥1},maxi=1,..,nmaxϱ∈¯ϱiχ{y∗n(ϱ):n≥1}}, | (19) |
Let ϱ∈(0,ϱ1]. Then,
y∗n(ϱ)=Ψ(0)Γ(μ)+(ψ(ϱ)−ψ(0))1−μΓ(ϑ)∫ϱ0(ψ(ϱ)−ψ(0))ϑ−1ψ′(s)fn(s)ds, |
which yields with (17)
χ{y∗n(ϱ):n≥1}≤(ψ(ϱ)−ψ(0))1−μΓ(ϑ)∫ϱ0(ψ(ϱ)−ψ(0))ϑ−1ψ′(s)χ{fn(s):n≥1}ds≤(ψ(b)−ψ(0))1−μΓ(ϑ)χH(Bn−1)∫ϱ0(ψ(ϱ)−ψ(0))ϑ−1ψ′(s)ζ(s)ds≤κ1Γ(ϑ) η||ς||L1q(I,R+)χH(Bn−1). | (20) |
Notice that
χ{y∗n(0):n≥1}=χ{limϱ→0+y∗n(ϱ):n≥1}=χ{Ψ(0)Γ(μ)}=0. | (21) |
Moreover, from the fact that ||xn||H≤n0,the set {xn(ϱ−i):n≥1} is bounded for every i=1,2,...,n and hence from the assumption (H1), we get
χ{gi(ϱ,xn(ϱ−i)),n≥1}=0,ϱ∈(ϱi,si],i=1,...,n. | (22) |
and
χ{gi(ϱi,xn(ϱ−i)),n≥1}=χ{gi(si,xn(ϱ−i)):n≥1}=0. | (23) |
Similarly,
χ{y∗n(si):n≥1}=χ{limϱ→s+i(ψ(ϱ)−ψ(si))1−μyn(ϱ),:n≥1}=χ{1Γ(μ)gi(si,xn(ϱ−i)):n≥1}=0,i=1,2,..,n. | (24) |
By arguing as in (20), one can show that
maxϱ∈¯Jkχ{y∗n(ϱ):n≥1}≤κ1Γ(ϑ) η||ς||L1q(I,R+)χH(Bn−1),k=1,2,..n. | (25) |
From the relations (7), (19)–(25), it follows that
χH(K)≤χH(K)2κ1Γ(ϑ) η||ς||L1q(I,R+) <χH(K). |
So, χPC1−μ(J,E)(K)=0, and hence K is relatively compact.
Step5. Φ maps compact sets into relatively compact sets.
Let B be a compact subset of Bn0, {yn,:n≥1}⊆Φ(B) Then, there is xn∈B,n≥1, such that yn ∈Φ(xn). So, there is fn∈S1F(.,¯xn(.)) such that (15) holds. We have to show that the set Z={yn:n≥1} is relatively compact in H. Note that, since B is compact in H, then from (F3) we get for a.e. s∈J0,
χE{fn(s):n≥1}≤(ψ(ϱ)−ψ(0))1−μς(ϱ)χ{F(s,τ(s)¯xn):n≥1}≤ζ(s)(ψ(ϱ)−ψ(0))1−μsupθ∈[−r,0] χ{¯xn(s+θ):n≥1}≤ζ(s)(ψ(ϱ)−ψ(0))1−μsupδ∈[−r,s] χ{¯xn(δ):n≥1}≤ζ(s)(ψ(ϱ)−ψ(0))1−μ supδ∈[0,s]χ{xn(δ):n≥1}=ζ(s)supδ∈[0,s]χ{(ψ(ϱ)−ψ(0))1−μxn(δ):n≥1}≤ζ(s)χPC1−μ,ψ(J,E){xn:n≥1}≤ζ(s)χH(B)=0. |
By the same reasons, one can show that for a.e. s∈Jk, k=1,2,..,n
χE{fn(s):n≥1}=0. |
By arguing as in the previous step one can show that Z is relatively compact, and hence Φ(B) is relatively compact.
Now, by applying Lemma 8, there is x∈H andf ∈S1F(.,τ(.)¯x) such that
x(ϱ)={0,ϱ∈[−r,0](ψ(ϱ)−ψ(0))μ−1Γ(μ)Ψ(0)+1Γ(ϑ)∫ϱ0(ψ(ϱ)−ψ(s))ϑ−1ψ′(s)f(s)ds,ϱ∈(0,ϱ1]gi(ϱ,x(ϱ−i)),ϱ∈(ϱi,si],i=1,......n,(ψ(ϱ)−ψ(0))μ−1Γ(μ)gi(si,x(ϱ−i))+1Γ(ϑ)∫ϱsi(ψ(ϱ)−ψ(s))ϑ−1ψ′(s)f(s),ϱ∈(siϱi+1],i=1,......n. |
Next, in view of (F1), there is h∈PC1−μ,ψ(J,E) such that f(ϱ)=Iν(1−ϑ),ψsi+h(ϱ),ϱ∈Jk,k=0,..,n, and hence, from Lemma 3.3,
Dν(1−ϑ),ψsi+f(ϱ)=Dν(1−ϑ),ψsi+Iν(1−ϑ)ψsi+h(ϱ)=h(ϱ),ϱ∈Jk,k=0,..,n. |
This yields that f∈PCν(1−ϑ)1−μ,ψ(J,E). Then, the function
¯x(ϱ)={Ψ(ϱ),ϱ∈[−r,0](ψ(ϱ)−ψ(0))μ−1Γ(μ)Ψ(0)+1Γ(ϑ)∫ϱ0(ψ(ϱ)−ψ(s))ϑ−1ψ′(s)f(s)ds,ϱ∈(0,ϱ1]gi(ϱ,x(ϱ−i)),ϱ∈(ϱi,si],i=1,......n,(ψ(ϱ)−ψ(si))μ−1Γ(μ)gi(si,x(ϱ−i))+1Γ(ϑ)∫ϱsi(ψ(ϱ)−ψ(s))ϑ−1ψ′(s),ϱ∈(siϱi+1],i=1,......n. | (26) |
belongs to H and in view of Lemma(6) it is a solution for (1). This completes the proof.
Remark 3.2. Theorem 1 remains true if condition (F2) is replaced by the following assumption:
(F2)∗There is a φ∈L1q(I,R+) such that for any z∈Θ
‖F(ϱ,z)‖≤{φ(ϱ) (1+(ψ(ϱ)−ψ(si))1−μ‖z(0)‖E),a.e.ϱ∈∪i=ni=0(siϱi+1],φ(ϱ) (1+‖z(0)‖E),a.e.ϱ∈∪i=ni=1(ϱi,si]. |
In fact, condition (F2) is used only to prove relation (11). We show this relation by using (F2)∗. Let fn∈ S1F(.,τ(.)¯xn), n≥1 Then, by (F2)∗, for almost ϱ∈∪i=ni=0(si ϱi+1] we get
||fn(ϱ)||≤φ(ϱ)(1+(ψ(ϱ)−ψ(si))1−μ||τ(ϱ)¯xn(0)||)≤φ(ϱ)(1+(ψ(ϱ)−ψ(si))1−μ||¯xn(ϱ)||)≤φ(ϱ)(1+(ψ(ϱ)−ψ(si))1−μ||xn(ϱ)||≤φ(ϱ)(1+‖xn‖H). |
and for almost ϱ∈∪i=ni=1(ϱi si]
||fn(ϱ)||≤φ(ϱ)(1+||τ(ϱ)¯xn(0)||E)≤φ(ϱ)(1+¯xn(ϱ))≤φ(ϱ)(1+||xn(ϱ)||)≤φ(ϱ)(1+‖xn‖H). |
In this section, we give another version for the existence of solutions and investigate the generalized ψ−generalized Ulam-Hyers stability of problem (1). For basic information about multivalued weakly Picard operators we refer to [38].
Definition 3. [39,40] A An increasing function ς:[0,∞)→[0,∞) is called comparison function if limn→∞ςn(s)=0, ∀s∈[0,∞),where ςn(s)=ςn−1(ς(s)).If in addition ∞∑n=1ςn(ϱ)<∞,∀s∈(0,∞), then ς:[0,∞)→[0,∞) is called strictly comparison.
Remark 4.1. [41]
1- If ς:[0,∞)→[0,∞) is comparisons, then ς(s)<s,∀s>0,ς(0)=0 and ς is continuous at 0.
2-the functions ς1(s)=cϱ;c∈[0,1) and ς2(s)=ss+1:s∈[0,∞) are strictly comparison
Let M:={y∈H:y|J∈PCμ1−μ,ψ(J,E), Dϑ,ν,ψsi+y(ϱ) exists for any ϱ∈(si ϱi+1], i=0,1,......n,}
Definition 4. [12] Problem (1) is called ψ−generalized Ulam- Hyers stable if there is a continuous function θ:[0,∞)→[0,∞) and θ(0)=0 such that for each ϵ>0 and each solution y∈M of the inequality
{d(Dϑ,ν,ψs+iy(ϱ),F(ϱ,τ(ϱ)y))≤ϵ,a.e.ϱ∈(si,ϱi+1],i=0,1,....,n,||y(ϱ)−gi(ϱ,y(ϱ−i))||≤ϵ,ϱ∈(ϱisi],i=1,......n,,||y(ϱ+i)−gi(ϱi,y(ϱ−i))||≤ϵ,i=1,......n, | (27) |
there is a solution x∈H for (1) with
supϱ∈J||x(ϱ)−y(ϱ)||≤θ(ϵ). | (28) |
In the following theorem, we establish the existence and generalized ψ− Ulam-Hyers stability of solutions Problem (1).
Theorem 4.1. Let F:J×Θ→Pck(E), Ψ∈Θ andgi:[ϱi,si]×E→E (i=1,2,....n)be such that:
(F4) For every z∈H, S1F(.,τ(.)z) is a non-empty subset of Iν(1−ϑ),ψ(PC1−μ,ψ(J,E)).
(F5) There is a function σ∈L1q(I,R+),0<q<ϑ and a strict comparison function ς:[0,∞)→[0,∞) such that
(ⅰ) For every x1,x2∈H
h(F(ϱ,τ(ϱ)¯x1),F(ϱ,τ(ϱ)¯x2))≤σ(ϱ)ς(‖x1−x2‖)H,∀ϱ∈J. |
(ⅱ)
Sup{ ||y|| : y∈F(ϱ,τ(ϱ)¯x0)}≤σ(ϱ),fora.e.ϱ∈J, |
where ¯x0(s)=Ψ(s);s∈[−r,0] and ¯x0(s)=0;s∈(0,b].
(H2)For anyi=1,2,....n, there ξi >0 such that for any ϱ∈[ϱi,si]
‖gi(ϱ,x)−gi(ϱ,y)‖≤ξi ς(‖x−y‖E),∀x,y∈E. |
Then, there is solution for Problem(1) provided that
ξΓ(μ)+ ηδq(ψ(b)−ψ(0))1−μΓ(ϑ)||σ||L1q(I,R+)<1, | (29) |
where ξ=∑i=ni=1ξi.Moreover, if, in addition there, is a c>1such that ς(cϱ)≤c ς(ϱ);ϱ∈[0,∞) and ϱ=0 is a point of uniform convergence for the series ∞∑n=1ςn(ϱ), then problem (1) is a ψ−generalized Ulam-Hyers stable.
Proof. Condition (F4) allows to define a multifunction Φ:H→P(H) as in (9) Note that by (F5), for every x∈H, and for a.e.ϱ∈J
||F(ϱ,τ(ϱ)¯x)||=h(F(ϱ,τ(ϱ)¯x),{0})≤h(F(ϱ,τ(ϱ)¯x),F(ϱ,τ(ϱ)¯x0))+h(F(ϱ,τ(ϱ)¯x0),{0})≤σ(ϱ)ς(‖x‖H)+σ(ϱ)≤σ(ϱ)(1+ς(‖x‖H)),for a.e.ϱ∈J. |
Then, as the reasons of Theorem 1, we can show that the values of Φ are closed. We show that Φ is ς−contraction. Let x2,x1∈H and y1∈Φ(x1). Then, there exists f∈S1F(.,τ(.)¯x1)such that
y1(ϱ)={0,ϱ∈[−r,0](ψ(ϱ)−ψ(0))μ−1Γ(μ)Ψ(0)+1Γ(ϑ)∫ϱ0(ψ(ϱ)−ψ(s))ϑ−1ψ′(s)f(s)ds,ϱ∈(0,ϱ1]gi(ϱ,x(ϱ−i)),ϱ∈(ϱi,si],i=1,......n,(ψ(ϱ)−ψ(si))μ−1Γ(μ)gi(si,x(ϱ−i))+1Γ(ϑ)∫ϱsi(ψ(ϱ)−ψ(s))ϑ−1ψ′(s)f(s),ϱ∈(siϱi+1],i=1,......n. |
Define Π:J→2E as
Π(ϱ)={u∈E:||f (ϱ)−u || ≤σ(ϱ)ς(‖x1−x2‖H)}. |
From (F5)(i), we have
h(F(ϱ,τ(ϱ)¯x1),F(ϱ,τ(ϱ)¯x2))≤σ(ϱ)ς(‖x1−x2‖H),ϱ∈J, |
Since the values of F are compact, there exists uϱ ∈F(ϱ,τ(ϱ)x2) such that
||uϱ−f(ϱ)|| ≤σ(ϱ)ς(‖x1−x2‖H), |
as a consequence Π(ϱ);ϱ∈J is not empty. Furthermore, because f,σ,x1,x2 are measurable and E is separable, it follows from [[42], Theorem Ⅲ.41], the multifunction s→ Π(s)∩F(s,τ(s)¯x2) is measurable and since its values are non-empty and compact, there is h∈S1F(ϱ,τ(ϱ)¯x2) such that
||h(ϱ)−f(ϱ)|| ≤ σ(ϱ)ς(‖x1−x2‖H),a.e.ϱ∈J. |
Set
y2(ϱ)={0,ϱ∈[−r,0](ψ(ϱ)−ψ(0))μ−1Γ(μ)Ψ(0)+1Γ(ϑ)∫ϱ0(ψ(ϱ)−ψ(s))ϑ−1ψ′(s)h(s)ds,ϱ∈(0,ϱ1]gi(ϱ,x2(ϱ−i)),ϱ∈(ϱi,si],i=1,......n,(ψ(ϱ)−ψ(si))μ−1Γ(μ)gi(si,x2(ϱ−i))+1Γ(ϑ)∫ϱsi(ψ(ϱ)−ψ(s))ϑ−1ψ′(s)h(s),ϱ∈(siϱi+1],i=1,......n. |
Notice that, y2∈ Φ(x2) and if ϱ∈J0,, we get from Holder's inequality
(ψ(ϱ)−ψ(0))1−μ||y2(ϱ)−y1(ϱ)|| ≤(ψ(ϱ)−ψ(0))1−μΓ(ϑ)∫ϱ0(ψ(ϱ)−ψ(s))ϑ−1ψ′(s)σ(s) ||h(s)−f(s)||ds≤(ψ(ϱ)−ψ(0))1−μΓ(ϑ)ς(‖x1−x2‖H)∫ϱ0(ψ(ϱ)−ψ(s))ϑ−1ψ′(s)σ(s)ds≤(ψ(ϱ)−ψ(0))1−μΓ(ϑ) ς(‖x1−x2‖H)||σ||L1q(I,R+)δq(ψ(b)−ψ(0))ϑ−q(ϑ−q1−q)1−q≤ηδq(ψ(b)−ψ(0))1−μΓ(ϑ)||σ||L1q(I,R+)ς(‖x1−x2‖H). | (30) |
Similarly, if ϱ∈(si ϱi+1],i=1,......,n, we get
(ψ(ϱ)−ψ(si))||y2(ϱ)−y1(ϱ)||≤1Γ(μ)||gi(si,x1(ϱ−i))−gi(si,x2(ϱ−i))||+ηδq(ψ(b)−ψ(0))1−μΓ(ϑ)||σ||L1q(I,R+)ς(‖x1−x2‖H).≤ξΓ(μ)ς(‖x1−x2‖H)+ηδq(ψ(b)−ψ(0))1−μΓ(ϑ)||σ||L1q(I,R+)ς(‖x1−x2‖H).≤ς(‖x1−x2‖H)[ξΓ(μ)+ ηδq(ψ(b)−ψ(0))1−μΓ(ϑ)||σ||L1q(I,R+)], | (31) |
where ξ=∑i=ni=1ξi. Furthermore, If ϱ∈(ϱi,si],i=1,......n, then by (H2)one obtains
||y2(ϱ)−y1(ϱ)||≤||gi(ϱ,x(ϱ−i))−gi(ϱ,x(ϱ−i))||≤ξς(||x1(ϱ−i)−x2(ϱ−i)||) ≤ξ ς(‖x1−x2‖H). | (32) |
It follows from (30)–(32),
hH(Φ(x2),Φ(x1))<ς(‖x1−x2‖H)[ξΓ(μ)+ ηδq(ψ(ϱ)−ψ(0))1−μΓ(ϑ)||σ||L1q(I,R+)], |
where, hH is the Hausdorff distance in H. This inequality and (29) imply that Φ is ς−contraction and thus by Theorem 3.1(i) in [41], Φ has a fixed point and as in Theorem 1, this fixed point is a solution for Problem (1).
Now, in order to demonstrate that Problem(1) is a ψ−generalized Ulam-Hyers stable. Let ϵ>0 and y∈M be a solution of inequality (27). Because the values of F are compact, there f∈S1F(.,τ(ϱ)y) such that
||Dϑ,ν,ψs+iy(ϱ)−f(ϱ)||=d(Dϑ,ν,ψs+iy(ϱ),F(ϱ,τ(ϱ)y)),a.e.ϱ∈(si,ϱi+1],i=0,1,....,n. |
Then, for almost everywhere ϱ∈∪i=ni=0(si,ϱi+1]
Dϑ,ν,ψs+iy(ϱ)=w(ϱ)+f(ϱ), |
where w∈PCν(1−ϑ)1−μ,ψ(J,E) and ||w(ϱ)|≤ϵ,∀ϱ∈J. Furthermore, according to (F4), there is h∈PC1−μ,ψ(J,E) such that f(ϱ)=Iν(1−ϑ),ψsi+h(ϱ),ϱ∈Jk, k=0,1,2,..,n, and hence, from Lemma 3.2,
Dν(1−ϑ),ψsi+f(ϱ)=Dν(1−ϑ),ψsi+Iν(1−ϑ)ψsi+h(ϱ)=h(ϱ),ϱ∈Jk, k=0,1,2,..,n. |
This yields that f∈PCν(1−ϑ)1−μ,ψ(J,E). Therefore, y(ϱ)=Ψ(ϱ),ϱ∈[−r,0] and
{Dϑ,ν,ψs+iy(ϱ)=f(ϱ)+ϑ(ϱ),a.e.ϱ∈(si,ϱi+1],i=0,1,....,n,y(ϱ)=gi(ϱ,y(ϱ−i))+ϵ,ϱ∈(ϱisi],i=1,......n,,y(ϱ+i)=gi(ϱi,y(ϱ−i))+ϵ,i=1,......n, | (33) |
In view of the second assertion of Lemma 3.3, relation (33) one obtains
y(ϱ)={Ψ(ϱ),ϱ∈[−r,0] (ψ(ϱ)−ψ(0))μ−1Γ(μ)Ψ(0)+1Γ(ϑ)∫ϱ0(ψ(ϱ)−ψ(s))ϑ−1f(s)ds+1Γ(ϑ)∫ϱ0(ψ(ϱ)−ψ(s))ϑ−1ϑ(ϱ)ds,ϱ∈(0,ϱ1],gi(ϱ,y(ϱ−i))+ϵ,ϱ∈(ϱisi],i=1,2,…,n,(ψ(ϱ)−ψ(si))μ−1Γ(μ)y(si)+1Γ(ϑ)∫ϱ0(ψ(ϱ)−ψ(s))ϑ−1ψ′(s)(f(s))ds+1Γ(ϑ)∫ϱ0(ψ(ϱ)−ψ(s))ϑ−1ψ′(s)ϑ(ϱ)ds,ϱ∈(siϱi+1],i=1,2,…,n. | (34) |
Let
z(ϱ)={0,ϱ∈[−r,0] (ψ(ϱ)−ψ(0))μ−1Γ(μ)Ψ(0)+1Γ(ϑ)∫ϱ0(ψ(ϱ)−ψ(s))ϑ−1ψ′(s)f(s)ds,ϱ∈(0,ϱ1],gi(ϱ,y(ϱ−i)),ϱ∈(ϱisi],i=1,2,…,n,(ψ(ϱ)−ψ(si))μ−1Γ(μ)gi(si,y(ϱ−i))+1Γ(ϑ)∫ϱ0(ψ(ϱ)−ψ(s))ϑ−1ψ′(s)(f(s))ds|+1Γ(ϑ)∫ϱ0(ψ(ϱ)−ψ(s))ϑ−1ψ′(s)ϑ(ϱ)ds,ϱ∈(siϱi+1],i=1,2,…,n. |
Obviously, z∈Φ(y∗), where
y∗(ϱ)={0,ϱ∈[−r,0] (ψ(ϱ)−ψ(0))μ−1Γ(μ)Ψ(0)+1Γ(ϑ)∫ϱ0(ψ(ϱ)−ψ(s))ϑ−1ψ′(s)f(s)ds+1Γ(ϑ)∫ϱ0(ψ(ϱ)−ψ(s))ϑ−1ϑ(ϱ)ds,ϱ∈(0,ϱ1],gi(ϱ,y(ϱ−i))+ϵ,ϱ∈(ϱisi],i=1,2,…,n,(ψ(ϱ)−ψ(si))μ−1Γ(μ)y(si)+1Γ(ϑ)∫ϱ0(ψ(ϱ)−ψ(s))ϑ−1ψ′(s)(f(s))ds|+1Γ(ϑ)∫ϱ0(ψ(ϱ)−ψ(s))ϑ−1ψ′(s)ϑ(ϱ)ds,ϱ∈(siϱi+1],i=1,2,…,n. |
Moreover, from (33) and (34) we get for ϱ∈(0,ϱ1]
||z(ϱ)−y∗(ϱ)||≤1Γ(ϑ)∫ϱ0(ψ(ϱ)−ψ(s))ϑ−1ψ′(s)||ϑ(ϱ)||ds≤ϵ(ψ(b)−ψ(0))ϑΓ(ϑ+1). |
Ifϱ∈(ϱi si], then||¯z(ϱ)−y(ϱ)||≤ϵ.If ϱ∈(si ϱi+1], then
||z(ϱ)−y∗(ϱ)||≤||(ψ(ϱ)−ψ(si))μ−1Γ(μ)gi(si,y(ϱ−i))(ψ(ϱ)−ψ(si))μ−1Γ(μ)y(si)||+1Γ(ϑ)∫ϱ0(ψ(ϱ)−ψ(s))ϑ−1ψ′(s)||ϑ(ϱ)||ds≤(ψ(ϱ)−ψ(si))μ−1Γ(μ)||gi(si,y(ϱ−i))−y(si)||+ϵ(ψ(b)−ψ(0))ϑΓ(ϑ+1)≤ϵ(ψ(b)−ψ(0))μ−1Γ(μ)+ϵ(ψ(b)−ψ(0))ϑΓ(ϑ+1)=ϵ[(ψ(b)−ψ(0))μ−1Γ(μ)+(ψ(b)−ψ(0))ϑΓ(ϑ+1)+1] |
As a consequence,
supϱ∈J||y∗(ϱ)−z(ϱ)||≤ϵ[(ψ(b)−ψ(0))μ−1Γ(μ)+(ψ(b)−ψ(0))ϑΓ(ϑ+1)+1]. | (35) |
On the other hand, from the facts that Φ is ς- contraction, ς is strictly comparison, ς(cϱ)≤c ς(ϱ)(c>1) for everyϱ∈[0,∞) and ϱ=0 is a point of uniform convergence for the series ∞∑n=1ςn(ϱ), it yields from Theorem3.1(ii)[41], that Φ is ρ−multivalued weakly Picard operator, where ρ(ϱ)=ϱ+∞∑n=1 ςn(ϱ). Then, the function Φ∞:Graph(Φ)→Fix(Φ) is well defined and
||y∗−Φ∞(y∗,z)||H≤ρ(||y∗−z||H). | (36) |
Put x=Φ∞(y∗,z). So, x∈Φ(x) and from (35) and (36) we get
dH(y,¯x)≤||y∗−Φ∞(y∗,z)||H≤ρ(||y∗−z||H)=ρ(ϵ[(ψ(b)−ψ(0))μ−1Γ(μ)+(ψ(b)−ψ(0))ϑΓ(ϑ+1)+1])=θ(ϵ), |
where θ(ϱ)=ρ(ϱ[(ψ(b)−ψ(0))μ−1Γ(μ)+(ψ(b)−ψ(0))ϑΓ(ϑ+1)+1]). Consequently, Problem (1) is ψ− generalized Ulam-Hyers stable.
In this section we give examples to clarify the possibility of applicability our assumptions.
EXAMPLE 1 Let E be a Hilbert space, K a non-empty convex compact subset of E,r=1,s0=0,ϱ1=1,s1=2,ϱ2=3, ϑ=12,ν=14and μ=ϑ+ν−ϑν=58. Let F:J×Θ→Pck(E) be defined as follows:
F(ϱ,z)={(ψ(ϱ)−ψ(0))38||z(0)||K, ϱ∈[0,1],||z(0)||K,ϱ∈(1,2],(ψ(ϱ)−ψ(2))38||z(0)||K,ϱ∈(2,3]. | (37) |
Clearly, for almost every ϱ∈J,z⟶F(ϱ,τ(ϱ)z) is upper semicontinuous. Set λ=sup{||x||:x∈K}. Then, for any (ϱ,z)∈J×Θ and any y∈F(ϱ,z) we have
||y(ϱ)||≤{(1+(ψ(ϱ)−ψ(0))38||z(0)||)λ, ϱ∈[0,1],||z(0)||λ,ϱ∈(1,2],(1+(ψ(ϱ)−ψ(2))38||z(0)||)λ,ϱ∈(2,3]. |
Then condition (F2)of Theorem1 is verified with φ(ϱ)=λ;ϱ∈J. Moreover, if f∈S1F(,.z);z∈Θ, then
f(ϱ)={(ψ(ϱ)−ψ(0))38||z(0)||y0, ϱ∈[0,1],y1,ϱ∈(1,2],(ψ(ϱ)−ψ(2))38y2,ϱ∈(2,3]. | (38) |
We define
hx1,x2,x3(ϱ)={Γ(118)Γ(108)(ψ(ϱ)−ψ(0))28x1, ϱ∈[0,1],x2,ϱ∈(1,2],Γ(118)Γ(108)(ψ(ϱ)−ψ(2))28x3,ϱ∈(2,3]. | (39) |
Obviously hx1,x2,x3∈PC1−μ,ψ(J,E) and from Lemma 2.2, we get Iν(1−ϑ),ψs+ihx1,x2,x3(ϱ)=f(ϱ);i=0,1.
Now, let D be a bounded subset of H, z1, z2∈D, ϱ∈(0,1], x∈F(ϱ.τ(ϱ)¯z1) and y∈F(ϱ.τ(ϱ)¯z2). Then x=(ψ(ϱ)−ψ(0))38||z1(ϱ)||x∗ and y=(ψ(ϱ)−ψ(0))38||z2(ϱ)||y∗, where x∗,y∗∈K. It follows that
||x−y||≤(ψ(ϱ)−ψ(0))38||x∗||z1(ϱ)||−||z2(ϱ)||y∗||, |
As a result,
infy∈F(ϱ.τ(ϱ)¯z2)||x−y||≤(ψ(ϱ)−ψ(0))38||x∗||z1(ϱ)||−||z2(ϱ)||x∗||=||x∗||(ψ(ϱ)−ψ(0))38|||z1(ϱ)||−||z2(ϱ)||≤λ((ψ(ϱ)−ψ(0))38||z1(ϱ)−z2(ϱ)||, |
which means that
h(F(ϱ.τ(ϱ)¯z1),F(ϱ.τ(ϱ)¯z2)≤λ((ψ(ϱ)−ψ(0))38||z1(ϱ)−z2(ϱ)||. | (40) |
Then, for ϱ∈(0,ϱ1]
χE(F(ϱ,{τ(ϱ)¯x:x∈D))≤λ(ψ(ϱ)−ψ(0))1−μς(ϱ)χE{x(ϱ):x∈D}≤λ(ψ(ϱ)−ψ(0))1−μς(ϱ)supθ∈[−r,0]χE{(τ(ϱ))¯x(θ):x∈D}, |
Similarly, one can show that if for ϱ∈(s2,ϱ3]
χE(F(ϱ,D))≤λ(ψ(ϱ)−ψ(s2))1−μsupθ∈[−r,0]χE{(τ(ϱ))¯x(θ):x∈D} |
and consequently, by choosing λ small enough such that the relation (7) becomes realized with ς(ϱ)=λ;ϱ∈J. Next, let g1:[1,2]×E→E such that
g(ϱ,x)=ρ(ψ(1)−ψ(0))38ϱ(x), | (41) |
where, ρ is a positive number andϱ is a linear, bounded and compact operator on E. So, condition (H1) is satisfied. As a consequence, from Theorem (1), the problem (1) has a solution where F and g are given by (47) and (50) and Ψ∈Θ provided that
3λκ1Γ(12)η+ρ+ρΓ(58)<1, | (42) |
where, η=δq(ψ(b)−ψ(0))ϑ−q(ϑ−q1−q)1−q, δ=maxs∈Jψ′(s), q=14 and κ1=(ψ(b)−ψ(0))1−μ. By choosing ρ small enough relation (42) becomes realized.
EXAMPLE 2 Let E, K, r, s0, ϱ1, s1, ϱ2, ϑ, ν, μ, F:J×Θ→Pck(E) and g1:[1,2]×E→E be as in EXAMPLE 1. Then assumption (F4) is satisfied. Moreover, in view of (40) for any z1,z1∈H
h(F(ϱ.τ(ϱ)¯z1),F(ϱ.τ(ϱ)¯z2)≤λ((ψ(ϱ)−ψ(0))38||z1(ϱ)−z2(ϱ)||≤λ||z1−z2||H=σ(ϱ)ς(||z1−z2||H), |
where, σ(ϱ)=2λ and ς(ϱ)=ϱ2;ϱ∈J. Observe that ς is strictly comparison and ς(cϱ)=c ς(ϱ), for every c>0, every ϱ∈[0,∞) and ϱ=0 is a point of uniform convergence for the series ∞∑n=1ςn(ϱ), and hence assumption (F5)is satisfied. Now, for any ϱ∈[1,2] and any x,y∈E
||g1(ϱ,x)−g1(ϱ,y)||≤ρ(ψ(1)−ψ(0))38||ϱ(x)−ϱ(y)||≤ρ||ϱ||(ψ(1)−ψ(0))38||x−y||E=ξς(||x−y||E), |
where ξ=2ρ||ϱ||(ψ(1)−ψ(0))38. It follows that (H2) is satisfied. By applying Theorem 2, Problem (1) has a solution and it is ψ−generalized Ulam-Hyers stable, where F and g are given by (37) and (40) and Ψ∈Θ provided that
2ρ||ϱ||(ψ(1)−ψ(0))38Γ(μ)+2ρ||ϱ||(ψ(1)−ψ(0))38Γ(μ)+η||σ||L1q(I,R+)<1. |
By choosing ρ small enough this inequality becomes realized.
A relation between a solution of the considered problem and the corresponding fractional integral equation is given, then two existence results of solutions for an ψ-Hilfer fractional non-instantaneous impulsive differential inclusion in the presence of delay in an infinite dimensional Banah spaces are obtained. Moreover, by using the multivalued weakly Picard operator theory, the stability of solutions for the considered problem in the sense of generalized Ulam-Hyers is studied. This work generalizes many recent results in the literature, for example [6,10,11,17,19,31]. Moreover, our technique can be used to study the existence and Ulam-Hyers stability of solutions or mild solutions for the problems considered in [3,4,15,16,20,21,22] to the case when, there are impulses and delay on the system, the right hand side is a multi-valued function and involving ψ−Hilfer fractional derivative. There are many directions for future work, for example: Generalize the obtained results in [3,4,15,16] when, the considered problems in these works involving ψ-Hilfer fractional derivative.
The authors acknowledge the Deanship of Scientific Research at King Faisal University for the financial support under Nasher Track (Grant No.186042).
The authors declare that they have no conflicts of interest.
All authors contributed equality and read and approved the final version of the manuscript.
[1] |
Y. Jia, S. Wang, J. Peng, Y. Gao, D. Hu, X. Zhao, Evaluation of pavement rutting based on driving safety of vehicles, Int. J. Pavement Res. Technol., 15 (2022), 457–469. https://doi.org/10.1007/s42947-021-00032-2 doi: 10.1007/s42947-021-00032-2
![]() |
[2] |
J. S. Tong, T. Ma, K. R. Shen, H. Y. Zhang, S. P. Wu, A criterion of asphalt pavement rutting based on the thermal-visco-elastic-plastic model, Int. J. Pavement Eng., 23 (2020), 4. https://doi.org/10.1080/10298436.2020.1792470 doi: 10.1080/10298436.2020.1792470
![]() |
[3] |
Y. Luan, Y. Ma, T. Ma, C. Wang, F. Xia, Research on the effects of asphalt performance on rutting resistance and its correlation with rutting performance indicators, Constr. Build. Mater., 400 (2023), 132773. https://doi.org/10.1016/j.conbuildmat.2023.132773 doi: 10.1016/j.conbuildmat.2023.132773
![]() |
[4] |
R. S. McDaniel, H. U. Bahia, Field evaluation of asphalt additives to control rutting and cracking, Transp. Res. Record, 1829 (2003), 47–54. https://doi.org/10.3141/1829-07 doi: 10.3141/1829-07
![]() |
[5] |
A. K. Ghuzlan, W. B. Al-Mistarehi, S. A. Al-Momani, Rutting performance of asphalt mixtures with gradations designed using Bailey and conventional Superpave methods, Constr. Build. Mater., 261 (2020), 119941. https://doi.org/10.1016/j.conbuildmat.2020.119941 doi: 10.1016/j.conbuildmat.2020.119941
![]() |
[6] |
J. Yang, J. Cui, J. Wan, X. Shi, Strategy of improving rutting resistance based on contribution rate of different layers, J. Southeast Univ. (Nat. Sci. Ed.), 37 (2007), 350–354. https://doi.org/10.3969/j.issn.1001-0505.2007.02.033 doi: 10.3969/j.issn.1001-0505.2007.02.033
![]() |
[7] | Y. L. Zhang, N. Li, Research on applicability of anti-rutting agent LY to different asphalt mixtures, in IOP Conference Series: Earth and Environmental Science, 587 (2020), 1. https://doi.org/10.1088/1755-1315/587/1/012029 |
[8] | W. J. Wang, Design of Key Parameters and Study on the Performance of Bituminous Mixtures Mixed with Anti-Rutting Additive, M.S. thesis, Chang'an University, 2011. |
[9] |
Z. Q. Zhang, Y. F. Luo, K. Zhang, Research on influence of different anti-rutting additives on high-temperature performance of asphalt mixture, J. Railw. Sci. Eng., 14 (2017), 964–970. https://doi.org/10.19713/j.cnki.43-1423/u.2017.05.012 doi: 10.19713/j.cnki.43-1423/u.2017.05.012
![]() |
[10] |
Y. S. Guan, H. Li, Z. X. Zhang, X. Y. Li, A new high rutting resistance asphalt mixture for intersection maintenance under heavy traffic, J. Test. Eval., 40 (2012), 1156–1161. https://doi.org/10.1520/JTE20120126 doi: 10.1520/JTE20120126
![]() |
[11] |
B. Wang, X. D. Hu, C. F. Cao, J. M. Wan, W. X. Gan, Z. W. Chen, et al., Development and characterization of permeative anti-rutting agent for asphalt mixture enhancement, Constr. Build. Mater., 364 (2023), 129937. https://doi.org/10.1016/j.conbuildmat.2022.129937 doi: 10.1016/j.conbuildmat.2022.129937
![]() |
[12] |
S. Han, W. Y. Ren, J. X. Jia, H. Yang, D. Y. Niu, Study on road performances of asphalt mixtures with granulated polymer anti-rutting additive, Int. J. Pavement Eng., 21 (2020), 257–265. https://doi.org/10.1080/10298436.2018.1461868 doi: 10.1080/10298436.2018.1461868
![]() |
[13] |
F. P. Xiao, D. F. Ma, J. Y. Wang, D. G. Cai, L. W. Lou, J. Yuan, Impacts of high modulus agent and anti-rutting agent on performances of airfield asphalt pavement, Constr. Build. Mater., 204 (2019), 1–9. https://doi.org/10.1016/j.conbuildmat.2019.01.138 doi: 10.1016/j.conbuildmat.2019.01.138
![]() |
[14] | JTG E20-2011, Standard Test Methods of Bitumen and Bituminous Mixtures for Highway Engineering, Ministry of Transport of the People's Republic of China, 2011. |
[15] | JTG E42-2005, Test Methods of Aggregate for Highway Engineering, Ministry of Transport of the People's Republic of China, 2005. |
[16] |
Z. Y. Wei, Y. S. Jia, S. Q. Wang, Z. R. Li, Y. S. Li, X. M. Wang, et al., Utilization of iron ore tailing as an alternative mineral filler in asphalt mastic: High-temperature performance and environmental aspects, J. Cleaner Prod., 335 (2022), 130318. https://doi.org/10.1016/j.jclepro.2021.130318 doi: 10.1016/j.jclepro.2021.130318
![]() |
[17] |
M. D. I. Domingos, A. L. Faxina, Susceptibility of asphalt binders to rutting: Literature review, J. Mater. Civ. Eng., 28 (2016), 04015134. https://doi.org/10.1061/(ASCE)MT.1943-5533.0001364 doi: 10.1061/(ASCE)MT.1943-5533.0001364
![]() |
1. | 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, 2023, 8, 2473-6988, 3469, 10.3934/math.2023177 | |
2. | Zainab Alsheekhhussain, Ahmed Gamal Ibrahim, Rabie A. Ramadan, Existence of $ S $-asymptotically $ \omega $-periodic solutions for non-instantaneous impulsive semilinear differential equations and inclusions of fractional order $ 1 < \alpha < 2 $, 2023, 8, 2473-6988, 76, 10.3934/math.2023004 | |
3. | Ramasamy Arul, Panjayan Karthikeyan, Kulandhaivel Karthikeyan, Palanisamy Geetha, Ymnah Alruwaily, Lamya Almaghamsi, El-sayed El-hady, On Ψ-Hilfer Fractional Integro-Differential Equations with Non-Instantaneous Impulsive Conditions, 2022, 6, 2504-3110, 732, 10.3390/fractalfract6120732 | |
4. | M. Latha Maheswari, K. S. Keerthana Shri, Karthik Muthusamy, Existence results for coupled sequential ψ-Hilfer fractional impulsive BVPs: topological degree theory approach, 2024, 2024, 1687-2770, 10.1186/s13661-024-01901-y | |
5. | Zainab Alsheekhhussain, Ahmad Gamal Ibrahim, Mohammed Mossa Al-Sawalha, Yousef Jawarneh, The Existence of Solutions for w-Weighted ψ-Hilfer Fractional Differential Inclusions of Order μ ∈ (1, 2) with Non-Instantaneous Impulses in Banach Spaces, 2024, 8, 2504-3110, 144, 10.3390/fractalfract8030144 | |
6. | Feryal Aladsani, Ahmed Gamal Ibrahim, Existence and Stability of Solutions for p-Proportional ω-Weighted κ-Hilfer Fractional Differential Inclusions in the Presence of Non-Instantaneous Impulses in Banach Spaces, 2024, 8, 2504-3110, 475, 10.3390/fractalfract8080475 | |
7. | Zainab Alsheekhhussain, Ahmed Gamal Ibrahim, M. Mossa Al-Sawalha, Khudhayr A. Rashedi, Mild Solutions for w-Weighted, Φ-Hilfer, Non-Instantaneous, Impulsive, w-Weighted, Fractional, Semilinear Differential Inclusions of Order μ ∈ (1, 2) in Banach Spaces, 2024, 8, 2504-3110, 289, 10.3390/fractalfract8050289 | |
8. | F. M. Ismaael, On a New Class of Impulsive η-Hilfer Fractional Volterra-Fredholm Integro-Differential Equations, 2023, 17, 1823-8343, 691, 10.47836/mjms.17.4.10 |