Skip to main content

Existence–uniqueness and monotone iteration of positive solutions to nonlinear tempered fractional differential equation with p-Laplacian operator

Abstract

In this paper, without requiring the complete continuity of integral operators and the existence of upper–lower solutions, by means of the sum-type mixed monotone operator fixed point theorem based on the cone \(P_{h}\), we investigate a kind of p-Laplacian differential equation Riemann–Stieltjes integral boundary value problem involving a tempered fractional derivative. Not only the existence and uniqueness of positive solutions are obtained, but also we can construct successively sequences for approximating the unique positive solution. As an application of our fundamental aims, we offer a realistic example to illustrate the effectiveness and practicability of the main results.

1 Introduction

In this paper, we devote our study to the kind of p-Laplacian differential equations Riemann–Stieltjes integral boundary value problems involving tempered fractional derivatives as follows:

$$ \textstyle\begin{cases} {{}_{0}^{R}\mathbb{D}}_{t}^{\alpha _{2},\lambda } (\varphi _{p} ({{}_{0}^{R} \mathbb{D}}_{t}^{\alpha _{1},\lambda } u(t) ) )=f (t,u(t),u(t) )+g (t,u(t) ),\quad t\in [0,1], \\ u(0)=u'(0)=\cdots =u^{(n-2)}(0)=0, \\ \varphi _{p} ({{}_{0}^{R}\mathbb{D}}_{t}^{\alpha _{1},\lambda }u )(0)=0, \\ u(1)=\beta \int _{0}^{1}e^{-\lambda (1-t)}u(t)\,dt, \\ {{}_{0}^{R}\mathbb{D}}_{t}^{\gamma _{1},\lambda } (\varphi _{p} ({{}_{0}^{R}\mathbb{D}}_{t}^{\alpha _{1},\lambda } u ) )(1)= \int _{0}^{\eta }a(t){{}_{0}^{R}\mathbb{D}}_{t}^{\gamma _{2},\lambda } [\varphi _{p} ({{}_{0}^{R}\mathbb{D}}_{t}^{\alpha _{1},\lambda } u(t) ) ]\,dA(t), \end{cases} $$
(1.1)

where \(n-1<\alpha _{1}\leq n\), \(1<\alpha _{2}\leq 2\), \(0<\gamma _{2}<\gamma _{1}<\alpha _{2}-1\), \(\beta <\alpha _{1}\) and \(\lambda >0\) is constant, \(\varphi _{p}\) is a p-Laplacian operator. \({}_{0}^{R}\mathbb{D}_{t}^{\alpha _{1},\lambda }\) are tempered fractional derivatives, which are defined by

$$ _{0}^{R}\mathbb{D}_{t}^{\alpha _{1},\lambda }u(t)=e^{-\lambda t}{{}_{0}^{R}}D_{t}^{ \alpha _{1}} \bigl(e^{\lambda t}u(t)\bigr). $$
(1.2)

Here, \({{}_{0}^{R}}D_{t}^{\alpha _{1}}\) denotes the standard Riemann–Liouville fractional derivative

$$ {{}_{0}^{R}}D_{t}^{\alpha _{1}}u(t)= \frac{d^{n}}{dt^{n}} \bigl({}_{0}I_{t}^{n- \alpha _{1}}u(t) \bigr), $$
(1.3)

where \({}_{0}I_{t}^{\upsilon }\) for \(\upsilon >0\) is the fractional integral operator of order υ defined by

$$ _{0}I_{t}^{\upsilon }\psi = \frac{1}{\varGamma (\upsilon )} \int _{0}^{t}(t-s)^{ \upsilon -1}\psi (s)\,ds. $$
(1.4)

A is a function of a bounded variation, \(\int _{0}^{\eta }a(t){{}_{0}^{R}\mathbb{D}}_{t}^{\gamma _{2},\lambda } [\varphi _{p} ({{}_{0}^{R}\mathbb{D}}_{t}^{\alpha _{1},\lambda } u(t) ) ]\,dA(t)\) denotes a Riemann–Stieltjes integral with respect to A. By using the sum-type mixed monotone fixed theorem based on the cone \(P_{h}\), we show the existence and uniqueness of positive solutions for the p-Laplacian differential system (1.1).

In recent years, many theories and experiments have shown that a large number of abnormal phenomena that occurs in the applied science and engineering can be well described by fractional calculus. Especially, fractional differential equations have been proved to be powerful tools in the modeling of various phenomena in various fields of science and engineering, for example fluid mechanics, physics and heat conduction; see for instance [16]. Meanwhile, it is well known that the p-Laplacian operator is also used in analyzing biology, physics, mechanics and the related fields of mathematical modeling; see [714]. In [7], for studying the turbulent flow in porous media, Leibenson introduced the p-Laplacian differential equation as follows:

