Next Article in Journal
Application of Mixed Sampling to Real Life Data: A Case Study on Socio-Economic Determinants by Using SEM and CFA Techniques
Next Article in Special Issue
Integral Representation for the Solutions of Autonomous Linear Neutral Fractional Systems with Distributed Delay
Previous Article in Journal
On the Stability with Respect to H-Manifolds for Cohen–Grossberg-Type Bidirectional Associative Memory Neural Networks with Variable Impulsive Perturbations and Time-Varying Delays
Previous Article in Special Issue
Certain Fractional Proportional Integral Inequalities via Convex Functions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Nonlinear Integro-Differential Equations Involving Mixed Right and Left Fractional Derivatives and Integrals with Nonlocal Boundary Data

1
Nonlinear Analysis and Applied Mathematics (NAAM)-Research Group, Department of Mathematics, Faculty of Science, King Abdulaziz University, P.O. Box 80203, Jeddah 21589, Saudi Arabia
2
Department of Mathematics, University of Ioannina, 451 10 Ioannina, Greece
*
Author to whom correspondence should be addressed.
Mathematics 2020, 8(3), 336; https://doi.org/10.3390/math8030336
Submission received: 2 February 2020 / Revised: 22 February 2020 / Accepted: 25 February 2020 / Published: 3 March 2020

Abstract

:
In this paper, we study the existence of solutions for a new nonlocal boundary value problem of integro-differential equations involving mixed left and right Caputo and Riemann–Liouville fractional derivatives and Riemann–Liouville fractional integrals of different orders. Our results rely on the standard tools of functional analysis. Examples are constructed to demonstrate the application of the derived results.

1. Introduction

In the last few decades, fractional-order single-valued and multivalued boundary value problems containing different fractional derivatives such as Caputo, Riemann–Liouville, Hadamard, etc., and classical, nonlocal, integral boundary conditions have been extensively studied, for example, see the articles [1,2,3,4,5,6,7,8,9,10,11,12] and the references cited therein.
In the study of variational principles, fractional differential equations involving both left and right fractional derivatives give rise to a special class of Euler–Lagrange equations, for details, see [13] and the references cited therein. Let us consider some works on mixed fractional-order boundary value problems. In [14], the authors discussed the existence of an extremal solution to a nonlinear system involving the right-handed Riemann–Liouville fractional derivative. In [15], a two-point nonlinear higher order fractional boundary value problem involving left Riemann–Liouville and right Caputo fractional derivatives was investigated, while a problem in terms of left Caputo and right Riemann–Liouville fractional derivatives was studied in [16]. A nonlinear fractional oscillator equation containing left Riemann–Liouville and right Caputo fractional derivatives was investigated in [17]. In a recent paper [18], the authors proved some existence results for nonlocal boundary value problems of differential equations and inclusions containing both left Caputo and right Riemann–Liouville fractional derivatives.
Integro-differential equations appear in the mathematical modeling of several real world problems such as, heat transfer phenomena [19,20], forced-convective flow over a heat-conducting plate [21], etc. In [22], the authors studied the steady heat-transfer in fractal media via the local fractional nonlinear Volterra integro-differential equations. Electromagnetic waves in a variety of dielectric media with susceptibility following a fractional power-law are described by the fractional integro-differential equations [23].
Motivated by aforementioned applications of integro-differential equations and [18], we introduce a new kind of integro-differential equation involving right-Caputo and left-Riemann–Liouville fractional derivatives of different orders and right-left Riemann–Liouville fractional integrals and solve it subject to nonlocal boundary conditions. In precise terms, we prove existence and uniqueness of solutions for the problem given by
C D 1 α R L D 0 + β y ( t ) + λ I 1 p I 0 + q h ( t , y ( t ) ) = f ( t , y ( t ) ) , t J : = [ 0 , 1 ] ,
y ( 0 ) = y ( ξ ) = 0 , y ( 1 ) = δ y ( μ ) , 0 < ξ < μ < 1 ,
where C D 1 α and R L D 0 + β denote the right Caputo fractional derivative of order α ( 1 , 2 ] and the left Riemann–Liouville fractional derivative of order β ( 0 , 1 ] , I 1 p and I 0 + q denote the right and left Riemann–Liouville fractional integrals of orders p , q > 0 respectively, f , h : [ 0 , 1 ] × R R are given continuous functions and δ , λ R . It is imperative to notice that the integro-differential equation in (1) and (2) contains mixed type (integral and nonintegral) nonlinearities.
We organize the rest of the paper as follows. Section 2 contains some preliminary concepts related to our work. In Section 3, we prove an auxiliary lemma for the linear variant of the problem (1) and (2). Then we derive the existence results for the problem (1) and (2) by applying a fixed point theorem due to Krasnoselski and Leray–Schauder nonlinear alternative, while the uniqueness result is established via Banach contraction mapping principle. Examples illustrating the main results are also presented.

2. Preliminaries

In this section, we recall some related definitions of fractional calculus [1].
Definition 1.
The left and right Riemann–Liouville fractional integrals of order β > 0 for an integrable function g : ( 0 , ) R are respectively defined by
I 0 + β g ( t ) = 0 t ( t s ) β 1 Γ ( β ) g ( s ) d s a n d I 1 β g ( t ) = t 1 ( s t ) β 1 Γ ( β ) g ( s ) d s .
Definition 2.
The left Riemann–Liouville fractional derivative and the right Caputo fractional derivative of order β ( n 1 , n ] , n N for a function g : ( 0 , ) R with g C n ( ( 0 , ) , R ) are respectively given by
D 0 + β g ( t ) = d n d t n 0 t ( t s ) n β 1 Γ ( n β ) g ( s ) d s a n d c D 1 β g ( t ) = ( 1 ) n t 1 ( s t ) n β 1 Γ ( n β ) g ( n ) ( s ) d s .
Lemma 1.
If p > 0 and q > 0 , then the following relations hold almost everywhere on [ a , b ] :
I 1 p I 1 q f ( x ) = I 1 p + q f ( x ) , I 0 + p I 0 + q f ( x ) = I 0 + p + q f ( x ) .

