812
Views
6
CrossRef citations to date
0
Altmetric
Research Articles

Existence and exact asymptotic behaviour of positive solutions for fractional boundary value problem with P-Laplacian operator

& ORCID Icon
Pages 370-376 | Received 10 Dec 2018, Accepted 03 Feb 2019, Published online: 17 Feb 2019

ABSTRACT

This paper deals with existence, uniqueness and global behaviour of a positive solution for the fractional boundary value problem Dβ(ψ(x)Φp(Dαu))=a(x)uσ in (0,1) with the condition limx0Dβ1(ψ(x)Φp(Dαu(x)))=limx1ψ(x)Φp(Dαu(x))=0andlimx0Dα1u(x)=u(1)=0, where β,α(1,2], Φp(t)=t|t|p2, p>1, σ(1p,p1), the differential operator is taken in the Riemann–Liouville sense and ψ,a :(0,1)R are non-negative and continuous functions that may are singular at x=0 or x=1 and satisfies some appropriate conditions.

1. Introduction

In this paper, we consider the following nonlinear boundary value problem of fractional differential equation with p-Laplacian (1) Dβ(ψ(x)Φp(Dαu))=a(x)uσ, x(0,1),limx0Dβ1(ψ(x)Φp(Dαu(x)))=limx1ψ(x)Φp(Dαu(x))=0limx0Dα1u(x)=u(1)=0,(1) where α,β(1,2] and for tR, Φp(t):=t|t|p2, (p>1) and σ(1p,p1). The operator Dα (and Dβ) is the Riemann–Liouville derivative of order α (of order β). The functions ψ and a are positive and continuous in (0,1) that may are singular at x=0 or x=1 and satisfying some conditions detailed below. Then, we will adress the question of existence, uniqueness and exact asymptotic behaviour of positive solutions to problem (Equation1).

We believe that there are few papers dedicated to the study of fractional differential equations with p-Laplacian operator, see for instance [Citation1–11]. These types of problems arise in applied field such as turbulent flow of a gas in a porous medium, biophysics, plasma physics and material science.

This work is motivated by recent advances in the study of p-Laplacian fractional differential equations involving singular or sublinear nonlinearities with different boundary conditions. Namely, in [Citation12], Liu considered the fractional differential equation (2) Dβ(Φp(ψ(x)Dαu(x)))=f(x,u(x)), x(0,1),(2) where 0<α,β1, ψC((0,1)) and f is a non-negative function on (0,1]×R allowed to be singular at t=0. The author proved the existence of positive solution with fractional non-local integral boundary conditions.

Recently, in [Citation13], Mâagli et al. considered the following problem (3) Dαu(x)=a(x)uσ, x(0,1),limx0+Dα1u(x)=0, u(1)=0,(3) where 1<α2, 1<σ<1 and the function a is required to satisfy some assumptions related to K, the set of all Karamata functions L defined on (0,η], by L(t):=cexptηz(s)sds, for some η>1, where c>0 and z is a continuous function on [0,η], with z(0)=0. To describe the result of [Citation13] in more details, we need some notations.

  • For two non-negative functions f and g defined on a set S, the notation f(x)g(x), xS means that there exists c>0 such that (1/c)f(x)g(x)cf(x), for all xS.

  • C2α([0,1])={f|tt2αf(t)C([0,1])}.

In [Citation13], Mâagli et al. studied problem (Equation3) where a verifies the assumption (H0):

(H0): aC((0,1)) satisfying for each x (0,1), (4) a(x)xλ(1x)μL(x)L~(1x),(4) where λ+(2α)σ1, μα and L,L~K such that (5) 0ηL(t)tλ+(2α)σdt<and0ηL~(t)tμα+1dt<.(5) Based on the Schauder fixed-point theorem, the authors showed in [Citation13] the following result.

Theorem 1.1

Assume that a satisfies (H0). Then problem (Equation3) has a unique positive solution uC2α([0,1]) satisfying for x(0,1), (6) u(x)xα2(1x)min(1,((αμ)/(1σ)))ΠL~,μ,σ,α(1x),(6) where ΠL~,μ,σ,α defined on (0,1) by (7) ΠL~,μ,σ,α(t):=1,if μ<σ+α1,tηL~(s)sds1/(1σ),if μ=σ+α1,L~(t)1/(1σ),ifσ+α1<μ<α,0tL~(s)sds1/(1σ),if μ=α.(7)

