Next Article in Journal
Existence and Uniqueness Results of Coupled Fractional-Order Differential Systems Involving Riemann–Liouville Derivative in the Space Wa+γ1,1(a,b)×Wa+γ2,1(a,b) with Perov’s Fixed Point Theorem
Next Article in Special Issue
Nonexistence of Global Solutions to Time-Fractional Damped Wave Inequalities in Bounded Domains with a Singular Potential on the Boundary
Previous Article in Journal
Certain Recurrence Relations of Two Parametric Mittag-Leffler Function and Their Application in Fractional Calculus
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mawhin’s Continuation Technique for a Nonlinear BVP of Variable Order at Resonance via Piecewise Constant Functions

by
Shahram Rezapour
1,2,†,
Mohammed Said Souid
3,†,
Sina Etemad
1,†,
Zoubida Bouazza
4,†,
Sotiris K. Ntouyas
5,6,†,
Suphawat Asawasamrit
7,† and
Jessada Tariboon
7,*,†
1
Department of Mathematics, Azarbaijan Shahid Madani University, Tabriz 53751-71379, Iran
2
Department of Medical Research, China Medical University Hospital, China Medical University, Taichung 40402, Taiwan
3
Department of Economic Sciences, University of Tiaret, Tiaret 14035, Algeria
4
Laboratory of Mathematics, Djillali Liabes University, Sidi Bel-Abbès 22000, Algeria
5
Department of Mathematics, University of Ioannina, 451 10 Ioannina, Greece
6
Nonlinear Analysis and Applied Mathematics (NAAM)-Research Group, Department of Mathematics, Faculty of Science, King Abdulaziz University, P.O. Box 80203, Jeddah 21589, Saudi Arabia
7
Intelligent and Nonlinear Dynamic Innovations Research Center, Department of Mathematics, Faculty of Applied Science, King Mongkut’s University of Technology North Bangkok, Bangkok 10800, Thailand
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Fractal Fract. 2021, 5(4), 216; https://doi.org/10.3390/fractalfract5040216
Submission received: 17 October 2021 / Revised: 9 November 2021 / Accepted: 10 November 2021 / Published: 12 November 2021

Abstract

:
In this paper, we establish the existence of solutions to a nonlinear boundary value problem (BVP) of variable order at resonance. The main theorem in this study is proved with the help of generalized intervals and piecewise constant functions, in which we convert the mentioned Caputo BVP of fractional variable order to an equivalent standard Caputo BVP at resonance of constant order. In fact, to use the Mawhin’s continuation technique, we have to transform the variable order BVP into a constant order BVP. We prove the existence of solutions based on the existing notions in the coincidence degree theory and Mawhin’s continuation theorem (MCTH). Finally, an example is provided according to the given variable order BVP to show the correctness of results.

1. Introduction

The initial idea of fractional calculus is taken from the powers of real or complex numbers in the order of differentiation and integration operators. In recent decades, fractional operators of variable order are appeared extensively in a vast domain of sciences including chaotic dynamical systems, fractal theory, rheology, signal processing, mathematical modeling, control theory, and biomedical applications. This range of applications is due to the fact that fractional derivatives provide a strong tool in the mathematics to describe the memory and hereditary properties of processes and various materials; see, for example [1,2,3].
Before the variable order systems, discussion of boundary value problems with fractional constant orders has attracted the attention of most researchers, and valuable findings have been established. Various researches have been conducted to study the behaviors of different fractional BVPs by means of some known methods such as fixed point theorems, numerical methods, monotone iterative methods, variational methods, and etc. [4,5,6,7,8,9,10,11,12].
Nevertheless, in addition to numerous published papers on fractional constant order problems, few studies on the existence theory have been done in relation to variable order problems [13,14,15,16,17,18,19]. Hence, investigation of this interesting and general topic makes all our findings worthy.
In 1970, Gaines and Mawhin [20] introduced the theory of coincidence degree for analysis of differential and functional equations. Mawhin has made important contributions since then, and the mentioned theory is also famous to the Mawhin’s coincidence theory. Coincidence theory is considered as a powerful technique, especially with regard to questions about the existence of solutions for nonlinear differential equations. Mawhin’s theory permits the use of a method based on the topological degree notion for some problems which can be written as an abstract operator equation of the form Θ x = W x , where Θ is a linear non-invertible operator and W is a nonlinear operator acting on a Banach space.
In 1972, Mawhin extended a technique to solve this operator equation in his famous paper [21]. He assumed that Θ is a Fredholm operator of index zero. Then, he developed a new theory of topological degree known as the degree of coincidence for ( Θ , W ) , that is also known as Mawhin’s coincidence degree theory in honor of him.
A given boundary value problem is said to be at resonance if the corresponding linear homogeneous BVP has a non-trivial solution. Many authors studied ordinary BVPs at resonance using Mawhin’s coincidence degree theory; we can cite some works done by Feng and Webb [22], Guezane-Lakoud and Frioui [23], Mawhin and Ward [24], Infante [25], and references therein.
Based on the aforementioned technique in relation to Mawhin’s method, in this paper, we shall investigate a nonlinear boundary value problem of variable order at resonance which takes the form as follows
c D 0 + u ( t ) ϕ ( t ) = g ( t , ϕ ( t ) ) , t A , ϕ ( 0 ) = ϕ ( T ) ,
where A = [ 0 , T ] , T ( 0 , ) , the function u ( t ) : A ( 0 , 1 ] is the order of the existing derivative in the above boundary problem, c D 0 + u ( t ) is the variable order Caputo derivative, and also g C ( A × R , R ) .
The important aim of this research is to investigate some qualitative properties of solutions of the given Caputo boundary value problem of variable order (1). The main novelty of this paper is that we use the Mawhin’s continuation technique for the first time for proving the existence of solutions of a Caputo boundary value problem at resonance equipped with variable order. Most papers apply this technique on the constant order systems, while we here try to derive the necessary conditions on a variable order system. In comparison to variable order partial systems, a linear analogue of this problem can be observed in the framework of partial differential equation [26] and this shows another version of such problems and specify our main contribution in this work. It is notable for young researchers that they can implement and investigate this methods and techniques on hidden-memory variable-order fractional problems introduced in [27,28] in the future.
The structure of the paper is organized as follows: Initially, some auxiliary definitions and remarks are collected for recalling the required notions in Section 2. Further, in Section 3, based on coincidence degree theory, a partition of the given interval A is applied, and by defining the relevant piecewise constant functions, the existence results are derived for an equivalent constant-order BVP at resonance and accordingly, for the given Caputo BVP of variable order (1). This proof is completed in some steps. In Section 4, we give an example to illustrate the theoretical existence theorems. The paper is completed with conclusions in Section 5.

2. Auxiliary Concepts

At first, some needed concepts about our study are collected from different sources. Here, the Banach space C ( A , R ) consisting of continuous functions like ϕ : A R is equipped with the sup-norm ϕ = sup { | ϕ ( t ) | : t A } .
Definition 1
([29,30]). The Riemann-Liouville fractional integral (RLFI) of variable order u ( t ) for the function ϕ is defined by
I 0 + u ( t ) ϕ ( t ) = 1 Γ ( u ( t ) ) 0 t ( t s ) u ( t ) 1 ϕ ( s ) d s , t A ,
where Γ ( z ) = 0 x z 1 e x d x , and the left Caputo fractional derivative (CFD) of variable order u ( t ) for ϕ ( t ) is defined by
c D 0 + u ( t ) ϕ ( t ) = 1 1 Γ ( u ( t ) ) 0 t ( t s ) u ( t ) ϕ ( s ) d s , t A .
Remark 1.
Notice that in (2), we have specified the variable order as the function u : A ( 0 , 1 ] , while for defining RLFI, we can consider it as a function with extended values like u : A ( 0 , ) .
Remark 2
([30]). When we define the variable order u as a constant-valued function in both (2) and (3), then the variable order RLFI and CFD operators are the same as the usual RLFI and CFD operators, respectively.
Remark 3
([29]). As we know, the semigroup property is satisfied for the standard RLFI operators equipped with constant order, while it is not valid for extended case of variable orders β 1 ( t ) and β 2 ( t ) . In other words, I 0 + β 1 ( t ) ( I 0 + β 2 ( t ) ) ϕ ( t ) I 0 + β 1 ( t ) + β 2 ( t ) ϕ ( t ) .
To see this problem, we give the following example.
Example 1.
Let A = [ 0 , 3 ] and ϕ ( t ) 1 , t A . The variable orders of RLFI operator can be taken as: β 1 ( t ) = t 2 and β 2 ( t ) = 1 , t [ 0 , 1 ] 2 , t [ 1 , 3 ] .
Then for all t A , and according to Definition (2), we compute
I 0 + β 1 ( t ) I 0 + β 2 ( t ) ϕ ( t ) = 0 t ( t s ) β 1 ( t ) 1 Γ ( β 1 ( t ) ) 0 s ( s τ ) β 2 ( s ) 1 Γ ( β 2 ( s ) ) ϕ ( τ ) d τ d s = 0 t ( t s ) β 1 ( t ) 1 Γ ( u ( t ) ) [ 0 1 ( s τ ) 0 Γ ( 1 ) d τ + 1 s ( s τ ) Γ ( 2 ) d τ ] d s = 0 t ( t s ) β 1 ( t ) 1 Γ ( β 1 ( t ) ) [ s 2 2 s + 3 2 ] d s ,
and
I 0 + β 1 ( t ) + β 2 ( t ) ϕ ( t ) = 0 t ( t s ) β 1 ( t ) + β 2 ( t ) 1 Γ ( b e t a 1 ( t ) + β 2 ( t ) ) ϕ ( s ) d s .
For t = 2 , it becomes
I 0 + β 1 ( t ) I 0 + β 2 ( t ) ϕ ( t ) | t = 2 = 0 2 ( 2 s ) 0 Γ ( 1 ) [ s 2 2 s + 3 2 ] d s = 0 2 ( s 2 2 s + 3 2 ) d s = 7 3 ,
and
I 0 + β 1 ( t ) + β 2 ( t ) ϕ ( t ) | t = 2 = 0 2 ( 2 s ) β 1 ( t ) + β 2 ( t ) 1 Γ ( β 1 ( t ) + β 2 ( t ) ) ϕ ( s ) d s = 0 1 ( 2 s ) 1 Γ ( 2 ) d s + 1 2 ( 2 s ) 2 Γ ( 3 ) d s = 3 2 + 1 6 = 5 3 .
Hence, it is simply seen that the mentioned property is not correct for the generalized RLFI operators with respect to variable orders.
The following expansion is key for our argument.
Lemma 1
([31]). Let a 1 , α 1 > 0 and n = 1 + [ α 1 ] . Then
I a 1 + α 1 ( c D a 1 + α 1 ϕ ( t ) ) = ϕ ( t ) k = 0 n 1 ϕ ( k ) ( a 1 ) k ! t k .
Lemma 2
([32]). Let α 1 , α 2 > 0 , ϕ , c D a 1 + α 1 ϕ L 1 ( a 1 , a 2 ) . Then, the differential equation
c D a 1 + α 1 ϕ ( t ) = 0 ,
has unique solution
ϕ ( t ) = r 0 + r 1 ( t a 1 ) + r 2 ( t a 1 ) 2 + + r n 1 ( t a 1 ) n 1 ,
and we have
I a 1 + α 1 ( c D a 1 + α 1 ) ϕ ( t ) = ϕ ( t ) + r 0 + r 1 ( t a 1 ) + r 2 ( t a 1 ) 2 + + r n 1 ( t a 1 ) n 1 ,
such that n 1 < α 1 n , r j R , j = 1 , 2 , , n .
Furthermore, we have
c D a 1 + α 1 ( I a 1 + α 1 ) ϕ ( t ) = ϕ ( t ) ,
and
I a 1 + α 1 ( I a 1 + α 2 ) ϕ ( t ) = I a 1 + α 2 ( I a 1 + α 1 ) ϕ ( t ) = I a 1 + α 1 + α 2 ϕ ( t ) .
We recall some properties of variable order RLFI operator formulated by (2) which will be used in the sequel.
Lemma 3
([33]). If u : A ( 0 , 1 ] has the continuity property, then for
h C δ ( A , R ) = { h ( t ) C ( A , R ) , t δ h ( t ) C ( A , R ) } , 0 δ 1 ,
the integral I 0 + u ( t ) h ( t ) admits a finite value for all t A .
Lemma 4
([33]). Assume that u : A ( 0 , 1 ] has the continuity property. Then
I 0 + u ( t ) h ( t ) C ( A , R )   for   h C ( A , R ) .
Definition 2
([34,35]). An interval J R is termed as a generalized interval if I is either an interval, or { a 1 } , or ∅. A finite set F is defined to be a partition of J if every x J belongs to exactly one and one generalized interval I in F . Finally, w : J R is piecewise constant w.r.t F as a partition of J, if for each I F , w is constant on I .
The next definitions and basic lemmas from coincidence degree theory are fundamental in the proof of theorems which we will establish them later.
Definition 3
([20,36]). Consider two normed spaces S 1 and S 2 . A Fredholm operator of index zero is a linear operator like Θ : D o m ( Θ ) S 1 S 2 satisfying:
(a) 
IMG ( Θ ) S 2 is closed;
(b) 
dim KER ( Θ ) = codim IMG ( Θ ) < + .
In view of Definition 3, it is followed the existence of continuous projections Ψ : S 1 S 1 and Φ : S 2 S 2 such that IMG ( Ψ ) = KER ( Θ ) , KER ( Φ ) = IMG ( Θ ) , S 1 = KER ( Θ ) KER ( Ψ ) , and S 2 = IMG ( Θ ) IMG ( Φ ) .
It is known that the restriction of Θ to D o m ( Θ ) KER ( Ψ ) , which we shall represent by Θ Ψ , will be an isomorphism onto its image [20,36].
Definition 4
([20,36]). Let Θ be a Fredholm operator of index zero and Ω S 1 be bounded with D o m ( Θ ) Ω . We say W : Ω ¯ S 2 has the Θ -compactness property in Ω ¯ whenever:
(H1) 
Φ W : Ω ¯ S 2 is continuous, and Φ W ( Ω ¯ ) S 2 is bounded,
(H2) 
( Θ Ψ ) 1 ( I Φ ) W : Ω ¯ S 1 is completely continuous.
The next theorem entitled Mawhin’s Continuation Theorem is our main criterion in the present study which proves the existence of solution.
Theorem 1
([37]). Assume that S 1 and S 2 are two Banach spaces and Ω S 1 is an open, bounded and symmetric set with 0 Ω . Also, assume that:
(A1) 
the Fredholm operator Θ : D o m ( Θ ) S 1 S 2 of index zero is such that
D o m ( Θ ) Ω ¯ ,
(A2) 
the operator W : S 1 S 2 is Θ -compact on Ω ¯ ,
(A3) 
x D o m ( Θ ) Ω and λ ( 0 , 1 ] ,
Θ x W x λ ( Θ x + W ( x ) ) ,
where Ω denotes the boundary of Ω w.r.t. S 1 .
Then the operator equation Θ x = W x has at least one solution on D o m ( Θ ) Ω ¯ .