3. Main Results

In the following lemma, we solve a linear variant of the problem (1) and (2).
Lemma 2.
Let H , F C [ 0 , 1 ] L ( 0 , 1 ) . Then the linear problem
C D 1 α R L D 0 + β y ( t ) + λ I 1 p I 0 + q H ( t ) = F ( t ) , t J : = [ 0 , 1 ] , y ( 0 ) = y ( ξ ) = 0 , y ( 1 ) = δ y ( μ ) ,
is equivalent to the fractional integral equation:
y ( t ) = 0 t ( t s ) β 1 Γ ( β ) I 1 α F ( s ) λ I 1 α + p I 0 + q H ( s ) d s + a 1 ( t ) { δ 0 μ ( μ s ) β 1 Γ ( β ) I 1 α F ( s ) λ I 1 α + p I 0 + q H ( s ) d s 0 1 ( 1 s ) β 1 Γ ( β ) I 1 α F ( s ) λ I 1 α + p I 0 + q H ( s ) d s } + a 2 ( t ) 0 ξ ( ξ s ) β 1 Γ ( β ) I 1 α F ( s ) λ I 1 α + p I 0 + q H ( s ) d s ,
where
a 1 ( t ) = 1 Λ ξ β + 1 t β ξ β t β + 1 , a 2 ( t ) = 1 Λ t β ( 1 δ μ β + 1 ) t β + 1 ( 1 δ μ β ) ,
and it is assumed that
Λ = ξ β + 1 ( 1 δ μ β ) ξ β ( 1 δ μ β + 1 ) 0 .
Proof. 
Applying the left and right fractional integrals I 1 α and I 0 + β successively to the integro-differential equation in (3), and then using Lemma 1, we get
y ( t ) = I 0 + β I 1 α F ( t ) λ I 0 + β I 1 α + p I 0 + q H ( t ) + c 0 t β Γ ( β + 1 ) + c 1 t β + 1 Γ ( β + 2 ) + c 2 t β 1 ,
where c 0 , c 1 and c 2 are unknown arbitrary constants.
In view of the condition y ( 0 ) = 0 , it follows from (7) that c 2 = 0 . Inserting c 2 = 0 in (7) and then using the nonlocal boundary conditions y ( ξ ) = 0 , y ( 1 ) = δ y ( μ ) in the resulting equation, we obtain a system of equations in c 0 and c 1 given by
c 0 1 δ μ β Γ ( β + 1 ) + c 1 1 δ μ β + 1 Γ ( β + 2 ) = δ A 1 A 2 , c 0 ξ β Γ ( β + 1 ) + c 1 ξ β + 1 Γ ( β + 2 ) = A 3 ,
where
A 1 = I 0 + β I 1 α F ( μ ) λ I 0 + β I 1 α + p I 0 + q H ( μ ) , A 2 = I 0 + β I 1 α F ( 1 ) λ I 0 + β I 1 α + p I 0 + q H ( 1 ) , A 3 = I 0 + β I 1 α F ( ξ ) λ I 0 + β I 1 α + p I 0 + q H ( ξ ) .
Solving the system (8), we find that
c 0 = Γ ( β + 1 ) Λ ξ β + 1 δ A 1 A 2 + ( 1 δ μ β + 1 ) A 3 , c 1 = Γ ( β + 2 ) Λ ξ β δ A 1 A 2 + ( 1 δ μ β ) A 3 ,
where Λ is defined by (6). Substituting the values of c 0 and c 1 together with the notations (5) in (7), we obtain the solution (4). The converse follows by direct computation. This completes the proof. □
Let X = C ( [ 0 , 1 ] , R ) denote the Banach space of all continuous functions from [ 0 , 1 ] R equipped with the norm y = sup { | y ( t ) | : t [ 0 , 1 ] } . By Lemma 2, we define an operator G : X X associated with the problem (1) and (2) as
G y ( t ) = 0 t ( t s ) β 1 Γ ( β ) I 1 α f ( s , y ( s ) ) λ I 1 α + p I 0 + q h ( s , y ( s ) ) d s + a 1 ( t ) [ δ 0 μ ( μ s ) β 1 Γ ( β ) I 1 α f ( s , y ( s ) ) λ I 1 α + p I 0 + q h ( s , y ( s ) ) d s 0 1 ( 1 s ) β 1 Γ ( β ) I 1 α f ( s , y ( s ) ) λ I 1 α + p I 0 + q h ( s , y ( s ) ) d s ] + a 2 ( t ) 0 ξ ( ξ s ) β 1 Γ ( β ) I 1 α f ( s , y ( s ) ) λ I 1 α + p I 0 + q h ( s , y ( s ) ) d s .
Notice that the fixed points of the operator G are solutions of the problem (1) and (2).
In the forthcoming analysis, we use the following estimates:
0 t ( t s ) β 1 Γ ( β ) I 1 α + p I 0 + q d s = 0 t ( t s ) β 1 Γ ( β ) s 1 ( u s ) α + p 1 Γ ( α + p ) 0 u ( u r ) q 1 Γ ( q ) d r d u d s t β Γ ( β + 1 ) Γ ( α + p + 1 ) Γ ( q + 1 ) , 0 t ( t s ) β 1 Γ ( β ) I 1 α d s = 0 t ( t s ) β 1 Γ ( β ) s 1 ( u s ) α 1 Γ ( α ) d u d s t β Γ ( α + 1 ) Γ ( β + 1 ) ,
where we have used u q 1 , ( 1 s ) α + p < 1 ; p , q > 0 , 1 < α 2 .
In the sequel, we set
Ω 1 = Δ Γ ( α + 1 ) , Ω 2 = | λ | Δ Γ ( α + p + 1 ) Γ ( q + 1 ) ,
where
Δ = 1 Γ ( β + 1 ) 1 + a ¯ 1 ( | δ | μ β + 1 ) + a ¯ 2 ξ β ,
a ¯ 1 = max t [ 0 , 1 ] | a 1 ( t ) | , a ¯ 2 = max t [ 0 , 1 ] | a 2 ( t ) | .