The main goal of this paper is to improve and extend the above results on the boundary behaviour of solutions to problem (Equation1). More precisely, we consider the following assumption: (H1): a and ψ are non-negative functions in C((0,1)) satisfying for each x (0,1), (8) a(x)xλ(1x)μL1(x)L2(1x)andψ(x)(1x)(p1)α+βμ,(8) where λ+(2α)σ1, β1<μ<β, 3<β+p and L1,L2K such that (9) 0ηL1(t)tλ+(2α)σdt<and0η(L2(t))1p1tdt<.(9)

Theorem 1.2

Assume (H1), then problem (Equation1) has a unique positive solution uC2α([0,1]) satisfying for x(0,1), (10) u(x)xα201x(L2(t))1/(p1)tdt(p1)/(p1σ).(10)

The rest of the paper is as follows. In Section 2, we state some already known results on functions in K and some definitions and lemmas from fractional calculus theory. In Section 3, we present some necessary conditions to existence result and we prove our main result.

2. Preliminaries: Karamata properties and fractional calculus

In what follows, we are quoting without proof some fundamental properties of functions belonging to the class K collected from [Citation14, Citation15]. We recall that a function L defined on (0,η] belongs to the class K, if (11) L(t):=cexptηz(s)sds,(11) for some η>1, c>0 and z is a continuous function on [0,η] with z(0)=0.

Proposition 2.1

  1. A function L is in K if and only if L is a positive function in C1((0,η]) such that (12) limt0+tL(t)L(t)=0.(12)

  2. Let L1, L2K and pR. Then we have (13) L1+L2K,L1L2KandL1pK.(13)

  3. Let LK and ε>0. Then we have (14) limt0+tεL(t)=0.(14)

As a standard example of functions belonging to the class K (see [Citation15]), we give

Example 2.2

Let mN and η>0. Let (μ1,μ2,,μm)Rm and w be a sufficiently large positif real number such that the function (15) L(t)=1imlogiwtμi(15) is defined and positive on (0,η], where logit=log.logt (i×). Then we have LK.

Lemma 2.3

  1. Let L be a function in K. Then we have (16) limt0+L(t)tηL(s)sds=0.(16) In particular, (17) ttηL(s)sdsK.(17)

  2. If 0ηL(s)sds converges, then (18) limt0+L(t)0tL(s)sds=0.(18) In particular, (19) t0tL(s)sdsK.(19)

Applying Karamata's theorem, we get the following.

Lemma 2.4

Let γR and L be a function in K defined on (0,η] for some η>1. We have

  1. If γ>1, then 0ηsγL(s)ds converges and (20) 0tsγL(s)dst0+t1+γL(t)1+γ.(20)

  2. If γ<1, then 0ηsγL(s)ds diverges and (21) tηsγL(s)dst0+t1+γL(t)1+γ.(21)

Next, we give some definitions and fundamental facts of fractional calculus theory, which can be found in [Citation16, Citation17].

Definition 2.5

Let γ>0, the Riemann–Liouville fractional integral of order γ of a measurable function f:(0,)R is given by (22) Iγf(x)=1Γ(γ)0x(xt)γ1f(t)dt,(22) provided that the right-hand side is pointwise defined on (0,).

Here Γ is the Euler Gamma function.

Definition 2.6

The Riemann–Liouville derivative of order γ>0 of a measurable function

f: (0,)R is given by Dγf(x)=1Γ(nγ)(ddx)n0x(xt)nγ1f(t)dt=(ddx)nInγf(x), provided that the right-hand side is pointwise defined on (0,).

Here n=[γ]+1, where [γ] means the integer part of the number γ.

Lemma 2.7

Let α,γ>0, n be the smallest integer greater than or equal to α and fL1((0,1]). Then, we have

  1. IαIγf(x)=Iα+γf(x) for x[0,1] and α+γ1.

  2. DαIαf(x)=f(x) for a.e. x[0,1].

  3. Dαf(x)=0 if and only if f(x)=j=1ncjxαj, where (c1,,cn)Rn.

