Next Article in Journal
Stabilization of a Non-linear Fractional Problem by a Non-Integer Frictional Damping and a Viscoelastic Term
Next Article in Special Issue
The Effect of Fractional Derivatives on Thermo-Mechanical Interaction in Biological Tissues during Hyperthermia Treatment Using Eigenvalues Approach
Previous Article in Journal
Dynamic Properties for a Second-Order Stochastic SEIR Model with Infectivity in Incubation Period and Homestead-Isolation of the Susceptible Population
Previous Article in Special Issue
Lie Symmetries and Third- and Fifth-Order Time-Fractional Polynomial Evolution Equations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Combined Liouville–Caputo Fractional Differential Equation

by
McSylvester Ejighikeme Omaba
*,
Hamdan Al Sulaimani
,
Soh Edwin Mukiawa
,
Cyril Dennis Enyi
and
Tijani Abdul-Aziz Apalara
Department of Mathematics, College of Science, University of Hafr Al Batin, P.O. Box 1803, Hafr Al Batin 31991, Saudi Arabia
*
Author to whom correspondence should be addressed.
Fractal Fract. 2023, 7(5), 366; https://doi.org/10.3390/fractalfract7050366
Submission received: 28 January 2023 / Revised: 4 March 2023 / Accepted: 29 March 2023 / Published: 28 April 2023
(This article belongs to the Special Issue Advances in Fractional Order Derivatives and Their Applications)

Abstract

:
This paper studies a fractional differential equation combined with a Liouville–Caputo fractional differential operator, namely, L C D η β , γ Q ( t ) = λ ϑ ( t , Q ( t ) ) , t [ c , d ] , β , γ ( 0 , 1 ] , η [ 0 , 1 ] , where Q ( c ) = q c is a bounded and non-negative initial value. The function ϑ : [ c , d ] × R R is Lipschitz continuous in the second variable, λ > 0 is a constant and the operator L C D η β , γ is a convex combination of the left and the right Liouville–Caputo fractional derivatives. We study the well-posedness using the fixed-point theorem, estimate the growth bounds of the solution and examine the asymptotic behaviours of the solutions. Our findings are illustrated with some analytical and numerical examples. Furthermore, we investigate the effect of noise on the growth behaviour of the solution to the combined Liouville–Caputo fractional differential equation.

1. Introduction

Fractional-order derivatives are known to give more accurate and realistic mathematical models when compared to the classical-order models. The Liouville–Caputo fractional derivative is unarguably a useful operator that mostly models nonlocal behaviours by fractional DEs [1]. The Liouville–Caputo fractional derivative and others alike have found applications in different science and engineering fields and have been used to model many real-life problems. For example, they are used in the mathematical modelling of a human brain tissue (HBT) constitutive model in the framework of anisotropic hyperelasticity [2]; modelling the growth of many economical processes, specifically memory effect on the economic growth model, that is, in the application of economic growth models with memory effect [3] in physics and the environment; studying the chaotic behaviour(s) of dynamical systems; and developing the fractional-order models of neurons [4] and porous media, among others. See also [5,6] for other fractional models. It is worthy of note that all of the above fractional differential equations and more studied in the literature make use of one-sided (left- or right-sided) fractional derivative operators.
The new operator is a convex combination of left and right Liouville–Caputo fractional operators. There is little that one can find in the literature regarding this combined fractional derivative operators. Importantly, this new combined fractional operator is more general than other fractional derivatives [7]. The new operators were studied and defined by Malinowska and Torres [7] as follows:
L C D η β , γ = η c L C D t β + ( 1 η ) t L C D d γ , β , γ ( 0 , 1 ] & η [ 0 , 1 ] ,
where
L C D 0 β , γ Q ( t ) = t L C D d γ Q ( t ) ,
and
L C D 1 β , γ Q ( t ) = c L C D t β Q ( t ) .
Another advantage of this new fractional derivative L C D η β , γ is that it can describe variational problems in a broad perspective [7]. By drawing inspirations from the diamond-alpha derivative on a time scale, which is a linear combination of the delta and nabla derivatives [8,9,10], the model (1) was birthed. Malinowska and Torres [8] showed that the approximation of exact derivatives by the diamond-alpha derivative was better that those of the delta and nabla derivatives.
Therefore, this research considered a special operator by studying the combined Liouville–Caputo fractional differential equation as follows:
L C D η β , γ Q ( t ) = λ ϑ ( t , Q ( t ) ) , β , γ ( 0 , 1 ] & η [ 0 , 1 ] ,
where Q ( c ) = q c is the initial condition, ϑ : [ c , d ] × R R is Lipschitz in the second variable, λ > 0 and L C D η β , γ is known as the combined Liouville–Caputo fractional derivative operator. The operator L C D η β , γ is a convex combination of left and right Liouville–Caputo fractional derivatives. Thus, our contribution and aim in this research paper was to examine the qualitative properties, the well-posedness, estimates of the growth bounds, and the asymptotic behaviour of the solution to a class of combined fractional differential equations.
The organization of the paper is as follows. Section 2 contains the preliminary concepts and definitions needed in the article; Section 3 contains the main results—existence, uniqueness, upper growth estimate, asymptotic behaviour of solution to the combined L–C fractional differential equation. Numerical and analytical illustration of our results are given in Section 4. In Section 5, we consider the stochastic (non-deterministic) case of our equation. Summary of the paper is given in Section 6.

2. Preliminary Concepts

Here, we give some definitions and basic materials. Readers can refer to [5] for further materials on fractional calculus.
Definition 1
([7]). For β , γ ( 0 , 1 ) and 0 η 1 , the combined Liouville–Caputo fractional derivative operator L C D η β , γ is a convex combination of left and right Liouville–Caputo fractional derivatives, defined by
L C D η β , γ = η c L C D t β + ( 1 η ) t L C D d γ .
Definition 2
([11]). Let 0 < c < d and f : [ c , d ] R be an integrable function. The left-sided Katugampola fractional integral of order β > 0 and parameter σ > 0 is given by
I c + β , σ u ( t ) = σ 1 β Γ ( β ) c t s σ 1 ( t σ s σ ) σ 1 u ( s ) d s ,
provided the integral converges.
For σ = 1 , we define the fractional integral as follows.
Definition 3
([11]). Let u : [ c , d ] R be an integrable function where c , d > 0 with c < d . Then,
c I t β u ( t ) = 1 Γ ( β ) c t ( t s ) β 1 u ( s ) d s ,
is the left-sided Riemann–Liouville fractional integral of u of order 0 < β < 1 , provided the integral converges.
Definition 4
([5]). Let u : [ c , d ] R , then
c D t β u ( t ) = 1 Γ ( 1 β ) d d t c t ( t s ) β u ( s ) d s ,
is the Riemann–Liouville fractional derivative of u of order 0 < β < 1 provided the integral converges.
Definition 5
([11]). Let c , d > 0 with c < d , σ > 0 , β R + and m N satisfying m 1 < β < m . Let u : [ c , d ] R be a C m -function. Then,
C D c + β , σ u ( t ) = I c + β , σ t 1 σ d d t m u ( t ) = σ 1 m + β Γ ( m β ) c t s σ 1 ( t σ s σ ) β 1 t 1 σ d d s m u ( s ) d s
is the left-sided Caputo–Katugampola fractional derivative of u of order β and parameter σ.
For σ = 1 and m = 1 , we define the Liouville–Caputo fractional derivative as follows.
Definition 6.
Given c , d > 0 with c < d , β R + so that 0 < β < 1 , and a C 1 -function u : [ c , d ] R . The left-sided Liouville–Caputo fractional derivative of order β is given by
c L C D t β u ( t ) = 1 Γ ( 1 β ) c t ( t s ) β u ( s ) d s ,
provided that the integral converges. Similarly, the right-sided Liouville–Caputo fractional derivative of order γ is given by
t L C D d γ u ( t ) = 1 Γ ( 1 β ) t d ( s t ) γ u ( s ) d s ,
provided that the integral converges.
Lemma 1
([12]). The Liouville–Caputo derivative is connected to the Riemann–Liouville derivatives by
c L C D t β u ( t ) = c D t β u ( t ) k = 0 m 1 ( t c ) k β Γ ( k β + 1 ) u ( k ) ( c ) , m 1 < β m .
In particular, when m = 1 , one obtains
c L C D t β u ( t ) = c D t β u ( t ) ( t c ) β Γ ( 1 β ) u ( c ) .
Suppose the function u : [ c , d ] R is defined as u ( t ) = t a . Then,
c D t β u ( t ) = Γ ( a + 1 ) Γ ( a β + 1 ) t a β .
Theorem 1
([11]). Suppose u C m [ c , d ] , then for m 1 < β m , 0 < σ 1
I c + β , σ C D c + β , σ u ( t ) = u ( t ) k = 0 m 1 σ k k ! ( t σ c σ ) k u ( k ) ( c ) .
Note that if u C 1 [ c , d ] and σ = 1 , then
I c + β L C D c + β u ( t ) = u ( t ) u ( c ) .
Thus,
Corollary 1.
Suppose u C 1 [ c , d ] , 0 < β 1 and σ = 1 , we have
c I t β c L C D t β u ( t ) = u ( t ) u ( c ) .
The next result is a formula for the generalized fractional integration by parts:
Theorem 2
([11]). Suppose u C [ c , d ] and v C m [ c , d ] . Then, for m 1 < β m , 0 < σ 1 ,
c d u ( t ) C D c + β , σ v ( t ) d t = c d t σ 1 v ( t ) D d β , σ ( t 1 σ ) u ( t ) d t + j = 0 m 1 t 1 σ d d t j I d m β , σ ( t 1 β u ( t ) ) v ( m j 1 ) ( t ) t = c t = d .
Particularly for m = 1 and σ = 1 , we have
c d u ( t ) L C D c + β v ( t ) d t = c d v ( t ) D d β u ( t ) d t .
Therefore,
Corollary 2.
Suppose u C [ c , d ] and v C 1 [ c , d ] , then for 0 < γ 1 and σ = 1 ,
c d u ( t ) t L C D d γ v ( t ) d t = c d v ( t ) c D t γ u ( t ) d t .

Formulation of the Solution

Here, we apply the properties or relationship between the fractional integral and fractional differential operators in Equation (4) to make sense of the solution to problem (2).
Lemma 2.
Let η ( 0 , 1 ] . Then, the solution to fractional differential Equation (2) is defined as
Q ( t ) = q c + 1 η η Γ ( β γ ) c t ( t s ) β γ 1 Q ( s ) d s + λ η Γ ( β ) c t ( t s ) β 1 ϑ ( s , Q ( s ) ) d s .
Proof. 
The application of the integral operator c I t β to Equation (2) on both sides gives
c I t β L C D η β , γ Q ( t ) = λ c I t β [ ϑ ( t , Q ( t ) ) ] .
That is,
c I t β η c L C D t β Q ( t ) + ( 1 η ) t L C D d γ Q ( t ) = λ c I t β [ ϑ ( t , Q ( t ) ) ] .
By the linearity of the operator c I t β , we have
η c I t β c L C D t β Q ( t ) + ( 1 η ) c I t β t L C D d γ Q ( t ) = λ c I t β [ ϑ ( t , Q ( t ) ) ] .
From Equation (4) in Corollary 1 and Equation (3), we get
η [ Q ( t ) q c ] + 1 η Γ ( β ) c t ( t s ) β 1 s L C D d γ Q ( s ) d s = λ Γ ( β ) c t ( t s ) β 1 ϑ ( s , Q ( s ) ) d s .
By Equation (5) in Corollary 2, we obtain
η [ Q ( t ) q c ] + 1 η Γ ( β ) c t Q ( s ) c D s γ [ ( t s ) β 1 ] d s = λ Γ ( β ) c t ( t s ) β 1 ϑ ( s , Q ( s ) ) d s .
The derivative c D s γ ( t s ) β 1 is evaluated as,
c D s γ ( t s ) β 1 = Γ ( β 1 + 1 ) Γ ( β 1 γ + 1 ) ( t s ) β γ 1 = Γ ( β ) Γ ( β γ ) ( t s ) β γ 1 .
From Equation (6), one gets
η [ Q ( t ) q c ] 1 η Γ ( β γ ) c t ( t s ) β γ 1 Q ( s ) d s = λ Γ ( β ) c t ( t s ) β 1 ϑ ( s , Q ( s ) ) d s .
Thus, for η ( 0 , 1 ] , one obtains
Q ( t ) = q c + 1 η η Γ ( β γ ) c t ( t s ) β γ 1 Q ( s ) d s + λ η Γ ( β ) c t ( t s ) β 1 ϑ ( s , Q ( s ) ) d s .
The following sup-norm on Q defined by
Q : = sup c t d | Q ( t ) | ,
is useful in the next section.

3. Main Results

The following global Lipschitz condition on the function ϑ ( . , Q ) with respect to the second variable is important for establishing our main result.
Condition 1.
Suppose 0 < Lip ϑ < , t [ c , d ] and Q , R R , we assume
| ϑ ( t , Q ) ϑ ( t , R ) | Lip ϑ | Q R | .
For convenience, we shall set  ϑ ( t , 0 ) = 0  in our computations.

3.1. Well-Posedness

Here, we apply the Banach’s fixed-point theorem to prove the existence and uniqueness of the solution to our problem (2). We begin by defining the operator
A Q ( t ) = q c + 1 η η Γ ( β γ ) c t ( t s ) β γ 1 Q ( s ) d s + λ η Γ ( β ) c t ( t s ) β 1 ϑ ( s , Q ( s ) ) d s ,
and show that A has a unique fixed point which gives the solution to problem (2).
Lemma 3.
Suppose Q is a solution to problem (2) and Condition 1 holds. Then, for 0 < γ < β < 1 and η ( 0 , 1 ] , we have
A Q q c + c 1 Q ,
with positive constant
c 1 : = 1 η η Γ ( β γ + 1 ) ( d c ) β γ + λ Lip ϑ η Γ ( β + 1 ) ( d c ) β < .
Proof. 
Taking the absolute value of Equation (8) leads to
| A Q ( t ) | q c + 1 η η Γ ( β γ ) c t ( t s ) β γ 1 | Q ( s ) | d s + λ η Γ ( β ) c t ( t s ) β 1 | ϑ ( s , Q ( s ) ) | d s .
Applying Equation (7) of Condition 1, we have
| A Q ( t ) | q c + 1 η η Γ ( β γ ) c t ( t s ) β γ 1 | Q ( s ) | d s + λ Lip ϑ η Γ ( β ) c t ( t s ) β 1 | Q ( s ) | d s q c + 1 η η Γ ( β γ ) Q c t ( t s ) β γ 1 d s + λ Lip ϑ η Γ ( β ) Q c t ( t s ) β 1 d s = q c + 1 η η Γ ( β γ ) Q ( t c ) β γ β γ + λ Lip ϑ η Γ ( β ) Q ( t c ) β β = q c + 1 η η Γ ( β γ + 1 ) ( t c ) β γ + λ Lip ϑ η Γ ( β + 1 ) ( t c ) β Q .
Next, we take the supremum over t [ c , d ] to get
A Q   q c + 1 η η Γ ( β γ + 1 ) ( d c ) β γ + λ Lip ϑ η Γ ( β + 1 ) ( d c ) β Q ,
and the result is obtained. □
Lemma 4.
Let Condition 1 hold. Suppose Q 1 and Q 2 are two solutions of (2); for 0 < γ < β < 1 and η ( 0 , 1 ] , we have
A Q 1 A Q 2 c 1 Q 1 Q 2 .
Proof. 
Proceeding with similar steps as in the proof of Lemma 3, we arrive at the desired result. □
Theorem 3.
Suppose Condition 1 holds, 0 < γ < β < 1 and η ( 0 , 1 ] . There exists a constant 0 < c 1 < 1 such that Equation (2) possesses a unique solution.
Proof. 
Applying the Banach fixed point theorem gives Q ( t ) = A Q ( t ) . It follows from Equation (9) in Lemma 3 that
Q = A Q q c + c 1 Q .
This gives Q 1 c 1 q c and therefore, Q < whenever c 1 < 1 .
Next, suppose for a contradiction that Q 1 Q 2 are solutions of Equation (2). Then, from Equation (10) of Lemma 4, we have
Q 1 Q 2 = A Q 1 A Q 2 c 1 Q 1 Q 2 .
Therefore, Q 1 Q 2 1 c 1 0 . However, 1 c 1 > 0 , therefore Q 1 Q 2 < 0 . This is a contradiction, hence Q 1 Q 2 = 0 . By the contraction principle, the existence and uniqueness result follows. □

3.2. Upper Growth Bound

Agarwal et al. [13] presented the following retarded Gronwall type inequality:
u ( t ) f ( t ) + i = 1 n r i ( t 0 ) r i ( t ) h i ( t , s ) w i ( u ( s ) ) d s , t 0 t < t 1 .
Theorem 4
(Theorem 2.1 of [13]). Assume the hypotheses of (Theorem 2.1 of [13]) hold and u ( t ) is a non-negative continuous function on [ t 0 , t 1 ) satisfying (11). Then,
u ( t ) X n 1 X n ( q n ( t ) ) + r n ( t 0 ) r n ( t ) max t 0 τ t h n ( τ , s ) d s , t 0 t T 1 ,
where q n ( t ) is determined recursively by
q 1 ( t ) : = f ( t 0 ) + t 0 t | f ( s ) | d s ,
q i + 1 : = X i 1 X i ( q i ( t ) ) + r i ( t 0 ) r i ( t ) max t 0 τ t h i ( τ , s ) d s , i = 1 , . . . , n 1 ,
and X i ( τ , τ i ) : = τ i τ d ω w i ( ω ) .
Remark 1.
For the case n = 2 , if
u ( t ) f ( t ) + r 1 ( t 0 ) r 1 ( t ) h 1 ( t , s ) w 1 ( u ( s ) ) d s + r 2 ( t 0 ) r 2 ( t ) h 2 ( t , s ) w 2 ( u ( s ) ) d s ,
then
u ( t ) X 2 1 X 2 ( q 2 ( t ) ) + r 2 ( t 0 ) r 2 ( t ) max t 0 τ t h 2 ( τ , s ) d s ,
with q 2 ( t ) = X 1 1 X 1 ( r 1 ( t ) ) + r 1 ( t 0 ) r 1 ( t ) max t 0 τ t h 1 ( τ , s ) d s .
Here, we take w i ( u ( s ) ) = u ( s ) , r i ( t 0 ) = t 0 = c , and r i ( t ) = t for i = 1 , 2 .
Consequently, we present the upper growth bound estimate for the solution.
Theorem 5.
Assume Condition 1 holds. Then, t [ c , d ] , 0 < c < d , and c 2 , c 3 > 0 , we have
| Q ( t ) | q c exp c 2 ( c t ) β γ + c 3 ( c t ) β ,
with c 2 = 1 η η Γ ( β γ + 1 ) , c 3 = λ Lip ϑ η Γ ( β + 1 ) , for η ( 0 , 1 ] , 0 < γ < β < 1 .
Proof. 
We already obtained from Lemma 3 that
| Q ( t ) | q c + 1 η η Γ ( β γ ) c t ( t s ) β γ 1 | Q ( s ) | d s + λ Lip ϑ η Γ ( β ) c t ( t s ) β 1 | Q ( s ) | d s .
Let g ( t ) : = | Q ( t ) | , t [ c , d ] ; it follows that
g ( t ) q c + 1 η η Γ ( β γ ) c t ( t s ) β γ 1 g ( s ) d s + λ Lip ϑ η Γ ( β ) c t ( t s ) β 1 g ( s ) d s .
Applying Theorem 4 to (12), we take w i ( ω ) = ω for i = 1 , 2 and it follows that
X 2 ( τ , τ 2 ) = τ 2 τ d ω ω = ln τ τ 2 .
If one takes τ 2 = 1 for convenience, then X 2 ( τ ) = ln τ with the inverse X 2 1 ( τ ) = e τ . Similarly, X 1 ( τ ) = ln τ has the inverse X 1 1 ( τ ) = e τ . Moreover, f ( t ) = q c and f ( t ) = 0 , hence q 1 ( t ) = q c . Now, we define the non-negative functions h 1 , h 2 : [ c , d ] × [ c , d ] R + as follows:
h 1 ( ζ , s ) : = 1 η η Γ ( β γ ) ( ζ s ) β γ 1 ,
and
h 2 ( ζ , s ) : = λ Lip ϑ η Γ ( β ) ( ζ s ) β 1 .
Let c s < ζ ; since β γ 1 < 0 , then h 1 is decreasing and continuous, thus
max c τ t h 1 ( ζ , s ) = 1 η η Γ ( β γ ) ( c s ) β γ 1 ,
and it follows that
q 2 ( t ) = exp ln ( q c ) + 1 η η Γ ( β γ ) c t ( c s ) β γ 1 d s = exp ln ( q c ) 1 η η Γ ( β γ ) ( c t ) β γ β γ = exp ln ( q c ) 1 η η Γ ( β γ + 1 ) ( c t ) β γ .
In addition, for c s < ζ , and for all β < 1 , h 2 is decreasing and continuous, and
max c ζ t h 2 ( ζ , s ) = λ Lip ϑ η Γ ( β ) ( c s ) β 1 .
Thus,
g ( t ) exp ln ( q 2 ( t ) ) + λ Lip ϑ η Γ ( β ) c t ( c s ) β 1 d s = exp ln ( q c ) 1 η η Γ ( β γ + 1 ) ( c t ) β γ λ Lip ϑ η Γ ( β ) ( c t ) β β = q c exp 1 η η Γ ( β γ + 1 ) ( c t ) β γ λ Lip ϑ η Γ ( β + 1 ) ( c t ) β ,
and this completes the proof. □

3.3. Asymptotic Behaviours

Here, we show that the solution exhibits some asymptotic properties. By the growth bound in Theorem 5, we have
| Q ( t ) | q c exp 1 η 1 η Γ ( β γ + 1 ) ( c t ) β γ + λ Lip ϑ Γ ( β + 1 ) ( c t ) β .
Now, taking the limit as η 0 , we obtain
lim η 0 | Q ( t ) | = 0 .
Next, taking limit as t c + , we get
lim t c + | Q ( t ) | q c .

4. Examples

The example below illustrates Theorem 3 as follows: For β = 9 10 , γ = 1 10 , η = 1 2 , and the nonlinear Lipschitz continuous function ϑ : [ 0.01 , 0.05 ] × R R defined by ϑ ( t , Q ( t ) ) = sin ( Q ( t ) ) having a Lipschitz constant Lip ϑ = 1 , the Combined Liouville–Caputo fractional differential equation
L C D 1 2 9 10 , 1 10 Q ( t ) = λ sin ( Q ( t ) ) , 0.01 < t 0.05 , Q ( 0.01 ) = q c ,
has a unique solution provided by
c 1 : = ( 1 1 2 ) 1 2 η ( 9 10 1 10 + 1 ) ( 0.04 ) 0.8 + λ 1 2 η ( 9 10 + 1 ) ( 0.04 ) 0.9 < 1 ,
if and only if c 1 = 0.081756 + 0.0551892 λ 0.480883 < 1 . That is, for all λ such that 0 < λ < 8.00099 .

Numerical Comparisons

We present numerical simulations and different plots of the upper growth bound functions q c exp 1 η 1 η Γ ( β γ + 1 ) ( c t ) β γ + λ Lip ϑ η ( β + 1 ) ( c t ) β and compare their behaviours for various values of parameters β , γ , η and λ over different time intervals. See Figure 1, Figure 2, Figure 3, Figure 4, Figure 5 and Figure 6 below.

5. Stochastic Combined Fractional Differential Equation

In this section, we study the effect of an external force (a noise term) on the growth behaviour of Equation (2). Thus, we perturb the combined L–C fractional differential equation with a multiplicative noise term w ˙ ( t ) known as the generalized derivative of a Wiener process w ( t ) (a Gaussian white noise process) and consider the following stochastic combined L–C fractional differential equation:
L C D η β , γ Φ ( t ) = λ ϑ ( t , Φ ( t ) ) w ˙ ( t ) , β , γ ( 0 , 1 ] & η [ 0 , 1 ] ,
where Φ ( c ) = ρ c is the initial condition, ϑ : [ c , d ] × R R is Lipschitz in the second variable, L C D η β , γ is known as the combined Liouville–Caputo fractional derivative operator, w ˙ ( t ) is the noise term and λ > 0 denotes the level of the noise term. For recent work on stochastic fractional differential equations, see [14,15] and their references.
Following the formulation of solution in Lemma 2, the solution of Equation (13) is given by
Φ ( t ) = ρ c + 1 η η Γ ( β γ ) c t ( t s ) β γ 1 Φ ( s ) d s + λ η Γ ( β ) c t ( t s ) β 1 ϑ ( s , Φ ( s ) ) d w ( s ) .

5.1. Well-Posedness of the Solution to Equation (13)

We define the norm of the solution (14) above by
Φ 2 2 : = sup c t d E | Φ ( t ) | 2
and define the operator as follows:
L Φ ( t ) = ρ c + 1 η η Γ ( β γ ) c t ( t s ) β γ 1 Φ ( s ) d s + λ η Γ ( β ) c t ( t s ) β 1 ϑ ( s , Φ ( s ) ) d w ( s ) .
We show that the fixed point of the operator L gives the solution to Equation (13).
Lemma 5.
Suppose Φ is a solution to problem (13) and Condition 1 holds. Then, for 0 < γ < β < 1 and η ( 0 , 1 ] , we have
L Φ 2 2 3 ρ c 2 + c 4 Φ 2 2 ,
with positive constant
c 4 : = 3 ( 1 η ) 2 η 2 Γ 2 ( β γ ) [ 2 ( β γ ) 1 ] ( d c ) 2 ( β γ ) + 3 λ 2 Lip ϑ 2 η 2 Γ 2 ( β ) [ 2 β 1 ] ( d c ) 2 β 1 < .
Proof. 
Take the second moment of Equation (15) to get
E | L Φ ( t ) | 2 3 ρ c 2 + 3 ( 1 η ) 2 η 2 Γ 2 ( β γ ) E | c t ( t s ) β γ 1 Φ ( s ) d s | 2 + 3 λ 2 η 2 Γ 2 ( β ) E | c t ( t s ) β 1 ϑ ( s , Φ ( s ) ) d w ( s ) | 2 .
Applying Holder’s inequality on the first integral and Itô isometry on the second integral, we obtain
E | L Φ ( t ) | 2 3 ρ c 2 + 3 ( 1 η ) 2 η 2 Γ 2 ( β γ ) E c t ( t s ) 2 ( β γ 1 ) d s 1 2 c t | Φ ( s ) | 2 d s 1 2 2 + 3 λ 2 η 2 Γ 2 ( β ) c t ( t s ) 2 ( β 1 ) E | ϑ ( s , Φ ( s ) ) | 2 d s .
Next, we apply Equation (7) of Lipschitz Condition 1 to arrive at
E | L Φ ( t ) | 2 3 ρ c 2 + 3 ( 1 η ) 2 η 2 Γ 2 ( β γ ) c t ( t s ) 2 ( β γ 1 ) d s c t E | Φ ( s ) | 2 d s + 3 λ 2 Lip ϑ η 2 Γ 2 ( β ) c t ( t s ) 2 ( β 1 ) E | Φ ( s ) | 2 d s 3 ρ c 2 + 3 ( 1 η ) 2 η 2 Γ 2 ( β γ ) c t ( t s ) 2 ( β γ 1 ) d s Φ 2 2 c t 1 d s + 3 λ 2 Lip ϑ η 2 Γ 2 ( β ) Φ 2 2 c t ( t s ) 2 ( β 1 ) d s = 3 ρ c 2 + 3 ( 1 η ) 2 η 2 Γ 2 ( β γ ) ( t c ) 2 ( β γ ) 1 2 ( β γ ) 1 Φ 2 2 ( t c ) + 3 λ 2 Lip ϑ η 2 Γ 2 ( β ) Φ 2 2 ( t c ) 2 β 1 2 β 1 = 3 ρ c 2 + 3 ( 1 η ) 2 η 2 Γ 2 ( β γ ) ( t c ) 2 ( β γ ) 2 ( β γ ) Φ 2 2 + 3 λ 2 Lip ϑ η 2 Γ 2 ( β ) Φ 2 2 ( t c ) 2 β 1 2 β 1 .
For 0 < γ < β < 1 such that β > 1 2 and β γ > 1 2 , we take the supremum over t [ c , d ] to arrive at
L Φ 2 2 3 ρ c 2 + 3 ( 1 η ) 2 η 2 Γ 2 ( β γ ) ( d c ) 2 ( β γ ) 2 ( β γ ) + 3 λ 2 Lip ϑ η 2 Γ 2 ( β ) ( d c ) 2 β 1 2 β 1 Φ 2 2 ,
and the result follows. □
Similarly, we obtain the following result:
Lemma 6.
Let Condition 1 hold. Suppose Φ and Θ are two solutions of (13); for 0 < γ < β < 1 and η ( 0 , 1 ] , we have
L Φ L Θ 2 2 c 4 Φ Θ 2 2 .
Applying Lemmas 5 and 6, we obtain the following existence and uniqueness result:
Theorem 6.
Suppose Condition 1 holds, 0 < γ < β < 1 and η ( 0 , 1 ] . There exists a constant 0 < c 4 < 1 such that Equation (13) possesses a unique solution.

5.2. Growth Moment Bound

Here, we state and prove the second moment growth estimate for the solution to Equation (13):
Theorem 7.
Assume Condition 1 holds. Then, t [ c , d ] , 0 < c < d and c 5 , c 6 > 0 , we have
E | Φ ( t ) | 2 3 ρ c 2 exp c 5 ( c t ) 2 ( β γ ) 1 + c 6 ( c t ) 2 β 1 ,
with c 5 = 3 ( 1 η ) 2 η 2 Γ 2 ( β γ ) 2 ( β γ ) 1 , c 6 = 3 λ 2 Lip ϑ 2 η 2 Γ 2 ( β ) ( 2 β 1 ) , for η ( 0 , 1 ] , 0 < γ < β < 1 .
Proof. 
Following the proof of Lemma 5, we obtain that
E | Φ ( t ) | 2 3 ρ c 2 + 3 ( 1 η ) 2 η 2 Γ 2 ( β γ ) E | c t ( t s ) β γ 1 Φ ( s ) d s | 2 + 3 λ 2 η 2 Γ 2 ( β ) E | c t ( t s ) β 1 ϑ ( s , Φ ( s ) ) d w ( s ) | 2 .
According to Holder’s inequality and Itô isometry on the first and second integrals, respectively, together with Equation (7) of Condition 1, one arrives at
E | Φ ( t ) | 2 3 ρ c 2 + 3 ( 1 η ) 2 η 2 Γ 2 ( β γ ) E c t ( t s ) 2 ( β γ 1 ) | Φ ( s ) | 2 d s 1 / 2 c t 1 d s 1 / 2 2 + 3 λ 2 Lip ϑ 2 η 2 Γ 2 ( β ) c t ( t s ) 2 ( β 1 ) E | Φ ( s ) | 2 d s = 3 ρ c 2 + 3 ( 1 η ) 2 η 2 Γ 2 ( β γ ) ( t c ) c t ( t s ) 2 ( β γ 1 ) E | Φ ( s ) | 2 d s + 3 λ 2 Lip ϑ 2 η 2 Γ 2 ( β ) c t ( t s ) 2 ( β 1 ) E | Φ ( s ) | 2 d s 3 ρ c 2 + 3 ( 1 η ) 2 η 2 Γ 2 ( β γ ) ( d c ) c t ( t s ) 2 ( β γ 1 ) E | Φ ( s ) | 2 d s + 3 λ 2 Lip ϑ 2 η 2 Γ 2 ( β ) c t ( t s ) 2 ( β 1 ) E | Φ ( s ) | 2 d s .
The last line follows since t [ c , d ] . Let g ( t ) : = E | Φ ( t ) | 2 , t [ c , d ] ; it follows that
g ( t ) 3 ρ c 2 + 3 ( 1 η ) 2 η 2 Γ 2 ( β γ ) ( d c ) c t ( t s ) 2 ( β γ 1 ) g ( s ) d s + 3 λ 2 Lip ϑ 2 η 2 Γ 2 ( β ) c t ( t s ) 2 ( β 1 ) g ( s ) d s .
Next, we apply the retarded Gronwall-type inequality of Theorem 4 on Equation (16) and follow the proof of Theorem 5. Consequently, the non-negative functions h 1 , h 2 : [ c , d ] × [ c , d ] R + are defined as follows:
h 1 ( ζ , s ) : = 3 ( 1 η ) 2 ( d c ) η 2 Γ 2 ( β γ ) ( ζ s ) 2 ( β γ 1 ) ,
and
h 2 ( ζ , s ) : = 3 λ 2 Lip ϑ 2 η 2 Γ 2 ( β ) ( ζ s ) 2 ( β 1 ) .
For c s < ζ and β γ 1 < 0 , it follows that h 1 is continuous and decreasing and one obtains
max c τ t h 1 ( ζ , s ) = 3 ( 1 η ) 2 ( d c ) η 2 Γ 2 ( β γ ) ( c s ) 2 ( β γ 1 ) .
This gives (for β γ > 1 2 ),
q 2 ( t ) = exp ln ( 3 ρ c 2 ) + 3 ( 1 η ) 2 ( d c ) η 2 Γ 2 ( β γ ) c t ( c s ) 2 ( β γ 1 ) d s = exp ln ( 3 ρ c 2 ) 3 ( 1 η ) 2 ( d c ) η 2 Γ 2 ( β γ ) ( c t ) 2 ( β γ ) 1 2 ( β γ ) 1 .
Similarly, for h 2 , we have that, given c s < ζ and β < 1 , h 2 is decreasing, continuous and
max c ζ t h 2 ( ζ , s ) = 3 λ 2 Lip ϑ 2 η 2 Γ 2 ( β ) ( c s ) 2 ( β 1 ) .
Then, it follows that (for β > 1 2 and β γ > 1 2 )
g ( t ) exp ln ( q 2 ( t ) ) + 3 λ 2 Lip ϑ 2 η 2 Γ 2 ( β ) c t ( c s ) 2 ( β 1 ) d s = exp ln ( 3 ρ c 2 ) 3 ( 1 η ) 2 ( d c ) η 2 Γ 2 ( β γ ) ( c t ) 2 ( β γ ) 1 2 ( β γ ) 1 3 λ 2 Lip ϑ 2 η 2 Γ 2 ( β ) ( c t ) 2 β 1 2 β 1 = 3 ρ c 2 exp 3 ( 1 η ) 2 ( d c ) η 2 Γ 2 ( β γ ) ( c t ) 2 ( β γ ) 1 2 ( β γ ) 1 3 λ 2 Lip ϑ 2 η 2 Γ 2 ( β ) ( c t ) 2 β 1 2 β 1 ,
and this proves the result. □
Remark 2.
The mild solution in Equation (14) also satisfies the same asymptotic properties of Section 3.3. Since
E | Φ ( t ) | 2 3 ρ c 2 exp 3 η 2 ( 1 η ) 2 ( d c ) Γ 2 ( β γ ) ( c t ) 2 ( β γ ) 1 2 ( β γ ) 1 + λ 2 Lip ϑ 2 Γ 2 ( β ) ( c t ) 2 β 1 2 β 1 ,
then, taking limits as η 0 and t c + , we obtain the following: lim η 0 E | Φ ( t ) | 2 = 0 and lim t c + E | Φ ( t ) | 2 3 ρ c 2 .

6. Conclusions

The combined Liouville–Caputo fractional derivative operator plays a major role in describing a more general class of variational problems. We used this new operator to investigate the behaviour of the solution to a class of combined Liouville–Caputo fractional differential equations. The existence and uniqueness of the solution was established through Banach’s fixed-point theorem. The solution’s growth bound and asymptotic behaviours were also established. Moreover, in Figure 1, Figure 2, Figure 3, Figure 4, Figure 5 and Figure 6, we presented a numerical simulation of the behaviour of the combined Liouville–Caputo fractional differential equations via the upper growth bound, with respect to different initial points t = c , different values β and γ of the fractional orders, different convex combination parameters η , as well as different values of λ . In general, the solution function Q of problem (2) decayed to zero; however, the speed of decay was determined by the values of the parameters η , β , γ and λ . We further investigated the effects of a noise term on the growth of the solution to the combined fractional differential equation and observed that the presence or introduction of the noise term does not affect the growth behaviour of the solution to the combined fractional differential equation. For future work, we will study the dependence of the solution(s) on the initial condition and will also estimate the lower bounds, etc.

Author Contributions

Conceptualization, M.E.O.; methodology, M.E.O., H.A.S., S.E.M., C.D.E. and T.A.-A.A.; software, C.D.E.; validation, M.E.O., H.A.S., S.E.M., C.D.E. and T.A.-A.A.; software, C.D.E.; formal analysis, M.E.O., H.A.S., S.E.M., C.D.E. and T.A.-A.A.; investigation, M.E.O., H.A.S., S.E.M., C.D.E. and T.A.-A.A.; writing—original draft preparation, M.E.O.; writing—review and editing, M.E.O., H.A.S., S.E.M., C.D.E. and T.A.-A.A.; project administration, M.E.O.; funding acquisition, M.E.O. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the University of Hafr Al Batin, Institutional Financial Program under project number IFP-A-2022-2-1-09.

Data Availability Statement

Not Applicable.

Acknowledgments

The authors gratefully acknowledge the technical and financial support from the Agency for Research and Innovation, Ministry of Education and University of Hafr Al Batin, Saudi Arabia.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Baleanu, D.; Fernadez, A.; Akgul, A. On a Fractional Combining proportional and Classical Differintegrals. Mathematics 2020, 8, 360. [Google Scholar] [CrossRef]
  2. Voyiadjis, G.Z.; Sumelka, W. Brain modelling in the framework of anisotropic hyperelasticity with time fractional damage evolution governed by the Caputo–Almeida fractional derivative. J. Mech. Behav. Biomed. Mater. 2019, 89, 209–216. [Google Scholar] [CrossRef] [PubMed]
  3. Johansyah, M.D.; Supriatna, A.K.; Rusyaman, E.; Saputra, J. Application of fractional differential equation in economic growth model: A Systematic review approach. Aims Math. 2021, 6, 10266–10280. [Google Scholar] [CrossRef]
  4. Mu’lla, M. Fractional calculus, Fractional Differential Equations and Applications. Open Access Libr. J. 2020, 7, 1–9. [Google Scholar] [CrossRef]
  5. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations, 1st ed.; North–Holland Mathematics Studies; Elsevier: Amsterdam, The Netherlands, 2006; Volume 204, pp. 1–523. [Google Scholar]
  6. Baleanu, D.; Machado, J.A.T.; Luo, A.C. Fractional Dynamics and Control; Springer Science and Business Mdeia: Berlin, Germany, 2011. [Google Scholar]
  7. Malinowska, A.B.; Torres, D.F.M. Fractional calculus of variations for a combined Caputo derivative. Fract. Calc. Appl. Anal. 2011, 14, 523–537. [Google Scholar] [CrossRef]
  8. Malinowska, A.B.; Torres, D.F.M. The diamond-alpha Riemann integral and mean value theorems on time scales. Dynam. Syst. Appl. 2009, 18, 469–482. [Google Scholar]
  9. Sidi Ammi, M.R.; Ferreira, R.A.C.; Torres, D.F.M. Diamond-α Jensen’s Inequality on time scales. J. Inequal. Appl. 2008, 576876. [Google Scholar] [CrossRef]
  10. Atasever, N.A. On Diamond-Alpha Dynamic Equations and Inequalities. Master’s Thesis, Georgia Southern University, Statesboro, GA, USA, 2011. [Google Scholar]
  11. Almeida, R. A Gronwall inequality for a general Caputo fractional operator. Math. Inequalities Appl. 2017, 20, 1089–1105. [Google Scholar] [CrossRef]
  12. Abdeljawad, T. On Riemann and Caputo fractional differences. Comput. Math. Appl. 2011, 62, 1602–1611. [Google Scholar] [CrossRef]
  13. Agarwal, R.P.; Deng, S.; Zhang, W. Generalization of a retarded Gronwall-like inequality and its applications. Appl. Math. Comput. 2005, 165, 599–612. [Google Scholar] [CrossRef]
  14. Omaba, M.E.; Enyi, C.D. Atangana–Baleanu time-fractional stochastic integro-differential equation. Partial. Differ. Equations Appl. Math. 2021, 4, 100100. [Google Scholar] [CrossRef]
  15. Omaba, M.E.; Nweze, E.R. A Nonlinear Fractional Langevin Equation of Two Fractional Orders with Multiplicative Noise. Fractal Fract. 2022, 6, 290. [Google Scholar] [CrossRef]
Figure 1. Behaviour of the upper growth bound for different η and λ values: (a) 0.01 t 10 , β = 0.2 , γ = 0.1 , λ = 1 ; (b) 0.01 t 10 .
Figure 1. Behaviour of the upper growth bound for different η and λ values: (a) 0.01 t 10 , β = 0.2 , γ = 0.1 , λ = 1 ; (b) 0.01 t 10 .
Fractalfract 07 00366 g001
Figure 2. Behaviour of the upper growth bound for different η and λ values: (a) 0.01 t 10 4 , β = 0.2 , γ = 0.1 , λ = 1 ; (b) 0.01 t 10 4 .
Figure 2. Behaviour of the upper growth bound for different η and λ values: (a) 0.01 t 10 4 , β = 0.2 , γ = 0.1 , λ = 1 ; (b) 0.01 t 10 4 .
Fractalfract 07 00366 g002
Figure 3. Behaviour of the upper growth bound for different η and λ values: (a) 0.01 t 10 10 , β = 0.2 , γ = 0.1 , λ = 1 ; (b) 0.01 t 10 10 .
Figure 3. Behaviour of the upper growth bound for different η and λ values: (a) 0.01 t 10 10 , β = 0.2 , γ = 0.1 , λ = 1 ; (b) 0.01 t 10 10 .
Fractalfract 07 00366 g003
Figure 4. Behaviour of the upper growth bound for different η , β and γ values: (a) 0.01 t 10 , λ = 1 ; (b) 0.01 t 10 , λ = 1 .
Figure 4. Behaviour of the upper growth bound for different η , β and γ values: (a) 0.01 t 10 , λ = 1 ; (b) 0.01 t 10 , λ = 1 .
Fractalfract 07 00366 g004
Figure 5. Behaviour of the upper growth bound for different initial points t = c : (a) c t 15 , λ = 1 ; (b) c t 15 , λ = 1 .
Figure 5. Behaviour of the upper growth bound for different initial points t = c : (a) c t 15 , λ = 1 ; (b) c t 15 , λ = 1 .
Fractalfract 07 00366 g005
Figure 6. Behaviour of the upper growth bound for different initial points t = c : (a) c t 10 10 , λ = 1 ; (b) c t 10 10 , λ = 1 .
Figure 6. Behaviour of the upper growth bound for different initial points t = c : (a) c t 10 10 , λ = 1 ; (b) c t 10 10 , λ = 1 .
Fractalfract 07 00366 g006
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

Omaba, M.E.; Al Sulaimani, H.; Mukiawa, S.E.; Enyi, C.D.; Apalara, T.A.-A. Combined Liouville–Caputo Fractional Differential Equation. Fractal Fract. 2023, 7, 366. https://doi.org/10.3390/fractalfract7050366

AMA Style

Omaba ME, Al Sulaimani H, Mukiawa SE, Enyi CD, Apalara TA-A. Combined Liouville–Caputo Fractional Differential Equation. Fractal and Fractional. 2023; 7(5):366. https://doi.org/10.3390/fractalfract7050366

Chicago/Turabian Style

Omaba, McSylvester Ejighikeme, Hamdan Al Sulaimani, Soh Edwin Mukiawa, Cyril Dennis Enyi, and Tijani Abdul-Aziz Apalara. 2023. "Combined Liouville–Caputo Fractional Differential Equation" Fractal and Fractional 7, no. 5: 366. https://doi.org/10.3390/fractalfract7050366

APA Style

Omaba, M. E., Al Sulaimani, H., Mukiawa, S. E., Enyi, C. D., & Apalara, T. A. -A. (2023). Combined Liouville–Caputo Fractional Differential Equation. Fractal and Fractional, 7(5), 366. https://doi.org/10.3390/fractalfract7050366

Article Metrics

Back to TopTop