Next Article in Journal
Harmony Search-Based Approach for Multi-Objective Software Architecture Reconstruction
Next Article in Special Issue
Quasilinear Dirichlet Problems with Degenerated p-Laplacian and Convection Term
Previous Article in Journal
Role of Economic Policy Uncertainty in the Connectedness of Cross-Country Stock Market Volatilities
Previous Article in Special Issue
Fractional-Order Integro-Differential Multivalued Problems with Fixed and Nonlocal Anti-Periodic Boundary Conditions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Boundary Value Problems for Hilfer Fractional Differential Inclusions with Nonlocal Integral Boundary Conditions

by
Athasit Wongcharoen
1,†,
Sotiris K. Ntouyas
2,3,† and
Jessada Tariboon
4,*,†
1
Department of Mechanical Engineering Technology, College of Industrial Engineering Technology, King Mongkut’s University of Technology North Bangkok, Bangkok 10800, Thailand
2
Department of Mathematics, University of Ioannina, 451 10 Ioannina, Greece
3
Nonlinear Analysis and Applied Mathematics (NAAM)-Research Group, Department of Mathematices, Faculty of Science, King Abdulaziz University, P.O. Box 80203, Jeddah 21589, Saudi Arabia
4
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.
These authors contributed equally to this work.
Submission received: 2 October 2020 / Revised: 26 October 2020 / Accepted: 28 October 2020 / Published: 31 October 2020
(This article belongs to the Special Issue Nonlinear Equations: Theory, Methods, and Applications)

Abstract

:
In this paper, we study boundary value problems for differential inclusions, involving Hilfer fractional derivatives and nonlocal integral boundary conditions. New existence results are obtained by using standard fixed point theorems for multivalued analysis. Examples illustrating our results are also presented.

1. Introduction

The theory of fractional differential equations received in recent years considerable interest both in pure mathematics and applications, see [1,2,3,4,5,6,7,8,9]. In the literature, there exist several different definitions of fractional integrals and derivatives, for example, the most popular of them are fractional derivatives in the sense of Riemann–Liouville and Caputo. Other known definitions are the Hadamard fractional derivative, the Erdeyl–Kober fractional derivative, and so on. We refer the interested in fractional calculus reader to the classical reference texts such as [1,2,3,4,5,6,7]. A generalization of derivatives of both Riemann–Liouville and Caputo was given by R. Hilfer in [10] when he studied fractional time evolution in physical phenomena. He named it as generalized fractional derivative of order α ( 0 , 1 ) and a type β [ 0 , 1 ] which can be reduced to the Riemann–Liouville and Caputo fractional derivatives when β = 0 and β = 1 , respectively. Such a derivative interpolates between the Riemann–Liouville and Caputo derivative. Some properties and applications of the Hilfer derivative are given in [11,12] and references cited therein.
Initial value problems involving Hilfer fractional derivatives were studied by several authors, see for example [13,14,15] and references therein. Recently, in [16], the authors initiated the study of boundary value problems for Hilfer fractional differential equations with nonlocal integral boundary conditions of the form
H D α , β x ( t ) = f ( t , x ( t ) ) , t J : = [ a , b ] , 1 < α < 2 , 0 β 1 ,
x ( a ) = 0 , x ( b ) = i = 1 m δ i I φ i x ( ξ i ) , φ i > 0 , δ i R , ξ i [ a , b ] ,
where H D α , β is the Hilfer fractional derivative of order α , 1 < α < 2 and parameter β , 0 β 1 , I φ i is the Riemann–Liouville fractional integral of order φ i > 0 , ξ i J , a 0 , δ i R , and f : J × R R is a continuous function. Several existence and uniqueness results were proved by using a variety of fixed point theorems, such as Banach’s fixed point theorem, Hölder’s inequality, and Boyd and Wong fixed point theorem for nonlinear contractions to obtain the uniqueness results, while nonlinear alternative of Leray–Schauder type and Krasnoselskii’s fixed point theorem were applied to obtain the existence results. For some other recent results, see [17,18].
The objective of the present work is to study the multivalued case of the problem (1) and (2); that is, the boundary value problem of Hilfer-type fractional differential inclusions with nonlocal integral boundary conditions
H D α , β x ( t ) F ( t , x ( t ) ) , t J , 1 < α < 2 , 0 β 1 ,
x ( a ) = 0 , x ( b ) = i = 1 m δ i I φ i x ( ξ i ) , φ i > 0 , δ i R , ξ i [ a , b ] ,
where F : J × R P ( R ) is a multivalued map, P ( R ) is the family of all nonempty subsets of R , and the other parameters are as in problem (1) and (2).
We obtain existence results for the problem (3) and (4) by applying the nonlinear alternative of Leray–Schauder type, and a fixed point theorem due to Nadler. Examples are also constructed for the illustration of main results.
The rest of the paper is arranged as follows. In Section 2, we recall some preliminary concepts of fractional calculus and multivalued analysis. Section 3 is devoted to the main results, while examples illustrating the obtained results are presented in Section 4.