3.1. Existence Results

In the following, we prove our first existence result for the problem (1) and (2), which relies on Krasnoselskii’s fixed point theorem [24].
Theorem 1.
Assumed that:
  • ( B 1 ) There exist L > 0 such that | f ( t , x ) f ( t , y ) | L | x y | , t [ 0 , 1 ] , x , y R ;
  • ( B 2 ) There exist K > 0 such that | h ( t , x ) h ( t , y ) | K | x y | , t [ 0 , 1 ] , x , y R ;
  • ( B 3 ) | f ( t , y ) | σ ( t ) and | h ( t , y ) | ρ ( t ) , where σ , ρ C ( [ 0 , 1 ] , R + ) .
Then the problem (1) and (2) has at least one solution on [ 0 , 1 ] if L γ 1 + K γ 2 < 1 , where
γ 1 = 1 Γ ( β + 1 ) Γ ( α + 1 ) , γ 2 = | λ | Γ ( β + 1 ) Γ ( α + p + 1 ) Γ ( q + 1 ) .
Proof. 
Introduce the ball B θ = { y X : y θ } , where σ = sup t [ 0 , 1 ] | σ ( t ) | , ρ = sup t [ 0 , 1 ] | ρ ( t ) | and
θ σ Ω 1 + ρ Ω 2 .
Let us split the operator G : X X on B θ as G = G 1 + G 2 , where
G 1 y ( t ) = 0 t ( t s ) β 1 Γ ( β ) I 1 α f ( s , y ( s ) ) d s λ 0 t ( t s ) β 1 Γ ( β ) I 1 α + p I 0 + q h ( s , y ( s ) ) d s , G 2 y ( t ) = a 1 ( t ) [ δ 0 μ ( μ s ) β 1 Γ ( β ) I 1 α f ( s , y ( s ) ) d s λ 0 μ ( μ s ) β 1 Γ ( β ) I 1 α + p I 0 + q h ( s , y ( s ) ) d s 0 1 ( 1 s ) β 1 Γ ( β ) I 1 α f ( s , y ( s ) ) d s λ 0 1 ( 1 s ) β 1 Γ ( β ) I 1 α + p I 0 + q h ( s , y ( s ) ) d s ] + a 2 ( t ) 0 ξ ( ξ s ) β 1 Γ ( β ) I 1 α f ( s , y ( s ) ) d s λ 0 ξ ( ξ s ) β 1 Γ ( β ) I 1 α + p I 0 + q h ( s , y ( s ) ) d s .
Now, we verify that the operators G 1 and G 2 satisfy the hypothesis of Krasnoselskii’s theorem [24] in three steps.
(i) For y , x B θ , we have
G 1 y + G 2 x sup t [ 0 , 1 ] { 0 t ( t s ) β 1 Γ ( β ) I 1 α | f ( s , y ( s ) ) | d s + | λ | 0 t ( t s ) β 1 Γ ( β ) I 1 α + p I 0 + q | h ( s , y ( s ) ) | d s + | a 1 ( t ) | { | δ | 0 μ ( μ s ) β 1 Γ ( β ) I 1 α | f ( s , x ( s ) ) | d s + | λ | 0 μ ( μ s ) β 1 Γ ( β ) I 1 α + p I 0 + q | h ( s , x ( s ) ) | d s + 0 1 ( 1 s ) β 1 Γ ( β ) I 1 α | f ( s , x ( s ) ) | d s + | λ | 0 1 ( 1 s ) β 1 Γ ( β ) I 1 α + p I 0 + q | h ( s , x ( s ) ) | d s } + | a 2 ( t ) | 0 ξ ( ξ s ) β 1 Γ ( β ) I 1 α | f ( s , x ( s ) ) | d s + | λ | 0 ξ ( ξ s ) β 1 Γ ( β ) I 1 α + p I 0 + q | h ( s , x ( s ) ) | d s } σ sup t [ 0 , 1 ] { 0 t ( t s ) β 1 Γ ( β ) I 1 α d s + | a 1 ( t ) | [ | δ | 0 μ ( μ s ) β 1 Γ ( β ) I 1 α d s + 0 1 ( 1 s ) β 1 Γ ( β ) I 1 α d s ] + | a 2 ( t ) | 0 ξ ( ξ s ) β 1 Γ ( β ) I 1 α d s } + ρ | λ | sup t [ 0 , 1 ] { 0 t ( t s ) β 1 Γ ( β ) I 1 α + p I 0 + q d s + | a 1 ( t ) | [ | δ | 0 μ ( μ s ) β 1 Γ ( β ) I 1 α + p I 0 + q d s + 0 1 ( 1 s ) β 1 Γ ( β ) I 1 α + p I 0 + q d s ] + | a 2 ( t ) | 0 ξ ( ξ s ) β 1 Γ ( β ) I 1 α + p I 0 + q d s } σ Γ ( α + 1 ) + ρ | λ | Γ ( α + p + 1 ) Γ ( q + 1 ) Δ = σ Ω 1 + ρ Ω 2 < θ ,
where we used (11). Thus G 1 y + G 2 x B θ .
(ii) Using ( B 1 ) and ( B 2 ), it is easy to show that
G 1 y G 1 x sup t [ 0 , 1 ] { 0 t ( t s ) β 1 Γ ( β ) I 1 α | f ( s , y ( s ) ) f ( s , x ( s ) ) | d s + | λ | 0 t ( t s ) β 1 Γ ( β ) I 1 α + p I 0 + q | h ( s , y ( s ) ) h ( s , x ( s ) ) | d s } ( L γ 1 + K γ 2 ) y x ,
which, in view of the condition: L γ 1 + K γ 2 < 1 , implies that the operator G 1 is a contraction.
(iii) Continuity of the functions f , h implies that the operator G 2 is continuous. In addition, G 2 is uniformly bounded on B θ as
G 2 y sup t [ 0 , 1 ] { | a 1 ( t ) | { | δ | 0 μ ( μ s ) β 1 Γ ( β ) I 1 α | f ( s , y ( s ) ) | + | λ | I 1 α + p I 0 + q | h ( s , y ( s ) ) | d s + 0 1 ( 1 s ) β 1 Γ ( β ) I 1 α | f ( s , y ( s ) ) | + | λ | I 1 α + p I 0 + q | h ( s , y ( s ) ) | d s } + | a 2 ( t ) | 0 ξ ( ξ s ) β 1 Γ ( β ) I 1 α | f ( s , y ( s ) ) | + | λ | I 1 α + p I 0 + q | h ( s , y ( s ) ) | d s } σ sup t [ 0 , 1 ] { | a 1 ( t ) | | δ | 0 μ ( μ s ) β 1 Γ ( β ) I 1 α d s + 0 1 ( 1 s ) β 1 Γ ( β ) I 1 α d s + | a 2 ( t ) | 0 ξ ( ξ s ) β 1 Γ ( β ) I 1 α d s } + ρ | λ | sup t [ 0 , 1 ] { | a 1 ( t ) | | δ | 0 μ ( μ s ) β 1 Γ ( β ) I 1 α + p I 0 + q d s + 0 1 ( 1 s ) β 1 Γ ( β ) I 1 α + p I 0 + q d s + | a 2 ( t ) | 0 ξ ( ξ s ) β 1 Γ ( β ) I 1 α + p I 0 + q d s } σ ( Ω 1 γ 1 ) + ρ ( Ω 2 γ 2 ) ,
where Ω i , and γ i ( i = 1 , 2 ) are defined by (9) and (10) respectively.
To show the compactness of G 2 , we fix sup ( t , y ) [ 0 , 1 ] × B θ | f ( t , y ) | = f ¯ , sup ( t , y ) [ 0 , 1 ] × B θ | h ( t , y ) | = h ¯ . Then, for 0 < t 1 < t 2 < 1 , we have
| ( G 2 y ) ( t 2 ) ( G 2 y ) ( t 1 ) | | a 1 ( t 2 ) a 1 ( t 1 ) | { | δ | 0 μ ( μ s ) β 1 Γ ( β ) I 1 α | f ( s , y ( s ) ) | + | λ | I 1 α + p I 0 + q | h ( s , y ( s ) ) | d s + 0 1 ( 1 s ) β 1 Γ ( β ) I 1 α | f ( s , y ( s ) ) | + | λ | I 1 α + p I 0 + q | h ( s , y ( s ) ) | d s } + | a 2 ( t 2 ) a 2 ( t 1 ) | 0 ξ ( ξ s ) β 1 Γ ( β ) I 1 α | f ( s , y ( s ) ) | + | λ | I 1 α + p I 0 + q | h ( s , y ( s ) ) | d s γ 1 f ¯ + γ 2 h ¯ { ξ β + 1 | t 2 β t 1 β | + ξ β | t 1 β + 1 t 2 β + 1 | ( | δ | μ β + 1 ) | Λ | + | 1 δ μ β + 1 | | t 2 β t 1 β | + | 1 δ μ β | | t 1 β + 1 t 2 β + 1 | ξ β | Λ | } ,
which tends to zero independent of y as t 2 t 1 . This shows that G 2 is equicontinuous. It is clear from the foregoing arguments that the operator G 2 is relatively compact on B θ . Hence, by the Arzelá-Ascoli theorem, G 2 is compact on B θ .
In view of the foregoing arguments (i)-(iii), the hypothesis of the Krasnoselskii’s fixed point theorem [24] holds true. Thus, the operator G 1 + G 2 = G has a fixed point, which implies that the problem (1) and (2) has at least one solution on [ 0 , 1 ] . The proof is finished. □
Remark 1.
If we interchange the roles of the operators G 1 and G 2 in the previous result, the condition L γ 1 + K γ 2 < 1 , is replaced with the following one:
L ( Ω 1 γ 1 ) + K ( Ω 2 γ 2 ) < 1 ,
where Ω 1 , Ω 2 and γ 1 , γ 2 are defined by (9), (10) respectively.
The following existence result relies on Leray–Schauder nonlinear alternative [25].
Theorem 2.
Suppose that the following conditions hold:
  • ( B 4 ) There exist continuous nondecreasing functions ϕ , ψ : [ 0 , ) ( 0 , ) such that ( t , y ) [ 0 , 1 ] × R , | f ( t , y ) | ω 1 ( t ) ϕ ( y ) and | h ( t , y ) | ω 2 ( t ) ψ ( y ) , where ω 1 , ω 2 C ( [ 0 , T ] , R + ) .
  • ( B 5 ) There exist a constant M > 0 such that
    M ω 1 ϕ ( M ) Ω 1 + ω 2 ψ ( M ) Ω 2 > 1 ,
