Next Article in Journal
Experiments on the Electron Impact Excitation of Hydrogen Molecules Indicate the Presence of the Second Flavor of Hydrogen Atoms
Next Article in Special Issue
Existence and Stability Results for Fractional Hybrid q-Difference Equations with q-Integro-Initial Condition
Previous Article in Journal
Longevity, Aging and Cancer: Thermodynamics and Complexity
Previous Article in Special Issue
Simpson’s Type Inequalities for s-Convex Functions via a Generalized Proportional Fractional Integral
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Nonlocal Boundary Value Problems for (k,ψ)-Hilfer Fractional Differential Equations and Inclusions

1
Department of Mathematics, University of Ioannina, 451 10 Ioannina, Greece
2
Nonlinear Analysis and Applied Mathematics (NAAM)-Research Group, Department of Mathematics, Faculty of Science, King Abdulaziz University, P.O. Box 80203, Jeddah 21589, Saudi Arabia
3
Intelligent and Nonlinear Dynamic Innovations Research Center, Department of Mathematics, Faculty of Applied Science, King Mongkut’s University of Technology North Bangkok, Bangkok 10800, Thailand
*
Author to whom correspondence should be addressed.
Foundations 2022, 2(3), 681-696; https://doi.org/10.3390/foundations2030046
Submission received: 24 June 2022 / Revised: 13 August 2022 / Accepted: 16 August 2022 / Published: 19 August 2022
(This article belongs to the Special Issue Recent Advances in Fractional Differential Equations and Inclusions)

Abstract

:
In the present research, single and multi-valued (k,ψ)-Hilfer type fractional boundary value problems of order in (1,2] involving nonlocal integral boundary conditions were studied. In the single-valued case, the Banach and Krasnosel’skiĭ fixed point theorems as well as the Leray–Schauder nonlinear alternative were used to establish the existence and uniqueness results. In the multi-valued case, when the right-hand side of the inclusion has convex values, we established an existence result via the Leray–Schauder nonlinear alternative method for multi-valued maps, while the second existence result, dealing with the non-convex valued right-hand side of the inclusion, was obtained by applying Covitz-Nadler fixed point theorem for multi-valued contractions. The obtained theoretical results are well illustrated by the numerical examples provided.

1. Introduction

Fractional calculus, dealing with integral and differential operators of non-integer order, has found interesting applications in many engineering and scientific disciplines such as physics, chemistry, mathematical biology, mechanics, and so forth, see the monographs [1,2,3,4,5,6,7,8,9]. Usually, fractional derivatives are defined in terms of fractional integral operators with different forms of the kernel function. Examples include Riemann–Liouville, Caputo, Hadamard, Katugampola and Hilfer fractional derivatives. Certain forms of fractional operators contain a number of different fractional operators. For example, the generalized fractional derivative of Katugampola [10,11] includes both Riemann–Liouville and Hadamard fractional derivatives. The Hilfer fractional derivative operator [12] contains Riemann–Liouville as well as Caputo fractional derivative operators. Another fractional derivative operator unifying Caputo, Caputo–Hadamard and Caputo–Erdélyi–Kober fractional derivative operators is the ψ -fractional derivative operator [13]. The ( k , ψ ) -Hilfer fractional derivative operator introduced in [14] generalizes many of the well-known fractional derivative operators, see [15].
Initial and boundary value problems involving the ( k , ψ ) -Hilfer fractional derivative operator recently received considerable attention. In [14], an existence and uniqueness result for a ( k , ψ ) -Hilfer type fractional initial value problem was derived. The authors discussed the existence of solutions for ( k , ψ ) -Hilfer fractional differential equations and inclusions supplemented with nonlocal boundary conditions in [15].
Motivated by the work presented in [14,15], in the present paper, we studied the existence of solutions for a ( k , ψ ) -Hilfer type fractional differential equation of order in ( 1 , 2 ] , equipped with nonlocal ( k , ψ ) -Riemann–Liouville fractional integral boundary conditions. In precise terms, we investigated the following ( k , ψ ) -Hilfer type fractional boundary value problem:
k , H D α , β ; ψ z ( t ) = f ( t , z ( t ) ) , t ( a , b ] , z ( a ) = 0 , z ( b ) = λ z ( ξ ) + μ k I v , ψ z ( σ ) ,
where k , H D α , β ; ψ denotes the ( k , ψ ) -Hilfer type fractional derivative operator of order α , 1 < α < 2 , 0 β 1 , k > 0 , f C ( [ a , b ] × R , R ) , k I v , ψ is the ( k , ψ ) -Riemann–Liouville fractional integral of order v > 0 , λ , μ R , and a < ξ ,   σ < b .
The corresponding multi-valued analogue of ( k , ψ ) -Hilfer type boundary value problem (1) provided by
k , H D α , β ; ψ z ( t ) F ( t , z ( t ) ) , t ( a , b ] , z ( a ) = 0 , z ( b ) = λ z ( ξ ) + μ k I v , ψ z ( σ ) ,
was also studied. In (2), F : [ a , b ] × R P ( R ) represents a multi-valued map ( P ( R ) denotes a family of all nonempty subsets of R ), while the other quantities are the same as described in the problem (1).
We applied Krasnosel’skiĭ’s fixed point theorem and the Leray–Schauder nonlinear alternative to prove the existence results for the problem (1), while the uniqueness of solutions for the problem (1) was established via Banach’s fixed point theorem. The existence results for the multivalued problem (2) for convex-valued and non-convex-valued cases were, respectively, obtained by means of the Leray–Schauder nonlinear alternative for multi-valued maps and the Covitz–Nadler fixed point theorem for multi-valued contractions. Concerning the advantages of the methods employed in the present study over other existing methods, we mention that the tools of the fixed point theory provide a suitable platform to establish the existence theory for boundary value problems once the problem at hand is converted into a fixed point problem.
Here, we recall that Hilfer fractional differential equations find useful applications in real world problems such as filtration processes [16,17], advection–diffusion phenomena [18], glass forming materials [19], etc. On the other hand, the nonlocal integral boundary conditions have potential applications in physical problems such as diffusion processes [20], blood flow problems [21], bacteria self-organization models [22], etc. We anticipated that the modeling of physical phenomena based on the Hilfer fractional derivative would be improved by using the ( k , ψ ) -Hilfer fractional derivative. Further, the ( k , ψ ) -Hilfer type boundary value problems considered in this paper correspond to a variety of fractional boundary value problems for different choices of ψ , for details, see [15]. In fact, the results obtained for the problems (1) and (2) are not only new in the given configuration but also correspond to several special cases for an appropriate choice of the values of ψ and the parameters involved in the given problems. Hence, the work established in this paper enriches the existing literature on the class of ( k , ψ ) -Hilfer boundary value problems.
The remainder of our paper is arranged as follows. In Section 2, we recall some fundamental concepts related to the study of the proposed problems. Section 3 contains an auxiliary result that plays a key role in converting the given problems into equivalent fixed point problems. Section 4 is devoted to the derivation of the main results for the single-valued problem (1), while the existence results for the multi-valued problem (2) are established in Section 5. Illustrative numerical examples demonstrating the applicability of the obtained theoretical results are presented in Section 6. The paper concludes with some interesting observations.

2. Preliminaries

