Boundary value problems for Hilfer type sequential fractional differential equations and inclusions involving Riemann–Stieltjes integral multi-strip boundary conditions

In this paper, we study boundary value problems for sequential fractional differential equations and inclusions involving Hilfer fractional derivatives, supplemented with Riemann–Stieltjes integral multi-strip boundary conditions. Existence and uniqueness results are obtained in the single-valued case by using the classical Banach and Krasnosel’skĭi fixed point theorems and the Leray–Schauder nonlinear alternative. In the multi-valued case an existence result is proved by using nonlinear alternative for contractive maps. Examples illustrating our results are also presented.


Introduction and preliminaries
Differential equations of fractional order describe many real world processes more accurately as compared to classical order differential equations. For this reason differential equations of fractional order with initial/boundary conditions have been studied by many researchers. Fractional differential equations arise in lots of engineering and clinical disciplines which include biology, physics, chemistry, economics, signal and image processing, control theory, and so on; see the monographs [1][2][3][4][5][6][7][8].
Several different definitions of fractional integrals and derivatives exist in the literature.
The most popular ones are the Riemann-Liouville fractional derivative of order α > 0 defined for a continuous function by RL D α u(t) := D n I n-α u(t) = 1 (nα) and the Caputo fractional derivative of order α > 0 defined by C D α u(t) := I n-α D n u(t) = 1 (nα) t a (ts) n-α-1 d ds n u(s) ds, n -1 < α < n.
In the above definitions I (·) is the Riemann-Liouville fractional integral of order α > 0 which is defined by (ts) α-1 u(s) ds, n -1 < α < n, provided the right-hand side is point-wise defined on (a, ∞). Other less-known definitions of fractional integrals and derivatives are the Hadamard fractional derivative, the Erdélyi-Kober fractional derivative, and so on. Hilfer in [9] generalized both Riemann-Liouville and Caputo fractional derivatives. This derivative, known as the Hilfer fractional derivative, is defined by where n -1 < α < n, 0 ≤ β ≤ 1, t > a, D = d dt . Fractional differential equations involving Hilfer derivative have many applications, we refer to [10] and the references cited therein. There are actual world occurrences with uncharacteristic dynamics such as atmospheric diffusion of pollution, signals transmissions through strong magnetic fields, the effect of the theory of the profitability of stocks in economic markets, the theoretical simulation of dielectric relaxation in glass forming materials, network traffic, and so on. See [11,12] and the references cited therein. We refer to [13,14] and the references cited therein for some properties of the Hilfer derivative, and to [15][16][17] for initial value problems involving Hilfer fractional derivatives.
In contrast to classical boundary conditions, nonlocal conditions help to formulate the changes happening on certain positions and segments of the given domain [18]. On the other hand, boundary conditions involving derivatives and integrals provide a decent approach to describe non-uniformities on segments of curved boundary structures. Examples include fluid problems [19], biomedical applications [20], engineering applications [21,22], etc.
In [23] the authors initiated the study of nonlocal boundary value problems for Hilfer fractional derivative by studying the boundary value problem of Hilfer-type fractional differential equations with nonlocal integral boundary conditions 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 ∈ [a, b], a ≥ 0, and δ i ∈ R. Several existence and uniqueness results were proved by using a variety of fixed point theorems.
In a series of papers  nonlocal boundary value problems involving Hilfer fractional derivatives were studied with a variety of boundary conditions. Thus, in [7] the authors studied Hilfer-Langevin three-point fractional boundary value problems, in [24] pantograph Hilfer fractional boundary value problems with nonlocal integral boundary conditions were studied, in [25] Hilfer fractional boundary value problems with nonlocal integral integro-multipoint boundary conditions were studied, in [26] Hilfer fractional boundary value problems with nonlocal multipoint, fractional derivative multi-order, and fractional integral multi-order boundary conditions were studied, and in [27] sequential Hilfer fractional boundary value problems with nonlocal integro-multipoint boundary conditions were studied.
In the present paper, motivated by the recent interest in studying boundary value problems involving Hilfer fractional derivative operators, we discuss a new class of sequential Hilfer-type boundary value problems for fractional differential equations involving Riemann-Stieltjes integral multi-strip boundary conditions of the form where H D α,β denotes the Hilfer fractional derivative operator of order α, 1 < α < 2, and parameter β, We establish existence and uniqueness results by using the basic tools from fixed point theory. The existence of a unique solution is proved by using the Banach contraction mapping principle, while in the two existence results we make use of the Krasnosel'skiȋ fixed point theorem and the nonlinear alternative of Leray-Schauder type. The main results are presented in Sect. 2. The obtained results are well illustrated by numerical examples. The results obtained in the present paper are new and significantly contribute to the existing literature on the topic.
We also cover the multi-valued case of problem (1.3) by considering the following sequential inclusion boundary value problem: where F : [a, b] × R → P(R) is a multi-valued function, (P(R) is the family of all nonempty subsets of R). An existence result is proved in Sect. 3 for the sequential Hilfer inclusion boundary value problem (1.4) via nonlinear alternative for contractive maps.