Then, the problem (1) and (2) has at least one solution on [ 0 , 1 ] .
Proof. 
First we show that the operator G is completely continuous. This will be established in several steps.
(i) G maps bounded sets into bounded sets in X .
Let y B r = { y X : y r } , where r is a fixed number. Then, using the strategy employed in the proof of Theorem 1, we obtain
G y ( t ) ω 1 ϕ ( r ) Γ ( α + 1 ) + ω 2 ψ ( r ) | λ | Γ ( α + p + 1 ) Γ ( q + 1 ) Δ = ω 1 ϕ ( r ) Ω 1 + ω 2 ψ ( r ) Ω 2 < .
(ii) G maps bounded sets into equicontinuous sets.
Let 0 < t 1 < t 2 < 1 and y B r , where B r is bounded set in X . Then we obtain
| G y ( t 2 ) G y ( t 1 ) | | 0 t 1 ( t 2 s ) β 1 ( t 1 s ) β 1 Γ ( β ) I 1 α | f ( s , y ( s ) ) | + | λ | I 1 α + p I 0 + q | h ( s , y ( s ) ) | d s | + | t 1 t 2 ( t 2 s ) β 1 Γ ( β ) I 1 α | f ( s , y ( s ) ) | + | λ | I 1 α + p I 0 + q | h ( s , y ( s ) ) | d s | + | a 1 ( t 2 ) a 1 ( t 1 ) | { | δ | 0 μ ( μ s ) β 1 Γ ( β ) I 1 α | f ( s , y ( s ) ) | + | λ | I 1 α + p I 0 + q | h ( s , y ( s ) ) | d s + 0 1 ( 1 s ) β 1 Γ ( β ) I 1 α | f ( s , y ( s ) ) | + | λ | I 1 α + p I 0 + q | h ( s , y ( s ) ) | d s } + | a 2 ( t 2 ) a 2 ( t 1 ) | 0 ξ ( ξ s ) β 1 Γ ( β ) I 1 α | f ( s , y ( s ) ) | + | λ | I 1 α + p I 0 + q | h ( s , y ( s ) ) | d s ω 1 ϕ ( r ) Γ ( β + 1 ) Γ ( α + 1 ) + ω 2 ψ ( r ) | λ | Γ ( β + 1 ) Γ ( α + p + 1 ) Γ ( q + 1 ) × { 2 ( t 2 t 1 ) β + | t 2 β t 1 β | + ( | δ | μ β + 1 ) | Λ | ξ β + 1 | t 2 β t 1 β | + ξ β | t 2 β + 1 t 1 β + 1 | + ξ β | Λ | | 1 δ μ β + 1 | | t 2 β t 1 β | + | 1 δ μ β | | t 2 β + 1 t 1 β + 1 | } .
Notice that the right-hand side of the above inequality tends to 0 as t 2 t 1 , independent of y B r . In view of the foregoing arguments, it follows by the Arzelá–Ascoli theorem that G : X X is completely continuous.
The conclusion of the Leray–Schauder nonlinear alternative [25] will be applicable once it is shown that there exists an open set U C ( [ 0 , 1 ] , R ) with y ν G y for ν ( 0 , 1 ) and y U . Let y C ( [ 0 , 1 ] , R ) such that y = ν G y for ν ( 0 , 1 ) . As argued in proving that the operator G is bounded, one can obtain that
| y ( t ) | = | ν G y ( t ) | | ω 1 ( t ) | ϕ ( y ) Ω 1 + | ω 2 ( t ) | ψ ( y ) Ω 2 ,
which can be written as
y ω 1 ϕ ( y ) Ω 1 + ω 2 ψ ( y ) Ω 2 1 .
On the other hand, we can find a positive number M such that y M by assumption ( B 5 ). Let us set
U = { y X : y < M } .
Clearly, U contains a solution only when y = M . In other words, there is no solution y U such that y = ν G y for some ν ( 0 , 1 ) . Therefore, G has a fixed point y U ¯ which is a solution of the problem (1) and (2). The proof is finished. □