Let us begin this section by introducing some preliminary concepts of fractional calculus.
Definition 1
([23]).Let k , α R + . The k-Riemann–Liouville fractional derivative of order α for the function h L 1 ( [ a , b ] , R ) is obtained by
k I a + α h ( t ) = 1 k Γ k ( α ) a t ( t u ) α k 1 h ( u ) d u ,
where Γ k is the k-Gamma function defined by [24]
Γ k ( z ) = 0 s z 1 e s k k d s , z C , ( z ) > 0 ,
which satisfies the following properties:
lim k 1 Γ k ( θ ) = Γ ( θ ) , Γ k ( θ ) = k θ k 1 Γ θ k and Γ k ( θ + k ) = θ Γ k ( θ ) .
Definition 2
([25]).The k-Riemann-Liouville fractional derivative of order α for the function h L 1 ( [ a , b ] , R ) with k , α R + is defined by
k , R L D a + α h ( t ) = k d d t n k I a + n k α h ( t ) , n = α k ,
where α k represents the ceiling function of α k .
Definition 3.
Let ψ : [ a , b ] R be an increasing function with ψ ( θ ) 0 for all θ [ a , b ] . Then, the ψ-Riemann–Liouville and ( k , ψ ) -Riemann–Liouville fractional integrals of order α for the function h L 1 ( [ a , b ] , R ) , respectively defined in [2,26], are given by
I α ; ψ h ( t ) = 1 Γ k ( α ) a t ψ ( u ) ( ψ ( t ) ψ ( u ) ) α 1 h ( u ) d u ,
k I a + α ; ψ h ( t ) = 1 k Γ k ( α ) a t ψ ( u ) ( ψ ( t ) ψ ( u ) ) α k 1 h ( u ) d u .
Definition 4
([27]).The ψ-Hilfer fractional derivative for the function h C ( [ a , b ] , R ) of order α ( n 1 , n ] , n N and type β [ 0 , 1 ] is defined by
H D α , β ; ψ h ( t ) = I a + β ( n α ) ; ψ 1 ψ ( t ) d d t n I a + ( 1 β ) ( n α ) ; ψ h ( t ) ,
where ψ : [ a , b ] R is an increasing function such that ψ C n ( [ a , b ] , R ) and ψ ( θ ) 0 ,   θ [ a , b ] . On the other hand, the ( k , ψ ) -Hilfer fractional derivative of order α and type β for the function h C n ( [ a , b ] , R ) defined in [14] is
k , H D α , β ; ψ h ( t ) = k I a + β ( n k α ) ; ψ k ψ ( t ) d d t n k I a + ( 1 β ) ( n k α ) ; ψ h ( t ) , n = α k .
Remark 1.
For β [ 0 , 1 ] and n 1 < α k n , we have n 1 < θ k k n , θ k = α + β ( n k α ) . Further, one can notice that
k , H D α , β ; ψ h ( t ) = k I a + θ k α ; ψ k ψ ( t ) d d t n k I a + n k θ k ; ψ h ( t ) = k I a + θ k α ; ψ k , R L D θ k ; ψ h ( t ) .
Now, we recall some useful results involving the ( k , ψ ) -Hilfer fractional derivative and integral operators.
Lemma 1
([14]).Let h C n ( [ a , b ] , R ) and k I a + n k μ ; ψ h C n ( [ a , b ] , R ) with μ , k R + = ( 0 , ) and n = μ k . Then,
k I μ ; ψ k , R L D μ ; ψ h ( t ) = h ( t ) j = 1 n ( ψ ( t ) ψ ( a ) ) μ k j Γ k ( μ j k + k ) k ψ ( t ) d d t n j k I a + n k μ ; ψ h ( t ) z = a .
Lemma 2
([14]).Let θ k = α + β ( k α ) with α , k R + = ( 0 , ) , α < k , and β [ 0 , 1 ] . Then,
k I θ k ; ψ k , R L D θ k ; ψ h ( t ) = k I α ; ψ k , H D α , β ; ψ h ( t ) , h C n ( [ a , b ] , R ) .

3. An Auxiliary Result

In this section, we provide an auxiliary result, which helps us in transforming the nonlinear ( k , ψ ) -Hilfer type fractional boundary value problem (1) into a fixed point problem. The following lemma concerns a linear variant of the ( k , ψ ) -Hilfer type fractional boundary value problem (1).
Lemma 3.
For g C ( a , b ) L 1 ( a , b ) and Λ 0 , the function z C ( [ a , b ] , R ) is a solution of the following problem
k , H D α , β ; ψ z ( t ) = g ( t ) , k > 0 , 1 < α 2 , β [ 0 , 1 ] , t ( a , b ] , z ( a ) = 0 , z ( b ) = λ z ( ξ ) + μ k I v , ψ z ( σ ) ,
if and only if
z ( t ) = k I α ; ψ g ( t ) + ( ψ ( t ) ψ ( a ) ) θ k k 1 Λ Γ k ( θ k ) λ k I α ; ψ g ( ξ ) + μ k I v + α , ψ g ( σ ) k I α ; ψ g ( b ) ,
where θ k = α + β ( 2 k α ) and
Λ : = ( ψ ( b ) ψ ( a ) ) θ k k 1 Γ k ( θ k ) λ ( ψ ( ξ ) ψ ( a ) ) θ k k 1 Γ k ( θ k ) μ ( ψ ( σ ) ψ ( a ) ) θ k + v k 1 Γ k ( θ k + v ) .
Proof. 
Assume that z is a solution of the problem (3). Operating on both sides of equation in (3), the fractional integral k I α ; ψ and using Lemmas 1 and 2, we obtain
k I α ; ψ k , H D α , β ; ψ z ( t ) = k I θ k ; ψ k , R L D θ k ; ψ z ( t ) = z ( t ) ( ψ ( t ) ψ ( a ) ) θ k k 1 Γ k ( θ k ) k ψ ( t ) d d t k I 2 k θ k ; ψ z ( t ) w = a ( ψ ( t ) ψ ( a ) ) θ k k 2 Γ k ( θ k k ) k I 2 k θ k ; ψ z ( t ) w = a .
Consequently
z ( t ) = k I α ; ψ g ( t ) + c 0 ( ψ ( t ) ψ ( a ) ) θ k k 1 Γ k ( θ k ) + c 1 ( ψ ( t ) ψ ( a ) ) θ k k 2 Γ k ( θ k k ) ,
where
c 0 = k ψ ( t ) d d t k I 2 k θ k ; ψ z ( t ) w = a , c 1 = k I 2 k θ k ; ψ z ( t ) w = a .
By the condition z ( a ) = 0 , we obtain c 1 = 0 , since θ k k 2 < 0 as shown by Remark 1. Now, using the nonlocal condition: z ( b ) = λ z ( ξ ) + μ k I v , ψ z ( σ ) and the following result from [14]:
k I v , ψ ( ψ ( t ) ψ ( a ) ) θ k k 1 = Γ k ( θ k ) Γ k ( θ k + v ) ( ψ ( t ) ψ ( a ) ) θ k + v k 1 ,
we find that
c 0 = 1 Λ λ k I α ; ψ g ( ξ ) + μ k I v + α , ψ g ( σ ) k I α ; ψ g ( b ) .
Substituting the above value of c 0 and c 1 = 0 in (6), we obtain the solution (4). By carrying out direct computation we can easily establish the converse of this lemma. The proof is completed. □

4. The Single-Valued Problem