Since our approach is based on potential theory, we recall in the following some basic tools. For γ(1,2], we denote by Gγ(x,t) the Green's function for the boundary value problem (Equation3). From [Citation13], we have (23) Gγ(x,t)=1Γ(γ)[xγ2(1t)γ1((xt)+)γ1],(23) where x+=max{x,0}.

Proposition 2.8

see [Citation13]

Let 1<γ2 and f be a non-negative measurable function on (0,1). Then we have

  1. For x,t(0,1), (24) Gγ(x,t)xγ2(1t)γ2(1max(x,t)).(24)

  2. For x(0,1), Gγf(x):=01Gγ(x,t)f(t)dt< if and only if 01(1t)γ1f(t)dt<.

  3. If the map t(1t)γ1f(t) is continuous and integrable on (0,1), then Gγf is the unique solution in C2γ([0,1]) of the boundary value problem (25) Dγu(x)=f(x), x(0,1),limx0+Dγ1u(x)=0,u(1)=0.(25)

Lemma 2.9

[Citation13]

Let L~1,L~2K and let for x(0,1) (26) b(t)=tλ1(1t)μ1L~1(t)L~2(1t),(26) with λ11 and μ1β. Assume that (27) 0ηtλ1L~1(t)dt<and0ηtβ1μ1L~2(t)dt<.(27) Then we have for x(0,1), (28) x2βGβb(x)(1x)min(1,βμ1)ΠL~2,μ1,0,β(1x).(28)

Remark 2.1

We need to verify condition 0ηtλ1L~1(t)dt< and 0ηtβ1μ1L~2(t)dt< in Lemma 6, only if λ1=1 and μ1=β. This is due to Lemma 2.

3. Proof of main result

We begin this section by announcing and proving some propositions and lemmas, that will play a crucial role in the proof of our main result.

Let p>1 and let q=p/(p1). Then Φq(t)=t|t|(2p)/(p1) is the function inverse of Φp. We say easily that Φq is increasing monotone and multiplicative on (0,).

Lemma 3.1

Let α,β(1,2]. Let f be a non-negative function in B((0,1)) such that x(1x)β1f(x) and x(1x)α1Φq((1/ψ)Gβf)(x) are continuous and integrable on (0,1), then the boundary value problem (29) Dβ(ψ(x)Φp(Dαu(x)))=f(x), x(0,1),limx0Dβ1(ψ(x)Φp(Dαu(x)))=limx1ψ(x)Φp(Dαu)(x)=0,limx0Dα1u(x)=u(1)=0,(29) has a unique positive solution in C2α([0,1]) given by (30) u(x)=Gα(Φq1ψGβf)(x):=01Gα(x,t)Φq1ψ(t)01Gβ(t,s)f(s)dsdt.(30)

Proof.

Since x(1x)β1f(x) is continuous and integrable on (0,1), we deduce by Proposition 2 (iii) that for x(0,1), we have (31) ψ(x)Φp(Dαu)(x)=Gβf(x).(31) Thus (32) Dαu(x)=Φq1ψ(x)Gβf(x).(32) In addition, using the fact that x(1x)α1Φq((1/ψ(x))Gβf(x)) is continuous and integrable on (0,1), we conclude again by Proposition 2 (iii) that problem (Equation29) has a unique solution u in C2α([0,1]) given by (33) u(x)=Gα(Φq1ψGβf)(x),0<x<1.(33)

Here, below we provide a crucial property concerning continuity.

Lemma 3.2

Let ϕ and ψ be tow non-negative functions in B((0,1)) such that (34) 01(1t)α1Φq1ψ(t)Gβϕ(t)dt<.(34) Then the family (35) F={Sf:xx2αGαΦq1ψGβf(x);|f|ϕ}(35) is relatively compact in C([0,1]).

Proof.

Let fB((0,1)) such that |f(x)|ϕ(x) for all x(0,1).

Let x(0,1), we have Sf(x)=x2αGαΦq1ψGβf(x)=x2α01Gα(x,t)Φq1ψ(t)Gβf(t)dt. Using Proposition 2 (i), we obtain that (36) |Sf(x)|1Γ(α)01(1t)α1Φq1ψ(t)Gβϕ(t)dt<.(36) Thus F is uniformly bounded. Now, let us prove that F is equicontinuous in [0,1].

