Next Article in Journal
Sugeno Integral Based on Overlap Function and Its Application to Fuzzy Quantifiers and Multi-Attribute Decision-Making
Next Article in Special Issue
Asynchronous Switching Control of Discrete Time Delay Linear Switched Systems Based on MDADT
Previous Article in Journal
An Equilibrium Strategy for Target Benefit Pension Plans with a Longevity Trend and Partial Information
Previous Article in Special Issue
Weighted Pseudo-θ-Almost Periodic Sequence and Finite-Time Guaranteed Cost Control for Discrete-Space and Discrete-Time Stochastic Genetic Regulatory Networks with Time Delays
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Solvability, Approximation and Stability of Periodic Boundary Value Problem for a Nonlinear Hadamard Fractional Differential Equation with p-Laplacian

Department of Mathematics, School of Electronics & Information Engineering, Taizhou University, Taizhou 318000, China
Submission received: 23 June 2023 / Revised: 23 July 2023 / Accepted: 26 July 2023 / Published: 27 July 2023
(This article belongs to the Special Issue Differential Equations in Applied Mathematics)

Abstract

:
The fractional order p -Laplacian differential equation model is a powerful tool for describing turbulent problems in porous viscoelastic media. The study of such models helps to reveal the dynamic behavior of turbulence. Therefore, this article is mainly concerned with the periodic boundary value problem (BVP) for a class of nonlinear Hadamard fractional differential equation with p -Laplacian operator. By virtue of an important fixed point theorem on a complete metric space with two distances, we study the solvability and approximation of this BVP. Based on nonlinear analysis methods, we further discuss the generalized Ulam-Hyers (GUH) stability of this problem. Eventually, we supply two example and simulations to verify the correctness and availability of our main results. Compared to many previous studies, our approach enables the solution of the system to exist in metric space rather than normed space. In summary, we obtain some sufficient conditions for the existence, uniqueness, and stability of solutions in the metric space.

1. Introduction

The p -Laplacian differential equation is one of the famous and important second-order nonlinear ordinary differential equations (ODEs). This equation first appeared in Leibenson’s study [1] of turbulence in porous media in 1983. The underlying form of p -Laplacian differential equation is written as
Φ p ( u ( t ) ) = f ( t , u ( t ) ) , t ( 0 , 1 ) ,
where Φ p : x | x | p 2 x ( p > 1 ) is called the p -Laplacian operator. Its inverse is Φ p 1 = Φ q with 1 p + 1 q = 1 . Due to its description of fundamental mechanical problems in turbulence, p -Laplacian differential equations have been extensively and deeply studied. In recent years, some scholars have begun to focus on the nonlinear fractional differential system with p -Laplacian. For example, the authors in [2] investigated the multiple positive solutions of a nonlinear high order Riemann-Liouville fractional p -Laplacian equation with integral boundary value conditions. In [3], the author explored the existence and GUH-stability of a nonlinear Caputo-Fabrizio fractional coupled Laplacian equations. In [4], based on the Guo-Krasnosel’skii fixed point theorem, the authors probed into the multiple positive solutions of a system of mixed Hadamard fractional BVP with ( p 1 , p 2 ) -Laplacian. In fact, some articles have been disposed of the BVP of p -Laplacian system involving Riemann-Liouville or Caputo fractional derivatives (see [5,6,7,8,9,10,11,12,13]).
Hadamard [14] raised a novel fractional integral and derivative in 1892, which was later named Hadamard-type fractional calculus. There are some obvious differences between Hadamard fractional calculus and Riemann-Liouville fractional calculus. For example, k H ( t , s ) = ( log t s ) γ 1 is the integral kernel corresponding to γ -order Hadamard fractional derivative, while k R ( t , s ) = ( t s ) γ 1 is the integral kernel corresponding to γ -order Riemann-Liouville fractional derivative. Furthermore, for any λ > 0 , k H ( λ t , λ s ) = k H ( t , s ) is different from k R ( λ t , λ s ) = λ γ 1 k R ( t , s ) k R ( t , s ) . The study on Hadamard fractional differential equations has attracted the attention of many scholars. There have been a series of fruitful achievements (see [15,16,17,18,19,20,21]). In 1940s, Ulam and Hyers [22,23] put forward a new stability that describes the stationarity of the exact and approximate solutions of system. Subsequently, extensive and in-depth research was conducted on the Ulam-Hyers stability of various systems. Especially, many excellent research results have emerged regarding the Ulam-Hyers stability of fractional order differential systems (see some of them [3,21,24,25,26,27,28,29,30,31,32]). Moreover, it is rare to combine the Hadamard fractional derivative with Laplacian operator. Therefore, it is novel and interesting to probe these problems.
Illuminated by the above arguments, this manuscript deals with the periodic BVP of a nonlinear Hadamard fractional differential equation with p -Laplacian operator as follows:
H D 1 + α Φ p ( H D 1 + β u ( t ) ) = f ( t , u ( t ) ) , t ( 1 , e ] , u ( 1 ) = u ( e ) , H D 1 + β u ( 1 ) = H D 1 + β u ( e ) ,
where 0 < α 1 , 1 < β 2 , p > 1 , H D 1 + * is the *-order Hadamard fractional derivative, Φ p = | x | p 2 x , and its inverse Φ p 1 = Φ q with 1 p + 1 q = 1 , f C ( [ 1 , e ] × R , R ) . In addition, our study has also been inspired by the latest achievements in fractional differential equations, such as numerical algorithms and simulations [33,34,35,36,37,38], as well as the application of some nonsingular fractional derivative models [34,37,38,39,40,41,42,43].
The paper aims to discuss the approximation and GUH-stability of BVP (1). The novelty of this paper is mainly reflected as follows: (a) Since there is no paper dealing with the approximation problem of nonlinear Hadamard fractional differential systems with Laplace operator, we first consider the system (1) to fill this gap. (b) By applying a fixed point theorem on complete metric space with two kinds of distance, we obtain some sufficient conditions to ensure that system (1) has a unique solution. In addition, we build the generalized Ulam-Hyers stability of system (1) based on nonlinear analysis methods and inequality techniques. (c) Many previous papers (see [2,3,4,5,6,7,8,9,10,11,12,13,17,18,19,24,25]) usually used some fixed-point theorems on Banach spaces to study the existence of solutions of fractional differential equations. However, we handle the existence of solutions to fractional order differential equations by defining two different distances on a complete distance space. This allows for the discussion of the existence of solutions in a broader space, and there are relatively few restrictions on the existence of solutions. Therefore, our research methods and results are novel and interesting.
The rest sections of this paper are organized as follows. In Section 2, we recollect the definition of Hadamard fractional integrals and derivatives and some necessary lemmas. In Section 3, we discuss the existence, uniqueness, and approximation of solutions to BVP (1) by constructing two different distances and applying an important fixed point theorem on metric space. Furthermore, we use nonlinear analysis methods and inequality techniques to establish the GUH-stability of BVP (1) in Section 4. Section 5 provides the numerical solutions and simulations for two examples by means of ODE113 toolbox in MATLAB. Finally, we have made a brief summary in Section 6.