2. Preliminaries

In this section, we introduce some notation and definitions of fractional calculus and present preliminary results needed in our proofs later [2,5].
Definition 1.
The Riemann–Liouville fractional integral of order α > 0 of a continuous function is defined by
I α u ( t ) = 1 Γ ( α ) a t ( t s ) α 1 u ( s ) d s , n 1 < α < n ,
where n = [ α ] + 1 , [ α ] denotes the integer part of real number α, provided the right-hand side is point-wise defined on ( a , ) , where Γ ( · ) is the Euler gamma function, which is defined by Γ ( r ) = 0 t r 1 e t d t .
Definition 2.
The Riemann–Liouville fractional derivative of order α > 0 of a continuous function is defined by
R L D α u ( t ) : = D n I n α u ( t ) = 1 Γ ( n α ) d d t n a t ( t s ) n α 1 u ( s ) d s , n 1 < α < n ,
where n = [ α ] + 1 , provided the right-hand side is point-wise defined on ( a , ) .
Definition 3.
The Caputo fractional derivative of order α > 0 of a continuous function is defined by
C D α u ( t ) : = I n α D n u ( t ) = 1 Γ ( n α ) a t ( t s ) n α 1 d d s n u ( s ) d s , n 1 < α < n ,
provided the right-hand side is point-wise defined on ( a , ) .
In [10] (see also [12]), another new definition of the fractional derivative was suggested. The generalized Riemann–Liouville fractional derivative defined as
Definition 4.
The generalized Riemann–Liouville fractional derivative or Hilfer fractional derivative of order α and parameter β of a function is defined by
H D α , β u ( t ) = I β ( n α ) D n I ( 1 β ) ( n α ) u ( t ) ,
where n 1 < α < n , 0 β 1 , t > a , D = d d t .
Remark 1.
When β = 0 , the Hilfer fractional derivative corresponds to the Riemann–Liouville fractional derivative
H D α , 0 u ( t ) = D n I n α u ( t ) ,
while when β = 1 the Hilfer fractional derivative corresponds to the Caputo fractional derivative
H D α , 1 u ( t ) = I n α D n u ( t ) .
In the following lemma we present the compositional property of Riemann–Liouville fractional integral operator with the Hilfer fractional derivative operator.
Lemma 1
([12]). Let f L ( a , b ) , n 1 < α n , n N , 0 β 1 , I ( n α ) ( 1 β ) f A C k [ a , b ] . Then
I α H D α , β f ( t ) = f ( t ) k = 0 n 1 ( t a ) k ( n α ) ( 1 β ) Γ ( k ( n α ) ( 1 β ) + 1 ) lim t a + I ( 1 β ) ( n α ) f ( t ) .
The following lemma deals with a linear variant of the boundary value problem (3) and (4) and is proved in [16].
Lemma 2.
Let φ i > 0 , ξ i [ a , b ] , a 0 , δ i R , 1 < α < 2 , γ = α + 2 β α β , Λ 0 , and h C ( [ a , b ] , R ) . Then the function x is a solution of the boundary value
H D α , β x ( t ) = h ( t ) , t [ a , b ] , 1 < α < 2 , 0 β 1 ,
x ( a ) = 0 , x ( b ) = i = 1 m δ i I φ i x ( ξ i ) , φ i > 0 , δ i R , ξ i [ a , b ] ,
if and only if
x ( t ) = ( t a ) γ 1 Λ Γ ( γ ) I α h ( b ) ) i = 1 m δ i I α + φ i h ( ξ i ) + I α h ( t ) ,
where
Λ = i = 1 m δ i ( ξ i a ) γ + φ i 1 Γ ( γ + φ i ) ( b a ) γ 1 Γ ( γ ) .

3. Existence Results

