Next Article in Journal
Adapting to PSTs’ Pedagogical Changes in Sustainable Mathematics Education through Flipped E-Learning: Ranking Its Criteria with MCDA/F-DEMATEL
Previous Article in Journal
Outliers Detection Models in Shewhart Control Charts; an Application in Photolithography: A Semiconductor Manufacturing Industry
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A New Proof of the Existence of Nonzero Weak Solutions of Impulsive Fractional Boundary Value Problems

by
Asma Alharbi
1,†,
Rafik Guefaifia
2,† and
Salah Boulaaras
1,3,*,†
1
Department of Mathematics, College of Sciences and Arts, ArRass, Qassim University, Buraidah 51452, Saudi Arabia
2
Department of Mathematics, Faculty of Exact Sciences, University Tebessa 12002, Tebessa 12002, Algeria
3
Laboratory of Fundamental and Applied Mathematics of Oran (LMFAO), University of Oran 1, Ahmed Benbella, Oran 31000, Algeria
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Mathematics 2020, 8(5), 856; https://doi.org/10.3390/math8050856
Submission received: 15 April 2020 / Revised: 6 May 2020 / Accepted: 20 May 2020 / Published: 25 May 2020

Abstract

:
The paper deals with the existence of at least two non zero weak solutions to a new class of impulsive fractional boundary value problems via Brezis and Nirenberg’s Linking Theorem. Finally, an example is presented to illustrate our results.

1. Introduction

The paper deals with the existence of weak solutions to the following boundary value problems for impulsive fractional Hamiltonian differential equations:
t D T α 0 c D t α u t = F t , u , t 0 , T , t D T α 1 0 c D t α u i t j = 1 2 I i j u i t j , i A , j B u 0 = u T = 0 ,
where 0 α < 1 , 0 c D t α and t D T α denote the left Caputo fractional derivative and the right Riemann–Liouville fractional derivative of order α , respectively,
T > 0 , , A : = 1 , , N , B : = 1 , , L , u t = u 1 t , u 2 t , , u N t , 0 = t 0 < t 1 < t 2 < < t L < t L + 1 = T ,
F t , x denotes the gradient of F t , x in x, and Δ the operator is defined as
t D T α 1 0 c D t α u i t j = t D T α 1 0 c D t α u i t j + t D T α 1 0 c D t α u i t j ,
where
t D T α 1 0 c D t α u i t j + = lim t t j + t D T α 1 0 c D t α u i t and t D T α 1 0 c D t α u i t j = lim t t j t D T α 1 0 c D t α u i t
and 0 c D t α is the left Caputo fractional derivatives of order α . I i j : R R , i A , j B are continuous with F : 0 , T × R N R satisfies the following condition:
(A) F t , u is measurable according to t for any u R N and continuously differentiable according u for a.e. t 0 , T . However, there exist a C R + , R + and b L 1 0 , T , R + , such that
F t , u a u b t , F t , u a u b t
for any u R N and a.e. t 0 , T .
Fractional differential equations have recently proved to be valuable tools for modeling many phenomena in various fields of science and engineering. Indeed, one can find many applications in viscoelasticity, electrochemistry, control, porous media, electromagnetic, etc., for example, see [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19]. Moreover, recently, the existence of solutions of boundary value problems for Fractional differential equations have widely been studied in many papers and we refer the reader to the papers [3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18] and the references therein. For instance, in [3], the authors created a variable structure and, using the critical point theory, investigated the existence of multiple solutions for a class of fractional advection–dispersion equations derived from symmetric mass flow transmission. In [20], Torres consider the fractional Hamiltonian system given by
t D T α 0 D t α u t = F t , u t , t 0 , T , u 0 = u T = 0 ,
where α 1 2 , 1 , u R N , and F : 0 , T × R N R satisfies some conditions. In addition by using a modified version of mountain pass theorem for functional bounded from below due to Bonanno [21], the author studied the existence of at least three different solutions for problem (2) (see Theorems 1 and 2). In addition, under certain conditions and using some critical point theorems In [22] Zhou et al. proved the following fractional Hamiltonian system with impulsive effects has at least one weak solution:
t D T α 0 c D t α u t = F t , u t , t 0 , T , u 0 = u T = 0 , t D T α 1 0 c D t α u i t j = I i j u i t j , i = 1 , , N , j = 1 , , p .
The outline of the paper is as follows: In Section 2, we lay down preliminaries and assumptions, some of which will be needed in the body of the paper. Then, in Section 3, the main result is obtained, which gives the result of the existence of at least two non zero weak solutions to problem (1) via Brezis and Nirenberg’s Linking theorem (see Lemma 2 in [23]).

2. Preliminaries and Assumptions