3.2. Uniqueness Result

Here we prove a uniqueness result for the problem (1) and (2) with the aid of Banach contraction mapping principle.
Theorem 3.
If the conditions ( B 1 ) and ( B 2 ) hold, then the problem (1) and (2) has a unique solution on [ 0 , 1 ] if
L Ω 1 + K Ω 2 < 1 ,
where Ω 1 and Ω 2 are defined by (9).
Proof. 
In the first step, we show that G B r B r , where B r = { y X : y r } with
r f 0 Ω 1 + h 0 Ω 2 1 ( L Ω 1 + K Ω 2 ) , f 0 = sup t [ 0 , 1 ] | f ( t , 0 ) | , h 0 = sup t [ 0 , 1 ] | h ( t , 0 ) | .
For y B r and using the condition ( B 1 ), we have
| f ( t , y ) | = | f ( t , y ) f ( t , 0 ) + f ( t , 0 ) | | f ( t , y ) f ( t , 0 ) | + | f ( t , 0 ) | L y + f 0 L r + f 0 .
Similarly, using ( B 2 ), we get
| h ( t , y ) | K r + h 0 .
In view of (13) and (14), we obtain
G y = sup t [ 0 , 1 ] | G y ( t ) | sup t [ 0 , 1 ] { 0 t ( t s ) β 1 Γ ( β ) I 1 α | f ( s , y ( s ) ) | + | λ | I 1 α + p I 0 + q | h ( s , y ( s ) ) | d s + | a 1 ( t ) | { | δ | 0 μ ( μ s ) β 1 Γ ( β ) I 1 α | f ( s , y ( s ) ) | + | λ | I 1 α + p I 0 + q | h ( s , y ( s ) ) | d s + 0 1 ( 1 s ) β 1 Γ ( β ) I 1 α | f ( s , y ( s ) ) | + | λ | I 1 α + p I 0 + q | h ( s , y ( s ) ) | d s } + | a 2 ( t ) | 0 ξ ( ξ s ) β 1 Γ ( β ) I 1 α | f ( s , y ( s ) ) | + | λ | I 1 α + p I 0 + q | h ( s , y ( s ) ) | d s } ( L r + f 0 ) sup t [ 0 , 1 ] { 0 t ( t s ) β 1 Γ ( β ) I 1 α d s + | a 1 ( t ) | [ | δ | 0 μ ( μ s ) β 1 Γ ( β ) I 1 α d s + 0 1 ( 1 s ) β 1 Γ ( β ) I 1 α d s ] + | a 2 ( t ) | 0 ξ ( ξ s ) β 1 Γ ( β ) I 1 α d s } + ( K r + h 0 ) | λ | sup t [ 0 , 1 ] { 0 t ( t s ) β 1 Γ ( β ) I 1 α + p I 0 + q d s + | a 1 ( t ) | | δ | 0 μ ( μ s ) β 1 Γ ( β ) I 1 α + p I 0 + q d s + 0 1 ( 1 s ) β 1 Γ ( β ) I 1 α + p I 0 + q d s + | a 2 ( t ) | 0 ξ ( ξ s ) β 1 Γ ( β ) I 1 α + p I 0 + q d s } ( L r + f 0 ) Γ ( α + 1 ) + ( K r + h 0 ) | λ | Γ ( α + p + 1 ) Γ ( q + 1 ) Δ = ( L r + f 0 ) Ω 1 + ( K r + h 0 ) Ω 2 < r ,
which implies that G y B r , for any y B r . Therefore, G B r B r . Next, we prove that G is a contraction. For that, let x , y X and t [ 0 , 1 ] . Then, by the conditions ( B 1 ) and ( B 2 ) , we obtain
G y G x = sup t [ 0 , 1 ] | ( G y ) ( t ) ( G x ) ( t ) | sup t [ 0 , 1 ] { 0 t ( t s ) β 1 Γ ( β ) I 1 α | f ( s , y ( s ) ) f ( s , x ( s ) ) | d s + | λ | 0 t ( t s ) β 1 Γ ( β ) I 1 α + p I 0 + q | h ( s , y ( s ) ) h ( s , x ( s ) ) | d s + | a 1 ( t ) | [ | δ | ( 0 μ ( μ s ) β 1 Γ ( β ) I 1 α | f ( s , y ( s ) ) f ( s , x ( s ) ) | d s + | λ | 0 μ ( μ s ) β 1 Γ ( β ) I 1 α + p I 0 + q | h ( s , y ( s ) ) h ( s , x ( s ) ) | d s ) + ( 0 1 ( 1 s ) β 1 Γ ( β ) I 1 α | f ( s , y ( s ) ) f ( s , x ( s ) ) | d s + | λ | 0 1 ( 1 s ) β 1 Γ ( β ) I 1 α + p I 0 + q | h ( s , y ( s ) ) h ( s , x ( s ) ) | d s ) ] + | a 2 ( t ) | [ 0 ξ ( ξ s ) β 1 Γ ( β ) I 1 α | f ( s , y ( s ) ) f ( s , x ( s ) ) | d s + | λ | 0 ξ ( ξ s ) β 1 Γ ( β ) I 1 α + p I 0 + q | h ( s , y ( s ) ) h ( s , x ( s ) ) | d r d u d s ] } L y x sup t [ 0 , 1 ] { 0 t ( t s ) β 1 Γ ( β ) I 1 α d s + | a 1 ( t ) | | δ | 0 μ ( μ s ) β 1 Γ ( β ) I 1 α d s + 0 1 ( 1 s ) β 1 Γ ( β ) I 1 α d s + | a 2 ( t ) | 0 ξ ( ξ s ) β 1 Γ ( β ) I 1 α d s } + K y x | λ | sup t [ 0 , 1 ] { 0 t ( t s ) β 1 Γ ( β ) I 1 α + p I 0 + q d s + | a 1 ( t ) | | δ | 0 μ ( μ s ) β 1 Γ ( β ) I 1 α + p I 0 + q d s + 0 1 ( 1 s ) β 1 Γ ( β ) I 1 α + p I 0 + q d s + | a 2 ( t ) | 0 ξ ( ξ s ) β 1 Γ ( β ) I 1 α + p I 0 + q d s } L Δ Γ ( α + 1 ) + K | λ | Δ Γ ( α + p + 1 ) Γ ( q + 1 ) y x = ( L Ω 1 + K Ω 2 ) y x .
From the above inequality, it follows by the assumption ( L Ω 1 + K Ω 2 ) < 1 that G is a contraction. Therefore, we deduce by Banach contraction mapping principle that there exists a unique fixed point for the operator G , which corresponds to a unique solution for the problem (1) and (2) on [ 0 , 1 ] . The proof is completed. □

