1.
Introduction
Fractional differential equations (FDEs) have great interest for many mathematicians. This is due to extensive applications of these equations in the mathematical modeling in various fields of both science and engineering such as: control theory, physics, biological phenomena, viscoelasticity, and signal processing (see [1,2,3]). Furthermore, integro-differential equations are prevalent in various physical phenomena, such as fluid dynamics, biological models, and chemical kinetics. These equations arise due to the complex interactions and behaviors observed in these systems, requiring a combination of differential and integral terms to accurately model their dynamics. For instance, aero-elastic coupling in structures like wings and wind turbine blades leads to integro-differential problems, where control techniques play a crucial role in preventing instabilities. Overall, the presence of integro-differential equations in physical phenomena underscores the need for advanced mathematical tools to understand and predict the behavior of complex systems (see [4,5,6]). Recently, authors used various fixed-point theorems to prove the existence and uniqueness for the fractional differential equations with initial and boundary conditions. For example, the existence and uniqueness of solutions of differential equations with a mixture of integer and fractional derivatives have been investigated in [7]. The authors in [8] established existence and uniqueness results of solutions for fractional differential equations with integral boundary conditions by means of the Banach contraction mapping principle under sufficient conditions. The existence of solutions of integro-fractional differential equation when δ∈(2,3] through fixed-point theorem have been studied in [9]. Researchers in [10,11] study the existence and uniqueness of solutions for certain differential equations by using boundary and initial conditions, along with various techniques based on fixed-point theorems. The existence theory concerning fractional-order three-dimensional differential systems at resonance is presented in [12]; for additional details see these manuscripts [13,14,15].
On the other hand, the properties of Lp-solutions received a large share of researchers focus. Arshad et al. [16] examined Lp-solutions of fractional integral equations involving the Riemann-Liouville integral operator using a compactness condition. In [17] the author estimated the existence of an integrable solution for the nonlinear fractional differential equations involving two Caputo's fractional derivatives by means of Hölder's inequality together with Banach contraction principle and Schaefer's fixed-point theorem. Also see [18,19,20,21]. The Ulam-Hyers (UH) stability analysis has been studied and obtain a great part from the work of audiences [22,23]. Murad and Ameen in [24] researched the existence and UH stability of nonlinear fractional differential equations of mixed Caputo-Riemann derivatives. Vu et al. [25] proved the UH stability for the nonlinear Volterra integro-differential equations. Caputo-Hadamard (CH) fractional differential equations have various applications in modeling complex systems with memory effects in uncertain environments. Some potential applications include: Describing physical systems with memory effects and uncertain parameters, analyzing energy harvesting systems with fractional order properties, modeling biological systems with uncertain dynamics; see [26,27,28]. A series of research papers investigated the Hadamard derivative and CH derivative to prove the existence and stability theorems. In [29], existence and uniqueness of solution for Hadamard fractional differential equations on an infinite interval with integral boundary value has been developed. The theoretical analysis of CH fractional boundary-value problems in Lp-spaces was introduced in [30]. The authors in [31,32], focus on the existence and Ulam stability of solutions for certain CH fractional differential equations. The study in [33] highlights the existence of a solution for the boundary value problem of a nonlinear CH fractional differential equation with integral and anti-periodic conditions. Among the immense number of papers dealing with Caputo-Hadamard and Hadamard fractional differential equations subject to a variety of boundary conditions using fixed-point theory; we refer to [34,35,36]. Muthaiah et al. [37] discussed existence and of solutions for Hadamard fractional differential equations with integral boundary conditions. In [38] the authors applied the Monch's fixed-point theorem to prove the existence result for the fractional boundary value problems with CH derivative. Subsequently, many authors discussed the subject of approximation solutions by the Adomian decomposition method (ADM) for various types of FDE, we allude to [39,40,41]. Abdulahad et al. in [42] proved the existence of Lp-solutions for the following boundary value problem
Benhamida et al. [43] studied the existence of a solution for the boundary value problem:
where CHDδ1+ is the CH derivative, (0<δ≤1) and a,b,c are constants with
Wang et al. [44] employed the existence and uniqueness of positive solutions for the following integral boundary value problem:
where Dδ is the Riemann-Liouville fractional derivative.
In this paper, first we study the following nonlinear fractional differential equation with boundary conditions:
where CHDδ1+ {CH} derivative, with 1<δ≤2 and
is continuous function.
Second, the following fractional integro-differential equations with boundary conditions are investigated:
where CHDδ1+ is the CH derivative and 0<δ≤1. Here λ1–λ3 are constants
with
are continuous functions.
This paper is arranged as follows: In Section 2, we present some preliminary results to be used later. In Section 3, the Krasnoselskii's fixed-point theorem and Banach contraction principle are applied to analyze the existence and uniqueness of solutions to the problems (1.1)–(1.4) in Lp-spaces. Moreover, we discuss the UH and UH-Rassias stability for the problems in Section 4. The ADM is implemented to find the approximate solutions for the given problems in Section 5. Finally, examples are also given to show the applicability of our results.
2.
Preliminaries
Let us give some definitions and lemmas that are basic and needed at various places in this work.
Definition 2.1. [45] The Hadamard fractional integral of order δ∈R for a continuous function V is defined as
provided the integral exists.
Definition 2.2. [45] The Hadamard derivative of fractional order δ∈R for a continuous function f is defined as
where
where [δ] denotes the integer part of the real number δ.
Definition 2.3. [45] The CH derivative of fractional order δ∈R for a continuous function V is defined as follows:
where
and [δ] denotes the integer part of the real number δ, and Γ is the gamma function.
Lemma 2.1. [45] Let
If ϕ∈ACnδ[a,b], then the differential equation
has solutions
and the following formula holds:
where ck∈R,k=1,2,⋯,n−1.
Definition 2.4. [46] The Eq (1.1) is UH stable if there exists a real number cf>0 such that for each ε>0 and for each solution z∈C1([a,b],R) of the inequality
there exists a solution ϕ∈C1([a,b],R) of Eq (1.1) with
Theorem 2.1. [47] (Krasnoselskii fixed-point theorem)
Let H be a closed, bounded, convex, and nonempty subset of a Banach space V. Let A and B be two operators such that
(1)Az1+Bz2∈H whenever z1,z2∈H;
(2) A is compact and continuous;
(3) B is a contraction mapping.
Then there exists z∈H such that
Lemma 2.2. [48] (Bochner integrable)
A measurable function
is Bochner integrable, if ||V|| is Lebesgue integrable.
Theorem 2.2. [49] (Kolmogorov compactness criterion)
Let ν⊆Lp[a,b], 1≤p<∞. If:
(i) ν is bounded in Lp[a,b];
(ii) xh→xash→0 uniformly with respect to x∈ν, then ν is relatively compact in Lp[a,b], where
Lemma 2.3. [50] (Hölder's inequality)
Let X be a measurable space, let p and q satisfy
and
If h∈Lp(X) and g∈Lq(X), then (hg) belongs to L(X) and satisfies
Lemma 2.4. [36] If
then
3.
Main results
3.1. Existence and uniqueness results for problems (1.1) and (1.2)
This section deals with the existence and uniqueness of a solution for the fractional differential Eq (1.1) with boundary condition (1.2). For measurable functions
define the norm
where Lp(J,R) is the Banach space of all Lebesgue measurable functions. Now, consider the following assumptions:
(F1) There exists a constant μ>0 such that
for each t∈J and for all ϕ∈R.
(F2) V(t,ϕ) is continuous and satisfies the Lipschitz condition, there exists a constant ω1>0 such that
for each ϕ1,ϕ2∈R.
For the sake of convenience, we set the notation:
Lemma 3.1. For any
then the boundary value problems (1.1) and (1.2) have a solution
Proof. Applying Lemma 2.1, we can reduce the problems (1.1) and (1.2) to an equivalent integral equation
to find c0 and c1, from the first boundary condition
we obtain
by using the condition
the result is
and
by using Fubini's theorem, the following is obtain
Hence
this implies that
This completes the poof.
□
The first result is based on Banach contraction principle.
Theorem 3.1. Assume that (F1) and (F2) hold. If
Then the boundary value problems (1.1) and (1.2) have a unique solution.
Proof. Define the operator ⊤ by
we need to prove that the operator ⊤ has a fixed-point on the set
For ϕ∈ϑu, we have
By Hölder's inequality and Lemma 2.4 we obtain
Now, by the same way we find that
and
Thus, Eqs (3.3)–(3.5) are Lebesgue integrable; by using Lemma 2.2, we conclude that (lnes)δ−1V(s,ϕ(s)), (lnes)δV(s,ϕ(s)), and (lnts)δ−1V(s,ϕ(s)) are Bochner integrable with respect to s∈[1,t]; for all t∈J, then the Eq (3.2) becomes
Then, by the condition (F1) implies that
Integrate by parts; the following is obtained:
and
which implies that ⊤ϑu⊆ϑu.
Hence, ⊤(ϕ)(t) is Lebesgue integrable and ⊤ maps ϑu into itself. We have to show that ⊤ is a contraction mapping. Let ϕ1,ϕ2∈Lp(J), we have
Using (F2) and H¨older's inequality, one has
Integrate by parts, leads to
and
If
then by the contraction mapping principle, the boundary value problems (1.1) and (1.2) have a unique solution. □
The following result is based on Krasnoselskii's fixed-point theorem.
Theorem 3.2. Assume that (F1) and (F2) hold. Then the boundary value problems (1.1) and (1.2) have at least one solution.
Proof. Let us define two operators, χ1 and χ2, from Eq (3.1) as
Consider the set
For x,y∈ðr, we have
By (F1) and Hölders inequality, Eq (3.6) becomes
It follows from integration by parts, that
and
Hence,
Now, to prove that χ2 is a contraction mapping on ðr, from (F2) and H¨older inequality, it is easy to see that
and
If ℵ2ω1<1, then χ2 is a contraction mapping.
We need to show that χ1 is compact and continuous, for any x∈ðr, we have
Hence, χ1 is uniformly bounded. To show that χ1 is completely continuous, we apply Theorem 2.2, the Kolmogorov compactness criterion. Let Ω be a bounded subset of ðr. Then χ1(Ω) is bounded in Lp(J), the condition (i) of Theorem 2.2 is applied. Next we will show that (χ1x)h→χ1x in Lp(J) as h→0, uniformly with respect to x∈Ω. We have the following estimation:
Since V∈Lp(J), we get that IδV∈Lp(J)
Hence
uniformly a h→0. Then, by Theorem 2.2, we deduce that χ1(Ω) is relatively compact; that is, χ1 is a compact operator. As a consequence of Krasnoselskiis fixed-point theorem, the boundary value problems (1.1) and (1.2) have at least one solution in ðr. □
3.2. Existence of solution for the fractional integro-differential equation
In this section, we prove the existence and uniqueness of solutions through Krasnoselskiis and Banach fixed-point theorems for the integro-differential equations
with the boundary condition
For measurable functions
define the norm
where Lp(I,R) is the Banach space of all Lebesgue measurable functions. In order to achieve the results, the following assumptions are required:
(P1) There exists positive constants η1 and η2 such that
and
for each t∈I and all ϕ∈R.
(P2) There exists a positive constants ϱ1,ϱ2>0, such that
for each ϕ1,ϕ2∈R.
For computational convenience, we set
Lemma 3.2. Let ϕ∈C(I,R) and 0<δ≤1, then the solution of the boundary value problems (3.7) and (3.8) is given by
Proof. By applying Lemma 2.1, we can reduce the problems (3.7) and (3.8) to an integral equation
from the boundary condition (3.8), we obtain
Then the solution is
The proof is completed.
□
Our first result is based on krasnoselskii's fixed-point theorem.
Theorem 3.3. Assume that (P1) and (P2) hold. Then the boundary value problems (3.7) and (3.8) have at least one solution.
Proof. Let us define two operators, ϖ1 and ϖ, from Eq (3.9) as
Consider the set
For x,ϕ∈δr, we have
By using (P1) and H¨older's inequality, the below is found:
Integrate by parts, the result are
and
Hence, ϖ1x+ϖϕ∈δr.
Now, to prove that ϖ is a contraction in Lp(I). Letting ϕ1,ϕ2∈Lp(I), we have
Then from (P2) and H¨older's inequality, Eq (3.10) becomes
and
If ℑ3<1, then ϖ is a contraction mapping principle. Moreover, continuity of x(t) implies that the operator ϖ1x is continuous
Hence, ϖ1 is uniformly bounded on δr.
Next to show that ϖ1 is completely continuous, we apply Theorem 2.2, the Kolmogorov compactness criterion. Let ζ be a bounded subset of δr. Then ϖ1(ζ) is bounded in Lp(I) and the condition (i) of Theorem 2.2 is applied. Next, to show that (ϖ1x)h→ϖ1x in Lp(I) as h→0, uniformly with respect to x∈ζ. Let
and
Since V∈Lp(I), we get that IδV∈Lp(I), so we have
Then by Theorem 2.2, we deduce that ϖ1(ζ) is relatively compact; this implies that ϖ1 is a compact operator. As a consequence of Krasnoselskiis fixed-point theorem the boundary value problems (3.7) and (3.8) have at least one solution. The proof is complete. □
Now, the uniqueness result for the problems (3.7) and (3.8) is based on the Banach contraction principle.
Theorem 3.4. Suppose that (P1) and (P2) holds. If
Then the boundary value problems (3.7) and (3.8) have a unique solution.
Proof. Define the operator
as follows:
We need to show that θσr⊆σr, where,
For ϕ∈σr, we have
It follows from condition (P1) and H¨older's inequality that
Hence, θ maps Lp(I) into itself. Now, to prove that θ is a contraction mapping. Let ϕ1,ϕ2∈Lp(I), we get
By repeating the same technique of Theorem 3.3, it immediately follows that
and
If
then θ is a contraction mapping. Therefore, by using Banach contraction mapping, θ has a unique fixed point, which is a unique solution of the boundary value problems (3.7) and (3.8). □
4.
UH stability
In this section, we will study the analysis of UH stability of the fractional differential Eq (1.1) with boundary condition (1.2) and for the problems (3.7) and (3.8).
4.1. UH stability for problems (1.1) and (1.2)
Theorem 4.1. If the hypothesis (F2) holds with
Then the boundary value problems (1.1) and (1.2) are UH stable.
Proof. For ϵ>0 and w be a solution that satisfies the following inequality
there exists a solution ϕ∈Lp(J) of the boundary value problems (1.1) and (1.2). Then ϕ(t) is given by
From the inequality (4.1) and for each t∈J, we have
On the other hand, for each t∈J, the below is found
Thus, by condition (F2) and Hölder inequality which implies that
Integrating by parts, we have
Hence
where
which implies that the boundary value problems (1.1) and (1.2) have UH stability. □
4.2. UH stability for problems (3.7) and (3.8)
Theorem 4.2. If the hypothesis (P2) holds with ℧3<1. Then the boundary value problems (3.7) and (3.8) are UH stable.
Proof. For ϵ>0 and each solution w∈Lp(I) of the inequality
and there exists a solution ϕ∈Lp(I) of the boundary value problems (3.7) and (3.8). Then ϕ(t) is given by
From the inequality (4.2) and for each t∈I, we obtain
for each t∈I, the below is found
by (P2), for each t∈I, we obtain
It follows from Hölder's inequality and integration by parts that
and
Hence
where
Thus, the solution of (3.7) and (3.8) is UH stable. □
5.
Examples
In this section, some examples are given to illustrate our main results.
Example 5.1. Consider the fractional boundary differential equation
where
and μ=0.0198, from condition (F2), we get
To prove the existence of a solution, Theorem 3.2 is applied as follows:
After taken r=1, one can has:
If p=2, then (℧1)1p(0.0198)r=0.4489.
If p=3, then (℧1)1p(0.0198)r=0.3485.
If p=4, then (℧1)1p(0.0198)r=0.3137.
The second step shows that χ2 is a contraction mapping
If p=2, then (ℵ2)12=11.0867,(ℵ2)12ω1=0.2194.
If p=3, then (ℵ2)13=8.7781,(ℵ2)13ω1=0.1737.
If p=4, then (ℵ2)14=7.9218,(ℵ2)14ω1=0.1568.
Hence, χ2 is a contraction mapping.
The third step shows that χ1 is compact and continuous, one can has
If p=2, then ||χ1x||p<0.5370∗0.0198=0.0106.
If p=3, then ||χ1x||p<0.3115∗0.0198=0.0062.
If p=4, then ||χ1x||p<0.2000∗0.0198=0.0040.
Hence, χ1 is uniformly bounded and relatively compact. All steps of Theorem 3.2 are satisfied; therefore, we deduce that the problem has at least one solution.
Next, to explain the uniqueness of the solution, and according to Theorem (3.1), the results are:
If p=2, then (℧1)12=22.2253, (℧1)12ω1=0.43988<1.
If p=3, then (℧1)13=17.5586, (℧1)13ω1=0.34752<1.
If p=4, then (℧1)14=15.8437, (℧1)14ω1=0.31357<1.
Then, the problem (5.1) has a unique solution.
Example 5.2. Consider the following boundary value problem:
where
and
by using the condition (F2), one has ω1=0.05011255. Moreover, from Theorem 3.1, we see that:
If p=2, then
If p=3, then
If p=4, then
By Theorem 3.1, the problem (5.2) has a unique solution.
Example 5.3. Consider the boundary value problem
where
and using the Lipshitz condition (P2), the out comes are
To estimate the problem has at least one solution, apply Theorem 3.3. For the first step, we have
Let
to get:
If p=2, ℑ4=0.9302<r.
If p=3, ℑ4=0.8025<r.
If p=4, ℑ4=0.7637<r.
Hence,
For the second step,
to prove that ϖ is a contraction mapping, the results are:
If p=2, ℑ3=0.0699<1.
If p=3, ℑ3=0.0473<1.
If p=4, ℑ3=0.0407<1.
Hence, ϖ is a contraction mapping.
For the third step, show that ϖ1 is compact, we have
If p=2 then, ||(ϖ1x)(t)||p<0.0762.
If p=3 then, ||(ϖ1x)(t)||p<0.0501.
If p=4 then, ||(ϖ1x)(t)||p<0.0437.
Hence, ϖ1 is uniformly bounded and relatively compact. All conditions of Krasnoselskiis fixed-point theorem are satisfied, then the problem has at least one solution.
Now, To exhibit there is only one solution, the Banach fixed-point Theorem 3.4 is applied as follows
If p=2 then (℧3)12=0.20659<1.
If p=3 then (℧3)13=0.12266<1.
If p=4 then (℧3)14=0.10042<1.
Then the problem (5.3) has a unique solution.
Example 5.4. Consider the fractional boundary value problem
where δ=0.6, λ1=4,λ2=2,λ3=0.5, and by using the condition (P2), we have ϱ1=0.129198 and ϱ1=0.4, from Theorem 3.4, the results are:
If p=2, then (℧3)12=0.403497<1.
If p=3, then (℧3)13=0.169727<1.
If p=4, then (℧3)14=0.132042<1.
Then by Theorem 3.4, the problem (5.4) has a unique solution.
6.
Numerical results
In this section, we deal with ADM to find the approximate solution of fractional differential and integro-differential equations; some numerical examples are presented to compare between the exact and approximate solutions.
6.1. The ADM
George Adomian established the ADM in the 1980s. The ADM has been paid much attention in the recent years in applied mathematics, and in the field of series solutions particular. Moreover, it is a fact that this method is powerful and effective, and it easily solves many types of linear or nonlinear ordinary or partial differential equations, and integral equations; see [51,52]. This method generates a solution in the form of a series whose terms are determined by a recursive relationship using these Adomian polynomials. A brief outline of the method follows. For every nonlinear differential equation, it can be decomposed into the following form:
where L is the highest order differential operator, R(ϕ) is the remainder of the linear part, N(ϕ) represents the nonlinear part and h is a given function. In general, the operator L is invertible. If we take L−1(integral operator) on both sides of Eq (6.1), an equivalent expression can be given
here g satisfies Lg=0 and the initial conditions. If L is the second-order derivative, L−1 is the two-fold definite integral. For the ADM, the solution u is expressed in terms of a series form
If we have a nonlinear term N(ϕ) it is represented by the Adomian polynomials Ak
Ak depends on ϕ0,ϕ1,⋯,ϕk and can be formulated by
Then Eq (6.2) can be written as
Example 6.1. Consider the boundary value problem
The exact solution is
Applying the inverse operator
on (6.3), we find that:
In order to obtain ϕ(t), we apply the Adomian iterative scheme
Now, to find, ϕ1,ϕ2,ϕ3,⋯, it follows
The approximate solution of problem (6.3) is:
Tables 1-3 show the approximate and exact solutions for Example 6.3.
One can observe on Figure 1 a decrease of the calculated errors towards zero, which confirms the results of convergences of the approximate solution to the exact solution in Example 6.1 when t∈[1,e] and for various fractional orders δ=1.2,1.3,⋯,1.9.
Example 6.2. Consider the boundary value problem
The exact solution is
Now applying the inverse operator
on (6.4), one obtains
where
By the same technique, we find
The solution is the given by
Tables 4-6, approximate and exact solutions for Eq (6.4).
The behavior of the absolute errors approaching zero, as displayed in Figure 2, supports the convergence results of the approximate solution toward the exact solution referenced in Example 6.2.
Example 6.3. Consider the fractional integro-differential equation with initial condition
The exact solution is
Applying the inverse operator
on (6.5), we find that
where
and
The analogous process gives
Tables 7-9 show the approximate and exact solutions for Eq (6.5).
The Figure 3 shows that the absolute errors between the approximate solution and the exact solution described in Example 6.3 are approaching zero.
This indicates that the ADM is a powerful and effective technique for obtaining accurate results.
7.
Conclusions
In this paper, we investigated the Lp-solutions for nonlinear fractional differential and integro-differential equations with boundary conditions in the sense of the CH derivative. By means of the Krasnoselskii fixed-point theorem and the Banach contraction principle, we have established sufficient conditions for the existence and uniqueness of solutions for a nonlinear problems. In addition, the UH stability of the solutions for the indicated problems is studied. We also employed the ADM to estimate the approximate solutions. Finally, we present examples to demonstrate the consistency of the theoretical findings. In future works, one can extend the given problems to more fractional derivatives such as the Hilfer derivative and Caputo-Fabrizio derivative.
Use of Generative-AI tools declaration
The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.
Conflict of interest
The authors declare that they have no conflicts of interest in this paper.