By C ( J , R ) , we denote the Banach space of all continuous functions from J into R with the norm
x : = sup { | x ( t ) | : t J } .
Additionally by L 1 ( J , R ) , we denote the space of functions x : J R such that x L 1 = 1 b | x ( t ) | d t .
For a normed space ( X , · ) , we define: P ( X ) = { Y X : Y } ; P c p ( X ) = { Y P ( X ) : Y   i s   c o m p a c t } ; P c , c p ( X ) = { Y P ( X ) : Y   i s   c o n v e x   a n d   c o m p a c t } ; P b , c l ( X ) = { Y P ( X ) : Y   i s   b o u n d e d   a n d   c l o s e d } and P c l ( X ) = { Y P ( X ) : Y   i s   c l o s e d } .
For the basic concepts of multivalued analysis we refer to [19,20]. See also [21].
For each y C ( J , R ) , define the set of selections of F by
S F , x : = { v L 1 ( J , R ) : v ( t ) F ( t , y ( t ) ) o n J } .
A multivalued map G : X P ( X ) has a fixed point if there is x X such that x G ( x ) . The fixed point set of the multivalued operator G will be denoted by Fix G .
Definition 5.
A function x C ( J , R ) is said to be a solution of the problem (3) and (4) if there exists a function v L 1 ( J , R ) with v ( t ) F ( t , x ) for a.e. t J such that x satisfies the differential equation H D α , β x ( t ) = v ( t ) on J and the boundary conditions x ( a ) = 0 , x ( b ) = i = 1 m δ i I φ i x ( ξ i ) .
Definition 6.
A multivalued map F : J × R P ( R ) is said to be Carathéodory if
(i)
t F ( t , x ) is measurable for each x R ;
(ii)
x F ( t , x ) is upper semicontinuous for almost all t J . Further a Carathéodory function F is called L 1 —Carathéodory if
(iii)
for each ρ > 0 , there exists φ ρ L 1 ( J , R + ) such that
F ( t , x ) = sup { | v | : v F ( t , x ) } φ ρ ( t )
for all x R with x ρ and for a.e. t J .
Our first existence result, dealing with the case when F has convex values, is based on nonlinear alternative of Leray–Schauder type [22] with the assumption that F is Carathéodory. For quick reference, we state the well-known nonlinear alternative of Leray–Schauder for multivalued maps.
Lemma 3.
(Nonlinear alternative for Kakutani maps [22]). Let E be a Banach space, C a closed convex subset of E , U an open subset of C and 0 U . Suppose that F : U ¯ P c p , c ( C ) is an upper semicontinuous compact map. Then either
(i)
F has a fixed point in U ¯ , or
(ii)
there is a u U and λ ( 0 , 1 ) with u λ F ( u ) .
Theorem 1.
Assume that:
(H1)
F : J × R P c p , c ( R ) is L 1 -Carathéodory;
(H2)
there exists a continuous nondecreasing function ψ : [ 0 , ) ( 0 , ) and a function p C ( J , R + ) such that
F ( t , x ) P : = sup { | x | : x F ( t , x ) } p ( t ) ψ ( x ) f o r   e a c h ( t , x ) J × R ;
(H3)
there exists a constant M > 0 such that
M p ψ ( M ) ( b a ) α Γ ( α + 1 ) + ( b a ) γ 1 | Λ | Γ ( γ ) ( b a ) α Γ ( α + 1 ) + i = 1 m | δ i | ( ξ i a ) α + φ i Γ ( α + φ i + 1 ) > 1 .
Then the boundary value problem (3) and (4) has at least one solution on [ a , b ] .
Proof. 
To transform the problem (3) and (4) into a fixed point problem, we define an operator F : C ( [ a , b ] , R ) P ( C ( [ a , b ] , R ) ) by
F ( x ) = h C ( [ a , b ] , R ) : h ( t ) = ( t a ) γ 1 Λ Γ ( γ ) I α v ( b ) i = 1 m δ i I α + φ i h ( ξ i ) + I α v ( t )
for t [ a , b ] and v S F , x . It is obvious that the fixed points of F are solutions of the boundary value problem (3) and (4).
We will show that the operator F satisfies the assumptions of Leray–Schauder nonlinear alternative [22] in a series of steps.
Step 1. F ( x ) is convex for each x C ( [ a , b ] , R ) .
This step is obvious since S F , x is convex (F has convex values), and therefore we omit the proof.
Step 2. F maps bounded sets (balls) into bounded sets in C ( [ a , b ] , R ) .
Let B r = { x C ( [ a , b ] , R ) : x r } be a bounded set in C ( [ a , b ] , R ) . Then, for each h B ( x ) , x B r , there exists v S F , x such that
h ( t ) = I α v ( s ) ( t ) + ( t a ) γ 1 Λ Γ ( γ ) I α v ( s ) ( b ) i = 1 m δ i I α + φ i v ( s ) ( ξ i ) .
Then, for t J , we have
| h ( t ) | I α | v ( s ) | ( b ) + ( b a ) γ 1 | Λ | Γ ( γ ) I α | v ( s ) | ( b ) + i = 1 m δ i I α + φ i | v ( s ) | ( ξ i ) p ψ ( x ) ( b a ) α Γ ( α + 1 ) + ( b a ) γ 1 | Λ | Γ ( γ ) ( b a ) α Γ ( α + 1 ) + i = 1 m | δ i | ( ξ i a ) α + φ i Γ ( α + φ i + 1 ) .
Thus,
h p ψ ( r ) ( b a ) α Γ ( α + 1 ) + ( b a ) γ 1 | Λ | Γ ( γ ) ( b a ) α Γ ( α + 1 ) + i = 1 m | δ i | ( ξ i a ) α + φ i Γ ( α + φ i + 1 ) .
Step 3. F maps bounded sets into equicontinuous sets of C ( [ a , b ] , R ) .
Let τ 1 , τ 2 [ a , b ] with τ 1 < τ 2 and x B r . Then, for each h B ( x ) , we obtain
| h ( t 2 ) h ( t 1 ) | ( τ 2 a ) γ 1 ( τ 1 a ) γ 1 | Λ | Γ ( γ ) I α | v ( s ) | ( b ) + ( τ 2 a ) γ 1 ( τ 1 a ) γ 1 | Λ | Γ ( γ ) i = 1 m | δ i | I α + φ i | v ( s ) | ( ξ i ) + 1 Γ ( α ) | a τ 1 [ ( τ 2 s ) α 1 ( τ 1 s ) α 1 ] v ( s ) d s + τ 1 τ 2 ( τ 2 s ) α 1 v ( s ) d s | ( τ 2 a ) γ 1 ( τ 1 a ) γ 1 | Λ | Γ ( γ ) ψ ( r ) I α p ( s ) ( b ) + ( τ 2 a ) γ 1 ( τ 1 a ) γ 1 | Λ | Γ ( γ ) ψ ( r ) i = 1 m | δ i | I α + φ i p ( s ) ( ξ i ) + ψ ( r ) Γ ( α ) | a τ 1 [ ( τ 2 s ) α 1 ( τ 1 s ) α 1 ] p ( s ) d s + τ 1 τ 2 ( τ 2 s ) α 1 p ( s ) d s | ( τ 2 a ) γ 1 ( τ 1 a ) γ 1 | Λ | Γ ( γ ) p ψ ( r ) ( b a ) α Γ ( α + 1 ) + i = 1 m | δ i | ( ξ i a ) α + φ i Γ ( α + φ i + 1 ) + p ψ ( r ) Γ ( α + 1 ) [ 2 ( τ 2 τ 1 ) α + | ( τ 2 a ) α ( τ 1 a ) α | ] .
Obviously the right hand side of the above inequality tends to zero independently of x B r as τ 2 τ 1 0 . Therefore, it follows by the Arzelá–Ascoli theorem that F : C ( [ a , b ] , R ) P ( C ( [ a , b ] , R ) ) is completely continuous.
Next we show that F is an upper semi-continuous multivalued mapping. It is known by [19] (Proposition 1.2) that F will be upper semicontinuous if we establish that it has a closed graph. Thus, in the next step we show that F has a closed graph.
Step 4. F has a closed graph.
Let x n x * , h n F ( x n ) and h n h * . Then we need to show that h * F ( x * ) . Associated with h n F ( x n ) , there exists v n S F , x n such that for each t [ a , b ] ,
h n ( t ) = I α v n ( s ) ( t ) + ( t a ) γ 1 Λ Γ ( γ ) I α v n ( s ) ( b ) i = 1 m δ i I α + φ i v n ( s ) ( ξ i ) .
Thus it suffices to show that there exists v * S F , x * such that for each t [ a , b ] ,
h * ( t ) = I α v ( s ) ( t ) + ( t a ) γ 1 Λ Γ ( γ ) I α v * ( s ) ( b ) i = 1 m δ i I α + φ i v * ( s ) ( ξ i ) .
Let us consider the linear operator Θ : L 1 ( J , R ) C ( J , R ) given by
v Θ ( v ) ( t ) = I α v ( s ) ( t ) + ( t a ) γ 1 Λ Γ ( γ ) I α v ( s ) ( b ) i = 1 m δ i I α + φ i v ( s ) ( ξ i ) .
Observe that h n ( t ) h * ( t ) 0 , as n . Therefore, it follows by a Lazota–Opial result [23], that Θ S F is a closed graph operator. Further, we have h n ( t ) Θ ( S F , x n ) . Since x n x * , we have
h * ( t ) = I α v * ( s ) ( t ) + ( t a ) γ 1 Λ Γ ( γ ) I α v * ( s ) ( b ) i = 1 m δ i I α + φ i v * ( s ) ( ξ i ) ,
for some v * S F , x * .
Step 5.We show there exists an open set U C ( [ a , b ] , R ) with x ν F ( x ) for any ν ( 0 , 1 ) and all y U .
Let ν ( 0 , 1 ) and x ν F ( x ) . Then there exists v L 1 ( [ a , b ] , R ) with v S F , x such that, for t [ a , b ] , we have
x ( t ) = ν I α v ( s ) ( t ) + ν ( t a ) γ 1 Λ Γ ( γ ) I α v ( s ) ( b ) i = 1 m δ i I α + φ i v ( s ) ( ξ i ) ,
As in the second step, it can be shown that
| x ( t ) | p ψ ( x ) ( b a ) α Γ ( α + 1 ) + ( b a ) γ 1 | Λ | Γ ( γ ) ( b a ) α Γ ( α + 1 ) + i = 1 m | δ i | ( ξ i a ) α + φ i Γ ( α + φ i + 1 ) .
which implies that
x p ψ ( x ) ( b a ) α Γ ( α + 1 ) + ( b a ) γ 1 | Λ | Γ ( γ ) ( b a ) α Γ ( α + 1 ) + i = 1 m | δ i | ( ξ i a ) α + φ i Γ ( α + φ i + 1 ) .
Consequently
x p ψ ( x ) ( b a ) α Γ ( α + 1 ) + ( b a ) γ 1 | Λ | Γ ( γ ) ( b a ) α Γ ( α + 1 ) + i = 1 m | δ i | ( ξ i a ) α + φ i Γ ( α + φ i + 1 ) 1 .
In view of ( H 3 ) , there exists M such that x M . Let us set
U = { x C ( [ a , b ] , R ) : x < M } .
Note that the operator F : U ¯ P ( C ( [ a , b ] , R ) ) is a compact multivalued map, u.s.c. with convex closed values. From the choice of U, there is no x U such that x ν F ( x ) for some ν ( 0 , 1 ) . In consequence, we deduce by the nonlinear alternative of Leray–Schauder type [22] that F has a fixed point y U ¯ which is a solution of the problem (3) and (4). This completes the proof. □
In our second result, we show the existence of solutions for the problem (3) and (4) when F is not necessary nonconvex valued by applying a fixed point theorem for multivalued maps due to Nadler [24]. Let us first describe the necessary material needed to establish this result.
Let ( X , d ) be a metric space induced from the normed space ( X ; · ) . Consider H d : P ( X ) × P ( X ) R { } given by
H d ( A , B ) = max { sup a A d ( a , B ) , sup b B d ( A , b ) } ,
where d ( A , b ) = inf a A d ( a ; b ) and d ( a , B ) = inf b B d ( a ; b ) . Then ( P c l , b ( X ) , H d ) is a metric space (see [25]).
Definition 7.
A multivalued operator N : X P c l ( X ) is called
(a) 
θ-Lipschitz if and only if there exists θ > 0 such that
H d ( N ( x ) , N ( y ) ) θ d ( x , y ) f o r   e a c h x , y X ;
(b) 
a contraction if and only if it is θ-Lipschitz with θ < 1 .
Lemma 4
([24]). Let ( X , d ) be a complete metric space. If N : X P b , c l ( X ) is a multivalued contraction then N has a fixed point.
Theorem 2.
Assume that the following conditions hold:
(A1)
F : [ a , b ] × R P c p ( R ) is such that F ( · , x ) : [ a , b ] P c p ( R ) is measurable for each x R .
(A2)
H d ( F ( t , x ) , F ( t , x ¯ ) ) m ( t ) | x x ¯ | for almost all t J and x , x ¯ R with m C ( [ a , b ] , R + ) and d ( 0 , F ( t , 0 ) ) m ( t ) for almost all t [ a , b ] .
Then the boundary value problem (3) and (4) has at least one solution on [ a , b ] if
δ : = ( b a ) α Γ ( α + 1 ) + ( b a ) γ 1 | Λ | Γ ( γ ) ( b a ) α Γ ( α + 1 ) + i = 1 m | δ i | ( ξ i a ) α + φ i Γ ( α + φ i + 1 ) m < 1 .
Proof. 
Consider the operator F defined by (9). Observe that the set S F , x is nonempty for each y C ( [ a , b ] , R ) by the assumption ( A 1 ) , so F has a measurable selection (see Theorem III.6 [26]). Now we show that the operator F satisfies the assumptions of Nadler theorem [24]. We show that F ( x ) P c l ( C ( [ a , b ] , R ) ) for each x C ( [ a , b ] , R ) . Let { u n } n 0 F ( x ) be such that u n u ( n ) in C ( [ a , b ] , R ) . Then u C ( [ a , b ] , R ) and there exists v n S F , x n such that, for each t [ a , b ] ,
u n ( t ) = I α v n ( s ) ( t ) + ( t a ) γ 1 Λ Γ ( γ ) I α v n ( s ) ( b ) i = 1 m δ i I α + φ i v n ( s ) ( ξ i )
As F has compact values, we pass onto a subsequence (if necessary) to obtain that v n converges to v in L 1 ( [ a , b ] , R ) . Thus, v S F , x and for each t J , we have
u n ( t ) u ( t ) = I α v ( s ) ( t ) + ( t a ) γ 1 Λ Γ ( γ ) I α v ( s ) ( b ) i = 1 m δ i I α + φ i v ( s ) ( ξ i ) .
Hence, u F ( x ) .
Next we show that there exists δ < 1 (defined by (10)), such that
H d ( F ( x ) , F ( x ¯ ) ) δ x x ¯ for each x , x ¯ C 2 ( [ a , b ] , R ) .
Let x , x ¯ C 2 ( [ a , b ] , R ) and h 1 F ( x ) . Then there exists v 1 ( t ) F ( t , x ( t ) ) such that, for each t [ a , b ] ,
h 1 ( t ) = I α v 1 ( s ) ( t ) + ( t a ) γ 1 Λ Γ ( γ ) I α v 1 ( s ) ( b ) i = 1 m δ i I α + φ i v 1 ( s ) ( ξ i ) .
By ( A 2 ) , we have
H d ( F ( t , x ) , F ( t , x ¯ ) ) m ( t ) | x ( t ) x ¯ ( t ) | .
Therefore, there exists w F ( t , x ¯ ( t ) ) , such that
| v 1 ( t ) w | m ( t ) | x ( t ) x ¯ ( t ) | , t [ a , b ] .
Define U : J P ( R ) by
U ( t ) = { w R : | v 1 ( t ) w | m ( t ) | x ( t ) x ¯ ( t ) | } .
Since the multivalued operator U F ( t , x ¯ ( t ) ) is measurable (Proposition III.4 [26]), there exists a function v 2 ( t ) which is a measurable selection for U. Therefore, v 2 ( t ) F ( t , x ¯ ( t ) ) and for each t [ a , b ] , we have | v 1 ( t ) v 2 ( t ) | m ( t ) | x ( t ) x ¯ ( t ) | .
For each t [ a , b ] , let us define
h 2 ( t ) = I α v 2 ( s ) ( t ) + ( t a ) γ 1 Λ Γ ( γ ) I α v 2 ( s ) ( b ) i = 1 m δ i I α + φ i v 2 ( s ) ( ξ i ) .
Thus,
| h 1 ( t ) h 2 ( t ) | I α ( | v 1 ( s ) v 2 ( s ) | ) ( t ) + ( b a ) γ 1 Λ Γ ( γ ) I α ( | v 1 ( s ) v 2 ( s ) | ) ( b ) + i = 1 m δ i I α + φ i ( | v 1 ( s ) v 2 ( s ) | ) ( ξ i ) ( b a ) α Γ ( α + 1 ) + ( b a ) γ 1 | Λ | Γ ( γ ) ( b a ) α Γ ( α + 1 ) + i = 1 m | δ i | ( ξ i a ) α + φ i Γ ( α + φ i + 1 ) m x x ¯ .
Hence,
h 1 h 2 ( b a ) α Γ ( α + 1 ) + ( b a ) γ 1 | Λ | Γ ( γ ) ( b a ) α Γ ( α + 1 ) + i = 1 m | δ i | ( ξ i a ) α + φ i Γ ( α + φ i + 1 ) m x x ¯ .
Analogously, interchanging the roles of x and x ¯ , we obtain
H d ( F ( x ) , F ( x ¯ ) ) ( b a ) α Γ ( α + 1 ) + ( b a ) γ 1 | Λ | Γ ( γ ) ( b a ) α Γ ( α + 1 ) + i = 1 m | δ i | ( ξ i a ) α + φ i Γ ( α + φ i + 1 ) m x x ¯ .
Therefore, F is a contraction. Therefore, it follows by Nadler theorem [24] that F has a fixed point x which is a solution of (3) and (4). This completes the proof. □