3.3. Examples

In this subsection, we construct examples to illustrate the existence and uniqueness results obtained in the last two subsections. Let us consider the following problem:
D 1 3 / 2 D 0 + 1 / 2 y ( t ) + 2 I 1 4 / 3 I 0 + 5 / 4 h ( t , y ( t ) ) = f ( t , y ( t ) ) , t J : = [ 0 , 1 ] , y ( 0 ) = y ( 2 / 3 ) = 0 , y ( 1 ) = 1 2 y ( 3 / 4 ) .
Here α = 3 / 2 , β = 1 / 2 , λ = 2 , p = 4 / 3 , q = 5 / 4 , μ = 3 / 4 , δ = 1 / 2 , ξ = 2 / 3 , and
f ( t , y ) = 1 ( t 2 + 8 ) tan 1 y + e t , h ( t , y ) = 1 2 t 2 + 9 | y | 1 + | y | + e t .
Using the given data, it is found that L = 1 / 8 , K = 1 / 6 ,
a ¯ 1 = max t [ 0 , 1 ] | a 1 ( t ) | = | a 1 ( t ) | t = 1 1.121394517474712 , a ¯ 2 = max t [ 0 , 1 ] | a 2 ( t ) | = | a 2 ( t ) | t = t a 2 1.168623082364286 ,
where
t a 2 = β ( 1 δ μ β + 1 ) ( 1 δ μ β ) ( β + 1 ) 0.396975661732535 .
In consequence, we get
Ω 1 3.022797441671726 , Ω 2 1.451691300771574 , | Λ | 0.242702744426469 ,
where Ω 1 , Ω 2 are defined by (9) and Λ is given by (6).
(i) For illustrating Theorem 1, we have
| f ( t , y ) | σ ( t ) = e t + ( π / 2 ) t 2 + 8 , | h ( t , y ) | ρ ( t ) = e t + 1 2 t 2 + 9 ,
and that
L γ 1 + K γ 2 0.174044436618777 < 1 ,
where γ 1 0.848826363156775 and γ 2 0.407646847345084 . Clearly, the hypothesis of Theorem 1 is satisfied and consequently its conclusion applies to the problem (15).
(ii) In order to explain Theorem 2, we take the following values (instead of (16)) in the problem (15):
f ( t , y ) = 1 t 2 + 25 y cos y + π / 2 , h ( t , y ) = 1 5 t 2 + 4 sin y + 1 / 4 ,
and note that ω 1 ( t ) = 1 t 2 + 25 , ω 1 = 1 / 5 , ω 2 ( t ) = 1 5 t 2 + 4 , ω 2 = 1 / 10 , ϕ ( y ) = y + π / 2 and ψ ( y ) = y + 1 / 4 . By the condition ( B 5 ), we find that M > 3.939452045479877 . Thus, all the conditions of Theorem 2 are satisfied and, hence the problem (15) with f ( t , y ) and h ( t , y ) given by (17) has at least one solution on [ 0 , 1 ] .
(iii) It is easy to show that f ( t , y ) and h ( t , y ) satisfy the conditions ( B 1 ) and ( B 2 ) respectively with L = 1 / 8 and K = 1 / 6 and that L Ω 1 + K Ω 2 0.619798230337561 < 1 . Thus, all the assumptions of Theorem 3 hold true and hence the problem (15) has a unique solution on [ 0 , 1 ] .

