Abstract

The existence of at least one positive solution is established for a class of semipositone fractional differential equations with Riemann-Stieltjes integral boundary condition. The technical approach is mainly based on the fixed-point theory in a cone.

1. Introduction

In this paper, we discuss the existence of positive solutions for the following singular semipositone fractional differential equation with nonlocal condition: -?????????(??)=????,??(??),-???????????(??)+??(??),???(0,1),??????(0)=??????+1??(0)=0,?????????(1)=10????????(??)????(??),(1.1) where 2<??=3, 0<??<1, and ??-??>2, ???? is the standard Riemann-Liouville derivative. ?10????????(??)????(??) denotes the Riemann-Stieltjes integral, and ?? is a function of bounded variation. ??:(0,1)×[0,+8)×(-8,0]?[0,+8) is continuous, ??:(0,1)?(-8,+8) is Lebesgue integrable.

Differential equations of fractional order have been recently proved to be valuable tools in the modeling of many phenomena in various fields of science and engineering. Indeed, we can find numerous applications in physics, engineering (like traffic, transportation, logistics, etc.), mechanics, chemistry, and so forth, (see [15]). There has been a significant development in the study of fractional differential equations in recent years, see the monographs of Kilbas et al. [6], Lakshmikantham et al. [7], Podlubny [4], Samko et al. [8], the survey by Agarwal et al. [9] and some recent results [1014].

On the other hand, the nonlocal condition given by a Riemann-Stieltjes integral is due to Webb and Infante in [1517] and gives a unified approach to many BVPs in [15, 16]. Motivated by [1517], Hao et al. [18] studied the existence of positive solutions for ??th-order singular nonlocal boundary value problem: ??(??)??(??)+??(??)??(??,??(??))=0,0<??<1,(??)?(0)=0,0=??=??-2,??(1)=10??(??)????(??),(1.2) where ?? can be singular at ??=0,1, ?? also can be singular at ??=0, but there is no singularity at ??=0,1. The existence of positive solutions of the BVP (1.2) is obtained by means of the fixed point index theory in cones.

More recently, Zhang [19] considered the following problem whose nonlinear term and boundary condition contain integer order derivatives of unknown functions ???????(??)+??(??)????,???,,??(??-2)???=0,0<??<1,??-1<??=??,(0)=???(0)=?=??(??-2)(0)=??(??-2)(1)=0,(1.3) where ???? is the standard Riemann-Liouville fractional derivative of order ??, ?? may be singular at ??=0 and ?? may be singular at ??=0, ??'=0,,??(??-2)=0. By using fixed point theorem of the mixed monotone operator, the unique existence result of positive solution to problem (1.3) was established. And then, Goodrich [20] was concerned with a partial extension of the problem (1.3) by extending boundary conditions -????????(??)=??(??,??(??)),0<??<1,??-1<??=??,(??)(0)=0,0=??=??-2,??????(1)=0,1=??=??-2,(1.4) and the author derived the Green’s function for the problem (1.4) and showed that it satisfies certain properties, then by using cone theoretic techniques, a general existence theorem for (1.4) was obtained when ??(??,??) satisfies some growth conditions.

Recently, Rehman and Khan [21] investigated the multipoint boundary value problems for fractional differential equations of the form: ?????????(??)=????,??(??),?????????(??),???(0,1),??(0)=0,????????(1)-??-2???=1??????????????????=??0,(1.5) where 1<??=2, 0<??<1, 0<????<1, ?????[0,+8) with ???-2??=1??????????-??-1<1. By using the Schauder fixed point theorem and the contraction mapping principle, the authors established the existence and uniqueness of nontrivial solutions for the BVP (1.5) provided that the nonlinear function ??:[0,1]×R×R is continuous and satisfies certain growth conditions. Since ?10????????(??)????(??) covers the multipoint BVP and integral BVP as special case, the fractional differential equations with the Riemann-Stieltjes integral condition also were extensively studied by many authors, see [22, 23]. In [23], Zhang and Han considered the existence of positive solution of the following singular fractional differential equation: ??????????(??)+??(??,??(??))=0,0<??<1,??-1<??=??,(??)?(0)=0,0=??=??-2,??(1)=10??(??)????(??),(1.6) where ??=2 and ????(??) can be a signed measure. Some growth conditions were adopted to guarantee that (1.6) has an unique positive solution, moreover, the authors also gave the iterative sequence of the solution, an error estimation, and the convergence rate of the positive solution.