2. Preliminaries

This portion mainly introduces some important concepts and lemmas.
Definition 1 
([44]). For a > 0 , the left-sided Hadamard fractional integral of order γ > 0 for a function ξ : [ a , ) R is defined by
H J a + γ ξ ( t ) = 1 Γ ( γ ) a t log t s γ 1 ξ ( s ) d s s ,
provided the integral exists, where Γ ( γ ) = 0 t γ 1 e t d t and log ( · ) = log e ( · ) .
Definition 2 
([44]). Let a , γ > 0 and ξ C m [ a , ) , the γ-order left-sided Hadamard fractional derivative is defined by
H D a + α ξ ( t ) = 1 Γ ( m γ ) t d d t m a t log t s m γ 1 ξ ( s ) d s s ,
where m 1 < α m , m = [ γ ] + 1 , and [ · ] is the Gaussian truncating integer function.
Lemma 1 
([44]). Let a , b , γ > 0 and ξ C m ( a , b ) L 1 ( a , b ) , then
H J a + γ ( H D a + γ ξ ( t ) ) = ξ ( t ) + i = 1 m c i log t a γ i ,
where c 1 , c 2 , , c m are some real constants, and m = [ γ ] + 1 .
Lemma 2. 
Let p > 1 . The p -Laplacian operator Φ p ( z ) = | z | p 2 z has the followings:
( i )
If z 0 , then Φ p ( z ) = z p 1 , and Φ p ( z ) is increasing with respect to z ;
( ii )
For all z , w R , Φ p ( z w ) = Φ p ( z ) Φ p ( w ) ;
( iii )
If 1 p + 1 q = 1 , then Φ q [ Φ p ( z ) ] = Φ p [ Φ q ( z ) ] = z , for all z R ;
( iv )
For all z , w 0 , z w Φ q ( z ) Φ q ( w ) ;
( v )
0 z Φ q 1 ( w ) 0 Φ q ( z ) w ;
( vi )
| Φ q ( z ) Φ q ( w ) | ( q 1 ) M ¯ q 2 | z w | , q 2 , 0 z , w M ¯ , ( q 1 ) M ̲ q 2 | z w | , 1 < q < 2 , z , w M ̲ 0 .
Now we introduce the following important fixed point theorem on a complete metric space involving two different distances, which will be used to prove the existence and uniqueness of solution to BVP (1).
Lemma 3 
([45]). Let ρ and ϱ be two different metrics on a nonempty set X , and define an operator T : X X . Assume that
( a 1 )
For all x , y X , there has a constant ι > 0 such that ϱ ( T x , T y ) ι ρ ( x , y ) ;
( a 2 )
( X , ϱ ) is a complete metric space;
( a 3 )
T : ( X , ϱ ) ( X , ϱ ) is continuous;
( a 4 )
For all x , y X , there has a constant 0 < κ < 1 such that ρ ( T x , T y ) κ ρ ( x , y ) .
Then there has a unique x * X such that T x * = x * , and lim k T k x 0 = x * for any x 0 X .
It is worth noting that the application techniques and related generalization of Lemma 3 can also be found in [46,47,48,49] and the references therein.

3. Solvability and Approximation