Let x,y(0,1), then we have |Sf(x)Sf(y)|=|x2αGαΦq1ψGβf(x)y2αGαΦq1ψGβf(y)|=1Γ(α)|01y2α((yt)+)α1x2α((xt)+)α1Φq1ψ(t)Gβf(t)dt|1Γ(α)01|y2α((yt)+)α1x2α((xt)+)α1|Φq1ψ(t)Gβϕ(t)dt. For every t(0,1), we have (37) |y2α((yt)+)α1x2α((xt)+)α1|0as|xy|0,(37) and (38) |y2α((yt)+)α1x2α((xt)+)α1|2(1t)α1.(38) Then we obtain by Lesbegue's theorem that (39) |Sf(x)Sf(y)|0 as |xy|0.(39) Now, let x(0,1), we have |Sf(x)1Γ(α)01(1t)α1Φq1ψ(t)Gβf(t)dt|=1Γ(α)|01x2α((xt)+)α1Φq1ψ(t)Gβf(t)dt|1Γ(α)01x2α((xt)+)α1Φq1ψ(t)Gβϕ(t)dt. Using the fact that for t(0,1), (40) x2α((xt)+)α10 as x0(40) and (41) 0x2α((xt)+)α1(1t)α1,(41) we get again by Lesbegue's theorem that (42) |Sf(x)1Γ(α)01(1t)α1Φq1ψ(t)Gβf(t)dt|0 as x0.(42) Furthermore, for x(0,1), we have (43) |Sf(x)|1Γ(α)01((1t)α1x2α((xt)+)α1)Φq1ψ(t)Gβϕ(t)dt.(43) By similar arguments as above, we deduce that (44) |Sf(x)|0 as x1.(44) Finally, we conclude that the familly F is equicontinuous in [0,1]. Hence, by Ascoli's theorem, we deduce that F is relatively compact in C([0,1]).

Proposition 3.3

Assume (H1) and suppose that there exist a non-negative function θ in C([0,1]) such that (45) 01(1t)β1w(t)dt<(45) and (46) x2αGαΦq1ψGβw(x)θ(1x),(46) where w(t)=a(t)t(α2)σθσ(1t) for t(0,1).

Then problem  (Equation1) has a unique solution uC2α([0,1]) satisfying for each x(0,1) (47) u(x)xα2θ(1x).(47)

Proof.

Let m1 and θ be a non-negative function satisfying for each x(0,1) (48) 1mθ(1x)x2αGαΦq1ψGβw(x)mθ(1x).(48) Put c0:=m(p1)/(p1|σ|). We consider the closed convex set given by (49) Y:={vC([0,1]);1c0θ(1x)v(x)c0θ(1x)}.(49) Using Lemma 3.2 and Proposition 2 (ii), we easily see that the function xx2αGα(Φq((1/ψ)Gβw))(x) belongs to C([0,1]) and satisfies (Equation48). So Y is not empty. In order to use the Schauder's fixed point theorem, we denote a~(x)=x(α2)σa(x) and we define the operator T on Y by (50) Tv(x)=x2αGαΦq1ψGβ(a~vσ)(x).(50) We need to check that the operator T has a fixed point v in Y. For this choice of c0, we will prove that T map Y into itself. Indeed, let vY, by using (Equation48), we have Tv(x)Φq(c0|σ|)x2αGαΦq1ψGβw(x)c0|σ|/p1mθ(1x)=c0θ(1x) and Tv(x)Φq(c0|σ|)x2αGαΦq1ψGβw(x)c0|σ|/p1mθ(1x)=1c0θ(1x). Furthermore, using (Equation48), we have for each x(0,1) (51) GαΦq1ψGβw(x)<.(51) This implies by Proposition 2 (ii) that (52) 01(1t)α1Φq1ψ(t)Gβw(t)dt<.(52) Hence, it follows from Lemma 3.2 that the family TY is relatively compact in C([0,1]). So Y is invariant under T.

Next, we shall prove the continuity of T. Let {vk} be a sequence in Y which converges uniformly to v in Y.

For x(0,1), we have |Tvk(x)Tv(x)|=x2α|GαΦq1ψGβ(a~vkσ)(x)GαΦq1ψGβ(a~vσ)(x)|x2α01Gα(x,t)|Φq1ψ(t)Gβ(a~vkσ)(t)Φq1ψ(t)Gβ(a~vσ)(t)|dt. For t(0,1), we have (53) |Gβ(a~vkσ)(t)Gβ(a~vσ)(t)|01Gβ(t,s)|(a~vkσ)(s)(a~vσ)(s)|ds(53) and for every s(0,1), (54) |(a~vkσ)(s)(a~vσ)(s)|2c0|σ|w(s).(54) Using Proposition 2 (ii) and (Equation45), we obtain by Lesbegue's theorem that (55) |Gβ(a~vkσ)(t)Gβ(a~vσ)(t)|0 as k.(55) Since Φq is continuous, we deduce that (56) |Φq1ψ(t)Gβ(a~vkσ)(t)Φq1ψ(t)Gβ(a~vσ)(t)|0 as k.(56) We have (57) |Φq1ψ(t)Gβ(a~vkσ)(t)Φq1ψ(t)Gβ(a~vσ(t))|2Φq(c0|σ|)Φq1ψ(t)Gβw(t).(57) Using (Equation48), we obtain by Lesbegue's theorem that for x[0,1] (58) Tvk(x)Tv(x),ask.(58) Since TY is a relatively compact in C([0,1]), we have the uniform convergence, namely (59) TvkTv0,ask.(59) Thus, we have proved that T is a compact mapping from Y into itself. It follows by the Schauder fixed point theorem that there exists vY such that Tv=v. Put (60) u(x)=xα2v(x),(60) then uC2α([0,1]) and u satisfies the equation (61) u(x)=GαΦq1ψGβ(auσ)(x).(61) Then due to Lemma 3.1, u is a positive continuous solution of problem (Equation1).

Finally, let us prove that u is the unique positive continuous solution satisfying (Equation47). To this aim, we assume that (1.1) has two positive solutions u and v satisfying (Equation47). Then, there exists a constant m>1 such that (62) 1mvumv.(62) This implies that the set (63) J:=t1,:1tvutv(63) is not empty. Now, put c:=infJ, then we aim to show that c=1. Suppose that c>1. Then by simple calculus, we obtain that (64) Dβψ(x)(Φp(Dα(c(|σ|/(p1))v))Φp(Dαu))=a(c|σ|vσuσ)0,limx0Dβ1ψ(x)(Φp(Dα(c(|σ|/(p1))v))Φp(Dαu))(x)=0,limx1ψ(x)Φp(Dα(c(|σ|/(p1))v)(x))Φp(Dαu(x))=0.(64) We conclude by Proposition 2 (iii) that (65) ψΦp(Dα(c|σ|/p1v))Φp(Dαu)=Gβ(a(c|σ|vσuσ))0.(65) Then we have (66) Φp(Dα(c|σ|/p1v))Φp(Dαu).(66) Which implies with the fact that Φp is increasing monotone that (67) Dα(c|σ|/p1vu)0,limx0Dα1(c|σ|/p1vu)(x)=0,c|σ|/p1v(1)u(1)=0.(67) Using again Proposition 2 (iii), we conclude that (68) c|σ|/p1vu0.(68) By symmetry, we obtain that vc|σ|/p1u. So c|σ|/p1J. Since (|σ|/(p1))<1 and c>1, we have c|σ|/p1<c. This yields to a contradiction with the fact that c=infJ. Hence c=1 and consequently u=v.

3.1. Proof of Theorem  2

Suppose that β1<μ<β and L2 satisfy (69) 0η(L2(t))(1/(p1))tdt<.(69) Let θ be the function defined on [0,1] by (70) θ(x)=0x(L2(t))(1/(p1))tdt(p1)/(p1σ).(70) Then, for x(0,1) w(x):=a(x)x(α2)σθσ(1x)xλ+(α2)σ(1x)μL1(x)L2(1x)01x(L2(t))(1/(p1))tdtσ(p1)/(p1σ). We conclude by Lemma 2.9 that for x(0,1) (71) Gβw(x)xβ2(1x)βμL2(1x)01x(L2(t))(1/(p1))tdtσ(p1)/(p1σ).(71) This implies by Proposition 2 (ii) that 01(1t)β1w(t)dt< and for x(0,1), we have (72) Φq1ψGβw(x)x((2β)/(p1))(1x)α(L2(1x))1/(p1)01x(L2(t))(1/(p1))tdtσ/(p1σ).(72) Using again (H1) and by Lemma 2.9, we deduce that for x(0,1) (73) x2αGαΦq1ψGβw(x)θ(1x).(73) Hence it follows from Proposition 3 that problem (Equation1) has a unique positive solution u in C2α([0,1]) satisfying for x(0,1), (74) u(x)xα2θ(1x).(74) As an application of our main result, we give the following example.