Fractional differential equations like (1.1), with nonlinearities which are allowed to change sign and boundary conditions which contain nonlocal condition given by a Riemann-Stieltjes integral with a signed measure, are rarely studied. This type problems are referred to as semipositone problems in the literature, which arise naturally in chemical reactor theory [24]. In the recent work [25], by constructing a modified function, Zhang and Liu studied the existence of positive solution of a class of semipositone singular second-order Dirichlet boundary value problem, and when ?? is superlinear, a sufficient condition for the existence of positive solution is obtained under the simple assumptions.

Motivated by the above work, in this paper, we establish the existence of positive solutions for the semipositone fractional differential equations (1.1) when ?? is superlinear and involves fractional derivatives of unknown functions.

2. Preliminaries and Lemmas

Definition 2.1 (see [4, 5]). The Riemann-Liouville fractional integral of order ??>0 of a function ??:(0,+8)?R is given by ????1??(??)=?G(??)??0(??-??)??-1??(??)????,(2.1) provided that the right-hand side is pointwise defined on (0,+8).

Definition 2.2 (see [4, 5]). The Riemann-Liouville fractional derivative of order ??>0 of a function ??:(0,+8)?R is given by ??????1??(??)=???G(??-??)??????????0(??-??)??-??-1??(??)????,(2.2) where ??=[??]+1, [??] denotes the integer part of number ??, provided that the right-hand side is pointwise defined on (0,+8).

Remark 2.3. If ??,??:(0,+8)?R with order ??>0, then ??????(??(??)+??(??))=????????(??)+????????(??).(2.3)

Proposition 2.4 (see [4, 5]). (1) If ?????1(0,1),??>??>0, then ??????????(??)=????+????(??),????????????(??)=????-????(??),????????????(??)=??(??).(2.4)
(2) If ??>0,??>0, then ??????????-1=G(??)??G(??-??)??-??-1.(2.5)

Proposition 2.5 (see [4, 5]). Let ??>0, and ??(??) is integrable, then ????????????(??)=??(??)+??1????-1+??2????-2+?+????????-??,(2.6) where ?????R(??=1,2,,??), ?? is the smallest integer greater than or equal to ??.
Let ??(??)=??????(??),??(??)???[0,1], by standard discuss, one easily reduces the BVP (1.1) to the following modified problems, -??????-?????(??)=????,???????(??),-??(??)+??(??),??(0)=????(0)=0,??(1)=10??(??)????(??),(2.7) and the BVP (2.7) is equivalent to the BVP (1.1).

Lemma 2.6 (see [26]). Given ?????1(0,1), then the problem, ??????-????(??)+??(??)=0,0<??<1,??(0)=???(0)=0,??(1)=0,(2.8) has the unique solution ???(??)=10??(??,??)??(??)????,(2.9) where ??(??,??) is the Green function of the BVP (2.8) and is given by 1??(??,??)=?[]G(??-??)??(1-??)??-??-1[],0=??=??=1,??(1-??)??-??-1-(??-??)??-??-1,0=??=??=1.(2.10)

Lemma 2.7 (see [26]). For any ??,???[0,1], ??(??,??) satisfies: ????-??-1(1-??)??(1-??)??-??-1G=(??-??)=??(??,??)??(1-??)??-??-1G???(??-??-1),or??-??-1(1-??)G?.(??-??-1)(2.11)

By Lemma 2.6, the unique solution of the problem, ????????(??)=0,0<??<1,??(0)=???(0)=0,??(1)=1,(2.12) is ????-??-1. Let ???=10????-??-1?????(??),B=10????-??-1(1-??)????(??)(2.13) and define ?????(??)=10??(??,??)????(??),(2.14) as in [25], one can get that the Green function for the nonlocal BVP (2.7) is given by ????(??,??)=??-??-1??1-????(??)+??(??,??).(2.15) Throughout paper one always assumes the following holds.

(????)?? is a increasing function of bounded variation such that ????(??)=0 for ???[0,1] and 0=??<1, where ?? is defined by (2.13).

Define ??+(??)=max{??(??),0},??-(??)=max{-??(??),0}.(2.16) One has the following Lemma.