4. Examples

Example 1.
Consider the nonlocal boundary value problem with Hilfer fractional differential inclusion of the form:
H D 3 2 , 1 2 x ( t ) F ( t , x ( t ) ) , t 1 4 , 3 2 , x 1 4 = 0 , x 3 2 = 2 3 I 1 2 x 1 2 + 3 4 I 3 2 x 3 4 + 4 5 I 5 2 x 1 + 5 6 I 7 2 x 5 4 .
Here, α = 3 / 2 , β = 1 / 2 , a = 1 / 4 , b = 3 / 2 , m = 4 , δ i = ( i + 1 ) / ( i + 2 ) , φ i = ( 2 i 1 ) / 2 and ξ i = ( i + 1 ) / 4 , i = 1 , 2 , 3 , 4 . From all information, we can find that γ = 7 / 4 , | Λ | 1.058845070 and
( b a ) α Γ ( α + 1 ) + ( b a ) γ 1 | Λ | Γ ( γ ) ( b a ) α Γ ( α + 1 ) + i = 1 m | δ i | ( ξ i a ) α + φ i Γ ( α + φ i + 1 ) 1.116843293 .
( i ) Let the set F ( t , x ) be given by
F ( t , x ) = 4 4 t + 7 x 2 + 2 | x | 1 + | x | + 2 3 , 4 4 t + 11 e | x | sin x 2 x 4 + 1 + 1 4 .
From (12), we have
F ( t , x ) P 4 4 t + 7 x + 5 3 ,
which can be set p ( t ) = 4 / ( 4 t + 7 ) and ψ ( x ) = x + ( 5 / 3 ) . Then we have p = 1 / 2 and there exists a constant M > 2.107672935 satisfying inequality in ( H 3 ) . Therefore, by the benefit of Theorem (1), the boundary value problem (11) with (12) has at least one solution on [ 1 / 4 , 3 / 2 ] .
( i i ) The set F ( t , x ) is defined by
F ( t , x ) = t + 1 4 tan 1 | x | + 1 5 , 1 + sin t 6 x 2 + | x | 2 ( 1 + | x | ) + 1 4 .
It is clear that F defined in (13) is measurable for all x R . Next, we see that
H d ( F ( t , x ) , F ( t , y ) ) 1 4 ( t + 1 ) | x y | ,
for all x , y R . Setting m ( t ) = ( t + 1 ) / 4 , we have m = 5 / 8 and d ( 0 , F ( t , 0 ) ) m ( t ) for all t [ 1 / 4 , 3 / 2 ] . Then we have
( b a ) α Γ ( α + 1 ) + ( b a ) γ 1 | Λ | Γ ( γ ) ( b a ) α Γ ( α + 1 ) + i = 1 m | δ i | ( ξ i a ) α + φ i Γ ( α + φ i + 1 ) m 0.6980270581 < 1 .
Hence, by applying Theorem 2, we deduce that the boundary value problem (11) with (13) has at least one solution on [ 1 / 4 , 3 / 2 ] .