In this section, we present some assumptions, definitions and basic properties of fractional computing, used later in this article. For omitted evidence, we refer the reader to ([24]) or other texts on the basic fractional calculation.
To state our main results, we set
ω : = 2 π T , ξ : = T α 1 2 Γ α 2 α 1 and t 0 , T ,
where Γ is the standard gamma function given by
Γ z = 0 + z α 1 e z d z .
Now, we assume that I i j s and F t , u satisfy the following assumptions:
h 1 There exists a constant γ > 0 , such that
lim u inf F t , u u 2 γ uniformly for a . e . t 0 , T .
h 2 There exist constants ρ > 0 , λ > 0 and k 1 verify
F t , u 1 2 λ k 2 ω 2 1 Γ 2 α 2 u 2 ,
for all u ρ and a.e. t 0 , T .
h 3 There exist constants γ i j [ 0.1 ) for any i A , j B , satisfying
γ i j 0 τ I i j s d s I i j τ τ , for every τ R .
h 4 I i j satisfy c i j : = inf τ = 1 0 τ I i j s d s > 0 .
h 5 There exist constants μ i j > 0 with an integer k 1 , which is defined in h 2 for all i A , j B , satisfying
0 < 0 z I i j s d s μ i j k 2 ω 2 z 2 , for every z R .
Theorem 1.
Suppose condition (A) is satisfied. Suppose that one of the following two conditions is true:
H 1 h 1 , h 3 and h 4 hold with γ > 1 2 ;
H 2 h 1 , h 3 and h 4 hold with
ξ 2 i = 1 N j = 1 L c i j 2 < γ 1 2 .
In addition, we assume also that the following condition holds:
H 3 h 2 , h 3 and h 5 hold with
λ 4 T 2 2 α 3 2 α + L N μ ξ 2 + 4 L k 2 ω 2 μ ξ 2 T 2 2 α 3 2 α ,
where μ = max i A , j B μ i j . Then, problem (1) has at least two nonzero weak solutions.
Definition 1.
([24]) Let u be a function defined on interval a , b . left and right Riemann-Liouville fractional integrals of order α > 0 for a function u denoted by a D t α u t and t D b α u t , respectively, are defined by
0 D t α u t = 1 Γ α 0 t t s α 1 u s d s , t a , b
and
t D b α u t = 1 Γ α t b s t α 1 u s d s , t a , b ,
Definition 2.
([24]) Let u be a function defined on a , b . The left and right Riemann–Liouville fractional derivatives of order α for a function u denoted by a D t α u t and t D b α u t are defined by
a D t α u t = d n d t n a D t α n u t = 1 Γ n α d n d t n a t t s n α 1 u s d s , t a , b
and
t D b α u t = 1 n d n d t n t D b α n u t = 1 Γ n α d n d t n t b s t n α 1 u s d s , t a , b ,
where n 1 α < 1 and n N . In particular, if 0 α < 1 , then
a D t α u t = d d t a D t α 1 u t = 1 Γ 1 α d d t a t t s α u s d s , t a , b
and
t D b α u t = d d t t D b α 1 u t = 1 Γ 1 α d d t t b s t α u s d s , t a , b
Definition 3.
([24]) Let α 0 and n N , we have:
i If α n 1 , 1 and u A C n a , b , R n , then the left and right Caputo fractional derivatives of order α for function u denoted by a c D t α u and t c D b α u respectively, are represented by
a c D t α u = 1 Γ n α a t t s n α 1 u n s d s , t a , b ,
and
t c D b α u = 1 n Γ n α t b t s n α 1 u n s d s , t a , b ,
respectively.
i i If α = n 1 and u A C n a , b , R n , then a c D t n 1 u t and t c D b n 1 u t are represented by a c D t n 1 u t = u n 1 t and t c D b n 1 u t = 1 n 1 u n 1 t .
Lemma 1.
([22], Proposition 3.2 ) Let 0 < α 1 , and 1 < p < . For all u L p 0 , T , R N . Then
0 D ξ α u L p 0 , t t α Γ α + 1 u L p 0 , t f o r ξ 0 , t , t 0 , T ,
where 0 D ξ α is left Riemann–Liouville fractional integral of order α .
Proposition 1.
([24]) We have the following property of fractional integration
0 T 0 D t α u t v t d t = 0 T t D b α v t u t d t , α > 0
provided that u L p 0 , T , R , v L q 0 , T , R and p 1 , 1 p + 1 q 1 + α or p 1 , q 1 , 1 p + 1 q = 1 + α .
Proposition 2.
([24]) Let n N and n 1 < α n . If u A C n a , b , R N or u C n a , b , R N , then
a D t α a c D t α u t = u t k = 0 n 1 u k a k ! t a k
t D b α t c D b α u t = u t k = 0 n 1 1 k u k b k ! b t k
and for t a , b . In particular, if 0 < α 1 and u A C a , b , R N or u C 1 a , b , R N , then
a D t α a c D t α u t = u t u a , a n d t D b α t c D b α u t = u t u b .
The Riemann - Liouville fractional derivative and the Caputo fractional derivative are related to one another by the following relationships.
Proposition 3.
([24]) Let n 1 α < n . for any n N . If u is a function defined on interval [ a , b ] for which the derivatives of Caputo fractional a c D t α u t and t c D b α u t of order α exist together with the Riemann–Liouville fractional derivatives a D t α u t and t D b α u t , thus
a c D t α u t = a D t α u t k = 0 n 1 u k a Γ k α + 1 t a k α
and
t c D b α u t = t D b α u t k = 0 n 1 u k b Γ k α + 1 b t k α ,
where t a , b . In particular, when 0 < α < 1 , we have
a c D t α u t = a D t α u t u a Γ 1 α t a α
and
t c D b α u t = t D b α u t u b Γ 1 α b t α .
Definition 4.
[25] Let 0 < α 1 and 1 < p < . The space of fractional derivative H 0 α , p is defined by the closure of C 0 0 , T , R N , that is
H 0 α , p = C 0 0 , T , R N ¯
with respect to the weighted norm
u α , p = 0 T u t p d t + 0 T 0 D t α u t p d t 1 p
for every u H 0 α , p .
Proposition 4.
([25], Proposition 3.1) Let 0 < α 1 , and 1 < p < . The fractional derivative space H 0 α , p is a reflexive and separable Banach space.
Proposition 5.
([25], Proposition 3.2) Let 0 < α 1 , and 1 < p < . For all u H 0 α , p , we have
u L p T α Γ α + 1 0 c D t α u t L p
Moreover, if α > 1 p , then
u T α 1 p Γ α α 1 q + 1 1 q 0 c D t α u t L p
Proposition 6.
([25], Proposition 3.3) Let 0 < α 1 and 1 < p < . Assume that α > 1 p and the sequence u k converges weakly to u in H 0 α , p , i.e., u k u . Then u k u in C 0 0 , T , R N , i.e., u u k 0 , as k .
In our study in this paper, we use Hilbert space H α = H 0 α , 2 with the inner product and the following corresponding norm
u , v = 0 T u t , v t d t + 0 T 0 c D t α u t , 0 c D t α v t d t , u , v H α
and
u α = u α , 2 = 0 T u t 2 d t + 0 T 0 c D t α u t 2 d t 1 2 , u H α
In view of (11), we have that, for t 0 , T and p = q = 2
u ξ 0 c D t α u L 2 ξ u α
Definition 5.
Let X be a Banach space with J : X R differentiable. It can be said that J satisfies the Palais–Smale P S -condition if for all sequence u n in X which J u n is bounded and J u n 0 as n possesses a convergent subsequence.
Definition 6.
Let X be a Banach space and J : X R differentiable and c R . It can be said that J satisfies the P S c -condition if the existence of a sequence u n in X such that
J u n c , J u n 0
as n , implies that c is a critical value of J .
Remark 1.
It can be remarked the P S -condition implies the P S c condition for each c R .
Lemma 2.
[23] Let X be a Banach space with a direct sum decomposition X = X 1 X 2 and k : = dim X 2 < . Let J C 1 X , R with J 0 = 0 , satisfying the P S -condition, we assume that, for ρ > 0 ,
J u 0 , u X 2 , u α ρ , J u 0 , u X 1 , u α ρ
In addition, we assume also that J is bounded below and inf X J < 0 . Thus J has at least two non zero critical points.

