On (k,ψ)-Hilfer Fractional Differential Equations and Inclusions with Mixed (k,ψ)-Derivative and Integral Boundary Conditions

Abstract

In this paper we study single-valued and multi-valued (k,ψ)-Hilfer-type boundary value problems of fractional order in (1,2], subject to nonlocal boundary conditions involving (k,ψ)-Hilfer-type derivative and integral operators. The results for single-valued case are established by using Banach and Krasnosel’skiĭ fixed point theorems as well as Leray–Schauder nonlinear alternative. In the multi-valued case, we establish an existence result for the convex valued right-hand side of the inclusion via Leray–Schauder nonlinear alternative for multi-valued maps, while the second one when the right-hand side has non-convex values is obtained by applying Covitz–Nadler fixed point theorem for multi-valued contractions. Numerical examples illustrating the obtained theoretical results are also presented.

1. Introduction

Fractional calculus is concerned with integral and derivative operators of non-integer order and arises in many engineering and scientific fields such as physics, chemistry, mathematical biology, mechanics, and so forth; see the monographs [1,2,3,4,5,6,7,8,9]. Usually, fractional integral operators are used to define fractional derivative operators. Many kinds of fractional derivative operators such as Rieman–Liouville, Caputo, Hadamard, Katugampola, Hilfer, etc., are proposed in the literature. Certain forms of fractional operators include definitions of other fractional operators. For example, the concept of generalized fractional derivatives and integrals introduced by Katugampola in [10,11] includes Riemann–Liouville and Hadamard fractional derivatives. The Hilfer fractional derivative operator [12] includes Rieman–Liouville and Caputo fractional derivative operators. The $\psi$-fractional derivative operator [13] unifies Caputo, Caputo–Hadamard and Caputo–Erdélyi-Kober fractional derivative operators. A wide class of fractional operators is covered by the $(k,\psi )$-Hilfer fractional derivative operator introduced in [14,15].

In [14], the authors, by applying Banach’s fixed point theorem, proved the existence of a unique solution for a nonlinear initial value problem involving $(k,\psi )$-Hilfer-type fractional derivative operator. Tariboon et al. [16] studied the existence and uniqueness of solutions for $(k,\psi )$-Hilfer fractional differential equations and inclusions with multi-point boundary conditions.

To enrich the literature in this new research topic, which is very limited at the moment, we continue in the present paper the study of boundary value problems involving $(k,\psi )$-Hilfer-type fractional derivatives in the order of $(1,2],$ supplemented with nonlocal boundary conditions involving $(k,\psi )$-Hilfer-type derivative and integral operators of the form