Lemma 2.8. Let 1<??-??=2 and (????) hold, then the unique solution ??(??) of the linear problem, -??????-????(??)=??-??(??),???(0,1),(0)=????(0)=0,??(1)=10??(??)????(??),(2.17) satisfies ??(??)=????(??),(2.18) where 1??=???G(??-??-1)1+B??10??-(??)????,??(??)=????-??-1B(1-??)+??1-????-??-1.(2.19)

Proof. By (2.11) and that ??(??) is a increasing function of bounded variation, we have ?????(??)=10???(??,??)????(??)=10????-??-1??G(??-??-1)????(??)=.G(??-??-1)(2.20) Consequently, ???(??)=10??(??,??)??-?(??)????=10?????-??-1??1-???????(??)+??(??,??)-=??(??)??????-??-1(1-??)?G(??-??-1)10??-(??)????+??????-??-1?G(??-??-1)(1-??)10??-=1(??)???????G(??-??-1)1+B??????-??-1B(1-??)+??1-????-??-1??10??-(??)????=????(??).(2.21)

Remark 2.9. (1) ??(??) satisfies (1-??)??(??)2=(1-??)??(??)??-??-1=1,(2.22)
(2) ????(??,??)=??(??)+1-????(1-??)??-??-1G(??-??-1)=F(??).(2.23)
In fact, since 1<??-??-1<2, the left side of (1) clearly holds. For right side of (1), from B=??, one gets B??(??)=1+??1-??=1+=11-??,1-??(2.24) thus we have (1-??)??(??)??-??-1=1.(2.25) (2) is obvious from (2.11).
Now define a function [·]* for any ?????[0,1] by []??(??)*=???(??),??(??)=0,0,??(??)<0(2.26) and consider the following approximate problem of the BVP (2.7): -??????-?????(??)=????,????[]??(??)-??(??)*[],-??(??)-??(??)*?+??+(??),??(0)=????(0)=0,??(1)=10??(??)????(??).(2.27)

Lemma 2.10. Suppose ?? is a positive solution of the problem (2.27) and satisfies ??(??)=??(??),???[0,1], then ??-?? is a positive solution of the problem (2.7), consequently, ????[??(??)-??(??)] also is a positive solution of the BVP (1.1).

Proof. In fact, if ?? is a positive solution of the BVP (2.27) such that ??(??)=??(??) for any ???[0,1], then, from (2.27) and the definition of [??(??)]*, we have -??????-?????(??)=????,?????(??(??)-??(??)),-(??(??)-??(??))+??+(??),??(0)=????(0)=0,??(1)=10??(??)????(??).(2.28)
Let ??=??-??, then we have ??(0)=????(0)=0,??(1)=10??(??)????(??),(2.29) and ??????-????(??)=??????-????(??)-??????-????(??), which implies that -??????-????(??)=-??????-????(??)+??-(??).(2.30) Since ??(??)=??+(??)-??-(??), then (2.27) is transformed to (2.7), that is, ??-?? is a positive solution of the BVP (2.27). By (2.7), ????[??(??)-??(??)] is a positive solution of the BVP (1.1).

It is well known that the BVP (2.27) is equivalent to the fixed points of the mapping ?? given by ?(????)(??)=10??????(??,??)??,????[]??(??)-??(??)*[],-??(??)-??(??)*?+??+?(??)????.(2.31)

The basic space used in this paper is ??=??([0,1];R), where R is a real number set. Obviously, the space ?? is a Banach space if it is endowed with the norm as follows: ????=max[]???0,1||||,??(??)(2.32) for any ?????. Let ???=?????:??(??)=(1-??)??(??)4?????,(2.33) then ?? is a cone of ??.

For the convenience in presentation, we now present some assumptions to be used in the rest of the paper.

(????) For any fixed ???(0,1) and for any ???(0,1], there exist constants ????>1,??=1,2,3,4, such that, for all (??,??,??)?(0,1)×[0,+8)×(-8,+0], ????1??(??,??,??)=??(??,????,??)=????3????(??,??,??),??2??(??,??,??)=??(??,??,????)=????4??(??,??,??).(2.34)