3. Existence of Solutions

To begin the desired analysis, we consider the following assumptions:
(AS1) 
Consider a sequence of finite many points { T k } k = 0 n so that 1 = T 0 < T k < T n = T , k N 1 n 1 . For k N 1 n , denote the subintervals A k as A k : = ( T k 1 , T k ] . Then P = k = 1 n A k is a partition of A.
(AS2) 
Let g C ( A j × R , R ) and there exists δ ( 0 , 1 ) such that t δ g C ( A j × R , R ) and there exists K > 0 with K < min 1 , Γ ( u j + 1 ) ( T j T j 1 ) u j so that t δ | g ( t , ϕ 1 ) g ( t , ϕ 2 ) | K | ϕ 1 ϕ 2 | , for any ϕ 1 , ϕ 2 R and t A j .
For each j N 1 n , the notation E j = C ( A j , R ) denotes the Banach space of continuous functions ϕ : A j R with the sup-norm ϕ E j = sup t A j | ϕ ( t ) | .
On the other side, consider the piecewise constant mapping u ( t ) : A ( 0 , 1 ] w.r.t. P , i.e.,
u ( t ) = j = 1 n u j I j ( t ) ,
where 0 < u j 1 are real numbers, and I j denotes the indicator of A j , j N 1 n ; that is, I j ( t ) = 1 if t A j and I j ( t ) = 0 otherwise. In this case, the left CFD of variable order u ( t ) for ϕ ( t ) C ( A , R ) , defined as (3), can be formulated as a sum of the left CFD operators of constant orders u k R which takes the form
c D 0 + u ( t ) ϕ ( t ) = k = 1 j 1 T k 1 T k ( t s ) u k Γ ( 1 u k ) ϕ ( s ) d s + T j 1 t ( t s ) u j Γ ( 1 u j ) ϕ ( s ) d s .
Thus, the given Caputo BVP of variable order (1) can be reformulated for each t A j , j N 1 n in the following structure
k = 1 j 1 T k 1 T k ( t s ) u k Γ ( 1 u k ) ϕ ( s ) d s + T j 1 t ( t s ) u j Γ ( 1 u j ) ϕ ( s ) d s = g ( t , ϕ ( t ) ) .
Let the function ϕ ˜ E j be so that ϕ ˜ ( t ) 0 on t [ 0 , T j 1 ] and it satisfies the above integral Equation (5). In such a situation, (5) is converted to the standard constant order fractional differential equation (FDE) as
c D T j 1 + u j ϕ ˜ ( t ) = g ( t , ϕ ˜ ( t ) ) , t A j .
In accordance with above equation, for each j N 1 n , we have the auxiliary FBVP equipped with Caputo constant order CFD operator
c D T j 1 + u j ϕ ( t ) = g ( t , ϕ ( t ) ) , t A j , ϕ ( T j 1 ) = ϕ ( T j ) .
A resonance problem is a boundary problem in which the corresponding homogeneous BVP has a non–trivial solution. Hence, we consider the homogeneous version of the given equivalent constant order FBVP (6) by
c D T j 1 + u j ϕ ( t ) = 0 , t A j , ϕ ( T j 1 ) = ϕ ( T j ) .
By Lemma 2, the homogeneous constant order FBVP (7) has nontrivial solution ϕ ( t ) = c which converts the equivalent constant order FBVP (6) to a resonance FBVP.
As well as, on the given subintervals, let S 1 = { ϕ E j : ϕ ( t ) = I T j 1 + u j v ( t ) : v E j , t A j } with the norm
ϕ S 1 = ϕ E j .
The linear operator Θ : D o m ( Θ ) S 1 E j along with the operator W : S 1 E j are defined as
Θ [ ϕ ( t ) ] : = c D T j 1 + u j ϕ ( t ) ,
and
W [ ϕ ( t ) ] : = g ( t , ϕ ( t ) ) , t A j ,
where
D o m ( Θ ) = { ϕ S 1 : c D T j 1 + u j ϕ E j and ϕ ( T j 1 ) = ϕ ( T j ) } .
Then the equivalent constant order resonance FBVP (6) can be reformulated by the equation Θ ϕ = W ϕ .
The first theorem on the existence of solutions for the equivalent constant order resonance FBVP (6) is established in this position.
Theorem 2.
If the condition (AS2) holds, then the equivalent constant order resonance FBVP (6) has at least one solution.
Proof. 
The proof will be followed in a sequence of claims.
Claim 1. We show that
KER ( Θ ) = { c : c R } ,
and
IMG ( Θ ) = ϕ E j : T j 1 T j ( T j s ) u j 1 ϕ ( s ) d s = 0 .
Let Θ (defined by (8)) be such that for t A j and by Lemma 2, the equation Θ [ ϕ ( t ) ] = c D T j 1 + u j ϕ ( t ) = 0 has the solution ϕ ( t ) = c , c R . Then
KER ( Θ ) = { ϕ ( t ) = c : c R } .
On the other hand, for v IMG ( Θ ) , there exits ϕ D o m ( Θ ) such that v = Θ ϕ E j . By Lemma 1, for any t A j , we have
ϕ ( t ) = ϕ ( T j 1 ) + 1 Γ ( u j ) T j 1 t ( t s ) u j 1 v ( s ) d s .
Since ϕ D o m ( Θ ) , v satisfies
1 Γ ( u j ) T j 1 T j ( T j s ) u j 1 v ( s ) d s = 0 .
Also, assume that v E j satisfies
T j 1 T j ( T j s ) u j 1 v ( s ) d s = 0 .
Let ϕ ( t ) = I T j 1 + u j v ( t ) . Then v ( t ) = c D T j 1 + u j ϕ ( t ) and so ϕ D o m ( Θ ) . Hence, v IMG ( Θ ) , so
IMG ( Θ ) = ϕ S 2 : T j 1 T j ( T j s ) u j 1 ϕ ( s ) d s = 0 .
Claim 2.  Θ is a Fredholm operator of index zero.
The linear continuous projector operators Ψ : S 1 S 1 and Φ : E j E j can be considered by the following forms
Ψ ϕ = ϕ ( T j 1 ) , Φ v ( t ) = u j ( T j T j 1 ) u j T j 1 T j ( T j s ) u j 1 v ( s ) d s .
Clearly, IMG ( Ψ ) = KER ( Θ ) and Ψ 2 = Ψ . It follows that for any ϕ S 1 ,
ϕ = ( ϕ Ψ ϕ ) + Ψ ϕ ,
i.e., S 1 = KER ( Ψ ) + KER ( Θ ) . A simple computation shows that KER ( Ψ ) KER ( Θ ) = 0 . Therefore, S 1 = KER ( Ψ ) KER ( Θ ) . A similar argument shows that for every v E j , Φ 2 v = Φ v and v = ( v Φ ( v ) ) + Φ ( v ) , where ( v Φ ( v ) ) KER ( Φ ) = IMG ( Θ ) .
From IMG ( Θ ) = KER ( Φ ) and Φ 2 = Φ , we have
IMG ( Φ ) IMG ( Θ ) = 0 .
Then, E j = IMG ( Θ ) IMG ( Φ ) .
In this case,
dim ( KER ( Θ ) = dim IMG ( Φ ) = codim IMG ( Θ ) .
The obtained result shows that Θ is a Fredholm operator of index zero.
Claim 3.  Θ Ψ 1 = ( Θ | D o m ( Θ ) KER ( Ψ ) ) 1 (the inverse of Θ | D o m ( Θ ) KER ( Ψ ) ).
Clearly, Θ Ψ 1 : IMG ( Θ ) S 1 KER ( Ψ ) satisfies
Θ Ψ 1 ( v ) ( t ) = I T j 1 + u j v ( t ) .
Let v IMG ( Θ ) . Then
Θ Θ Ψ 1 ( v ) = c D T j 1 + u j ( I T j 1 + u j v ) = v .
Furthermore, for ϕ D o m ( Θ ) KER ( Ψ ) , we get
Θ Ψ 1 ( Θ ( ϕ ( t ) ) ) = I T j 1 + u j ( c D T j 1 + u j ϕ ( t ) ) = ϕ ( t ) ϕ ( T j 1 ) .
Since ϕ D o m ( Θ ) KER ( Ψ ) , we know that ϕ ( T j 1 ) = 0 . Thus
Θ Ψ 1 ( Θ ( ϕ ( t ) ) ) = ϕ ( t ) .
Combining (10) and (11) shows that Θ Ψ 1 = ( Θ | D o m ( Θ ) KER ( Ψ ) ) 1 .
Claim 4. On every bounded and open set Ω S 1 , W is Θ -compact. Define Ω = { ϕ S 1 : ϕ S 1 < M } as a bounded and open set, where M > 0 .
The proof of this claim will be done in three steps.
Step 1.  Φ W is continuous.
This property for Φ W is derived due to the imposed conditions on the nonlinear function g and the Lebesgue dominated convergence criterion, immediately.
Step 2.  Φ W ( Ω ¯ ) is bounded.
Now, for each ϕ Ω ¯ and for all t A j , we have
| Φ W ( ϕ ) ( t ) | u j ( T j T j 1 ) u j T j 1 T j ( T j s ) u j 1 | g ( s , ϕ ( s ) ) | d s u j ( T j T j 1 ) u j T j 1 T j ( T j s ) u j 1 | g ( s , ϕ ( s ) ) g ( s , 0 ) | d s + u j ( T j T j 1 ) u j T j 1 T j ( T j s ) u j 1 | g ( s , 0 ) | d s g * + u j ( T j T j 1 ) u j T j 1 T j ( T j s ) u j 1 s δ ( K | ϕ ( s ) | ) d s g * + M K T j 1 δ ,
by assuming g * = sup t A j | g ( t , 0 ) | . Thus,
Φ W ( ϕ ) E j g * + M K T j 1 δ : = R > 0 .
This shows that Φ W ( Ω ¯ ) E j is bounded.
Step 3.  Θ Ψ 1 ( I Φ ) W : Ω ¯ S 1 is completely continuous.
By considering the existing hypotheses in relation to Ascoli-Arzelà theorem, it is necessary that we prove two properties of the boundedness and equi-continuity for Θ Ψ 1 ( I Φ ) W ( Ω ¯ ) S 1 . At first, for each ϕ Ω ¯ and for all t A j , we have
Θ Ψ 1 ( I Φ ) W ϕ ( t ) = Θ Ψ 1 ( W ϕ ( t ) Φ W ϕ ( t ) ) = I T j 1 + u j g ( t , ϕ ( t ) ) u j ( T j T j 1 ) j u T j 1 T j ( T j s ) u j 1 g ( s , ϕ ( s ) ) d s = 1 Γ ( u j ) T j 1 t ( t s ) u j 1 g ( s , ϕ ( s ) ) d s t u j ( T j T j 1 ) u j Γ ( u j ) T j 1 T j ( T j s ) u j 1 g ( s , ϕ ( s ) ) d s .
Further, for each ϕ Ω ¯ and for all t A j , we get
| Θ Ψ 1 ( I Φ ) W ϕ ( t ) | 2 Γ ( u j ) T j 1 T j ( T j s ) u j 1 | g ( s , ϕ ( s ) ) g ( t , 0 ) | d s + 2 Γ ( u j ) T j 1 T j ( T j s ) u j 1 | g ( t , 0 ) | d s [ g * + M K T j 1 δ ] 2 ( T j T j 1 ) u j Γ ( u j + 1 ) : = B 1 .
so
Θ Ψ 1 ( I Φ ) W ϕ E j B 1 ,
which gives the uniform boundedness of Θ Ψ 1 ( I Φ ) W ( Ω ¯ ) in S 1 .
To prove the equi-continuity of Θ Ψ 1 ( I Φ ) W ( Ω ¯ ) , notice that for T j 1 t 1 t 2 T j and ϕ Ω ¯ , we get
| Θ Ψ 1 ( I Φ ) W ϕ ( t 2 ) Θ Ψ 1 ( I Φ ) W ϕ ( t 1 ) | g * + T j 1 δ M K Γ ( u j ) [ t 1 t 2 ( t 2 s ) u j 1 d s + T j 1 t 1 | ( t 2 s ) u j 1 ( t 1 s ) u j 1 | d s ] + T j 1 δ M K + g * Γ ( u j + 1 ) ( t 2 u j t 1 u j ) .
The right-hand side of above inequality tends to zero as t 1 t 2 . Thus, Θ Ψ 1 ( I Φ ) W ( Ω ¯ ) is equicontinuous in S 1 . On the basis of the Ascoli-Arzelà theorem, L Ψ 1 ( I Φ ) W ( Ω ¯ ) is relatively compact. In accordance with the steps 1 to 3, we can follow that W is Θ -compact in Ω ¯ .
Claim 5. There exists ϵ > 0 (not depending on λ ) so that if
Θ ( ϕ ) W ( ϕ ) = λ [ Θ ( ϕ ) + W ( ϕ ) ] , λ ( 0 , 1 ] ,
then ϕ S 1 ϵ . By the condition (AS2) and for each ϕ S 1 satisfying (12), we get
Θ ( ϕ ) W ( ϕ ) = λ Θ ( ϕ ) λ W ( ϕ ) .
So
Θ ( ϕ ) = 1 1 + λ W ( ϕ ) λ 1 + λ W ( ϕ ) .
By (13), and for all t A j , we get
ϕ ( t ) = 1 1 + λ Θ Ψ 1 W ϕ ( t ) λ 1 + λ Θ Ψ 1 W ( ϕ ( t ) ) ,
and so we estimate
| ϕ ( t ) | 1 ( 1 + λ ) Γ ( u j ) T j 1 t ( t s ) u j 1 g ( s , ϕ ( s ) ) g ( s , 0 ) d s + λ ( 1 + λ ) Γ ( u j ) T j 1 t ( t s ) u j 1 g ( s , ϕ ( s ) ) g ( s , 0 ) d s + g * ( T j T j 1 ) u j ( 1 + λ ) Γ ( u j + 1 ) + λ g * ( T j T j 1 ) u j ( 1 + λ ) Γ ( u j + 1 ) 1 1 + λ + λ 1 + λ T j 1 δ ( T j T j 1 ) u j Γ ( u j + 1 ) ( K ϕ E j ) + 1 1 + λ + λ 1 + λ g * ( T j T j 1 ) u j Γ ( u j + 1 ) = K T j 1 δ ( T j T j 1 ) u j Γ ( u j + 1 ) ϕ E j + g * ( T j T j 1 ) u j Γ ( u j + 1 ) .
Hence,
ϕ E j g * + K T j 1 δ ϕ E j ( T j T j 1 ) u j Γ ( u j + 1 ) ,
and so
ϕ S 1 g * Γ ( u j + 1 ) ( T j T j 1 ) u j K T j 1 δ : = ϵ .
Claim 6. There exists a bounded and open set Ω S 1 such that
Θ ( ϕ ) W ( ϕ ) λ [ Θ ( ϕ ) + W ( ϕ ) ] ,
for all ϕ Ω and all λ ( 0 , 1 ] .
By the condition (AS2) and Claim 5, there exits ϵ > 0 (independent of λ ) such that if ϕ solves
Θ ( ϕ ) W ( ϕ ) = λ [ Θ ( ϕ ) + W ( ϕ ) ] , λ ( 0 , 1 ] ,
then ϕ S 1 ϵ . Consequently, if
Ω = { ϕ S 1 : ϕ S 1 < B } ,
then from the condition (AS2), it is immediately obtained that the set Ω introduced by (15), is symmetric, 0 Ω , and S 1 Ω ¯ = Ω ¯ .
Furthermore, it is obtained that
Θ ( ϕ ) W ( ϕ ) λ [ Θ ( ϕ ) W ( ϕ ) ] ,
for all ϕ Ω = { ϕ S 1 : ϕ S 1 = B } and for all λ ( 0 , 1 ] , where B > ϵ . This together with Theorem 1 imply that the equivalent constant order resonance FBVP (6) has at least one solution, and this completes the proof. □
Now, we complete our deduction on the existence property for solutions of the given Caputo FBVP of variable order (1).
Theorem 3.
Let the conditions (AS1) and (AS2) be satisfied for all j N 1 n . Then, the given Caputo FBVP of variable order (1) has at least a solution in C ( A , R ) .
Proof. 
We know that for all j N 1 n , and according to Theorem 2, the equivalent constant order resonance FBVP (6) has at least one solution ϕ j ˜ E j . For each j N 1 n , and on the existing subintervals, define
ϕ j = 0 , t [ 0 , T j 1 ] , ϕ ˜ j , t A j .
In such a case, ϕ j C ( [ 0 , T j ] , R ) satisfies the integral equation (5) for t A j , which means that ϕ j ( 0 ) = 0 , ϕ j ( T j ) = ϕ ˜ j ( T j ) = 0 and satisfies (5) for t A j , j N 1 n . Therefore, the piecewise function
ϕ ( t ) = ϕ 1 ( t ) , t A 1 , ϕ 2 ( t ) , t A 2 , ϕ n ( t ) , t A n = [ 0 , T ] ,
is a solution to the given Caputo FBVP of variable order (1) in C ( A , R ) . □

4. Example

Example 2.
Consider the following FBVP (based on the given Caputo FBVP of variable order (1)) as follows
c D 0 + u ( t ) ϕ ( t ) = sin ϕ ( t ) ϕ ( t ) cos t 5 1 + t , t A : = [ 0 , 2 ] , ϕ ( 0 ) = ϕ ( 2 ) .
Let
g ( t , ϕ ) = sin ϕ ϕ cos t 5 1 + t , ( t , ϕ ) [ 0 , 2 ] × [ 0 , + ) ,
and
u ( t ) = 7 5 , t A 1 : = [ 0 , 1 ] , 3 2 , t A 2 : = [ 1 , 2 ] .
In this case,
t 1 2 | g ( t , ϕ 1 ) g 1 ( t , ϕ 2 ) | = t 1 2 ( sin ϕ 1 ϕ 1 cos t ) 5 1 + t t 1 2 ( sin ϕ 2 ϕ 2 cos t ) 5 1 + t 1 5 t 1 + t | sin ϕ 1 sin ϕ 2 | + | cos t | | ϕ 1 ϕ 2 | 2 5 | ϕ 1 ϕ 2 | .
By (17) and (6), on every subintervals A 1 and A 2 , two auxiliary constant order resonance FBVPs are considered as
c D 0 + 7 5 ϕ ( t ) = sin ϕ ( t ) ϕ ( t ) cos t 5 1 + t , t A 1 , ϕ ( 0 ) = ϕ ( 1 ) ,
and
c D 1 + 3 2 ϕ ( t ) = sin ϕ ( t ) ϕ ( t ) cos t 5 1 + t , t A 2 , ϕ ( 1 ) = ϕ ( 2 ) .
Evidently, the condition (AS2) is satisfied for j = 1 with δ = 1 2 and K = 2 5 , and
0 < K = 2 5 < min 1 , Γ ( u j + 1 ) ( T j T j 1 ) u j = min 1 , Γ ( 12 5 ) = 1 .
According to Theorem 2, the constant order resonance FBVP (18) has a solution like ϕ ˜ 1 E 1 .
Next, the condition (AS2) is also valid for j = 2 with δ = 1 2 and K = 2 5 , and
0 < K = 2 5 < min 1 , Γ ( u j + 1 ) ( T j T j 1 ) u j = min 1 , Γ ( 5 2 ) = 1 .
According to Theorem 2, the constant order resonance FBVP (19) has a solution like ϕ ˜ 2 E 2 .
Then, by Theorem 3, the given Caputo FBVP of variable order (16) has a solution as
ϕ ( t ) = ϕ ˜ 1 ( t ) , t A 1 , ϕ 2 ( t ) , t A 2 ,
where
ϕ 2 ( t ) = 0 , t A 1 , ϕ ˜ 2 ( t ) , t A 2 ,
and this shows the correctness of our results.

5. Conclusions

In this paper, a theoretical study was done for the given Caputo BVP of variable order (1) at resonance. To conduct this research, we defined some generalized subintervals as a partition of the main interval, and then on each subinterval, the piecewise constant functions were defined. With the help of these notions, we converted the given variable-order system to a constant-order system at resonance. In this case, we implemented the conditions of the Mawhin’s continuation theorem for proving the existence criterion for solutions of the corresponding BVP. Finally, an example was simulated numerically to show the correctness of our results. This technique on a variable-order BVP is new and determines the novelty of this work compared with other limited published papers in the form of variable orders. In relation to next studies, we aim to work on hidden-memory variable order systems and analyze the qualitative behaviors of their solutions such as existence, stability, and numerical solutions.

Author Contributions

Conceptualization, M.S.S. and Z.B.; Formal analysis, S.R., M.S.S., S.E., Z.B., S.K.N., S.A. and J.T.; Funding acquisition, J.T.; Methodology, S.R., M.S.S., S.E., Z.B., S.K.N. and S.A.; Software, S.E. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by King Mongkut’s University of Technology North Bangkok. Contract No. KMUTNB-62-KNOW-29.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data sharing not applicable to this article as no datasets were generated or analyzed during the current study.

Acknowledgments

The first and third authors would like to thank Azarbaijan Shahid Madani University.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
BVPBoundary Value Problem

References

  1. Gomez-Aguilar, J.F. Analytical and numerical solutions of a nonlinear alcoholism model via variable-order fractional differential equations. Phys. A Stat. Mech. Its Appl. 2018, 494, 52–75. [Google Scholar] [CrossRef]
  2. Almeida, R.; Tavares, D.; Torres, D.F.M. The Variable-Order Fractional Calculus of Variations; Springer: Cham, Switzerland, 2019. [Google Scholar]
  3. Sun, H.G.; Chang, A.; Zhang, Y.; Chen, W. A review on variable-order fractional differential equations: Mathematical foundations, Physical models, numerical methods and applications. Fract. Calc. Appl. Anal. 2019, 22, 27–59. [Google Scholar] [CrossRef] [Green Version]
  4. Abdeljawad, A.; Agarwal, R.P.; Karapinar, E.; Kumari, P.S. Solutions of the nonlinear integral equation and fractional differential equation using the technique of a fixed point with a numerical experiment in extended b-Metric Space. Symmetry 2019, 11, 686. [Google Scholar] [CrossRef] [Green Version]
  5. Abbas, M.I.; Ragusa, M.A. On the hybrid fractional differential equations with fractional proportional derivatives of a function with respect to a certain function. Symmetry 2021, 13, 264. [Google Scholar] [CrossRef]
  6. Boutiara, A.; Adjimi, N.; Benbachir, M.; Abdo, M.S. Analysis of a fractional boundary value problem involving Riesz-Caputo fractional derivative. Adv. Theory Nonlinear Anal. Its Appl. 2022, 6, 14–27. [Google Scholar]
  7. Ntouyas, S.K.; Etemad, S.; Tariboon, J. Existence results for multi-term fractional differential inclusions. Adv. Differ. Equ. 2015, 2015, 140. [Google Scholar] [CrossRef] [Green Version]
  8. Rezapour, S.; Ntouyas, S.K.; Iqbal, M.Q.; Hussain, A.; Etemad, S.; Tariboon, J. An analytical survey on the solutions of the generalized double-order φ-integrodifferential equation. J. Funct. Spaces 2021, 2021, 6667757. [Google Scholar] [CrossRef]
  9. Adiguzel, R.S.; Aksoy, U.; Karapinar, E.; Erhan, I.M. On the solution of a boundary value problem associated with a fractional differential equation. Math. Methods Appl. Sci. 2020. [Google Scholar] [CrossRef]
  10. Afshari, H.; Karapinar, E. A discussion on the existence of positive solutions of the boundary value problems via ψ-Hilfer fractional derivative on b-metric spaces. Adv. Differ. Equ. 2020, 2020, 616. [Google Scholar] [CrossRef]
  11. Etemad, S.; Rezapour, S.; Samei, M.E. α-ψ-contractions and solutions of a q-fractional differential inclusion with three-point boundary value conditions via computational results. Adv. Differ. Equ. 2020, 2020, 218. [Google Scholar] [CrossRef]
  12. Boutiara, A.; Etemad, S.; Hussain, A.; Rezapour, S. The generalized U-H and U-H stability and existence analysis of a coupled hybrid system of integro-differential IVPs involving φ-Caputo fractional operators. Adv. Differ. Equ. 2021, 2021, 95. [Google Scholar] [CrossRef]
  13. Zhang, S.; Hu, L. The existence of solutions and generalized Lyapunov-type inequalities to boundary value problems of differential equations of variable order. AIMS Math. 2020, 5, 2923–2943. [Google Scholar] [CrossRef]
  14. Sousa, J.V.d.C.; de Oliveira, E.C. Two new fractional derivatives of variable order with non-singular kernel and fractional differential equation. Comput. Appl. Math. 2018, 37, 5375–5394. [Google Scholar] [CrossRef]
  15. Benkerrouche, A.; Souid, M.S.; Etemad, S.; Hakem, A.; Agarwal, P.; Rezapour, S.; Ntouyas, S.K.; Tariboon, J. Qualitative study on solutions of a Hadamard variable order boundary problem via the Ulam-Hyers-Rassias stability. Fractal Fract. 2021, 5, 108. [Google Scholar] [CrossRef]
  16. Kaabar, M.K.A.; Refice, A.; Souid, M.S.; Martínez, F.; Etemad, S.; Siri, Z.; Rezapour, S. Existence and U-H-R stability of solutions to the implicit nonlinear FBVP in the variable order settings. Mathematics 2021, 9, 1693. [Google Scholar] [CrossRef]
  17. Zhang, S.; Hu, L. Unique existence result of approximate solution to initial value problem for fractional differential equation of variable order involving the derivative arguments on the half-axis. Mathematics 2019, 7, 286. [Google Scholar] [CrossRef] [Green Version]
  18. Zhang, S.; Li, S.; Hu, L. The existence and uniqueness result of solutions to initial value problems of nonlinear diffusion equations involving with the conformable variable derivative. Revista de la Real Academia de Ciencias Exactas Físicas y Naturales Serie A Matemáticas 2019, 113, 1601–1623. [Google Scholar] [CrossRef]
  19. Bouazza, Z.; Etemad, S.; Souid, M.S.; Rezapour, S.; Martínez, F.; Kaabar, M.K.A. A study on the solutions of a multiterm FBVP of variable order. J. Funct. Spaces 2021, 2021, 9939147. [Google Scholar] [CrossRef]
  20. Gaines, R.E.; Mawhin, J.L. Coincidence Degree and Nonlinear Differential Equations; Lecture Notes in Mathematics; Springer: Berlin/Heidelberg, Germany, 1977. [Google Scholar]
  21. Mawhin, J. Topological degree and boundary value problems for nonlinear differential equations. In Topological Methods for Ordinary Differential Equations; Furi, M., Zecca, P., Eds.; Lecture Notes in Mathematics; Springer: Berlin/Heidelberg, Germany, 1993; Volume 1537. [Google Scholar]
  22. Feng, W.; Webb, J.R.L. Solvability of three point boundary value problems at resonance. Nonlinear Anal. Theory Methods Appl. 1997, 30, 3227–3238. [Google Scholar] [CrossRef]
  23. Guezane-Lakoud, A.; Frioui, A. Third order boundary value problem with integral condition at resonance. Theory Appl. Math. Comput. Sci. 2013, 3, 56–64. [Google Scholar]
  24. Mawhin, J.; Ward, J.R., Jr. Periodic solutions of some forced Lienard differential equations at resonance. Archiv. Math. 1983, 41, 337–351. [Google Scholar] [CrossRef]
  25. Infante, G.; Zima, M. Positive solutions of multi-point boundary value problems at resonance. Nonlinear Anal. Theory Methods Appl. 2008, 69, 2458–2465. [Google Scholar] [CrossRef]
  26. Umarov, S.; Steinberg, S. Variable order differential equations with piecewise constant order-function and diffusion with changing modes. Z. Anal. Ihre Anwendungen 2009, 28, 431–450. [Google Scholar] [CrossRef]
  27. Zheng, X.; Wang, H. A hidden-memory variable-order time-fractional optimal control model: Analysis and approximation. SIAM J. Control Optim. 2021, 59, 1851–1880. [Google Scholar] [CrossRef]
  28. Lorenzo, C.F.; Hartley, T.T. Initialization, conceptualization, and application in the generalized (fractional) calculus. Crit. Rev. Biomed. Eng. 2007, 35, 447–553. [Google Scholar] [CrossRef] [PubMed]
  29. Samko, S.G. Fractional integration and differentiation of variable order. Anal. Math. 1995, 21, 213–236. [Google Scholar] [CrossRef]
  30. Samko, S.G.; Ross, B. Integration and differentiation to a variable fractional order. Integral Transform. Spec. Funct. 1993, 1, 277–300. [Google Scholar] [CrossRef]
  31. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differenatial Equations; North-Holland Mathematics Studies, 204; Elsevier Science B.V.: Amsterdam, The Netherlands, 2006. [Google Scholar]
  32. Podlubny, I. Fractional Differential Equations; Academic Press: San Diego, CA, USA, 1999; Volume 198. [Google Scholar]
  33. Zhang, S.; Sun, S.; Hu, L. Approximate solutions to initial value problem for differential equation of variable order. J. Fract. Calc. Appl. 2018, 9, 93–112. [Google Scholar]
  34. An, J.; Chen, P. Uniqueness of solutions to initial value problem of fractional differential equations of variable-order. Dyn. Syst. Appl. 2019, 28, 607–623. [Google Scholar]
  35. Zhang, S. The uniqueness result of solutions to initial value problems of differential equations of variable-order. Revista de la Real Academia de Ciencias Exactas Físicas y Naturales Serie A Matemáticas 2018, 112, 407–423. [Google Scholar] [CrossRef]
  36. Mawhin, J. Topological Degree Methods in Nonlinear Boundary Value Problems; CBMS Regional Conference Series in Mathematics; American Mathematical Society: Providence, RI, USA, 1979; Volume 40. [Google Scholar]
  37. O’Regan, D.; Cho, Y.J.; Chen, Y.Q. Topological Degree Theory and Application; Chapman and Hall/CRC: Boca Raton, FL, USA, 2006; Volume 10. [Google Scholar]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Rezapour, S.; Souid, M.S.; Etemad, S.; Bouazza, Z.; Ntouyas, S.K.; Asawasamrit, S.; Tariboon, J. Mawhin’s Continuation Technique for a Nonlinear BVP of Variable Order at Resonance via Piecewise Constant Functions. Fractal Fract. 2021, 5, 216. https://doi.org/10.3390/fractalfract5040216

AMA Style

Rezapour S, Souid MS, Etemad S, Bouazza Z, Ntouyas SK, Asawasamrit S, Tariboon J. Mawhin’s Continuation Technique for a Nonlinear BVP of Variable Order at Resonance via Piecewise Constant Functions. Fractal and Fractional. 2021; 5(4):216. https://doi.org/10.3390/fractalfract5040216

Chicago/Turabian Style

Rezapour, Shahram, Mohammed Said Souid, Sina Etemad, Zoubida Bouazza, Sotiris K. Ntouyas, Suphawat Asawasamrit, and Jessada Tariboon. 2021. "Mawhin’s Continuation Technique for a Nonlinear BVP of Variable Order at Resonance via Piecewise Constant Functions" Fractal and Fractional 5, no. 4: 216. https://doi.org/10.3390/fractalfract5040216

APA Style

Rezapour, S., Souid, M. S., Etemad, S., Bouazza, Z., Ntouyas, S. K., Asawasamrit, S., & Tariboon, J. (2021). Mawhin’s Continuation Technique for a Nonlinear BVP of Variable Order at Resonance via Piecewise Constant Functions. Fractal and Fractional, 5(4), 216. https://doi.org/10.3390/fractalfract5040216

Article Metrics

Back to TopTop