In this portion, we will prove the existence of a unique solution for system (1) based on Lemma 3. To this end, we need the following important lemma.
Lemma 4. 
Assume that 0 < α 1 , 1 < β 2 and p > 1 are some constants, f C ( [ 1 , e ] × R , R ) . Then BVP (1) is equivalent to the following integral equation
u ( t ) = A u ( e ) ( log t ) β 1 + A u ( t ) ,
where A u ( t ) = H J 1 + β Φ q H J 1 + α f ( t , u ( t ) ) .
Proof. 
If u ( t ) C ( ( 1 , e ] , R ) is a solution of system (1), then it follows from Lemma 1 that
Φ p ( H D 1 + β u ( t ) ) = c ( log t ) α 1 + H J 1 + α f ( t , u ( t ) ) ,
which implies that
H D 1 + β u ( t ) = Φ q [ c ( log t ) α 1 + H J 1 + α f ( t , u ( t ) ) ] .
According to the existence of H D 1 + β u ( 1 ) and (4), we know that c = 0 and H D 1 + β u ( 1 ) = 0 . By Lemma 1 and (4), we have
u ( t ) = d 1 ( log t ) β 1 + d 2 ( log t ) β 2 + H J 1 + β Φ q H J 1 + α f ( t , u ( t ) ) .
Similarly, we drive from u ( 1 ) = u ( e ) and (5) that d 2 = 0 and
d 1 = H J 1 + β Φ q H J 1 + α f ( t , u ( t ) ) | t = e = A u ( e ) .
In view of (5) and (6), we have
u ( t ) = A u ( e ) ( log t ) β 1 + A u ( t ) .
Thus, u ( t ) C ( ( 1 , e ] , R ) is a solution of system (2). Vice versa, if u ( t ) C ( ( 1 , e ] , R ) is a solution of (2), then it is also a solution of (1) because the above derivation is completely reversible. The proof is completed. □
Let X = C ( [ 1 , e ] , R ) , two different distances ρ , ϱ : X X are respectively defined by
ρ ( u ( t ) , v ( t ) ) = sup t [ 1 , e ] | u ( t ) v ( t ) | , ϱ ( u ( t ) , v ( t ) ) = 1 e | u ( t ) v ( t ) | d t ,
for all u ( t ) , v ( t ) X . It is easy to prove that ( X , ρ ) and ( X , ϱ ) are all complete metric spaces. In addition, we need the following underlying assumptions in the whole paper.
( H 1 )
0 < α 1 , 1 < β 2 and 1 < p 2 are some constants, f C ( [ 1 , e ] × R , R ) .
( H 2 )
There has a constant M > 0 such that
0 f ( t , u ) M , t [ 1 , e ] , u R .
( H 3 )
There has a function 0 l ( t ) C [ 1 , e ] such that, for all t [ 1 , e ] and u , v R ,
| f ( t , u ) f ( t , v ) | l ( t ) | u v | .
( H 4 )
κ = 2 ( q 1 ) Γ ( α + β + 1 ) M Γ ( α + 1 ) q 2 l e < 1 , where l e = max 1 t e { l ( t ) } .
Theorem 1. 
If ( H 1 ) ( H 4 ) are fulfilled, then BVP (1) has a unique solution u * ( t ) X .
Proof. 
In what follows, we will apply Lemma 3 to prove Theorem 1. Two different distances ρ , ϱ : X X are defined as (8), then ( X , ρ ) and ( X , ϱ ) are all complete metric spaces, which indicates that the condition ( a 2 )  holds. According to Lemma 4, for all u ( t ) X , an operator T : X X is defined by
T ( u ( t ) ) = A u ( e ) ( log t ) β 1 + A u ( t ) ,
where A u ( t ) = H J 1 + β Φ q H J 1 + α f ( t , u ( t ) ) .
From the continuity of Φ q and Hadamard fractional integral, we know that T : ( X , ϱ ) ( X , ϱ ) is continuous, which means that the condition ( a 3 )  holds. By ( H 2 ) , we have
0 H J 1 + α f ( t , u ( t ) , u ( t τ ( t ) ) ) M Γ ( α + 1 ) ( log t ) α M Γ ( α + 1 ) , t [ 1 , e ] .
For all u , v X , t [ 1 , e ] , we derive from ( vi ) in Lemma 2, ( H 3 ) and (10) that
| A u ( t ) A v ( t ) | = | H J 1 + β Φ q H J 1 + α f ( t , u ( t ) ) Φ q H J 1 + α f ( t , v ( t ) ) | H J 1 + β | Φ q H J 1 + α f ( t , u ( t ) ) Φ q H J 1 + α f ( t , v ( t ) ) | ( q 1 ) M Γ ( α + 1 ) q 2 H J 1 + β | H J 1 + α f ( t , u ( t ) ) H J 1 + α f ( t , v ( t ) ) | ( q 1 ) M Γ ( α + 1 ) q 2 H J 1 + β H J 1 + α | f ( t , u ( t ) ) f ( t , v ( t ) ) | ( q 1 ) M Γ ( α + 1 ) q 2 H J 1 + α + β [ l ( t ) | u ( t ) v ( t ) | ] ( q 1 ) M Γ ( α + 1 ) q 2 1 Γ ( α + β + 1 ) l e · ρ ( u , v ) .
It follows from (9) and (11) that
| T ( u ( t ) ) T ( v ( t ) ) | = | ( A u ( e ) A v ( e ) ) ( log t ) β 1 + ( A u ( t ) A v ( t ) ) | | A u ( e ) A v ( e ) | ( log t ) β 1 + | ( A u ( t ) A v ( t ) ) | 2 ( q 1 ) Γ ( α + β + 1 ) M Γ ( α + 1 ) q 2 l e · ρ ( u , v ) = κ ρ ( u , v ) .
In light of (12), we have
ρ ( T ( u ) , T ( v ) ) κ ρ ( u , v ) , u , v X , t [ 1 , e ] .
According to ( H 4 ) and (13), we know that ( a 4 ) in Lemma 3 holds.
Similar to (11), noticing that log t s α + β 1 and 1 s are monotonically decreasing with respect to s in [ 1 , e ] , we have
| A u ( t ) A v ( t ) | ( q 1 ) M Γ ( α + 1 ) q 2 H J 1 + α + β [ l ( t ) | u ( t ) v ( t ) | ] = ( q 1 ) M Γ ( α + 1 ) q 2 1 Γ ( α + β ) l e 1 t log t s α + β 1 | u ( s ) v ( s ) | d s s ( q 1 ) M Γ ( α + 1 ) q 2 1 Γ ( α + β ) l e · log t α + β 1 1 t | u ( s ) v ( s ) | d s q 1 Γ ( α + β ) M Γ ( α + 1 ) q 2 l e · log e α + β 1 1 e | u ( s ) v ( s ) | d s = q 1 Γ ( α + β ) M Γ ( α + 1 ) q 2 l e · ϱ ( u , v ) .
Similar to (12), we derive from (14) that
| T ( u ( t ) ) T ( v ( t ) ) | 2 ( q 1 ) Γ ( α + β ) M Γ ( α + 1 ) q 2 l e · ϱ ( u , v ) .
From (15), we get
ϱ ( T ( u ) , T ( v ) ) = 1 e | T ( u ( t ) ) T ( v ( t ) ) | d t 2 ( q 1 ) e Γ ( α + β ) M Γ ( α + 1 ) q 2 l e · ϱ ( u , v ) .
Equation (16) indicates that ( a 1 ) in Lemma 3 also holds. Thus, it follows from Lemma 3 that T exists a unique fixed point u * ( t ) X , which is the unique solution of (1). The proof is completed. □
Next, we shall discuss the approximation of solution for system (1). In fact, from Lemma 3, we conclude that the unique solution u * ( t ) X of (1) satisfies u * ( t ) = lim n T n u 0 for any u 0 X . Denote u n ( t ) = T n u 0 , then { u n ( t ) } is a approximation sequence of solution to system (1). Based on (9), u n ( t ) can be represented as
u n ( t ) = A u n 1 ( e ) ( log t ) β 1 + A u n 1 ( t ) ,
where A u n 1 ( t ) = H J 1 + β Φ q H J 1 + α f ( t , u n 1 ( t ) ) .
Similar to (12), we derive from (17) that
| u n + 1 ( t ) u n ( t ) | κ | u n ( t ) u n 1 ( t ) | ,
which implies that
ρ ( u n + 1 , u n ) κ ρ ( u n , u n 1 ) .
By virtue of ( H 4 ) and (18), we know that { u n ( t ) } converges exponentially on ( X , ρ ) .