3. Mains Result

We establish the existence of at least two non zero weak solutions to problem (1) via Brezis and Nirenberg’s Linking Theorem.
Definition 7.
A function u H α is said to be a weak solution of problem (1) if the identity
0 T 0 c D t α u t 0 c D t α v t d t + 1 2 i = 1 N j = 1 L I i j u i t j v i t j d t + 0 T F t , u t , v t d t = 0
holds for any v H α .
Consider the functional J : H α R defined by
J u = 1 2 0 T 0 c D t α u t 2 d t + 0 T F t , u t d t + 1 2 i = 1 N j = 1 L 0 u i t j I i j s d s = J 1 u + J 2 u
where
J 1 u = 1 2 0 T 0 c D t α u t 2 d t + 0 T F t , u t d t
and
J 2 u = 1 2 i = 1 N j = 1 L 0 u i t j I i j s d t
Proposition 7.
The functional J is continuously differentiable on H α and
J u , v = 0 T 0 c D t α u t 0 c D t α v t d t + 1 2 i = 1 N j = 1 L I i j u i t j v i t j d t + 0 T F t , u t , v t d t
Proof. 
The proof of this proposition is very simple, we omit it. □
It is clear that, the critical points of J are weak solutions of (1). Now, we give the proof of our main results.
Proof of Theorem 1.
We apply Lemma 1 to J. Knowing that H α is a Banach space and J C 1 H α , R (see Proposition 5). By (15), it can be easily checked that functional J satisfies J 0 = 0 . We decompose the proof of the theorem into the following three steps.
Step1 We get h 1 or h 2 implies that
lim u α J u = + ,
and J u is bounded below on H α . To this end, for any 0 < ε < γ , h 1 implies that there exists δ > 0 such that
F t , u t γ ε u 2 , u R N with u δ and a . e . t 0 , T .
Let a δ = max u δ a u ; in view of (A), one can get
F t , u a u b t a δ b t + γ ε u 2 γ ε δ 2 , u R N with u δ and a . e . t 0 , T .
Then it follows from (20) and (21) that
F t , u γ ε u 2 γ ε δ 2 a δ b t , u R N with u δ and a . e . t 0 , T .
Thus, by (16) and (22), we have
J 1 u = 1 2 0 T 0 c D t α u t 2 d t + 0 T F t , u t d t 1 2 0 T 0 c D t α u t 2 d t + γ ε 0 T u t 2 d t γ ε δ 2 T a δ 0 T b t d t min 1 2 , γ ε u α 2 γ ε δ 2 T a δ 0 T b t d t
for all u H 0 α . For any z R , we set
K τ : = 0 τ 1 z I i j s d s τ γ i j , for all τ 1 .
By h 3
K τ : = τ γ i j 2 z I i j τ 1 z + γ i j 0 τ 1 z I i j s d s τ γ i j 1 = τ γ i j 1 γ i j 0 τ 1 z I i j s d s τ 1 z I i j τ 1 z 0
so that K is non-increasing in [ 1 , ) . Therefore, for any z 1 , we have K 1 K z , that is,
0 z I i j s d s 0 z 1 z I i j s d s z γ i j c i j z γ i j
where c i j = inf z = 1 0 z I i j s d s > 0 by h 4 Thus, it follows from (15), (16), (18), (23) and (24) that
J u min 1 2 , γ ε u α 2 γ ε δ 2 T a δ 0 T b t d t + 1 2 i = 1 N j = 1 L 0 u i t j I i j s d t min 1 2 , γ ε u α 2 γ ε δ 2 T a δ 0 T b t d t 1 2 i = 1 N j = 1 L c i j u i t j γ i j min 1 2 , γ ε u α 2 1 2 i = 1 N j = 1 L c i j u i t j γ i j γ ε δ 2 T a δ 0 T b t d t min 1 2 , γ ε u α 2 1 2 i = 1 N j = 1 L c i j ξ γ i j u i t j α γ i j γ ε δ 2 T a δ 0 T b t d t
Then in view of (25), for any γ > 0 , choosing 0 < ε < γ . The following two cases may occur.
Case1 γ > 1 2 choosing ε = γ 1 2 2 , we can get γ ε > 1 2 , so
min 1 2 , γ ε = 1 2 .
Thus, it follows from (25) that
J u 1 2 u α 2 1 2 i = 1 N j = 1 L c i j ξ γ i j u i t j α γ i j 1 2 γ + 1 2 δ 2 T a δ 0 T b t d t .
So ( H 1 ) yields that (19) holds, i.e.,
J u , a s u α .
Case2 γ 1 2 choosing ε = γ ξ 2 i = 1 N j = 1 L c i j 2 2 It follows from ( H 2 ) that
γ ε = 1 2 γ + 1 2 ξ 2 i = 1 N j = 1 L c i j > 1 2 ξ 2 i = 1 N j = 1 L c i j > 0
and ε > 0 . Since γ ε 1 2 ε < 1 2 , we have
min 1 2 , γ ε = γ ε .
Thus, it follows from (19) and (20) that
J u γ ε u α 2 1 2 i = 1 N j = 1 L c i j ξ γ i j u i t j α γ i j γ ε δ 2 T a δ 0 T b t d t γ ε u α 2 1 2 i = 1 N j = 1 L c i j ξ γ i j u i t j α γ i j 1 2 γ + 1 2 ξ 2 i = 1 N j = 1 L c i j δ 2 T a δ 0 T b t d t ,
which combined with ( H 2 ) yields that (19) holds, i.e.,
J u , a s u α .
Therefore, ( H 1 ) or ( H 2 ) implies (19). Thus, J u is bounded below on H α .
Step 2 We prove that ( H 1 ) or ( H 2 ) implies that J u satisfies the P S condition. Suppose that u n is a sequence in H α such that J u n is bounded and J u n 0 as n . Then u n is bounded on H 0 α . In fact, if u n is an unbounded sequence, without loss of generality we assume that u n α as n . By Step 1, we know that ( H 1 ) or ( H 2 ) implies (19). Thus J u n , which contradicts the boundedness of J u n . Since u n H α is bounded and H α is a reflexive Banach space and so by passing to a subsequence (for simplicity denoted again by u n ) if necessary, by Proposition 6, we may assume that
u n u , weakly in H α , u n u , strongly in C 0 , T , R N
By (18), we have
J u n J u , u n u = 0 T 0 c D t α u n t 0 c D t α u t 2 d t + 1 2 i = 1 N j = 1 L I i j u n i t j I i j u i t j u n i t j u i t j d t + 0 T F t , u n t F t , u t , u n t u t d t , 0 T 0 c D t α u n t 0 c D t α u t 2 d t + 1 2 i = 1 N j = 1 L I i j u n i t j I i j u i t j u n i t j u i t j d t + 0 T F t , u n t F t , u t , u n t u t d t ,
One has
0 T 0 c D t α u n t 0 c D t α u t 2 d t J u n J u , u n u 1 2 i = 1 N j = 1 L I i j u n i t j I i j u i t j u n i t j u i t j d t 0 T F t , u n t F t , u t , u n t u t d t
By (27), we know that
0 T u n t u t 2 d t 0 , as n
and for any i A , j B , we have that u n i t j u i t j , as n . In fact
u n i t j u i t j u n t j u t j for any i A , j B ,
Thus, it follows from the continuity of all I i j that
i = 1 N j = 1 L I i j u n i t j I i j u i t j u n i t j u i t j 0 , as n .
In view of ( A ) and (27), we have
0 T F t , u n t F t , u t , u n t u t d t 0 , as n .
Since J u n 0 , then using (28), (30) and (31), we have that
0 T 0 c D t α u n t 0 c D t α u t 2 d t 0 , as n .
By (29) and (32), one has u n u α 0 , which means that J u satisfies the ( P S ) condition.
Step 3 We prove that H 3 implies that (14) holds for some ρ > 0 . To this end, it follows from (13), (17) and h 5 , that
0 J 2 u = 1 2 i = 1 N j = 1 L 0 u i t j I i j s d t
1 2 i = 1 N j = 1 L k 2 ω 2 μ i j u i t j 2
1 2 i = 1 N j = 1 L k 2 ω 2 μ u 2
1 2 i = 1 N j = 1 L k 2 ω 2 μ ξ 2 u α 2
= 1 2 N L k 2 ω 2 μ ξ 2 u α 2
where μ = max i A , j B μ i j . Let
X 2 = m = 0 k a m sin m ω t ; a m R N ,
where a m = a m 1 , a m 2 , , a m N with a m i = 1 for i = 1 , , N . and X 1 be the orthogonal complement of X 2 in H α . Then dim X 2 < and H α = X 1 X 2 .
If ( H 3 ) holds, we will consider X 2 with k 1 . When u X 2 , we have
0 T u t 2 d t = 0 T i = 1 N u i t 2 d t = 0 T i = 1 N m = 0 k a m i sin m ω t 2 d t = 0 T i = 1 N m = 0 k a m i 2 sin 2 m ω t + m = 0 k j = 0 , j m k a m i a j i sin m ω t sin j ω t d t
On the other hand, one has
0 T sin m ω t sin j ω t d t = T 2 , m = j 0 , other wise .
So, in view of (34) and (35), we have
0 T u t 2 d t = T 2 i = 1 N m = 0 k a m i 2 = N k + 1 T 2
Also, For all u X 2 , and by Proposition 3, one has 0 c D t α u t = 0 D t α u t and
0 c D t α u i t = 0 c D t α u i t = 1 Γ 1 α d d t 0 T t s α m = 0 k a m i sin m ω s d s 1 Γ 1 α d d t 0 T t s α m = 0 k a m i m ω s d s = 1 Γ 1 α d d t m = 0 k a m i m ω t 2 α 1 α 2 α = 1 1 α Γ 1 α m = 0 k a m i m ω t 1 α = 1 Γ 2 α m = 0 k a m i m ω t 1 α
which implies
0 T 0 c D t α u t 2 d t 1 Γ 2 α 2 0 T i = 1 N m = 0 k a m i 2 m 2 ω 2 + m = 0 k j = 0 , j m k m j ω 2 a m i a j i t 2 1 α d t = 1 Γ 2 α 2 i = 1 N m = 0 k a m i 2 m 2 ω 2 + m = 0 k j = 0 , j m k m j ω 2 a m i a j i T 3 2 α 3 2 α 1 Γ 2 α 2 2 k + 1 k 2 ω 2 T 3 2 α 3 2 α
By (16) and h 2 , for all u α ρ , we have
J 1 u = 1 2 0 T 0 c D t α u t 2 d t + 0 T F t , u t d t 1 2 0 T 0 c D t α u t 2 d t 1 2 λ k 2 ω 2 1 Γ 2 α 2 0 T u t 2 d t .
By (15), (33), (36)–(38), we have
J u 1 2 0 T 0 c D t α u t 2 d t 1 2 λ k 2 ω 2 1 Γ 2 α 2 0 T u t 2 d t + 1 2 N L k 2 ω 2 μ ξ 2 u α 2 , 1 2 1 Γ 2 α 2 2 k + 1 k 2 ω 2 T 3 2 α 3 2 α 1 4 λ k 2 ω 2 N k + 1 T 1 Γ 2 α 2 + 1 2 N L k 2 ω 2 μ ξ 2 u α 2 ,
for all u X 2 . Therefore, by the fact 1 < 1 Γ 2 α 2 for α ( 1 2 , 1 ] , 2 k + 1 < 2 k + 2 , N > 1 , (36), (37) and (39), we have
J u 1 Γ 2 α 2 2 k + 1 k 2 ω 2 T 3 2 α 3 2 α 1 4 λ k 2 ω 2 N k + 1 T 1 Γ 2 α 2 + 1 4 L N 2 k + 1 k 2 ω 2 ξ 2 T + 1 Γ 2 α L N 2 k + 1 k 4 ω 4 μ ξ 2 T 3 2 α 3 2 α 1 Γ 2 α 2 2 k + 2 k 2 ω 2 T 3 2 α 3 2 α 1 4 λ k 2 ω 2 N k + 1 T 1 Γ 2 α 2 + 1 4 L N 2 k + 1 k 2 ω 2 μ ξ 2 T 1 Γ 2 α 2 + 1 Γ 2 α 2 L N 2 k + 2 k 4 ω 4 μ ξ 2 T 3 2 α 3 2 α 2 k 2 k + 1 π 2 T 1 Γ 2 α 2 × 2 T 2 2 α 3 2 α 1 2 λ N + 1 2 L N 2 μ ξ 2 + 2 L N k 2 ω 2 μ ξ 2 T 2 2 α 3 2 α
Thus
J u k 2 k + 1 N π 2 T 1 Γ 2 α 2 × 4 T 2 2 α 3 2 α λ + L N μ ξ 2 + 4 L k 2 ω 2 μ ξ 2 T 2 2 α 3 2 α
for all u X 2 with u α ρ . Thus, H 3 implies that J u 0 for all u X 2 with u α ρ .
In view of (22) by choosing ε = γ 2 , we get
J 1 u = 1 2 0 T 0 c D t α u t 2 d t + 0 T F t , u t d t 1 2 0 T 0 c D t α u t 2 d t + γ 2 0 T u t 2 d t γ 2 δ 2 T a δ 0 T b t d t
min 1 2 , γ 2 u α 2 γ 2 δ 2 T a δ 0 T b t d t ,
for all u H α . Since b L 1 0 , T , R + , then J u 0 for all u X 1 with u α ρ for ρ large enough. Thus, H 3 implies that (14) holds for ρ large enough.
Moreover, it follows from h 2 that F t , 0 = 0 for a.e. t 0 , T . Thus, h 2 implies J 0 = 0 . Now if inf H α J 0 , by Step 3 we have that all u X 2 with u α ρ are minima of J, which implies that J has infinitely many critical points. If inf H α J 0 , then it follows from Lemma 2 that J has at least two nonzero critical points. Therefore, problem (1) has at least two nonzero weak solutions in H α . The proof of Theorem 1 is complete. □
Example 1.
Let α = 3 4 , T = 2 π , N = 1 and k = 1 . Consider the boundary value problem of the fractional differential equation with impulsive effects
t D T 3 4 0 c D t 3 4 u t = F t , u , t 0.2 π , t D T 1 4 0 c D T 3 4 u i t j = 1 2 I i j u i t j , i = 1 , j = 1 u 0 = u 2 π = 0 ,
where I i j s = 1 5 0 . 0005 s . , Direct computation shows that ξ = 2 π 1 4 2 Γ 3 4 1.8 , ω = 1 , h 3 , h 4 and h 5 hold with γ i j = 1 2 and μ = 1 2000 .
Case 1Assume that F of (42) is
F t , u = 1 2 t + 2 u 2 u , u > 1 , 5.71187 t + 1 u 4 u 2 , u 1 ,
where t 0.2 π . In this case, we have that condition(A)holds with a u = max u 4 + u 2 , 2 u 2 + 2 u , 4 u 3 + 2 u , 4 u + 2 and b t = 7 t + 2 . Direct computation shows that h 1 holds with γ = 1 since, we have
lim u inf F t , u u 2 = 1 2 t + 2 γ > 1 2
Moreover, direct computation shows that h 2 holds with λ = 9.46 and
λ 9 . 45 + 24 . 89 μ 9 . 46 ,
where Γ 2 α = Γ 5 4 0.91 . Since lim u 0 F t , u u 2 = 5.71187 t + 1 , then by choosing ε = 1 100 such that F t , u 5.71187 u 2 .
Thus, by F t , 0 = 0 , H 1 and H 3 holds. According to Theorem 1, the problem (42) has at least two non zero weak solutions.
Case 2Assume that F of (42) is
F t , u = t + 1 3 u 2 u , u > 1 , 9.77 t + 0.605 u 4 u 2 , u 1 ,
where t 0.2 π . In this case, we have that condition(A)holds with a u = max u 4 + u 2 , 2 u 2 + 2 u , 4 u 3 + 2 u , 4 u + 2 and b t = 10 t + 1 . Direct computation shows that h 1 holds with γ = 1 3 , since
lim u inf F t , u u 2 = t + 1 3 γ ,
and
0.06593 ξ 2 i 1 j 1 c i j 2 < γ 1 2 .
Moreover, direct computation shows that ( h 2 ) holds. Since
lim u 0 F t , u u 2 = 9.77 t + 0.605 .
Therefore, by choosing ε = 132 1000 , we have F t , u 5.71187 u 2 . Then, by F t , 0 = 0 , H 2 and H 3 holds. According to Theorem 1, the problem (42) admits at least two nonzero weak solutions.

4. Conclusions

Fractional differential equations have recently proved to be valuable tools for modeling many phenomena in various fields of science and engineering. Indeed, one can find many applications in viscoelasticity, electrochemistry, control, porous media, electromagnetic, etc., for example, see [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19]. Moreover, recently, the existence of solutions of boundary value problems for fractional differential equations have widely been studied in many papers and we refer the reader to the papers [3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18]. In this work, we can extend the previous mentioned works for proving the existence of at least two non zero weak solutions to a new class of impulsive fractional boundary value problems via Brezis and Nirenberg’s Linking Theorem. Finally, an example is presented to illustrate our results. In the next work, we will try to prove the existence of three different weak solutions of the p-Laplacian fractional for an overdetermined nonlinear fractional partial Fredholm–Volterra integro-differential system by using variational methods combined with a critical point theorem due to Bonanno and Marano.

Author Contributions

Conceptualization, methodology and writing-original manuscript, A.A., R.G. and S.B.; supervision, S.B.; writing review and editing, S.B. and R.G., formal analysis, S.B. and A.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research received NO external or internal funding.

Acknowledgments

The authors would like to thank the both anonymous referees for their careful reading and for relevant remarks/suggestions which helped them to improve the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Agarwal, P.; Karimov, E.; Mamchuev, M.; Ruzhansky, M. On boundary-value problems for a partial differential equation with Caputo and Bessel operators. In Recent Applications of Harmonic Analysis to Function Spaces, Differential Equations, and Data Science; Birkhauser: Cham, Switzerland, 2017; Volume 152, pp. 707–718. [Google Scholar]
  2. Atangana, A.; Hammouch, Z. Fractional calculus with power law: The cradle of our ancestors. Eur. Phys. J. Plus 2017, 134, 429. [Google Scholar] [CrossRef]
  3. Bai, C. Existence of solutions for a nonlinear fractional boundary value problem via a local minimum theorem. Elect. J. Differ. Equat. 2012, 176, 1–9. [Google Scholar]
  4. Benson, D.; Wheatcraft, S.; Meerschaert, M. Application of a fractional advection-dispersion equation. Water Resour. Res. 2000, 36, 1403–1412. [Google Scholar] [CrossRef] [Green Version]
  5. Benson, D.; Wheatcraft, S.; Meerschaert, M. The fractional-order governing equation of Levy motion. Water Resour. Res. 2000, 36, 1413–1423. [Google Scholar] [CrossRef]
  6. Boulaaras, S.; Bouizem, Y.; Guefaifia, R. Existence of positive solutions of (p(x), q(x))-Laplacian parabolic systems with right hand side defined as a multiplication of two separate functions. Math. Methods Appl. Sci. 2020, 43, 2615–2625. [Google Scholar] [CrossRef]
  7. Boulaaras, S.; Allahem, A. Existence of Positive Solutions of Nonlocal p(x)-Kirchhoff Evolutionary Systems via Sub-Super Solutions Concept. Symmetry 2019, 11, 253. [Google Scholar] [CrossRef] [Green Version]
  8. Boulaaras, S.; Guefaifia, R. Existence of positive weak solutions for a class of Kirrchoff elliptic systems with multiple parameters. Math. Methods Appl. Sci. 2018, 41, 5203–5210. [Google Scholar] [CrossRef]
  9. Diethelm, K.; Freed, A. On the solution of nonlinear fractional order differential equations used in the modeling of viscoelasticity. In Scientific Computing in Chemical Engineering II-Computational Fluid Dynamics, Reaction Engineering and Molecular Properties; Keil, F., Mackens, W., Voss, H., Werther, J., Eds.; Springer: Berlin/Heidelberg, Germany, 1999; pp. 217–224. [Google Scholar]
  10. Devi, J.; Lakshmikantham, V. Nonsmooth analysis and fractional differential equations. Nonlinear Anal. TMA 2009, 70, 4151–4157. [Google Scholar] [CrossRef]
  11. Erwin, V.; Roop, J. Variational formulation for the stationary fractional advection dispersion equation. Numer. Methods Partial Differ. Equ. 2006, 22, 58–76. [Google Scholar]
  12. Fix, G.; Roop, J. Least squares finite-element solution of a fractional order two-point boundary value problem. Comput. Math. Appl. 2004, 48, 1017–1033. [Google Scholar] [CrossRef] [Green Version]
  13. Glockle, W.; Nonnenmacher, T. A fractional calculus approach of self-similar protein dynamics. Biophys. J. 1995, 68, 46–53. [Google Scholar] [CrossRef] [Green Version]
  14. Guefaifia, R.; Boulaaras, S. Sub-super solutions method for elliptic systems involving (p1,…,pm) Laplacian operator. Math. Methods Appl. Sci. 2020, 43, 4191–4199. [Google Scholar]
  15. Hilfer, R. Applications of Fractional Calculus in Physics; World Scientific: Singapore, 2000. [Google Scholar]
  16. Jarad, F.; Abdeljawad, T.; Hammouch, Z. On a class of ordinary differential equations in the frame of Atangana–Baleanu fractional derivative. Chaos Solitons Fractals 2018, 117, 16–20. [Google Scholar] [CrossRef]
  17. Kamache, F.; Guefaifia, R.; Boulaaras, S.; Alharbi, A. Existence of weak solutions for a new class of fractional p-Laplacian boundary value systems. Mathematics 2020, 8, 475. [Google Scholar] [CrossRef] [Green Version]
  18. Kirchner, J.; Feng, X.; Neal, C. Fractal streamchemistry and its implications for contaminant transport in catchments. Nature 2000, 403, 524–526. [Google Scholar] [CrossRef]
  19. Kıymaz, I.O.; Cetinkaya, A.; Agarwal, P. An extension of Caputo fractional derivative operator and its applications. J.Nonlinear Sci. Appl. 2016, 9, 3611–3621. [Google Scholar] [CrossRef]
  20. Torres, C.; Oliverio, D. Existence of three solution for fractional Hamiltonian system. J. Homepage 2017, 4, 51–58. [Google Scholar] [CrossRef]
  21. Zhang, S. Positive solutions to singular boundary value problem for nonlinear fractional differential equation. Comput. Math. Appl. 2010, 59, 1300–1309. [Google Scholar] [CrossRef] [Green Version]
  22. Zhou, J.; Wang, Y.; Li, Y. Existence of solutions for a class Of Fractional Hamiltonian Systems with impulsive effects. Fract. Differ. Calc. 2018, 8, 233–253. [Google Scholar] [CrossRef] [Green Version]
  23. Brézis, H.; Nirenberg, L. Remarks on finding critical points. Commun. Pure Appl. Math. 1991, 44, 939–963. [Google Scholar] [CrossRef]
  24. Podlubny, I. Fractional Differential Equations; Academic Press: San Diego, CA, USA, 1999. [Google Scholar]
  25. Rodriguez-Lopez, R.; Tersian, S. Multiple solutions to boundary value problem for impulsive fractional differential equations. Fract. Calc. Appl. Anal. 2014, 17, 1016–1038. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Alharbi, A.; Guefaifia, R.; Boulaaras, S. A New Proof of the Existence of Nonzero Weak Solutions of Impulsive Fractional Boundary Value Problems. Mathematics 2020, 8, 856. https://doi.org/10.3390/math8050856

AMA Style

Alharbi A, Guefaifia R, Boulaaras S. A New Proof of the Existence of Nonzero Weak Solutions of Impulsive Fractional Boundary Value Problems. Mathematics. 2020; 8(5):856. https://doi.org/10.3390/math8050856

Chicago/Turabian Style

Alharbi, Asma, Rafik Guefaifia, and Salah Boulaaras. 2020. "A New Proof of the Existence of Nonzero Weak Solutions of Impulsive Fractional Boundary Value Problems" Mathematics 8, no. 5: 856. https://doi.org/10.3390/math8050856

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