4. Conclusions

We considered a fractional differential equation involving left Caputo and right Riemann–Liouville fractional derivatives of different orders and a pair of nonlinearities: I 1 p I 0 + q h ( t , y ( t ) ) = t 1 ( s t ) p 1 Γ ( p ) 0 s ( s v ) q 1 Γ ( q ) h ( v , y ( v ) ) d v d s (integral type) and f ( t , y ( t ) ) , equipped with four-point nonlocal boundary conditions. Different criteria ensuring the existence of solutions for the given problem are presented in Theorems 1 and 2, while the uniqueness of solutions is shown in Theorem 3. An interesting and scientific feature of the fractional differential Equation (1) is that the integral type of nonlinearity can describe composition of a physical quantity (like density) over two different arbitrary subsegments of the given domain. In the case of p = q = 1 , this composition takes the form t 1 0 s h ( v , y ( v ) ) d v d s . As pointed out in the introduction, fractional differential equations containing mixed (left Caputo and right Riemann–Liouville) fractional derivatives appear as Euler–Lagrange equations in the study of variational principles. So, such equations in the presence of the integral type of nonlinearity of the form introduced in (1) enhances the scope of Euler–Lagrange equations studied in [26]. Moreover, the fractional integro-differential Equation (1) can improve the description of the electromagnetic waves in dielectric media considered in [23]. As a special case, our results correspond to a three-point nonlocal mixed fractional order boundary value problem by letting δ = 0 , which is indeed new in the given configuration.

Author Contributions

Conceptualization, B.A.; Formal analysis, B.A., A.B., A.A. and S.K.N.; Funding acquisition, A.A.; Methodology, B.A., A.B., A.A. and S.K.N. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Deanship of Scientific Research (DSR), King Abdulaziz University, Jeddah, Saudi Arabia under grant number KEP-MSc-23-130-40.

Acknowledgments

