Boundary value problems for fractional differential equations with ...

3 downloads 0 Views 176KB Size Report
Jun 28, 2012 - There has been a significant development in the theory of initial and .... Definition 2.2 The Riemann-Liouville fractional integral of order q for a ...
Sudsutad and Tariboon Advances in Difference Equations 2012, 2012:93 http://www.advancesindifferenceequations.com/content/2012/1/93

RESEARCH

Open Access

Boundary value problems for fractional differential equations with three-point fractional integral boundary conditions Weerawat Sudsutad1,2 and Jessada Tariboon1,2* *

Correspondence: [email protected] 1 Department of Mathematics, Faculty of Applied Science, King Mongkut’s University of Technology North Bangkok, Bangkok, 10800, Thailand 2 Centre of Excellence in Mathematics, CHE, Sri Ayutthaya Road, Bangkok, 10400, Thailand

Abstract This article studies a boundary value problem of nonlinear fractional differential equations with three-point fractional integral boundary conditions. Some new existence results are obtained by applying standard fixed point theorems. As an application, we give two examples that illustrate our results. MSC: 26A33; 34B15 Keywords: existence; Caputo fractional derivative; Riemann-Liouville fractional integral; boundary value problem

1 Introduction Fractional differential equations have recently proved to be valuable tools in the modelling of many phenomena in various field of science and applications, such as physics, mechanics, chemistry, biology, economics, control theory, aerodynamics, engineering, etc. See [–]. There has been a significant development in the theory of initial and boundary value problems for nonlinear fractional differential equations; see, for example, [–]. Ahmad and co-authors have studied the existence and uniqueness of solutions of nonlinear fractional differential and integro-differential equations for a variety of boundary conditions using standard fixed-point theorems and Leray-Schauder degree theory. Ahmad et al. [] discusses the existence and uniqueness of solutions of fractional integrodifferential equations for fractional nonlocal integral boundary conditions. Ahmad et al. [] and references therein give details of recent work on the properties of solutions of sequential fractional differential equations. Ahmad et al. [] considers solutions of fractional differential equations with non-separated type integral boundary conditions. In Ahmad et al. [], the Krasnoselskii fixed point theorem and the contraction mapping principle are used to prove the existence of solutions of the nonlinear Langevin equation with two fractional orders for a number of different intervals. Ahmad et al. [] discusses the existence and uniqueness of solutions of nonlinear fractional differential equations with three-point integral boundary conditions. Cabada et al. [] have also studied properties of solutions of nonlinear fractional differential equations. They used the properties of the associated Green’s function and the Guo-Krasnosellskii fixed-point theorem to investigate the existence of positive solutions of nonlinear fractional differential equations with integral boundary-value conditions. © 2012 Sudsutad and Tariboon; licensee Springer. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Sudsutad and Tariboon Advances in Difference Equations 2012, 2012:93 http://www.advancesindifferenceequations.com/content/2012/1/93

Page 2 of 10

Motivated by the papers [] and [], this article is concerned with the existence and uniqueness of solutions for a boundary value problem of nonlinear fractional differential equations with three-point fractional integral boundary conditions given by   Dq x(t) = f t, x(t) , t ∈ [, ], q ∈ (, ],   x() = , α I p x (η) = x(),  < η < ,

c

(.) (.)

where c Dq denotes the Caputo fractional derivative of order q, I p is the Riemann-Liouville fractional integral of order p, f : [, ] × R → R is a continuous function and α ∈ R is such that α = (p + )/ηp+ . By C([, ], R) we denote the Banach space of all continuous functions from [, ] into R with the norm    x = sup x(t); t ∈ [, ] . We note that if p = , then condition (.) reduces to the usual three-point integral condition. In such a case, the boundary condition corresponds to the area under the curve of solutions x(t) from t =  to t = η.

2 Preliminaries In this section, we introduce notations, definitions of fractional calculus and prove a lemma before stating our main results. Definition . For a continuous function f : [, ∞) → R, the Caputo derivative of fractional order q is defined as c

Dq f (t) =

 (n – q)



t

(t – s)n–q– f (n) (s) ds,

n –  < q < n, n = [q] + ,



provided that f (n) (t) exists, where [q] denotes the integer part of the real number q. Definition . The Riemann-Liouville fractional integral of order q for a continuous function f (t) is defined as I q f (t) =

 (q)



t

(t – s)q– f (s) ds,

q > ,



provided that such integral exists. Definition . The Riemann-Liouville fractional derivative of order q for a continuous function f (t) is defined by Dq f (t) =

n t d  (t – s)n–q– f (s) ds, (n – q) dt 