(1)$\left\{\begin{array}{c}{\displaystyle {}^{k,H}{D}^{\alpha ,\beta ;\psi }u\left(t\right)=\mathfrak{f}(t,u\left(t\right)),t\in (a,b],}\hfill \\ {\displaystyle u\left(a\right)=0,u\left(b\right)=\lambda {}^{k,H}{D}^{p,q;\psi }u\left(\eta \right)+\mu {}^k{\mathfrak{I}}^{v,\psi }u\left(\sigma \right),}\hfill \end{array}\right.$

where ${}^{k,H}{D}^{\alpha ,\beta ;\psi }$ denotes the $(k,\psi )$-Hilfer-type fractional derivative of order $\alpha$, $1<\alpha <2$ and parameter $\beta$, $0\le \beta \le 1$, $k>0$, $\mathfrak{f}:[a,b]\times \mathbb{R}\to \mathbb{R}$ is a continuous function, ${}^{k,H}{D}^{p,q;\psi }$ denotes the $(k,\psi )$-Hilfer-type fractional derivative of order p, $1<p<2$ and parameter q, $0\le q\le 1$, $p<\alpha$, ${}^k{\mathfrak{I}}^{v,\psi }$ is the $(k,\psi )$-Riemann–Liouville fractional integral of order $v>0,\lambda ,\mu \in \mathbb{R},$ and $a<\xi ,\sigma <b.$ Existence and uniqueness will be established, by using Banach’s and Krasnosel’skiĭ’s fixed point theorems, as well as the Laray–Schauder nonlinear alternative.

We also study the corresponding multi-valued case of the problem (1) given by

(2)$\left\{\begin{array}{c}{\displaystyle {}^{k,H}{D}^{\alpha ,\beta ;\psi }u\left(t\right)\in \mathfrak{F}(t,u\left(t\right)),t\in (a,b],}\hfill \\ {\displaystyle u\left(a\right)=0,u\left(b\right)=\lambda {}^{k,H}{D}^{p,q;\psi }u\left(\eta \right)+\mu {}^k{\mathfrak{I}}^{v,\psi }u\left(\sigma \right),}\hfill \end{array}\right.$

in which $\mathfrak{F}:[a,b]\times \mathbb{R}\to \mathcal{P}\left(\mathbb{R}\right)$ is a multivalued map ($\mathcal{P}\left(\mathbb{R}\right)$ denotes the family of all nonempty subsets of $\mathbb{R}$) while the other parameters are the same as defined in problem (1). The existence results for the problem (2) associated with convex and non-convex cases of the multi-valued map $\mathfrak{F}(t,u(t\left)\right)$ will be obtained by using the Laray–Schauder nonlinear alternative for multi-valued maps and the Covitz–Nadler fixed point theorem for multi-valued contractions, respectively.

The fixed point theory provides an excellent approach to establish the existence theory for initial and boundary value problems. For some recent publications on this branch of mathematical analysis, we refer the reader to some recent books [17,18,19].

The content in the rest of paper is arranged as follows. Section 2 contains some necessary definitions and lemmas, while Section 3 is concerned with an auxiliary lemma which enables us to transform the nonlinear $(k,\psi )$-Hilfer type boundary value problem (1) into an equivalent fixed point problem. The main results for the problems (1) and (2) are presented in Section 4 and Section 5, respectively. Finally, in Section 6, numerical examples illustrating the obtained theoretical results are presented.

2. Preliminaries

In this section we introduce some definitions and lemmas that will be used throughout the paper.

Definition 1

([20]). Let $\mathfrak{h}\in L^1([a,b],\mathbb{R})$ and $k,\alpha \in {\mathbb{R}}^{+}$. Then, the k-Riemann–Liouville fractional derivative of order α of the function $\mathfrak{h}$ is given by

(3)${}^k{\mathfrak{I}}_{a+}^{\alpha }\mathfrak{h}\left(t\right)=\frac{1}{k{\Gamma }_k\left(\alpha \right)}{\int }_a^t{(t-w)}^{\frac{\alpha }{k}-1}\mathfrak{h}\left(w\right)dw,$

where ${\Gamma }_k$ is the k-Gamma function for $z\in \mathbb{C}$ with $\Re \left(z\right)>0$ and $k\in \mathbb{R},k>0$ which is defined in [21] by

${\Gamma }_k\left(z\right)={\int }_0^{+\infty }{s}^{z-1}{e}^{-\frac{s^k}{k}}ds.$

The following relations are well known.

$\Gamma \left(\theta \right)=\underset{k\to 1}{lim}{\Gamma }_k\left(\theta \right),{\Gamma }_k\left(\theta \right)=k^{\frac{\theta }{k}-1}\Gamma \left(\frac{\theta }{k}\right)\mathrm{and}{\Gamma }_k(\theta +k)=\theta {\Gamma }_k\left(\theta \right).$

Definition 2

([22]). Let $\mathfrak{h}\in L^1([a,b],\mathbb{R})$ and $k,\alpha \in {\mathbb{R}}^{+}.$ Then the k-Riemann–Liouville fractional derivative of order α of the function $\mathfrak{h}$ is given by

(4)${}^{k,RL}{D}_{a+}^{\alpha }\mathfrak{h}\left(t\right)={\left(k\frac{d}{dt}\right)}^n{}^k{\mathfrak{I}}_{a+}^{nk-\alpha }\mathfrak{h}\left(t\right),n=⌈\frac{\alpha }{k}⌉,$

where $⌈\frac{\alpha }{k}⌉$ is the ceiling function of $\frac{\alpha }{k}.$

Definition 3

([2]). Let $\mathfrak{h}\in L^1([a,b],\mathbb{R})$ and an increasing function $\psi :[a,b]\to \mathbb{R}$ with ${\psi }^{\prime }\left(t\right)\ne 0$ for all $t\in [a,b].$ Then the ψ-Riemann–Liouville fractional integral of the function $\mathfrak{h}$ is given by

(5)${\mathfrak{I}}^{\alpha ;\psi }\mathfrak{h}\left(t\right)=\frac{1}{{\Gamma }_k\left(\alpha \right)}{\int }_a^t{\psi }^{\prime }\left(u\right){(\psi \left(t\right)-\psi \left(u\right))}^{\alpha -1}\mathfrak{h}\left(u\right)du.$

Definition 4.

Let $n-1<\alpha \le n,\psi \in C^n([a,b],\mathbb{R})$ is an increasing function with ${\psi }^{\prime }\left(t\right)\ne 0,t\in [a,b],$ and $\mathfrak{h}\in C\left(\right[a,b],\mathbb{R}).$

• (a) 

  The ψ-Riemann–Liouville fractional derivative of the function $\mathfrak{h}$ of order α is given in [2] as

  (6) ${}^{RL}{D}^{\alpha ;\psi }\mathfrak{h}\left(t\right)={\left(\frac{1}{{\psi }^{\prime }\left(t\right)}\frac{d}{dt}\right)}^n{\mathfrak{I}}_{a+}^{n-\alpha ;\psi }\mathfrak{h}\left(t\right).$

• (b) 

  The ψ-Caputo fractional derivative of the function $\mathfrak{h}$ of order α is defined in [13] as

  (7) ${}^C{D}^{\alpha ;\psi }\mathfrak{h}\left(t\right)={\mathfrak{I}}_{a+}^{n-\alpha ;\psi }{\left(\frac{1}{{\psi }^{\prime }\left(t\right)}\frac{d}{dt}\right)}^n\mathfrak{h}\left(t\right).$

• (c) 

  The ψ-Hilfer fractional derivative of the function $\mathfrak{h}\in C\left(\right[a,b],\mathbb{R})$ of order $\alpha \in (n-1,n]$ and type $\beta \in [0,1]$ is defined in [23] as

  (8) ${}^H{D}^{\alpha ,\beta ;\psi }\mathfrak{h}\left(t\right)={\mathfrak{I}}_{a+}^{\beta (n-\alpha );\psi }{\left(\frac{1}{{\psi }^{\prime }\left(t\right)}\frac{d}{dt}\right)}^n{\mathfrak{I}}_{a+}^{(1-\beta )(n-\alpha );\psi }\mathfrak{h}\left(t\right).$

Definition 5

([24]). Let $\psi :[a,b]\to \mathbb{R}$ be an increasing function with ${\psi }^{\prime }\left(t\right)\ne 0$ for all $t\in [a,b]$. Then the $(k,\psi )$-Riemann–Liouville fractional integral of order $\alpha >0$ ($\alpha \in \mathbb{R}$) of a function $\mathfrak{h}\in L^1([a,b],\mathbb{R})$ is given by

(9)${}^k{\mathfrak{I}}_{a+}^{\alpha ;\psi }\mathfrak{h}\left(t\right)=\frac{1}{k{\Gamma }_k\left(\alpha \right)}{\int }_a^t{\psi }^{\prime }\left(u\right){(\psi \left(t\right)-\psi \left(u\right))}^{\frac{\alpha }{k}-1}\mathfrak{h}\left(u\right)du,k>0.$

Definition 6

([14]). Let $\alpha ,k\in {\mathbb{R}}^{+}=(0,+\infty ),\beta \in [0,1],\psi \in C^n([a,b],\mathbb{R})$ is an increasing function with ${\psi }^{\prime }\left(t\right)\ne 0,t\in [a,b]$ and $\mathfrak{h}\in C^n([a,b],\mathbb{R}).$ Then the $(k,\psi )$-Hilfer fractional derivative of the function $\mathfrak{h}$ of order α and type $\beta ,$ is defined by

(10)${}^{k,H}{D}^{\alpha ,\beta ;\psi }\mathfrak{h}\left(t\right)={}^k{\mathfrak{I}}_{a+}^{\beta (nk-\alpha );\psi }{\left(\frac{k}{{\psi }^{\prime }\left(t\right)}\frac{d}{dt}\right)}^n{}^k{\mathfrak{I}}_{a+}^{(1-\beta )(nk-\alpha );\psi }\mathfrak{h}\left(t\right),n=⌈\frac{\alpha }{k}⌉.$

Remark 1.

Observe that the $(k,\psi )$-Hilfer-type fractional derivative can be expressed in terms of $(k,\psi )$-Riemann–Liouville fractional derivative as

$\begin{array}{ccc}\hfill {}^{k,H}{D}^{\alpha ,\beta ;\psi }\mathfrak{h}\left(t\right)& =& {}^k{\mathfrak{I}}_{a+}^{{\theta }_k-\alpha ;\psi }{\left(\frac{k}{{\psi }^{\prime }\left(t\right)}\frac{d}{dt}\right)}^n{}^k{\mathfrak{I}}_{a+}^{nk-{\theta }_k;\psi }\mathfrak{h}\left(t\right)\hfill \\ & =& {}^k{\mathfrak{I}}_{a+}^{{\theta }_k-\alpha ;\psi }\left({}^{k,RL}{D}^{{\theta }_k;\psi }\mathfrak{h}\right)\left(t\right),\hfill \end{array}$

where ${\theta }_k=\alpha +\beta (nk-\alpha ),\beta (nk-\alpha )={\theta }_k-\alpha$ and $(1-\beta )(nk-\alpha )=nk-{\theta }_k,\beta \in [0,1].$ Note that $n-1<\frac{{\theta }_k}{k}\le n$ when $n-1<\frac{\alpha }{k}\le n.$

We recall now some useful lemmas.

Lemma 1

([14]). Assume that $\mathfrak{h}\in C^n([a,b],\mathbb{R})$ and ${}^k{\mathfrak{I}}_{a+}^{nk-\mu ;\psi }\mathfrak{h}\in C^n([a,b],\mathbb{R})$ with $\mu ,k\in {\mathbb{R}}^{+}=(0,+\infty )$ and $n=⌈\frac{\mu }{k}⌉.$ Then

${}^k{\mathfrak{I}}^{\mu ;\psi }\left({}^{k,RL}{D}^{\mu ;\psi }\mathfrak{h}\left(t\right)\right)=\mathfrak{h}\left(t\right)-\sum _{j=1}^n\frac{{(\psi \left(t\right)-\psi \left(a\right))}^{\frac{\mu }{k}-j}}{{\Gamma }_k(\mu -jk+k)}{\left[{\left(\frac{k}{{\psi }^{\prime }\left(t\right)}\frac{d}{dt}\right)}^{n-j}{}^k{\mathfrak{I}}_{a+}^{nk-\mu ;\psi }\mathfrak{h}\left(t\right)\right]}_{z=a}.$

Lemma 2

([14]). Let ${\theta }_k=\alpha +\beta (k-\alpha )$ with $\alpha ,k\in {\mathbb{R}}^{+}=(0,+\infty ),\alpha <k,\beta \in [0,1].$ Then

${}^k{\mathfrak{I}}^{{\theta }_k;\psi }\left({}^{k,RL}{D}^{{\theta }_k;\psi }\mathfrak{h}\right)\left(t\right)={}^k{\mathfrak{I}}^{\alpha ;\psi }\left({}^{k,H}{D}^{\alpha ,\beta ;\psi }\mathfrak{h}\right)\left(t\right),\mathfrak{h}\in C^n([a,b],\mathbb{R}).$

Lemma 3

([14]). Let $\zeta ,k\in {\mathbb{R}}^{+}=(0,+\infty )$ and $\eta \in \mathbb{R}$ such that $\frac{\eta }{k}>-1.$ Then

$\left(i\right).{}^k{\mathcal{I}}^{\zeta ,\psi }{(\psi \left(t\right)-\psi \left(a\right))}^{\frac{\eta }{k}}=\frac{{\Gamma }_k(\eta +k)}{{\Gamma }_k(\eta +k+\zeta )}{(\psi \left(t\right)-\psi \left(a\right))}^{\frac{\eta +\zeta }{k}}.$

$\left(ii\right).{}^k{\mathcal{D}}^{\zeta ,\psi }{(\psi \left(t\right)-\psi \left(a\right))}^{\frac{\eta }{k}}=\frac{{\Gamma }_k(\eta +k)}{{\Gamma }_k(\eta +k-\zeta )}{(\psi \left(t\right)-\psi \left(a\right))}^{\frac{\eta -\zeta }{k}}.$

3. An Auxiliary Result

In this section, an auxiliary result dealing with the linear variant the problem (1) is presented.

Lemma 4.

Let $g\in C(a,b)\cap L^1(a,b)$ (see [25,26]) and

(11) $\begin{array}{ccc}\hfill \Omega :& =& \frac{{(\psi \left(b\right)-\psi \left(a\right))}^{\frac{{\theta }_k}{k}-1}}{{\Gamma }_k\left({\theta }_k\right)}-\lambda {\Gamma }_k\left({\theta }_k\right)\frac{{(\psi \left(\eta \right)-\psi \left(a\right))}^{\frac{{\theta }_k-p}{k}-1}}{{\Gamma }_k({\theta }_k-p)}\hfill \\ & & -\mu {\Gamma }_k\left({\theta }_k\right)\frac{{(\psi \left(\sigma \right)-\psi \left(a\right))}^{\frac{{\theta }_k+v}{k}-1}}{{\Gamma }_k({\theta }_k+v)}\ne 0.\hfill \end{array}$

Then, the function $u\in C^2([a,b],\mathbb{R})$ is a solution of the $(k,\psi )$-Hilfer-type nonlocal fractional boundary value problem

(12) $\left\{\begin{array}{c}{\displaystyle {}^{k,H}{D}^{\alpha ,\beta ;\psi }u\left(t\right)=g\left(t\right),t\in (a,b],k>0,1<\alpha \le 2,\beta \in [0,1],}\hfill \\ {\displaystyle u\left(a\right)=0,u\left(b\right)=\lambda {}^{k,H}{D}^{p,q;\psi }u\left(\eta \right)+\mu {}^k{\mathfrak{I}}^{v,\psi }u\left(\sigma \right),}\hfill \end{array}\right.$

if and only if

(13) $u\left(t\right)={}^k{\mathfrak{I}}^{\alpha ;\psi }g\left(t\right)+\frac{{(\psi \left(t\right)-\psi \left(a\right))}^{\frac{{\theta }_k}{k}-1}}{\Omega {\Gamma }_k\left({\theta }_k\right)}\left[\lambda {}^k{\mathfrak{I}}^{\alpha -p;\psi }g\left(\eta \right)+\mu {}^k{\mathfrak{I}}^{v+\alpha ,\psi }g\left(\sigma \right)-{}^k{\mathfrak{I}}^{\alpha ;\psi }g\left(b\right)\right],$

where ${\theta }_k=\alpha +\beta (2k-\alpha ).$

Proof. 

Assume that u is a solution of the $(k,\psi )$-Hilfer-type nonlocal fractional boundary value problem (12). Operating on both sides of equation in (12), the fractional integral ${}^k{\mathfrak{I}}^{\alpha ;\psi }$ and using Lemmas 2 and 1, we obtain

$\begin{array}{ccc}\hfill {}^k{\mathfrak{I}}^{\alpha ;\psi }\left({}^{k,H}{D}^{\alpha ,\beta ;\psi }u\right)\left(t\right)& =& {}^k{\mathfrak{I}}^{{\theta }_k;\psi }\left({}^{k,RL}{D}^{{\theta }_k;\psi }u\right)\left(t\right)\hfill \\ & =& u\left(t\right)-\frac{{(\psi \left(t\right)-\psi \left(a\right))}^{\frac{{\theta }_k}{k}-1}}{{\Gamma }_k\left({\theta }_k\right)}{\left[\left(\frac{k}{{\psi }^{\prime }\left(t\right)}\frac{d}{dt}\right){}^k{\mathfrak{I}}^{2k-{\theta }_k;\psi }u\left(t\right)\right]}_{w=a}\hfill \\ & & -\frac{{(\psi \left(t\right)-\psi \left(a\right))}^{\frac{{\theta }_k}{k}-2}}{{\Gamma }_k({\theta }_k-k)}{\left[{}^k{\mathfrak{I}}^{2k-{\theta }_k;\psi }u\left(t\right)\right]}_{w=a}.\hfill \end{array}$

Consequently

(14)$u\left(t\right)={}^k{\mathfrak{I}}^{\alpha ;\psi }g\left(t\right)+c_0\frac{{(\psi \left(t\right)-\psi \left(a\right))}^{\frac{{\theta }_k}{k}-1}}{{\Gamma }_k\left({\theta }_k\right)}+c_1\frac{{(\psi \left(t\right)-\psi \left(a\right))}^{\frac{{\theta }_k}{k}-2}}{{\Gamma }_k({\theta }_k-k)},$

where

$c_0={\left[\left(\frac{k}{{\psi }^{\prime }\left(t\right)}\frac{d}{dt}\right){}^k{\mathfrak{I}}^{2k-{\theta }_k;\psi }u\left(t\right)\right]}_{w=a},c_1={\left[{}^k{\mathfrak{I}}^{2k-{\theta }_k;\psi }u\left(t\right)\right]}_{w=a}.$

By the condition $u\left(a\right)=0$, we find that $c_1=0$ as $\frac{{\theta }_k}{k}-2<0$ by Remark 1. By using Lemma 3, we obtain

(15)${}^k{D}^{p,q;\psi }{(\psi \left(t\right)-\psi \left(a\right))}^{\frac{{\theta }_k}{k}-1}=\frac{{\Gamma }_k\left({\theta }_k\right)}{{\Gamma }_k({\theta }_k-p)}{(\psi \left(t\right)-\psi \left(a\right))}^{\frac{{\theta }_k-p}{k}-1}.$

and

(16)${}^k{\mathfrak{I}}^{v,\psi }{(\psi \left(t\right)-\psi \left(a\right))}^{\frac{{\theta }_k}{k}-1}=\frac{{\Gamma }_k\left({\theta }_k\right)}{{\Gamma }_k({\theta }_k+v)}{(\psi \left(t\right)-\psi \left(a\right))}^{\frac{{\theta }_k+v}{k}-1}.$

From (15), (16) and the boundary condition: $u\left(b\right)=\lambda u\left(\xi \right)+\mu {}^k{\mathfrak{I}}^{v,\psi }u\left(\sigma \right)$, we obtain

$c_0=\frac{1}{\Omega }\left[\lambda {}^k{\mathfrak{I}}^{\alpha -p;\psi }g\left(\eta \right)+\mu {}^k{\mathfrak{I}}^{v+\alpha ,\psi }g\left(\sigma \right)-{}^k{\mathfrak{I}}^{\alpha ;\psi }g\left(b\right)\right].$

Replacing $c_0$ and $c_1$ in (14) by their above values, we obtain the solution (13). The converse can be proved easily by direct computation. This finishes the proof. □

4. The Single Valued Problem

Let us begin this section by defining the solution of problem (1).

Definition 7.

A function $u\in C^2([a,b],\mathbb{R})$ possessing the $(k,\psi )$-Hilfer fractional derivative in the sense of Definition 6 is said to be a solution of the $(k,\psi )$-Hilfer-type nonlocal fractional boundary value problem (1) if it satisfies the equations ${\displaystyle {}^{k,H}{D}^{\alpha ,\beta ;\psi }u\left(t\right)=\mathfrak{f}(t,u\left(t\right)),t\in (a,b]}$ with $f\in C\left(\right([a,b],\mathbb{R}),\mathbb{R})$ and the boundary conditions $u\left(a\right)=0,u\left(b\right)=\lambda {}^{k,H}{D}^{p,q;\psi }u\left(\eta \right)+\mu {}^k{\mathfrak{I}}^{v,\psi }u\left(\sigma \right)$.

In view of Lemma 4, we define an operator $\mathcal{A}:C\left(\right[a,b],\mathbb{R})\to C\left(\right[a,b],\mathbb{R})$ by

(17)$\begin{array}{ccc}\hfill \left(\mathcal{A}u\right)\left(t\right)& =& \frac{{(\psi \left(t\right)-\psi \left(a\right))}^{\frac{{\theta }_k}{k}-1}}{\Omega {\Gamma }_k\left({\theta }_k\right)}[\lambda {}^k{\mathfrak{I}}^{\alpha -p;\psi }\mathfrak{f}(\eta ,u\left(\eta \right))+\mu {}^k{\mathfrak{I}}^{v+\alpha ,\psi }\mathfrak{f}(\sigma ,u\left(\sigma \right))\hfill \\ & & -{}^k{\mathfrak{I}}^{\alpha ;\psi }\mathfrak{f}(b,u\left(b\right))]+{}^k{\mathfrak{I}}^{\alpha ;\psi }\mathfrak{f}(t,u\left(t\right)),t\in [a,b],\hfill \end{array}$

where $C\left(\right[a,b],\mathbb{R})$ denotes the Banach space of all continuous real valued functions defined on $[a,b]$ equipped with the sup-norm $\parallel u\parallel ={sup}_{t\in [a,b]}\left|u\left(t\right)\right|.$

Observe that the solutions of the $(k,\psi )$-Hilfer-type nonlocal fractional boundary value problem (1) must be sought among the fixed points of $\mathcal{A}$.

For computational convenience we put:

(18)$\begin{array}{ccc}\hfill \mathbb{G}& =& \frac{{(\psi \left(b\right)-\psi \left(a\right))}^{\frac{\alpha }{k}}}{{\Gamma }_k(\alpha +k)}+\frac{{(\psi \left(b\right)-\psi \left(a\right))}^{\frac{{\theta }_k}{k}-1}}{|\Omega |{\Gamma }_k\left({\theta }_k\right)}[\left|\lambda \right|\frac{{(\psi \left(\eta \right)-\psi \left(a\right))}^{\frac{\alpha -p}{k}}}{{\Gamma }_k(\alpha -p+k)}\hfill \\ & & +\left|\mu \right|\frac{{(\psi \left(\sigma \right)-\psi \left(a\right))}^{\frac{\alpha +v}{k}}}{{\Gamma }_k(\alpha +v+k)}+\frac{{(\psi \left(b\right)-\psi \left(a\right))}^{\frac{\alpha }{k}}}{{\Gamma }_k(\alpha +k)}].\hfill \end{array}$

4.1. Existence of a Unique Solution

In this subsection, we make use of Banach’s fixed point theorem [27] to prove the existence of a unique solution to the problem (1).

Theorem 1.

Let the following condition hold:

• $\left(H_1\right)$

  $\left|\mathfrak{f}\right(t,u)-\mathfrak{f}(t,y\left)\right|\le \mathfrak{L}|u-y|,\mathfrak{L}>0for eacht\in [a,b]andu,y\in \mathbb{R}.$

Then there exists a unique solution to the $(k,\psi )$-Hilfer-type nonlocal fractional boundary value problem (1) on $[a,b],$ provided that

(19) $\mathfrak{L}\mathbb{G}<1,$

where $\mathbb{G}$ is defined by $\left(18\right).$

Proof. 

We verify that the operator $\mathcal{A}$ defined in (17) satisfies the hypothesis of Banach’s fixed point theorem. Letting ${sup}_{t\in [a,b]}\left|\mathfrak{f}(t,0)\right|=\mathbb{M}<+\infty$, we define $B_r=\{u\in C([a,b],\mathbb{R}):\parallel u\parallel \le r\}$ with

(20)$r\ge \frac{{\displaystyle \mathbb{M}\mathbb{G}}}{1-\mathfrak{L}\mathbb{G}}.$

We will first show that $\mathcal{A}{B}_r\subset B_r.$ By the assumptions $\left(H_1\right),$ we have

$\begin{array}{ccc}\hfill \left|\mathfrak{f}\right(t,u\left(t\right)\left)\right|& \le & \left|\mathfrak{f}\right(t,u\left(t\right))-\mathfrak{f}(t,0\left)\right|+\left|\mathfrak{f}\right(t,0\left)\right|\hfill \\ & \le & \mathfrak{L}\left|u\right(t\left)\right|+\mathbb{M}\le \mathfrak{L}\parallel u\parallel +\mathbb{M}\le \mathfrak{L}r+\mathbb{M}.\hfill \end{array}$

Then, for any $u\in B_r$, we obtain

$\begin{array}{ccc}\hfill \left|\right(\mathcal{A}u\left)\right(t\left)\right|& \le & \underset{t\in [a,b]}{sup}\{\frac{{(\psi \left(t\right)-\psi \left(a\right))}^{\frac{{\theta }_k}{k}-1}}{|\Omega |{\Gamma }_k\left({\theta }_k\right)}\left[\right|\lambda |{}^k{\mathfrak{I}}^{\alpha -p;\psi }\left|\mathfrak{f}(\eta ,u\left(\eta \right))\right|+\left|\mu \right|{}^k{\mathfrak{I}}^{v+\alpha ,\psi }\left|\mathfrak{f}(\sigma ,u\left(\sigma \right))\right|\hfill \\ & & +{}^k{\mathfrak{I}}^{\alpha ;\psi }\left|\mathfrak{f}(b,u\left(b\right))\right|]+{}^k{\mathfrak{I}}^{\alpha ;\psi }\left|\mathfrak{f}(t,u\left(t\right))\right|\}\hfill \\ & \le & {}^k{\mathfrak{I}}^{\alpha ;\psi }\left(\right|\mathfrak{f}(t,u\left(t\right))-\mathfrak{f}(t,0)|+|\mathfrak{f}(t,0)\left|\right)\hfill \\ & & +\frac{{(\psi \left(b\right)-\psi \left(a\right))}^{\frac{{\theta }_k}{k}-1}}{|\Omega |{\Gamma }_k\left({\theta }_k\right)}\left(\right|\lambda |{}^k{\mathfrak{I}}^{\alpha -p;\psi }\left(\right|\mathfrak{f}(\eta ,u\left(\eta \right))-\mathfrak{f}(\eta ,0)|+|\mathfrak{f}(\eta ,0)\left|\right)\hfill \\ & & +\left|\mu \right|{}^k{\mathfrak{I}}^{v+\alpha ,\psi }\left(\right|\mathfrak{f}(\sigma ,u\left(\sigma \right))-\mathfrak{f}(\sigma ,0)|+|\mathfrak{f}(\sigma ,0)\left|\right)\hfill \\ & & {+}^k{\mathfrak{I}}^{\alpha ;\psi }\left(\right|\mathfrak{f}(b,u\left(b\right))-\mathfrak{f}(b,0)|+|\mathfrak{f}(b,0)\left|\right))\hfill \\ & \le & \{\frac{{(\psi \left(b\right)-\psi \left(a\right))}^{\frac{\alpha }{k}}}{{\Gamma }_k(\alpha +k)}+\frac{{(\psi \left(b\right)-\psi \left(a\right))}^{\frac{{\theta }_k}{k}-1}}{|\Omega |{\Gamma }_k\left({\theta }_k\right)}[\left|\lambda \right|\frac{{(\psi \left(\eta \right)-\psi \left(a\right))}^{\frac{\alpha -p}{k}}}{{\Gamma }_k(\alpha -p+k)}\hfill \\ & & +\left|\mu \right|\frac{{(\psi \left(\sigma \right)-\psi \left(a\right))}^{\frac{\alpha +v}{k}}}{{\Gamma }_k(\alpha +v+k)}+\frac{{(\psi \left(b\right)-\psi \left(a\right))}^{\frac{\alpha }{k}}}{{\Gamma }_k(\alpha +k)}\left]\right\}(\mathfrak{L}\parallel u\parallel +\mathbb{M})\hfill \\ & \le & (\mathfrak{L}r+\mathbb{M})\mathbb{G}\le r.\hfill \end{array}$

Hence, $\parallel \mathcal{A}u\parallel \le r,$ which means that $\mathcal{A}{B}_r\subset B_r$ as $u\in B_r$ is an arbitrary element.

In the second step, it will be established that $\mathcal{A}$ is a contraction. For $u,y\in C\left(\right[a,b],\mathbb{R})$ and $t\in [a,b]$, we obtain

$\begin{array}{ccc}& & \left|\right(\mathcal{A}u\left)\right(t)-(\mathcal{A}y\left)\right(t\left)\right|\hfill \\ & \le & {}^k{\mathfrak{I}}^{\alpha ;\psi }|\mathfrak{f}(t,u\left(t\right))-\mathfrak{f}(t,y\left(t\right))|\hfill \\ & & +\frac{{(\psi \left(b\right)-\psi \left(a\right))}^{\frac{{\theta }_k}{k}-1}}{|\Omega |{\Gamma }_k\left({\theta }_k\right)}{\left(\right|\lambda |}^k{\mathfrak{I}}^{\alpha -p;\psi }|\mathfrak{f}(\eta ,u\left(\eta \right))-\mathfrak{f}(\eta ,y\left(\eta \right)|\hfill \\ & & +\left|\mu \right|{}^k{\mathfrak{I}}^{\alpha +v;\psi }|\mathfrak{f}(\sigma ,u\left(\sigma \right))-\mathfrak{f}(\sigma ,y\left(\sigma \right)|\hfill \\ & & {+}^k{\mathfrak{I}}^{\alpha ;\psi }\left(\right|\mathfrak{f}(b,u\left(b\right))-\mathfrak{f}(b,y\left(b\right))\left|\right)\hfill \\ & \le & \{\frac{{(\psi \left(b\right)-\psi \left(a\right))}^{\frac{\alpha }{k}}}{{\Gamma }_k(\alpha +k)}+\frac{{(\psi \left(b\right)-\psi \left(a\right))}^{\frac{{\theta }_k}{k}-1}}{|\Omega |{\Gamma }_k\left({\theta }_k\right)}[\left|\lambda \right|\frac{{(\psi \left(\eta \right)-\psi \left(a\right))}^{\frac{\alpha -p}{k}}}{{\Gamma }_k(\alpha -p+k)}\hfill \\ & & +\left|\mu \right|\frac{{(\psi \left(\sigma \right)-\psi \left(a\right))}^{\frac{\alpha +v}{k}}}{{\Gamma }_k(\alpha +v+k)}+\frac{{(\psi \left(b\right)-\psi \left(a\right))}^{\frac{\alpha }{k}}}{{\Gamma }_k(\alpha +k)}\left]\right\}\mathfrak{L}\parallel u-y\parallel \hfill \\ & =& \mathfrak{L}\mathbb{G}\parallel u-y\parallel ,\hfill \end{array}$

which, on taking the norm for $t\in [a,b]$, yields $\parallel \mathcal{A}u-\mathcal{A}y\parallel \le \mathfrak{L}\mathbb{G}\parallel u-y\parallel .$ Since $\mathfrak{L}\mathbb{G}<1,$ therefore $\mathcal{A}$ is a contraction. As the hypothesis of the Banach’s fixed point theorem is verified, we conclude that the operator $\mathcal{A}$ has a unique fixed point, which is indeed a unique solution of the problem (1). This finishes the proof. □

4.2. Existence Results

Here we present two existence results for the problem (1) by applying Krasnosel’skiĭ’s fixed point theorem [28] and nonlinear alternative of Leray–Schauder type [29].

Theorem 2.

Suppose that $\left(H_1\right)$ and th following condition hold:

• $\left(H_2\right)$

  $\forall (t,u)\in [a,b]\times \mathbb{R}$, there exists $\theta \in C([a,b],{\mathbb{R}}^{+})$ such that $\left|\mathfrak{f}\right(t,u\left)\right|\le \theta \left(t\right).$

Then, if ${\mathbb{G}}_1\mathfrak{L}<1$, where

(21) $\begin{array}{ccc}\hfill {\mathbb{G}}_1:& =& \frac{{(\psi \left(b\right)-\psi \left(a\right))}^{\frac{{\theta }_k}{k}-1}}{|\Omega |{\Gamma }_k\left({\theta }_k\right)}\left[\right|\lambda |\frac{{(\psi \left(\eta \right)-\psi \left(a\right))}^{\frac{\alpha -p}{k}}}{{\Gamma }_k(\alpha -p+k)}+\left|\mu \right|\frac{{(\psi \left(\sigma \right)-\psi \left(a\right))}^{\frac{\alpha +v}{k}}}{{\Gamma }_k(\alpha +v+k)}\hfill \\ & & +\frac{{(\psi \left(b\right)-\psi \left(a\right))}^{\frac{\alpha }{k}}}{{\Gamma }_k(\alpha +k)}],\hfill \end{array}$

the problem (1) has at least one solution on $[a,b].$

Proof. 

Set ${sup}_{t\in [a,b]}\theta \left(t\right)=\parallel \theta \parallel$ and consider the ball $B_{\rho }=\{u\in C([a,b],\mathbb{R}):\parallel u\parallel \le \rho \},$ with $\rho \ge \parallel \theta \parallel \mathbb{G}.$ Introduce the operators ${\mathcal{A}}_1$ and ${\mathcal{A}}_2$ on $B_{\rho }$ to $\mathbb{R}$ as

$\begin{array}{c}\hfill {\mathcal{A}}_{1}u\left(t\right)={}^k{\mathfrak{I}}^{\alpha ;\psi }\mathfrak{f}(t,u\left(t\right)),t\in [a,b],\end{array}$

$\begin{array}{ccc}\hfill \left({\mathcal{A}}_{2}u\right)\left(t\right)& =& \frac{{(\psi \left(t\right)-\psi \left(a\right))}^{\frac{{\theta }_k}{k}-1}}{\Omega {\Gamma }_k\left({\theta }_k\right)}[\lambda {}^k{\mathfrak{I}}^{\alpha -p;\psi }\mathfrak{f}(\eta ,u\left(\eta \right))+\mu {}^k{\mathfrak{I}}^{v+\alpha ,\psi }\mathfrak{f}(\sigma ,u\left(\sigma \right))\hfill \\ & & -{}^k{\mathfrak{I}}^{\alpha ;\psi }\mathfrak{f}(b,u\left(b\right))],t\in [a,b].\hfill \end{array}$

For any $u,y\in B_{\rho }$, we have

$\begin{array}{ccc}& & |\left({\mathcal{A}}_{1}u\right)\left(t\right)+\left({\mathcal{A}}_{2}y\right)\left(t\right)|\hfill \\ & \le & \underset{t\in [a,b]}{sup}\{\frac{{(\psi \left(t\right)-\psi \left(a\right))}^{\frac{{\theta }_k}{k}-1}}{|\Omega |{\Gamma }_k\left({\theta }_k\right)}\left[\right|\lambda |{}^k{\mathfrak{I}}^{\alpha -p;\psi }\left|\mathfrak{f}(\eta ,u\left(\eta \right))\right|+\left|\mu \right|{}^k{\mathfrak{I}}^{v+\alpha ,\psi }\left|\mathfrak{f}(\sigma ,y\left(\sigma \right))\right|\hfill \\ & & +{}^k{\mathfrak{I}}^{\alpha ;\psi }\left|\mathfrak{f}(b,y\left(b\right))\right|]+{}^k{\mathfrak{I}}^{\alpha ;\psi }\left|\mathfrak{f}(t,u\left(t\right))\right|\}\hfill \\ & \le & \{\frac{{(\psi \left(b\right)-\psi \left(a\right))}^{\frac{\alpha }{k}}}{{\Gamma }_k(\alpha +k)}+\frac{{(\psi \left(b\right)-\psi \left(a\right))}^{\frac{{\theta }_k}{k}-1}}{|\Omega |{\Gamma }_k\left({\theta }_k\right)}[\left|\lambda \right|\frac{{(\psi \left(\eta \right)-\psi \left(a\right))}^{\frac{\alpha -p}{k}}}{{\Gamma }_k(\alpha -p+k)}\hfill \\ & & +\left|\mu \right|\frac{{(\psi \left(\sigma \right)-\psi \left(a\right))}^{\frac{\alpha +v}{k}}}{{\Gamma }_k(\alpha +v+k)}+\frac{{(\psi \left(b\right)-\psi \left(a\right))}^{\frac{\alpha }{k}}}{{\Gamma }_k(\alpha +k)}\left]\right\}\parallel \theta \parallel \hfill \\ & =& \mathbb{G}\parallel \theta \parallel \le \rho .\hfill \end{array}$

Therefore $\parallel \left({\mathcal{A}}_{1}u\right)+\left({\mathcal{A}}_{2}y\right)\parallel \le \rho ,$ which shows that ${\mathcal{A}}_{1}u+{\mathcal{A}}_{2}y\in B_{\rho }.$ It is easy to show, by using (21), that the operator ${\mathcal{A}}_2$ is a contraction.

Note that continuity of $\mathfrak{f}$ implies the continuity of the operator ${\mathcal{A}}_1.$ Furthermore, the operator ${\mathcal{A}}_1$ is uniformly bounded on $B_{\rho }$ as

$\parallel {\mathcal{A}}_{1}u\parallel \le \frac{{(\psi \left(b\right)-\psi \left(a\right))}^{\frac{\alpha }{k}}}{{\Gamma }_k(\alpha +k)}\parallel \theta \parallel .$

Now we establish that the operator ${\mathcal{A}}_1$ is compact. For $t_1<t_2$ ($t_1,t_2\in [a,b]$), we have

$\begin{array}{ccc}& & |\left({\mathcal{A}}_{1}u\right)\left(t_2\right)-\left({\mathcal{A}}_{1}u\right)\left(t_1\right)|\hfill \\ & \le & \frac{1}{{\Gamma }_k\left(\alpha \right)}|{\int }_a^{t_1}{\psi }^{\prime }\left(s\right)[{(\psi \left(t_2\right)-\psi \left(s\right))}^{\frac{\alpha }{k}-1}-{(\psi \left(t_1\right)-\psi \left(s\right))}^{\frac{\alpha }{k}-1}]\mathfrak{f}(s,u\left(s\right))ds\hfill \\ & & +{\int }_{t_1}^{t_2}{\psi }^{\prime }\left(s\right){(\psi \left(t_2\right)-\psi \left(s\right))}^{\frac{\alpha }{k}-1}\mathfrak{f}(s,u\left(s\right))ds|\hfill \\ & \le & \frac{\parallel \theta \parallel }{{\Gamma }_k(\alpha +k)}[2{(\psi \left(t_2\right)-\psi \left(t_1\right))}^{\frac{\alpha }{k}}+|{(\psi \left(t_2\right)-\psi \left(a\right))}^{\frac{\alpha }{k}}-{(\psi \left(t_1\right)-\psi \left(a\right))}^{\frac{\alpha }{k}}\left|\right]\hfill \\ & ⟶& 0\mathrm{as}{t}_2-t_1\to 0,\hfill \end{array}$

independently of $u.$ Hence, the operator ${\mathcal{A}}_1$ is equicontinuous, and consequently, it follows by the Arzelá–Ascoli theorem that the operator ${\mathcal{A}}_1$ is completely continuous. Thus, the conclusion of Krasnosel’skiĭ’s fixed point theorem applies, and the problem (1) has at least one solution on $[a,b].$ The proof is completed. □

Theorem 3.

Suppose that the following conditions hold:

• ($H_3$)

  ∃ a continuous, nondecreasing function $\chi :[0,+\infty )\to (0,+\infty )$ and a function $\gamma \in C([a,b],{\mathbb{R}}^{+})$ such that, $\forall (t,u)\in [a,b]\times \mathbb{R},\left|\mathfrak{f}\right(t,u\left)\right|\le \gamma \left(t\right)\chi \left(\right|u\left|\right)$;

• ($H_4$)

  ∃ a constant $\mathcal{K}>0$ such that

  $\frac{\mathcal{K}}{\chi \left(\mathcal{K}\right)\parallel \gamma \parallel \mathbb{G}}>1.$

Then there exists at least one solution for the problem (1) on $[a,b]$.

Proof. 

Consider the operator $\mathcal{A}$ is defined by (17). We will split the proof into several steps. It will be shown in the first step that the operator $\mathcal{A}$ maps bounded sets into bounded set in $C\left(\right[a,b],\mathbb{R})$. For $r>0$, let $B_r=\{u\in C([a,b],\mathbb{R}):\parallel u\parallel \le r\}.$ Then, for $t\in [a,b]$, we obtain

$\begin{array}{ccc}\hfill \left|\right(\mathcal{A}u\left)\right(t\left)\right|& \le & \underset{t\in [a,b]}{sup}\{\frac{{(\psi \left(t\right)-\psi \left(a\right))}^{\frac{{\theta }_k}{k}-1}}{|\Omega |{\Gamma }_k\left({\theta }_k\right)}\left[\right|\lambda |{}^k{\mathfrak{I}}^{\alpha -p;\psi }\left|\mathfrak{f}(\eta ,u\left(\eta \right))\right|+\left|\mu \right|{}^k{\mathfrak{I}}^{v+\alpha ,\psi }\left|\mathfrak{f}(\sigma ,u\left(\sigma \right))\right|\hfill \\ & & +{}^k{\mathfrak{I}}^{\alpha ;\psi }\left|\mathfrak{f}(b,u\left(b\right))\right|]+{}^k{\mathfrak{I}}^{\alpha ;\psi }\left|\mathfrak{f}(t,u\left(t\right))\right|\}\hfill \\ & \le & \{\frac{{(\psi \left(b\right)-\psi \left(a\right))}^{\frac{\alpha }{k}}}{{\Gamma }_k(\alpha +k)}+\frac{{(\psi \left(b\right)-\psi \left(a\right))}^{\frac{{\theta }_k}{k}-1}}{|\Omega |{\Gamma }_k\left({\theta }_k\right)}[\left|\lambda \right|\frac{{(\psi \left(\eta \right)-\psi \left(a\right))}^{\frac{\alpha -p}{k}}}{{\Gamma }_k(\alpha -p+k)}\hfill \\ & & +\left|\mu \right|\frac{{(\psi \left(\sigma \right)-\psi \left(a\right))}^{\frac{\alpha +v}{k}}}{{\Gamma }_k(\alpha +v+k)}+\frac{{(\psi \left(b\right)-\psi \left(a\right))}^{\frac{\alpha }{k}}}{{\Gamma }_k(\alpha +k)}\left]\right\}\parallel \gamma \parallel \chi (\parallel u\parallel ),\hfill \end{array}$

which implies that

$\parallel \mathcal{A}x\parallel \le \chi \left(r\right)\parallel \gamma \parallel \mathbb{G}.$

Now, we establish that $\mathcal{A}$ maps bounded sets into equicontinuous sets of $C\left(\right[a,b],\mathbb{R})$. Consider $u\in B_r$ and $t_1,t_2\in [a,b]$ such that $t_1<t_2.$ Then we obtain

$\begin{array}{ccc}& & |\left(\mathcal{A}u\right)\left(t_2\right)-\left(\mathcal{A}u\right)\left(t_1\right)|\hfill \\ & \le & \frac{1}{{\Gamma }_k\left(\alpha \right)}|{\int }_a^{t_1}{\psi }^{\prime }\left(s\right)[{(\psi \left(t_2\right)-\psi \left(s\right))}^{\frac{\alpha }{k}-1}-{(\psi \left(t_1\right)-\psi \left(s\right))}^{\frac{\alpha }{k}-1}]\mathfrak{f}(s,u\left(s\right))ds\hfill \\ & & +{\int }_{t_1}^{t_2}{\psi }^{\prime }\left(s\right){(\psi \left(t_2\right)-\psi \left(s\right))}^{\frac{\alpha }{k}-1}\mathfrak{f}(s,u\left(s\right))ds|\hfill \\ & & +\frac{{(\psi \left(t_2\right)-\psi \left(a\right))}^{\frac{{\theta }_k}{k}-1}-{(\psi \left(t_1\right)-\psi \left(a\right))}^{\frac{{\theta }_k}{k}-1}}{|\Omega |{\Gamma }_k\left({\theta }_k\right)}\left[\right|\lambda |{}^k{\mathfrak{I}}^{\alpha -p;\psi }\left|\mathfrak{f}(\eta ,u\left(\eta \right))\right|\hfill \\ & & +\left|\mu \right|{}^k{\mathfrak{I}}^{v+\alpha ,\psi }\left|\mathfrak{f}(\sigma ,u\left(\sigma \right))\right|+{}^k{\mathfrak{I}}^{\alpha ;\psi }\left|\mathfrak{f}(b,u\left(b\right))\right|]\hfill \\ & \le & \frac{\parallel \gamma \parallel \chi \left(r\right)}{{\Gamma }_k(\alpha +k)}[2{(\psi \left(t_2\right)-\psi \left(t_1\right))}^{\frac{\alpha }{k}}+|{(\psi \left(t_2\right)-\psi \left(a\right))}^{\frac{\alpha }{k}}-{(\psi \left(t_1\right)-\psi \left(a\right))}^{\frac{\alpha }{k}}\left|\right],\hfill \\ & & +\frac{{(\psi \left(t_2\right)-\psi \left(a\right))}^{\frac{{\theta }_k}{k}-1}-{(\psi \left(t_1\right)-\psi \left(a\right))}^{\frac{{\theta }_k}{k}-1}}{|\Omega |{\Gamma }_k\left({\theta }_k\right)}[\left|\lambda \right|\frac{{(\psi \left(\eta \right)-\psi \left(a\right))}^{\frac{\alpha -p}{k}}}{{\Gamma }_k(\alpha -p+k)}\hfill \\ & & +\left|\mu \right|\frac{{(\psi \left(\sigma \right)-\psi \left(a\right))}^{\frac{\alpha +v}{k}}}{{\Gamma }_k(\alpha +v+k)}+\frac{{(\psi \left(b\right)-\psi \left(a\right))}^{\frac{\alpha }{k}}}{{\Gamma }_k(\alpha +k)}]\parallel \gamma \parallel \chi \left(r\right)\hfill \\ & & ⟶0\mathrm{as}{t}_2-t_1\to 0,\hfill \end{array}$

independently of $u\in B_r.$ Thus, an immediate consequence of the Arzelá–Ascoli theorem implies that the operator $\mathcal{A}:C\left(\right[a,b],\mathbb{R})\to C\left(\right[a,b],\mathbb{R})$ is completely continuous.

In the final step, we establish the boundedness of the set of all solutions to the equation $u=\omega \mathcal{A}u$ for $\omega \in (0,1)$.

As in the first step, we obtain

$\left|u\right(t\left)\right|\le \chi (\parallel u\parallel )\parallel \gamma \parallel \mathbb{G},$

or

$\frac{\parallel u\parallel }{\chi (\parallel u\parallel )\parallel \gamma \parallel \mathbb{G}}\le 1.$

By $\left(H_4\right)$, we can find $\mathcal{K}$ such that $\parallel u\parallel \ne \mathcal{K}.$ Consider the set

$U=\{u\in C([a,b],\mathbb{R}):\parallel u\parallel <\mathcal{K}\}.$

Observe that the operator $\mathcal{A}:\overline{U}\to C([a,b],\mathbb{R})$ is continuous and completely continuous, where $\overline{U}$ denotes the closure of U. By the given choice of $U,$ we cannot find any $u\in \partial U$ ($\partial U$ denotes the boundary of U) such that $u=\omega \mathcal{A}u$ for some $\omega \in (0,1)$ in view of the assumption $\left(H_4\right)$. In consequence, we deduce by the nonlinear alternative of Leray–Schauder type the operator $\mathcal{A}$ has a fixed point $u\in \overline{U}$, which is a solution of the problem (1). This finishes the proof. □

5. The Multivalued Problem

Definition 8.

A function $u\in C^2([a,b],\mathbb{R})$ possessing the $(k, \psi )$-Hilfer fractional derivative in the sense of Definition 6 is said to be a solution of the $(k, \psi )$-Hilfer-type nonlocal fractional multi-valued boundary value problem (2) if there exists a function $f\in L^1([a,b],\mathbb{R})$ with $f\left(t\right)\in \mathfrak{F}(t,u)$ for a.e. $t\in [a,b]$ such that u satisfies the differential equation ${}^{k,H}{D}^{\alpha ,\beta ;\psi }u\left(t\right)=f\left(t\right)$ on $[a,b]$ and the boundary conditions $u\left(a\right)=0,u\left(b\right)=\lambda {}^{k,H}{D}^{p,q;\psi }u\left(\eta \right)+\mu {}^k{\mathfrak{I}}^{v,\psi }u\left(\sigma \right)$.

For each $u\in C\left(\right[a,b],\mathbb{R}),$ we define the set of selections of $\mathfrak{F}$ as

$S_{\mathfrak{F},u}:=\{f\in L^1([a,b],\mathbb{R}):f\left(t\right)\in \mathfrak{F}(t,u\left(t\right))\mathrm{on}[a,b]\}.$

Our first result for the multi-valued problem (2) is concerned with the case when the multi-valued map $\mathfrak{F}$ has convex values, and relies on the nonlinear alternative of Leray–Schauder type for multi-valued maps [29].

Theorem 4.

Suppose that:

• $\left(G_1\right)$

  $\mathfrak{F}:[a,b]\times \mathbb{R}\to {\mathcal{P}}_{cp,c}\left(\mathbb{R}\right)$ is $L^1$-Carathéodory, where ${\mathcal{P}}_{cp,c}\left(\mathbb{R}\right)=\{\mathfrak{R}\in \mathcal{P}\left(\mathbb{R}\right):\mathfrak{R} is compact and convex\};$

• $\left(G_2\right)$

  ∃ a continuous nondecreasing function $z:[0,+\infty )\to (0,+\infty )$ and a positive continuous real valued function q such that, $\forall (t,u)\in [a,b]\times \mathbb{R}$,

  ${\parallel \mathfrak{F}(t,u)\parallel }_{\mathcal{P}}:=sup\left\{\right|f|:f\in \mathfrak{F}(t,u)\}\le q\left(t\right)z(\parallel u\parallel );$

• $\left(G_3\right)$

  ∃ a constant $\mathfrak{M}>0$ such that

  $\frac{{\displaystyle \mathfrak{M}}}{{\displaystyle \parallel q\parallel z\left(\mathfrak{M}\right)\mathbb{G}}}>1,$

  where $\mathbb{G}$ is defined by (18).

Then the multi-valued problem (2) has at least one solution on $[a,b].$

Proof. 

Define an operator $\mathcal{F}:C\left(\right[a,b],\mathbb{R})⟶\mathcal{P}\left(C\right([a,b],\mathbb{R}\left)\right)$ as

$\mathcal{F}\left(u\right)=\left\{\begin{array}{cc}& h\in C\left(\right[a,b],\mathbb{R}):\hfill \\ & h\left(t\right)=\left\{\begin{array}{c}{\displaystyle \frac{{(\psi \left(t\right)-\psi \left(a\right))}^{\frac{{\theta }_k}{k}-1}}{\Omega {\Gamma }_k\left({\theta }_k\right)}\left[\lambda {}^k{\mathfrak{I}}^{\alpha -p;\psi }f\left(\eta \right)+\mu {}^k{\mathfrak{I}}^{\alpha +v;\psi }f\left(\sigma \right)-{}^k{\mathfrak{I}}^{\alpha ;\psi }f\left(b\right)\right]}\hfill \\ {\displaystyle +{}^k{\mathfrak{I}}^{\alpha ;\psi }f\left(t\right),}\hfill \end{array}\right.\hfill \end{array}\right\}$

for $t\in [a,b]$ and $f\in S_{\mathfrak{F},x}.$ Observe that the fixed points of the operator $\mathcal{F}$ are solutions to the multi-valued problem (2).

We will split the proof in different steps.

Step 1. $\mathcal{F}\left(u\right)$ is convex for each $u\in C\left(\right[a,b],\mathbb{R}.$

Since $S_{F,u}$ is convex, this step is obvious.

Step 2. Bounded sets are mapped into bounded sets in $C\left(\right[a,b],\mathbb{R})$ by $\mathcal{F}.$

Let $B_r=\{u\in C([a,b],\mathbb{R}):\parallel u\parallel \le r\},r>0.$ Then, for each $h\in \mathcal{F}\left(u\right),u\in B_r$, there exists $f\in S_{\mathfrak{F},x}$ such that

$\begin{array}{c}\hfill h\left(t\right)=\frac{{(\psi \left(t\right)-\psi \left(a\right))}^{\frac{{\theta }_k}{k}-1}}{\Omega {\Gamma }_k\left({\theta }_k\right)}\left[\lambda {}^k{\mathfrak{I}}^{\alpha -p;\psi }f\left(\eta \right)+\mu {}^k{\mathfrak{I}}^{\alpha +v;\psi }f\left(\sigma \right)-{}^k{\mathfrak{I}}^{\alpha ;\psi }f\left(b\right)\right]+{}^k{\mathfrak{I}}^{\alpha ;\psi }f\left(t\right).\end{array}$

Then, for $t\in [a,b],$ we have

$\begin{array}{ccc}\hfill \left|h\right(t\left)\right|& \le & \underset{t\in [a,b]}{sup}\{\frac{{(\psi \left(t\right)-\psi \left(a\right))}^{\frac{{\theta }_k}{k}-1}}{|\Omega |{\Gamma }_k\left({\theta }_k\right)}\left[\left|\lambda \right|{}^k{\mathfrak{I}}^{\alpha -p;\psi }\left|f\left(\eta \right)\right|+\left|\mu \right|{}^k{\mathfrak{I}}^{\alpha +v;\psi }\left|f\left(\sigma \right)\right|+{}^k{\mathfrak{I}}^{\alpha ;\psi }\left|f\left(b\right)\right|\right]\hfill \\ & & +{}^k{\mathfrak{I}}^{\alpha ;\psi }\left|f\left(t\right)\right|\}\hfill \\ & \le & \{\frac{{(\psi \left(b\right)-\psi \left(a\right))}^{\frac{\alpha }{k}}}{{\Gamma }_k(\alpha +k)}+\frac{{(\psi \left(b\right)-\psi \left(a\right))}^{\frac{{\theta }_k}{k}-1}}{|\Omega |{\Gamma }_k\left({\theta }_k\right)}[\left|\lambda \right|\frac{{(\psi \left(\eta \right)-\psi \left(a\right))}^{\frac{\alpha -p}{k}}}{{\Gamma }_k(\alpha -p+k)}\hfill \\ & & +\left|\mu \right|\frac{{(\psi \left(\sigma \right)-\psi \left(a\right))}^{\frac{\alpha +v}{k}}}{{\Gamma }_k(\alpha +v+k)}+\frac{{(\psi \left(b\right)-\psi \left(a\right))}^{\frac{\alpha }{k}}}{{\Gamma }_k(\alpha +k)}\left]\right\}\parallel q\parallel z(\parallel u\parallel ),\hfill \end{array}$

which implies that $\parallel h\parallel \le z\left(r\right)\parallel q\parallel \mathbb{G}.$

Step 3. Bounded sets are mapped into equicontinuous sets of $C\left(\right[a,b],\mathbb{R})$ by $\mathcal{F}$.

Consider $t_1,t_2\in [a,b],t_1<t_2$ and $u\in B_r.$ Then, for each $h\in \mathcal{F}\left(u\right),$ we have

$\begin{array}{ccc}& & |h\left(t_2\right)-h\left(t_1\right)|\hfill \\ & \le & \frac{1}{{\Gamma }_k\left(\alpha \right)}|{\int }_a^{t_1}{\psi }^{\prime }\left(s\right)[{(\psi \left(t_2\right)-\psi \left(s\right))}^{\frac{\alpha }{k}-1}-{(\psi \left(t_1\right)-\psi \left(s\right))}^{\frac{\alpha }{k}-1}]f\left(s\right)ds\hfill \\ & & +{\int }_{t_1}^{t_2}{\psi }^{\prime }\left(s\right){(\psi \left(t_2\right)-\psi \left(s\right))}^{\frac{\alpha }{k}-1}f\left(s\right)ds|\hfill \\ & & +\frac{{(\psi \left(t_2\right)-\psi \left(a\right))}^{\frac{{\theta }_k}{k}-1}-{(\psi \left(t_1\right)-\psi \left(a\right))}^{\frac{{\theta }_k}{k}-1}}{|\Omega |{\Gamma }_k\left({\theta }_k\right)}\left[\right|\lambda |{}^k{\mathfrak{I}}^{\alpha -p;\psi }\left|f\left(\eta \right)\right|\hfill \\ & & +\left|\mu \right|{}^k{\mathfrak{I}}^{\alpha +v;\psi }\left|f\left(\sigma \right)\right|+{}^k{\mathfrak{I}}^{\alpha ;\psi }\left|f\left(b\right)\right|]\hfill \\ & \le & \frac{\parallel q\parallel z\left(r\right)}{{\Gamma }_k(\alpha +k)}[2{(\psi \left(t_2\right)-\psi \left(t_1\right))}^{\frac{\alpha }{k}}+|{(\psi \left(t_2\right)-\psi \left(a\right))}^{\frac{\alpha }{k}}-{(\psi \left(t_1\right)-\psi \left(a\right))}^{\frac{\alpha }{k}}\left|\right],\hfill \\ & & +\frac{{(\psi \left(t_2\right)-\psi \left(a\right))}^{\frac{{\theta }_k}{k}-1}-{(\psi \left(t_1\right)-\psi \left(a\right))}^{\frac{{\theta }_k}{k}-1}}{|\Omega |{\Gamma }_k\left({\theta }_k\right)}[\left|\lambda \right|\frac{{(\psi \left(\eta \right)-\psi \left(a\right))}^{\frac{\alpha -p}{k}}}{{\Gamma }_k(\alpha -p+k)}\hfill \\ & & +\left|\mu \right|\frac{{(\psi \left(\sigma \right)-\psi \left(a\right))}^{\frac{\alpha +v}{k}}}{{\Gamma }_k(\alpha +v+k)}+\frac{{(\psi \left(b\right)-\psi \left(a\right))}^{\frac{\alpha }{k}}}{{\Gamma }_k(\alpha +k)}]\parallel q\parallel z\left(r\right)\hfill \\ & & ⟶0\mathrm{as}{t}_2-t_1\to 0,\hfill \end{array}$

independently of $u\in B_r$. Hence it follows by the Arzelá–Ascoli theorem that $\mathcal{F}:C\left(\right[a,b],\mathbb{R})\to \mathcal{P}\left(C\right([a,b],\mathbb{R}\left)\right)$ is completely continuous.

To prove that $\mathcal{F}$ is upper semi-continuous multivalued mapping, it is enough to prove that the $\mathcal{F}$ has a closed graph, by Proposition 1.2 of [30].

Step 4. $\mathcal{F}$ has a closed graph.

Assume that $u_n\to u_{*},h_n\in \mathcal{F}\left(u_n\right)$ and $h_n\to h_{*}.$ Then we must to show that $h_{*}\in \mathcal{F}\left(u_{*}\right).$ Since $h_n\in \mathcal{F}\left(u_n\right),$ there exists $v_n\in S_{\mathfrak{F},u_n}$ such that for each $t\in [a,b],$

$\begin{array}{c}\hfill h_n\left(t\right)=\frac{{(\psi \left(t\right)-\psi \left(a\right))}^{\frac{{\theta }_k}{k}-1}}{\Omega {\Gamma }_k\left({\theta }_k\right)}\left[\lambda {}^k{\mathfrak{I}}^{\alpha -p;\psi }{f}_n\left(\eta \right)+\mu {}^k{\mathfrak{I}}^{\alpha +v;\psi }{f}_n\left(\sigma \right)-{}^k{\mathfrak{I}}^{\alpha ;\psi }{f}_n\left(b\right)\right]+{}^k{\mathfrak{I}}^{\alpha ;\psi }{f}_n\left(t\right).\end{array}$

We must show that there exists $v_{*}\in S_{\mathfrak{F},u_{*}}$ such that for each $t\in [a,b],$

$\begin{array}{c}\hfill h_{*}\left(t\right)=\frac{{(\psi \left(t\right)-\psi \left(a\right))}^{\frac{{\theta }_k}{k}-1}}{\Omega {\Gamma }_k\left({\theta }_k\right)}\left[\lambda {}^k{\mathfrak{I}}^{\alpha -p;\psi }{f}_{*}\left(\eta \right)+\mu {}^k{\mathfrak{I}}^{\alpha +v;\psi }{f}_{*}\left(\sigma \right)-{}^k{\mathfrak{I}}^{\alpha ;\psi }{f}_{*}\left(b\right)\right]+{}^k{\mathfrak{I}}^{\alpha ;\psi }{f}_{*}\left(t\right).\end{array}$

Consider the linear operator $\Theta :L^1([a,b],\mathbb{R})\to C([a,b],\mathbb{R})$ given by

$\begin{array}{c}\hfill v\mapsto \Theta \left(v\right)\left(t\right)=\frac{{(\psi \left(t\right)-\psi \left(a\right))}^{\frac{{\theta }_k}{k}-1}}{\Omega {\Gamma }_k\left({\theta }_k\right)}\left[\lambda {}^k{\mathfrak{I}}^{\alpha -p;\psi }f\left(\eta \right)+\mu {}^k{\mathfrak{I}}^{\alpha +v;\psi }f\left(\sigma \right)-{}^k{\mathfrak{I}}^{\alpha ;\psi }f\left(b\right)\right]+{}^k{\mathfrak{I}}^{\alpha ;\psi }f\left(t\right).\end{array}$

Observe that $\parallel h_n\left(t\right)-h_{*}\left(t\right)\parallel \to 0,$ as $n\to +\infty .$ From a result due to Lazota–Opial [31], we deduce that $\Theta \circ S_{\mathfrak{F}}$ is a closed graph operator, and moreover we have $h_n\left(t\right)\in \Theta \left(S_{\mathfrak{F},u_n}\right).$ We have, since $u_n\to u_{*},$

$\begin{array}{c}\hfill h_{*}\left(t\right)=\frac{{(\psi \left(t\right)-\psi \left(a\right))}^{\frac{{\theta }_k}{k}-1}}{\Omega {\Gamma }_k\left({\theta }_k\right)}\left[\lambda {}^k{\mathfrak{I}}^{\alpha -p;\psi }{f}_{*}\left(\eta \right)+\mu {}^k{\mathfrak{I}}^{\alpha +v;\psi }{f}_{*}\left(\sigma \right)-{}^k{\mathfrak{I}}^{\alpha ;\psi }{f}_{*}\left(b\right)\right]+{}^k{\mathfrak{I}}^{\alpha ;\psi }{f}_{*}\left(t\right),\end{array}$

for some $v_{*}\in S_{\mathfrak{F},u_{*}}$.

Step 5. An open set $\mathcal{U}\subseteq C\left(\right[a,b],\mathbb{R})$ exists with $u\notin \nu \mathcal{F}\left(u\right)$ for any $\nu \in (0,1)$ and all $u\in \partial \mathcal{U}.$

Consider $\nu \in (0,1)$ and $u\in \nu \mathcal{F}\left(u\right).$ Then there exists $f\in L^1([a,b],\mathbb{R})$ with $f\in S_{\mathfrak{F},u}$ such that, for $t\in [a,b]$, we have

$\begin{array}{c}\hfill u\left(t\right)=\nu \frac{{(\psi \left(t\right)-\psi \left(a\right))}^{\frac{{\theta }_k}{k}-1}}{\Omega {\Gamma }_k\left({\theta }_k\right)}\left[\lambda {}^k{\mathfrak{I}}^{\alpha -p;\psi }f\left(\eta \right)+\mu {}^k{\mathfrak{I}}^{\alpha +v;\psi }f\left(\sigma \right)-{}^k{\mathfrak{I}}^{\alpha ;\psi }f\left(b\right)\right]+\nu {}^k{\mathfrak{I}}^{\alpha ;\psi }f\left(t\right).\end{array}$

As in second step, we obtain

$\begin{array}{c}\hfill \left|u\right(t\left)\right|\le \parallel q\parallel z(\parallel u\parallel )\mathbb{G}.\end{array}$

Consequently

$\begin{array}{c}\hfill \parallel u\parallel \le \parallel q\parallel z(\parallel u\parallel )\mathbb{G},\end{array}$

or

$\frac{{\displaystyle \parallel u\parallel }}{{\displaystyle \parallel q\parallel z(\parallel u\parallel )\mathbb{G}}}\le 1.$

By $\left(H_3\right)$, we can find $\mathfrak{M}$ satisfying $\parallel u\parallel \ne \mathfrak{M}$. Consider the set

$\mathcal{U}=\{u\in C([a,b],\mathbb{R}):\parallel u\parallel <\mathfrak{M}\}.$

From the preceding arguments, $\mathcal{F}:\overline{U}\to \mathcal{P}\left(C([a,b],\mathbb{R})\right)$ is a compact and upper semi-continuous multivalued map with convex closed values. By definition $\mathcal{U},$ there does not exist any $u\in \partial \mathcal{U}$ such that $u\in \nu \mathcal{F}\left(u\right)$ for some $\nu \in (0,1)$. Hence, it follows by the nonlinear alternative of Leray–Schauder type for multi-valued maps [29] that $\mathcal{F}$ has a fixed point $u\in \overline{\mathcal{U}}$, which is indeed a solution of the multi-valued problem (2). The proof is complete. □

Now we consider the case when F is not necessarily convex valued and show that there exists a solution to the problem (2) with the aid of a fixed point theorem for multivalued contractive maps due to Covitz and Nadler [32].

Theorem 5.

Suppose that:

• $\left(A_1\right)$

  $\mathfrak{F}:[a,b]\times \mathbb{R}\to {\mathcal{P}}_{cp}\left(\mathbb{R}\right)$ is such that $\mathfrak{f}(·,u):[a,b]\to {\mathcal{P}}_{cp}\left(\mathbb{R}\right)$ is measurable for each $u\in \mathbb{R}$, where ${\mathcal{P}}_{cp}\left(\mathbb{R}\right)=\{\mathfrak{K}\in \mathcal{P}\left(\mathbb{R}\right):\mathfrak{K}is compact\};$

• $\left(A_2\right)$

  ∃ a function $m\in C([a,b],{\mathbb{R}}^{+})$ such that

  $H_d(\mathfrak{F}(t,u),\mathfrak{F}(t,\overline{u}))\le m\left(t\right)|u-\overline{u}|,$

  with $d(0,\mathfrak{F}(t,0\left)\right)\le m\left(t\right)$ for almost all $t\in [a,b]$ and $u,\overline{u}\in \mathbb{R}.$

Then the problem (2) has at least one solution on $[a,b],$ provided that

(22) $\delta :=\mathbb{G}\parallel m\parallel <1,$

where $\mathbb{G}$ is given by (18).

Proof. 

Observe that the set $S_{\mathfrak{F},u}$ is nonempty for each $u\in C\left(\right[a,b],\mathbb{R}),$ by the assumption $\left(A_1\right).$ Thus, by Theorem III.6 [33], $\mathfrak{F}$ has a measurable selection. Now for each $u\in C\left(\right[a,b],\mathbb{R}),$ it will be shown that $\mathcal{F}\left(u\right)\in {\mathcal{P}}_{cl}\left(C([a,b],\mathbb{R})\right),$ where ${\mathcal{P}}_{cl}\left(\mathbb{R}\right)=\{\mathfrak{K}\in \mathcal{P}\left(\mathbb{R}\right):\mathfrak{K}\mathrm{is} \mathrm{closed}\}.$ Assume that ${\left\{u_n\right\}}_{n\ge 0}\in \mathcal{F}\left(u\right)$ with $u_n\to u(n\to +\infty )$ in $C\left(\right[a,b],\mathbb{R}).$ Then we have $u\in C\left(\right[a,b],\mathbb{R})$ and $v_n\in S_{\mathfrak{F},u_n}$ such that, for each $t\in [a,b]$,

$\begin{array}{c}\hfill u_n\left(t\right)=\frac{{(\psi \left(t\right)-\psi \left(a\right))}^{\frac{{\theta }_k}{k}-1}}{\Omega {\Gamma }_k\left({\theta }_k\right)}\left[\lambda {}^k{\mathfrak{I}}^{\alpha -p;\psi }{f}_n\left(\eta \right)+\mu {}^k{\mathfrak{I}}^{\alpha +v;\psi }{f}_n\left(\sigma \right)-{}^k{\mathfrak{I}}^{\alpha ;\psi }{f}_n\left(b\right)\right]+{}^k{\mathfrak{I}}^{\alpha ;\psi }{f}_n\left(t\right).\end{array}$

Since $\mathfrak{F}$ has compact values, one can pass onto a subsequence (if necessary) to get $v_n$ converges to v in $L^1([a,b],\mathbb{R}).$ In consequence, $v\in S_{\mathfrak{F},x}$ and we have

$\begin{array}{ccc}\hfill u_n\left(t\right)\to u\left(t\right)& =& \frac{{(\psi \left(t\right)-\psi \left(a\right))}^{\frac{{\theta }_k}{k}-1}}{\Omega {\Gamma }_k\left({\theta }_k\right)}\left[\lambda {}^k{\mathfrak{I}}^{\alpha -p;\psi }f\left(\eta \right)+\mu {}^k{\mathfrak{I}}^{\alpha +v;\psi }f\left(\sigma \right)-{}^k{\mathfrak{I}}^{\alpha ;\psi }f\left(b\right)\right]\hfill \\ & & +{}^k{\mathfrak{I}}^{\alpha ;\psi }f\left(t\right),\mathrm{for}\mathrm{each}t\in [a,b].\hfill \end{array}$

Therefore, $u\in \mathcal{F}\left(u\right).$

Finally we show that

$H_d(\mathcal{F}\left(u\right),\mathcal{F}\left(\overline{u}\right))\le \delta \parallel u-\overline{u}\parallel ,\delta <1,\mathrm{for} \mathrm{each}u,\overline{u}\in C^2([a,b],\mathbb{R}).$

Assume that $u,\overline{u}\in C^2([a,b],\mathbb{R})$ and $h_1\in \mathcal{F}\left(u\right)$. Then we can find $v_1\left(t\right)\in \mathfrak{F}(t,u\left(t\right))$ satisfying

$\begin{array}{c}\hfill h_1\left(t\right)=\frac{{(\psi \left(t\right)-\psi \left(a\right))}^{\frac{{\theta }_k}{k}-1}}{\Omega {\Gamma }_k\left({\theta }_k\right)}\left[\lambda {}^k{\mathfrak{I}}^{\alpha -p;\psi }{f}_1\left(\eta \right)+\mu {}^k{\mathfrak{I}}^{\alpha +v;\psi }{f}_1\left(\sigma \right)-{}^k{\mathfrak{I}}^{\alpha ;\psi }{f}_1\left(b\right)\right]+{}^k{\mathfrak{I}}^{\alpha ;\psi }{f}_1\left(t\right),\mathrm{for}\mathrm{each}t\in [a,b].\end{array}$

By $\left(A_2\right)$, we have

$H_d(\mathfrak{F}(t,u),\mathfrak{F}(t,\overline{u}))\le m\left(t\right)|u\left(t\right)-\overline{u}\left(t\right)|.$

So, there exists an element $\upsilon \in \mathfrak{F}(t,\overline{x}\left(t\right))$ such that

$|v_1\left(t\right)-\upsilon |\le m\left(t\right)|u\left(t\right)-\overline{u}\left(t\right)|,t\in [a,b].$

Define $V:[a,b]\to \mathcal{P}\left(\mathbb{R}\right)$ by

$V\left(t\right)=\{\upsilon \in \mathbb{R}:|v_1\left(t\right)-\upsilon |\le m\left(t\right)|u\left(t\right)-\overline{u}\left(t\right)\left|\right\}.$

By Proposition III.4 [33], the multivalued operator $V\left(t\right)\cap \mathfrak{F}(t,\overline{u}\left(t\right))$ is measurable, and thus we can find a measurable selection $v_2\left(t\right)$ for V. So $v_2\left(t\right)\in \mathfrak{F}(t,\overline{u}\left(t\right))$ and $|v_1\left(t\right)-v_2\left(t\right)|\le m\left(t\right)|u\left(t\right)-\overline{u}\left(t\right)|$, for each $t\in [a,b]$. Let us define

$\begin{array}{c}\hfill h_2\left(t\right)=\frac{{(\psi \left(t\right)-\psi \left(a\right))}^{\frac{{\theta }_k}{k}-1}}{\Omega {\Gamma }_k\left({\theta }_k\right)}\left[\lambda {}^k{\mathfrak{I}}^{\alpha -p;\psi }{f}_2\left(\eta \right)+\mu {}^k{\mathfrak{I}}^{\alpha +v;\psi }{f}_2\left(\sigma \right)-{}^k{\mathfrak{I}}^{\alpha ;\psi }{f}_2\left(b\right)\right]+{}^k{\mathfrak{I}}^{\alpha ;\psi }{f}_2\left(t\right),\end{array}$

for each $t\in [a,b]$. Consequently, we have

$\begin{array}{ccc}& & |h_1\left(t\right)-h_2\left(t\right)|\hfill \\ & \le & \frac{{(\psi \left(t\right)-\psi \left(a\right))}^{\frac{{\theta }_k}{k}-1}}{\Omega {\Gamma }_k\left({\theta }_k\right)}\left[\right|\lambda |{}^k{\mathfrak{I}}^{\alpha -p;\psi }\left(\right|f_1\left(s\right)-f_2\left(s\right)\left|\right)\left(\eta \right)\hfill \\ & & +\left|\mu \right|{}^k{\mathfrak{I}}^{\alpha +v;\psi }\left(\right|f_1\left(s\right)-f_2\left(s\right)\left|\right)\left(\sigma \right)+{}^k{\mathfrak{I}}^{\alpha ;\psi }\left(\right|f_1\left(s\right)-f_2\left(s\right)\left|\right)\left(b\right)]\hfill \\ & & +{}^k{\mathfrak{I}}^{\alpha ;\psi }\left(\right|f_1\left(s\right)-f_2\left(s\right)\left|\right)\left(t\right)\hfill \\ & \le & \{\frac{{(\psi \left(b\right)-\psi \left(a\right))}^{\frac{\alpha }{k}}}{{\Gamma }_k(\alpha +k)}+\frac{{(\psi \left(b\right)-\psi \left(a\right))}^{\frac{{\theta }_k}{k}-1}}{|\Omega |{\Gamma }_k\left({\theta }_k\right)}[\left|\lambda \right|\frac{{(\psi \left(\eta \right)-\psi \left(a\right))}^{\frac{\alpha -p}{k}}}{{\Gamma }_k(\alpha -p+k)}\hfill \\ & & +\left|\mu \right|\frac{{(\psi \left(\sigma \right)-\psi \left(a\right))}^{\frac{\alpha +v}{k}}}{{\Gamma }_k(\alpha +v+k)}+\frac{{(\psi \left(b\right)-\psi \left(a\right))}^{\frac{\alpha }{k}}}{{\Gamma }_k(\alpha +k)}\left]\right\}\parallel m\parallel \parallel u-\overline{u}\parallel \hfill \\ & =& \mathbb{G}\parallel m\parallel \parallel u-\overline{u}\parallel ,\hfill \end{array}$

which yields

$\parallel h_1-h_2\parallel \le \mathbb{G}\parallel m\parallel \parallel u-\overline{u}\parallel .$

Similarly, switching the roles of u and $\overline{u}$, we can obtain

$\begin{array}{c}\hfill H_d(\mathcal{F}\left(u\right),\mathcal{F}\left(\overline{u}\right))\le \mathbb{G}\parallel m\parallel \parallel u-\overline{u}\parallel .\end{array}$

Hence, $\mathcal{F}$ is a contraction and we deduce by Covitz and Nadler fixed point theorem [32] that $\mathcal{F}$ has a fixed point $u,$ which is indeed a solution of the multi-valued problem (2). This ends the proof. □

6. Examples

In this section, numerical examples illustrating the applicability of our theoretical results are presented.

Example 1.

Consider the following $(k,\psi )$-Hilfer-type fractional differential equations and inclusions with mixed $(k,\psi )$-derivative and integral boundary conditions of the form

(23) $\left\{\begin{array}{c}{\displaystyle {}^{\frac{3}{2},H}{D}^{\frac{7}{4},\frac{1}{4};\frac{t+1}{t+2}}u\left(t\right)=\mathfrak{f}(t,u\left(t\right)),\frac{1}{5}<t<\frac{8}{5},}\hfill \\ {\displaystyle u\left(\frac{1}{5}\right)=0,u\left(\frac{8}{5}\right)=\frac{2}{31}{}^{\frac{3}{2},H}{D}^{\frac{3}{4},\frac{1}{2};\frac{t+1}{t+2}}u\left(\frac{3}{5}\right)+\frac{3}{51}{}^{\frac{3}{2}}{\mathfrak{I}}^{\frac{15}{16},\frac{t+1}{t+2}}u\left(\frac{6}{5}\right).}\hfill \end{array}\right.$

Here, we choose $k=3/2$, $\alpha =7/4$, $\beta =1/4$, $\psi \left(t\right)=(t+1)/(t+2)$, $a=1/5$, $b=8/5$, $\lambda =2/31$, $p=3/4$, $q=1/2$, $\eta =3/5$, $\mu =3/51$, $v=15/16$, $\sigma =6/5$. The computational yields ${\theta }_{\frac{3}{2}}=33/16$, $\Omega \approx 0.3485295041$, $\mathbb{G}\approx 0.2007468250$, ${\mathbb{G}}_1\approx 0.1245096496$.

(i) Let the nonlinear unbounded Lipschitzian function $\mathfrak{f}:[1/5,8/5]\times \mathbb{R}\to \mathbb{R}$ be given by

(24)$\mathfrak{f}(t,u)=\left(\frac{u^2+4\left|u\right|}{1+\left|u\right|}\right){sin}^2\pi t+\frac{1}{3}{t}^2+\frac{1}{4}t+\frac{1}{5}.$

Now, we see that $\mathfrak{f}$ satisfies the Lipschitzian condition as

$|\mathfrak{f}(t,u_1)-\mathfrak{f}(t,u_2)|\le 4|u_1-u_2|,$

for each $u_1,u_2\in \mathbb{R}$ and $t\in [1/5,8/5]$ with $\mathfrak{L}=4$. Therefore, we obtain $\mathfrak{L}\mathbb{G}\approx 0.802987300<1$, which means that the condition (19) is satisfied. Hence, by Theorem 1, the $(k,\psi )$-Hilfer-type fractional differential equation with mixed $(k,\psi )$-derivative and integral boundary conditions (23) with the function f given by (24), has a unique solution on the interval $[1/5,8/5]$.

(ii) Suppose that a nonlinear bounded Lipschitzian function $\mathfrak{f}:[1/5,8/5]\times \mathbb{R}\to \mathbb{R}$ is

(25)$\mathfrak{f}(t,u)=\left(\frac{8\left|u\right|}{1+\left|u\right|}\right){cos}^4\pi t+3t+2,$

which is bounded as

$\left|\mathfrak{f}(t,u)\right|\le 8{cos}^4\pi t+3t+2:=\theta \left(t\right),$

for all $t\in [1/5,8/5]$. Observe that $\mathfrak{f}$ satisfies the Lipschitz condition: $|\mathfrak{f}(t,u_1)-\mathfrak{f}(t,u_2)|\le 8|u_1-u_2|,$ with Lipschitz constant $\mathfrak{L}=8$. Hence, we obtain $\mathfrak{L}{\mathfrak{G}}_{\mathfrak{1}}\approx 0.9960771968<1$. Then, the result in Theorem 2 yields that the $(k,\psi )$-Hilfer-type boundary value problem (23) with $\mathfrak{f}$ presented by (25), has at least one solution on $[1/5,8/5]$. As $\mathfrak{L}\mathbb{G}\approx 1.605974600>1$, the uniqueness result (Theorem 1) does not apply in this situation.

(iii) Assume that

(26)$\mathfrak{f}(t,u)=\theta \left(t\right)\left(A{g}_1\left(u\right)+B\right),$

where $\theta :[1/5,8/5]\to {\mathbb{R}}^{+}$, $g_1:\mathbb{R}\to \mathbb{R}$ with $|g_1\left(u\right)|\le |u|$, for example, $g_1\left(u\right)=u{sin}^{8}u$. In addition, $0\le A<1/(\parallel \theta \parallel \mathbb{G})$ and $B>0$. Then we have

$\left|\mathfrak{f}(t,u)\right|\le \parallel \theta \parallel \left(A\left|u\right|+B\right).$

Writing $\chi \left(u\right)=A\left|u\right|+B$, there exists a constant $\mathfrak{K}$ satisfying condition $\left(H_4\right)$ in Theorem 3:

$\mathfrak{K}>\frac{B\parallel \theta \parallel \mathbb{G}}{1-A\parallel \theta \parallel \mathbb{G}}.$

Therefore, by applying Theorem 3, we deduce that the boundary value problem of $(k,\psi )$-Hilfer-type fractional differential equation (23) with $\mathfrak{f}$ given in (26) has at least one solution on $[1/5,8/5].$

(iv) Suppose that

(27)$\mathfrak{f}(t,u)=\theta \left(t\right)\left(A{g}_2\left(u\right)+B\right),$

where $\theta :[1/5,8/5]\to {\mathbb{R}}^{+}$, $g_2:\mathbb{R}\to \mathbb{R}$ with $|g_2\left(u\right)|\le u^2$, for example, $g_2\left(u\right)=u^{18}/(1+u^{16})$. In addition, two positive constants A and B are such that $AB<1/(4\parallel \theta {\parallel }^2{\mathbb{G}}^2)$. Then, we have

$\left|\mathfrak{f}(t,u)\right|\le \parallel \theta \parallel \left(A{u}^2+B\right).$

Setting a function $\chi \left(u\right)=A{u}^2+B$, there exists a constant

$\mathfrak{K}\in \left(\frac{1-\sqrt{1-{4AB\parallel \theta \parallel }^2{\mathbb{G}}^2}}{2A\parallel \theta \parallel \mathbb{G}},\frac{1+\sqrt{1-{4AB\parallel \theta \parallel }^2{\mathbb{G}}^2}}{2A\parallel \theta \parallel \mathbb{G}}\right),$

satisfying condition $\left(H_4\right)$ in Theorem 3. Therefore, all conditions of Theorem 3 are satisfied. Hence $(k,\psi )$-Hilfer-type nonlocal fractional boundary value problem (23), with $\mathfrak{f}$ given in (27) has at least one solution on $[1/5,8/5]$.

(v) Suppose that the first equation of (23) is presented by

(28)${}^{\frac{3}{2},H}{D}^{\frac{7}{4},\frac{1}{4};\frac{t+1}{t+2}}u\left(t\right)\in \mathfrak{F}(t,u\left(t\right)),\frac{1}{5}<t<\frac{8}{5},$

where

$\mathfrak{F}(t,u)=\left[0,\frac{1}{5t+1}\left(\frac{u^2+6\left|u\right|}{1+\left|u\right|}+\frac{1}{4}{sin}^{2}t+\frac{3}{4}{e}^{-|5t-1|}\right)\right].$

Obviously $\mathfrak{F}(t,u)$ is a measurable set. Moreover,

$H_d\left(\mathfrak{F}(t,u),\mathfrak{F}(t,\overline{u})\right)\le \frac{6}{5t+1}|u-\overline{u}|.$

Let us choose $m\left(t\right)=6/(5t+1)$ and observe that $d(0,\mathfrak{F}(t,0\left)\right)\le 1/(5t+1)<6/(5t+1)=m\left(t\right)$ for almost all $t\in [1/8,5/8]$. Since $\delta =\mathbb{G}\parallel m\parallel \approx 0.6022404750<1$, the $(k,\psi )$-Hilfer-type nonlocal fractional inclusion (28) with mixed $(k,\psi )$-derivative and integral boundary conditions as in (23), has at least one solution on $[1/5,8/5]$.

7. Conclusions

This research is devoted to the analysis of single-valued and multi-valued $(k,\psi )$-Hilfer-type nonlocal fractional boundary value problems involving $(k,\psi )$-Hilfer fractional derivative and integral operators in boundary conditions. We established existence and uniqueness results for the single-valued case after transforming the given problem into a fixed point problem, with the help of Banach contraction mapping principle, Krasnosel’skiĭ fixed point theorem and the Leray–Schauder nonlinear alternative. Two existence results for the multi-valued problem are obtained by applying Leray–Schauder nonlinear alternative for multivalued maps and Covitz–Nadler fixed point theorem for contractive multivalued maps, respectively, for the cases of convex and non-convex multivalued map involved in the problem. All the obtained theoretical results are well-illustrated by numerical examples. The results are new and enrich the new research area on $(k,\psi )$-Hilfer nonlocal fractional boundary value problems in the order of $(1,2].$ In a future study, we plan to extend the results of this paper to cover $(k,\psi )$-Hilfer nonlocal fractional systems of order in $(1,2].$

Author Contributions

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

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.

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

References

1. Diethelm, K. The Analysis of Fractional Differential Equations; Lecture Notes in Mathematics Springer: New York, NY, USA, 2010.

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.

3. Lakshmikantham, V.; Leela, S.; Devi, J.V. Theory of Fractional Dynamic Systems; Cambridge Scientific Publishers: Cambridge, UK, 2009.

4. Miller, K.S.; Ross, B. An Introduction to the Fractional Calculus and Differential Equations; John Wiley: NewYork, NY, USA, 1993.

5. Podlubny, I. Fractional Differential Equations; Academic Press: New York, NY, USA, 1999.

6. Samko, S.G.; Kilbas, A.A.; Marichev, O.I. Fractional Integrals and Derivatives; Gordon and Breach Science: Yverdon, Switzerland, 1993.

7. Ahmad, B.; Alsaedi, A.; Ntouyas, S.K.; Tariboon, J. Hadamard-Type Fractional Differential Equations, Inclusions and Inequalities; Springer: Cham, Switzerland, 2017.

8. Ahmad, B.; Ntouyas, S.K. Nonlocal Nonlinear Fractional-Order Boundary Value Problems; World Scientific: Singapore, 2021.

9. Zhou, Y. Basic Theory of Fractional Differential Equations; World Scientific: Singapore, 2014.

10. Katugampola, U.N. New approach to a generalized fractional integral. Appl. Math. Comput.; 2011; 218, pp. 860-865. [DOI: https://dx.doi.org/10.1016/j.amc.2011.03.062]

11. Katugampola, U.N. A new approach to generalized fractional derivatives. Bull. Math. Anal. Appl.; 2014; 6, pp. 1-154.

12. Hilfer, R. Applications of Fractional Calculus in Physics; World Scientific: Singapore, 2000.

13. Almeida, R. A Caputo fractional derivative of a function with respect to another function. Commun. Nonlinear Sci. Numer. Simul.; 2017; 44, pp. 460-481. [DOI: https://dx.doi.org/10.1016/j.cnsns.2016.09.006]

14. Kucche, K.D.; Mali, A.D. On the nonlinear (k,ψ)-Hilfer fractional differential equations. Chaos Solitons Fractals; 2021; 152, 111335. [DOI: https://dx.doi.org/10.1016/j.chaos.2021.111335]

15. Torres Ledesma, C.E.; Nyamoradi, N. (k,ψ)-Hilfer variational problem. J. Elliptic Parabol. Equ.; 2022; [DOI: https://dx.doi.org/10.1007/s41808-022-00173-w]

16. Tariboon, J.; Samadi, A.; Ntouyas, S.K. Multi-point boundary value problems for (k,ψ)-Hilfer fractional differential equations and inclusions. Axioms; 2022; 11, 110. [DOI: https://dx.doi.org/10.3390/axioms11030110]

17. Debnath, P.; Konwar, N.; Radenovic, S. Metric Fixed Point Theory—Applications in Science, Engineering and Behavioural Sciences; Forum for Interdisciplinary Mathematics Springer: Singapore, 2021.

18. Todorcevic, V. Harmonic Quasiconformal Mappings and Hyperbolic Type Metrics; Springer: Cham, Switzerland, 2019.

19. Cho, Y.J.; Jleli, M.; Mursaleen, M.; Samet, B.; Vetro, C. Advances in Metric Fixed Point Theory and Applications; Springer: Singapore, 2021.

20. Mubeen, S.; Habibullah, G.M. k–fractional integrals and applications. Int. J. Contemp. Math. Sci.; 2012; 7, pp. 89-94.

21. Diaz, R.; Pariguan, E. On hypergeometric functions and Pochhammer k-symbol. Divulg. Mat.; 2007; 2, pp. 179-192.

22. Dorrego, G.A. An alternative definition for the k-Riemann-Liouville fractional derivative. Appl. Math. Sci.; 2015; 9, pp. 481-491. [DOI: https://dx.doi.org/10.12988/ams.2015.411893]

23. Sousa, J.V.D.C.; de Oliveira, E.C. On the ψ-Hilfer fractional derivative. Commun. Nonlinear Sci. Numer. Simul.; 2018; 60, pp. 72-91. [DOI: https://dx.doi.org/10.1016/j.cnsns.2018.01.005]

24. 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, pp. 64946-64953. [DOI: https://dx.doi.org/10.1109/ACCESS.2018.2878266]

25. Zhang, S.Q. Existence of solution for a boundary value problem of fractional order. Acta Math. Sci.; 2006; 26, pp. 220-228. [DOI: https://dx.doi.org/10.1016/S0252-9602(06)60044-1]

26. Luca, R. On a system of Riemann-Liouville fractional differential equations with coupled nonlocal boundary conditions. Adv. Differ. Equ.; 2021; 2021, 134. [DOI: https://dx.doi.org/10.1186/s13662-021-03303-1]

27. Deimling, K. Nonlinear Functional Analysis; Springer: New York, NY, USA, 1985.

28. Krasnosel’skiĭ, M.A. Two remarks on the method of successive approximations. Uspekhi Mat. Nauk; 1955; 10, pp. 123-127.

29. Granas, A.; Dugundji, J. Fixed Point Theory; Springer: New York, NY, USA, 2005.

30. Deimling, K. Multivalued Differential Equations; Walter De Gruyter: Berlin, Germany, New York, NY, USA, 1992.

31. 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, pp. 781-786.

32. Covitz, H.; Nadler, S.B., Jr. Multivalued contraction mappings in generalized metric spaces. Israel J. Math.; 1970; 8, pp. 5-11. [DOI: https://dx.doi.org/10.1007/BF02771543]

33. Castaing, C.; Valadier, M. Convex Analysis and Measurable Multifunctions; Lecture Notes in Mathematics 580 Springer: Berlin/Heidelberg, Germany, New York, NY, USA, 1977.