This project was funded by the Deanship of Scientific Research (DSR), King Abdulaziz University, Jeddah, Saudi Arabia under grant no. (KEP-MSc-23-130-40). The authors, therefore, acknowledge with thanks DSR technical and financial support. We also thank the reviewers for their useful remarks on our work.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; North-Holland Mathematics Studies, 204; Elsevier Science B.V.: Amsterdam, The Netherlands, 2006. [Google Scholar]
  2. Henderson, J.; Luca, R.; Tudorache, A. On a system of fractional differential equations with coupled integral boundary conditions. Fract. Calc. Appl. Anal. 2015, 18, 361–386. [Google Scholar] [CrossRef]
  3. Peng, L.; Zhou, Y. Bifurcation from interval and positive solutions of the three-point boundary value problem for fractional differential equations. Appl. Math. Comput. 2015, 257, 458–466. [Google Scholar] [CrossRef]
  4. Ahmad, B.; Alsaedi, A.; Ntouyas, S.K.; Tariboon, J. Hadamard-Type Fractional Differential Equations, Inclusions and Inequalities; Springer: Cham, Switzerland, 2017. [Google Scholar]
  5. Ahmad, B.; Ntouyas, S.K. Nonlocal initial value problems for Hadamard-type fractional differential equations and inclusions. Rocky Mt. J. Math. 2018, 48, 1043–1068. [Google Scholar] [CrossRef]
  6. Cui, Y.; Ma, W.; Sun, Q.; Su, X. New uniqueness results for boundary value problem of fractional differential equation. Nonlinear Anal. Model. Control 2018, 23, 31–39. [Google Scholar] [CrossRef]
  7. Ahmad, B.; Luca, R. Existence of solutions for sequential fractional integro-differential equations and inclusions with nonlocal boundary conditions. Appl. Math. Comput. 2018, 339, 516–534. [Google Scholar] [CrossRef]
  8. Alsaedi, A.; Ahmad, B.; Alghanmi, M.; Ntouyas, S.K. On a generalized Langevin type nonlocal fractional integral multivalued problem. Mathematics 2019, 7, 1015. [Google Scholar] [CrossRef] [Green Version]
  9. Ahmad, B.; Alghamdi, N.; Alsaedi, A.; Ntouyas, S.K. A system of coupled multi-term fractional differential equations with three-point coupled boundary conditions. Fract. Calc. Appl. Anal. 2019, 22, 601–618. [Google Scholar] [CrossRef]
  10. Alsaedi, A.; Ahmad, B.; Alghanmi, M. Extremal solutions for generalized Caputo fractional differential equations with Steiltjes-type fractional integro-initial conditions. Appl. Math. Lett. 2019, 91, 113120. [Google Scholar] [CrossRef]
  11. Ahmad, B.; Alsaedi, A.; Alruwaily, Y.; Ntouyas, S.K. Nonlinear multi-term fractional differential equations with Riemann-Stieltjes integro-multipoint boundary conditions. AIMS Math. 2020, 5, 1446–1461. [Google Scholar] [CrossRef]
  12. Liang, S.; Wang, L.; Yin, G. Fractional differential equation approach for convex optimization with convergence rate analysis. Optim. Lett. 2020, 14, 145–155. [Google Scholar] [CrossRef]
  13. Atanackovic, T.M.; Stankovic, B. On a differential equation with left and right fractional derivatives. Fract. Calc. Appl. Anal. 2007, 10, 139–150. [Google Scholar] [CrossRef]
  14. Zhang, L.; Ahmad, B.; Wang, G. The existence of an extremal solution to a nonlinear system with the right-handed Riemann-Liouville fractional derivative. Appl. Math. Lett. 2014, 31, 1–6. [Google Scholar] [CrossRef]
  15. Khaldi, R.; Guezane-Lakoud, A. Higher order fractional boundary value problems for mixed type derivatives. J. Nonlinear Funct. Anal. 2017, 30, 9. [Google Scholar]
  16. Lakoud, A.G.; Khaldi, R.; Kilicman, A. Existence of solutions for a mixed fractional boundary value problem. Adv. Differ. Equ. 2017, 2017, 164. [Google Scholar] [CrossRef]
  17. Guezane-Lakoud, A.; Khaldi, R.; Torres, D.F.M. On a fractional oscillator equation with natural boundary conditions. Progr. Fract. Differ. Appl. 2017, 3, 191–197. [Google Scholar] [CrossRef] [Green Version]
  18. Ahmad, B.; Ntouyas, S.K.; Alsaedi, A. Existence theory for nonlocal boundary value problems involving mixed fractional derivatives. Nonlinear Anal. Model. Control 2019, 24, 937–957. [Google Scholar] [CrossRef]
  19. Laitinen, M.; Tiihonen, T. Heat transfer in conducting and radiating bodies. Appl. Math. Lett. 1997, 10, 5–8. [Google Scholar] [CrossRef] [Green Version]
  20. Laitinen, M.; Tiihonen, T. Integro-differential equation modelling heat transfer in conducting, radiating and semitransparent materials. Math. Methods Appl. Sci. 1998, 21, 375–392. [Google Scholar] [CrossRef]
  21. Hajmohammadi, M.R.; Nourazar, S.S.; Manesh, A.H. Semi-analytical treatments of conjugate heat transfer. J. Mech. Eng. Sci. 2013, 227, 492–503. [Google Scholar] [CrossRef]
  22. Yang, A.M.; Han, Y.; Zhang, Y.Z.; Wang, L.T.; Zhang, D.; Yang, X.J. On nonlocal fractional Volterra integro-differential equations in fractional steady heat transfer. Therm. Sci. 2016, 20, S789–S793. [Google Scholar] [CrossRef] [Green Version]
  23. Tarasov, V.E. Fractional integro-differential equations for electromagnetic waves in dielectric media. Theor. Math. Phys. 2009, 158, 355–359. [Google Scholar] [CrossRef] [Green Version]
  24. Krasnoselskii, M.A. Two remarks on the method of successive approximations. Uspekhi Mat. Nauk. 1955, 10, 123–127. [Google Scholar]
  25. Granas, A.; Dugundji, J. Fixed Point Theory; Springer: New York, NY, USA, 2005. [Google Scholar]
  26. Agrawal, O.P. Formulation of Euler-Lagrange equations for fractional variational problems. J. Math. Anal. Appl. 2002, 272, 368–379. [Google Scholar] [CrossRef] [Green Version]

Share and Cite

MDPI and ACS Style

Ahmad, B.; Broom, A.; Alsaedi, A.; Ntouyas, S.K. Nonlinear Integro-Differential Equations Involving Mixed Right and Left Fractional Derivatives and Integrals with Nonlocal Boundary Data. Mathematics 2020, 8, 336. https://doi.org/10.3390/math8030336

AMA Style

Ahmad B, Broom A, Alsaedi A, Ntouyas SK. Nonlinear Integro-Differential Equations Involving Mixed Right and Left Fractional Derivatives and Integrals with Nonlocal Boundary Data. Mathematics. 2020; 8(3):336. https://doi.org/10.3390/math8030336

Chicago/Turabian Style

Ahmad, Bashir, Abrar Broom, Ahmed Alsaedi, and Sotiris K. Ntouyas. 2020. "Nonlinear Integro-Differential Equations Involving Mixed Right and Left Fractional Derivatives and Integrals with Nonlocal Boundary Data" Mathematics 8, no. 3: 336. https://doi.org/10.3390/math8030336

APA Style

Ahmad, B., Broom, A., Alsaedi, A., & Ntouyas, S. K. (2020). Nonlinear Integro-Differential Equations Involving Mixed Right and Left Fractional Derivatives and Integrals with Nonlocal Boundary Data. Mathematics, 8(3), 336. https://doi.org/10.3390/math8030336

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