Existence result for a semipositone fractional boundary value problem

Author(s): Serife Müge Ege1, Fatma Serap Topal1
1Department of Mathematics, Ege University, Bornova, Izmir 35100 Turkey.
Copyright © Serife Müge Ege, Fatma Serap Topal. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

This work deals with a boundary value problem for a nonlinear semipositone multi-point fractional differential equation. By using the Schauder fixed point theorem, we show the existence of one solution for this problem. Our result extend some recent works in the literature.

Keywords: Caputo fractional derivative; Boundary value problem; Changing sign nonlinearity; Fixed point theorem.

1. Introduction

Fractional calculus is the extension of integer derivatives to real or complex number derivatives. Therefore, it has gained an important place in terms of examining the sensitivity of many phenomena in natural sciences. Most of the problems in scientific fields such as physics, biology, astronomy, medical sciences, optical fibers, chemical biology, radiology can be modeled by differential equations [1,2,3,4,5,6]. The fractional orders of the differential equations enable the investigations to be more comprehensive and obtain optimal solutions. That’s why fractional differential equations have been the focus of various researchers over the past decades. Articles and books have been written to demonstrate the existence of solutions to such equations or to develop new solution methods [7,8,9].

Many people have studied the existence and positive solutions or multiplicity of solutions for nonlinear fractional boundary value problems using fixed point theorems such as the Schauder fixed-point theorem, the Guo-Krasnosel’skii fixed-point theorem, and the Leggett-Williams fixed-point theorem. In recent studies, boundary value problems involving multi point boundary condition and fractional differential equations have attracted attention.

In order to guarantee the existence of positive solutions for boundary value problems, the nonlinearity is usually nonnegative. If the nonlinearity changing sign, it will bring much more difficulties to the study of the problem. Because of this, there is only a few study on such problems [10,11,12,13,14].

In [15], Ahmad consider the two point Liouville-Caputo boundary value problem

\begin{eqnarray*}\label{1} \left. \begin{array}{ll} ^{c}D^\beta x(t)=-f(t,x(t)), \quad t \in (a,b)\\ x(a)=\delta_1, \quad x(b)=\delta_2, \quad \delta_i \in \mathbb{R}, i=1,2 \end{array} \right. \end{eqnarray*} where \(^{c}D^\beta\) denotes the Caputo fractional derivative of order \(\beta\) with \(1< \beta< 2\) and \(f\) is a continuous function.

In this paper, we investigate the existence result for the following m-point fractional boundary value problem with changing sign nonlinearity,

\begin{eqnarray}\label{11} D_{a}^\beta u(t)+f(t,u(t))=0, \quad t \in (a,b) \end{eqnarray}
(1)
\begin{eqnarray}\label{12} u(a)=\delta_1, \quad u(b)=\sum_{i=1}^{m-2}\xi_iu(\eta_i)+\delta_2, \quad \delta_i \in \mathbb{R}, i=1,2, \end{eqnarray}
(2)
where \(D_a^\beta\) is the standard Caputo fractional derivatives of fractional order \(\beta\) with \(1< \beta\leq 2\) , \(\xi_i \quad (1\leq i\leq m-2)\) are positive real constants with \( 0< \sum_{i=1}^{m-2} \xi_i< 1\), \(\eta_i \in(a,b)\) and a continuous function \(f:[a,b]\times \mathbb{R}\rightarrow \mathbb{R}\) may change sign.

The paper is organized as follows: In Section 2, present some background materials and preliminaries. In Section 3, we give an existence result and also an example is given to exemplify main result.

2. Preliminaries

In this section we gather some preliminary definitions, theorems and fundamental results of fractional calculus theory that will be used in subsequent section. Basic properties of the fractional calculus and its applications can be found in Kilbas, Srivastava and Trujillo [4] and Podlubny [6].

Definition 1. For a function \(f\) given on the interval \([a,b]\), the Caputo derivative of fractional order \(r\) is defined as \begin{equation*} D^rf(t)=\frac{1}{\Gamma(n-r)}\int_0^t(t-s)^{n-r-1}f^{(n)}(s)ds, \quad n=[r]+1, \end{equation*} where \([r]\) denotes the integer part of \(r\).

Definition 2. The Riemann-Liouville fractional integral of order \(r\) for a function \(f\) is defined as \begin{equation*} I^rf(t)=\frac{1}{\Gamma(r)}\int_0^t (t-s)^{r-1}f(s)ds, \quad r>0, \end{equation*} where \([r]\) denotes the integer part of \(r\).

Lemma 1. Let \(r>0\). Then the differential equation \(D^rx(t)=0\) has solutions \begin{equation*} x(t)=c_0+c_1t+c_2t^2+\cdots+c_{n-1}t^{n-1}, \end{equation*} where \(c_i \in \mathbb{R}\), \(i=0,1,2, . . .,n\), \(n=[r]+1\).

Lemma 2. Let \(r>0\). Then \begin{equation*} I^r(D^rx)(t)=x(t)+c_0+c_1t+c_2t^2+ . . . +c_{n-1}t^{n-1}, \end{equation*} where \(c_i \in \mathbb{R}\), \(i=0,1,2, . . ., n\), \(n=[r]+1\).

For finding a solution of the problem (1)-(2) , we first consider the following fractional differential equation
\begin{eqnarray}\label{21} \left. \begin{array}{ll} -D_a^\beta u(t)=h(t), \end{array} \right. \end{eqnarray}
(3)
with the boundary condition (2).

Let we define \( \Delta:=b-a-\sum_{i=1}^{m-2}\xi_{i}(\eta_i-a)\).

Lemma 3. Let \(\beta \in (1, 2]\) and \(t\in [a,b]\). The boundary value problem (3)-(2) has a unique solution \(u\) in the form \begin{eqnarray*} u(t)&=&\left[1+\frac{ \sum_{i=1}^{m-2} \xi_{i}-1}{\Delta}(t-a)\right]\delta_1+\frac{t-a}{\Delta}\delta_2 +\int_a^b G(t,s)h(s)ds, \end{eqnarray*} where

\begin{eqnarray}\label{22} G(t,s)=\frac{1}{\Gamma(\beta)}\left\{ \begin{array}{ll} -(t-s)^{\beta-1}+\frac{t-a}{\Delta}\left[(b-s)^{\beta -1}-\sum_{j=i}^{m-2}\xi_j (\eta_j -s)^{\beta -1}\right], & \hbox{\(s\leq t\),&nbsp \(\eta_{i-1}< s\leq \eta_i\);} \\ \frac{t-a}{\Delta}\left[(b-s)^{\beta-1}-\sum_{j=i}^{m-2}\xi_j(\eta_j-s)^{\beta-1}\right], & \hbox{\(t\leq s\), &nbsp \(\eta_{i-1}< s\leq \eta_i\),} \end{array} \right. \end{eqnarray}
(4)
\(i=1,2,\ldots m-2.\)

Proof. The equation \(D^\beta u(t)=-h(t)\) has a unique solution

\begin{equation}\label{23} u(t)=-\frac{1}{\Gamma(\beta)}\int_a^t (t-s)^{\beta-1}h(s)ds+c_0+c_1(t-a), \end{equation}
(5)
where \(c_0, c_1 \in \mathbb{R}\).

By \(u(a)=\delta_1\) and \( u(b)=\sum\limits_{i=1}^{m-2}\xi u(\eta_i)+\delta_2\), we have \(c_0=\delta_1\) and

\[ c_1=\frac{1}{\Delta}\left(-\frac{1}{\Gamma(\beta)}\sum_{i=1}^{m-2} \xi_{i} \int_a^{\eta_i} (\eta_i-s)^{\beta-1}h(s)ds+\frac{1}{\Gamma(\beta)}\int_a^b(b-s)^{\beta-1}h(s)ds+\delta_1 \left[\sum_{i=1}^{m-2}\xi_{i}-1\right] +\delta_2 \right).\] Substituting \(c_0, c_1\) into equation (5) we find, \begin{align*} u(t)=&\left[1+\frac{ \sum_{i=1}^{m-2}\xi_i -1}{\Delta}(t-a)\right]\delta_1+\frac{t-a}{\Delta}\delta_2 -\frac{1}{\Gamma(\beta)}\left(\int_a^t (t-s)^{\beta -1}h(s)ds+\frac{t-a}{\Delta} \sum_{i=1}^{m-2} \xi_i \int_a^{\eta_i} (\eta_i-s)^{\beta-1}h(s)ds\right.\\ &\left.-\frac{t-a}{\Delta}\int_a^b (b-s)^{\beta-1} h(s)ds\right)\\ =& \left[1+\frac{ \sum_{i=1}^{m-2}\xi_i -1}{\Delta}(t-a)\right]\delta_1+\frac{t-a}{\Delta}\delta_2+\int_a^b G(t,s)h(s)ds, \end{align*} where \(G(t,s)\) is defined by (4). The proof is completed.

Lemma 4. If \( 0< \sum_{i=1}^{m-2}\xi_i < 1\), then

  • i) \(\Delta>0\),
  • ii) \( (b-s)^{\beta-1}- \sum\limits_{j=i}^{m-2}\xi_j(\eta_j – s)^{\beta -1 }>0\).

Proof.

  • i) We can easily see that \begin{align*} &\eta_i < b ,\\ &\xi_i(\eta_i-a)-\sum_{i=1}^{m-2}\xi_i(b-a),\\ &b-a-\sum_{i=1}^{m-2}\xi_i(\eta_i-a)>b-a-\sum_{i=1}^{m-2}\xi_i(b-a)=(b-a)\left[1-\sum_{i=1}^{m-2}\xi_{i}\right]. \end{align*} If \( 1-\sum_{i=1}^{m-2}\xi_i>0\), then \( b-a- \sum_{i=1}^{m-2}\xi_i (\eta_i -a)>0\). So we have \(\Delta>0\).
  • ii) Since \( 0b-s>(b-s)^{\beta-1}\) and \(\eta_i – a>\eta_i -s > (\eta_i -s)^{\beta -1}\). Also we know that \((\eta_j-s)^{\beta-1}< (b-s)^{\beta -1}\). Thus we have \begin{eqnarray*} \sum_{j=i}^{m-2}\xi_j(\eta_j-s)^{\beta-1}< \sum_{j=i}^{m-2}\xi_j(b-s)^{\beta-1} \leq (b-s)^{\beta-1}\sum_{i=1}^{m-2}\xi_i 0. \)

Remark 1. Let \(G(t,s)\) be the Green’s function for the problem (1)-(2). It is easy to find that

\begin{align} \label{new} \int_a^b \vert G(t,s)\vert ds &\leq \frac{1}{\Gamma(\beta)}\int_a^t (t-s)^{\beta-1}ds+\frac{t-a}{\Gamma(\beta)\Delta}\sum_{i=1}^{m-2}\xi_i \int_a^{\eta_i}(\eta_i-s)^{\beta -1}ds+ \frac{t-a}{\Delta \Gamma(\beta)}\int_a^b (b-s)^{\beta -1}ds\notag\\ &= \frac{(t-a)^\beta}{\Gamma(\beta +1)}+\frac{t-a}{\Delta \Gamma(\beta +1)}\sum_{i=1}^{m-2}\xi_i (\eta_i -a)^\beta +\frac{t-a}{\Delta \Gamma(\beta +1)}(b-a)^\beta\notag\\ &\leq \frac{(b-a)^\beta}{\Gamma(\beta +1)}+\frac{b-a}{\Delta \Gamma(\beta +1)}\sum_{i=1}^{m-2}\xi_i (\eta_i -a)^\beta+\frac{(b-a)^{\beta +1}}{\Delta \Gamma(\beta +1)}:=M. \end{align}
(6)

Remark 2. Let \(p(t)\in L^1[a,b]\) and \(w(t)\) be a solution of the problem

\begin{eqnarray}\label{25} \begin{cases} \begin{array}{ll} D_a^\beta u(t)+p(t)=0\\ u(a)=0, \quad u(b)= \sum_{i=1}^{m-2}\xi_i u(\eta_i) \end{array} \end{cases} \end{eqnarray}
(7)
then \( w(t)=\int_a^b G(t,s)p(s)ds\).

The following fixed point theorem is fundamental and important to the proof of our main result.

Theorem 1. [14][ Schauder-Tychonoff Fixed Point Theorem] Let \( X \) be a Banach space. Assume that \( K \) is a closed, bounded, convex subset of \( X \). If \( T:K\longrightarrow K \) is compact, then \( T \) has a fixed point in \(K.\)

3. Existence Result

We also make the following assumption throughout this paper.
  • (H1) There exists a nonnegative function \(p \in L^1[a,b]\) and \( \int_a^b p(t)dt>0\) such that \(f(t,u)\geq – p(t)\) for all \((t,u)\in [a,b]\times \mathbb{R}\).
  • (H2) \(f(t,u)\neq 0\), for \((t,u)\in[a,b]\times \mathbb{R}\).
Let \(B=\mathbb{C}([a,b],\mathbb{R})\) denote the Banach space of all continuous function from \([a,b]\) into \(\mathbb{R}\) endowed with the usual norm \(\Vert u \Vert = sup\{\vert u(t) \vert: t\in[a,b]\}\).

First we shall show that the following fractional equation

\begin{equation}\label{31} D_a^\beta u(t)+F(t,u^*(t))=0, \quad t\in [a,b] \end{equation}
(8)
with the boundary condition (2) has a solution, where \(F:[a,b]\times \mathbb{R}\rightarrow \mathbb{R}\)
\begin{eqnarray}\label{32} F(t,z)=\left\{ \begin{array}{ll} f(t,z)+p(t), & \hbox{\(z\geq 0\),} \\ f(t,0)+p(t), & \hbox{\(z\leq 0\),} \end{array} \right. \end{eqnarray}
(9)
and \(u^*(t)=max\{(u-w)(t),0\}\) such that \(w\) is the unique solution of the problem (7). We define an operator \(T:B \rightarrow B\) associated with the problem (8)-(2) as
\begin{equation} \label{33} (Tu)(t)=\left[1+\frac{ \sum_{i=1}^{m-2}\xi_i -1}{\Delta}(t-a)\right]\delta_1+\frac{t-a}{\Delta}\delta_2+\int_a^b G(t,s)F(s,u^*(s))ds \end{equation}
(10)
where \(G(t,s)\) is given by (4). The existence of a fixed point for the operator \(T\) implies the existence of a solution for the problem (8)-(2).

Theorem 2. Assume that (H1)-(H2) are satisfied. If \( K > 0 \) satisfies \(\left[1+\frac{ \sum_{i=1}^{m-2}\xi_i -1}{\Delta}(b-a)\right]\delta_1+ \frac{b-a}{\Delta}\delta_2 + L M \leq K \) where \(L\geq max\{\left|F(t,u) \right|:t \in [a,b], \ \left|u \right|\leq K \} \) and \(M\) is given in (6) then the problem (8)-(9) has a solution \(u(t)\).

Proof. Let we define \( P:=\{u\in B:\left\|u \right\|\leq K \}\). It can be easily seen that \( P \) is a closed, bounded and convex subset of \( B \) to which Schauder fixed point theorem is applicable. Define \( T:P\longrightarrow B\) by (10). It is easily seen that \( T:P\longrightarrow B\) is continuous. Claim \(T:P \longrightarrow P\). Let \( u \in P \). Consider \(u^*(t) \leq u(t) \leq K\) for all \(t \in [a, b]\). Then \begin{align*} \left|Tu(t) \right|&=\left|\left[1+\frac{ \sum_{i=1}^{m-2}\xi_i -1}{\Delta}(t-a)\right]\delta_1+\frac{t-a}{\Delta}\delta_2+\int_a^b G(t,s)F(s,u^*(s))ds \right|\\ &\leq \left[1+\frac{ \sum_{i=1}^{m-2}\xi_i -1}{\Delta}(b-a)\right]\delta_1+\frac{b-a}{\Delta}\delta_2+ L M \leq K \end{align*} for all \( t \in [a,b].\) This implies that \( \left\|Tu\right\|\leq K \). Hence \( T:K\longrightarrow K.\) Using the Arzela-Ascoli theorem it can be shown that \( T:K\longrightarrow K \) is a compact operator. Hence \( T \) has a fixed point \( u \) in \(P\) by the Schauder-Tychonov theorem. This implies that \(u\) is a solution of the problem (8)-(2).

Lemma 5. \(u^*(t)\) is the solution of the boundary value problem (1)-(2) with \(u(t)>w(t)\) for all \(t\in[a,b]\) if and only if \(u=u^*+w\) is the positive solution of the boundary value problem (8)-(2).

Proof. Let \(u(t)\) is the solution of the boundary value problem (8)-(2). Then \begin{align*} u(t)&=\left[1+\frac{ \sum_{i=1}^{m-2}\xi_i-1}{\Delta}(t-a)\right]\delta_1+\frac{t-a}{\Delta}\delta_2 + \frac{1}{\Gamma(\beta)}\int_a^b G(t,s)F(s,u^*(s))ds\\ &=\left[1+\frac{ \sum_{i=1}^{m-2}\xi_i-1}{\Delta}(t-a)\right]\delta_1+\frac{t-a}{\Delta}\delta_2 + \frac{1}{\Gamma(\beta)}\int_a^b G(t,s)(f(s,u^*(s))+p(s))ds\\ &=\left[1+\frac{ \sum_{i=1}^{m-2}\xi_i-1}{\Delta}(t-a)\right]\delta_1+\frac{t-a}{\Delta}\delta_2 + \frac{1}{\Gamma(\beta)}\int_a^b G(t,s)f(s,(u-w)(s))ds+\frac{1}{\Gamma(\beta)}\int_a^b G(t,s)p(s)ds\\ &=\left[1+\frac{ \sum_{i=1}^{m-2}\xi_i-1}{\Delta}(t-a)\right]\delta_1+\frac{t-a}{\Delta}\delta_2 + \frac{1}{\Gamma(\beta)}\int_a^b G(t,s)f(s,(u-w)(s))ds+w(t) \end{align*} or \begin{eqnarray*} u(t)-w(t)=&\left[1+\frac{ \sum_{i=1}^{m-2}\xi_i-1}{\Delta}(t-a)\right]\delta_1+\frac{t-a}{\Delta}\delta_2 + \frac{1}{\Gamma(\beta)}\int_a^b G(t,s)f(s,(u-w)(s))ds \end{eqnarray*} and hence we get \begin{eqnarray*} u^*(t)&= \left[1+\frac{ \sum_{i=1}^{m-2}\xi_i-1}{\Delta}(t-a)\right]\delta_1+\frac{t-a}{\Delta}\delta_2+\frac{1}{\Gamma(\beta)}\int_a^b G(t,s)f(s,u^*(s))ds. \end{eqnarray*} In other words, if \(u^*\) is a solution of the boundary value problem (1)-(2) then we get \begin{eqnarray*} D_a^\beta (u^*(t)+w(t))=D_a^\beta u^*(t)+ D_a^\beta w(t) = -f(t,u^*(t))-p(t)= -[f(t,u^*(t))+p(t)]= – F(t,u^*(t)), \end{eqnarray*} which implies that \begin{equation*} D_a^\beta u(t)=-F(t,u^*(t)). \end{equation*} Also from the boundary conditions, we easily see that \begin{equation*} u^*(a)=u(a)-w(a)=u(a)-0=\delta_1, \end{equation*} i.e., \( u(0)=\delta_1\) and \begin{align*} &u^*(b)=\sum_{i=1}^{m-2}\xi_i u^*(\eta_i)+\delta_2\\ &u(b)-w(b)= \sum_{i=1}^{m-2}\xi_i (u(\eta_i)-w(\eta_i))+\delta_2=\sum_{i=1}^{m-2}\xi_i u(\eta_i)-\sum_{i=1}^{m-2}\xi_i w(\eta_i)+\delta_2, \end{align*} i.e., \[u(b)= \sum_{i=1}^{m-2}\xi_i u(\eta_i)+\delta_2.\] Therefore \(u(t)\) is a solution of the boundary value problem (8)-(9).

Example 1. Consider the following fractional boundary value problem

\begin{equation}\label{3.4} D^\frac{4}{3}u(t)+f(t,u(t))=0,\ \ \ t\in(0, 1), \end{equation}
(11)
\begin{equation}\label{3.5} \begin{array}{ll} u(0)=1,\ \ \ \ u(1)=\frac{1}{4}u(\frac{1}{2})+\frac{1}{2}u(\frac{1}{4})+\frac{1}{8}u(\frac{1}{5})-1 \end{array} \end{equation}
(12)
with the function \(f(t,u(t)) = \frac{e^t}{1+t^2}\sin(u(t))\).

Choosing \(p(t)= e^t\) we get \(\int_{0}^{1} e^t dt =e-1 > 0\), so it is easy to check that the assumptions (H1)-(H2) hold. Calculating \(\Delta = \frac{29}{40}, M\cong2.9\) and seeing \(|F(t,u)| < 2 e = L\) such that \(|u|\leq K\) where \(K=20\), we can easily verify that

\[\left[1+\frac{ \sum_{i=1}^{m-2}\xi_i -1}{\Delta}(b-a)\right]\delta_1+ \frac{b-a}{\Delta}\delta_2 + L M \cong 15,3 \leq 20.\] Then applying Theorem 2 the problem (11)-(12) has a solution \(u(t)\).

Author Contributions

All authors contributed equally to the writing of this paper. All authors read and approved the final manuscript.

Conflicts of Interest

”The authors declare no conflict of interest.”

References:

  1. Ahmad, B. (2010). Existence of solutions for irregular boundary value problems of nonlinear fractional differential equations. Applied Mathematics Letters, 23(4), 390-394. [Google Scholor]
  2. Bai, Z. (2010). On positive solutions of a nonlocal fractional boundary value problem. Nonlinear Analysis: Theory, Methods & Applications, 72(2), 916-924. [Google Scholor]
  3. Bai, Z., & Lü, H. (2005). Positive solutions for boundary value problem of nonlinear fractional differential equation. Journal of Mathematical Analysis and Applications, 311(2), 495-505. [Google Scholor]
  4. Kilbas, A. A., Srivastava, H. M., & Trujillo, J. J. (2006). Theory and Applications of Fractional Differential Equations(Vol. 204). elsevier. [Google Scholor]
  5. Miller, K. S., & Ross, B. (1993). An Introduction to the Fractional Calculus and Fractional Differential Equations. Wiley, New York. [Google Scholor]
  6. Podlubny, I. (1999). Fractional Differential Equations, Mathematics in Science and Engineering. Academic Press New York. [Google Scholor]
  7. Liang, S., & Zhang, J. (2009). Positive solutions for boundary value problems of nonlinear fractional differential equation. Nonlinear Analysis: Theory, Methods & Applications, 71(11), 5545-5550. [Google Scholor]
  8. Zhang, S. (2006). Positive solutions for boundary-value problems of nonlinear fractional differential equations. Electronic Journal of Differential Equations, 36(2), 1-12. [Google Scholor]
  9. Zhong, W., & Lin, W. (2010). Nonlocal and multiple-point boundary value problem for fractional differential equations. Computers & Mathematics with Applications, 59(3), 1345-1351. [Google Scholor]
  10. Li, F., Zhang, Y., & Li, Y. (2008). Sign-changing solutions on a kind of fourth-order Neumann boundary value problem. Journal of Mathematical Analysis and Applications, 344(1), 417-428. [Google Scholor]
  11. Li, Y., & Li, F. (2008). Sign-changing solutions to second-order integral boundary value problems. Nonlinear Analysis: Theory, Methods & Applications, 69(4), 1179-1187. [Google Scholor]
  12. Liu, Z., Ding, Y., Liu, C., & Zhao, C. (2020). Existence and uniqueness of solutions for singular fractional differential equation boundary value problem with p-Laplacian. Advances in Difference Equations, 2020, Article No. 83. https://doi.org/10.1186/s13662-019-2482-9. [Google Scholor]
  13. Xu, X. (2004). Multiple sign-changing solutions for some m-point boundary-value problems. Electronic Journal of Differential Equations, 89, 1-14. [Google Scholor]
  14. Tychonoff, A. (1935). Ein fixpunktsatz. Mathematische Annalen, 111(1), 767-776. [Google Scholor]
  15. Ahmad, B. (2017). Sharp estimates for the unique solution of two-point fractional-order boundary value problems. Applied Mathematics Letters, 65, 77-82. [Google Scholor]