Next Article in Journal
Almost Boyd-Wong Type Contractions under Binary Relations with Applications to Boundary Value Problems
Next Article in Special Issue
Solution of Two-Dimensional Solute Transport Model for Heterogeneous Porous Medium Using Fractional Reduced Differential Transform Method
Previous Article in Journal
A Least Squares Estimator for Gradual Change-Point in Time Series with m-Asymptotically Almost Negatively Associated Errors
Previous Article in Special Issue
Two Novel Computational Techniques for Solving Nonlinear Time-Fractional Lax’s Korteweg-de Vries Equation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

New Results Achieved for Fractional Differential Equations with Riemann–Liouville Derivatives of Nonlinear Variable Order

by
Hallouz Abdelhamid
1,†,
Gani Stamov
2,*,†,
Mohammed Said Souid
3,† and
Ivanka Stamova
2,†
1
Laboratory of Informatics and Mathematics, Department of Mathematics, University of Tiaret, Tiaret 14035, Algeria
2
Department of Mathematics, University of Texas at San Antonio, San Antonio, TX 78249, USA
3
Department of Economic Sciences, University of Tiaret, Tiaret 14035, Algeria
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Axioms 2023, 12(9), 895; https://doi.org/10.3390/axioms12090895
Submission received: 12 August 2023 / Revised: 13 September 2023 / Accepted: 18 September 2023 / Published: 20 September 2023
(This article belongs to the Special Issue Fractional Calculus and the Applied Analysis)

Abstract

:
This paper proposes new existence and uniqueness results for an initial value problem (IVP) of fractional differential equations of nonlinear variable order. Riemann–Liouville-type fractional derivatives are considered in the problem. The new fundamental results achieved in this work are obtained by using the inequalities technique and the fixed point theory. In addition, uniform stability criteria for the solutions are derived. The accomplished results are new and complement the scientific research in the field. A numerical example is composed to show the efficacy and potency of the proposed criteria.

1. Introduction

Differential equations having fractional-order derivatives can be used as models of processes that generally possess boundless memory, which finds an advantage over integer-order differential equations. Indeed, it is known that the heredity properties of processes and materials are better represented by fractional-order derivatives. In addition, fractional-order derivatives and integrals provide more degrees of freedom to the model represented. As such, differential equations in terms of fractional derivatives have been extensively investigated during the last several decades as a technique for precisely describing real systems investigated in numerous fields of science, engineering and medicine [1,2,3,4,5,6,7].
Note that, historically, the first introduced derivatives and integrals in noninteger order were those of the Riemann–Liouville type [8,9]. It is also known that although a Caputo-type derivative of zero is zero, this is not true for the Riemann–Liouville-type fractional derivatives. The last fact leads to some complications in the investigation of fractional differential equations in terms of Riemann–Liouville-type derivatives. However, the excellent suggestions of the physical meaning of initial conditions for Riemann–Liouville fractional differential equations (see, for example, [10]) put these equations into a preferred modeling tool in various fields of science [2,5,11,12,13].
As an extension of the classical fractional differential systems, the type of differential systems with variable order of fractional derivatives has been also studied by numerous researchers using a variety of analytical and numerical methodologies. The fundamental idea that led to such an extension is to replace the constant β by a function β ( · ) . The basic notion of variable order derivatives has been first proposed by Samko and Ross in [14]. Since then, the investigations of variable order fractional operators in different forms and their applications have attracted more and more attention. For some excellent contributions to the theory of fractional differentiable systems with variable order derivatives we will refer to [15,16,17,18] and the corresponding references. The fractional variable order generalization makes the fractional differential systems in terms of variable order derivatives a more flexible apparatus in modeling various processes and natural phenomena. Hence, there has been an increasing research activity in the theory of such equations [19,20,21,22,23,24], including recently studied applications [25,26,27,28,29] which demonstrated the flexibility of this modeling approach. In fact, fractional differential equations in terms of variable order fractional derivatives have proven to be suitable in modeling numerous phenomena such as anomalous diffusion [28,30], tumor modeling [29], petroleum engineering [31], viscoelastic mechanics [32] and many others [17,25,26,27]. In addition, it has been shown [33] that variable order fractional calculus is a potential candidate to provide an accurate mathematical framework for efficient characterization of complex physical processes and systems.
Even though the substitution of a constant fractional order with a variable fractional order in the fractional derivatives seems simple, there are numerous difficulties in the mathematical investigations of differential systems in terms of variable order fractional derivatives. One of them is related to the absence of the semigroup property [19]. Another mathematical characterization of fractional variable order calculus which is problematic is the circumstance that a fractional derivative of variable order is not necessarily a left-inverse of the corresponding integral [21]. All these complications are reasons for an incomplete evolution of the theory of variable-order fractional differential systems.
Recently, Souid et al. have intensively contributed to this field [34,35,36,37,38,39,40,41,42]. A series of papers is concerned with the questions of the existence and uniqueness of solutions of different classes of differential systems with fractional derivatives of variable order [34,36,37,38,39,40,41,42]. Some of our research papers are devoted to the qualitative analysis of such problems. See, for example, ref. [35] and some of the references therein. In our studies, we apply different techniques such as fixed point theorem, a measure of noncompactness, upper–lower solution methods, piecewise constant functions and some others. The variety of problems investigated includes fractional differential systems with variable order derivatives of Caputo type, Hadamard type and Riemann–Liouville type as well as multiterm fractional boundary value problems of variable order.
In this paper, motivated by the above related works in this regard, we investigate the the following initial value problem (IVP) for fractional differential equations (FDE) with nonlinear variable order (NVO) derivatives of Riemann–Liouville type, defined as follows:
D 0 + β ( t , y ( t ) ) y ( t ) = ψ ( t , y ( t ) ) , t Δ : = ( 0 , M ] , 0 < M < , ( A ) ( IVPFDENVO ) y ( 0 ) = 0 , ( B )
where D 0 + β ( t , y ( t ) ) stands for the Riemann–Liouville fractional derivative of the variable order β ( t , y ( t ) ) , ψ is a given function and β satisfies 0 < β * β ( t , y ( t ) ) β * < 1 .
Note that, although differential equations with nonlinear variable-order fractional derivatives are already applied as models of numerous problems investigated in science and engineering, the results related to such equations are limited. Hence, the development of the area requires the establishment of new existence and stability results.
The main contributions of our paper are stated as follows:
(1) an IVPFDENVO is defined, which extends some existing problems for fractional-order systems with Riemann–Liouville fractional derivatives of variable order investigated in the literature;
(2) new inequalities are proved for fractional integrals and derivatives of nonlinear variable order;
(3) new criteria for the existence and uniqueness of the solutions to the introduced problem are proposed;
(4) we consider two different Banach spaces of functions to which we apply the inequalities technique and fixed points theorems;
(5) novel uniform stability results are established via an inequalities technique.
The body of the manuscript is organized in the following manner. In Section 2 we present notations, definitions and lemmas that will be necessary to carry out our study. Section 3 is devoted to new existence results for solutions of the IVPFDENVO. A Banach-type fixed point theorem is applied to two different spaces of functions. In Section 4, uniform stability criteria are derived. Section 5 represents numerical applications and simulations. Finally, conclusion remarks are presented in Section 6.