Existence and uniqueness results for problem (1.3)
The following auxiliary lemma, concerning a linear variant of the boundary value problem (1.3), plays a fundamental role in establishing the existence and uniqueness results for the given nonlinear problem.
, R), and Then the function x is a solution of the sequential boundary value problem

2)
if and only if Proof Let x be a solution of the nonlocal sequential boundary value problem (2.2). By Lemma 2.1, applying on both sides of equation (2.2) the fractional integral I α , we have where c 0 and c 1 are arbitrary real constants. From x(a) = 0 we have c 0 = 0. Then we get Substituting the value of c 1 in (2.4), we obtain solution (2.3). The converse can be proven by direct computation. The proof is finished.
It is obvious that the sequential boundary value problem (1.3) has solutions if and only if the operator A has fixed points.
For computational convenience we use the notations: By using well-known fixed point theorems, we prove existence, as well as existence and uniqueness results, for the sequential Hilfer-type boundary value problem (1.3) in the following sections.

Existence and uniqueness result
In our first result we prove the existence of a unique solution for the sequential Hilfer boundary value problem (1.3) via Banach's contraction mapping principle.

Theorem 2.1 Assume that:
( , and x, y ∈ R. If Set sup t∈ [a,b] |f (t, 0)| = M < ∞, and choose r > 0 such that For any x ∈ B r , we have Hence Ax -Ay ≤ (L + 1 ) xy . Consequently, A is a contraction as, by assumption, L + 1 < 1. We deduce, by the Banach contraction mapping principle, that A has a fixed point, which means that the sequential Hilfer boundary value problem (1.3) has a unique solution. This completes the proof.

Existence results
Two existence results are presented in this subsection. The first one is based on the wellknown Krasnosel'skiȋ fixed point theorem [28]. Proof We set sup t∈ [a,b]  For x, y ∈ B ρ , we have Therefore A 1 x + A 2 y ≤ ρ, which shows that A 1 x + A 2 y ∈ B ρ . It is easy to see, using the condition 1 < 1, that A 2 is a contraction mapping. The operator A 1 is continuous since f is continuous. Moreover, the uniform boundedness on B ρ of the operator A 1 follows from the relation Now we prove the equicontinuity of the operator A 1 . Let t 1 , t 2 ∈ [a, b] with t 1 < t 2 . Then we have as t 2 -t 1 → 0 and is independent of x. Therefore, the operator A 1 is equicontinuous. By the Arzelá-Ascoli theorem, A 1 is compact on B ρ . Consequently all the assumptions of Krasnosel'skiiȋ's fixed point theorem [28] are satisfied. Hence the sequential Hilfer boundary value problem (1.3) has at least one solution on [a, b]. The proof is completed.
In our second existence result we apply the Leray-Schauder nonlinear alternative [29].

Theorem 2.3 Let 1 < 1.
In addition, we assume that: and consequently, In the next step we will show that bounded sets are mapped by the operator A into equicontinuous sets of C([a, b], R). Let t 1 , t 2 ∈ [a, b] with t 1 < t 2 and x ∈ B r . Then we have as t 2t 1 → 0 and is independent of x ∈ B r . By the Arzelá-Ascoli theorem, the operator A : C([a, b], R) → C ([a, b], R) is completely continuous. Finally we will prove that the set of all solutions to equation x = θ Ax for θ ∈ (0, 1) is bounded.
(i) If the nonlinear function f (t, x) in (2.11) is given by

Existence results for problem (1.4)
In the following we use the notation P cp,c (X) = {Y ∈ P(X) : Y is compact and convex} for a normed space (X, · ). For details on multi-valued analysis, we refer to the monographs [7,[30][31][32]. F(t, x(t)) a.e. on [a, b] and Theorem 3.1 Assume that 1 < 1. In addition, we suppose that: for all x ∈ R with x ≤ α and for a.e. t ∈ [a, b];  (C([a, b], R)) by Next we introduce the single-valued operator A : C([a, b], R) − → C ([a, b], R) and the multi-valued operator B : C([a, b], R) − → P (C([a, b], R)) by and . (3.4) Observe that N = A + B. To prove our existence result, we make use of the nonlinear alternative for contractive maps [33,Corollary 3.8].
Step 1: First we show that A is a contraction on C([a, b], R). We have, for x, y ∈ C([a, b], R), Hence, since 1 < 1, the operator A is a contraction.
Step 2: for almost all t ∈ [a, b]. Let 0 ≤ θ ≤ 1. Then we have The fact that F has convex values implies that S F,x is convex, and thus θ v 1 (s) + (1θ )v 2 (s) ∈ S F,x . Hence θ z 1 + (1θ )z 2 ∈ B(x), which proves that the operator B is convexvalued.
Step 3: Now we will prove that the operator B is compact and upper semicontinuous. First, we prove that bounded sets are mapped by B into bounded sets in C([a, b], R). Let We have, for t ∈ [a, b], Thus, Secondly we prove that bounded sets are mapped by the operator B into equicontinuous sets. Let t 1 , t 2 ∈ [a, b] with t 1 < t 2 and x ∈ B r . Then, for each h ∈ B(x), we obtain as t 2t 1 → 0 and is independent of x ∈ B r . By the Arzelá-Ascoli theorem the operator B : C([a, b], R) → P (C([a, b], R)) is completely continuous. By [30, Proposition 1.2] we know that a completely continuous operator is upper semicontinuous if it has a closed graph. Thus in the following we will prove that the operator B has a closed graph.
Let x n → x * , h n ∈ B(x n ) and h n → h * . We shall show that h * ∈ B(x * ). Since h n ∈ B(x n ), there exists v n ∈ S F,x n such that, for each t ∈ [a, b], It is suffices to show that there exists v * ∈ S F,x * such that, Let us consider the linear operator : Observe that h n (t)h * (t) = I α v n (t)v * (t) + (ta) γ -1 |λ| b a I α v n (t)v * (t) dH(s) as n → ∞. Thus, it follows by a closed graph lemma [34] that • S F is a closed graph operator. Moreover, h n (t) ∈ (S F,x n ). Since x n → x * , we have that for some v * ∈ S F,x * . Thus the operator B has a closed graph. Thus the operators A and B satisfy all the conditions of the nonlinear alternative for contractive maps [33,Corollary 3.8]. It implies that either (i) N has a fixed point in [a, b] or (ii) there is a point x ∈ ∂B M = {x ∈ C([a, b], R) : x ≤ M} and θ ∈ (0, 1) with x = θN (x).

Conclusion
This paper discussed a new class of boundary value problems for sequential fractional differential equations and inclusions involving Hilfer fractional derivatives, supplemented with Riemann-Stieltjes integral multi-strip boundary conditions. Existence and uniqueness results are established in the single-valued case by using the classical Banach and Krasnosel'skiȋ fixed point theorems and the nonlinear alternative of Leray-Schauder type. In the multi-valued case, an existence result is proved by using nonlinear alternative for contractive maps. By suitable numerical examples we verified the derived analysis. The results of the present paper are new and significantly contribute to the existing literature on the topic.