Let C ( [ a , b ] , R ) denotes the Banach space of all continuous functions from [ a , b ] to R endowed with the norm z = sup t [ a , b ] | z ( t ) | . By using Lemma 3, an operator A : C ( [ a , b ] , R ) C ( [ a , b ] , R ) can be defined as
( A z ) ( t ) = ( ψ ( t ) ψ ( a ) ) θ k k 1 Λ Γ k ( θ k ) λ k I α ; ψ f ( ξ i , z ( ξ i ) ) + μ k I v + α , ψ f ( σ , z ( σ ) ) k I α ; ψ f ( b , z ( b ) ) + k I α ; ψ f ( t , z ( t ) ) , t [ a , b ] .
Notice that the fixed points of the operator A will be solutions of the nonlocal ( k , ψ ) -Hilfer type fractional boundary value problem (1).
For computational convenience, we set the following notation:
G = ( ψ ( b ) ψ ( a ) ) α k Γ k ( α + k ) + ( ψ ( b ) ψ ( a ) ) θ k k 1 | Λ | Γ k ( θ k ) [ | λ | ( ψ ( ξ ) ψ ( a ) ) α k Γ k ( α + k ) + | μ | ( ψ ( b ) ψ ( a ) ) α + v k Γ k ( α + v + k ) + ( ψ ( b ) ψ ( a ) ) α k Γ k ( α + k ) ] .

4.1. Existence and Uniqueness Result

In this subsection, we make use of Banach’s fixed point theorem [28] to prove a uniqueness result for the ( k , ψ ) -Hilfer type fractional boundary value problem (1).
Theorem 1.
Suppose that
(H1) | f ( t , z ) f ( t , y ) | L | z y | , L > 0 f o r e a c h t [ a , b ] a n d z , y R .
Then, there exists a unique solution for the problem (1) on [ a , b ] provided that
L G < 1 ,
where G is defined by (9)
Proof. 
Transform the problem (1) into a fixed point problem A z = z , where the operator A : C ( [ a , b ] , R ) C ( [ a , b ] , R ) is defined in (8). We verify the hypothesis of Banach’s fixed point theorem to deduce that the operator A has a unique fixed point.
Let us first show that A B r B r , where B r = { z C ( [ a , b ] , R ) : z r } with
r Q G 1 L G , sup t [ a , b ] | f ( t , 0 ) | = Q < .
By ( H 1 ) , we have
| f ( t , z ( t ) ) | | f ( t , z ( t ) ) f ( t , 0 ) | + | f ( t , 0 ) | L | z ( t ) | + Q L z + Q L r + Q .
For any z B r , we have
| ( A z ) ( t ) | sup t [ a , b ] { ( ψ ( t ) ψ ( a ) ) θ k k 1 | Λ | Γ k ( θ k ) [ | λ | k I α ; ψ | f ( ξ , z ( ξ ) ) | + | μ | k I v + α , ψ | f ( σ , z ( σ ) ) | + k I α ; ψ | f ( b , z ( b ) ) | ] + k I α ; ψ | f ( t , z ( t ) ) | } k I α ; ψ ( | f ( t , z ( t ) ) f ( t , 0 ) | + | f ( t , 0 ) | ) + ( ψ ( b ) ψ ( a ) ) θ k k 1 | Λ | Γ k ( θ k ) ( | λ | k I α ; ψ ( | f ( ξ , z ( ξ ) ) f ( ξ , 0 ) | + | f ( ξ , 0 ) | ) + | μ | k I v + α , ψ ( | f ( σ , z ( σ ) ) f ( σ , 0 ) | + | f ( σ , 0 ) | ) + k I α ; ψ ( | f ( b , z ( b ) ) f ( b , 0 ) | + | f ( b , 0 ) | ) ) { ( ψ ( b ) ψ ( a ) ) α k Γ k ( α + k ) + ( ψ ( b ) ψ ( a ) ) θ k k 1 | Λ | Γ k ( θ k ) [ | λ | ( ψ ( ξ ) ψ ( a ) ) α k Γ k ( α + k ) + | μ | ( ψ ( b ) ψ ( a ) ) α + v k Γ k ( α + v + k ) + ( ψ ( b ) ψ ( a ) ) α k Γ k ( α + k ) ] } ( L z + Q ) ( L r + Q ) G r ,
where (11) has been applied. Consequently, A z r , which means that A B r B r .
In order to show that A is a contraction, let z , y C ( [ a , b ] , R ) . Then, for t [ a , b ] , we obtain
| ( A z ) ( t ) ( A y ) ( t ) | k I α ; ψ | f ( t , z ( t ) ) f ( t , y ( t ) ) | + ( ψ ( b ) ψ ( a ) ) θ k k 1 | Λ | Γ k ( θ k ) ( | λ | k I α ; ψ | f ( ξ , z ( ξ ) ) f ( ξ , y ( ξ ) | + | μ | k I α + v ; ψ | f ( σ , z ( σ ) ) f ( σ , y ( σ ) | + k I α ; ψ ( | f ( b , z ( b ) ) f ( b , y ( b ) ) | ) { ( ψ ( b ) ψ ( a ) ) α k Γ k ( α + k ) + ( ψ ( b ) ψ ( a ) ) θ k k 1 | Λ | Γ k ( θ k ) [ | λ | ( ψ ( ξ ) ψ ( a ) ) α k Γ k ( α + k ) + | μ | ( ψ ( b ) ψ ( a ) ) α + v k Γ k ( α + v + k ) + ( ψ ( b ) ψ ( a ) ) α k Γ k ( α + k ) ] } L x y = L G x y .
Thus, A x A y L G x y , which shows that the operator A is a contraction in view of the condition (10). Hence, by Banach’s fixed point theorem, the operator A has a unique fixed point, which is indeed a unique solution of the problem (1) on [ a , b ] . This finishes the proof. □

4.2. Existence Results

In this subsection, we present two existence results for problem (1), which are proved with the aid of Krasnosel’skiĭ’s fixed point theorem [29] and nonlinear alternative of the Leray–Schauder type [30].
Theorem 2.
Let the continuous function f : [ a , b ] × R R satisfy assumption ( H 1 ) . In addition, we suppose that
(H2) | f ( t , z ) | φ ( t ) , ( t , z ) [ a , b ] × R , and φ C ( [ a , b ] , R + ) .
Then, there exists at least one solution for the problem (1) on [ a , b ] , if G 1 L < 1 , where
G 1 : = ( ψ ( b ) ψ ( a ) ) θ k k 1 | Λ | Γ k ( θ k ) | λ | ( ψ ( ξ ) ψ ( a ) ) α k Γ k ( α + k ) + | μ | ( ψ ( b ) ψ ( a ) ) α + v k Γ k ( α + v + k ) + ( ψ ( b ) ψ ( a ) ) α k Γ k ( α + k ) .
.
Proof. 
Set sup t [ a , b ] φ ( t ) = φ and B ρ = { z C ( [ a , b ] , R ) : z ρ } , with ρ φ G . We define on B ρ two operators A 1 , A 2 by
A 1 z ( t ) = k I α ; ψ f ( t , z ( t ) ) , t [ a , b ] ,
( A 2 z ) ( t ) = ( ψ ( t ) ψ ( a ) ) θ k k 1 Λ Γ k ( θ k ) [ λ k I α ; ψ f ( ξ , z ( ξ ) ) + μ k I v + α , ψ f ( σ , z ( σ ) ) k I α ; ψ f ( b , z ( b ) ) ] , t [ a , b ] .
For any z , y B ρ , we have
| ( A 1 z ) ( t ) + ( A 2 y ) ( t ) | sup t [ a , b ] { ( ψ ( t ) ψ ( a ) ) θ k k 1 | Λ | Γ k ( θ k ) [ | λ | k I α ; ψ | f ( ξ , y ( ξ ) ) | + | μ | k I v + α , ψ | f ( σ , y ( σ ) ) | + k I α ; ψ | f ( b , y ( b ) ) | ] + k I α ; ψ | f ( t , z ( t ) ) | } { ( ψ ( b ) ψ ( a ) ) α k Γ k ( α + k ) + ( ψ ( b ) ψ ( a ) ) θ k k 1 | Λ | Γ k ( θ k ) [ | λ | ( ψ ( ξ ) ψ ( a ) ) α k Γ k ( α + k ) + | μ | ( ψ ( b ) ψ ( a ) ) α + v k Γ k ( α + v + k ) + ( ψ ( b ) ψ ( a ) ) α k Γ k ( α + k ) ] } φ = G φ ρ .
Therefore, ( A 1 z ) + ( A 2 y ) ρ , which shows that A 1 z + A 2 y B ρ . In the next step, by using (12), we can easily show that the operator A 2 is a contraction mapping.
Observe that the continuity of f implies that of the operator A 1 . Additionally, A 1 is uniformly bounded on B ρ as
A 1 z ( ψ ( b ) ψ ( a ) ) α k Γ k ( α + k ) φ .
Next, we establish equicontinuity of the operator A 1 . For t 1 , t 2 [ a , b ] with t 1 < t 2 , it is easy show that
| ( A 1 z ) ( t 2 ) ( A 1 z ) ( t 1 ) | φ Γ k ( α + k ) [ 2 ( ψ ( t 2 ) ψ ( t 1 ) ) α k + | ( ψ ( t 2 ) ψ ( a ) ) α k ( ψ ( t 1 ) ψ ( a ) ) α k | ] ,
which, in the limit t 2 t 1 0 , tends to zero independently of z B ρ . Thus, A 1 is equicontinuous. Consequently, using the Arzelá-Ascoli theorem, A 1 is completely continuous. Thus, the hypothesis of Krasnosel’skiĭ’s fixed point theorem is verified, which guarantees that the problem (1) has at least one solution on [ a , b ] . The proof is complete. □
Theorem 3.
Assume that
(H3)
a continuous nondecreasing function χ : [ 0 , ) ( 0 , ) and a function σ C ( [ a , b ] , R + ) satisfying | f ( t , z ) | σ ( t ) χ ( | z | ) , ( t , z ) [ a , b ] × R ;
(H4)
a constant K > 0 such that K χ ( K ) σ G > 1 .
Then, there exists at least one solution on [ a , b ] for the problem (1).
Proof. 
Let us first show that the operator A defined by (8) maps bounded sets into C ( [ a , b ] , R ) . For r > 0 , let B r = { z C ( [ a , b ] , R ) : z r } . Then, for t [ a , b ] , we obtain
| ( A z ) ( t ) | sup t [ a , b ] { ( ψ ( t ) ψ ( a ) ) θ k k 1 | Λ | Γ k ( θ k ) [ | λ | k I α ; ψ | f ( ξ , z ( ξ ) ) | + | μ | k I v + α , ψ | f ( σ , z ( σ ) ) | + k I α ; ψ | f ( b , z ( b ) ) | ] + k I α ; ψ | f ( t , z ( t ) ) | } { ( ψ ( b ) ψ ( a ) ) α k Γ k ( α + k ) + ( ψ ( b ) ψ ( a ) ) θ k k 1 | Λ | Γ k ( θ k ) [ | λ | ( ψ ( ξ ) ψ ( a ) ) α k Γ k ( α + k ) + | μ | ( ψ ( b ) ψ ( a ) ) α + v k Γ k ( α + v + k ) + ( ψ ( b ) ψ ( a ) ) α k Γ k ( α + k ) ] } σ χ ( z ) ,
which implies that
A z χ ( r ) σ G .
Next, it is shown that the operator A maps bounded sets into equicontinuous sets of C ( [ a , b ] , R ) . For t 1 , t 2 [ a , b ] with t 1 < t 2 and z B r , we obtain
| ( A z ) ( t 2 ) ( A z ) ( t 1 ) | 1 Γ k ( α ) | a t 1 ψ ( s ) [ ( ψ ( t 2 ) ψ ( s ) ) α k 1 ( ψ ( t 1 ) ψ ( s ) ) α k 1 ] f ( s , z ( s ) ) d s + t 1 t 2 ψ ( s ) ( ψ ( t 2 ) ψ ( s ) ) α k 1 f ( s , z ( s ) ) d s | + ( ψ ( t 2 ) ψ ( a ) ) θ k k 1 ( ψ ( t 1 ) ψ ( a ) ) θ k k 1 | Λ | Γ k ( θ k ) [ | λ | k I α ; ψ | f ( ξ , z ( ξ ) ) | + | μ | k I v + α , ψ | f ( σ , z ( σ ) ) | + k I α ; ψ | f ( b , z ( b ) ) | ] σ χ ( r ) Γ k ( α + k ) [ 2 ( ψ ( t 2 ) ψ ( t 1 ) ) α k + | ( ψ ( t 2 ) ψ ( a ) ) α k ( ψ ( t 1 ) ψ ( a ) ) α k | ] , + ( ψ ( t 2 ) ψ ( a ) ) θ k k 1 ( ψ ( t 1 ) ψ ( a ) ) θ k k 1 | Λ | Γ k ( θ k ) [ | λ | ( ψ ( ξ ) ψ ( a ) ) α k Γ k ( α + k ) + | μ | ( ψ ( b ) ψ ( a ) ) α + v k Γ k ( α + v + k ) + ( ψ ( b ) ψ ( a ) ) α k Γ k ( α + k ) ] σ χ ( r ) ,
which, as t 2 t 1 0 , tends to zero independently of z B r . As a consequence, we deduce by the Arzelá-Ascoli theorem that the operator A : C ( [ a , b ] , R ) C ( [ a , b ] , R ) is completely continuous.
Finally, we prove that the set of all solutions to the equation z = λ A z is bounded for λ ( 0 , 1 ) .
As in the first step, one can obtain
z χ ( z ) σ G 1 .
By ( H 4 ) , we can find K > 0 such that z K . Consider the set
M = { z C ( [ a , b ] , R ) : z < K } .
Observe that the operator A : M ¯ C ( [ a , b ] , R ) is continuous and completely continuous. Thus, there does not exist any z M satisfying z = λ A z for some λ ( 0 , 1 ) by the given choice of M . So, A has a fixed point z M ¯ by the application of the nonlinear alternative of the Leray–Schauder type, which means that there exists at least one solution for the problem (1) on [ 0 , 1 ] . This finishes the proof. □

5. The Multi-Valued Problem

For each z C ( [ a , b ] , R ) , we define the set of selections of F as
S F , z : = { f L 1 ( [ a , b ] , R ) : f ( t ) F ( t , z ( t ) ) o n [ a , b ] } .
Definition 5.
A continuous function z is said to be a solution of the ( k , ψ ) -Hilfer type nonlocal integral fractional boundary value problem (2), if it satisfies the boundary conditions z ( a ) = 0 , z ( b ) = λ z ( ξ ) + μ k I v , ψ z ( σ ) , and there exists an integrable function f with f ( t ) F ( t , z ) for a.e. t [ a , b ] such that z satisfies the differential equation k , H D α , β ; ψ z ( t ) = f ( t ) on [ a , b ] .
Our first result for the multi-valued problem (2) is concerned with the case when the multi-valued map F has convex values, and relies on the nonlinear alternative of the Leray–Schauder type for multi-valued maps [30].
Theorem 4.
Suppose that:
( G 1 )
F : [ a , b ] × R P c p , c ( R ) is L 1 -Carathéodory, where P c p , c ( R ) = { R P ( R ) : R   i s   c o m p a c t   a n d   c o n v e x } ;
( G 2 )
a continuous nondecreasing function H : [ 0 , + ) ( 0 , + ) and a positive continuous real valued function q such that, ∀ ( t , z ) [ a , b ] × R ,
F ( t , z ) P : = sup { | f | : f F ( t , z ) } q ( t ) H ( z ) ;
( G 3 )
a constant K > 0 such that
K q H ( K ) G > 1 ,
where G is defined by (9).
Then, the multi-valued problem (2) has at least one solution on [ a , b ] .
Proof. 
We define an operator F : C ( [ a , b ] , R ) P ( C ( [ a , b ] , R ) ) by
F ( z ) = h C ( [ a , b ] , R ) : h ( t ) = ( ψ ( t ) ψ ( a ) ) θ k k 1 Λ Γ k ( θ k ) λ k I α ; ψ f ( ξ ) + μ k I α + v ; ψ f ( σ ) k I α ; ψ f ( b ) + k I α ; ψ f ( t ) , t [ a , b ] , f S F , z .
Notice that the fixed points of F are solutions to the problem (2).
We split the proof into several steps.
Step 1. F ( z ) is convex, for each z C ( [ a , b ] , R ,
Since S F , z is convex, this step is obvious, and so the proof is omitted.
Step 2. Bounded sets are mapped by F into bounded sets in C ( [ a , b ] , R ) .
Let B r = { z C ( [ a , b ] , R ) : z r } , r > 0 . Then, for each h F ( z ) , z B r , there exists f S F , z such that
h ( t ) = ( ψ ( t ) ψ ( a ) ) θ k k 1 Λ Γ k ( θ k ) λ k I α ; ψ f ( ξ ) + μ k I α + v ; ψ f ( σ ) k I α ; ψ f ( b ) + k I α ; ψ f ( t ) .
Further, for t [ a , b ] , we have
| h ( t ) | sup t [ a , b ] { ( ψ ( t ) ψ ( a ) ) θ k k 1 | Λ | Γ k ( θ k ) | λ | k I α ; ψ | f ( ξ ) | + | μ | k I α + v ; ψ | f ( σ ) | + k I α ; ψ | f ( b ) | + k I α ; ψ | f ( t ) | } { ( ψ ( b ) ψ ( a ) ) α k Γ k ( α + k ) + ( ψ ( b ) ψ ( a ) ) θ k k 1 | Λ | Γ k ( θ k ) [ | λ | ( ψ ( ξ ) ψ ( a ) ) α k Γ k ( α + k ) + | μ | ( ψ ( b ) ψ ( a ) ) α + v k Γ k ( α + v + k ) + ( ψ ( b ) ψ ( a ) ) α k Γ k ( α + k ) ] } q H ( z ) ,
which implies that
h H ( r ) q G .
Step 3. F maps bounded sets into equicontinuous sets of C ( [ a , b ] , R ) .
Let t 1 , t 2 [ a , b ] with t 1 < t 2 and z B r . Then, for each h F ( z ) , we find that
| h ( t 2 ) h ( t 1 ) | 1 Γ k ( α ) | a t 1 ψ ( s ) [ ( ψ ( t 2 ) ψ ( s ) ) α k 1 ( ψ ( t 1 ) ψ ( s ) ) α k 1 ] f ( s ) d s + t 1 t 2 ψ ( s ) ( ψ ( t 2 ) ψ ( s ) ) α k 1 f ( s ) d s | + ( ψ ( t 2 ) ψ ( a ) ) θ k k 1 ( ψ ( t 1 ) ψ ( a ) ) θ k k 1 | Λ | Γ k ( θ k ) [ | λ | k I α ; ψ | f ( ξ ) | + | μ | k I α + v ; ψ | f ( σ ) | + k I α ; ψ | f ( b ) | ] q H ( r ) Γ k ( α + k ) [ 2 ( ψ ( t 2 ) ψ ( t 1 ) ) α k + | ( ψ ( t 2 ) ψ ( a ) ) α k ( ψ ( t 1 ) ψ ( a ) ) α k | ] , + ( ψ ( t 2 ) ψ ( a ) ) θ k k 1 ( ψ ( t 1 ) ψ ( a ) ) θ k k 1 | Λ | Γ k ( θ k ) [ | λ | ( ψ ( ξ ) ψ ( a ) ) α k Γ k ( α + k ) + | μ | ( ψ ( b ) ψ ( a ) ) α + v k Γ k ( α + v + k ) + ( ψ ( b ) ψ ( a ) ) α k Γ k ( α + k ) ] q z ( r ) ,
which shows that | h ( t 2 ) h ( t 1 ) | 0 independently of z B r as t 2 t 1 0 . Therefore, F : C ( [ a , b ] , R ) P ( C ( [ a , b ] , R ) ) is completely continuous by the application of the Arzelá–Ascoli theorem.
In the next step, we show that F has a closed graph, which is equivalent to the fact that F is a upper semi-continuous multivalued map by Proposition 1.2 in [31].
Step 4. F has a closed graph.
Let z n z * , h n F ( z n ) and h n h * . Then, we need to show that h * F ( z * ) . Since h n F ( z n ) , there exists v n S F , z n such that, for each t [ a , b ] ,
h n ( t ) = ( ψ ( t ) ψ ( a ) ) θ k k 1 Λ Γ k ( θ k ) λ k I α ; ψ f n ( ξ ) + μ k I α + v ; ψ f n ( σ ) k I α ; ψ f n ( b ) + k I α ; ψ f n ( t ) .
Thus, we must show that there exists v * S F , z * such that, for each t [ a , b ] ,
h * ( t ) = ( ψ ( t ) ψ ( a ) ) θ k k 1 Λ Γ k ( θ k ) λ k I α ; ψ f * ( ξ ) + μ k I α + v ; ψ f * ( σ ) k I α ; ψ f * ( b ) + k I α ; ψ f * ( t ) .
Let us consider the linear operator Ω : L 1 ( [ a , b ] , R ) C ( [ a , b ] , R ) given by
f Ω ( f ) ( t ) = ( ψ ( t ) ψ ( a ) ) θ k k 1 Λ Γ k ( θ k ) λ k I α ; ψ f ( ξ ) + μ k I α + v ; ψ f ( σ ) k I α ; ψ f ( b ) + k I α ; ψ f ( t ) .
Observe that h n ( t ) h * ( t ) 0 , as n . Consequently, Ω S F is a closed graph operator, by a Lazota–Opial result [32]. Further, we obtain h n ( t ) Ω ( S F , z n ) . Since z n z * , we obtain
h * ( t ) = ( ψ ( t ) ψ ( a ) ) θ k k 1 Λ Γ k ( θ k ) λ k I α ; ψ f * ( ξ ) + μ k I α + v ; ψ f * ( σ ) k I α ; ψ f * ( b ) + k I α ; ψ f * ( t ) ,
for some v * S F , z * .
Step 5. We prove that there exists an open set U C ( [ a , b ] , R ) such that, for any ν ( 0 , 1 ) and all z U , we have z ν F ( z ) .
Assume that z ν F ( z ) for ν ( 0 , 1 ) . Then, there exists f L 1 ( [ a , b ] , R ) with f S F , z such that, for t [ a , b ] , we have
z ( t ) = ν ( ψ ( t ) ψ ( a ) ) θ k k 1 Λ Γ k ( θ k ) λ k I α ; ψ f ( ξ ) + μ k I α + v ; ψ f ( σ ) k I α ; ψ f ( b ) + ν k I α ; ψ f ( t ) .
Then, as in the second step, one can obtain
z q H ( z ) G 1 .
By ( H 3 ) , there exists K such that z K . Let us set
U = { z C ( [ a , b ] , R ) : z < K } .
From the preceding arguments, F : U ¯ P ( C ( [ a , b ] , R ) ) is a compact and upper semi-continuous multivalued map with convex closed values. By definition of U , there does not exist any z U such that z ν F ( z ) for some ν ( 0 , 1 ) . Hence, it follows by the nonlinear alternative of the Leray–Schauder type for multi-valued maps [30] that F has a fixed point z U ¯ , which is indeed a solution to the multi-valued problem (2). The proof is complete. □
Now, we apply the fixed point theorem for multivalued contractive maps suggested by Covitz and Nadler [33] to show that there exists a solution to the problem (2) when F is not necessarily a convex valued map.
Theorem 5.
Assume that
( A 1 )
F : [ a , b ] × R P c p ( R ) is such that f ( · , z ) : [ a , b ] P c p ( R ) is measurable for each z R , where P c p ( R ) = { V P ( R ) : V   i s   c o m p a c t } ;
( A 2 )
a function m C ( [ a , b ] , R + ) such that
H d ( F ( t , z ) , F ( t , z ¯ ) ) m ( t ) | z z ¯ | ,
with d ( 0 , F ( t , 0 ) ) m ( t ) for almost all t [ a , b ] and z , z ¯ R .
Then, the problem (2) has at least one solution on [ a , b ] , provided that
δ : = G m < 1 ,
where G is given by (9).
Proof. 
By the assumption ( A 1 ) , the set S F , z is nonempty for each z C ( [ a , b ] , R ) . Hence, by implementing Theorem III.6 [34], F has a measurable selection. We now prove that F ( z ) P c l ( C ( [ a , b ] , R ) ) for each z C ( [ a , b ] , R ) . Consider { z n } n 0 F ( z ) such that z n z ( n ) in C ( [ a , b ] , R ) . Then, we have z C ( [ a , b ] , R ) and there exists v n S F , z n such that, for each t [ a , b ] ,
z n ( t ) = ( ψ ( t ) ψ ( a ) ) θ k k 1 Λ Γ k ( θ k ) λ k I α ; ψ f n ( ξ ) + μ k I α + v ; ψ f n ( σ ) k I α ; ψ f n ( b ) + k I α ; ψ f n ( t ) .
Since F has compact values, there exists v n which converges to v in L 1 ( [ a , b ] , R ) . Hence, v S F , z and for each t [ a , b ] , we have
z n ( t ) z ( t ) = ( ψ ( t ) ψ ( a ) ) θ k k 1 Λ Γ k ( θ k ) λ k I α ; ψ f ( ξ ) + μ k I α + v ; ψ f ( σ ) k I α ; ψ f ( b ) + k I α ; ψ f ( t ) .
Thus, z F ( z ) .
Next, we show that
H d ( F ( z ) , F ( z ¯ ) ) δ z z ¯ , δ < 1 , for each z , z ¯ C 2 ( [ a , b ] , R ) .
Let z , z ¯ C 2 ( [ a , b ] , R ) and h 1 F ( x ) . Then ∃ v 1 ( t ) F ( t , z ( t ) ) such that, for each t [ a , b ] ,
h 1 ( t ) = ( ψ ( t ) ψ ( a ) ) θ k k 1 Λ Γ k ( θ k ) λ k I α ; ψ f 1 ( ξ ) + μ k I α + v ; ψ f 1 ( σ ) k I α ; ψ f 1 ( b ) + k I α ; ψ f 1 ( t ) .
Using ( A 2 ) , we obtain
H d ( F ( t , z ) , F ( t , z ¯ ) ) m ( t ) | z ( t ) z ¯ ( t ) | .
So, ∃ w F ( t , x ¯ ( t ) ) such that
| v 1 ( t ) w | m ( t ) | z ( t ) z ¯ ( t ) | , t [ a , b ] .
Define V : [ a , b ] P ( R ) by
V ( t ) = { w R : | v 1 ( t ) w | m ( t ) | z ( t ) z ¯ ( t ) | } .
According to Proposition III.4 [34], the multivalued operator V ( t ) F ( t , z ¯ ( t ) ) is measurable, and thus ∃ a function v 2 ( t ) which is a measurable selection for V . So v 2 ( t ) F ( t , z ¯ ( t ) ) and for each t [ a , b ] , we have | v 1 ( t ) v 2 ( t ) | m ( t ) | z ( t ) z ¯ ( t ) | .
For each t [ a , b ] , let us define
h 2 ( t ) = ( ψ ( t ) ψ ( a ) ) θ k k 1 Λ Γ k ( θ k ) λ k I α ; ψ f 2 ( ξ ) + μ k I α + v ; ψ f 2 ( σ ) k I α ; ψ f 2 ( b ) + k I α ; ψ f 2 ( t ) .
Then, we have
| h 1 ( t ) h 2 ( t ) | ( ψ ( t ) ψ ( a ) ) θ k k 1 | Λ | Γ k ( θ k ) [ | λ | k I α ; ψ ( | f 1 ( s ) f 2 ( s ) | ) ( ξ ) + | μ | k I α + v ; ψ ( | f 1 ( s ) f 2 ( s ) | ) ( σ ) + k I α ; ψ ( | f 1 ( s ) f 2 ( s ) | ) ( b ) ] + k I α ; ψ ( | f 1 ( s ) f 2 ( s ) | ) ( t ) { ( ψ ( b ) ψ ( a ) ) α k Γ k ( α + k ) + ( ψ ( b ) ψ ( a ) ) θ k k 1 | Λ | Γ k ( θ k ) [ | λ | ( ψ ( ξ ) ψ ( a ) ) α k Γ k ( α + k ) + | μ | ( ψ ( b ) ψ ( a ) ) α + v k Γ k ( α + v + k ) + ( ψ ( b ) ψ ( a ) ) α k Γ k ( α + k ) ] } m z z ¯ = G m z z ¯ .
Hence
h 1 h 2 G m z z ¯ .
Analogously, by interchanging the roles of z and z ¯ , we obtain
H d ( F ( z ) , F ( z ¯ ) ) G m z z ¯ .
So, F is a contraction, and thus F has a fixed point z by application of the Covitz and Nadler theorem [33]. Consequently, there exists at least one solution on [ a , b ] to the problem (2). The proof is finished. □

6. Examples

In this section, some examples, illustrating the obtained theoretical results in the previous section, are presented.
Consider the following ( k , ψ ) -Hilfer type nonlocal integral fractional boundary value problem
13 10 , H D 8 5 , 3 4 ; t + 1 z ( t ) = f ( t , z ( t ) ) , 1 3 < t < 7 3 , z 1 3 = 0 , z 7 3 = 1 π z 2 3 + 6 17 13 10 I 4 3 , w + 1 z 5 3 .
Here, k = 13 / 10 , α = 8 / 5 , β = 3 / 4 , ψ ( t ) = t + 1 , a = 1 / 3 , b = 7 / 3 , λ = 1 / π , ξ = 2 / 3 , μ = 6 / 17 , v = 4 / 3 , σ = 5 / 3 , and we can compute that θ 13 10 = 47 / 20 , Λ 0.9298349818 , G 1.182496327 , G 1 0.5758310296 .
Example 1.
Let f be a nonlinear unbounded Lipschitz function defined by
f ( t , z ) = z 2 + 10 | z | 2 ( 6 + | z | ) cos 2 π w + 1 3 w 4 + 2 w 3 + 1 .
Clearly, f satisfies the Lipschitz condition with constant L = 5 / 6 as
| f ( t , z 1 ) f ( t , z 2 ) | 5 6 | z 1 z 2 | ,
for all z 1 , z 2 R and t [ 1 / 3 , 7 / 3 ] . Moreover, L G 0.9854136058 < 1 . Hence, by Theorem 1, the problem (14) with f given in (15), has a unique solution on [ 1 / 3 , 7 / 3 ] .
Example 2.
Let a nonlinear bounded Lipschitz function f be given by
f ( t , z ) = 10 | z | 6 + | z | e | 3 w 1 | + 1 4 w 3 + 1 3 w 2 + 4 5 ,
Note that
| f ( t , z ) | 10 e | 3 w 1 | + 1 4 w 3 + 1 3 w 2 + 4 5 : = φ ( t ) ,
for all t [ 1 / 3 , 7 / 3 ] . Observe that f satisfies the Lipschitz condition | f ( t , z 1 ) f ( t , z 2 ) | ( 10 / 6 ) | z 1 z 2 | , with the Lipschitz constant L = 5 / 3 . Using the given data, we obtain L G 1 0.9597183827 < 1 . As a consequence, the conclusion of Theorem 2 applies and hence, the problem (14) with f given by (16) has at least one solution on [ 1 / 3 , 7 / 3 ] . It is imperative to notice that the uniqueness of the solution for this problem cannot be guaranteed since L G 1.970827212 > 1 .
Example 3.
Consider the function
f ( t , z ) = σ ( t ) A g 1 ( z ) + B ,
where σ : [ 1 / 3 , 7 / 3 ] R + , g 1 : R R with | g 1 ( z ) | | z | , 0 A < 1 / ( σ G ) and B > 0 . Then, we have
| f ( t , z ) | σ A | z | + B .
Setting χ ( u ) = A | u | + B , we can find a constant K satisfying the condition ( H 4 ) of Theorem 3 as
K > B σ G 1 A σ G .
By applying Theorem 3, we deduce that the problem (14) with f given in (17) has at least one solution on [ 1 / 3 , 7 / 3 ] .
Example 4.
Let f be defined by
f ( t , z ) = σ ( t ) A g 2 ( z ) + B ,
where σ : [ 1 / 3 , 7 / 3 ] R + , g 2 : R R with | g 2 ( z ) | z 2 , and A , B are positive constants with A B < 1 / ( 4 σ 2 G 2 ) . Then, we obtain
| f ( t , z ) | σ A z 2 + B .
Choosing a function χ ( u ) = A u 2 + B , we find a constant
K 1 1 4 A B σ 2 G 2 2 A σ G , 1 + 1 4 A B σ 2 G 2 2 A σ G ,
satisfying the condition ( H 4 ) of Theorem 3. Thus, all the assumptions of Theorem 3 are satisfied. Hence, the problem (14) with f given by (18) has at least one solution on [ 1 / 3 , 7 / 3 ] .
Example 5.
Assume that the first equation of (14) is replaced by
13 10 , H D 8 5 , 3 4 ; t + 1 z ( t ) F ( t , z ( t ) ) , 1 3 < t < 7 3 ,
where
F ( t , z ) = 0 , 1 ( 3 t + 2 ) z 2 + 2 | z | 1 + | z | + sin 2 t .
Observe that F ( t , z ) is a measurable set. Additionally,
H d F ( t , z ) , F ( t , z ¯ ) 2 ( 3 t + 2 ) | z z ¯ | .
Now, we set m ( t ) = 2 / ( 3 t + 2 ) such that d ( 0 , F ( t , 0 ) ) 1 / ( 3 t + 2 ) < 2 / ( 3 t + 2 ) = m ( t ) for almost all t [ 1 / 3 , 7 / 3 ] . From δ = G m 0.7883308847 < 1 , we deduce that the ( k , ψ ) -Hilfer type fractional inclusion (19) with nonlocal integral boundary conditions given in (14), has at least one solution on [ 1 / 3 , 7 / 3 ] .

7. Conclusions

We presented the existence criteria for solutions to the ( k , ψ ) -Hilfer type fractional differential equations and inclusions of order in [ 1 , 2 ] complemented with nonlocal integral boundary conditions. We first transformed the nonlinear ( k , ψ ) -Hilfer type fractional boundary value problem into a fixed point problem. For the single-valued case, we established existence and uniqueness results by applying the Banach contraction mapping principle, Krasnosel’skiĭ fixed point theorem and the Leray–Schauder nonlinear alternative. Our first existence result dealing with the convex-valued multi-valued map involved in the inclusion was established by applying the Leray–Schauder nonlinear alternative for multivalued maps, while the existence result for the non-convex valued multivalued map in the inclusion was obtained by applying the Covitz–Nadler fixed point theorem for contractive multivalued maps. It is worthwhile to mention that the work established for ( k , ψ ) -Hilfer fractional differential equations supplemented with nonlocal ( k , ψ ) -Riemann–Liouville fractional integral boundary conditions is more general and significant as the ( k , ψ ) -Riemann–Liouville and ( k , ψ ) -Caputo fractional derivatives are special cases of the ( k , ψ ) -Hilfer fractional derivative. Moreover, the ( k , ψ ) -Riemann–Liouville fractional integral operator used in the boundary conditions is of a more general nature.

Author Contributions

Conceptualization, S.K.N.; methodology, S.K.N., B.A. and J.T.; validation, S.K.N., B.A. and J.T.; formal analysis, S.K.N., B.A. and J.T.; writing—original draft preparation, S.K.N., B.A. and J.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Diethelm, K. The Analysis of Fractional Differential Equations; Lecture Notes in Mathematics; Springer: New York, NY, USA, 2010. [Google Scholar]
  2. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of the Fractional Differential Equations; North-Holland Mathematics Studies; Elsevier: Amsterdam, The Netherlands, 2006; Volume 204. [Google Scholar]
  3. Lakshmikantham, V.; Leela, S.; Devi, J.V. Theory of Fractional Dynamic Systems; Cambridge Scientific Publishers: Cambridge, UK, 2009. [Google Scholar]
  4. Miller, K.S.; Ross, B. An Introduction to the Fractional Calculus and Differential Equations; John Wiley: New York, NY, USA, 1993. [Google Scholar]
  5. Podlubny, I. Fractional Differential Equations; Academic Press: New York, NY, USA, 1999. [Google Scholar]
  6. Samko, S.G.; Kilbas, A.A.; Marichev, O.I. Fractional Integrals and Derivatives; Gordon and Breach Science: Yverdon, Switzerland, 1993. [Google Scholar]
  7. Ahmad, B.; Alsaedi, A.; Ntouyas, S.K.; Tariboon, J. Hadamard-Type Fractional Differential Equations, Inclusions and Inequalities; Springer: Cham, Switzerland, 2017. [Google Scholar]
  8. Ahmad, B.; Ntouyas, S.K. Nonlocal Nonlinear Fractional-Order Boundary Value Problems; World Scientific: Singapore, 2021. [Google Scholar]
  9. Zhou, Y. Basic Theory of Fractional Differential Equations; World Scientific: Singapore, 2014. [Google Scholar]
  10. Katugampola, U.N. New approach to a generalized fractional integral. Appl. Math. Comput. 2011, 218, 860–865. [Google Scholar] [CrossRef]
  11. Katugampola, U.N. A new approach to generalized fractional derivatives. Bull. Math. Anal. Appl. 2014, 6, 1–154. [Google Scholar]
  12. Hilfer, R. (Ed.) Applications of Fractional Calculus in Physics; World Scientific: Singapore, 2000. [Google Scholar]
  13. Almeida, R. A Caputo fractional derivative of a function with respect to another function. Commun. Nonlinear Sci. Numer. Simul. 2017, 44, 460–481. [Google Scholar] [CrossRef]
  14. Kucche, K.D.; Mali, A.D. On the nonlinear (k,ψ)-Hilfer fractional differential equations. Chaos Solitons Fractals 2021, 152, 111335. [Google Scholar] [CrossRef]
  15. Tariboon, J.; Samadi, A.; Ntouyas, S.K. Multi-point boundary value problems for (k,ψ)-Hilfer fractional differential equations and inclusions. Axioms 2022, 11, 110. [Google Scholar] [CrossRef]
  16. Bulavatsky, V.M. Mathematical modeling of fractional differential filtration dynamics based on models with Hilfer–Prabhakar derivative. Cybern. Syst. Anal. 2017, 53, 204–216. [Google Scholar] [CrossRef]
  17. Bulavatsky, V.M. Mathematical models and problems of fractional-differential dynamics of some relaxation filtration processes. Cybern. Syst. Anal. 2018, 54, 727–736. [Google Scholar] [CrossRef]
  18. Ali, I.; Malik, N. Hilfer fractional advection-diffusion equations with power-law initial condition; a numerical study using variational iteration method. Comput. Math. Appl. 2014, 68, 1161–1179. [Google Scholar] [CrossRef]
  19. Hilfer, R. Experimental evidence for fractional time evolution in glass forming materials. Chem. Phys. 2002, 284, 399–408. [Google Scholar] [CrossRef]
  20. Andreu-Vaillo, F.; Mazon, J.M.; Rossi, J.D.; Toledo-Melero, J.J. Nonlocal Diffusion Problems. Mathematical Surveys and Monographs; American Mathematical Society: Providence, RI, USA, 2010; Volume 165. [Google Scholar]
  21. Alotta, G.; Paola, M.D.; Pinnola, F.P.; Zingales, M. A fractional nonlocal approach to nonlinear blood flow in small-lumen arterial vessels. Meccanica 2020, 55, 891–906. [Google Scholar] [CrossRef]
  22. Ciegis, R.; Bugajev, A. Numerical approximation of one model of bacterial self-organization. Nonlinear Anal. Model. Control 2012, 17, 253–270. [Google Scholar] [CrossRef]
  23. Mubeen, S.; Habibullah, G.M. k–fractional integrals and applications. Int. J. Contemp. Math. Sci. 2012, 7, 89–94. [Google Scholar]
  24. Diaz, R.; Pariguan, E. On hypergeometric functions and Pochhammer k-symbol. Divulg. Mat. 2007, 2, 179–192. [Google Scholar]
  25. Dorrego, G.A. An alternative definition for the k-Riemann-Liouville fractional derivative. Appl. Math. Sci. 2015, 9, 481–491. [Google Scholar] [CrossRef]
  26. Kwun, Y.C.; Farid, G.; Nazeer, W.; Ullah, S.; Kang, S.M. Generalized Riemann-Liouville k-fractional integrals associated with Ostrowski type inequalities and error bounds of Hadamard inequalities. IEEE Access 2018, 6, 64946–64953. [Google Scholar] [CrossRef]
  27. Vanterler, J.; de Oliveira, E.C. On the ψ-Hilfer fractional derivative. Commun. Nonlinear Sci. Numer. Simul. 2018, 60, 72–91. [Google Scholar] [CrossRef]
  28. Deimling, K. Nonlinear Functional Analysis; Springer: New York, NY, USA, 1985. [Google Scholar]
  29. Krasnosel’skiĭ, M.A. Two remarks on the method of successive approximations. Uspekhi Mat. Nauk 1955, 10, 123–127. [Google Scholar]
  30. Granas, A.; Dugundji, J. Fixed Point Theory; Springer: New York, NY, USA, 2005. [Google Scholar]
  31. Deimling, K. Multivalued Differential Equations; Walter De Gruyter: Berlin, Germany, 1992. [Google Scholar]
  32. Lasota, A.; Opial, Z. An application of the Kakutani-Ky Fan theorem in the theory of ordinary differential equations. Bull. Acad. Polon. Sci. Ser. Sci. Math. Astronom. Phys. 1965, 13, 781–786. [Google Scholar]
  33. Covitz, H.; Nadler, S.B., Jr. Multivalued contraction mappings in generalized metric spaces. Israel J. Math. 1970, 8, 5–11. [Google Scholar] [CrossRef]
  34. Castaing, C.; Valadier, M. Convex Analysis and Measurable Multifunctions; Lecture Notes in Mathematics 580; Springer: Berlin/Heidelberg, Germany, 1977. [Google Scholar]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ntouyas, S.K.; Ahmad, B.; Tariboon, J. Nonlocal Boundary Value Problems for (k,ψ)-Hilfer Fractional Differential Equations and Inclusions. Foundations 2022, 2, 681-696. https://doi.org/10.3390/foundations2030046

AMA Style

Ntouyas SK, Ahmad B, Tariboon J. Nonlocal Boundary Value Problems for (k,ψ)-Hilfer Fractional Differential Equations and Inclusions. Foundations. 2022; 2(3):681-696. https://doi.org/10.3390/foundations2030046

Chicago/Turabian Style

Ntouyas, Sotiris K., Bashir Ahmad, and Jessada Tariboon. 2022. "Nonlocal Boundary Value Problems for (k,ψ)-Hilfer Fractional Differential Equations and Inclusions" Foundations 2, no. 3: 681-696. https://doi.org/10.3390/foundations2030046

APA Style

Ntouyas, S. K., Ahmad, B., & Tariboon, J. (2022). Nonlocal Boundary Value Problems for (k,ψ)-Hilfer Fractional Differential Equations and Inclusions. Foundations, 2(3), 681-696. https://doi.org/10.3390/foundations2030046

Article Metrics

Back to TopTop