(????)?10??-(??)????>0, for any ???(0,1),??(??,1,-1)>0, and ?10?F(??)??(??,1,-1)+??+???(??)????=2(1+??)??1+??2,(2.35) where ??=8??,1-??(2.36) and ??,??,F are defined by (2.13), (2.19), and (2.23), respectively.

Remark 2.11. If ??>1, the reversed inequality of (2.34) holds, that is, for all (??,??,??)?(0,1)×[0,+8)×(-8,+0], one has ????3??(??,??,??)=??(??,????,??)=????1??(??,??,??),????4??(??,??,??)=??(??,??,????)=????2??(??,??,??).(2.37)

Lemma 2.12. Suppose (????)-(????) holds, for any fixed ???(0,1), ??(??,??,??) is nondecreasing in ?? on [0,+8) and nonincreasing in ?? on (-8,0]; and for any [??,??]?(0,1), lim(??,??)?(+8,-8)min[]?????,????(??,??,??)|??||??|=+8.(2.38)

Proof. Let 0=??1=??2. If ??2=0, obviously ??(??,??1,??)=??(??,??2,??) holds. If ??2?0, let ??0=??1/??2, then 0<??0=1. It follows from (????) that ?????,??1??,??=????,??0??2?,??=????30?????,??2??,??=????,??2?.,??(2.39) Thus ??(??,??,??) is nondecreasing in ?? on [0,8).
On the other hand, for any ??1=??2=0, if ??1=0, obviously ??(??,??,??1)=??(??,??,??2) holds. If ??1?0, let ??1=??2/??1, then 0<??1=1. By (????), we have ?????,??,??2??=????,??,??1??1?=????41?????,??,??1??=????,??,??1?.(2.40) Thus ??(??,??,??) is nonincreasing in ?? on (-8,0].
Now choose ??>1 and ??<-1, then by Remark 2.11, we have ??(??,??,??)=??(??,|??|,-|??|)=|??|??3??(??,1,-|??|)=|??|??3|??|??4??(??,1,-1).(2.41) Thus for any [??,??]?(0,1), and any ???[??,??], we have min[]?????,????(??,??,??)|??||??|=|??|??3-1|??|??4-1min[]?????,????(??,1,-1)>0.(2.42) Therefore lim(??,??)?(+8,-8)min[]?????,????(??,??,??)|??||??|=+8.(2.43)

Lemma 2.13. Assume that (????)(????) holds. Then ??:????? is well defined. Furthermore, ??:????? is a completely continuous operator.

Proof. For any fixed ?????, there exists a constant ??>0 such that ????=??. And then, []0=??(??)-??(??)*=??(??)=????=??,0=????[]??(??)-??(??)*=???0(??-??)??-1[]??(??)-??(??)*??G(??)????=.G(??)(2.44) By (2.44) and (????)-(????), we have ?(????)(??)=10??????(??,??)??,????[]??(??)-??(??)*[],-??(??)-??(??)*?+??+?=?(??)????10????????(??,??)??,?G(??)+1,-(??+1)+??+?=???(??)?????G(??)+1??1(??+1)??2?10?F(??)??(??,1,-1)+??+?(??)????<+8,(2.45) which implies that the operator ??:????? is well defined.
Next let ??*(??)=??(??,????[??(??)-??(??)]*,-[??(??)-??(??)]*)+??+(??), for any ?????, by (2.11), we have ???????=10??(1-??)??-??-1??G(??-??-1)*??(??)????+??-??-1?1-??10????(??)??*=1(??)?????1-??10??(1-??)??-??-1??G(??-??-1)*??(??)????+??-??-1?1-??10????(??)??*(??)????.(2.46) On the other hand, by (2.11), (2.22), and (2.46), we also have (????)(??)=????-??-1?(1-??)10??(1-??)??-??-1??G(??-??)*??(??)????+??-??-1?1-??10????(??)??*=1(??)????2????-??-1?(1-??)10??(1-??)??-??-1??G(??-??)*??(??)????+??-??-1?1-??10????(??)??*=1(??)????2?????-??-1?(1-??)10??(1-??)??-??-1??G(??-??)*??(??)????+??-??-1?1-??10????(??)??*?+??(??)??????-??-1?2(1-??)10????(??)??*=1(??)????2?????-??-1(?1-??)10??(1-??)??-??-1??G(??-??)*(??)????+B????-??-1?1-??10??(1-??)??-??-1??G(??-??)*(?+????)??????-??-1?2(1-??)10????(??)??*=(??)????(1-??)??(??)2×1(??-??-1)?1-??10??(1-??)??-??-1G??(??-??-1)*??(??)????+??-??-12?(1-??)10????(??)??*=(??)????(1-??)??(??)?12(??-??-1)?1-??10??(1-??)??-??-1??G(??-??-1)*??(??)????+??-??-1?1-??10????(??)??*?=(??)????(1-??)??(??)4?1?1-??10??(1-??)??-??-1??G(??-??-1)*??(??)????+??-??-1?1-??10????(??)??*?.(??)????(2.47) So we have (????)(??)=(1-??)??(??)4[].??????,???0,1(2.48) which yields that ??(??)???.
At the end, using standard arguments, according to the Ascoli-Arzela Theorem, one can show that ??:????? is continuous. Thus ??:????? is a completely continuous operator.