$$ \bigl(\varphi _{p} \bigl(u'(t) \bigr) \bigr)'=f \bigl(t,u(t),u'(t) \bigr),\quad t\in (0,1), $$
(1.5)

where \(\varphi _{p}(s)=|s|^{p-2}s\), \(p>1\). Motivated by Leibenson’s work, Guo et al. [8] studied the existence of a solution for an ordinary differential equation m-point boundary value problem with p-Laplacian operator. Lu et al. [9] investigated a fractional differential equation for a two points boundary value problem involving the p-Laplacian operator as follows:

$$ \textstyle\begin{cases} D_{0+}^{\beta } (\varphi _{p} (D_{0+}^{\alpha }u(t) ) )=f (t,u(t) ), \quad 0\leq t\leq 1; \\ u(0)=u'(0)=u'(1)=0; \\ D_{0+}^{\alpha }u(0)=D_{0+}^{\alpha }u(1)=0, \end{cases} $$
(1.6)

where \(2<\alpha \leq 3\), \(1<\beta \leq 2\) and \(\varphi _{p}(s)=|s|^{p-2}s\). \(D_{0+}^{\alpha }\), \(D_{0+}^{\beta }\) are standard Riemann–Liouville fractional derivatives. By employing the Guo–Krasnosel’skii fixed-point theorem and upper–lower solutions method, the existence of positive solutions was obtained.

In [10], Ren, Li and Zhang studied the existence of maximum and minimum solutions for the following nonlocal p-Laplacian fractional differential system:

$$ \textstyle\begin{cases} -D_{t}^{\beta _{1}} (\varphi _{p_{1}} (-D_{t}^{\alpha _{1}} x_{1} ) )(t)=f_{1} (x_{1}(t),x_{2}(t) ), \\ -D_{t}^{\beta _{2}} (\varphi _{p_{2}} (-D_{t}^{\alpha _{2}} x_{2} ) )(t)=f_{2} (x_{1}(t),x_{2}(t) ), \\ x_{1}(0)=0,\qquad D_{t}^{\alpha _{1}}x_{1}(0)=D_{t}^{\alpha _{1}}x_{1}(1)=0,\qquad x_{1}(1)= \int _{0}^{1}x_{1}(t)\,dA_{1}(t), \\ x_{2}(0)=0,\qquad D_{t}^{\alpha _{2}}x_{2}(0)=D_{t}^{\alpha _{2}}x_{2}(1)=0,\qquad x_{2}(1)= \int _{0}^{1}x_{2}(t)\,dA_{2}(t), \end{cases} $$
(1.7)

where \(D_{t}^{\alpha _{i}}\), \(D_{t}^{\beta _{i}}\) are the standard Riemann–Liouville derivatives satisfying \(1<\alpha _{i}\), \(\beta _{i}<2\), \(\int _{0}^{1}x_{i}(t)\,dA_{i}(t)\) denotes a Riemann–Stieltjes integral and \(A_{i}\) is a function of bounded variation, \(\varphi _{p_{i}}\) is a p-Laplacian operator. By using the monotone iterative technique, some new results as regards the existence of maximal and minimal solutions were established, and the estimation of the lower and upper bounds of the maximum and minimum solutions was also derived.

Recently, in [15], we investigated the conformable differential equation with p-Laplacian operator as follows:

$$ \textstyle\begin{cases} T_{\alpha }^{0+} (\varphi _{p} (T_{\alpha }^{0+}u(t) ) )=f (t,u(t),T_{\alpha }^{0+}u(t) ), \\ u^{(i)}(0)=0,\qquad [\varphi _{p}(T_{\alpha }^{0+}u)]^{(i)}(0)=0, \\ [T_{\beta }^{0+}u(t)]_{t=1}=0, \qquad [T_{\beta }^{0+} (\varphi _{p} (T_{\alpha }^{0+}u(t) ) ) ]_{t=1}=0, \end{cases} $$
(1.8)

where \(n-1\leq \alpha < n\) and \(T_{\alpha }^{0^{+}}\) is a new fractional derivative called “the conformable fractional derivative”. By using the Guo–Krasnosel’skii fixed point theorem, some new existence conclusions of positive solutions were obtained to the boundary value problem (1.8).

In [16], we continued to investigate the existence of multiple positive solutions for high order Riemann–Liouville fractional differential equation involving the p-Laplacian operator as follows:

$$ \textstyle\begin{cases} {{}_{0}^{R}D}_{t}^{\alpha } (\varphi _{p} ({{}_{0}^{R}D}_{t}^{\alpha }u(t) ) )=f (t,u(t),{{}_{0}^{R}D}_{t}^{\alpha }u(t) ),\quad 0\leq t \leq 1; \\ u^{(i)}(0)=0,\qquad [\varphi _{p}({{}_{0}^{R}D}_{t}^{\alpha }u)]^{(i)}(0)=0,\quad i=0,1,2,\ldots,n-2; \\ [{{}_{0}^{R}D}_{t}^{\beta }u(t)]_{t=1}=0, \quad 0< \beta \leq \alpha -1; \\ [{{}_{0}^{R}D}_{t}^{\beta } (\varphi _{p} ({{}_{0}^{R}D}_{t}^{\alpha }u(t) ) ) ]_{t=1}=0; \end{cases} $$
(1.9)

where \(n-1<\alpha \leq n\), \({{}_{0}^{R}D}_{t}^{\alpha }\) is the standard Riemann–Liouville fractional derivative, \(\varphi _{p}\) is the p-Laplacian operator. By means of the Leggett–Williams fixed point theorem and a functional-type cone expansion-compression fixed point theorem, not only the existence of two positive solutions was obtained, but also some sufficient conditions for the existence of at least three positive solutions was established.

In addition, Zhang et al. [17] investigated the eigenvalue problem for a kind of singular fractional differential equation Riemann–Stieltjes integral boundary value problem involving the p-Laplacian operator as follows:

$$ \textstyle\begin{cases} -D_{t}^{\beta } (\varphi _{p} (D_{t}^{\alpha }x(t) ) )= \lambda f (t,x(t) ),\quad 0\leq t\leq 1, \\ x(0)=0,\qquad D_{t}^{\alpha }x(0)=0, \\ x(1)=\int _{0}^{1}x(s)\,dA(s), \end{cases} $$
(1.10)

where \(D_{t}^{\beta }\) and \(D_{t}^{\alpha }\) are standard Riemann–Liouville fractional derivatives with \(0<\beta \leq 1\), \(1<\alpha \leq 2\), \(\int _{0}^{1}x(s)\,dA(s)\) is the standard Riemann–Stieltjes integral and A is a function of the bounded variation. By using the Schauder fixed point theorem and upper and lower solution methods, some new theorems on existence were obtained.

Inspired by the above work, in this paper, we investigate the existence and uniqueness of positive solutions for a p-Laplacian differential equation Riemann–Stieltjes integral boundary value problem involving a tempered fractional derivative (1.1). To the best of our knowledge, this kind of integral boundary value problem involving a tempered fractional derivative has seldom been researched up to now. Compared with other references, the present article has the following characteristics. Firstly, the tempered fractional derivative \({{}_{0}^{R}\mathbb{D}}_{t}^{\alpha ,\lambda }\) is more general than the standard fractional derivative \({{}_{0}^{R}D}_{t}^{\alpha }\). For example, letting \(\lambda =0\), it is easy to see that \({{}_{0}^{R}\mathbb{D}}_{t}^{\alpha ,\lambda }\) is equivalent to \({{}_{0}^{R}D}_{t}^{\alpha }\). Secondly, the Riemann–Stieltjes integral boundary conditions involving a tempered fractional derivative are more general cases, which cover the common integral boundary conditions as special cases. Thirdly, compared with the p-Laplacian differential system (1.8) and (1.9), in this paper, the integral operator need not be completely continuous or compact. Fourthly, in this paper, by employing the sum-type mixed monotone operators fixed points theorem, our conclusions cannot only guarantee the existence of a unique positive solution, but also construct successively sequences for approximating the unique positive solution. Finally, it is worth mentioning that some important properties of two different kernel functions rely on the parameter λ.

The rest of this paper is organized as follows. In Sect. 2, we briefly introduce some necessary basic definitions and preliminary results which will be used to prove our main results. In Sect. 3, we study the existence and uniqueness and monotone iteration of a positive solution to the p-Laplacian differential system (1.1) by means of sum-type mixed monotone fixed points theorems based on the cone \(P_{h}\). At last, in Sect. 4, we demonstrate the effectiveness and feasibility of the main results by an example.

2 Preliminaries

In the section, we first list some basic notations, concepts in ordered Banach spaces. For convenience, we refer the reader to [18, 19] for details.

Suppose that \((E,\| \cdot \| )\) is a real Banach space which is partially ordered by a cone \(P\subset E\), that is, \(x\leq y\) if and only if \(y-x\in P\). If \(x\leq y\) and \(x\neq y\), then we denote \(x< y\) or \(y>x\). By θ we denote the zero element of E. A nonempty closed convex set \(P\subset E\) is a cone if it satisfies: (i) \(x\in P\), \(\lambda \geq 0 \Rightarrow \lambda x\in P\); (ii) \(x\in P\), \(-x\in P \Rightarrow x=\theta \).

Definition 2.1

([18])

P is called normal if there exists \(M>0\) such that, for all \(x,y\in E\), \(\theta \leq x\leq y\) implies \(\|x\|\leq \|y\|\); in this case M is the infimum of such a constant, it is called the normality constant of P.

In addition, for a given \(h>\theta \), we denote by \(P_{h}\) the set \(P_{h}=\{x\in E\mid x\sim h\}\), in which is an equivalence relation, i.e., \(x\sim y\) means that there exist \(\lambda >0\) and \(\mu >0\) such that \(\lambda x\geq y\geq \mu x\) for all \(x,y\in E\).

Definition 2.2

([20])

An operator \(A:P\times P\rightarrow P\) is said to be a mixed monotone operator if \(A(x,y)\) is increasing in x and decreasing in y, i.e., \(u_{i}, v_{i}(i=1,2)\in P\), \(u_{1}< u_{2}\), \(v_{1}>v_{2}\) imply \(A(u_{1},v_{1})\leq A(u_{2},v_{2})\). An element \(x\in P\) is called a fixed point of A if \(A(x,x)=x\).

Definition 2.3

([9])

Let \(p>1\), the p-Laplacian operator is given by

$$ \varphi _{p}(x)= \vert x \vert ^{p-2}x\quad \mbox{and}\quad \varphi _{p}^{-1}=\varphi _{q},\quad \frac{1}{p}+\frac{1}{q}=1. $$

Definition 2.4

([21])

\(A:P\rightarrow P\) is said to be sub-homogeneous if it satisfies

$$ A(tx)\geq tAx,\quad \forall t\in (0,1), x\in P. $$

Lemma 2.1

([16])

Let\(h(t)\in C[0,1]\cap L^{1}[0,1]\), \(\alpha >0\), then

$$ _{0}I_{t}^{\alpha }{{}_{0}^{R}}D_{t}^{\alpha }h(t)=h(t)+c_{1}t^{\alpha -1}+c_{2}t^{ \alpha -2}+ \cdots +c_{n}t^{\alpha -n}, $$

where\(c_{i}\in R\), \(i=1, 2, 3, \ldots, n\) (\(n=[\alpha ]+1\)).

Lemma 2.2

([16])

  1. (1)

    If\(u\in L^{1}(0,1)\), \(\alpha >\beta >0\), then

    $$ {{}_{0}I}_{t}^{\alpha }{{}_{0}I}_{t}^{\beta }u(t)={{}_{0}I}_{t}^{\alpha +\beta }u(t),\qquad {{}_{0}^{R}D}_{t}^{\beta }{{}_{0}I}_{t}^{\alpha }u(t)=_{0}I_{t}^{\alpha -\beta }u(t),\qquad {{}_{0}^{R}D}_{t}^{\beta }{{}_{0}I}_{t}^{\beta }u(t)=u(t). $$
  2. (2)

    If\(\rho >0\), \(\mu >0\), then

    $$ {{}_{0}^{R}D}_{t}^{\rho }t^{\mu -1}= \frac{\varGamma (\mu )}{\varGamma (\mu -\rho )}t^{\mu -\rho -1}. $$

Lemma 2.3

Let\(g(t)\in C[0,1]\), then the unique solution of the linear problem

$$ \textstyle\begin{cases} {{}_{0}^{R}}\mathbb{D}_{t}^{\alpha _{1},\lambda }u(t)+g(t)=0,\quad n-1< \alpha \leq n; \\ u(0)=u'(0)=\cdots =u^{(n-2)}(0)=0; \\ u(1)=\beta \int _{0}^{1}e^{-\lambda (1-t)}u(t)\,dt, \quad \beta < \alpha _{1}; \end{cases} $$
(2.1)

is given by

$$ u(t)= \int _{0}^{1}H(t,s)g(s)\,ds, $$
(2.2)

where we have the Green function

$$ H(t,s)= \textstyle\begin{cases} \frac{\alpha _{1}(1-s)^{\alpha _{1}-1}(\alpha _{1}-\beta +\beta s)e^{\lambda s}t^{\alpha _{1}-1}-\alpha _{1}(\alpha _{1}-\beta )e^{\lambda s}(t-s)^{\alpha _{1}-1}}{(\alpha _{1}-\beta )\varGamma (\alpha _{1}+1)}e^{- \lambda t},& 0\leq s\leq t\leq 1; \\ \frac{\alpha _{1}(1-s)^{\alpha _{1}-1}(\alpha _{1}-\beta +\beta s)e^{\lambda s}}{(\alpha _{1}-\beta )\varGamma (\alpha _{1}+1)}e^{- \lambda t}t^{\alpha _{1}-1}, &0\leq t\leq s\leq 1. \end{cases} $$
(2.3)

Proof

For the system (2.1), by using Lemma 2.1, we get

$$ e^{\lambda t}u(t)=- \int _{0}^{t} \frac{(t-s)^{\alpha _{1}-1}}{\varGamma (\alpha _{1})}e^{\lambda s}g(s)\,ds+c_{1}t^{ \alpha _{1}-1}+c_{2}t^{\alpha _{1}-2}+ \cdots+c_{n}t^{\alpha _{1}-n}. $$

Furthermore, the boundary conditions \(u(0)=u'(0)=\cdots =u^{(n-2)}(0)=0\) imply that \(c_{n}=c_{n-1}=c_{n-2}=\cdots=c_{3}=c_{2}=0\). Thus, we have

$$ e^{\lambda t}u(t)=- \int _{0}^{t} \frac{(t-s)^{\alpha _{1}-1}}{\varGamma (\alpha _{1})}e^{\lambda s}g(s)\,ds+c_{1}t^{ \alpha _{1}-1}. $$
(2.4)

Integrating both sides of Eq. (2.4) from 0 to 1, we see that

$$ \begin{aligned}[b] \int _{0}^{1}e^{\lambda t} u(t)\,dt&=- \int _{0}^{1}\biggl( \int _{0}^{t} \frac{(t-s)^{\alpha _{1}-1}}{\varGamma (\alpha _{1})}e^{\lambda s}g(s)\,ds \biggr)\,dt+c_{1} \int _{0}^{1}t^{\alpha _{1}-1}\,dt \\ &=- \int _{0}^{1}\frac{e^{\lambda s}}{\varGamma (\alpha _{1})}g(s)\,ds \int _{s}^{1}(t-s)^{ \alpha _{1}-1}\,dt+ \frac{c_{1}}{\alpha _{1}} \\ &=\frac{c_{1}}{\alpha _{1}}- \int _{0}^{1} \frac{(1-s)^{\alpha _{1}}}{\varGamma (\alpha _{1}+1)}e^{\lambda s}g(s)\,ds. \end{aligned} $$
(2.5)

Letting \(t=1\) in (2.4), we obtain

$$ e^{\lambda }u(1)=- \int _{0}^{1} \frac{(1-s)^{\alpha _{1}-1}}{\varGamma (\alpha _{1})}e^{\lambda s}g(s)\,ds+c_{1}. $$
(2.6)

Combining the integral boundary value condition \(u(1)=\beta \int _{0}^{1}e^{-\lambda (1-t)}u(t)\,dt\), (2.6) and (2.5), we can clearly see that

$$ c_{1}= \int _{0}^{1} \frac{\alpha _{1}^{2}(1-s)^{\alpha _{1}-1}-\alpha _{1}\beta (1-s)^{\alpha _{1}}}{(\alpha _{1}-\beta )\varGamma (\alpha _{1}+1)}e^{ \lambda s}g(s)\,ds. $$
(2.7)

Finally, by simply substituting (2.7) into (2.4),

$$ \begin{aligned} u(t)&= \int _{0}^{1} \frac{[\alpha _{1}^{2}(1-s)^{\alpha _{1}-1}-\alpha _{1}\beta (1-s)^{\alpha _{1}}]t^{\alpha _{1}-1}}{(\alpha _{1}-\beta )\varGamma (\alpha _{1}+1)}e^{- \lambda t}e^{\lambda s}g(s)\,ds \\ &\quad {}- \int _{0}^{t} \frac{\alpha _{1}(\alpha _{1}-\beta )(t-s)^{\alpha _{1}-1}}{(\alpha _{1}-\beta )\varGamma (\alpha _{1}+1)}e^{- \lambda t}e^{\lambda s}g(s)\,ds \\ &= \int _{0}^{1}H(t,s)g(s), \end{aligned} $$

where the Green function \(H(t,s)\) is defined as (2.3). □

Lemma 2.4

If\(\widetilde{g}\in C[0,1]\)is given, then thep-Laplacian tempered fractional differential equation integral boundary value problem

$$ \textstyle\begin{cases} {{}_{0}^{R}\mathbb{D}}_{t}^{\alpha _{2},\lambda } (\varphi _{p} ({{}_{0}^{R} \mathbb{D}}_{t}^{\alpha _{1},\lambda } u(t) ) )=\widetilde{g}(t),\quad t \in [0,1], \\ u(0)=u'(0)=\cdots =u^{(n-2)}(0)=0, \\ \varphi _{p} ({{}_{0}^{R}\mathbb{D}}_{t}^{\alpha _{1},\lambda }u )(0)=0, \\ u(1)=\beta \int _{0}^{1}e^{-\lambda (1-t)}u(t)\,dt, \\ {{}_{0}^{R}\mathbb{D}}_{t}^{\gamma _{1},\lambda } (\varphi _{p} ({{}_{0}^{R}\mathbb{D}}_{t}^{\alpha _{1},\lambda } u ) )(1)= \int _{0}^{\eta }a(s){{}_{0}^{R}\mathbb{D}}_{t}^{\gamma _{2},\lambda } [\varphi _{p} ({{}_{0}^{R}\mathbb{D}}_{t}^{\alpha _{1},\lambda } u(s) ) ]\,dA(s), \end{cases} $$
(2.8)

has a unique integral formal solution

$$ u(t)= \int _{0}^{1}H(t,s)\varphi _{q} \biggl( \int _{0}^{1}G(s,\tau ) \widetilde{g}(\tau )\,d\tau \biggr)\,ds, $$
(2.9)

where\(H(t,s)\)is given as (2.3), \(G(t,s)\)is a Green function and

$$ G(t,s)=G_{1}(t,s)+ \frac{t^{\alpha _{2}-1}e^{-\lambda t}}{\Delta \varGamma (\alpha _{2}-\gamma _{2})} \int _{0}^{\eta }a(t)G_{2}(t,s)\,dA(t), $$
(2.10)

in which

$$\begin{aligned}& G_{1}(t,s)=\frac{e^{\lambda (s-t)}}{\varGamma (\alpha _{2})} \textstyle\begin{cases} (1-s)^{\alpha _{2}-\gamma _{1}-1}t^{\alpha _{2}-1}-(t-s)^{\alpha _{2}-1},& 0 \leq s\leq t\leq 1, \\ (1-s)^{\alpha _{2}-\gamma _{1}-1}t^{\alpha _{2}-1},& 0\leq t\leq s \leq 1, \end{cases}\displaystyle \\& G_{2}(t,s)=\frac{e^{\lambda (s-t)}}{\varGamma (\alpha )} \textstyle\begin{cases} (1-s)^{\alpha _{2}-\gamma _{1}-1}t^{\alpha _{2}-\gamma _{2}-1}-(t-s)^{ \alpha _{2}-\gamma _{2}-1},& 0\leq s\leq t\leq 1, \\ (1-s)^{\alpha _{2}-\gamma _{1}-1}t^{\alpha _{2}-\gamma _{2}-1},& 0 \leq t\leq s\leq 1, \end{cases}\displaystyle \end{aligned}$$

and

$$ \Delta =\frac{e^{-\lambda }}{\varGamma (\alpha _{2}-\gamma _{1})}- \frac{\delta }{\varGamma (\alpha _{2}-\gamma _{2})},\quad \delta = \int _{0}^{\eta }e^{-\lambda s}s^{\alpha _{2}-\gamma _{2}-1}a(s)\,dA(s). $$

Proof

From Lemma 2.1, integrating both sides of the first equation of (2.8), we obtain

$$ \begin{aligned} e^{\lambda t}\varphi _{p} \bigl({{}_{0}^{R}\mathbb{D}}_{t}^{ \alpha _{1},\lambda } u(t) \bigr)&={{}_{0}I}_{t}^{\alpha _{2}} \bigl(e^{ \lambda t}\widetilde{g}(t) \bigr)+d_{1}t^{\alpha _{2}-1}+d_{2}t^{\alpha _{2}-2} \\ &= \int _{0}^{t}\frac{(t-s)^{\alpha _{1}-1}}{\varGamma (\alpha )}e^{ \lambda s} \widetilde{g}(s)\,ds+d_{1}t^{\alpha _{2}-1}+d_{2}t^{\alpha _{2}-2}. \end{aligned} $$

Since \(\varphi _{p} ({{}_{0}^{R}\mathbb{D}}_{t}^{\alpha _{1},\lambda } u(0) )=0\), we see that \(d_{2}=0\), that is,

$$ \varphi _{p} \bigl({{}_{0}^{R} \mathbb{D}}_{t}^{\alpha _{1},\lambda } u(t) \bigr)=e^{-\lambda t}{{}_{0}I}_{t}^{\alpha _{2}} \bigl(e^{\lambda t} \widetilde{g}(t) \bigr)+d_{1}e^{-\lambda t}t^{\alpha _{2}-1}. $$
(2.11)

Furthermore, applying the tempered fractional derivative operators \({{}_{0}^{R}\mathbb{D}}_{t}^{\gamma _{i},\lambda }(i=1,2)\) on both sides of Eq. (2.11), we have

$$ \begin{aligned}[b] &{{}_{0}^{R} \mathbb{D}}_{t}^{\gamma _{i},\lambda } \bigl( \varphi _{p} \bigl({{}_{0}^{R}\mathbb{D}}_{t}^{\alpha _{1},\lambda } u(t) \bigr) \bigr) \\ &\quad ={{}_{0}^{R}\mathbb{D}}_{t}^{\gamma _{i},\lambda } \bigl(e^{- \lambda t}{{}_{0}I}_{t}^{\alpha _{2}} \bigl(e^{\lambda t}\widetilde{g}(t) \bigr) \bigr)+d_{1}{{}_{0}^{R} \mathbb{D}}_{t}^{\gamma _{i},\lambda }\bigl(e^{- \lambda t}t^{\alpha _{2}-1} \bigr) \\ &\quad =e^{-\lambda t}{{}_{0}I}_{t}^{\alpha _{2}-\gamma _{i}} \bigl(e^{\lambda t} \widetilde{g}(t) \bigr)+d_{1}e^{-\lambda t}{{}_{0}^{R}D}_{t}^{\gamma _{i}} \bigl(t^{ \alpha _{2}-1}\bigr) \\ &\quad = \int _{0}^{t} \frac{(t-s)^{\alpha _{2}-\gamma _{i}-1}e^{\lambda (s-t)}}{\varGamma (\alpha _{2}-\gamma _{i})} \widetilde{g}(s)\,ds+d_{1} \frac{\varGamma (\alpha _{2})}{\varGamma (\alpha _{2}-\gamma _{i})}e^{- \lambda t}t^{\alpha _{2}-1-\gamma _{i}}. \end{aligned} $$
(2.12)

From (2.12), we have

$$ \textstyle\begin{cases} {{}_{0}^{R}\mathbb{D}}_{t}^{\gamma _{1},\lambda } [\varphi _{p} ({{}_{0}^{R}\mathbb{D}}_{t}^{\alpha _{1},\lambda } u ) ](1)= \int _{0}^{1} \frac{(1-s)^{\alpha _{2}-\gamma _{1}-1}e^{\lambda (s-1)}}{\varGamma (\alpha _{2}-\gamma _{1})} \tilde{g}(s)\,ds+d_{1} \frac{\varGamma (\alpha _{2})e^{-\lambda }}{\varGamma (\alpha _{2}-\gamma _{1})}, \\ {{}_{0}^{R}\mathbb{D}}_{t}^{\gamma _{2},\lambda } [\varphi _{p} ({{}_{0}^{R}\mathbb{D}}_{t}^{\alpha _{1},\lambda } u ) ](t)= \int _{0}^{t} \frac{(t-s)^{\alpha _{2}-\gamma _{2}-1}e^{\lambda (s-t)}}{\varGamma (\alpha _{2}-\gamma _{2})} \tilde{g}(s)\,ds+d_{1} \frac{\varGamma (\alpha _{2})e^{-\lambda t}}{\varGamma (\alpha _{2}-\gamma _{2})}t^{ \alpha _{2}-1-\gamma _{2}}. \end{cases} $$
(2.13)

Substituting (2.13) into the integral boundary value condition \({{}_{0}^{R}\mathbb{D}}_{t}^{\gamma _{1},\lambda } (\varphi _{p} ({{}_{0}^{R}\mathbb{D}}_{t}^{\alpha _{1},\lambda } u ) )(1)= \int _{0}^{\eta }a(s){{}_{0}^{R}\mathbb{D}}_{t}^{\gamma _{2},\lambda } [\varphi _{p} ({{}_{0}^{R}\mathbb{D}}_{t}^{\alpha _{1},\lambda } u(s) ) ]\,dA(s)\), we obtain

$$ \begin{aligned}[b] d_{1}&= \frac{-1}{\varGamma (\alpha _{2})\Delta } \biggl\{ \int _{0}^{1} \frac{(1-s)^{\alpha _{2}-\gamma _{1}-1}e^{\lambda (s-1)}}{\varGamma (\alpha _{2}-\gamma _{1})} \widetilde{g}(s)\,ds \\ &\quad {}- \int _{0}^{\eta }a(t)\,dA(t) \int _{0}^{t} \frac{(t-s)^{\alpha _{2}-\gamma _{2}-1}e^{\lambda (s-t)}}{\varGamma (\alpha _{2}-\gamma _{2})} \widetilde{g}(s)\,ds \biggr\} . \end{aligned} $$
(2.14)

Substituting (2.14) into (2.11), we get

$$ \begin{aligned} &\varphi _{p} \bigl({{}_{0}^{R} \mathbb{D}}_{t}^{\alpha _{1}, \lambda } u(t) \bigr) \\ &\quad =e^{-\lambda t} \int _{0}^{t} \frac{(t-s)^{\alpha _{2}-1}e^{\lambda s}}{\varGamma (\alpha _{2})} \widetilde{g}(s)\,ds- \frac{e^{-\lambda t}t^{\alpha _{2}-1}}{\varGamma (\alpha _{2})\Delta } \biggl\{ \int _{0}^{1} \frac{(1-s)^{\alpha _{2}-\gamma _{1}-1}e^{-\lambda }}{\varGamma (\alpha _{2}-\gamma _{1})}e^{ \lambda s} \widetilde{g}(s)\,ds \\ &\qquad {}- \int _{0}^{\eta }a(t)\,dA(t) \int _{0}^{t} \frac{(t-s)^{\alpha _{2}-\gamma _{2}-1}e^{\lambda (s-t)}}{\varGamma (\alpha _{2}-\gamma _{2})} \widetilde{g}(s)\,ds \biggr\} \\ &\quad = \int _{0}^{t} \frac{(t-s)^{\alpha _{2}-1}e^{\lambda (s-t)}}{\varGamma (\alpha _{2})} \widetilde{g}(s)\,ds - \frac{e^{-\lambda t}t^{\alpha _{2}-1}}{\varGamma (\alpha _{2})} \int _{0}^{1} (1-s)^{\alpha _{2}-\gamma _{1}-1}e^{\lambda s} \widetilde{g}(s)\,ds \\ &\qquad {}- \frac{e^{-\lambda t}t^{\alpha _{2}-1}\delta }{\varGamma (\alpha _{2})\varGamma (\alpha _{2}-\gamma _{2})\Delta } \int _{0}^{1} (1-s)^{\alpha _{2}-\gamma _{1}-1}e^{\lambda s} \widetilde{g}(s)\,ds \\ &\qquad {}+ \frac{e^{-\lambda t}t^{\alpha _{2}-1}}{\varGamma (\alpha _{2})\varGamma (\alpha _{2}-\gamma _{2})\Delta } \int _{0}^{\eta }a(t)\,dA(t) \int _{0}^{t} (t-s)^{\alpha _{2}-\gamma _{2}-1}e^{ \lambda (s-t)} \widetilde{g}(s)\,ds \\ &\quad = \int _{0}^{t} \frac{(t-s)^{\alpha _{2}-1}e^{-\lambda t}e^{\lambda s}}{\varGamma (\alpha _{2})} \widetilde{g}(s)\,ds- \int _{0}^{1} \frac{(1-s)^{\alpha _{2}-\gamma _{1}-1}t^{\alpha _{2}-1}e^{\lambda (s-t)}}{\varGamma (\alpha _{2})} \widetilde{g}(s)\,ds \\ &\qquad {}- \frac{t^{\alpha _{2}-1}e^{-\lambda t}}{\varGamma (\alpha _{2}-\gamma _{2})\Delta } \int _{0}^{\eta }a(t)\,dA(t) \int _{0}^{1} \frac{(1-s)^{\alpha _{2}-\gamma _{1}-1}t^{\alpha _{2}-\gamma _{2}-1}e^{\lambda (s-t)}}{\varGamma (\alpha _{2})} \widetilde{g}(s)\,ds \\ &\qquad {}+ \frac{t^{\alpha _{2}-1}e^{-\lambda t}}{\varGamma (\alpha _{2})\varGamma (\alpha _{2}-\gamma _{2})\Delta } \int _{0}^{\eta }a(t)\,dA(t) \int _{0}^{t} \frac{(t-s)^{\alpha _{2}-\gamma _{2}-1}e^{\lambda (s-t)}}{\varGamma (\alpha _{2})} \widetilde{g}(s)\,ds \\ &\quad =- \int _{0}^{1}G_{1}(t,s) \widetilde{g}(s)\,ds- \frac{t^{\alpha _{2}-1}e^{-\lambda t}}{\varGamma (\alpha _{2}-\gamma _{2})\Delta } \int _{0}^{1}\widetilde{g}(s)\,ds \int _{0}^{\eta }G_{2}(t,s)a(t)\,dA(t) \\ &\quad =- \int _{0}^{1}G(t,s)\widetilde{g}(s)\,ds. \end{aligned} $$

By employing the p-Laplacian operator \(\varphi _{q}\) on both sides of the above equation, we have

$$ {{}_{0}^{R}\mathbb{D}}_{t}^{\alpha _{1},\lambda } u(t)+\varphi _{q} \biggl( \int _{0}^{1}G(t,s)\widetilde{g}(s)\,ds \biggr)=0. $$
(2.15)

Setting \({g}(t):\triangleq \varphi _{q} (\int _{0}^{1}G(t,s)\widetilde{g}(s)\,ds )\), thus, the p-Laplacian tempered fractional differential system (2.8) is equivalent to the integral boundary value problems as follows:

$$ \textstyle\begin{cases} {{}_{0}^{R}}\mathbb{D}_{t}^{\alpha _{1},\lambda }u(t)+g(t)=0,\quad n-1< \alpha _{1}\leq n; \\ u(0)=u'(0)=\cdots =u^{(n-2)}(0)=0; \\ u(1)=\beta \int _{0}^{1}e^{-\lambda (1-t)}u(t)\,dt, \quad \beta < \alpha _{1}. \end{cases} $$
(2.16)

By means of Lemma 2.3, we see that the integral boundary value problem (2.16) has a unique integral solution

$$ \begin{aligned} u(t)&= \int _{0}^{1}H(t,s)g(s)\,ds \\ &= \int _{0}^{1}H(t,s)\varphi _{q} \biggl( \int _{0}^{1}G(s,\tau ) \widetilde{g}(\tau )\,d\tau \biggr)\,ds, \end{aligned} $$

where the Green function \(G(t,s)\) and \(H(t,s)\) are given by (2.10) and (2.3), respectively. This constitutes the complete proof. □

Lemma 2.5

For\(\forall (s,t)\in [0,1]\times [0,1]\), the Green function\(H(t,s)\)given by (2.3) has the following properties:

\((A_{1})\):

\(H(t,s)\)is continuous and\(H(t,s)\geq 0\);

\((A_{2})\):

\(m_{1}(s)e^{-\lambda t}t^{\alpha _{1}-1}\leq H(t,s)\leq M_{1}(s)e^{- \lambda t}t^{\alpha _{1}-1}\), where

$$ M_{1}(s)= \frac{\alpha _{1}(1-s)^{\alpha _{1}-1}(\alpha _{1}-\beta +\beta s)e^{\lambda s}}{(\alpha _{1}-\beta )\varGamma (\alpha _{1}+1)},\qquad m_{1}(s)= \frac{\alpha _{1}\beta s(1-s)^{\alpha _{1}-1} e^{\lambda s} }{(\alpha _{1}-\beta )\varGamma (\alpha _{1}+1)}. $$

Proof

Evidently, \(H(t,s)\) is continuous and \(H(t,s)\leq M_{1}(s)e^{-\lambda t}t^{\alpha _{1}-1}\) holds. So, we only need to prove the inequality \(H(t,s)\geq m_{1}(s)e^{-\lambda t}t^{\alpha _{1}-1}\) and \(H(t,s)\geq 0\).

If \(0\leq s\leq t\leq 1\), then we have \(0\leq t-s\leq t-ts =t(1-s)\), and thus \((t-s)^{\alpha _{1}-1}\leq t^{\alpha _{1}-1}(1-s)^{\alpha _{1}-1}\). Hence, we get

$$ \begin{aligned} H(t,s)&= \frac{[\alpha _{1}^{2}(1-s)^{\alpha _{1}-1}-\alpha _{1}\beta (1-s)^{\alpha _{1}}]t^{\alpha _{1}-1}-\alpha _{1}(\alpha _{1}-\beta )(t-s)^{\alpha _{1}-1}}{(\alpha _{1}-\beta )\varGamma (\alpha _{1}+1)}e^{- \lambda t}e^{\lambda s} \\ &\geq \frac{[\alpha _{1}^{2}(1-s)^{\alpha _{1}-1}-\alpha _{1}\beta (1-s)^{\alpha _{1}}]t^{\alpha _{1}-1}-\alpha _{1}(\alpha _{1}-\beta )t^{\alpha _{1}-1}(1-s)^{\alpha _{1}-1}}{(\alpha _{1}-\beta )\varGamma (\alpha _{1}+1)}e^{- \lambda t}e^{\lambda s} \\ &= \frac{\alpha _{1}\beta s(1-s)^{\alpha _{1}-1} e^{\lambda s} }{(\alpha _{1}-\beta )\varGamma (\alpha _{1}+1)}e^{- \lambda _{1} t}t^{\alpha _{1}-1}\geq 0. \end{aligned} $$

If \(0\leq s\leq t\leq 1\), clearly, we can see that

$$ \begin{aligned} H(t,s)&= \frac{[\alpha _{1}^{2}(1-s)^{\alpha _{1}-1}-\alpha _{1}\beta (1-s)^{\alpha _{1}}]t^{\alpha _{1}-1}}{(\alpha _{1}-\beta )\varGamma (\alpha _{1}+1)}e^{- \lambda t}e^{\lambda s} \\ &\geq \frac{[\alpha _{1}^{2}(1-s)^{\alpha _{1}-1}-\alpha _{1}\beta (1-s)^{\alpha _{1}}]t^{\alpha _{1}-1}-\alpha _{1}(\alpha _{1}-\beta )t^{\alpha _{1}-1}(1-s)^{\alpha _{1}-1}}{(\alpha _{1}-\beta )\varGamma (\alpha _{1}+1)}e^{- \lambda t}e^{\lambda s} \\ &\geq \frac{\alpha _{1}\beta s(1-s)^{\alpha _{1}-1} e^{\lambda s} }{(\alpha _{1}-\beta )\varGamma (\alpha _{1}+1)}e^{- \lambda t}t^{\alpha _{1}-1}\geq 0. \end{aligned} $$

Hence, the proof is complete. □

Lemma 2.6

Suppose that

$$ (H)\quad e^{-\lambda }\varGamma (\alpha _{2}-\gamma _{2})>\varGamma (\alpha _{2}- \gamma _{1}) \int _{0}^{\eta }e^{-\lambda s}s^{\alpha _{2}-\gamma _{2}-1}a(s)\,dA(s), $$

then, for all\((t,s)\in [0,1]\times [0,1]\), the Green function\(G(t,s)\)is continuous and satisfies:

\((B_{1})\):

\(G_{1}(t,s)\geq 0\), \(G_{2}(t,s)\geq 0\), and\(G(t,s)\geq 0\);

\((B_{2})\):

\(\frac{e^{\lambda s}[(1-s)^{\alpha _{2}-\gamma _{1}-1}-(1-s)^{\alpha _{2}-1}]}{\varGamma (\alpha _{2})}e^{- \lambda t}t^{\alpha _{2}-1}\leq G_{1}(t,s)\leq \frac{e^{\lambda s}(1-s)^{\alpha _{2}-\gamma _{1}-1}}{\varGamma (\alpha _{2})}e^{- \lambda t}t^{\alpha _{2}-1}\);

\((B_{3})\):

\(\frac{e^{\lambda s}[(1-s)^{\alpha _{2}-\gamma _{1}-1}-(1-s)^{\alpha _{2}-\gamma _{2}-1}]}{\varGamma (\alpha _{2})}e^{- \lambda t}t^{\alpha _{2}-\gamma _{2}-1}\leq G_{2}(t,s)\leq \frac{e^{\lambda s}(1-s)^{\alpha _{2}-\gamma _{1}-1}}{\varGamma (\alpha _{2})}e^{- \lambda t}t^{\alpha _{2}-\gamma _{2}-1}\);

\((B_{4})\):

\(m_{2}(s)e^{-\lambda t}t^{\alpha _{2}-1}\leq G(t,s)\leq M_{2}(s)e^{- \lambda t}t^{\alpha _{2}-1}\), where

$$ \textstyle\begin{cases} M_{2}(s)=[\frac{1}{\varGamma (\alpha _{2})}+ \frac{\delta }{\Delta \varGamma (\alpha _{2})\varGamma (\alpha _{2}-\gamma _{2})} ]e^{\lambda s}(1-s)^{\alpha _{2}-\gamma _{1}-1}, \\ m_{2}(s)= \frac{e^{\lambda s}[(1-s)^{\alpha _{2}-\gamma _{1}-1}-(1-s)^{\alpha _{2}-1}]}{\varGamma (\alpha _{2})} + \frac{\delta e^{\lambda s}[(1-s)^{\alpha _{2}-\gamma _{1}-1}-(1-s)^{\alpha _{2}-\gamma _{2}-1}]}{\Delta \varGamma (\alpha _{2})\varGamma (\alpha _{2}-\gamma _{2})}. \end{cases} $$

Proof

Firstly, for \((t,s)\in [0,1]\times [0,1]\), it is evident that \(G(t,s)\) and \(G_{i}(t,s)(i=1,2)\) are continuous.

Secondly, for \((B_{2})\) and \((B_{3})\), it is easy to see that the right sides of the inequalities hold, so we only need to prove the left sides of the inequalities. If \(0\leq s\leq t\leq 1\), we have \(0\leq t-s\leq t-ts=(1-s)t\), and thus \((t-s)^{\alpha _{2}-1}\leq (1-s)^{\alpha _{2}-1}t^{\alpha _{2}-1}\). Hence, we have

$$ \begin{aligned} G_{1}(t,s)&= \frac{e^{\lambda (s-t)}}{\varGamma (\alpha _{2})} \bigl[t^{\alpha _{2}-1}(1-s)^{ \alpha _{2}-\gamma _{1}-1}-(t-s)^{\alpha _{2}-1}\bigr] \\ &\geq \frac{e^{\lambda (s-t)}}{\varGamma (\alpha _{2})}\bigl[t^{\alpha _{2}-1}(1-s)^{ \alpha _{2}-\gamma _{1}-1}-(1-s)^{\alpha _{2}-1}t^{\alpha _{2}-1} \bigr] \\ &= \frac{e^{\lambda s}[(1-s)^{\alpha _{2}-\gamma _{1}-1}-(1-s)^{\alpha _{2}-1}]}{\varGamma (\alpha _{2})}e^{- \lambda t}t^{\alpha _{2}-1}. \end{aligned} $$

If \(0\leq t\leq s\leq 1\),

$$ \begin{aligned} G_{1}(t,s)&= \frac{e^{\lambda (s-t)}}{\varGamma (\alpha _{2})}t^{\alpha _{2}-1}(1-s)^{ \alpha _{2}-\gamma _{1}-1} \\ &\geq \frac{e^{\lambda s}[(1-s)^{\alpha _{2}-\gamma _{1}-1}-(1-s)^{\alpha _{2}-1}]}{\varGamma (\alpha _{2})}e^{- \lambda t}t^{\alpha _{2}-1}. \end{aligned} $$

Furthermore, from \((1-s)^{\alpha _{2}-\gamma _{1}-1}>(1-s)^{\alpha _{2}-1}\), we get \(G_{1}(t,s)\geq 0\) for \(\forall (t,s)\in [0,1]\times [0,1]\). In the same way, similar conclusions can be obtained for \(G_{2}(t,s)\).

Finally, from \((B_{2})\) and \((B_{3})\), we can know that \(m_{2}(s)e^{-\lambda t}t^{\alpha -1}\leq G(t,s)\leq M_{2}(s)e^{- \lambda t}t^{\alpha -1}\). Since the condition \((H)\) holds, it is easy to see that \(\Delta >0\). Combining \((1-s)^{\alpha _{2}-\gamma _{1}-1}>(1-s)^{\alpha _{2}-1}\) with \(\Delta >0\), we obtain \(m_{2}(s)\geq 0\). Then \(G(s,t)\geq 0\) for \(\forall (t,s)\in [0,1]\times [0,1]\). Therefore, our justification for the proof is complete. □

Lemma 2.7

([20])

Let\(\xi \in (0,1)\), \(A:P\times P\rightarrow P\)be a mixed monotone operator that satisfies

$$ A\bigl(tx,t^{-1}y\bigr)\geq t^{\xi }A(x,y), \quad t\in (0,1), x,y\in P. $$
(2.17)

\(B:P\rightarrow P\)is an increasing sub-homogeneous operator. Assume that

  1. (I)

    there is\(h_{0}\in P_{h}\)such that\(A(h_{0},h_{0})\in P_{h}\)and\(Bh_{0}\in P_{h}\);

  2. (II)

    there exists a constant\(\delta _{0}>0\)such that\(A(x,y)\geq \delta _{0} Bx\), \(\forall x,y\in P\).

Then:

  1. (1)

    \(A:P_{h}\times P_{h}\rightarrow P_{h}\), \(B:P_{h}\rightarrow P_{h}\);

  2. (2)

    there exist\(u_{0},v_{0}\in P_{h}\)and\(r\in (0,1)\)such that

    $$ rv_{0}\leq u_{0}< v_{0},\qquad u_{0}\leq A(u_{0},v_{0})+Bu_{0} \leq A(v_{0},u_{0})+B(v_{0}) \leq v_{0}; $$
  3. (3)

    the operator equation\(A(x,x)+Bx=x\)has a unique solution\(x^{*}\)in\(P_{h}\);

  4. (4)

    for any initial values\(x_{0},y_{0}\in P_{h}\), constructing successively the sequences

    $$ x_{n}=A(x_{n-1},y_{n-1})+Bx_{n-1},\qquad y_{n}=A(y_{n-1},x_{n-1})+By_{n-1},\quad n=1,2, \ldots , $$

    we have\(x_{n}\rightarrow x^{*}\)and\(y_{n}\rightarrow x^{*}\)as\(n\rightarrow \infty \).

3 Main results

In this section, we will work in the Banach space \(C[0,1]\), the space of all continuous functions on \([0,1]\). It is obvious that this space can be equipped with a partial order

$$ x,y\in C[0,1],\quad x\leq y\quad \Longleftrightarrow\quad x(t)\leq y(t)\quad \mbox{for } t\in [0,1]. $$

Setting \(P=\{x\in C[0,1]\mid x(t)\geq 0, t\in [0,1]\}\) and \(h(t)=e^{-\lambda t}t^{\alpha _{1}-1}\), then we see that P is a normal cone in \(C[0,1]\).

From Lemma 2.4, we can recognize that the p-Laplacian differential equation integral boundary value problem (1.1) is equivalent to the integral formulation given by

$$ u(t)= \int _{0}^{1}H(t,s)\varphi _{q} \biggl( \int _{0}^{1}G(s,\tau ) \bigl[f \bigl(\tau ,u( \tau ),u(\tau ) \bigr)+g \bigl(\tau ,u(\tau ) \bigr) \bigr]\,d\tau \biggr) \,ds. $$

For convenience, we define an operator T by

$$ T(u,v) (t)= \int _{0}^{1}H(t,s)\varphi _{q} \biggl( \int _{0}^{1}G(s,\tau ) \bigl[f \bigl(\tau ,u( \tau ),v(\tau ) \bigr)+g \bigl(\tau ,u(\tau ) \bigr) \bigr]\,d\tau \biggr) \,ds. $$
(3.1)

It is evident that \(u^{*}\) is a solution of p-Laplacian differential equation integral boundary value problem (1.1) if and only if \(T(u^{*},u^{*})=u^{*}\).

Theorem 3.1

Assume that the condition\((H)\)holds, and the following conditions are satisfied:

\((H_{1})\):

\(f(t,u,v):[0,1]\times [0,+\infty )\times [0,+\infty )\rightarrow [0,+ \infty )\)is continuous, \(g(t,u):[0,1]\times [0,+\infty )\rightarrow [0,+\infty )\)is continuous with\(g(t,u)\not \equiv 0\)and\(a(t):[0,1]\rightarrow R^{+}\)is continuous; for fixed\(t\in [0,1]\), \(f(t,u,v)\)is increasing in\(u\in [0,+\infty )\)and decreasing in\(v\in [0,+\infty )\), \(g(t,u)\)is increasing in\(u\in [0,+\infty )\).

\((H_{2})\):

For\(\forall t\in [0,1]\), \(\gamma \in (0,1)\), \(u,v\in [0,+\infty )\), there exists a constant\(\xi \in (0,1)\)such that

$$\begin{aligned}& f\bigl(t,\gamma u,\gamma ^{-1} v\bigr)\geq \varphi _{p}^{\xi }(\gamma ) f(t,u,v), \end{aligned}$$
(3.2)
$$\begin{aligned}& g(t,\gamma u)\geq \varphi _{p}(\gamma ) g(t,u). \end{aligned}$$
(3.3)
\((H_{3})\):

For\(\forall t\in [0,1]\)and\(u,v\in [0,+\infty )\), there exists a constant\(\delta _{0}>0\)such that

$$ f(t,u,v)\geq \varphi _{p} (\delta _{0}) g(t,u). $$
(3.4)

Then we have:

  1. (I)

    thep-Laplacian differential equation integral boundary value problem involving tempered fractional derivative (1.1) has a unique positive solution\(u^{*}\in P_{h}\), where\(h(t)=e^{-\lambda t}t^{\alpha _{1}-1}\), \(t\in [0,1]\);

  2. (II)

    for\(\forall t\in [0,1]\), there exist\(u_{0},v_{0}\in P_{h}\)and\(r\in (0,1)\)such that\(rv_{0}\leq u_{0}< v_{0}\)and

    $$\begin{aligned}& u_{0}(t)\leq \int _{0}^{1}H(t,s)\varphi _{q} \biggl( \int _{0}^{1}G(s, \tau ) \bigl[f \bigl(\tau ,u_{0}(\tau ),v_{0}(\tau ) \bigr)+g \bigl(\tau ,u_{0}( \tau ) \bigr) \bigr]\,d\tau \biggr) \,ds, \\& v_{0}(t)\geq \int _{0}^{1}H(t,s)\varphi _{q} \biggl( \int _{0}^{1}G(s, \tau ) \bigl[f \bigl(\tau ,v_{0}(\tau ),u_{0}(\tau ) \bigr)+g \bigl(\tau ,v_{0}( \tau ) \bigr) \bigr]\,d\tau \biggr) \,ds; \end{aligned}$$
  3. (III)

    for any initial values\(x_{0},y_{0}\in P_{h}\), making successively the sequences

    $$\begin{aligned}& x_{n}= \int _{0}^{1}H(t,s)\varphi _{q} \biggl( \int _{0}^{1}G(s,\tau ) \bigl[f \bigl(\tau ,x_{n-1}(\tau ),y_{n-1}(\tau ) \bigr)+g \bigl(\tau ,x_{n-1}( \tau ) \bigr) \bigr]\,d\tau \biggr) \,ds, \\& y_{n}= \int _{0}^{1}H(t,s)\varphi _{q} \biggl( \int _{0}^{1}G(s,\tau ) \bigl[f \bigl(\tau ,y_{n-1}(\tau ),x_{n-1}(\tau ) \bigr)+g \bigl(\tau ,y_{n-1}( \tau ) \bigr) \bigr]\,d\tau \biggr) \,ds, \\& \quad n=0,1,2,\ldots, \end{aligned}$$

    we obtain\(x_{n} \rightarrow u^{*}\)and\(y_{n}\rightarrow u^{*}\)as\(n\rightarrow \infty \).

Proof

Firstly, we define two operators \(A:P\times P\rightarrow E\) and \(B:P\rightarrow E\) by

$$\begin{aligned}& A(u,v) (t)= \int _{0}^{1}H(t,s)\varphi _{q} \biggl( \int _{0}^{1}G(s,\tau )f \bigl(\tau ,u(\tau ),v( \tau ) \bigr)\,d\tau \biggr) \,ds, \end{aligned}$$
(3.5)
$$\begin{aligned}& B(u) (t)= \int _{0}^{1}H(t,s)\varphi _{q} \biggl( \int _{0}^{1}G(s,\tau )g \bigl(\tau ,u(\tau ) \bigr)\,d\tau \biggr) \,ds. \end{aligned}$$
(3.6)

From (3.1), we have \(T(u,v)=A(u,v)+B(u)\) and u is a solution of the p-Laplacian differential system (1.1) if and only if \(T(u,u)=u\). We show that the operator A satisfies the condition (2.17) in Lemma 2.7 and the operator B is a sub-homogeneous operator.

From \((H_{1})\), Lemma 2.5 and Lemma 2.6, we know that \(A:P\times P\rightarrow P\) and \(B:P\rightarrow P\). In addition, it follows from \((H_{1})\) and \((H_{2})\) that A is a mixed monotone operator and B is an increasing operator. For \(\forall \gamma \in (0,1)\) and \(u,v\in P\), from (3.2), we obtain

$$ \begin{aligned}[b] A\bigl(\gamma u,\gamma ^{-1}v\bigr) (t)&= \int _{0}^{1}H(t,s) \varphi _{q} \biggl( \int _{0}^{1}G(s,\tau )f \bigl(\tau ,\gamma u( \tau ), \gamma v(\tau ) \bigr)\,d\tau \biggr) \,ds \\ &\geq \int _{0}^{1}H(t,s)\varphi _{q} \biggl(\varphi _{p}^{\xi }(\gamma ) \int _{0}^{1}G(s,\tau )f \bigl(\tau ,u(\tau ),v( \tau ) \bigr)\,d\tau \biggr) \,ds \\ &=\gamma ^{\xi }A(u,v) (t). \end{aligned} $$
(3.7)

That is, \(A(\gamma u,\gamma ^{-1}v)\geq \gamma ^{\xi }A(u,v)\) for \(\forall \gamma \in (0,1)\), \(u,v\in P\). Furthermore, for \(\forall \gamma \in (0,1)\) and \(u\in P\), from (3.3), we have

$$ \begin{aligned}[b] B(\gamma u) (t)&= \int _{0}^{1}H(t,s)\varphi _{q} \biggl( \int _{0}^{1}G(s,\tau )g \bigl(\tau ,\gamma u( \tau ) \bigr)\,d\tau \biggr) \,ds \\ &\geq \varphi _{q} \bigl(\varphi _{p}(\gamma ) \bigr) \int _{0}^{1}H(t,s) \varphi _{q} \biggl( \int _{0}^{1}G(s,\tau )g \bigl(\tau ,u(\tau ) \bigr)\,d\tau \biggr) \,ds \\ &=\gamma B(u) (t). \end{aligned} $$
(3.8)

That is, the operator B is a sub-homogeneous operator.

Secondly, we show that \(A(h,h)\in P_{h}\) and \(Bh\in P_{h}\). From Lemma 2.5 and Lemma 2.6, we have

$$\begin{aligned} A(h,h) (t)&= \int _{0}^{1}H(t,s)\varphi _{q} \biggl( \int _{0}^{1}G(s, \tau )f \bigl(\tau ,h(\tau ),h(\tau ) \bigr)\,d\tau \biggr) \,ds \\ &\leq \int _{0}^{1}H(t,s)\varphi _{q} \biggl( \int _{0}^{1}M_{2}(\tau )e^{- \lambda s}s^{\alpha _{2}-1}f \bigl(\tau ,h(\tau ),h(\tau ) \bigr)\,d\tau \biggr) \,ds \\ &\leq \int _{0}^{1}M_{1}(s)e^{-\lambda t}t^{\alpha _{1}-1} \varphi _{q} \biggl( \int _{0}^{1}M_{2}(\tau )e^{-\lambda s}s^{\alpha _{2}-1}f \bigl( \tau ,h(\tau ),h(\tau ) \bigr)\,d\tau \biggr) \,ds \\ &\leq \biggl\{ \int _{0}^{1} \frac{M_{1}(s)s^{(\alpha _{2}-1)(q-1)}}{e^{\lambda s(q-1)}}\varphi _{q} \biggl( \int _{0}^{1}M_{2}(\tau )f (\tau ,h_{\max },0 )\,d\tau \biggr) \,ds \biggr\} e^{-\lambda t}t^{\alpha _{1}-1} \end{aligned}$$

and

$$ \begin{aligned} A(h,h) (t)&= \int _{0}^{1}H(t,s)\varphi _{q} \biggl( \int _{0}^{1}G(s, \tau )f \bigl(\tau ,h(\tau ),h(\tau ) \bigr)\,d\tau \biggr) \,ds \\ &\geq \int _{0}^{1}H(t,s)\varphi _{q} \biggl( \int _{0}^{1}m_{2}(\tau )e^{- \lambda s}s^{\alpha _{2}-1}f \bigl(\tau ,h(\tau ),h(\tau ) \bigr)\,d\tau \biggr) \,ds \\ &\geq \int _{0}^{1}m_{1}(s)e^{-\lambda t}t^{\alpha _{1}-1} \varphi _{q} \biggl( \int _{0}^{1}m_{2}(\tau )e^{-\lambda s}s^{\alpha _{2}-1}f \bigl( \tau ,h(\tau ),h(\tau ) \bigr)\,d\tau \biggr) \,ds \\ &\geq \biggl\{ \int _{0}^{1} \frac{m_{1}(s)s^{(\alpha _{2}-1)(q-1)}}{e^{\lambda s(q-1)}}\varphi _{q} \biggl( \int _{0}^{1}m_{2}(\tau )f (\tau ,0,h_{\max } )\,d\tau \biggr) \,ds \biggr\} e^{-\lambda t}t^{\alpha _{1}-1}, \end{aligned} $$

where \(h_{\max }=\max \{h(t):t\in [0,1]\}\). Setting

$$\begin{aligned}& L_{1}= \int _{0}^{1} \frac{M_{1}(s)s^{(\alpha _{2}-1)(q-1)}}{e^{\lambda s(q-1)}}\varphi _{q} \biggl( \int _{0}^{1}M_{2}(\tau )f (\tau ,h_{\max },0 )\,d\tau \biggr) \,ds, \\& l_{1}= \int _{0}^{1} \frac{m_{1}(s)s^{(\alpha _{2}-1)(q-1)}}{e^{\lambda s(q-1)}}\varphi _{q} \biggl( \int _{0}^{1}m_{2}(\tau )f (\tau ,0,h_{\max } )\,d\tau \biggr) \,ds, \end{aligned}$$

it is easy to see that \(L_{1}>l_{1}>0\). Hence, we get \(l_{1}h(t)\leq A(h,h)\leq L_{1}h(t)\). That is, \(A(h,h)\in P_{h}\). Similarly,

$$ \begin{aligned} B(h) (t)&= \int _{0}^{1}H(t,s)\varphi _{q} \biggl( \int _{0}^{1}G(s, \tau )g \bigl(\tau ,h(\tau ) \bigr)\,d\tau \biggr) \,ds \\ &\leq \biggl\{ \int _{0}^{1} \frac{M_{1}(s)s^{(\alpha _{2}-1)(q-1)}}{e^{\lambda s(q-1)}}\varphi _{q} \biggl( \int _{0}^{1}M_{2}(\tau )g (\tau ,h_{\max } )\,d\tau \biggr) \,ds \biggr\} e^{-\lambda t}t^{\alpha _{1}-1} \end{aligned} $$

and

$$ \begin{aligned} B(h) (t)&= \int _{0}^{1}H(t,s)\varphi _{q} \biggl( \int _{0}^{1}G(s, \tau )g \bigl(\tau ,h(\tau ) \bigr)\,d\tau \biggr) \,ds \\ &\geq \biggl\{ \int _{0}^{1} \frac{m_{1}(s)s^{(\alpha _{2}-1)(q-1)}}{e^{\lambda s(q-1)}}\varphi _{q} \biggl( \int _{0}^{1}m_{2}(\tau )g (\tau ,0 )\,d\tau \biggr) \,ds \biggr\} e^{-\lambda t}t^{\alpha _{1}-1}. \end{aligned} $$

Set

$$\begin{aligned}& L_{2}= \int _{0}^{1} \frac{M_{1}(s)s^{(\alpha _{2}-1)(q-1)}}{e^{\lambda s(q-1)}}\varphi _{q} \biggl( \int _{0}^{1}M_{2}(\tau )g (\tau ,h_{\max } )\,d\tau \biggr) \,ds, \\& l_{2}= \int _{0}^{1} \frac{m_{1}(s)s^{(\alpha _{2}-1)(q-1)}}{e^{\lambda s(q-1)}}\varphi _{q} \biggl( \int _{0}^{1}m_{2}(\tau )g (\tau ,0 )\,d\tau \biggr) \,ds. \end{aligned}$$

From \(L_{2}>l_{2}>0\) and \(l_{2}h\leq B(h)\leq L_{2}h\), we get \(Bh\in P_{h}\). Since \(h\in P_{h}\), letting \(h_{0}=h\), we see that the condition \((I_{1})\) of Lemma 2.7 is satisfied.

Finally, we show that the condition (II) of Lemma 2.7 is also satisfied. For \(u,v\in P\), from (3.4), we get

$$ \begin{aligned}[b] A(u,v) (t)&= \int _{0}^{1}H(t,s)\varphi _{q} \biggl( \int _{0}^{1}G(s, \tau )f \bigl(\tau ,u(\tau ),v(\tau ) \bigr)\,d\tau \biggr) \,ds \\ &\geq \int _{0}^{1}H(t,s)\varphi _{q} \biggl( \int _{0}^{1}\varphi _{p}( \delta _{0})G(s,\tau )g \bigl(\tau ,u(\tau ) \bigr)\,d\tau \biggr) \,ds \\ &=\delta _{0}B(u) (t). \end{aligned} $$
(3.9)

Now, all conditions of Lemma 2.7 are satisfied. Hence, the conclusions of Theorem 3.1 follow from Lemma 2.7. □

Corollary 3.1

Assume that the condition\((H)\)holds and

\((H_{1}')\):

\(f(t,u,v):[0,1]\times [0,+\infty )\times [0,+\infty )\rightarrow [0,+ \infty )\)is continuous, \(g(t,u)\equiv 0\)for\(\forall t\in [0,1]\)and\(u\in [0,+\infty )\), \(a(t):[0,1]\rightarrow R^{+}\)is continuous;

\((H_{2}')\):

\(f(t,u,v)\)is increasing in\(u\in [0,+\infty )\)for fixed\(t\in [0,1]\)and\(v\in [0,+\infty )\), decreasing in\(v\in [0,+\infty )\)for fixed\(t\in [0,1]\)and\(u\in [0,+\infty )\);

\((H_{3}')\):

for\(\forall t\in [0,1]\), \(\gamma \in (0,1)\), \(u,v\in [0,+\infty )\), there exists a constant\(\xi \in (0,1)\)such that

$$ f\bigl(t,\gamma u,\gamma ^{-1} v\bigr)\geq \varphi _{p}^{\xi }(\gamma ) f(t,u,v). $$

Then we have:

(I):

Thep-Laplacian differential equation integral boundary value problem

$$ \textstyle\begin{cases} {{}_{0}^{R}\mathbb{D}}_{t}^{\alpha _{2},\lambda } (\varphi _{p} ({{}_{0}^{R} \mathbb{D}}_{t}^{\alpha _{1},\lambda } u(t) ) )=f (t,u(t),u(t) ),\quad t\in [0,1], \\ u(0)=u'(0)=\cdots =u^{(n-2)}(0)=0, \\ \varphi _{p} ({{}_{0}^{R}\mathbb{D}}_{t}^{\alpha _{1},\lambda }u )(0)=0, \\ u(1)=\beta \int _{0}^{1}e^{-\lambda (1-t)}u(t)\,dt, \\ {{}_{0}^{R}\mathbb{D}}_{t}^{\gamma _{1},\lambda } (\varphi _{p} ({{}_{0}^{R}\mathbb{D}}_{t}^{\alpha _{1},\lambda } u ) )(1)= \int _{0}^{\eta }a(s){{}_{0}^{R}\mathbb{D}}_{t}^{\gamma _{2},\lambda } [\varphi _{p} ({{}_{0}^{R}\mathbb{D}}_{t}^{\alpha _{1},\lambda } u(s) ) ]\,dA(s), \end{cases} $$

has a unique positive solution\(u^{*}\in P_{h}\), where\(h(t)=e^{-\lambda t}t^{\alpha _{1}-1}\).

(II):

For\(\forall t\in [0,1]\), there exist\(u_{0},v_{0}\in P_{h}\)and\(r\in (0,1)\)such that\(rv_{0}\leq u_{0}< v_{0}\)and

$$\begin{aligned}& u_{0}(t)\leq \int _{0}^{1}H(t,s)\varphi _{q} \biggl( \int _{0}^{1}G(s, \tau )f \bigl(\tau ,u_{0}(\tau ),v_{0}(\tau ) \bigr)\,d\tau \biggr) \,ds, \\& v_{0}(t)\geq \int _{0}^{1}H(t,s)\varphi _{q} \biggl( \int _{0}^{1}G(s, \tau )f \bigl(\tau ,v_{0}(\tau ),u_{0}(\tau ) \bigr)\,d\tau \biggr) \,ds. \end{aligned}$$
(III):

For any initial values\(x_{0},y_{0}\in P_{h}\), making successively the sequences

$$\begin{aligned}& x_{n}= \int _{0}^{1}H(t,s)\varphi _{q} \biggl( \int _{0}^{1}G(s,\tau )f \bigl(\tau ,x_{n-1}(\tau ),y_{n-1}(\tau ) \bigr)\,d\tau \biggr) \,ds, \\& y_{n}= \int _{0}^{1}H(t,s)\varphi _{q} \biggl( \int _{0}^{1}G(s,\tau )f \bigl(\tau ,y_{n-1}(\tau ),x_{n-1}(\tau ) \bigr)\,d\tau \biggr) \,ds, \\& \quad n=0,1,2,\ldots, \end{aligned}$$

we obtain\(x_{n} \rightarrow u^{*}\)and\(y_{n}\rightarrow u^{*}\)as\(n\rightarrow \infty \).

Proof

Letting \(g(t,u(t))\equiv 0\), from Theorem 3.1, we get the conclusions. □

4 Applications

Example

We consider the p-Laplacian differential equation integral boundary value problem involving a tempered fractional derivative as follows:

$$ \textstyle\begin{cases} {{}_{0}^{R}\mathbb{D}}_{t}^{\frac{3}{2},1} (\varphi _{3} ({{}_{0}^{R} \mathbb{D}}_{t}^{\frac{5}{2},1} u(t) ) )=f (t,u(t),u(t) )+g (t,u(t) ), \quad 0\leq t\leq 1; \\ u(0)=u'(0)=0; \\ \varphi _{3} ({{}_{0}^{R}\mathbb{D}}_{t}^{\frac{5}{2},1} u )(0)=0; \\ u(1)=\int _{0}^{1}e^{t-1}u(t)\,dt, \\ {{}_{0}^{R}\mathbb{D}}_{t}^{\frac{3}{8},1} (\varphi _{3} ({{}_{0}^{R} \mathbb{D}}_{t}^{\frac{5}{2},1} u ) )(1)=\int _{0}^{1} {{}_{0}^{R} \mathbb{D}}_{t}^{\frac{2}{8},1} [\varphi _{3} ({{}_{0}^{R} \mathbb{D}}_{t}^{\frac{5}{2},1} u(t) ) ]\, d(\frac{t}{2}); \end{cases} $$
(4.1)

where \(f(t,u,v)=(1-t)^{-\frac{1}{3}}t^{-\frac{2}{3}}u^{\frac{1}{3}}+v^{- \frac{1}{5}}\), \(g(t,u)=(1-t)^{-\frac{1}{8}}t^{-\frac{1}{6}}u^{\frac{1}{3}}\), \(p=3\), \(\lambda =1>0\), \(\eta =1\) and \(A(t)=\frac{t}{2}\). For any \(t\in (0,1)\), \(u>0\) and \(v>0\) and we see that \(\alpha _{1}=\frac{5}{2}\), \(\alpha _{2}=\frac{3}{2} \gamma _{1}= \frac{3}{8}\), \(\gamma _{2}=\frac{2}{8}\), \(\beta =1\), \(a(t)\equiv 1\) in the systems (4.1).

Let us investigate if all the conditions required in Theorem 3.1 are satisfied.

  1. (1)

    From \(\delta =\int _{0}^{\eta }e^{-\lambda s}s^{\alpha _{2}-\gamma _{2}-1}a(s)\,dA(s)=0.2385\), it is easy to see that \(\varGamma (\alpha _{2}-\gamma _{1})\int _{0}^{\eta }e^{-\lambda s}s^{ \alpha _{2}-\gamma _{2}-1}a(s)\,dA(s)=0.2246\) and \(e^{-\lambda }\varGamma (\alpha _{2}-\gamma _{2})=0.3334\), clearly, \(e^{-\lambda }\varGamma (\alpha _{2}-\gamma _{2})>\varGamma (\alpha _{2}- \gamma _{1})\int _{0}^{\eta }e^{-\lambda s}s^{\alpha _{2}-\gamma _{2}-1}a(s)\,dA(s)\). Then the condition \((H)\) is satisfied.

  2. (2)

    It is obvious that \(f(t,u,v):(0,1)\times R^{+}\times R^{+}\rightarrow R^{+}\) and \(g(t,u):(0,1)\times R^{+}\rightarrow R^{+}\) are continuous. In addition, \(f(t,u,v)\) is increasing in u for fixed \(t\in (0,1)\) and \(v\in R^{+}\), decreasing in v for fixed \(t\in (0,1)\) and \(u\in R^{+}\); furthermore, for fixed \(t\in (0,1)\), \(g(t,u)\) is increasing in u.

  3. (3)

    For any \(\gamma \in (0,1)\), \(t\in (0,1)\), \(u,v>0\), taking \(\xi =\frac{1}{2}\in (0,1)\), we have

    $$ \begin{aligned} f\bigl(t,\gamma u,\gamma ^{-1} v \bigr)&=(1-t)^{-\frac{1}{3}}t^{- \frac{2}{3}}(\gamma u)^{\frac{1}{3}}+\bigl( \gamma ^{-1}v\bigr)^{-\frac{1}{5}} \\ &\geq \gamma ^{\frac{1}{2}}\bigl[(1-t)^{-\frac{1}{3}}t^{-\frac{2}{3}}u^{ \frac{1}{3}}+v^{-\frac{1}{5}} \bigr] \\ &\geq \gamma \bigl[(1-t)^{-\frac{1}{3}}t^{-\frac{2}{3}}u^{\frac{1}{3}}+v^{- \frac{1}{5}} \bigr] \\ &=\varphi _{p}^{\xi }(\gamma )f(t,u,v) \end{aligned} $$

    and

    $$ \begin{aligned} g(t,\gamma u)&=(1-t)^{-\frac{1}{8}}t^{-\frac{1}{6}}( \gamma u)^{\frac{1}{3}} \\ &\geq \gamma ^{2} \bigl[(1-t)^{-\frac{1}{8}}t^{-\frac{1}{6}}u^{\frac{1}{3}} \bigr] \\ &=\varphi _{p}(\gamma ) g(t,u). \end{aligned} $$
  4. (4)

    Taking \(\delta _{0}=\frac{1}{2}\), for \(\forall t\in (0,1)\) and \(u,v\in [0,+\infty )\), we have

    $$ \begin{aligned} f(t,u,v)&=(1-t)^{-\frac{1}{3}}t^{-\frac{2}{3}} u^{ \frac{1}{3}}+v^{-\frac{1}{5}} \\ &\geq \frac{1}{4}\bigl[(1-t)^{-\frac{1}{8}}t^{-\frac{1}{6}}u^{\frac{1}{3}} \bigr] \\ &=\varphi _{p}(\delta _{0})g(t,u). \end{aligned} $$

From the above conditions, we can see that all the assumptions of Theorem 3.1 are satisfied. Hence, Theorem 3.1 implies that the p-Laplacian differential system (4.1) has a unique positive solution \(u^{*}\in P_{h}\), where \(h(t)=e^{-t}t^{\frac{3}{2}}\). Furthermore, for any initial values \(x_{0},y_{0}\in P_{h}\), making successively the sequences \(x_{n}=T(x_{n-1},y_{n-1})\), \(y_{n}=T(y_{n-1},x_{n-1})\), \(n=0,1,2,\ldots\) , we obtain \(x_{n} \rightarrow u^{*}\) and \(y_{n}\rightarrow u^{*}\) as \(n\rightarrow \infty \).

References

  1. Kilbas, A.A., Srivastava, H.H., Trujillo, J.J.: Theory and Applications of Fractional Differential Equation. Elsevier, Amsterdam (2006)

    MATH  Google Scholar 

  2. Zhang, X., Liu, L., Wiwatanapataphee, B.: Nontrivial solutions for a fractional advection dispersion equation in anomalous diffusion. Appl. Math. Lett. 66, 1–8 (2017)

    Article  MathSciNet  Google Scholar 

  3. Agrawal, O.P.: Formulation of Euler–Lagrange equations for fractional variational problems. J. Math. Anal. Appl. 272, 368–379 (2002)

    Article  MathSciNet  Google Scholar 

  4. Machado, J.T., Kiryakova, V., Mainardi, F.: Recent history of fractional calculus. Commun. Nonlinear Sci. Numer. Simul. 16, 1140–1153 (2011)

    Article  MathSciNet  Google Scholar 

  5. He, J., Zhang, X., Liu, L., Wu, Y., Cui, Y.: A singular fractional Kelvin–Voigt model involving a nonlinear operator and their convergence properties. Bound. Value Probl. 2019, 112 (2019)

    Article  MathSciNet  Google Scholar 

  6. Zhang, X., Caccetta, L., Wu, Y.: Non local fractional order differential equations with changing-sign singular perturbation. Appl. Math. Model. 21(39), 6543–6552 (2015)

    Article  Google Scholar 

  7. Leibenson, L.S.: General problem of the movement of a compressible fluid in a porous medium. Izv. Akad. Nauk SSSR 9, 7–10 (1945)

    MathSciNet  Google Scholar 

  8. Guo, Y., Ji, Y., Liu, X.: Multiple positive solutions for some multi-point boundary value problems with p-Laplacian. J. Comput. Appl. Math. 216, 144–156 (2008)

    Article  MathSciNet  Google Scholar 

  9. Lu, H.L., Han, Z.L.: Existence on positive solutions for boundary value problems of nonlinear fractional differential equations with p-Laplacian. Adv. Differ. Equ. 2013, 30 (2013)

    Article  MathSciNet  Google Scholar 

  10. Ren, T., Li, S., Zhang, X., Liu, L.: Maximum and minimum solutions for a nonlocal p-Laplacian fractional differential system from eco-economical processes. Bound. Value Probl. 2017, 118 (2017)

    Article  MathSciNet  Google Scholar 

  11. Dong, X., Bai, Z., Zhang, S.: Positive solutions to boundary value problems of p-Laplacian with fractional derivative. Bound. Value Probl. 2017, 5 (2017)

    Article  MathSciNet  Google Scholar 

  12. Bai, C.: Existence and uniqueness of solutions for fractional boundary value problems with p-Laplacian operator. Adv. Differ. Equ. 2018, 4 (2018)

    Article  MathSciNet  Google Scholar 

  13. He, J., Zhang, X., Liu, L., Wu, Y., Cui, Y.: Existence and asymptotic analysis of positive solutions for a singular fractional differential equation with nonlocal boundary conditions. Bound. Value Probl. 2018, 189 (2018)

    Article  MathSciNet  Google Scholar 

  14. Zhang, X., Liu, L., Wu, Y.: Multiple positive solutions of a singular fractional differential equation with negatively perturbed term. Math. Comput. Model. 55, 1263–1274 (2012)

    Article  MathSciNet  Google Scholar 

  15. Zhou, B., Zhang, L.: Existence of positive solutions of boundary value problems for high-order nonlinear conformable differential equations with p-Laplacian operator. Adv. Differ. Equ. 2019, 351 (2019)

    Article  MathSciNet  Google Scholar 

  16. Zhou, B., Zhang, L., Addai, E., Zhang, N.: Multiple positive solutions for nonlinear high-order Riemann–Liouville fractional differential equations boundary value problems with p-Laplacian operator. Bound. Value Probl. 2020, 26 (2020)

    Article  MathSciNet  Google Scholar 

  17. Zhang, X., Liu, L.: The eigenvalue for a class of singular p-Laplacian fractional differential equations involving the Riemann–Stieltjes integral boundary condition. Appl. Math. Comput. 235, 412–422 (2014)

    MathSciNet  MATH  Google Scholar 

  18. Guo, D., Lakshmikantham, V.: Nolinear Problems in Abstact Cones. Academic Press, New York (1998)

    Google Scholar 

  19. Krasnoselskii, M.A.: Positive Solutions of Operator Equations. Noordhoff, Groningen (1964)

    Google Scholar 

  20. Zhai, C., Hao, M.: Fixed point theorems for mixed monotone operators with perturbation and applications to fractional differential equation boundary value problems. Nonlinear Anal. 75, 2542–2551 (2012)

    Article  MathSciNet  Google Scholar 

  21. Zhai, C., Yan, W., Yang, C.: A sum operator method for the existence and uniqueness of positive solutions to Riemann–Liouville fractional differential equation boundary value problems. Commun. Nonlinear Sci. Numer. Simul. 18, 858–866 (2013)

    Article  MathSciNet  Google Scholar 

Download references

Acknowledgements

We are thankful to the editor and the anonymous reviewers for many valuable suggestions to improve this paper.

Availability of data and materials

Data sharing is not applicable to this article as no datasets were generated or analyzed during the current study.

Funding

This paper is supported by the opening project of State Key Laboratory of Explosion Science and Technology (Beijing Institute of Technology). The opening project number is KFJJ19-06M. And it is also supported by the Key R&D Program of Shanxi Province (International Cooperation, 201903D421042).

Author information

Authors and Affiliations

Authors

Contributions

All authors contributed equally to the manuscript and typed, read, and approved the final manuscript.

Corresponding author

Correspondence to Lingling Zhang.

Ethics declarations

Competing interests

The authors declare that they have no competing interests.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhou, B., Zhang, L., Xing, G. et al. Existence–uniqueness and monotone iteration of positive solutions to nonlinear tempered fractional differential equation with p-Laplacian operator. Bound Value Probl 2020, 117 (2020). https://doi.org/10.1186/s13661-020-01414-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13661-020-01414-4

MSC

Keywords