5. Conclusions

We initiated in this work the study of boundary value problems for Hilfer fractional differential inclusions. We obtain existence results by considering the cases when the multivalued map has convex or nonconvex values. In the case of convex multivalued map the existence result is proved via nonlinear alternative of Leray–Schauder type, while in the case of nonconvex multivalued map the existence result is proved by using a fixed point theorem due to Nadler and concern multivalued contractions. The obtained results are well illustrated by numerical examples. We emphasize that our results are new in the context of Hilfer fractional differential inclusions and contribute significantly to the existing literature on the topic.

Author Contributions

Conceptualization, S.K.N. and J.T.; methodology, A.W., S.K.N. and J.T.; formal analysis, A.W., S.K.N. and J.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by King Mongkut’s University of Technology North Bangkok. Contract no. KMUTNB-62-KNOW-42.

Acknowledgments

The authors thanks the referees for their helpful suggestions which improved final version of this paper.

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: Amsterdam, The Netherlands, 2006. [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. Zhou, Y. Basic Theory of Fractional Differential Equations; World Scientific: Singapore, 2014. [Google Scholar]
  8. Martinez-Garcia, M.; Kalawsky, R.S.; Gordon, T.; Smith, T.; Flemisch, Q.M.F. Communication and Interaction with semiautonomous ground vehicles by force control steering. IEEE Trans. Cybern. 2020. [Google Scholar] [CrossRef] [PubMed]
  9. Martinez-Garcia, M.; Zhang, Y.; Gordon, T. Memory pattern identification for feedback tracking control in human-machine systems. Hum. Factors 2019. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Hilfer, R. (Ed.) Applications of Fractional Calculus in Physics; World Scientific: Singapore, 2000. [Google Scholar]
  11. Hilfer, R. Experimental evidence for fractional time evolution in glass forming materials. J. Chem. Phys. 2002, 284, 399–408. [Google Scholar] [CrossRef]
  12. Hilfer, R.; Luchko, Y.; Tomovski, Z. Operational method for the solution of fractional differential equations with generalized Riemann-Liouvill fractional derivatives. Frac. Calc. Appl. Anal. 2009, 12, 299–318. [Google Scholar]
  13. Furati, K.M.; Kassim, N.D.; Tatar, N.E. Existence and uniqueness for a problem involving Hilfer fractional derivative. Comput. Math. Appl. 2012, 64, 1616–1626. [Google Scholar] [CrossRef] [Green Version]
  14. Gu, H.; Trujillo, J.J. Existence of mild solution for evolution equation with Hilfer fractional derivative. Appl. Math. Comput. 2015, 257, 344–354. [Google Scholar] [CrossRef]
  15. Wang, J.; Zhang, Y. Nonlocal initial value problems for differential equations with Hilfer fractional derivative. Appl. Math. Comput. 2015, 266, 850–859. [Google Scholar] [CrossRef]
  16. Asawasamrit, S.; Kijjathanakorn, A.; Ntouyas, S.K.; Tariboon, J. Nonlocal boundary value problems for Hilfer fractional differential equations. Bull. Korean Math. Soc. 2018, 55, 1639–1657. [Google Scholar]
  17. Wongcharoen, A.; Ntouyas, S.K.; Tariboon, J. Nonlocal boundary value problems for Hilfer type pantograph fractional differential equations and inclusions. Adv. Differ. Equ. 2020, 2020, 279. [Google Scholar] [CrossRef]
  18. Wongcharoen, A.; Ntouyas, S.K.; Tariboon, J. On coupled systems for Hilfer fractional differential equations with nonlocal integral boundary conditions. J. Math. 2020, 2020, 2875152. [Google Scholar] [CrossRef]
  19. Deimling, K. Multivalued Differential Equations; Walter De Gruyter: Berlin, Germany; New York, NY, USA, 1992. [Google Scholar]
  20. Hu, S.; Papageorgiou, N. Handbook of Multivalued Analysis, Theory I; Kluwer: Dordrecht, The Netherlands, 1997. [Google Scholar]
  21. Ahmad, B.; Alsaedi, A.; Ntouyas, S.K.; Tariboon, J. Hadamard-Type Fractional Differential Equations, Inclusions and Inequalities; Springer: Cham, Switzerland, 2017. [Google Scholar]
  22. Granas, A.; Dugundji, J. Fixed Point Theory; Springer: New York, NY, USA, 2005. [Google Scholar]
  23. 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]
  24. Nadler, S.B., Jr. Multivalued contraction mappings. Pasific J. Math. 1969, 30, 475–488. [Google Scholar] [CrossRef]
  25. Kisielewicz, M. Differential Inclusions and Optimal Control; Kluwer: Dordrecht, The Netherlands, 1991. [Google Scholar]
  26. Castaing, C.; Valadier, M. Convex Analysis and Measurable Multifunctions; Lecture Notes in Mathematics 580; Springer: Berlin/Heidelberg, Germany; New York, NY, USA, 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

Wongcharoen, A.; Ntouyas, S.K.; Tariboon, J. Boundary Value Problems for Hilfer Fractional Differential Inclusions with Nonlocal Integral Boundary Conditions. Mathematics 2020, 8, 1905. https://0-doi-org.brum.beds.ac.uk/10.3390/math8111905

AMA Style

Wongcharoen A, Ntouyas SK, Tariboon J. Boundary Value Problems for Hilfer Fractional Differential Inclusions with Nonlocal Integral Boundary Conditions. Mathematics. 2020; 8(11):1905. https://0-doi-org.brum.beds.ac.uk/10.3390/math8111905

Chicago/Turabian Style

Wongcharoen, Athasit, Sotiris K. Ntouyas, and Jessada Tariboon. 2020. "Boundary Value Problems for Hilfer Fractional Differential Inclusions with Nonlocal Integral Boundary Conditions" Mathematics 8, no. 11: 1905. https://0-doi-org.brum.beds.ac.uk/10.3390/math8111905

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop