Next Article in Journal
Advancing Model Generalization in Continuous Cyclic Test-Time Adaptation with Matrix Perturbation Noise
Next Article in Special Issue
Pre-Compactness of Sets and Compactness of Commutators for Riesz Potential in Global Morrey-Type Spaces
Previous Article in Journal
Singularity, Observability, and Independence: Unveiling Lorenz’s Cryptographic Potential
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Unique Solutions for Caputo Fractional Differential Equations with Several Delays Using Progressive Contractions

by
Cemil Tunç
1,† and
Fahir Talay Akyildiz
2,*,†
1
Department of Mathematics, Faculty of Sciences, Van Yuzuncu Yıl University, 65080 Van, Turkey
2
Department of Mathematics and Statistics, Faculty of Sciences, Imam Mohammad Ibn Saud Islamic University (IMSIU), Riyadh 11564, Saudi Arabia
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Mathematics 2024, 12(18), 2799; https://doi.org/10.3390/math12182799
Submission received: 31 July 2024 / Revised: 8 September 2024 / Accepted: 9 September 2024 / Published: 10 September 2024
(This article belongs to the Special Issue Advances in Mathematics: Equations, Algebra, and Discrete Mathematics)

Abstract

:
We take into account a nonlinear Caputo fractional-order differential equation including several variable delays. We examine whether the solutions to the Caputo fractional-order differential equation taken under consideration, which has numerous variable delays, are unique. In the present study, first, we will apply the method of progressive contractions, which belongs to T.A. Burton, to Caputo fractional-order differential equation, including multiple variable delays, which has not yet appeared in the relevant literature by this time. The significant point of the method of progressive contractions consists of a very flexible idea to discuss the uniqueness of solutions for various mathematical models. Lastly, we provide two examples to demonstrate how this paper’s primary outcome can be applied.

1. Introduction

In the majority of the current literature, authors prove the existence and uniqueness of solutions of differential equations by writing them as integral equations and applying some type of fixed-point theorem that can be tedious and challenging, often patching together solutions on short intervals after making complicated translations. Examples of direct fixed-point mappings can be seen in Anderson and Avery [1], Burton [2], Burton and Zhang [3], and Mureşan and Nica [4]. In each case, there are excellent reasons for not first converting DEs to IEs. Across these difficulties, in 2016, Burton [5] developed a fascinating method known as progressive contractions to investigate the global existence and uniqueness of solutions of an integro-differential equation. As a result of the technique of progressive contractions, researchers avoided the process of establishing existence on a little interval, translating the equation to a later starting time, and then connecting a solution on the new interval to the preceding one. Later, the technique of progressive contractions was extended and applied to delay IEs, scalar nonlinear IEs, fractional differential equations of the Riemann–Liouville type, and nonlinear IEs in a Banach space by Burton and Purnaras [6], Burton [7], Burton [8], Ilea and Otrocol [9], respectively. Indeed, this technique is quite effective and adaptable, and easy to investigate the global existence and uniqueness of solutions that of IEs, IDEs, FDEs, etc.; see Burton [5,7,8], Burton and Purnaras [6], Ilea and Otrocol [9], Tunç et al. [10,11], and the references of these papers. In fact, it is highly desirable for researchers to have an appropriate technique, like the progressive contractions technique, for analyzing the uniqueness of solutions in an infinite interval case. We also direct readers to the construction of a fractional non-polynomial spline to approximate solutions of the KdV equation (Yousif and Hamasalh [12]).
We will now briefly review a few works related to different mathematical models where progressive contractions were used as the main techniques in the proofs of the outcomes of those works.
In 2016, Burton [5] presented a straightforward proof regarding the global existence and uniqueness of a solution of the following IDE:
x ( t ) = g ( t , x ( t ) ) + 0 t A ( t s ) f ( s , x ( s ) ) d s .
The proofs of Burton [5], Theorem 2.1, Theorem 2.2 employed the concept of progressive contractions, which is attributed to T. A. Burton. It is an extensive fixed-point theorem for differential equations.
Subsequently, in 2017, Burton and Purnaras [6] expanded the use of progressive contractions to obtain global unique solutions of delay IEs of the following type:
x ( t ) = L ( t ) + g ( t , x ( t ) ) + 0 t A ( t s ) f ( s , x ( s ) ) + f ( s , x ( s r ( s ) ) ) d s .
In Burton and Purnaras [6], Lemma 2.1 is the key result, which extends progressive contractions to delay IEs. Therefore, in Burton and Purnaras [6], considering Lemma 2.1, subsequently utilizing the method of progressive contractions, sufficient conditions were constructed regarding global unique solutions of this delay IE.
In 2017, the following scalar FDE of the Riemann–Liouville type was also taken into consideration by Burton [8]:
D q x ( t ) = f ( t , x ( t ) ) with lim t t 1 q x ( t ) = x 0 .
Burton [8] presented a simple proof of the global existence and uniqueness of a solution of this FDE and then parlayed this into a solution on ( 0 , ) . The proof of [8] utilized the concept of progressive contractions. This technique stays away from the traditional techniques that imply endless repeated translations.
Following that, in 2019, Burton [7] showed a very easy and direct method of obtaining a unique solution on [ 0 , ) for the scalar nonlinear IE that is as follows:
x ( t ) = g ( t , x ( t ) ) + 0 t A ( t s ) f ( s , x ( s ) ) d s .
In Burton [7], once more, the author used the progressive contraction technique to obtain new qualitative outcomes. Burton [7] concluded that there is only one function ξ ( t ) that can satisfy this IE on [ 0 , ) without resorting to any of the classical translations and extensions of solutions, which in actuality must invoke Zorn’s lemma and can encounter challenges.
The progressive contractions method of Burton in [5] was expanded upon in 2020 by Ilea and Otrocol [9] regarding the subsequent IEs in a Banach space:
x ( t ) = 0 t K ( t , s , x ( s ) ) d s
and
x ( t ) = g ( t , x ( t ) ) + 0 t f ( t , s , x ( s ) ) d s .
Ilea and Otrocol [9] obtained sufficient conditions to guarantee aforementioned IEs have a unique solution in a Banach space.
In 2023, Tunç et al. [13] examined the subsequent nonlinear and Hammerstein-type FIE in the Banach space:
x ( t ) = F ( x ( t ) ) + G ( t , x ( t ) ) + H ( t , x ( t ) ) 0 t K ( t , s , x ( s ) ) + Q ( t , s , q ( x ( s ) ) ) d s .
In [13], utilizing the complete metric and the Chebyshev norm, Tunç et al. enhanced Burton’s progressive contractions technique to the aforementioned nonlinear Hammerstein-type FIE in a Banach space. Hence, the authors acquired and enhanced new outcomes for this nonlinear Hammerstein-type FIE.
Motivated by Burton [6,7], Ilea and Otrocol [9] and that given above, Graef et al. [14] considered the Hammerstein-type IE of the form:
x ( t ) = p ( t ) + g ( t , x ( t ) ) + h ( t , x ( t ) ) 0 t A ( t s ) f ( t , s , x ( s ) ) d s .
Graef et al. [14] applied Burton’s method to this Hammerstein-type IE and obtained two new results regarding the existence of solutions.
Recently, Tunç et al. [10] also considered the following nonlinear IE with several variable time delays
x ( t ) = q ( t ) + r ( x ( t ) ) + h ( t , x ( t ) ) + g ( t , x ( t ) ) 0 t A ( t s ) f ( s , x ( s ) ) d s + i = 1 N 0 t A i ( t s ) f i ( s , x ( s ) , x ( s τ i ( s ) ) ) d s
and the nonlinear IDE without delay
x ( t ) = r ( x ( t ) ) + g ( t , x ( t ) ) + h ( t , x ( t ) ) 0 t A ( t s ) f ( s , x ( s ) ) d s .
Tunç et al. [10] studied the global existence of solutions of this IE with several variable time delays and the nonlinear IDE without delay by the fixed-point method using progressive contractions. The results of Tunç et al. [10] expanded upon and enhanced a few associated outcomes from the pertinent literature.
On the other hand, we also direct readers to a thorough discussion for some qualitative behaviors of solutions of various mathematical models as follows: existence and stability of implicit FDEs and IEs (Abbas et al. [15], Benchohra et al. [16]), existence of solutions to boundary value problems (Anderson and Avery [1]), theory and applications of FDEs (Balachandran [17], Diethelm [18], Kilbas et al. [19], Miller and Ross [20], Podlubny [21], Zhou [22]), existence of solutions of nonlinear FDEs of Riemann–Liouville type (Becker et al. [23]), fixed-point theory and applications (Burton [2]), periodicity in DDEs by fixed-point mappings (Burton and Zhang [3]), existence of solutions of IEs (Chauhan et al. [24], Deep et al. [25], Jung [26]), theory of functional differential equations (Hale [27]), existence of solutions of hybrid FDEs (Khan et al. [28,29]), the UHML stability of FDEs with delay (Tunç [30]), solution estimates for CFDEs with delay (Tunç and Tunç [11]), exponential discrete form for Caputo–Fabrizio fuzzy BAM neural networks (Zhang and Li [31]) and the references cited therein.
Inspired by the aforementioned works, especially by Burton [8], we will examine the subsequent CFDE with multiple variable time delays:
D t q 0 C ϑ ( t ) = i = 1 N G i ( t , ϑ ( t ) , ϑ ( θ i ( t ) ) ) ,
x ( t ) = ψ ( t ) , t [ ω , 0 ] ,
where t + , + = [ 0 , ) , ϑ ( t ) , D t q C ϑ ( t ) is the Caputo derivative of q order, q ( 0 , 1 ) . We assume that G i C + × 2 , , θ i C ( + , + ) , θ i ( t ) t with 0 θ i ( t ) ω i , ω i , ω i > 0 , ω = max ( ω i ) , i = 1 , 2 , , N .
Regarding the reason behind the motivation of this study, it is well known that classical calculus is a special case of fractional calculus. For example, any ordinary differential equation describing an electrical circuit has exactly one time-dependent current function as the solution. Nevertheless, depending on the order of the fractional derivative, we obtain many time-dependent current functions when solving the associated differential equation of fractional order. This fact regarding fractional derivatives offers a substantial advantage in real-world applications of physics and electronics, where achieving the desired current function is crucial. Next, there is extensive literature on the existence of solutions of delay systems; however there are extensively open problems for delay systems of fractional order, especially in the case of multiple variable time delays, see the books of Benchohra et al. [16], Balachandran [17], Diethelm [18], Kilbas et al. [19], Miller and Ross [20], Podlubny [21] and the sources in the references of this study. To the best of our knowledge, the database of pertinent literature does not contain any paper or work on the uniqueness of solutions of CFDEs with or without delays, where the technique of progressive contractions is a basic tool. Therefore, it deserves to study the uniqueness of solutions regarding the CFDE (1) with different variable time delays.
The primary goal of this work is to advance and contribute the theory of the existence and uniqueness of FDEs with and without delay, as well as to improve the results of Burton [5,7,8], Burton and Purnaras [6], Graef et al. [14], Ilea and Otrocol [9], and Tunç et al. [10,11] for CFDE (1). Indeed, we aim to acquire new insights in connection with the existence and uniqueness theory of FDEs. Therefore, this study aims to obtain a new result, i.e., Theorem 1, on the topic of CFDE (1) and to provide two new examples in relation to show the applications of this theorem by Examples 1 and 2. These data represent the primary contributions, novelty, and originality of this paper.
This paper’s remaining sections are organized as follows: Section 2 allows some fundamental definitions and information in connection with the theory and application of fractional calculus. Section 3 presents the primary and novel outcome of the study concerning the existence of solutions of the nonlinear CFDE (1) through the application of the progressive contractions technique. Section 4 provides two examples of numerical applications. Section 5 summarizes the outcomes and contributions of the paper. Finally, the paper’s conclusion is given in Section 6.

2. Preliminaries

Within this section, some common terminologies and background knowledge of fractional calculus are provided, as they are necessary for this study.
Definition 1
(Kilbaset al. [19], Podlubny [21]). The Riemann–Liouville fractional integral of order q with the lower limit zero for a function ϑ is defined as
I 0 + q ϑ ( t ) = 1 Γ ( q ) 0 t ϑ ( s ) ( t s ) 1 q d s , t > 0 , q > 0 ,
provided that the right side is point-wise defined on + , + = [ 0 , ) , where Γ ( . ) is the gamma function.
Definition 2
(Kilbas et al. [19], Podlubny [21]). The Caputo derivative of fractional-order q for a function ϑ is defined as
D 0 + q C ϑ ( t ) = 1 Γ ( n q ) 0 t 1 ( t s ) 1 + q n ϑ ( n ) ( s ) d s , t > 0 , q > 0 ,
where n = q + 1 , q represents the integer part of q.

3. Progressive Contractions and Existence of Solutions

We shall demonstrate the uniqueness of solutions of the CFDE (1), including different variable time delays in Theorem 1.
Theorem 1.
Let the following conditions hold:
(Hp1) 0 t E , E > 0 , E , G i C ( [ 0 , E ] × 2 , ) , θ i C ( [ 0 , E ] , + ) , θ i ( t ) t , 0 θ i ( t ) ω i ,
ω i , ω i > 0 , ω = max ( ω i ) , i = 1 , 2 , , N ;
(Hp2) there exist positive constants L G i such that
G i ( t , ψ 1 , ψ 1 ( θ i ( t ) ) ) G i ( t , ψ 2 , ψ 2 ( θ i ( t ) ) ) L G i ψ 1 ψ 2 + ψ 1 ( θ i ( t ) ) ψ 2 ( θ i ( t ) )
for all
t [ 0 , E ] , ψ 1 , ψ 2 , i = 1 , 2 , , N ;
(Hp3)
2 i = 1 N L G i Γ ( q ) × T T 0 q q < β , β ( 0 , 1 ) .
Then, the CFDE (1) has a unique solution.
Remark 1.
Here, T 0 , E is divided into equal portions, each having the endpoints T 0 , T 1 , , T n = E such that E > T 0 and the each part having also the length S < T T 0 , T > 0 . Let q ( 0 , 1 ) and β ( 0 , 1 ) . Therefore, it may be concluded that
2 i = 1 N L G i Γ ( q ) T 0 T T s q 1 d s = 2 i = 1 N L G i Γ ( q ) × T T 0 q q < β , β ( 0 , 1 ) .
We will also need a form in connection with (3). Let T i s = u , ( i = 1 , 2 , , n ) . Then, since T i T i 1 = S , we have
T i 1 T i T i s q 1 d s = 0 T i T i 1 u q 1 d u = S q q < T T 0 q q < Γ ( q ) β 2 i = 1 N L G i .
Remark 2.
Since the functions G i ( t , ϑ ( t ) , ϑ ( θ i ( t ) ) , i = 1 , 2 , , N are continuous and satisfy the Lipschitz condition of (Hp2), then, for each q ( 0 , 1 ) , there exists a T 0 ( 0 , E ) such that Equation (5) has a unique continuous solution ψ ( t ) = ζ on ( 0 , T 0 ] with
lim t 0 0 t ( t s ) q 1 i = 1 N G i s , ϑ ( s ) , ϑ ( θ i ( s ) ) d s = 0 and lim t 0 ϑ ( t ) = ψ ( 0 )
(see also Burton [8]).
Proof. 
It is simple to observe that the CFDE (1) with (2) is equivalent to the singular integral system that follows:
ϑ ( t ) = ψ ( t ) , t [ ω , 0 ] ψ ( 0 ) + 1 Γ ( q ) 0 t ( t s ) q 1 i = 1 N G i ( s , ϑ ( s ) , ϑ ( θ i ( s ) ) ) d s , t [ 0 , E ] .
Step 1. Let 1 , . be a complete metric space with the supremum norm such that it encompasses the functions ψ C ( 0 , T 1 ] , with ψ ( t ) = ζ of Remark 2 on ( 0 , T 0 ] . We describe an operator
P 1 : 1 1 with ψ 1 ,
which allows that
P 1 ϑ ( t ) = ψ ( t ) , t [ ω , 0 ] ψ ( 0 ) + 1 Γ ( q ) 0 t ( t s ) q 1 i = 1 N G i s , ϑ ( s ) , ϑ ( θ i ( s ) ) d s , t ( 0 , T 1 ] .
For t [ ω , 0 ] , from Equation (6) we have
P 1 ψ 1 ( t ) P 1 ψ 2 ( t ) = 0 , ψ 1 , ψ 2 1 .
Next, since ζ satisfies Equation (5) on ( 0 , T 0 ] , i.e., for 0 < t T 0 ,then
P 1 ψ ( t ) = ζ ( t ) = ψ ( 0 ) + 1 Γ ( q ) 0 t ( t s ) q 1 i = 1 N G i s , ζ ( s ) , ζ ( θ i ) d s
for t [ 0 , T 1 ] . Therefore, P 1 of Equation (7) does the map 1 1 .
Let ψ 1 , ψ 2 1 . Then, in the light of Equation (6), it follows that
P 1 ψ 1 ( t ) P 1 ψ 2 ( t ) = 1 Γ ( q ) 0 t ( t s ) q 1 i = 1 N G i s , ψ 1 ( s ) , ψ 1 ( θ i ( s ) ) G i s , ψ 2 ( s ) , ψ 2 ( θ i ( s ) ) d s
(since ψ 1 = ψ 2 = ζ on ( 0 , T 0 ] , using (Hp2) and let t > T 0 )
1 Γ ( q ) T 0 t ( t s ) q 1 i = 1 N L G i ψ 1 ( s ) ψ 2 ( s ) + ψ 1 ( θ i ( s ) ) ψ 2 ( θ i ( s ) ) d s 2 i = 1 N L G i Γ ( q ) ψ 1 ψ 2 T 0 , T 1 T 0 T 1 ( T 1 s ) q 1 d s = 2 i = 1 N L G i Γ ( q ) × ( T 1 T 0 ) q q ψ 1 ψ 2 T 0 , T 1 β ψ 1 ψ 2 ,
where ψ 1 ψ 2 T 0 , T 1 denotes the supremum on the given interval. Therefore, P 1 is a contraction on the set 1 with the unique fixed point ζ 1 , where ζ 1 agrees with ζ on ( 0 , T 0 ] since ζ and ζ 1 are unique on the interval ( 0 , T 0 ] .
Step 2. Let 2 , . be a complete metric space, which encompasses the functions ψ C ( 0 , T 2 ] , such that ψ ( t ) = ζ 1 ( t ) of Step 1 on the interval ( 0 , T 1 ] together with the supremum norm.
We describe a map
P 2 : 2 2 with ψ 2 ,
which allows that
P 2 ϑ ( t ) = ψ ( t ) , t [ ω , 0 ] ψ ( 0 ) + 1 Γ ( q ) 0 t ( t s ) q 1 i = 1 N G i s , ϑ ( s ) , ϑ ( θ i ( s ) ) d s , t ( 0 , T 2 ] .
Hence, for t [ ω , 0 ] , from Equation (8) we have
P 2 ψ 1 ( t ) P 2 ψ 2 ( t ) = 0 , ψ 1 , ψ 2 2 .
Since ζ 1 satisfies Equation (6) on the interval ( 0 , T 1 ] , i.e., for 0 < t T 1 and in 2 , each ψ ( t ) = ζ 1 ( t ) , then
P 2 ψ ( t ) = ζ 1 ( t ) = ψ ( 0 ) + 1 Γ ( q ) 0 t ( t s ) q 1 i = 1 N G i ( s , ζ 1 ( s ) , ζ 1 ( θ i ( s ) ) ) d s , t ( 0 , T 2 ] .
Therefore, P 2 of Equation (9) does the map 2 2 .
Let ψ 1 , ψ 2 2 . Then, according to the conditions of Theorem 1, we derive that
P 2 ψ 1 ( t ) P 2 ψ 2 ( t ) = 1 Γ ( q ) 0 t ( t s ) q 1 i = 1 N G i ( s , ψ 1 ( s ) , ψ 1 ( θ i ( s ) ) ) G i ( s , ψ 2 ( s ) , ψ 2 ( θ i ( s ) ) ) d s
(since ψ 1 = ψ 2 = ζ 1 on ( 0 , T 1 ] , using (Hp2) and let t > T 1 )
1 Γ ( q ) T 1 t ( t s ) q 1 i = 1 N L G i ψ 1 ( s ) ψ 2 ( s ) + ψ 1 ( θ i ( s ) ) ψ 2 ( θ i ( s ) ) d s
2 i = 1 N L G i Γ ( q ) ψ 1 ψ 2 T 1 , T 2 T 1 T 2 ( T 2 s ) q 1 d s
(using (3), (4) and let S = T 2 T 1 )
= 2 i = 1 N L G i Γ ( q ) ψ 1 ψ 2 T 1 , T 2 T 0 T 1 ( T 1 s ) q 1 d s = 2 i = 1 N L G i Γ ( q ) × ( T 1 T 0 ) q q ψ 1 ψ 2 T 0 , T 1 β ψ 1 ψ 2 ,
where ψ 1 ψ 2 T 1 , T 2 denotes the supremum on the given interval. Therefore, P 2 is a contraction on the set 2 with the unique fixed point ζ 2 , where ζ 2 agrees with ζ 1 on ( 0 , T 1 ] since ζ 1 and ζ 2 are unique on the interval ( 0 , T 1 ] .
Step 3. Let 3 , . be a complete metric space, which encompasses the functions ψ C ( 0 , T 3 ] , such that ψ ( t ) = ζ 2 ( t ) of Step 2 on the interval ( 0 , T 2 ] together with the supremum norm.
We describe a map
P 3 : 3 3 with ψ 3 ,
which allows that
P 3 ϑ ( t ) = ψ ( t ) , t [ ω , 0 ] ψ ( 0 ) + 1 Γ ( q ) 0 t ( t s ) q 1 i = 1 N G i s , ϑ ( s ) , ϑ ( θ i ( s ) ) d s , t ( 0 , T 3 ] .
For t [ ω , 0 ] , from Equation (10) we have
P 3 ψ 1 ( t ) P 3 ψ 2 ( t ) = 0 , ψ 1 , ψ 2 3 .
Since ζ 2 satisfies Equation (9) on the interval ( 0 , T 2 ] , i.e., for 0 < t T 2 and in 3 , each ψ ( t ) = ζ 2 ( t ) , then
P 3 ψ ( t ) = ζ 2 ( t ) = ψ ( 0 ) + 1 Γ ( q ) 0 t ( t s ) q 1 i = 1 N G i ( s , ζ 2 ( s ) , ζ 2 ( θ i ( s ) ) ) d s , t ( 0 , T 3 ] .
Therefore, P 3 of Equation (11) does the map 3 3 .
Let ψ 1 , ψ 2 3 . Then, in the light of the conditions of Theorem 1, we derive that
P 3 ψ 1 ( t ) P 3 ψ 2 ( t ) = 1 Γ ( q ) 0 t ( t s ) q 1 i = 1 N G i s , ψ 1 ( s ) , ψ 1 ( θ i ( s ) ) G i s , ψ 2 ( s ) , ψ 2 ( θ i ( s ) ) d s
(since ψ 1 = ψ 2 = ζ 2 on ( 0 , T 2 ] , using (Hp3) and let t > T 2 )
1 Γ ( q ) T 2 t ( t s ) q 1 i = 1 N L G i ψ 1 ( s ) ψ 2 ( s ) + ψ 1 ( θ i ( s ) ) ψ 2 ( θ i ( s ) ) d s 2 i = 1 N L G i Γ ( q ) ψ 1 ψ 2 T 2 , T 3 T 2 T 3 ( T 3 s ) α 1 d s
(using (3) and (4) also let S = T 3 T 2 )
= 2 i = 1 N L G i Γ ( q ) ψ 1 ψ 2 T 2 , T 3 T 0 T 1 ( T 1 s ) q 1 d s = 2 i = 1 N L G i Γ ( q ) × ( T 1 T 0 ) q q ψ 1 ψ 2 T 0 , T 1 β ψ 1 ψ 2 ,
where ψ 1 ψ 2 T 2 , T 3 denotes the supremum on the given interval. Therefore, P 3 is a contraction on the set 3 with the unique fixed point ζ 3 , where ζ 3 agrees with ζ 2 on the interval ( 0 , T 2 ] since ζ 2 and ζ 3 are unique on the interval ( 0 , T 2 ] .
Step i. Finally, for the i t h step, by applying the mathematical induction method, we will be able to derive a unique solution on [ 0 , E ] . Showing this case is sufficient for full comprehension; however, the details on the induction are available in the following.
It is now clear that
T i 1 T i ( T i s ) q 1 d s = T 0 T 1 ( T 1 u ) q 1 d u Γ ( q ) β 2 i = 1 N L G i .
Let i , . be a complete metric space, which encompasses the functions ψ C ( 0 , T i ] , such that ψ ( t ) = ζ i 1 ( t ) of the previous step on the interval ( 0 , T i 1 ] together with the supremum norm.
We describe the transformation
P i : i i with ψ i .
This transformation implies that
P i ϑ ( t ) = ψ ( t ) , t [ ω , 0 ] ψ ( 0 ) + 1 Γ ( q ) 0 t ( t s ) q 1 i = 1 N G i ( s , ϑ ( s ) , ϑ ( θ i ( s ) ) ) d s , t ( 0 , T i ] .
For t [ ω , 0 ] , from Equation (12) we have
P i ψ 1 ( t ) P i ψ 2 ( t ) = 0 , ψ 1 , ψ 2 i .
Similar as previous steps, ζ i 1 ( t ) also satisfies Equation (12) on the interval ( 0 , T i 1 ] , i.e., for 0 < t T i 1 and in i , ψ ( t ) = ζ i 1 ( t ) , afterward we have
P i ψ ( t ) = ζ i 1 ( t ) = ψ ( 0 ) + 1 Γ ( q ) 0 t ( t s ) q 1 i = 1 N G i ( s , ζ i 1 ( s ) , ζ i 1 ( θ i ( s ) ) ) d s , t ( 0 , T i ] .
Therefore, P i of Equation (13) does a map i i .
Let ψ 1 , ψ 2 i . Consequently, according to Equation (12), it follows that
P i ψ 1 ( t ) P i ψ 2 ( t ) = 1 Γ ( q ) 0 t ( t s ) q 1 i = 1 N G i ( s , ψ 1 ( s ) , ψ 1 ( θ i ( s ) ) ) G i ( s , ψ 2 ( s ) , ψ 2 ( θ i ( s ) ) ) d s
(since ψ 1 = ψ 2 = ζ i 1 on ( 0 , T i 1 ] , using (Hp2) and let t > T i 1 )
1 Γ ( q ) T i 1 t ( t s ) q 1 i = 1 N L G i ψ 1 ( s ) ψ 2 ( s ) + ψ 1 ( θ i ( s ) ) ψ 2 ( θ i ( s ) ) d s 2 i = 1 N L G i Γ ( q ) ψ 1 ψ 2 T i 1 , T i T i 1 T i ( T i s ) q 1 d s
(using (3) and (4) also let S = T i T i 1 )
= 2 i = 1 N L G i Γ ( q ) ψ 1 ψ 2 T i 1 , T i T 0 T 1 ( T 1 s ) q 1 d s = 2 i = 1 N L G i Γ ( q ) × ( T 1 T 0 ) q q ψ 1 ψ 2 T 0 , T 1 β ψ 1 ψ 2 ,
where ψ 1 ψ 2 T i 1 , T i denotes the supremum on the given interval. Therefore, P i is a contraction on the set i with the unique fixed point ζ i , where ζ i agrees with ζ i 1 on the interval ( 0 , T i 1 ] since ζ i 1 and ζ i are unique on the interval ( 0 , T i 1 ] . Therefore, after completing the n-steps, we obtain a solution ζ n on the interval ( 0 , E ] . Additionally, to obtain a solution on the interval ( 0 , ) , we obtain a sequence of solutions { ζ i } on the intervals ( 0 , i ] and construct the set of the functions { ζ i * } , which are solutions from the earlier sequence, but continued for t > i at the constant value ζ i ( i ) . This sequence converges uniformly on compacts sets to a continuous solution ϑ ( t ) on the interval ( 0 , ) since at every value of t and for t > i , ϑ ( t ) aligns with ζ i ( t ) . This ends the proof of the theorem.  □

4. Numerical Applications of Theorem 1

In the section on numerical applications, we present two examples that satisfy the conditions of Theorem 1 in connection with two specific cases of the CFDE (1).
Example 1.
Think about the following CpFrDE:
D t 1 2 0 C ϑ ( t ) = 1 50 sin ( t ) + 1 100 ϑ ( t ) + 1 100 ϑ 2 ( t ) 1 + ϑ 2 ( t ) + 1 100 cos ϑ ( t 5 1 ) , t [ 0 , E ] ϑ ( t ) = t , t [ 5 1 , 0 ] .
We note that the CFDE (14) has the same form as the CFDE (1) with the data as follows:
q = 1 2 , [ ω , 0 ] = [ 5 1 , 0 ] , 0 t E , E > 0 ,
0 < θ 1 ( t ) = 5 1 = ω , ω , ω > 0 ,
G 1 ( t , ϑ , ϑ ( θ 1 ( t ) ) ) = 1 50 sin ( t ) + 1 100 ϑ ( t ) + 1 100 ϑ 2 ( t ) 1 + ϑ 2 ( t ) + 1 100 cos ϑ ( t 5 1 ) .
We will now demonstrate that the conditions (Hp1), (Hp2) and (Hp3) of Theorem 1 hold. It is now evident that the function G 1 is continuous. Therefore, the condition (Hp1) holds.
Let L G 1 = 1 100 , T 0 = 0 and T < 25 π . Then, we obtain
| G 1 ( t , ϑ 1 , ϑ 1 ( θ 1 ( t ) ) ) G 1 ( t , ϑ 2 , ϑ 2 ( θ 1 ( t ) ) ) | 1 100 ϑ 1 ϑ 2 + 1 100 ϑ 1 2 1 + ϑ 1 2 ϑ 2 2 1 + ϑ 2 2 + 1 100 cos ϑ 1 ( t 5 1 ) cos ϑ 2 ( t 5 1 ) = 1 100 ϑ 1 ϑ 2 + 1 100 ϑ 1 + ϑ 2 ϑ 1 ϑ 2 ( 1 + ϑ 1 2 ) ( 1 + ϑ 2 2 ) + 1 50 sin ϑ 1 ( t 5 1 ) + ϑ 2 ( t 5 1 ) 2 sin ϑ 1 ( t 5 1 ) ϑ 2 ( t 5 1 ) 2 = 1 100 ϑ 1 ϑ 2 + 1 100 ϑ 1 + ϑ 2 ϑ 1 ϑ 2 ( 1 + ϑ 1 2 ) ( 1 + ϑ 2 2 ) + 1 50 sin ϑ 1 ( t 5 1 ) + ϑ 2 ( t 5 1 ) 2 × sin ϑ 1 ( t 5 1 ) ϑ 2 ( t 5 1 ) 2 1 100 ϑ 1 ϑ 2 + 1 100 ϑ 1 ϑ 2 + 1 100 ϑ 1 ( t 5 1 ) ϑ 2 ( t 5 1 ) 1 50 ϑ 1 ϑ 2 + 1 100 ϑ 1 ( t 5 1 ) ϑ 2 ( t 5 1 )
and
2 i = 1 N L G i Γ ( q ) × T T 0 q q = 2 L G 1 Γ ( q ) × 2 T = 2 × 1 100 Γ 1 2 × 2 T = T 25 π = β < 1 .
Therefore, the conditions (Hp2) and (Hp3) of Theorem 1 are also met. Consequently, the CFDE (14) admits a unique solution.
Example 2.
Let us consider the CFDE as follows:
D t 1 3 0 C ϑ ( t ) = 1 200 exp ( t ) 1 1 + ϑ ( t ) + ϑ t 10 1 + 1 2 + ϑ ( t ) + ϑ t 20 1 sin ( t ) , t [ 0 , E ] , ϑ ( t ) = t , t [ 10 1 , 0 ] .
We see that the CFDE (15) and the CFDE (1) have the same form with the data as follows:
q = 1 3 , [ ω , 0 ] = [ 10 1 , 0 ] , 0 t E , E > 0 ,
0 < θ 1 ( t ) = 10 1 , 0 θ 2 ( t ) = 20 1 sin ( t ) 20 1 , ω = 10 1 , ω , ω > 0 ,
G 1 ( t , ϑ , ϑ ( θ 1 ( t ) ) ) = 1 200 exp ( t ) 1 + ϑ + ϑ t 10 1 ,
G 2 ( t , ϑ , ϑ ( θ 2 ( t ) ) ) = 1 200 exp ( t ) 2 + ϑ + ϑ t 20 1 sin ( t ) .
We will now show that the conditions (Hp1), (Hp2) and (Hp3) of Theorem 1 hold. Subsequently, it is clear that the functions G 1 and G 2 are continuous. Therefore, the condition (Hp1) holds.
Let
L G 1 = 1 200 , L G 2 = 1 200 , T 0 = 0 and T 3 < 89.298 .
Next, while considering the above data, we obtain
| G 1 ( t , ϑ 1 , ϑ 1 ( θ 1 ( t ) ) ) G 1 ( t , ϑ 2 , ϑ 2 ( θ 1 ( t ) ) ) | = 1 200 exp ( t ) 1 1 + ϑ 1 + ϑ 1 t 10 1 1 1 + ϑ 2 + ϑ 2 t 10 1 = 1 200 exp ( t ) 1 + ϑ 2 + ϑ 2 t 10 1 1 + ϑ 1 + ϑ 1 t 10 1 1 + ϑ 1 + ϑ 1 t 10 1 × 1 + ϑ 2 + ϑ 2 t 10 1 1 200 1 + ϑ 2 + ϑ 2 t 10 1 1 + ϑ 1 + ϑ 1 t 10 1 1 200 ϑ 1 ϑ 2 + 1 200 ϑ 1 ( t 10 1 ) ϑ 2 ( t 10 1 ) ,
| G 2 ( t , ϑ 1 , ϑ 1 ( θ 1 ( t ) ) ) G 2 ( t , ϑ 2 , ϑ 2 ( θ 1 ( t ) ) ) | = 1 200 exp ( t ) 1 2 + ϑ 1 + ϑ 1 t 20 1 sin ( t ) 1 2 + ϑ 2 + ϑ 2 t 20 1 sin ( t ) = 1 200 exp ( t ) 2 + ϑ 2 + ϑ 2 t 20 1 sin ( t ) 2 + ϑ 1 + ϑ 1 t 20 1 sin ( t ) 2 + ϑ 1 + ϑ 1 t 20 1 sin ( t ) × 2 + ϑ 2 + ϑ 2 t 20 1 sin ( t ) 1 200 2 + ϑ 2 + ϑ 2 t 20 1 sin ( t ) 2 + ϑ 1 + ϑ 1 t 20 1 sin ( t ) 1 200 ϑ 1 ϑ 2 + ϑ 1 t 20 1 sin ( t ) ϑ 2 t 20 1 sin ( t )
and
2 i = 1 N L G i Γ ( q ) × T T 0 q q = 2 L G 1 + L G 2 Γ ( q ) × 3 T 3 = 2 × 1 200 Γ 1 3 × 3 T 3 = 3 100 × 2.67894 × T 3 = 3 267.894 × T 3 = β < 1 ,
where Γ 1 3 2.67894 .
Therefore, the conditions (Hp2) and (Hp3) of Theorem 1 hold. Consequently, the CFDE (15) admits a unique solution.

5. The Outcomes and Contributions

We will now succinctly summarize the findings and contributions of this study.
-
In view of the literature review that is provided in this study’s introduction and available in the pertinent literature in order to investigate the existence and uniqueness of solutions for a variety of mathematical models, including IDEs, delay IEs, a scalar FDE of the Riemann–Liouville type, scalar nonlinear IEs, Hammerstein-type FIEs, nonlinear IEs including several variable time delays and nonlinear IDEs without delay on the intervals [ 0 , E ] and [ 0 , ) , the method of progressive contractions, which belongs to T.A. Burton, has effectively been applied and very interesting results have been obtained in the pertinent literature by this time. Nevertheless, to the best of our information from the relevant literature, there is no result on the same topic for CFDEs, including time delay(s) and without time delay(s), where the method of progressive contractions is used as a basic technique on the subject. Therefore, the aim of this study is to achieve the application of the progressive contractions technique to the existence and uniqueness of solutions for CFDEs, including several variable time delays. This fact represents the paper’s novelty and originality.
-
Burton [8], Theorem 1.2 dealt with a scalar FDE of the Riemann–Liouville type without delay. To the best of our knowledge, the fractional differential equation considered in Burton [8] has a simple form and does not include any time delay. Moreover, Burton [8] provided no numerical application to validate the primary outcome of [8], Theorem 1.2. Despite this case, in this paper, we examine a distinct kind of fractional differential equation, namely CFDE (1) with time delays in multiple variables, and also offer two new two examples as numerical applications of our primary finding, i.e., Theorem 1.
-
Regarding the benefits of the technique called progressive contractions, which was used in the proof of this paper, the existence and the uniqueness of solutions of different kinds of nonlinear mathematical models, as well as FDEs of various kinds, for example, FDEs of the Riemann–Liouville type with time delay and without time delay, CFDEs with time delay and without time delay, Hilfer type FDEs with time delay and without time delay, etc. can be discussed via the progressive contractions throughout very simple and short steps for the finite and infinite interval cases.

6. Conclusions

The uniqueness of solutions to a nonlinear CFDE with various variable time delays is the main problem of this paper. In the current work, we have deduced novel sufficient conditions for the uniqueness of solutions to the considered nonlinear CFDE with multiple variable time delays, which is called a qualitative characteristic. The fixed-point method, together with the supremum norm, is the basic tool in the proof of the main result of the paper, which is applied throughout the technique called progressive contractions.
To the best of our knowledge, Theorem 1 is the first result in the pertinent literature on the uniqueness of solutions for nonlinear delay CFDEs. The unique solution on the interval [ 0 , ) is obtained using the progressive contractions together supremum norm. The findings of this paper complement those that can be found in the pertinent literature and offer new insights into the qualitative theory of FDEs with and without time delays. In this study, it is a desirable case to provide Example 1 and Example 2 as numerical applications of Theorem 1.
Considering the present new result at hand, as future recommendations, the uniqueness of solutions to the following proper mathematical models can be investigated using the progressive contractions technique: CFIDEs with or without time delays, Riemann–Liouville FDEs with time delays, Riemann–Liouville IDEs with or without time delays, Hilfer FDEs with or without time delays, and Hilfer fractional IDEs with or without time delays.

Author Contributions

Conceptualization, C.T. and F.T.A.; Data curation, C.T. and F.T.A.; Formal analysis, C.T. and F.T.A.; Funding acquisition, F.T.A.; Methodology, C.T. and F.T.A.; Project administration, F.T.A.; Supervision, C.T. and F.T.A.; Validation, C.T. and F.T.A.; Visualization, C.T.; Writing-original draft, C.T. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported and funded by the Deanship of Scientific Research at Imam Mohammad Ibn Saud Islamic University (IMSIU) (grant number IMSIU-RP23069).

Data Availability Statement

The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding author.

Acknowledgments

The authors are very grateful to the anonymous referees and the handling editor whose comments and suggestions significantly contributed to the current presentation of this paper.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations will be used throughout this paper:
DEDifferential equation
DEsDifferential equations
DDEsDelay differential equations
CFDECaputo fractional-order differential equation
CFDEsCaputo fractional-order differential equations
CFIDEsCaputo fractional-order integro-differential equations
IEIntegral equation
IEsIntegral equations
IDEIntegro-differential equation
IDEsIntegro-differential equations
FDEFractional differential equation
FDEsFractional differential equations
FIEFunctional integral equation
FIEsFunctional integral equations
UHML stabilityUlam–Hyers–Mittag–Leffler stability

References

  1. Anderson, D.R.; Avery, R.I. Existence of a solution to a conjugate boundary value problem applying a corollary of the omitted ray fixed point theorem. Commun. Appl. Anal. 2015, 19, 403–412. [Google Scholar]
  2. Burton, T.A. Fixed points, differential equations, and proper mappings. In Proceedings of the International Conference on Nonlinear Operators, Differential Equations and Applications, Cluj-Napoca, Romania, 12–15 September 2001. [Google Scholar]
  3. Burton, T.A.; Zhang, B. Periodicity in delay equations by direct fixed point mapping. Differential Equations Dynam. Systems 1998, 6, 413–424. [Google Scholar]
  4. Mureşan, S.; Nica, O.M. Comparing two integral operators in the approach of Cauchy problems. An. Univ. Oradea Fasc. Mat. 2012, 19, 197–205. [Google Scholar]
  5. Burton, T.A. Existence and uniqueness results by progressive contractions for integro-differential equations. Nonlinear Dyn. Syst. Theory 2016, 16, 366–371. [Google Scholar]
  6. Burton, T.A.; Purnaras, I.K. Global existence and uniqueness of solutions of integral equations with delay: Progressive contractions. Electron. J. Qual. Theory Differ. Equ. 2017, 49. [Google Scholar]
  7. Burton, T.A. A note on existence and uniqueness for integral equations with sum of two operators: Progressive contractions. Fixed Point Theory 2019, 20, 107–111. [Google Scholar] [CrossRef]
  8. Burton, T.A. An existence theorem for a fractional differential equation using progressive contractions. J. Fract. Calc. Appl. 2017, 8, 168–172. [Google Scholar]
  9. Ilea, V.; Otrocol, D. On the Burton method of progressive contractions for Volterra integral equations. Fixed Point Theory 2020, 21, 585–593. [Google Scholar] [CrossRef]
  10. Tunç, O.; Tunç, C.; Yao, J.C. Global existence and uniqueness of solutions of integral equations with multiple variable delays and integro-differential equations: Progressive contractions. Mathematics 2024, 12, 171. [Google Scholar] [CrossRef]
  11. Tunç, O.; Tunç, C. Solution estimates to Caputo proportional fractional derivative delay integro-differential equations. Rev. Real Acad. Cienc. Exactas Fís. Nat. Ser. A Mat. RACSAM 2023, 117, 12. [Google Scholar] [CrossRef]
  12. Yousif, M.A.; Hamasalh, F.K. The fractional non-polynomial spline method: Precision and modeling improvements. Math. Comput. Simul. 2024, 218, 512–525. [Google Scholar] [CrossRef]
  13. Tunç, C.; Alshammari, F.S.; Akyıldız, F.T. Existence and uniqueness of solutions of Hammerstein type functional integral equations. Symmetry 2023, 15, 2205. [Google Scholar] [CrossRef]
  14. Graef, J.R.; Tunç, O.; Tunç, C. On Unique Solutions of Integral Equations by Progressive Contractions. Comput. Methods Differ. Equ. 2024, in press. [Google Scholar] [CrossRef]
  15. Abbas, S.; Benchohra, M.; Graef, J.R.; Henderson, J. Implicit Fractional Differential and Integral Equations. Existence and Stability; De Gruyter Series in Nonlinear Analysis and Applications, 26; De Gruyter: Berlin, Germany, 2018. [Google Scholar]
  16. Benchohra, M.; Karapınar, E.; Lazreg, J.E.; Salim, A. Advanced Topics in Fractional Differential Equations—A Fixed Point Approach; Synthesis Lectures on Mathematics and Statistics; Springer: Cham, Switzerland, 2023. [Google Scholar]
  17. Balachandran, K. An Introduction to Fractional Differential Equations; Industrial and Applied Mathematics; Springer: Singapore, 2023. [Google Scholar]
  18. Diethelm, K. The Analysis of Fractional Differential Equations. An Application-Oriented Exposition Using Differential Operators of Caputo Type; Lecture Notes in Mathematics; Springer: Berlin/Heidelberg, Germany, 2010. [Google Scholar]
  19. 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]
  20. Miller, K.S.; Ross, B. An Introduction to the Fractional Calculus and Fractional Differential Equations; A Wiley-Interscience Publication; John Wiley & Sons, Inc.: New York, NY, USA, 1993. [Google Scholar]
  21. Podlubny, I. Fractional Differential Equations. An Introduction to Fractional Derivatives, Fractional Differential Equations, to Methods of Their Solution and Some of Their Applications; Mathematics in Science and Engineering, 198; Academic Press, Inc.: San Diego, CA, USA, 1999. [Google Scholar]
  22. Zhou, Y. Basic Theory of Fractional Differential Equations, 3rd ed.; World Scientific Publishing Co. Pte. Ltd.: Hackensack, NJ, USA, 2024. [Google Scholar]
  23. Becker, L.C.; Burton, T.A.; Purnaras, I.K. Existence of solutions of nonlinear fractional differential equations of Riemann-Liouville type. J. Fract. Calc. Appl. 2016, 7, 20–39. [Google Scholar]
  24. Chauhan, H.V.S.; Singh, B.S.; Tunç, C.; Tunç, O. On the existence of solutions of non-linear 2D Volterra integral equations in a Banach space. Rev. Real Acad. Cienc. Exactas Fís. Nat. Ser. A Mat. RACSAM 2022, 116, 101. [Google Scholar] [CrossRef]
  25. Deep, A.; Deepmala, T.C. On the existence of solutions of some non-linear functional integral equations in Banach algebra with applications. Arab. J. Basic Appl. Sci. 2020, 27, 279–286. [Google Scholar] [CrossRef]
  26. Jung, S.M. A fixed point approach to the stability of a Volterra integral equation. Fixed Point Theory Appl. 2007, 2007, 57064. [Google Scholar] [CrossRef]
  27. Hale, J. Theory of Functional Differential Equations, 2nd ed.; Applied Mathematical Sciences; Springer: New York, NY, USA; Berlin/Heidelberg, Germany, 1977; Volume 3. [Google Scholar]
  28. Khan, H.; Tunç, C.; Chen, W.; Khan, A. Existence theorems and Hyers-Ulam stability for a class of Hybrid fractional differential equations with p-Laplacian operator. J. Appl. Anal. Comput. 2018, 8, 1211–1226. [Google Scholar]
  29. Khan, H.; Tunç, C.; Khan, A. Stability results and existence theorems for nonlinear delay-fractional differential equations with Φp-operator. J. Appl. Anal. Comput. 2020, 10, 584–597. [Google Scholar]
  30. Tunç, O. New Results on the Ulam–Hyers–Mittag–Leffler Stability of Caputo Fractional-Order Delay Differential Equations. Mathematics 2024, 12, 1342. [Google Scholar] [CrossRef]
  31. Zhang, T.; Li, Y. Global exponential stability of discrete-time almost automorphic Caputo–Fabrizio BAM fuzzy neural networks via exponential Euler technique. Knowl.-Based Syst. 2022, 246, 108675. [Google Scholar] [CrossRef]
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

Tunç, C.; Akyildiz, F.T. Unique Solutions for Caputo Fractional Differential Equations with Several Delays Using Progressive Contractions. Mathematics 2024, 12, 2799. https://doi.org/10.3390/math12182799

AMA Style

Tunç C, Akyildiz FT. Unique Solutions for Caputo Fractional Differential Equations with Several Delays Using Progressive Contractions. Mathematics. 2024; 12(18):2799. https://doi.org/10.3390/math12182799

Chicago/Turabian Style

Tunç, Cemil, and Fahir Talay Akyildiz. 2024. "Unique Solutions for Caputo Fractional Differential Equations with Several Delays Using Progressive Contractions" Mathematics 12, no. 18: 2799. https://doi.org/10.3390/math12182799

APA Style

Tunç, C., & Akyildiz, F. T. (2024). Unique Solutions for Caputo Fractional Differential Equations with Several Delays Using Progressive Contractions. Mathematics, 12(18), 2799. https://doi.org/10.3390/math12182799

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