n = [q] + ,

provided that the right-hand side is pointwise defined on (, ∞). Furthermore, we note that the Riemann-Liouville fractional derivative of a constant is usually nonzero which can cause serious problems in real would applications. Actually,

Sudsutad and Tariboon Advances in Difference Equations 2012, 2012:93 http://www.advancesindifferenceequations.com/content/2012/1/93

Page 3 of 10

the relationship between the two-types of fractional derivative is as follows c

Dq f (t) =

 (n – q)



t 

f (n) (s) ds (t – s)q+–n

n–

f (k) () t k–q (k – q + ) k=

 n– (k) f () k q t , t > , n –  < q < n. = D f (t) – k! = Dq f (t) –

k=

So, we prefer to use Caputo’s definition which gives better results than those of RiemannLiouville. Lemma . [] Let q > , then the fractional differential equation c

Dq u(t) = 

has solution u(t) = k + k t + k t  + · · · + kn– t n– ,

ki ∈ R, i = , , , . . . , n – ,

where n is the smallest integer greater than or equal to q. Lemma . [] Let q > , then I q c Dq u(t) = u(t) + k + k t + k t  + · · · + kn– t n– , for some ki ∈ R, i = , , , . . . , n –  where n is the smallest integer greater than or equal to q. Lemma . Let α = c

(p+) , ηp+

 < q ≤ . Then for h ∈ C([, ], R), the problem

Dq x(t) = h(t),

 < t < ,   α I p x (η) = x(),

x() = ,

(.) (.)

has a unique solution  (p + )t (t – s) h(s) ds – ( – s)q– h(s) ds (q)((p + ) – αηp+ )   η s αp(p + )t (η – s)p– (s – r)q– h(r) dr ds. + (q)((p + ) – αηp+ )  

 x(t) = (q)



t

q–

(.)

Proof By applying Lemma ., we may reduce (.) to an equivalent integral equation x(t) =

 (q)



t

(t – s)q– h(s) ds – c – c t, 

for some c , c ∈ R.

(.)

Sudsutad and Tariboon Advances in Difference Equations 2012, 2012:93 http://www.advancesindifferenceequations.com/content/2012/1/93

Page 4 of 10

From x() = , it follows c = . Using the Riemann-Liouville integral of order p for (.), we have I p x(t) = =

 (p)





t

(t – s)p– 

 (p)(q)

t

 (q)



s

 (s – r)q– h(r) dr – c s ds



s

(t – s)p– (s – r)q– h(r) dr ds – c 



t p+ . (p + )

The second condition of (.) implies that α (p)(q) =

 (q)



η s



(η – s)p– (s – r)q– h(r) dr ds – c



αηp+ (p + )



( – s)q– h(s) ds – c . 

Thus, (p + ) c = (q)((p + ) – αηp+ ) –





( – s)q– h(s) ds 

αp(p + ) (q)((p + ) – αηp+ )



η s



(η – s)p– (s – r)q– h(r) dr ds.



Substituting the values of c and c in (.), we obtain the solution (.).



In the following, for the sake of convenience, set =

 (p + ) αηp+q (p + ) + + . p+ (q + ) (q + )|(p + ) – αη | (p + q + )|(p + ) – αηp+ |

(.)

3 Main results Now we are in the position to establish the main results. Theorem . Assume that there exists a constant L >  such that (H ) |f (t, x) – f (t, y)| ≤ L|x – y|, for each t ∈ [, ], and all x, y ∈ R. If L < , where  is defined by (.), then the BVP (.)-(.) has a unique solution on [, ]. Proof Transform the BVP (.)-(.) into a fixed point problem. In view of Lemma ., we consider the operator F : C([, ], R) → C([, ], R) defined by F(x)(t) =

 (q) – +



t

  (t – s)q– f s, x(s) ds



(p + )t (q)((p + ) – αηp+ ) αp(p + )t (q)((p + ) – αηp+ )









  ( – s)q– f s, x(s) ds

η s

  (η – s)p– (s – r)q– f r, x(r) dr ds.



Obviously, the fixed points of the operator F are solution of the problem (.)-(.). We shall use the Banach fixed point theorem to prove that F has a fixed point. We will show that F is a contraction.

Sudsutad and Tariboon Advances in Difference Equations 2012, 2012:93 http://www.advancesindifferenceequations.com/content/2012/1/93

Page 5 of 10

Let x, y ∈ C([, ], R). Then, for each t ∈ [, ] we have   F(x)(t) – F(y)(t) t       (t – s)q– f s, x(s) – f s, y(s)  ds ≤ (q)        (p + ) q–  f s, x(s) – f s, y(s)  ds ( – s) + (q)|(p + ) – αηp+ |  η s      αp(p + ) + (η – s)p– (s – r)q– f r, x(r) – f r, y(r)  dr ds p+ (q)|(p + ) – αη |   Lx – y t ≤ (t – s)q– ds (q)   L(p + )x – y ( – s)q– ds + (q)|(p + ) – αηp+ |  η s αp(p + )Lx – y (η – s)p– (s – r)q– dr ds + (q)|(p + ) – αηp+ |   =

L(p + )x – y Lx – y + (q + ) (q + )|(p + ) – αηp+ | +

αp(p + )Lx – yηp+q B(q + , p) , (q + )|(p + ) – αηp+ |

where B is the beta function.

By using the property of beta function, B(q + , p) =

(q+)(p) , (p+q+)

we have

  L(p + )x – y F(x)(t) – F(y)(t) ≤ Lx – y + (q + ) (q + )|(p + ) – αηp+ | +

αηp+q (p + )Lx – y . (p + q + )|(p + ) – αηp+ |

Thus   F(x) – F(y) ≤ Lx – y. Therefore, F is a contraction. Hence, by Banach fixed point theorem, we get that F has a fixed point which is a solution of the problem (.)-(.).  The following result is based on Schaefer’s fixed point theorem. Theorem . Assume that: (H ) The function f : [, ] × R → R is continuous. (H ) There exists a constant M >  such that |f (t, x)| ≤ M for each t ∈ J and all x ∈ R. Then the BVP (.)-(.) has at least one solution on [, ]. Proof We shall use Schaefer’s fixed point theorem to prove that F has a fixed point. We divide the proof into four steps.

Sudsutad and Tariboon Advances in Difference Equations 2012, 2012:93 http://www.advancesindifferenceequations.com/content/2012/1/93

Page 6 of 10

Step I. Continuity of F. Let {xn } be a sequence such that xn → x in C([, ], R). Then for each t ∈ [, ]   F(xn )(t) – F(x)(t) t       ≤ (t – s)q– f s, xn (s) – f s, x(s)  ds (q)        (p + ) ( – s)q– f s, xn (s) – f s, x(s)  ds + p+ (q)|(p + ) – αη |  η s      αp(p + ) + (η – s)p– (s – r)q– f r, xn (r) – f r, x(r)  dr ds p+ (q)|(p + ) – αη |   t       ≤ (t – s)q– sup f s, xn (s) – f s, x(s)  ds (q)  s∈[,]       (p + ) + ( – s)q– sup f s, xn (s) – f s, x(s)  ds p+ (q)|(p + ) – αη |  s∈[,] η s αp(p + ) + (η – s)p– (s – r)q– (q)|(p + ) – αηp+ |        × sup f r, xn (r) – f r, x(r)  dr ds. r∈[,]

Since f is continuous function, then F(xn ) – F(x) →  as n → ∞. This means that F is continuous. Step II. F maps bounded sets into bounded sets in C([, ], R). So, let us prove that for any r > , there exists a positive constant l such that for each x ∈ Br = {x ∈ C([, ], R) : x ≤ r}, we have F(x) ≤ l. Indeed, we have for any x ∈ Br   F(x)(t) ≤

 (q)



t

   (t – s)q– f s, x(s)  ds



(p + ) + (q)|(p + ) – αηp+ | +

αp(p + ) (q)|(p + ) – αηp+ |







   ( – s)q– f s, x(s)  ds

η s



   (η – s)p– (s – r)q– f r, x(r)  dr ds



which in view of (H ) gives   F(x)(t) ≤ M (q) + +



t

(t – s)q– ds 

(p + )M (q)|(p + ) – αηp+ | αp(p + )M (q)|(p + ) – αηp+ |





( – s)q– ds 





η s

(η – s)p– (s – r)q– dr ds.



Hence, we deduce that   F(x)(t) ≤

(p + )M M + (q + ) (q + )|(p + ) – αηp+ | +

α(p + )Mηp+q . (p + q + )|(p + ) – αηp+ |

Sudsutad and Tariboon Advances in Difference Equations 2012, 2012:93 http://www.advancesindifferenceequations.com/content/2012/1/93

Thus,   F(x) ≤ M := l. Step III. We prove that F(Br ) is equicontinuous with Br defined as in Step II. Let  ≤ t < t ≤  and x ∈ Br , then   F(x)(t ) – F(x)(t )  t       (t – s)q– – (t – s)q– f s, x(s) ds =  (q)  t    + (t – s)q– f s, x(s) ds (q) t    (t – t )(p + ) ( – s)q– f s, x(s) ds – p+ (q)((p + ) – αη )   η s    αp(p + )(t – t ) p– q– + (η – s) (s – r) f r, x(r) dr ds (q)((p + ) – αηp+ )   t      (t – s)q– – (t – s)q– f s, x(s)  ds ≤ (q)  t     (t – s)q– f s, x(s)  ds + (q) t     (t – t )(p + ) ( – s)q– f s, x(s)  ds + p+ (q)|(p + ) – αη |  η s    αp(p + )(t – t ) (η – s)p– (s – r)q– f r, x(r)  dr ds + (q)|(p + ) – αηp+ |   t   M (t – s)q– – (t – s)q– ds ≤ (q)  t M + (t – s)q– ds (q) t  (t – t )(p + )M ( – s)q– ds + (q)|(p + ) – αηp+ |  η s αp(p + )(t – t )M (η – s)p– (s – r)q– dr ds + (q)|(p + ) – αηp+ |   ≤

 q q  M(t – t )q M  (t – t )q + t – t + (q + ) (q + ) +

(t – t )(p + )M (q + )|(p + ) – αηp+ |

+

αηp+q (t – t )(p + )M . (p + q + )|(p + ) – αηp+ |

Actually, as t → t , the right-hand side of the above inequality tends to zero. As a consequence of Steps I to III together with the Arzela-Ascoli theorem, we get that F : C([, ], R) → C([, ], R) is completely continuous. Step IV. A priori bounds.

Page 7 of 10

Sudsutad and Tariboon Advances in Difference Equations 2012, 2012:93 http://www.advancesindifferenceequations.com/content/2012/1/93

Page 8 of 10

We show that the set     E = x ∈ C [, ], R : x = λF(x) for some  < λ <  is bounded. Let x ∈ E. Then x(t) = λF(x)(t) for some  < λ < . Thus, for each t ∈ [, ] we have λ (q)

F(x)(t) =



t

  (t – s)q– f s, x(s) ds



λ(p + )t – (q)((p + ) – αηp+ ) αλp(p + )t + (q)((p + ) – αηp+ )







  ( – s)q– f s, x(s) ds

η s



  (η – s)p– (s – r)q– f r, x(r) dr ds.



This implies by (H ) that for all t ∈ [, ], we get   F(x)(t) ≤ λM (q)



t

(t – s)q– ds 

λ(p + )M + (q)|(p + ) – αηp+ | + ≤

αλp(p + )M (q)|(p + ) – αηp+ |





( – s)q– ds 





η s

(η – s)p– (s – r)q– dr ds



(p + )M M + (q + ) (q + )|(p + ) – αηp+ | +

α(p + )Mηp+q . (p + q + )|(p + ) – αηp+ |

Hence, we deduce that   F(x) ≤ M := R. This shows that the set E is bounded. As a consequence of Schaefer’s fixed point theorem, we conclude that F has a fixed point which is a solution of the problem (.)-(.). 

4 Examples In this section, in order to illustrate our results, we consider two examples. Example . Consider the following three-point fractional integral boundary value problem 

e– sin t |x| , t ∈ [, ], · ( + t) |x| + 

√      = x(). x() = , I x  c



D  x(t) =

(.) (.)

√  Set η = /, q = /, p = /, α =  = (p + )/ηp+ = (/)/(/)  and f (t, x) =  (e– sin t /( + t) )(|x|/( + |x|)). Since f (t, x) – f (t, y) ≤ (/)x – y, then, (H ) is satisfied

Sudsutad and Tariboon Advances in Difference Equations 2012, 2012:93 http://www.advancesindifferenceequations.com/content/2012/1/93

Page 9 of 10

with L = /. We can show that 

(p + ) αηp+q (p + )  + + p+ (q + ) (q + )|(p + ) – αη | (p + q + )|(p + ) – αηp+ | √    ,  π  = + ≈ . < . √ + √ √   π , π –  , π – 



L = L

Hence, by Theorem ., the boundary value problem (.)-(.) has a unique solution on [, ]. Example . Consider the following three-point fractional integral boundary value problem 

t  e–t ,  ≤ t ≤ ,  + tx

    x() = , Ix = x().  

c



D  x(t) =

(.) (.) 

Set η = /, q = /, p = /, α = / = (p + )/ηp+ = (/)/(/)  and f (t, x) =   t  e–t /( + tx ). Clearly |f (t, x)| ≤ |t  e–t /( + tx )| ≤  = M. Hence, all the conditions of Theorem . are satisfied and consequently the problem (.)-(.) has at least one solution.