Lemma 2.14 (see [27]). Let ?? be a real Banach space, ????? be a cone. Assume O1,O2 are two bounded open subsets of ?? with ???O1, O1?O2, and let ??:??n(O2?O1)??? be a completely continuous operator such that either(1)??????=????, ?????n??O1, and ??????=????, ?????n??O2 or (2)??????=????, ?????n??O1, and ??????=????, ?????n??O2.
Then ?? has a fixed point in ???(O2?O1).

3. Main Results

Theorem 3.1. Suppose (????)(????) hold. Then the BVP (1.1) has at least one positive solutions ??(??), and ??(??) satisfies ??(??)=??G(??)G(??-??)??BG(??)1+???1-????-1-??-?????????[],???0,1.(3.1)

Proof. Let ??=8??,1-??(3.2) and O1={?????:????<??}. Then, for any ?????O1,???[0,1], notice that 0<??<1, we have []??(??)-??(??)*??=??(??)=????=??,??[]??(??)-??(??)*=1G?(??)??0(??-??)??-1[]??(??)-??(??)*??????=G(??)=??.(3.3) By (????) and (????), one has ??????=max0=??=1?10??????(??,??)??,????[]??(??)-??(??)*[],-??(??)-??(??)*?+??+?=?(??)????10F??????(??)??,?G(??),-??+??+?=?(??)????10?F(??)??(??,??+1,-(??+1))+??+?(??)????=2(??+1)??1+??2?10?F(??)??(??,1,-1)+??+?(??)????=??=????.(3.4) Therefore, ??????=????,?????n??O1.(3.5) On the other hand, choose a real number ??>0 such that ????-??(1-??)??-??×G(??-??)??(1-??)2???????64G(??+1)??-??-1B(1-??)+??1-????-??-1?2=1.(3.6)
From (2.38), there exists ??>?? such that, for any ???[??,??], ??(??,??,??)=??|??||??|,|??|,|??|=??.(3.7) Take ???????>max1,8?????(1-??)??-??-1B(1-??)+??1-????-??-1?,8??G(??+1)(1-??)?????????-??-1B(1-??)+??1-????-??-1??????,,??(3.8) then ??>??. Let O2={?????:????<??}, for any ?????n??O2 and for any ???[??,??], we have 4??(??)-??(??)=??(??)-????(??)=??(??)-1(1-??)????(??)=2=??(??)(1-??)??(??)8(??=1-??)8?????-??-1(B1-??)+??1-????-??-1?????=??>0,??[]=??(??)-??(??)(1-??)8?????-??-1(B1-??)+??1-????-??-1??????0(??-??)??-1=G(??)????(1-??)???????8G(??+1)??-??-1B(1-??)+??1-????-??-1???=??>0.(3.9) So for any ?????n??O2,???[??,??], by (3.7)-(3.9), we have ???????=10??????(??,??)??,????[]??(??)-??(??)*[],-??(??)-??(??)*?+??+?=?(??)????10???(??,??)????,????[??](??)-??(??)*[??],-(??)-??(??)*?=??????????-??-1(1-??)?G(??-??)??????(1-??)??-??-1?????,????[??](??)-??(??)*[??],-(??)-??(??)*?=????????-??-1(1-??)?G(??-??)??????(1-??)??-??-1??||????[]||||||=????(??)-??(??)??(??)-??(??)??????-??(1-??)??-??×??G(??-??)(1-??)2???????64G(??+1)??-??-1B(1-??)+??1-????-??-1?2??2=??2=??=????.(3.10) Thus ??????=????,?????n??O2.(3.11) By Lemma 2.14, ?? has at least one fixed points ?? such that ??=||??||=??.
In the end, ??(??)-??(??)=(1-??)??(??)4????-????(??)=(1-??)??(??)4??-????(??)=????(??)>0,???(0,1).(3.12) Equation (3.12) implies that ??(??)>??(??),???(0,1), and ??(??)=????(??(??)-??(??))=?????????(??)=????0(??-??)??-1=??(??)??????G(??)G(??-??)??BG(??)1+???1-????-1-??-?????????>0,???(0,1).(3.13) So by Lemma 2.10, the BVP (1.1) has at least one positive solution ??, and ?? satisfies (3.13).