3.2. Example 1

Let β,α(1,2], p>1 and σ(1p,p1). Let a and ψ be two positive continuous functions on (0,1) such that (75) a1(x)(1x)μlog31xp(75) and (76) ψ(x)(1x)(p1)α+βμ,(76) where β1<μ<β and β+p>3. Therefore by Theorem 2, problem (Equation1) has a unique positive solution uC2α([0,1]) satisfying for each x(0,1) (77) u(x)xα2log31x1/(p1σ).(77)

Acknowledgments

The authors would like to thank the anonymous reviewers for their valuable comments and suggestions to improve the manuscript.

Disclosure statement

No potential conflict of interest was reported by the authors.

Additional information

Funding

This work was supported by Taibah University.

References

  • Chai G. Positive solutions for boundary value problem of fractional differential equation with p-Laplacian operator. Bound Value Probl. 2012;2012(18):20.
  • Chen T, Liu W, Hu Z. A boundary value problem for fractional differential equation with p-Laplacian operator at resonance. Nonlinear Anal. 2012;75(6):3210–3217. doi: 10.1016/j.na.2011.12.020
  • Liu X, Jia M, Ge W. Multiple solutions of a p-Laplacian model involving a fractional derivative. Adv Differ Equ. 2013;2013(126):12. doi: 10.1186/1687-1847-2013-12
  • Liu X, Jia M, Geb W. The method of lower and upper solutions for mixed fractional four-point boundary value problem with p-Laplacian operator. Appl Math Lett. 2017;65:56–62. doi: 10.1016/j.aml.2016.10.001
  • Liu X, Jia M, Xiang X. On the solvability of a fractional differential equation model involving the p-Laplacian operator. Comput Math Appl. 2012;64:3267–3275. doi: 10.1016/j.camwa.2012.03.001
  • Lu H, Han Z, Sun S, et al. Existence of positive solutions for boundary value broblems of nonlinear fractional differential equations with p-Laplacian. Adv Differ Equ. 2013;2013(30):16.
  • Mahmudov N, Unul S. Existence of solutions of fractional boundary value problems with p-Laplacian operator. Bound Value Probl. 2015;2015(99):16.
  • Prasad KR, Krushna BMB. Existence of multiple positive solutions for p-Laplacian fractional order boundary value problems. Int J Anal Appl. 2014;6(1):63–81.
  • Wang J, Xiang H, Liu Z. Positive solutions for three-point boundary value problems of nonlinear fractional differential equations with p-Laplacian operator. Far East J Appl Math. 2009;37:33–47.
  • Wang J, Xiang H. Uper and lower solutions method for a class of singular fractional boundary value problems with p-Laplacian. Abstr Appl Anal Article ID. 2010;971824:12.
  • Bahaa GM, Hamiaz A. Optimal control problem for coupled time-fractional diffusion systems with final observations. J Taibah University Sci. 2019;13(1):124–135. doi: 10.1080/16583655.2018.1545560
  • Liu Y. Existence of solutions of periodic-type boundary value problems for multi-term fractional differential equations. Math Method Appl Sci. 2013;36:2187–2207. doi: 10.1002/mma.2746
  • Mâagli H, Mhadhebi N, Zeddini N. Existence and exact asymptotic behavior of positive solutions for a fractional boundary value problem. Abstr Appl Anal Article ID. 2013;420514:6.
  • Chemmam R, Mâagli H, Masmoudi S, et al. Combined effects in nonlinear singular elliptic problems in a bounded domain. Adv Nonlinear Anal. 2012;1:301–318.
  • Seneta R. Regular varying functions, Lecture notes in math. Vol. 508. Berlin: Springer-Verlag; 1976.
  • Kilbas AA, Srivastava HM, Trujillo JJ. Theory and applications of fractional differential equations. Amsterdam: Elsevier; 2006.
  • Podlubny I. Geometric and physical interpretation of fractional integration and fractional differentiation. Fract Calc Appl Anal. 2002;5:367–386.