4. Generalized Ulam-Hyers Stability

This section centres on the GUH-stability of BVP (1). We first provide the concept of GUH-stability for BVP (1).
For all δ > 0 , consider the following fractional differential inequality
H D 1 + α Φ p ( H D 1 + β u ( t ) ) f ( t , u ( t ) ) δ , t ( 1 , e ] , u ( 1 ) = u ( e ) , H D 1 + β u ( 1 ) = H D 1 + β u ( e ) .
Definition 3. 
BVP (1) is said to be generalized Ulam-Hyers (GUH) stable on the metric space ( X , ρ ) , provided that, for all δ > 0 and any solution u X of (19), there have an ω C ( R , R + ) with ω ( 0 ) = 0 and a unique solution u * X of (1) such that
ρ ( u , u * ) ω ( δ ) .
Remark 1. 
u ( t ) X solves the inequality 19) iff there has a continuous function  φ ( t ) such that
| φ ( t ) | δ , t ( 1 , e ] , H D 1 + α Φ p ( H D 1 + β u ( t ) ) = f ( t , u ( t ) ) + φ ( t ) , t ( 1 , e ] , u ( 1 ) = u ( e ) , H D 1 + β u ( 1 ) = H D 1 + β u ( e ) .
Theorem 2. 
If ( H 1 ) ( H 4 ) are satisfied, then BVP (1) is GUH-stable.
Proof. 
On the basis of Lemma 4 and Remark 1, the solution u ( t ) of inequality (19) is written by
u ( t ) = A u φ ( e ) ( log t ) β 1 + A u φ ( t ) ,
where A u φ ( t ) = H J 1 + β Φ q H J 1 + α [ f ( t , u ( t ) ) + φ ( t ) ] . In the light of Theorem 1 and Lemma 4, the unique solution u * ( t ) of BVP (1) is read as
u * ( t ) = A u * ( e ) ( log t ) β 1 + A u * ( t ) ,
where A u * ( t ) = H J 1 + β Φ q H J 1 + α f ( t , u * ( t ) ) .
Similar to (11), we obtain
| A u φ ( t ) A u * ( t ) | = | H J 1 + β Φ q H J 1 + α [ f ( t , u ( t ) ) + φ ( t ) ] Φ q H J 1 + α f ( t , u * ( t ) ) | H J 1 + β | Φ q H J 1 + α [ f ( t , u ( t ) ) + φ ( t ) ] Φ q H J 1 + α f ( t , u * ( t ) ) | ( q 1 ) M + δ Γ ( α + 1 ) q 2 H J 1 + β | H J 1 + α [ f ( t , u ( t ) ) + φ ( t ) ] H J 1 + α f ( t , u * ( t ) ) | ( q 1 ) M + δ Γ ( α + 1 ) q 2 H J 1 + β H J 1 + α [ | f ( t , u ( t ) ) f ( t , u * ( t ) ) | + | φ ( t ) | ] ( q 1 ) M + δ Γ ( α + 1 ) q 2 H J 1 + α + β [ l ( t ) | u ( t ) u * ( t ) | + | φ ( t ) | ] ( q 1 ) M + δ Γ ( α + 1 ) q 2 1 Γ ( α + β + 1 ) [ l e · ρ ( u , u * ) + δ ] .
By the same manner of (12), we derive from (21)–(23) that
| u ( t ) u * ( t ) | = | ( A u φ ( e ) A u * ( e ) ) ( log t ) β 1 + ( A u φ ( t ) A u * ( t ) ) | | A u φ ( e ) A u * ( e ) | ( log t ) β 1 + | ( A u φ ( t ) A u * ( t ) ) | 2 ( q 1 ) Γ ( α + β + 1 ) M + δ Γ ( α + 1 ) q 2 [ l e · ρ ( u , u * ) + δ ] = κ ( δ ) ρ ( u , u * ) + λ ( δ ) ,
where κ ( δ ) = 2 ( q 1 ) Γ ( α + β + 1 ) M + δ Γ ( α + 1 ) q 2 l e , λ ( δ ) = 2 ( q 1 ) δ Γ ( α + β + 1 ) M + δ Γ ( α + 1 ) q 2 .
For any sufficiently small δ > 0 , the condition ( H 4 ) ensures that 0 < κ ( δ ) < 1 . Thus, we know from (24) that
ρ ( u , u * ) λ ( δ ) 1 κ ( δ ) = ω ( δ ) .
Obviously, κ ( 0 ) = κ < 1 , λ ( 0 ) = 0 and ω ( δ ) = λ ( δ ) 1 κ ( δ ) > 0 with ω ( 0 ) = 0 . By virtue of Definition 3, (25) shows that BVP (1) is GUH-stable. The proof is completed. □