2. Preliminaries

Some preliminary results will be presented in this section together with the related notations and definitions.
We consider the Banach space of all real-valued continuous functions x : Δ R , R = ( , ) , with the norm
x = sup | x ( t ) | : t Δ
which we will denoted by C ( Δ , R ) .
The Banach space of all functions x : Δ R such that
0 < γ < 1 and t γ x ( t ) C ( Δ , R )
with the norm
x γ = sup { t γ | x ( t ) | : t Δ }
will be denoted by the symbol C γ ( Δ , R ) .
The symbol L p ( Δ , R ) represents the Banach space of all functions x : Δ R which are Lebesgue measurable such that
p 1 and 0 M | x ( s ) | p d s <
with the norm
x p = : 0 M | x ( s ) | p d s 1 p .
Remark 1.
The following observations are made to make our study easy in the sequel:
(1)
if  0 < M 1  then  M β ( s , y ( s ) ) 1 M β * 1 .
(2)
if  1 < M  then  M β ( s , y ( s ) ) 1 1 .
Set
Λ * = max { 1 , M β * 1 } .
(3)
The function Γ ( β ( t , f ( t ) ) ) is continuous as a composition of two continuous functions, hence we can set:
M f = max t [ 0 , M ] 1 Γ ( β ( t , f ( t ) ) .
Let us consider two continuous functions β : Δ × R ( 0 , β * ] and α : Δ × R ( 0 , α * ] . Let f ( t ) C ( Δ , R ) .
Definition 1
([14,16,18]). The left Riemann–Liouville type integral of fractional variable order β, β = β ( t , f ( t ) ) for f ( t ) is
I a 1 + β ( t , f ( t ) ) f ( t ) = a 1 t ( t s ) β ( s , f ( s ) ) 1 Γ ( β ( s , f ( s ) ) ) f ( s ) d s , t > a 1 ,
where Γ ( · ) denotes the Gamma function and a 1 Δ .
If β ( t , y ( t ) ) is a constant β, then (3) will be reduced to the standard Riemann–Liouville fractional integral of a constant fractional order given by [2,3,5]
I a 1 + β f ( t ) = 1 Γ ( β ) β * t ( t s ) β 1 f ( s ) d s , t > a 1 .
Definition 2
([14,16,18]). The left Riemann–Liouville type derivative of fractional variable order α, α = α ( t , f ( t ) ) for a function f ( t ) is given by
D a 1 + α ( t , f ( t ) ) f ( t ) = d d t I a 1 + 1 α ( t , f ( t ) ) f ( t ) = d d t a 1 t ( t s ) α ( s , f ( s ) ) Γ ( 1 α ( s , f ( s ) ) ) f ( s ) d s , t > a 1 .
If α ( t , y ( t ) ) is a constant α, then (5) will become
D a 1 + α f ( t ) = d d t I a 1 + 1 α f ( t ) = d d t a 1 t ( t s ) α Γ ( 1 α ) f ( s ) d s , t > a 1 ,
which represents the Riemann–Liouville derivative of a constant fractional order α [2,3,5].
For more characteristics of the integrals and derivatives of fractional constant orders we refer to [2,3,5], and about the integrals and derivatives of fractional variable orders, see [14,16,19].
Remark 2.
It is observed that the semigroup property is not always fulfilled for general functions β ( t , y ( t ) ) , α ( t , y ( t ) ) ; i.e., it is possible
I a 1 + β ( t , f ( t ) ) I a 1 + α ( t , f ( t ) ) f ( t ) I a 1 + β ( t , f ( t ) ) + α ( t , f ( t ) ) f ( t ) .
For more details, see [23,38,43].
We will present some specifications of the fractional integrals of the Riemann–Liouville type of constant orders which will be used in the coming lemmas.
Lemma 1
([3]). If γ R , then the Riemann–Liouville type integral of fractional constant order is bounded in C γ ( Δ , R ) and we have for f C γ ( Δ , R )
I 0 + α f γ M α Γ ( 1 γ ) Γ ( 1 + α γ ) f γ , α > 0
Lemma 2
([3]). If γ R , then the Riemann–Liouville type fractional integral is bounded in L p ( Δ , R ) and we have for f L p ( Δ , R )
I 0 + α f p M α α Γ ( α ) f p , α > 0 .
On the base of Lemmas 1 and 2, we will prove similar inequalities for fractional integrals of nonlinear variable order of Riemann–Liouville type.
Lemma 3.
If β : Δ × R ( 0 , 1 ] is a continuous function, such that 0 < β * β ( t , f ( t ) ) β * < 1 , then I 0 + β ( t , y ( t ) ) f ( t ) C γ ( Δ , R ) for f C γ ( Δ , R ) . Moreover, we have:
(i)
I 0 + β ( t , f ( t ) ) f γ M Γ ( 1 γ ) Γ ( β * ) M f Λ * Γ ( 1 + β * γ ) f γ .
(ii) For f , g C γ ( Δ , R ) , we have
I 0 + β ( t , f ( t ) ) f I 0 + β ( t , g ( t ) ) g γ 4 M B Γ ( β * ) Λ * Γ ( 1 γ ) Γ ( 1 + β * γ ) f g γ ,
where B = max M f , M g , Λ * is defined by Equation (1) and M f and M g are defined by Equation (2).
Proof. 
(i) Let f C γ ( Δ , R ) . From Equation (3) we have
| I 0 + β ( t , f ( t ) ) f ( t ) | M f 0 t ( t s ) β ( s , f ( s ) ) 1 | f ( s ) | d s M f 0 t M β ( s , f ( s ) ) 1 t s M β ( s , f ( s ) ) 1 | f ( s ) | d s M f Λ * M β * 1 0 t ( t s ) β * 1 | f ( s ) | d s Γ ( β * ) M f Λ * M β * 1 I 0 + β * | f ( t ) | .
The above estimate implies
I 0 + β ( t , f ( t ) ) f γ Γ ( β * ) M f Λ * M β * 1 I 0 + β * | f | γ .
We apply (7) to obtain
I 0 + β ( t , f ( t ) ) f γ M Γ ( 1 γ ) Γ ( β * ) M f Λ * Γ ( 1 + β * γ ) f γ .
(ii) For f , g C γ ( Δ , R ) , we have
I 0 + β ( t , f ( t ) ) f ( t ) I 0 + β ( t , g ( t ) ) g ( t ) = 0 t ( t s ) β ( s , f ( s ) ) 1 Γ ( β ( s , f ( s ) ) ) f ( s ) d s ( t s ) β ( s , g ( s ) ) 1 Γ ( β ( s , g ( s ) ) ) g ( s ) d s 2 B Λ * 0 t t s M β ( s , f ( s ) ) 1 + t s M β ( s , g ( s ) ) 1 × | f ( s ) g ( s ) | d s 4 B Λ * M β * 1 0 t ( t s ) β * 1 | f ( s ) g ( s ) | d s 4 B Γ ( β * ) Λ * M β * 1 I 0 + β * | f ( s ) g ( s ) | ,
which, after the application of (7), implies
I 0 + β ( t , f ( t ) ) f I 0 + β ( t , g ( t ) ) g γ 4 M B Γ ( β * ) Λ * Γ ( 1 γ ) Γ ( 1 + β * γ ) f g γ .
Lemma 4.
If β : Δ × R ( 0 , 1 ] is a continuous function, such that 0 < β * β ( t , f ( t ) ) β * < 1 , then I 0 + β ( t , f ( t ) ) f ( t ) L p ( Δ , R ) for f L p ( Δ , R ) . Moreover, we have:
(i)
I 0 + β ( t , f ( t ) ) f p M M f Λ * β * f p .
(ii) For f , g L p ( Δ , R ) we have
I 0 + β ( t , f ( t ) ) f I 0 + β ( t , g ( t ) ) g p 4 M B Λ * β * f g p .
Proof. 
(i) Using (9), we obtain
I 0 + β ( t , f ( t ) ) f p Γ ( β * ) M f Λ * M β * 1 I 0 + β * | f | p .
Now, we apply Equation (8) to obtain
I 0 + β ( t , f ( t ) ) f p M M f Λ * β * f p .
(ii) We combine Equation (12) with the Hölder’s inequality to obtain
| I 0 + β ( t , f ( t ) ) f ( t ) I 0 + β ( t , g ( t ) ) g ( t ) | p 4 B Γ ( β * ) Λ * M β * 1 I 0 + β * | f ( s ) g ( s ) | p 4 B Γ ( β * ) Λ * M β * 1 p I 0 + β * | f ( s ) g ( s ) | p .
Integrating both sides of (15) on [ 0 , M ] and take 1 p -root on both sides, we obtain
I 0 + β ( t , f ( t ) ) f I 0 + β ( t , g ( t ) ) g p 4 M B Λ * β * f g p .
The proof of (12) is completed. □
The following lemma will also be useful.
Lemma 5
([44]). Let Θ be a nonempty, bounded Banach space and F : Θ Θ be a mapping such that for some n N , F n is a contraction, where F n = F F F n times. Then F has a unique fixed point in Θ.

3. Achieved Existence Results

Definition 3.
A function y C γ ( Δ , R ) or y L p ( Δ , R ) is said to be a solution for (IVPFDENVO) if and only if it verifies (IVPFDENVO(A)) and (IVPFDENVO(B)), simultaneously.
In order to present our new existence results in the Banach spaces C γ ( Δ , R ) and L p ( Δ , R ) , we will analyze an equivalent integral form of the IVPFDENVO(A).
Lemma 6.
Let y be an element of C γ ( Δ , R ) or L p ( Δ , R ) . Then, equation (IVPFDENVO(A)) is equivalent to
I 0 + 1 β ( t , y ( t ) ) y ( t ) = 0 t ( t s ) β ( s , y ( s ) ) Γ ( 1 β ( s , y ( s ) ) ) y ( s ) d s = 0 t ψ ( s , y ( s ) ) d s , t Δ .
Proof. 
Let y C γ ( Δ , R ) or y L p ( Δ , R ) . Then, for equation (IVPFDENVO(A)) we have the following representation
D 0 + β ( t , y ( t ) ) y ( t ) = d d t 0 t ( t s ) β ( s , y ( s ) ) Γ ( 1 β ( s , y ( s ) ) ) y ( s ) d s = ψ ( t , y ( t ) ) .
Then, both sides of (17) can be integrated from [ 0 , t ] , to obtain
0 t ( t s ) β ( s , y ( s ) ) Γ ( 1 β ( s , y ( s ) ) ) y ( s ) d s = c 0 + 0 t ψ ( s , y ( s ) ) d s .
Evaluating (18) at t = 0 gives us c 0 = 0 .
Conversely, differentiating both sides of (16) to reach
d d t 0 t ( t s ) β ( s , y ( s ) ) Γ ( 1 β ( s , y ( s ) ) ) y ( s ) d s = ψ ( t , y ( t ) ) ,
from which we obtain (IVPFDENVO(A)).
The proof is concluded. □
The following assumptions will be essential in our analysis.
(A1)
The function β : Δ × R ( 0 , β * ] is continuous on its domain.
(A2)
The function ψ : Δ × R R is a continuous function with respect to its first variable and satisfies:
ψ ( s , χ 1 ) ψ ( s , χ 2 ) k | χ 1 χ 2 | , χ 1 , χ 2 R
for t Δ and k > 0 .

3.1. Existence Result in C γ ( Δ , R )

The first obtained result is based on Lemma 5.
Theorem 1.
Assume that (A1) and (A2) are satisfied. Then the (IVPFDENVO) has a unique solution in C 1 β * ( Δ , R ) .
Proof. 
Let us consider γ = 1 β * and the set of elements Θ in the space C γ ( Δ , R ) such that y ( 0 ) = 0 . Define the following operator
Π : Θ Θ ,
where
( Π y ) ( t ) = y ( t ) + I 0 + 1 β ( t , y ( t ) ) y ( t ) 0 t ψ ( s , y ( s ) ) d s .
First, for two x , y : Δ R using (A2), we have
0 t ψ ( s , x ( s ) ) ψ ( s , y ( s ) ) d s k t γ | | y x | | γ .
Then, from (18) we can obtain the following estimation
| ( Π x ) ( t ) ( Π y ) ( t ) | | y ( t ) x ( t ) | + | I 0 + 1 β ( t , x ( t ) ) x ( t ) I 0 + 1 β ( t , y ( t ) ) y ( t ) | + 0 t ψ ( s , x ( s ) ) ψ ( s , y ( s ) ) d s t γ | | y x | | γ + | I 0 + 1 β ( t , x ( t ) ) x ( t ) I 0 + 1 β ( t , y ( t ) ) y ( t ) | + k t γ | | y x | | γ .
We multiply both sides of Equation (21) with t γ and take the sup of both sides to obtain
| | Π x Π y | | γ | | y x | | γ + | | I 0 + 1 β ( t , x ( t ) x I 0 + 1 β ( t , y ( t ) y | | γ + k | | y x | | γ .
Using Equation (10) we obtain
| | Π x Π y | | γ | | y x | | γ + 4 M B Γ ( β * ) 2 Λ * Γ ( 2 β * ) y x γ + k | | y x | | γ .
If we set ζ = 1 + 4 M B Γ ( β * ) 2 Λ * Γ ( 2 β * ) + k , then we have
| | Π x Π y | | γ ζ | | y x | | γ .
By induction, it is trivial to prove that
| | Π n x Π n y | | γ ζ n n ! | | y x | | γ ,
where Π n = Π Π Π n times.
Since ζ n n ! is the general term of the convergent exponential series e ζ , it approaches zero as n approaches infinity, and so for n sufficiently large we have
ζ n n ! < 1 .
Lemma 5 asserts that the operator Π has a unique fixed point in Θ .
This implies that
I 0 + 1 β ( t , y ( t ) ) y ( t ) = 0 t ψ ( s , y ( s ) ) d s
with y ( 0 ) = 0
Finally, from Lemma 6, we obtain
D 0 + β ( t , y ( t ) ) y ( t ) = ψ ( t , y ( t ) ) w i t h y ( 0 ) = 0 .
This concludes our proof. □

3.2. Existence Result in L p ( Δ , R )

Theorem 2.
Under the assumptions (A1) and (A2), the (IVPFDENVO) has a unique solution in the Banach space L p ( Δ , R ) .
Proof. 
We consider the set Θ as an element in L p ( Δ , R ) such that y ( 0 ) = 0 , and the operator
Π : Θ Θ ,
where
( Π y ) ( t ) = y ( t ) + I 0 + 1 β ( t , y ( t ) ) y ( t ) 0 t ψ ( s , y ( s ) ) d s .
Then, we have from (A2) that for x , y : Δ R ,
0 t ψ ( s , x ( s ) ) ψ ( s , y ( s ) ) d s k M 1 / p | | y x | | p
and
| ( Π x ) ( t ) ( Π y ) ( t ) | p 2 p | y ( t ) x ( t ) | p + | I 0 + 1 β ( t , x ( t ) ) x ( t ) I 0 + 1 β ( t , y ( t ) ) y ( t ) ) | p + 0 t ψ ( s , x ( s ) ) ψ ( s , y ( s ) ) d s p 2 p | y ( t ) x ( t ) | p + | I 0 + 1 β ( t , x ( t ) ) x ( t ) I 0 + 1 β ( t , y ( t ) ) y ( t ) | p + k p M | | y x | | p p .
Integrating Equation (24) on [ 0 , M ] , we obtain
| | Π x Π y | | p p 2 p | | y x | | p p + | | I 0 + 1 β ( t , x ( t ) ) x I 0 + 1 β ( t , y ( t ) ) y | | p p + k p M 2 | | y x | | p p .
Using Equation (14) from Lemma 4, we obtain
| | Π x Π y | | p p 2 p ( | | y x | | p p + 4 M B Λ * 1 β * p y x p p + k p M 2 | | y x | | p p ) 2 p 1 + 4 M B Λ * 1 β * p + k p M 2 | | y x | | p p .
If we denote ζ = 2 1 + 4 M B Λ * 1 β * p + k p M 2 1 p , then the assertion of the theorem can be proved analogously to the final part of the proof of Theorem 1.
This concludes our proof. □
Remark 3.
Since the operators involved in the description of fractional differential equations of variable nonlinear fractional order have complex properties, the research results in this direction are still limited [23,34,36,37,38,39,40,41,42,43]. With the proposed new criteria in Theorems 1 and 2, we complement and extend the existence of theoretical results for such initial value problems. The delivered results are obtained by using fixed point theory and are presented in two different Banach spaces.
Remark 4.
The criteria presented in this section are also extensions and generalizations of some announced results that are considered initial value problems for fractional constant-order differential equations to the variable order case [45,46,47,48]. In fact, the consideration of nonlinear variable orders leads to the definition of more complex and generalized problems that can be used in the applications.
Remark 5.
Different from the existing results for differential systems with fractional derivatives of variable order, in this study, we consider fractional derivatives of Riemann–Liouville types of order β : Δ × R ( 0 , 1 ] . Instead of the approaches introduced in [23,43], such as piecewise continuous functions and the Picard scheme, we apply the operator approach and a Banach-type fixed point theorem, which we consider as more appropriate for the considered problem from the applied perspective. The proposed strategy can be applied to similar problems considering delays and impulsive factors, which can motivate future research.

4. Uniform Stability

In this Section, the newly achieved existence and uniqueness results will be applied to derive uniform stability criteria for the solution of the (IVPFDENVO).
Consider a solution x ( t ) of the equation IVPFDENVO(A) corresponding to an initial condition y ( 0 ) = x 0 . Let x ¯ ( t ) be another solution of equation IVPFDENVO(A) corresponding to an initial condition x ¯ ( 0 ) = x ¯ 0 .
Definition 4
([49]). The solution x ( t ) of the IVPFDENVO(A) is uniformly stable if for any ϵ > 0 , there exists υ ( ϵ ) > 0 such that | x 0 x ¯ 0 | υ ( ϵ ) implies
x x ˜ ϵ .
Note that we will apply Definition 4 for the (IVPFDENVO), i.e., we will consider the unique solution y ( t ) which satisfies the zero initial condition IVPFDENVO(B). Also, in the proof, we will use functions from the class C 1 β * ( Δ , R ) .
Theorem 3.
Assume that (A1)–(A2) are satisfied and, in addition,
4 M B Γ ( β * ) 2 Λ * Γ ( 2 β * ) + k < 1 .
Then, the unique solution of IVPFDENVO is uniformly stable.
Proof. 
Theorem 1 guarantees that the IVPFDENVO has a unique solution y ( t ) . Consider a solution y ¯ ( t ) of the IVPFDENVO(A) which corresponds to an initial condition y ¯ ( 0 ) = y ¯ 0 .
For the solution y ( t ) , we have from Lemma 6,
I 0 + 1 β ( t , y ( t ) ) y ( t ) = 0 t ψ ( s , y ( s ) ) d s , t Δ .
Similarly, for y ¯ ( t ) , from (16) we obtain
I 0 + 1 β ( t , y ¯ ( t ) ) y ¯ ( t ) = y ¯ 0 + 0 t ψ ( s , y ¯ ( s ) ) d s , t Δ .
Hence,
| y ( t ) y ¯ ( t ) | | y ¯ 0 | + I 0 + 1 β ( t , y ( t ) ) y I 0 + 1 β ( t , y ¯ ( t ) ) y ¯ + 0 t ψ ( s , y ( s ) ) ψ ( s , y ¯ ( s ) ) d s | y ¯ 0 | + I 0 + 1 β ( t , y ( t ) ) y I 0 + 1 β ( t , y ¯ ( t ) ) y ¯ + k t γ | | y y ¯ | | γ .
Using Equation (10) for γ = 1 β * , we obtain
| | y y ¯ | | γ M 1 β * | y ¯ 0 | + 4 M B Γ ( β * ) 2 Λ * Γ ( 2 β * ) y y ¯ γ + k | | y y ¯ | | γ .
Thus
| | y y ¯ | | γ M 1 β * 1 4 M B Γ ( β * ) 2 Λ * Γ ( 2 β * ) k 1 | y ¯ 0 | .
Therefore, if | y ¯ 0 | = | y ( 0 ) y ¯ 0 | < υ ( ϵ ) , then | | y y ¯ | | < ϵ , which completes the proof of the theorem. □
Remark 6.
Theorem 3 presents criteria for uniform stability of in C 1 β * ( Δ , R ) . A similar uniform stability result can be proved using the norm in L p ( Δ , R ) .
Remark 7.
Stability results for fractional differential systems with derivatives of variable order have been considered in the existent literature [20,35,38,40]. Hence, the proposed new stability criteria are a contribution to the development of the stability theory of such equations. Different from all existing studies which mainly considered Ulam–Hyers stability, we establish a uniform stability result. Also, the obtained qualitative result shows the applicability of the derived fundamental results in the previous section.

5. Approximate Numerical Applications

Example 1.
Let us consider the following fractional initial value problem
D 0 + β ( t , y ( t ) ) y ( t ) = ψ ( t , y ( t ) ) , t Δ : = [ 0 , 1 ] , ( A ) ( IVPNFDEVO ) y ( 0 ) = 0 , ( B )
with
β ( t , y ) = 1 4 t + 1 4 ( 1 + 2 y 3 )
and
ψ ( t , y ( t ) ) = e t y ( t ) ( 4 + e 2 t ) ( 1 + y ( t ) ) .
We have that β is a continuous function on Δ × R and 0 < β ( t , s ) < 1 .
Also,
| ψ ( t , x ) ψ ( t , y ) | e t ( 4 + e 2 t ) x 1 + x y 1 + y e t | x y | ( 4 + e 2 t ) ( 1 + x ) ( 1 + y ) e t ( 4 + e 2 t ) | x y | 1 5 | x y | ,
It is easy to check that for the given choice of nonlinear functions β and ψ Assumptions (A1) and (A2) are satisfied. Therefore, by Theorems 1 and 2, the problem (IVPNFDEVO) has a unique solution.
Example 2.
Let us consider the following fractional initial value problem
D 0 + β ( t , y ( t ) ) y ( t ) = ψ ( t , y ( t ) ) , t Δ : = [ 0 , 1 ] , ( A ) ( IVPNFDEVO ) y ( 0 ) = 0 , ( B )
with
β ( t , y ) = t 3 3 + t 3 + 1 3 ( y 2 + 1 )
and
ψ ( t , y ( t ) ) = y ( t ) 4 + 1 6 + t , t Δ .
We have that β is a continuous function on Δ × R and 0 < β ( t , y ) < 1 .
Also,
| ψ ( t , x ) ψ ( t , y ) | x 4 + 1 6 + t y 4 + 1 6 + t x 4 y 4 1 4 | x y | , x , y R .
Hence, for the given choice of nonlinear functions β and ψ Assumptions (A1) and (A2) are satisfied. Therefore, by Theorems 1 and 2, the problem (IVPNFDEVO)) has a unique solution.
In the remaining part, some numerical applications are demonstrated.
The approximate solution y ( t ) for β ( t , y ) = t 3 3 + t 3 + 1 3 ( y 2 + 1 ) with t [ 0 , 1 ] is represented in Figure 1.
On the other hand, Figure 2 and Figure 3 present the graphs of the functions β ( t , y ) and ψ ( t , y ) .
In Table 1, we present our β ( t , y ) and y ( t ) with different value of t [ 0 , 1 ] .

6. Conclusions

In this study, we introduce an initial value problem for a class of nonlinear differential equations with fractional Riemann–Liouville-type derivative of variable nonlinear order. The existence and uniqueness of the solution are investigated in two different Banach spaces and new criteria are achieved. A Banach-type fixed point theorem is applied as a proof technique. In addition, a uniform stability result is established for the solution of the investigated problem which shows the efficiency of the existence criteria. The newly achieved outcomes complement results for different classes of variable-order fractional differential equations and provide an extension of the theory of such problems. Numerical applications are also elaborated. The introduced problem and the achieved results can be developed. Some directions of the future expansion of the topic include considering delay terms, reaction-diffusion terms and impulsive effects.

Author Contributions

Conceptualization, M.S.S. and I.S.; methodology, H.A., G.S., M.S.S. and I.S.; formal analysis, H.A., G.S., M.S.S. and I.S.; investigation, H.A., G.S., M.S.S. and I.S.; writing—original draft preparation, I.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Baleanu, D.; Diethelm, K.; Scalas, E.; Trujillo, J.J. Fractional Calculus: Models and Numerical Methods, 1st ed.; World Scientific: Singapore, 2012; ISBN 978-981-4355-20-9. [Google Scholar]
  2. Magin, R. Fractional Calculus in Bioengineering, 1st ed.; Begell House: Redding, CA, USA, 2006; ISBN 978-1567002157. [Google Scholar]
  3. Kilbas, A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations, 1st ed.; Elsevier: New York, NY, USA, 2006; ISBN 9780444518323. [Google Scholar]
  4. Petráš, I. Fractional-Order Nonlinear Systems, 1st ed.; Springer: Heidelberg, Germany; Dordrecht, The Netherlands; London, UK; New York, NY, USA, 2011; ISBN 978-3-642-18101-6. [Google Scholar]
  5. Samko, S.G.; Kilbas, A.A.; Marichev, O.I. Fractional Integrals and Derivatives: Theory and Applications, 1st ed.; Gordon and Breach: Yverdon, Switzerland, 1993; ISBN 9782881248641. [Google Scholar]
  6. Stamova, I.M.; Stamov, G.T. Functional and Impulsive Differential Equations of Fractional Order: Qualitative Analysis and Applications, 1st ed.; Taylor & Francis Group: Boca Raton, FL, USA, 2017; ISBN 9781498764834. [Google Scholar]
  7. Tarasov, V.E. Fractional Dynamics: Application of Fractional Calculus to Dynamics of Particles, Fields and Media, 1st ed.; Springer: Beijing, China, 2015; ISBN 978-3-642-14003-7. [Google Scholar]
  8. Liouville, J. Mémoire sur quelques questions de géométrie et de mécanique, et sur un nouveau genre de calcul pour résoudre ces questions. J. Ec. Polytech. 1832, 13, 1–69. [Google Scholar]
  9. Liouville, J. Mémoire sur le calcul des différentielles à indices quelconques. J. Ec. Polytech. 1832, 13, 71–162. [Google Scholar]
  10. Heymans, N.; Podlubny, I. Physical interpretation of initial conditions for fractional differential equations with Riemann–Liouville fractional derivatives. Rheol. Acta 2006, 45, 765–771. [Google Scholar] [CrossRef]
  11. Heymans, N.; Kitagawa, M. Modelling “unusual” behaviour after strain reversal with hierarchical fractional models. Rheol. Acta 2004, 43, 383–389. [Google Scholar] [CrossRef]
  12. Koeller, R.C. Applications of fractional calculus to the theory of viscoelasticity. J. Appl. Mech. 1984, 51, 299–307. [Google Scholar] [CrossRef]
  13. Schiessel, H.; Blumen, A. Hierarchical analogues to fractional relaxation equations. J. Phys. A 1993, 26, 5057–5069. [Google Scholar] [CrossRef]
  14. Samko, S.; Ross, B. Integration and differentiation to a variable fractional order. Integral Transform. Spec. Funct. 1993, 1, 277–300. [Google Scholar] [CrossRef]
  15. Lorenzo, C.F.; Hartley, T.T. Initialization, conceptualization, and application in the generalized fractional calculus. Crit. Rev. Biomed. Eng. 2007, 35, 477–553. [Google Scholar] [CrossRef]
  16. Samko, S. Fractional integration and differentiation of variable order: An overview. Nonlinear Dyn. 2013, 71, 653–662. [Google Scholar] [CrossRef]
  17. Sun, H.; Chang, A.; Zhang, Y.; Chen, W. A review on variable-order fractional differential equations: Mathematical foundations, physical models, numerical methods and applications. J. Fract. Calc. Appl. 2019, 22, 27–59. [Google Scholar] [CrossRef]
  18. Valerio, D.; Costa, J.S. Variable-order fractional derivatives and their numerical approximations. Signal Process. 2011, 91, 470–483. [Google Scholar] [CrossRef]
  19. Almeida, R.; Tavares, D.; Torres, D.F.M. The Variable-Order Fractional Calculus of Variations, 1st ed.; Springer: Cham, Switzerland, 2019; ISBN 978-3-319-94005-2. [Google Scholar]
  20. Benkerrouche, A.; Souid, M.S.; Stamov, G.; Stamova, I. On the solutions of a quadratic integral equation of the Urysohn type of fractional variable order. Entropy 2022, 24, 886. [Google Scholar] [CrossRef]
  21. Garrappa, R.; Giusti, A.; Mainardi, F. Variable-order fractional calculus: A change of perspective. Commun. Nonlinear Sci. Numer. Simul. 2021, 102, 105904. [Google Scholar] [CrossRef]
  22. Odzijewicz, T.; Malinowska, A.B.; Torres, D.F.M. Fractional variational calculus of variable order. In Advances in Harmonic Analysis and Operator Theory, 1st ed.; Almeida, A., Castro, L., Speck, F.O., Eds.; Birkhäuser: Basel, Switzerland, 2013; Volume 229, pp. 291–301. [Google Scholar]
  23. Zhang, S.; Hu, L. The existence of solutions and generalized Lyapunov-type inequalities to boundary value problems of differential equations of variable order. AIMS Math. 2020, 5, 2923–2943. [Google Scholar] [CrossRef]
  24. Zhang, S.; Sun, S.; Hu, L. Approximate solutions to initial value problem for differential equation of variable order. J. Fract. Calc. Appl. 2018, 9, 93–112. [Google Scholar]
  25. Ali, U.; Ahmad, H.; Abu-Zindah, H. Soliton solutions for nonlinear variable-order fractional Korteweg–de Vries (KdV) equation arising in shallow water waves. J. Ocean Eng. Sci. 2022. [Google Scholar] [CrossRef]
  26. Patnaik, S.; Hollkamp, J.P.; Semperlotti, F. Applications of variable-order fractional operators: A review. Proc. R. Soc. A 2020, 476, 20190498. [Google Scholar] [CrossRef]
  27. Lu, X.; Li, H.; Chen, N. An indicator for the electrode aging of lithium-ion batteries using a fractional variable order model. Electrochim. Acta 2019, 299, 378–387. [Google Scholar] [CrossRef]
  28. Sun, H.G.; Chen, W.; Chen, Y.Q. Variable-order fractional differential operators in anomalous diffusion modeling. Physica A 2009, 388, 4586–4592. [Google Scholar] [CrossRef]
  29. Sweilam, N.H.; AL-Mekhlafi, S.M.; Alshomrani, A.S.; Baleanu, D. Comparative study for optimal control nonlinear variable-order fractional tumor model. Chaos Solitons Fract. 2020, 136, 1–10. [Google Scholar] [CrossRef]
  30. Zhuang, P.; Liu, F.; Anh, V.; Turner, I. Numerical methods for the variable-order fractional advection-diffusion equation with a nonlinear source term. SIAM J. Numer. Anal. 2009, 47, 1760–1781. [Google Scholar] [CrossRef]
  31. Obembe, A.D.; Hossain, M.E.; Abu-Khamsin, S.A. Variable-order derivative time fractional diffusion model for heterogeneous porous media. J. Petrol. Sci. Eng. 2017, 152, 391–405. [Google Scholar] [CrossRef]
  32. Diaz, G.; Coimbra, C.F.M. Nonlinear dynamics and control of a variable order oscillator with application to the van der Pol equation. Nonlinear Dynam. 2009, 56, 145–157. [Google Scholar] [CrossRef]
  33. Zhang, H.; Liu, F.; Phanikumar, M.S.; Meerschaert, M.M. A novel numerical method for the time variable fractional order mobile-immobile advection-dispersion model. Comput. Math. Appl. 2013, 66, 693–701. [Google Scholar] [CrossRef]
  34. Benkerrouche, A.; Baleanu, D.; Souid, M.S.; Hakem, A.; Inc, M. Boundary value problem for nonlinear fractional differential equations of variable order via Kuratowski MNC technique. Adv. Differ. Equ. 2021, 365, 1–19. [Google Scholar] [CrossRef]
  35. Benkerrouche, A.; Souid, M.S.; Etemad, S.; Hakem, A.; Agarwal, P.; Rezapour, S.; Ntouyas, S.K.; Tariboon, J. Qualitative study on solutions of a Hadamard variable order boundary problem via the Ulam-Hyers-Rassias stability. Fractal Fract. 2021, 5, 108. [Google Scholar] [CrossRef]
  36. Benkerrouche, A.; Souid, M.S.; Karapinar, E.; Hakem, A. On the boundary value problems of Hadamard fractional differential equations of variable order. Math. Meth. Appl. Sci. 2022, 46, 3187–3203. [Google Scholar] [CrossRef]
  37. Benkerrouche, A.; Souid, M.S.; Sitthithakerngkiet, K.; Hakem, A. Implicit nonlinear fractional differential equations of variable order. Bound. Value Probl. 2021, 2021, 64. [Google Scholar] [CrossRef]
  38. Benkerrouche, A.; Souid, M.S.; Stamov, G.; Stamova, I. Multiterm impulsive Caputo–Hadamard type differential equations of fractional variable order. Axioms 2022, 11, 634. [Google Scholar] [CrossRef]
  39. Bouazza, Z.; Souid, M.S.; Hatira Günerhan, H. Multiterm boundary value problem of Caputo fractional differential equations of variable order. Adv. Differ. Equ. 2021, 2021, 1–17. [Google Scholar] [CrossRef]
  40. Bouazza, Z.; Souid, M.S.; Etemad, S.; Kaabar, M.K.A. Darbo fixed point criterion on solutions of a Hadamard nonlinear variable order problem and Ulam-Hyers-Rassias stability. J. Funct. Spaces 2022, 2022, 1769359. [Google Scholar]
  41. Refice, A.; Souid, M.S.; Stamova, I. On the boundary value problems of Hadamard fractional differential equations of variable order via Kuratowski MNC technique. Mathematics 2021, 9, 1134. [Google Scholar] [CrossRef]
  42. Rezapour, S.; Souid, M.S.; Etemad, S.; Bouazza, Z.; Ntouyas, S.K.; Asawasamrit, S.; Tariboon, J. Mawhin continuation technique for a nonlinear BVP of variable order at resonance via piece-wise constant functions. Fractal Fract. 2021, 5, 216–230. [Google Scholar] [CrossRef]
  43. Zhang, S. Existence of solutions for two-point boundary-value problems with singular differential equations of variable order. Electron. J. Differ. Equ. 2013, 2013, 1–16. [Google Scholar]
  44. Kesavan, S. Functional Analysis, 1st ed.; Hindustan Book Agency: New Delhi, India, 2009; ISBN 978-93-86279-42-2. [Google Scholar]
  45. Alsaedi, A.; Ahmad, B.; Alghamdi, B.; Ntouyas, S.K. On a nonlinear system of Riemann–Liouville fractional differential equations with semi-coupled integro-multipoint boundary conditions. Open Math. 2021, 19, 760–772. [Google Scholar] [CrossRef]
  46. Abbas, S.; Benchohra, M.; Graef, J.R.; Henderson, J. Implicit Fractional Differential and Integral Equations: Existence and Stability, 1st ed.; Walter de Gruyter: Berlin, Germany, 2018; ISBN 9783110553130. [Google Scholar]
  47. Benchohra, M.; Bouriah, S.; Nieto, J.J. Existence and Ulam stability for nonlinear implicit differential equations with Riemann–Liouville fractional derivative. Demonstr. Math. 2019, 52, 437–450. [Google Scholar] [CrossRef]
  48. Luca, R. On a system of Riemann–Liouville fractional differential equations with coupled nonlocal boundary conditions. Adv. Differ. Equ. 2021, 2021, 134. [Google Scholar] [CrossRef]
  49. Corduneanu, C. Principles of Differential and Integral Equations, 2nd ed.; AMS Chelsea Publishing: Providence, RI, USA, 2008; ISBN 978-0821846223. [Google Scholar]
Figure 1. The approximate solution y ( t ) in [ 0 , 1 ] with β ( t , y ) = t 3 3 + t 3 + 1 3 ( y 2 + 1 ) .
Figure 1. The approximate solution y ( t ) in [ 0 , 1 ] with β ( t , y ) = t 3 3 + t 3 + 1 3 ( y 2 + 1 ) .
Axioms 12 00895 g001
Figure 2. The function β ( t , y ) for t [ 0 , 1 ] and y [ 0 , 20 ] .
Figure 2. The function β ( t , y ) for t [ 0 , 1 ] and y [ 0 , 20 ] .
Axioms 12 00895 g002
Figure 3. The function ψ ( t , y ) for t [ 0 , 1 ] and y [ 0 , 20 ] .
Figure 3. The function ψ ( t , y ) for t [ 0 , 1 ] and y [ 0 , 20 ] .
Axioms 12 00895 g003
Table 1. Some values of y ( t ) and β ( t , y ) for t [ 0 , 1 ] .
Table 1. Some values of y ( t ) and β ( t , y ) for t [ 0 , 1 ] .
t0.10.20.30.40.50.60.70.80.91.0
y ( t ) 0.230.721.63.15.157.710.513.3215.918
β ( t , y ( t ) ) 0.350.30.200.20.220.30.350.440.540.7
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Abdelhamid, H.; Stamov, G.; Souid, M.S.; Stamova, I. New Results Achieved for Fractional Differential Equations with Riemann–Liouville Derivatives of Nonlinear Variable Order. Axioms 2023, 12, 895. https://doi.org/10.3390/axioms12090895

AMA Style

Abdelhamid H, Stamov G, Souid MS, Stamova I. New Results Achieved for Fractional Differential Equations with Riemann–Liouville Derivatives of Nonlinear Variable Order. Axioms. 2023; 12(9):895. https://doi.org/10.3390/axioms12090895

Chicago/Turabian Style

Abdelhamid, Hallouz, Gani Stamov, Mohammed Said Souid, and Ivanka Stamova. 2023. "New Results Achieved for Fractional Differential Equations with Riemann–Liouville Derivatives of Nonlinear Variable Order" Axioms 12, no. 9: 895. https://doi.org/10.3390/axioms12090895

APA Style

Abdelhamid, H., Stamov, G., Souid, M. S., & Stamova, I. (2023). New Results Achieved for Fractional Differential Equations with Riemann–Liouville Derivatives of Nonlinear Variable Order. Axioms, 12(9), 895. https://doi.org/10.3390/axioms12090895

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