Example 3.2. Consider the following semipositone boundary value problem with fractional order ??=21/8: -??21/8????(??)=2???(??)+1/8???(??)31012×4v-1(1-??)v????,0<??<1,1/8??(0)=??9/8??(0)=0,??1/8???(1)=10??1/8??(??)????(??).(3.14) Let ????????1??(??)=0,???0,4?,?11,???4,34?,?32,???4?,,1(3.15) then the BVP (3.14) becomes the 4-Point BVP with coefficients -??21/8????(??)=2???(??)+1/8???(??)31012×4v-1(1-??)v????,0<??<1,1/8??(0)=??9/8??(0)=0,??1/8??(1)=??1/8???14?+??1/8???34?.(3.16) Then the BVP (3.14) has at least one positive solution ??(??), and ??(??) satisfies ???(??)=62.42891.1357??5/2-10??2121/8?[].,???0,1(3.17)

Proof. Obviously, ??=21/8,??=1/8, and ?0=??=10??3/2?1????(??)=4?3/2+?34?3/2?˜0.7745<1,B=10??3/23(1-??)????(??)=4?14?3/2+14?34?3/2˜0.2561.(3.18) On the other hand, we have ?????????(??,??)=1[](??,??)=??(1-??)3/2G??(5/2),0=??=??=1,2[??](??,??)=(1-??)3/2-(??-??)3/2??G(5/2),0=??=??=1,???????????(??)=2?14?,??+2??2?34?1,??,0=??<4,??1?14?,??+2??2?34?,1,??43=??<4,??1?14?,??+2??1?34?,3,??4=??=1.(3.19) Thus ????(??)=0, 0=??<1 and ??(??) is increasing. So (????) holds.
Take ????(??,??,??)=2-??31012×4v([]11-??),(??,??,??)?(0,1)×0,+8)×(-8,0,??(??)=-v??,(3.20) then ?????,??(??),-????1/8???=??(??)2???(??)+1/8???(??)31012×4v(1-??),??+(??)=0,??-1(??)=v??.(3.21) For any 0<??=1, we also have ??3??(??,??,??)=??(??,????,??)=??3/2????(??,??,??),2??(??,??,??)=??(??,????,??)=??4??(??,??,??).(3.22) Then (????) holds.
In the end, we notice ??F(??)=??(??)+1-????(1-??)??-??-1G=??(??-??-1)G+1(??-??-1)(1-??)G=1(??-??-1)G(??-??-1)(1-??)˜5.0041,(3.23) then ?10?F(??)??(??,1,-1)+??+??(??)????=5.00411021012×4v(1-??)????˜1.33443×1011.(3.24) Moreover, 1??=???G(??-??-1)1+B??10??-(??)????˜9.0817,??=8????1-??=322.1889,2(1+??)??1+??2=322.18892(1+322.1889)5˜4.5688×10-11,(3.25) which implies that (????) holds.
According to Theorem 3.1, the BVP (3.14) has at least one positive solution ??(??), and ??(??) satisfies ???(??)=62.42891.1357??5/2-10??2121/8?[].,???0,1(3.26)

Acknowledgments

The paper is Funded by the Humanities and Social Sciences Foundation of Ministry of Education (no. 09YJC630056), the National Natural Science Foundation of China (no. 51009060, no. 50909042), the Fundamental Research Funds for the Central Universities (no. 2009B13414), and the Priority Academic Program Development of Jiangsu Higher Education Institutions(Coastal Development Conservancy).