5. Two Examples and Simulations

This section provides two examples and simulations to inspect the correctness and validity of our main results. Consider the following nonlinear Hadamard fractional differential equation with p -Laplacian operator
H D 1 + α Φ p ( H D 1 + β u ( t ) ) = f ( t , u ( t ) ) , t ( 1 , e ] , u ( 1 ) = u ( e ) , H D 1 + β u ( 1 ) = H D 1 + β u ( e ) ,
Example 1. 
In (26), we take p = 5 4 , α = 0.2 , β = 1.4 , f ( t , u ) = 2 + sin ( 3 t ) 20 [ 3 π 4 + arctan ( u ) ] , then a simple computation gives that q = 5 > 2 , and
π 80 f ( t , u ) 3 π 16 , | f ( t , u ) f ( t , v ) | 2 + sin ( 3 t ) 20 | u v | .
In consequent, the conditions ( H 1 ) ( H 3 ) are fulfilled. In addition, M = 3 π 16 , l ( t ) = 2 + sin ( 3 t ) 20 , l e = 3 20 , and
κ = 2 ( q 1 ) Γ ( α + β + 1 ) M Γ ( α + 1 ) q 2 l e 0.1716 < 1 .
Thus, ( H 4 ) holds. From Theorem 1 and Theorem 2, we claim that Example 1 has a unique solution, which is GUH-stable.
Remark 2. 
In Example 1, p , α , β are all rational number. α = 0.2 is close to 0, and β = 1.4 is close to 1.5. To further verify the correctness of our results and the sensitivity of numerical simulation to parameters, we choose p , α , β as irrational number satisfying α close to 1 and β close to 2 in the following example.
Example 2. 
In (26), Choose p = 15 2 , α = 0.9 , β = 3.9 , and f ( t , u ) be same as Example 1. Then q = 2.0678 > 2 , and the conditions ( H 1 ) ( H 3 ) also hold. M, l ( t ) and l e are same as Example 1. In addition,
κ = 2 ( q 1 ) Γ ( α + β + 1 ) M Γ ( α + 1 ) q 2 l e 0.0567 < 1 .
Thus, ( H 4 ) is also true. From Theorems 1 and 2, we claim that Example 2 also has a unique GUH-stable solution.
To perform the numerical simulation on Examples 1 and 2, we need to give a concise algorithm below. Let v ( t ) = H D 1 + β u ( t ) , then the Equation (2) can be rewritten as
u ( t ) = 1 Γ ( β ) 1 e log e s β 1 v ( s ) d s s · ( log t ) β 1 + 1 Γ ( β ) 1 t log t s β 1 v ( s ) d s s , v ( t ) = 1 Γ ( α ) 1 t log t s α 1 f ( s , u ( s ) ) d s s q 1 .
Taking the derivative at both sides of (27), we get
d u ( t ) d t = β 1 Γ ( β ) 1 e log e s β 1 v ( s ) d s s · ( log t ) β 2 t + β 1 t Γ ( β ) 1 t log t s β 2 v ( s ) d s s , d v ( t ) d t = ( q 1 ) 1 Γ ( α ) 1 t log t s α 1 f ( s , u ( s ) ) d s s q 2 × α 1 t Γ ( α ) 1 t log t s α 2 f ( s , u ( s ) ) d s s .
For (28), we can apply the appropriate ODE toolbox in MATLAB to perform numerical solutions and simulations.
Based on the above algorithm, we employ the ODE113 toolbox in MATLAB R2019b on two examples to give their numerical solutions and simulations. Example 1 is shown as Table 1 and Table 2, Figure 1 and Figure 2. Example 2 is shown as Table 3 and Table 4, Figure 3 and Figure 4. Figure 2 and Table 1 and Table 2 show that Example 1 is GUH-stable. Figure 4 and Table 3 and Table 4 show that Example 2 is GUH-stable.

6. Summaries

It is well known that the p -Laplacian differential equation arises from the turbulence problem in porous medium. In viscoelastic mechanics, some studies have shown that fractional order differential equation models are more accurate than integer order differential equation models. Therefore, the fractional p -Laplacian differential model has greater advantages in the study of viscoelastic porous medium turbulence. In this article, we study BVP (1) of a nonlinear p -Laplacian Hadamard fractional differential equation. Unlike many published papers, we have established the existence, uniqueness, stability, and sequence approximation of solutions for fractional order differential equations on a wide range of complete metric spaces rather than Banach spaces. We have obtained some concise and easily verifiable sufficient criteria. Examples 1 and 2 and simulations demonstrate that our main results are correct and available. Meanwhile, Figure 1 and Figure 2 also indicate that the solution of BVP (1) is sensitive and dependent on parameters p , α and β . Our results provide theoretical support for revealing the mechanical problems of viscoelastic porous medium turbulence. The mathematical theories and methods used in the article have certain generality in solving similar problems. In addition, based on our recent research findings [50,51,52,53], we plan to study some ecosystems involving fractional derivatives or reaction diffusion effects in the future.