Competing interests The authors declare that they have no competing interests. Authors’ contributions All authors contributed equally and significantly in writing this article. All authors read and approved the final manuscript. Acknowledgements The authors thank the referees for several useful remarks and interesting comments. This research is supported by the Centre of Excellence in Mathematics, the Commission on Higher Education, Thailand. Received: 9 February 2012 Accepted: 8 June 2012 Published: 28 June 2012 References 1. Podlubny, I: Fractional Differential Equations. Academic Press, San Diego (1999) 2. Hilfer, R: Fractional Calculus in Physics. World Scientific, Singapore (2000) 3. Kilbas, AA, Srivastava, HM, Trujillo, JJ: Theory and Applications of Fractional Differential Equations. Elsevier, Amsterdam (2006) 4. Mainardi, F: Fractional Calculus and Waves in Linear Viscoelasticity: An Introduction to Mathematical Models. Imperial College Press, Singapore (2010) 5. Ortigueira, MD: Fractional Calculus for Scientists and Engineers. Heidelberg, Springer (2011) 6. Petras, I: Fractional-Order Nonlinear Systems: Modeling, Analysis and Simulation. Springer, Berlin (2011) 7. Ahmad, B: Existence of solutions for irregular boundary value problems of nonlinear fractional differential equations. Appl. Math. Lett. 23, 390-394 (2010) 8. Ahmad, B, Sivasundaram, S: On four-point nonlocal boundary value problems of nonlinear integro-differential equations of fractional order. Appl. Math. Comput. 217, 480-487 (2010) 9. Bai, Z, Qiu, T: Existence of positive solution for singular fractional differential equation. Appl. Math. Comput. 215, 2761-2767 (2009) 10. Delbosco, D, Rodino, L: Existence and uniqueness for a nonlinear fractional differential equation. J. Math. Anal. Appl. 204, 609-625 (1996) 11. Deng, J, Ma, L: Existence and uniqueness of solutions of initial value problems for nonlinear fractional differential equations. Appl. Math. Lett. 23, 676-680 (2010) 12. Han, X, Wang, T: The existence of solutions for a nonlinear fractional multi-point boundary value problem at resonance. Int. J. Differ. Equ. 2011, Article ID 401803 (2011) 13. Jiang, W: The existence of solutions to boundary value problems of fractional differential equations at resonance. Nonlinear Anal. 74, 1987-1994 (2011)

Sudsutad and Tariboon Advances in Difference Equations 2012, 2012:93 http://www.advancesindifferenceequations.com/content/2012/1/93

14. Rehman, MU, Khan, RA: Existence and uniqueness of solutions for multi-point boundary value problems for fractional differential equations. Appl. Math. Lett. 23, 1038-1044 (2010) 15. Wang, G, Liu, W, Yang, J, Zhu, S, Zheng, T: The existence of solutions for a fractional 2m-point boundary value problems. J. Appl. Math. 2012, Article ID 841349 (2012) 16. Ahmad, B, Nieto, JJ: Riemann-Liouville fractional integro-differential equations with fractional nonlocal integral conditions. Bound. Value Probl. 2011, 36 (2011) 17. Ahmad, B, Nieto, JJ: Sequential fractional differential equations with three-point boundary conditions. Comput. Math. Appl. (in press). http://dx.doi.org/10.1016/j.camwa.2012.02.036 18. Ahmad, B, Nieto, JJ, Alsaedi, A: Existence and uniqueness of solutions for nonlinear fractional differential equations with non-separated type integral boundary conditions. Acta Math. Sci., Ser. B 31(6), 2122-2130 (2011) 19. Ahmad, B, Nieto, JJ, Alsaedi, A, El-Shahed, M: A study of nonlinear Langevin equation involving two fractional orders in different intervals. Nonlinear Anal., Real World Appl. 13, 599-606 (2012) 20. Ahmad, B, Ntouyas, SK, Alsaedi, A: New existence results for nonlinear fractional differential equations with three-point integral boundary conditions. Adv. Differ. Equ. 2011, Article ID 107384 (2011) 21. Cabada, A, Wang, G: Positive solutions of nonlinear fractional differential equations with integral boundary value conditions. J. Math. Anal. Appl. 389, 403-411 (2012)

doi:10.1186/1687-1847-2012-93 Cite this article as: Sudsutad and Tariboon: Boundary value problems for fractional differential equations with three-point fractional integral boundary conditions. Advances in Difference Equations 2012 2012:93.

Page 10 of 10