Funding

The APC was funded by research start-up funds for high-level talents of Taizhou University.

Data Availability Statement

Not applicable.

Acknowledgments

The author would like to express his heartfelt gratitude to the reviewers for their constructive comments.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Leibenson, L. General problem of the movement of a compressible uid in a porous medium. Izv. Akad. Nauk Kirg. SSR Ser. Biol. Nauk 1983, 9, 7–10. [Google Scholar]
  2. Alsaedi, A.; Luca, R.; Ahmad, B. Existence of positive solutions for a system of singular fractional boundary value problems with p -Laplacian operators. Mathematics 2020, 8, 1890. [Google Scholar] [CrossRef]
  3. Zhao, K.H. Solvability and GUH-stability of a nonlinear CF-fractional coupled Laplacian equations. AIMS Math. 2023, 8, 13351–13367. [Google Scholar]
  4. Rao, S.; Ahmadini, A. Multiple positive solutions for system of mixed Hadamard fractional boundary value problems with (p1,p2)-Laplacian operator. AIMS Math. 2023, 8, 14767–14791. [Google Scholar] [CrossRef]
  5. Alsaedi, A.; Alghanmi, M.; Ahmad, B.; Alharbi, B. Uniqueness of solutions for a ψ-Hilfer fractional integral boundary value problem with the p -Laplacian operator. Demonstr. Math. 2023, 56, 20220195. [Google Scholar]
  6. Sun, B.; Zhang, S.; Jiang, W. Solvability of fractional functional boundary-value problems with p -Laplacian operator on a half-line at resonance. J. Appl. Anal. Comput. 2023, 13, 11–33. [Google Scholar]
  7. Ahmadkhanlu, A. On the existence and multiplicity of positive solutions for a p -Laplacian fractional boundary value problem with an integral boundary condition. Filomat 2023, 37, 235–250. [Google Scholar]
  8. Chabane, F.; Benbachir, M.; Hachama, M.; Samei, M.E. Existence of positive solutions for p -Laplacian boundary value problems of fractional differential equations. Bound. Value Probl. 2022, 2022, 65. [Google Scholar] [CrossRef]
  9. Rezapour, S.; Abbas, M.; Etemad, S.; Minh Dien, N. On a multi-point p -Laplacian fractional differential equation with generalized fractional derivatives. Math. Method. Appl. Sci. 2023, 46, 8390–8407. [Google Scholar] [CrossRef]
  10. Boutiara, A.; Abdo, M.; Almalahi, M.; Shah, K.; Abdalla, B.; Abdeljawad, T. Study of Sturm-Liouville boundary value problems with p -Laplacian by using generalized form of fractional order derivative. AIMS Math. 2022, 7, 18360–18376. [Google Scholar] [CrossRef]
  11. Salem, A.; Almaghamsi, L.; Alzahrani, F. An infinite system of fractional order with p -Laplacian operator in a tempered sequence space via measure of noncompactness technique. Fractal Fract. 2021, 5, 182. [Google Scholar] [CrossRef]
  12. Jong, K.; Choi, H.; Kim, M.; Kim, K.; Jo, S.; Ri, O. On the solvability and approximate solution of a one-dimensional singular problem for a p -Laplacian fractional differential equation. Chaos Soliton Fract. 2021, 147, 110948. [Google Scholar]
  13. Li, S.; Zhang, Z.; Jiang, W. Multiple positive solutions for four-point boundary value problem of fractional delay differential equations with p -Laplacian operator. Appl. Numer. Math. 2021, 165, 348–356. [Google Scholar] [CrossRef]
  14. Hadamard, J. Essai sur l’étude des fonctions données par leur développment de Taylor. J. Math. Pures Appl. 1892, 8, 101–186. [Google Scholar]
  15. Kilbas, A.A. Hadamard-type fractional calculus. J. Korean Math. Soc. 2001, 38, 1191–1204. [Google Scholar]
  16. Butzer, P.L.; Kilbas, A.A.; Trujillo, J.J. Compositions of Hadamard-type fractional integration operators and the semigroup property. J. Math. Anal. Appl. 2002, 269, 387–400. [Google Scholar] [CrossRef] [Green Version]
  17. Ahmad, B.; Ntouyas, S.K. On Hadamard fractional integro-differential boundary value problems. J. Appl. Math. Comput. 2015, 47, 119–131. [Google Scholar] [CrossRef]
  18. Aljoudi, S.; Ahmad, B.; Nieto, J.J.; Alsaedi, A. A coupled system of Hadamard type sequential fractional differential equations with coupled strip conditions. Chaos Soliton Fract. 2016, 91, 39–46. [Google Scholar]
  19. Benchohra, M.; Bouriah, S.; Graef, J.R. Boundary value problems for nonlinear implicit Caputo-Hadamard-type fractional differential equations with impulses. Mediterr. J. Math. 2017, 14, 206. [Google Scholar] [CrossRef]
  20. Huang, H.; Zhao, K.H.; Liu, X.D. On solvability of BVP for a coupled Hadamard fractional systems involving fractional derivative impulses. AIMS Math. 2022, 7, 19221–19236. [Google Scholar] [CrossRef]
  21. Zhao, K.H. Existence and UH-stability of integral boundary problem for a class of nonlinear higher-order Hadamard fractional Langevin equation via Mittag-Leffler functions. Filomat 2023, 37, 1053–1063. [Google Scholar] [CrossRef]
  22. Ulam, S. A Collection of Mathematical Problems-Interscience Tracts in Pure and Applied Mathmatics; Interscience: New York, NY, USA, 1906. [Google Scholar]
  23. Hyers, D. On the stability of the linear functional equation. Proc. Natl. Acad. Sci. USA 1941, 27, 2222–2240. [Google Scholar] [CrossRef] [Green Version]
  24. Zada, A.; Waheed, H.; Alzabut, J.; Wang, X. Existence and stability of impulsive coupled system of fractional integrodifferential equations. Demonstr. Math. 2019, 52, 296–335. [Google Scholar] [CrossRef]
  25. Yu, X. Existence and β-Ulam-Hyers stability for a class of fractional differential equations with non-instantaneous impulses. Adv. Differ. Equ. 2015, 2015, 104. [Google Scholar] [CrossRef] [Green Version]
  26. Zhao, K.H. Stability of a nonlinear fractional Langevin system with nonsingular exponential kernel and delay control. Discret. Dyn. Nat. Soc. 2022, 2022, 9169185. [Google Scholar] [CrossRef]
  27. Chen, C.; Li, M. Existence and Ulam type stability for impulsive fractional differential systems with pure delay. Fractal Fract. 2022, 6, 742. [Google Scholar] [CrossRef]
  28. Zhao, K.H. Existence, stability and simulation of a class of nonlinear fractional Langevin equations involving nonsingular Mittag-Leffler kernel. Fractal Fract. 2022, 6, 469. [Google Scholar] [CrossRef]
  29. Zhao, K.H. Stability of a nonlinear Langevin system of ML-type fractional derivative affected by time-varying delays and differential feedback control. Fractal Fract. 2022, 6, 725. [Google Scholar] [CrossRef]
  30. Mehmood, M.; Abbas, A.; Akgul, A.; Asjad, M.I.; Eldin, S.M.; Abd El-Rahman, M.; Baleanu, D. Existence and stability results for coupled system of fractional differential equations involving AB-caputo derivative. Fractals 2023, 31, 2340023. [Google Scholar] [CrossRef]
  31. Yaghoubi, H.; Zare, A.; Rasouli, M.; Alizadehsani, R. Novel frequency-based approach to analyze the stability of polynomial fractional differential equations. Axioms 2023, 12, 147. [Google Scholar] [CrossRef]
  32. Alqahtani, R.; Ntonga, J.; Ngondiep, E. Stability analysis and convergence rate of a two-step predictor-corrector approach for shallow water equations with source terms. AIMS Math. 2023, 8, 9265–9289. [Google Scholar]
  33. Salama, F.; Ali, N.; Hamid, N. Fast O(N) hybrid Laplace transform-finite difference method in solving 2D time fractional diffusion equation. J. Math. Comput. Sci. 2021, 23, 110–123. [Google Scholar] [CrossRef]
  34. Jassim, H.; Hussain, M. On approximate solutions for fractional system of differential equations with Caputo-Fabrizio fractional operator. J. Math. Comput. Sci. 2021, 23, 58–66. [Google Scholar] [CrossRef]
  35. Can, N.; Nikan, O.; Rasoulizadeh, M.; Jafari, H.; Gasimov, Y.S. Numerical computation of the time non-linear fractional generalized equal width model arising in shallow water channel. Therm. Sci. 2020, 24, 49–58. [Google Scholar] [CrossRef]
  36. Akram, T.; Abbas, M.; Ali, A. A numerical study on time fractional Fisher equation using an extended cubic B-spline approximation. J. Math. Comput. Sci. 2021, 22, 85–96. [Google Scholar] [CrossRef]
  37. Ahmad, I.; Ahmad, H.; Thounthong, P.; Chu, Y.M.; Cesarano, C. Solution of multi-term time-fractional PDE models arising in mathematical biology and physics by local meshless method. Symmetry 2020, 12, 1195. [Google Scholar] [CrossRef]
  38. Murtaza, S.; Ahmad, Z.; Ali, I.; Chu, Y.M.; Cesarano, C. Analysis and numerical simulation of fractal-fractional order non-linear couple stress nanofluid with cadmium telluride nanoparticles. J. King Saud Univ. Sci. 2023, 35, 102618. [Google Scholar] [CrossRef]
  39. Zhao, K.H. Generalized UH-stability of a nonlinear fractional coupling (p1,p2)-Laplacian system concerned with nonsingular Atangana-Baleanu fractional calculus. J. Inequal. Appl. 2023; accepted. [Google Scholar]
  40. Ahmad, Z.; Bonanomi, G.; di Serafino, D.; Giannino, F. Transmission dynamics and sensitivity analysis of pine wilt disease with asymptomatic carriers via fractal-fractional differential operator of Mittag-Leffler kernel. Appl. Numer. Math. 2023, 185, 446–465. [Google Scholar] [CrossRef]
  41. Ahmad, Z.; El-Kafrawy, S.; Alandijany, T.; Mirza, A.A.; El-Daly, M.M.; Faizo, A.A.; Bajrai, L.H.; Kamal, M.A.; Azhar, E.I. A global report on the dynamics of COVID-19 with quarantine and hospitalization: A fractional order model with non-local kernel. Comput. Biol. Chem. 2022, 98, 107645. [Google Scholar] [CrossRef]
  42. Ahmad, Z.; Ali, F.; Khan, N.; Khan, I. Dynamics of fractal-fractional model of a new chaotic system of integrated circuit with Mittag-Leffler kernel. Chaos Soliton. Fract. 2021, 153, 111602. [Google Scholar] [CrossRef]
  43. Wang, F.; Asjad, M.; Zahid, M.; Iqbal, A.; Ahmad, H.; Alsulami, M.D. Unsteady thermal transport flow of Casson nanofluids with generalized Mittag-Leffler kernel of Prabhakar’s type. J. Mater. Res. Technol. 2021, 14, 1292–1300. [Google Scholar] [CrossRef]
  44. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; Elsevier: Amsterdam, The Netherlands, 2006. [Google Scholar]
  45. Rus, I.A. On a fixed point theorem of Maia. Stud. Univ. Babes-Bolyai Math. 1977, 22, 40–42. [Google Scholar]
  46. Wardowski, D. Fixed points of a new type of contractive mappings in complete metric spaces. Fixed Point Theory Appl. 2012, 2012, 94. [Google Scholar] [CrossRef] [Green Version]
  47. Agarwal, R.; O’Regan, D. Fixed point theory for generalized contractions on spaces with two metrics. J. Math. Anal. Appl. 2000, 248, 402–414. [Google Scholar] [CrossRef] [Green Version]
  48. Karapinar, E.; Fulga, A.; Agarwal, R. A survey: F-contractions with related fixed point results. J. Fixed Point Theory Appl. 2020, 22, 69. [Google Scholar] [CrossRef]
  49. Mehmood, N.; Khan, I.; Nawaz, M.; Ahmad, N. Existence results for ABC-fractional BVP via new fixed point results of F-Lipschitzian mappings. Demonstr. Math. 2022, 55, 452–469. [Google Scholar] [CrossRef]
  50. Zhao, K.H. Local exponential stability of several almost periodic positive solutions for a classical controlled GA-predation ecosystem possessed distributed delays. Appl. Math. Comput. 2023, 437, 127540. [Google Scholar] [CrossRef]
  51. Zhao, K.H. Existence and stability of a nonlinear distributed delayed periodic AG-ecosystem with competition on time scales. Axioms 2023, 12, 315. [Google Scholar] [CrossRef]
  52. Zhao, K.H. Global stability of a novel nonlinear diffusion online game addiction model with unsustainable control. AIMS Math. 2022, 7, 20752–20766. [Google Scholar] [CrossRef]
  53. Zhao, K.H. Attractor of a nonlinear hybrid reaction-diffusion model of neuroendocrine transdifferentiation of human prostate cancer cells with time-lags. AIMS Math. 2023, 8, 14426–14448. [Google Scholar] [CrossRef]
Figure 1. Simulation of solutions to Example 1.
Figure 1. Simulation of solutions to Example 1.
Axioms 12 00733 g001
Figure 2. Evolution of the GUH-stability of Example 1.
Figure 2. Evolution of the GUH-stability of Example 1.
Axioms 12 00733 g002
Figure 3. Simulation of solutions to Example 2.
Figure 3. Simulation of solutions to Example 2.
Axioms 12 00733 g003
Figure 4. Evolution of the GUH-stability of Example 2.
Figure 4. Evolution of the GUH-stability of Example 2.
Axioms 12 00733 g004
Table 1. The numerical solution u ( t ) to Example 1 which needs to multiply by 10 10 .
Table 1. The numerical solution u ( t ) to Example 1 which needs to multiply by 10 10 .
t1.21.41.61.82.02.22.42.62.7183
u
δ = 0 0.00790.03140.04970.06450.07860.09630.12100.14480.1501
δ = 0.001 0.00800.03180.05010.06500.07900.09710.12180.14470.1498
Table 2. The numerical solution v ( t ) to Example 1 which needs to multiply by 10 10 .
Table 2. The numerical solution v ( t ) to Example 1 which needs to multiply by 10 10 .
t1.21.41.61.82.02.22.42.62.7183
u
δ = 0 0.16530.26900.33440.43020.73191.73074.03107.58439.7746
δ = 0.001 0.16710.27200.33550.43270.73191.77493.90197.26219.3765
Table 3. The numerical solution u ( t ) to Example 2.
Table 3. The numerical solution u ( t ) to Example 2.
t1.21.41.61.82.02.22.42.62.7183
u
δ = 0 0.00180.01210.02520.03840.05070.06140.06990.07490.0757
δ = 0.001 0.00190.01220.02540.03870.05110.06190.07040.07540.0763
Table 4. The numerical solution v ( t ) to Example 2 which needs to multiply by 1 .
Table 4. The numerical solution v ( t ) to Example 2 which needs to multiply by 1 .
t1.21.41.61.82.02.22.42.62.7183
u
δ = 0 0.08180.19470.27570.35300.44810.57210.71890.87410.9639
δ = 0.001 0.08240.19600.27790.35560.45160.57670.72440.88280.9723
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

Zhao, K. Solvability, Approximation and Stability of Periodic Boundary Value Problem for a Nonlinear Hadamard Fractional Differential Equation with p-Laplacian. Axioms 2023, 12, 733. https://0-doi-org.brum.beds.ac.uk/10.3390/axioms12080733

AMA Style

Zhao K. Solvability, Approximation and Stability of Periodic Boundary Value Problem for a Nonlinear Hadamard Fractional Differential Equation with p-Laplacian. Axioms. 2023; 12(8):733. https://0-doi-org.brum.beds.ac.uk/10.3390/axioms12080733

Chicago/Turabian Style

Zhao, Kaihong. 2023. "Solvability, Approximation and Stability of Periodic Boundary Value Problem for a Nonlinear Hadamard Fractional Differential Equation with p-Laplacian" Axioms 12, no. 8: 733. https://0-doi-org.brum.beds.ac.uk/10.3390/axioms12080733

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