Skip to main content

Discrete fractional order two-point boundary value problem with some relevant physical applications

Abstract

The results reported in this paper are concerned with the existence and uniqueness of solutions of discrete fractional order two-point boundary value problem. The results are developed by employing the properties of Caputo and Riemann–Liouville fractional difference operators, the contraction mapping principle and the Brouwer fixed point theorem. Furthermore, the conditions for Hyers–Ulam stability and Hyers–Ulam–Rassias stability of the proposed discrete fractional boundary value problem are established. The applicability of the theoretical findings has been demonstrated with relevant practical examples. The analysis of the considered mathematical models is illustrated by figures and presented in tabular forms. The results are compared and the occurrence of overlapping/non-overlapping has been discussed.

1 Introduction

Partial differential equations are invariably important in almost all fields of applied mathematics and science [13]. Particularly, one can observe that partial differential equations have been utilized in few places to help in the use of ordinary differential equations such as in the study of waves in liquids, propagation of sound, gravitational attraction and vibrations of strings [4]. On the other hand, partial fractional differential equations have presented adequate interpretations for many physical problems in areas such as fluid mechanics, biological populations, viscoelasticity, advection–diffusion, nuclear science and signals processing [5, 6]. In many cases, like the heat equation, wave equation, Poisson equation and Laplace equation, the problems remained unsolved due to the nonlinearity property of these equations. Owing to this limitation, various techniques like numerical methods for approximating solutions are used to problems modeled by nonlinear partial differential equations involving initial and boundary conditions [712].

Very recently, fractional differential equations (FDEs) have become intensively rich theory and found applications in various fields. These equations, which involve derivatives or integrals of fractional order, have resulted in a great interest for many researchers due to their effective applications in physics, chemistry, chaotic dynamical systems and random walks with memory in different fields of applied mathematics and engineering. Particularly emphasis has been put to the topics on existence, uniqueness and stability of solutions of differential equations of fractional order; see [1325] and the references cited therein. The corresponding discrete counter part, fractional order difference equations (FODEs), have appeared as a new research area for mathematicians and scientists. The study of discrete fractional calculus was initiated by Miller and Ross [26] and then developed by several other researchers [2741]. In the meantime, researchers have adopted the fact that dealing with FODEs provides a more accurate description than FDEs and the use of FODEs facilitates applications that require computational and simulation analysis.

The organization of the remaining part of the paper is outlined as follows: Fundamental definitions and concepts are introduced in Sect. 2. Section 3 is devoted to the discussion on existence and uniqueness results for a discrete FBVP (3.1). The main results of this section are obtained by using the contraction mapping principle and the Brouwer fixed point theorem. In Sect. 4, we develop conditions for Hyers–Ulam and Hyers–Ulam–Rassias stability of the discrete FBVP. The applications are discussed in Sect. 5 which is followed by our conclusion.

2 Auxiliary preliminaries

Now we present some fundamental definitions and essential lemmas of discrete fractional calculus that are to be used throughout this paper.

Definition 2.1

(see [31, 32])

Let \(\alpha > 0\). The αth fractional sum of a function Ψ is defined as

$$\begin{aligned} \Delta ^{-\alpha } \Psi (x)= \frac{1}{\Gamma (\alpha )} \sum _{\ell =a}^{x- \alpha } (x-\ell -1)^{(\alpha -1)} \Psi (\ell ), \end{aligned}$$

for all \(x\in \{a+\alpha , a+\alpha +1, \ldots \}:= {\mathbb{N}_{a+\alpha }}\) and \(x^{(\alpha )}:=\frac{\Gamma (x+1)}{\Gamma (x+1-\alpha )} \).

Definition 2.2

(see [3032])

Let \(\alpha >0\) and set \(\mu =n-\alpha \). The αth fractional Caputo difference operator is defined as

$$\begin{aligned} {}^{C}_{0}\Delta _{x}^{\alpha } \Psi (x)= \Delta ^{-\mu } \bigl(\Delta ^{n} \Psi (x) \bigr) = \frac{1}{\Gamma (\mu )} \sum_{\ell =a}^{x-\mu } (x- \ell -1)^{(\mu -1)}\Delta ^{n} \Psi (\ell ), \end{aligned}$$

for all \(x \in {\mathbb{N}_{a+\mu }}\) and \(n-1<\alpha \leq n\), where \(n= \lceil \alpha \rceil \) and \(\lceil . \rceil \) is the ceiling of a number.

Lemma 2.3

(see [28, 39])

Let x and α be any numbers for which \(x^{(\alpha )}\)and \(x^{(\alpha -1)}\)are defined. Then \(\Delta x^{(\alpha )}=\alpha x^{(\alpha -1)}\).

Lemma 2.4

(see [38, 39])

Let \(0\leq N-1<\alpha \leq N\). Then

$$\begin{aligned} \Delta ^{-\alpha } {}^{RL}_{0}\Delta _{x}^{\alpha }\Psi (x) =&\Psi (x)- \sum _{j=0}^{N-1}\frac{x^{(\alpha -N+j)}}{\Gamma (\alpha +j-1)} \Delta ^{j} \bigl[\Delta ^{\alpha -N}\Psi (a) \bigr], \\ =&\Psi (x)+B_{1} x^{(\alpha -1)}+B_{2} x^{(\alpha -2)} + \cdots +B_{N} x^{(\alpha -N)}, \end{aligned}$$

for \(B_{i} \in \mathbb{R}\), where \(i=1,2,\ldots ,N\).

Lemma 2.5

(see [27, 31, 32])

Suppose that \(\alpha >0\)and Ψ is defined on \({\mathbb{N}_{a}}\). Then

$$\begin{aligned} \Delta ^{-\alpha } {}^{C}_{0}\Delta _{x}^{\alpha }\Psi (x) =&\Psi (x)- \sum _{j=0}^{n-1}\frac{(x-a)^{(j)}}{j!}\Delta ^{j} \Psi (a) \\ =&\Psi (x)+C_{0}+C_{1} x+\cdots+C_{n-1} x^{(n-1)}, \end{aligned}$$

for \(C_{i} \in {\mathbb{R}}\), where \(i=0,1,2,\ldots ,n-1\).

Lemma 2.6

(see [30, 37])

If α and x are any numbers, then

  1. 1

    \(\sum_{\ell =0}^{x-\alpha } (x-\ell -1 )^{(\alpha -1)}= \frac{\Gamma (x+1)}{\alpha \Gamma (x-\alpha +1)} \).

  2. 2

    \(\sum_{\ell =0}^{L} (\alpha +L-\ell -1 )^{( \alpha -1)}=\frac{1}{\alpha }\frac{\Gamma (\alpha +L+1)}{\Gamma (L+1)}\).

Lemma 2.7

(see [29])

Let \(\mu \in \mathbb{R}\backslash \{\ldots ,-2,-1\}\). Then

$$ \Delta ^{-\alpha }x^{(\mu )}= \frac{\Gamma (\mu +1)}{\Gamma (\mu +\alpha +1)}x^{(\alpha +\mu )}. $$

3 Existence and uniqueness of solutions

In this section, we will discuss the existence and uniqueness of solutions to a discrete fractional boundary value problem (FBVP) of the form

$$ \textstyle\begin{cases} {}^{C}_{0}\Delta _{x}^{\alpha }w(x)=\Psi (x+\alpha -1, w(x+\alpha -1)), \quad 1< \alpha \leq 2, \\ \Delta w(\alpha -2)=A, \qquad w(\alpha +L)=B, \end{cases} $$
(3.1)

for \(x\in [0,L]_{\mathbb{N}_{0}}= \{ 0,1,\ldots ,L \} \), A, B are some real constants, \(\Psi : [\alpha -2, \alpha +L]_{\mathbb{N}_{\alpha -2}} \times \mathbb{R} \rightarrow \mathbb{R}\) is continuous function, \({}^{C}_{0}\Delta _{x}^{\alpha }\) is the Caputo fractional difference operator (CFDO) and \(L\in \mathbb{N}_{1}\). Now, we state and prove an important theorem which will be helpful to obtain a form of the solution of (3.1), provided that the solution exists.

Theorem 3.1

Let \(1<\alpha \leq 2\)and \(\Psi : [\alpha -2, \alpha +L]_{\mathbb{N}_{\alpha -2}} \rightarrow \mathbb{R}\)be given. Then a function w is a solution to the discrete FBVP

$$ \textstyle\begin{cases} {}^{C}_{0}\Delta _{x}^{\alpha }w(x)=\Psi (x+\alpha -1),\quad x\in [0,L]_{ \mathbb{N}_{0}}, \\ \Delta w(\alpha -2)=A,\qquad w(\alpha +L)=B, \end{cases} $$
(3.2)

if and only if \(w(x)\), for \(x \in [\alpha -2, \alpha +L]_{\mathbb{N}_{\alpha -2}}\)is a solution to the following fractional Taylor’s difference formula:

$$ \begin{aligned} w(x)&= p(x)+\frac{1}{\Gamma (\alpha )}\sum _{\ell =0}^{x- \alpha }(x-\ell -1)^{(\alpha -1)}\Psi (\ell +\alpha -1) \\ &\quad{}- \frac{1}{\Gamma (\alpha )}\sum_{\ell =0}^{L}( \alpha +L-\ell -1)^{( \alpha -1)}\Psi (\ell +\alpha -1), \end{aligned} $$
(3.3)

where p is the unique solution to the discrete FBVP

$$ \textstyle\begin{cases} {}^{C}_{0}\Delta _{x}^{\alpha }p(x)=0, \\ \Delta p(\alpha -2)=A, \qquad p(\alpha +L)=B. \end{cases} $$
(3.4)

Proof

Suppose that p is a solution to (3.4). Using Definition 2.1 together with Lemma 2.5 shows that

$$ p(x)=C_{0}+C_{1} x, \quad x \in [\alpha -2, \alpha +L]_{\mathbb{N}_{ \alpha -2}}, $$
(3.5)

where \(C_{0}, C_{1} \in \mathbb{R}\). Applying the operator Δ to both sides in (3.5), we get

$$ \Delta p(x)= C_{1}. $$
(3.6)

Using the boundary conditions \(\Delta p(\alpha -2)=A\) and \(p(\alpha +L)=B\) in (3.5) and (3.6), then it turns out that \(C_{0}=B-A(\alpha +L)\) and \(C_{1}=A\). Using \(C_{0}\) and \(C_{1}\) in \(p(x)\), we are left with

$$ p(x)= A\bigl[x-(\alpha +L)\bigr]+B. $$
(3.7)

Let \(w(x)\) be a solution to (3.2). In view of Lemma 2.5, we obtain a general solution to (3.2) in the form

$$ w(x)=\Delta ^{-\alpha } \Psi (x+\alpha -1)+C_{2}+C_{3} x, x \in [ \alpha -2, \alpha +L]_{\mathbb{N}_{\alpha -2}}, $$

where \(C_{2}, C_{3} \in \mathbb{R}\). Whence, by Definition 2.1, we have

$$ w(x)= \frac{1}{\Gamma (\alpha )}\sum_{\ell =0}^{x-\alpha }(x- \ell -1)^{( \alpha -1)}\Psi (\ell +\alpha -1)+C_{2}+C_{3} x. $$
(3.8)

Applying the operator Δ on both sides in (3.8), we get

$$ \Delta w(x)= \frac{1}{\Gamma (\alpha -1)}\sum_{\ell =0}^{x-\alpha +1}(x- \ell -1)^{(\alpha -2)}\Psi (\ell +\alpha -1)+ C_{3}. $$
(3.9)

By the boundary conditions \(\Delta w(\alpha -2)=A\) and \(w(\alpha +L)=B\), we obtain \(C_{2}=B-A(\alpha +L)-\frac{1}{\Gamma (\alpha )}\sum_{\ell =0}^{L}( \alpha +L-\ell -1)^{(\alpha -1)}\Psi (\ell +\alpha -1)\) and \(C_{3}=A\). Using the values of \(C_{2}\), \(C_{3}\) and \(p(x)\) in \(w(x)\), we deduce that

$$ \begin{aligned} w(x)&= p(x)+\frac{1}{\Gamma (\alpha )}\sum _{\ell =0}^{x- \alpha }(x-\ell -1)^{(\alpha -1)}\Psi (\ell +\alpha -1) \\ &\quad{}- \frac{1}{\Gamma (\alpha )}\sum_{\ell =0}^{L}( \alpha +L-\ell -1)^{( \alpha -1)}\Psi (\ell +\alpha -1). \end{aligned} $$
(3.10)

Conversely, it is easy to show that the solution (3.10) satisfies the discrete FBVP (3.2). The proof of the theorem is complete. □

For applications using the contraction mapping principle and the Brouwer fixed point theorems, the following operator is defined:

$$ \begin{aligned} (Tw) (x)&= p(x)+\frac{1}{\Gamma (\alpha )}\sum _{ \ell =0}^{x-\alpha }(x-\ell -1)^{(\alpha -1)}\Psi \bigl(\ell +\alpha -1, w( \ell +\alpha -1)\bigr) \\ &\quad{}- \frac{1}{\Gamma (\alpha )}\sum_{\ell =0}^{L}( \alpha +L-\ell -1)^{( \alpha -1)}\Psi \bigl(\ell +\alpha -1, w(\ell +\alpha -1) \bigr), \end{aligned} $$
(3.11)

for \(x\in [\alpha -2, \alpha +L]_{\mathbb{N}_{\alpha -2}}\). Obviously, \(w(x)\) is a solution to (3.1) if it is a fixed point of the operator T. For our convenience, we consider the Banach space \(\mathbb{E}\) with norm \(\Vert w \Vert =\max \vert w(x) \vert \) for \(x\in [\alpha -2, \alpha +L]_{\mathbb{N}_{\alpha -2}}\).

Theorem 3.2

Assume the following.

(\(H_{1}\)):

There exists a constant \(K>0\)such that \(\vert \Psi (x,w)-\Psi (x,w_{1}) \vert \leq K \vert w-w_{1} \vert \)for each \(x \in [\alpha -2, \alpha +L]_{\mathbb{N}_{\alpha -2}}\)and all \(w,w_{1}\in \mathbb{E}\).

Then the discrete FBVP (3.1) has a unique solution on \(\mathbb{E}\)provided that

$$ \begin{aligned} \frac{2\Gamma (\alpha +L+1)}{\Gamma (\alpha +1)\Gamma (L+1)}< \frac{1}{K}. \end{aligned} $$
(3.12)

Proof

Let \(w,w_{1}\in \mathbb{E}\). Then, for each \(x\in [\alpha -2, \alpha +L]_{\mathbb{N}_{\alpha -2}}\), we have

$$ \begin{aligned} & \bigl\vert (Tw) (x)-(Tw_{1}) (x) \bigr\vert \\ &\quad \leq \frac{1}{\Gamma (\alpha )}\sum_{\ell =0}^{x-\alpha }(x- \ell -1)^{( \alpha -1)} \\ & \qquad {}\times \bigl\vert \Psi \bigl(\ell +\alpha -1, w(\ell +\alpha -1) \bigr)-\Psi \bigl(\ell + \alpha -1, w_{1}(\ell +\alpha -1)\bigr) \bigr\vert \\ &\quad {}+\frac{1}{\Gamma (\alpha )}\sum_{\ell =0}^{L}( \alpha +L-\ell -1)^{( \alpha -1)} \\ & \qquad {}\times \bigl\vert \Psi \bigl(\ell +\alpha -1, w(\ell +\alpha -1) \bigr)-\Psi \bigl(\ell + \alpha -1, w_{1}(\ell +\alpha -1)\bigr) \bigr\vert \\ &\quad \leq \frac{K}{\Gamma (\alpha )}\sum_{\ell =0}^{x-\alpha }(x- \ell -1)^{( \alpha -1)} \bigl\vert w(\ell +\alpha -1)- w_{1}(\ell + \alpha -1) \bigr\vert \\ &\qquad {}+\frac{K}{\Gamma (\alpha )}\sum_{\ell =0}^{L}( \alpha +L-\ell -1)^{( \alpha -1)} \bigl\vert w(\ell +\alpha -1)-w_{1}( \ell +\alpha -1) \bigr\vert . \end{aligned} $$

It follows that

$$ \begin{aligned} & \Vert Tw-Tw_{1} \Vert \\ &\quad \leq \frac{K \Vert w-w_{1} \Vert }{\Gamma (\alpha )}\sum_{\ell =0}^{x- \alpha }(x- \ell -1)^{(\alpha -1)}+ \frac{K \Vert w-w_{1} \Vert }{\Gamma (\alpha )}\sum_{\ell =0}^{L}( \alpha +L-\ell -1)^{(\alpha -1)} \\ &\quad \leq \frac{K \Vert w-w_{1} \Vert }{\Gamma (\alpha )}\cdot \frac{\Gamma (x+1)}{\alpha \Gamma (x-\alpha +1)}+ \frac{K \Vert w-w_{1} \Vert }{\Gamma (\alpha )} \frac{\Gamma (\alpha +L+1)}{\alpha \Gamma (L+1)} \\ &\quad \leq \biggl[ \frac{2\Gamma (\alpha +L+1)}{\Gamma (\alpha +1)\Gamma (L+1)} \biggr] K \Vert w-w_{1} \Vert , \end{aligned} $$

which implies that T is a contraction. By the contraction mapping principle, T has a unique fixed point which is a unique solution to the discrete FBVP (3.1). The proof is complete. □

Theorem 3.3

Assume that \(\Psi :[\alpha -2, \alpha +L]_{\mathbb{N}_{\alpha -2}}\times \mathbb{R}\rightarrow \mathbb{R}\)is continuous and \(M \geq \max_{x\in [\alpha -2, \alpha +L]} \vert p(x) \vert \), where p is the unique solution of the discrete FBVP (3.4). Let \(Q= \max \{ \vert \Psi (x,w) \vert : x\in [\alpha -2, \alpha +L]_{\mathbb{N}_{\alpha -2}}, w\in \mathbb{E}, \vert w \vert \leq 2M \} \). Then the discrete FBVP (3.1) has a solution provided

$$ \begin{aligned} Q\leq \frac{M\Gamma (\alpha +1)\Gamma (L+1)}{2\Gamma (\alpha +L+1)}. \end{aligned} $$
(3.13)

Proof

Let \(M>0\) and we define the set \(S= \{ w(x)\in \mathbb{E} : \Vert w \Vert \leq 2M \} \). To prove this theorem, we only need to show that T maps S in S. For \(w(x)\in S\), we have

$$\begin{aligned}& \begin{aligned} \bigl\vert (Tw) (x) \bigr\vert &= \Biggl\vert p(x)+ \frac{1}{\Gamma (\alpha )}\sum_{\ell =0}^{x-\alpha }(x- \ell -1)^{( \alpha -1)}\Psi \bigl(\ell +\alpha -1, w(\ell +\alpha -1)\bigr) \\ & \quad {}-\frac{1}{\Gamma (\alpha )}\sum_{\ell =0}^{L}( \alpha +L- \ell -1)^{(\alpha -1)}\Psi \bigl(\ell +\alpha -1, w(\ell +\alpha -1) \bigr) \Biggr\vert \\ &\leq \bigl\vert p(x) \bigr\vert +\frac{1}{\Gamma (\alpha )}\sum _{\ell =0}^{x- \alpha }(x-\ell -1)^{(\alpha -1)} \bigl\vert \Psi \bigl(\ell +\alpha -1, w(\ell + \alpha -1)\bigr) \bigr\vert \\ &\quad {}+\frac{1}{\Gamma (\alpha )}\sum_{\ell =0}^{L}( \alpha +L-\ell -1)^{( \alpha -1)} \bigl\vert \Psi \bigl(\ell +\alpha -1, w(\ell +\alpha -1)\bigr) \bigr\vert \\ &\leq M + \frac{Q }{\Gamma (\alpha )} \Biggl[\sum_{\ell =0}^{x- \alpha }(x- \ell -1)^{(\alpha -1)}+\sum_{\ell =0}^{L}(\alpha +L-\ell -1)^{( \alpha -1)} \Biggr] \\ &\leq M + \frac{Q }{\Gamma (\alpha +1)} \biggl[ \frac{\Gamma (x+1)}{\Gamma (x+1-\alpha )}+ \frac{\Gamma (\alpha +L+1)}{\Gamma (L+1)} \biggr], \end{aligned} \\& \Vert Tw \Vert \leq M + 2Q \frac{\Gamma (\alpha +L+1)}{\Gamma (\alpha +1)\Gamma (L+1)}. \end{aligned}$$

From (3.13), we have \(\Vert Tw \Vert \leq 2M\), which implies that T maps S in S. Thus, T has at least one fixed point which is a solution to the BVP (3.1) according to the Brouwer fixed point theorem. □

4 Stability analysis

In this section, the stability analysis is presented for the following discrete FBVP:

$$ \textstyle\begin{cases} {}^{RL}_{0}\Delta _{x}^{\alpha }w(x)=\Psi (x+\alpha -1, w(x+\alpha -1)),\quad 1< \alpha \leq 2, \\ \Delta w(\alpha -2)=A,\qquad w(\alpha +L)=B, \end{cases} $$
(4.1)

for \(x\in [0,L]_{\mathbb{N}_{0}}\), where \({}^{RL}_{0}\Delta _{x}^{\alpha }\) is the Riemann–Liouville fractional difference operator (RLFDO). We now investigate the solution of (4.1), provided that the solution exists.

Theorem 4.1

Let \(1<\alpha \leq 2\)and \(\Psi : [\alpha -2, \alpha +L]_{\mathbb{N}_{\alpha -2}}\rightarrow \mathbb{R}\)be given. A solution to the discrete FBVP

$$ \textstyle\begin{cases} {}^{RL}_{0}\Delta _{x}^{\alpha }w(x)=\Psi (x+\alpha -1),\quad x\in [0, L]_{ \mathbb{N}_{0}}, \\ \Delta w(\alpha -2)=A, \qquad w(\alpha +L)=B , \end{cases} $$
(4.2)

has the form

$$ \begin{aligned} w(x)&= q(x)+ \frac{1}{\Gamma (\alpha )}\sum _{\ell =0}^{x- \alpha }(x-\ell -1)^{(\alpha -1)}\Psi (\ell +\alpha -1) \\ &\quad{}+\frac{u(x)}{\Gamma (\alpha )}\sum_{\ell =0}^{L}( \alpha +L- \ell -1)^{(\alpha -1)} \Psi (\ell +\alpha -1), \end{aligned} $$
(4.3)

where \(u(x)=\frac{\Gamma (L+2)}{\Gamma (\alpha +L+1)} [\beta \Gamma ( \alpha ) (x^{(\alpha -2)}-\frac{x^{(\alpha -1)}}{L+2} )-x^{( \alpha -1)} ]\)such that \(\beta = \frac{L+2}{(\alpha -2)(L+2)\Gamma (\alpha -1)-\Gamma (\alpha )}\)and q is the unique solution to the discrete FBVP

$$ \textstyle\begin{cases} {}^{RL}_{0}\Delta _{x}^{\alpha }q(x)=0, \\ \Delta q(\alpha -2)=A, \qquad q(\alpha +L)=B, \end{cases} $$
(4.4)

where \({}^{RL}_{0}\Delta _{x}^{\alpha }\)is the RLFDO.

Proof

Let q be a solution to (4.4) defined on \([\alpha -2, \alpha +L]_{\mathbb{N}_{\alpha -2}}\). Using Definition 2.1 and Lemma 2.4, we get

$$ \begin{aligned} q(x)= C_{4} x^{(\alpha -1)} +C_{5} x^{(\alpha -2)}, \end{aligned} $$
(4.5)

for some \(C_{4}\), \(C_{5} \in {\mathbb{R}}\). Applying the operator Δ on both sides in (4.5), we get

$$ \begin{aligned} \Delta u(x)= C_{4} (\alpha -1) x^{(\alpha -2)} +C_{5} ( \alpha -2) x^{(\alpha -3)}. \end{aligned} $$

Using the boundary conditions \(\Delta q(\alpha -2)=A\) and \(q(\alpha +L)=B\), we deduce that

$$ C_{4}=\frac{B\Gamma (L+2)}{\Gamma (\alpha +L+1)} \biggl[1+ \frac{\beta \Gamma (\alpha )}{L+2} \biggr]- \frac{A \beta }{L+2},\quad \text{and}\quad C_{5}=\beta A- \frac{\beta B \Gamma (\alpha )\Gamma (L+2)}{\Gamma (\alpha +L+1)}. $$

Substituting the values of \(C_{4}\) and \(C_{5}\) in q, we obtain

$$\begin{aligned} q(x)=\frac{B\Gamma (L+2)}{\Gamma (\alpha +L+1)}x^{(\alpha -1)}+ \beta \biggl(A- \frac{\Gamma (\alpha )\Gamma (L+2)}{\Gamma (\alpha +L+1)}B \biggr) \biggl[x^{(\alpha -2)}-\frac{x^{(\alpha -1)}}{L+2} \biggr]. \end{aligned}$$
(4.6)

Assume that w is a solution to (4.2). From Lemma 2.4, we obtain a general solution for (4.2) as

$$ \begin{aligned} w(x)=\frac{1}{\Gamma (\alpha )}\sum _{\ell =0}^{x- \alpha }(x-\ell -1)^{(\alpha -1)}\Psi (\ell +\alpha -1)+ C_{6} x^{( \alpha -1)} +C_{7} x^{(\alpha -2)}, \end{aligned} $$
(4.7)

for some \(C_{6}\), \(C_{7} \in {\mathbb{R}}\). Applying the operator Δ on both sides in (4.7), we get

$$ \begin{aligned} \Delta w(x)&= \frac{1}{\Gamma (\alpha -1)}\sum _{ \ell =0}^{x-\alpha +1}(x-\ell -1)^{(\alpha -2)}\Psi (\ell + \alpha -1) \\ &\quad {}+ C_{6} (\alpha -1) x^{(\alpha -2)} +C_{7} (\alpha -2) x^{(\alpha -3)}. \end{aligned} $$

In view of \(\Delta w(\alpha -2)=A\) and \(w(\alpha +L)=B\), we get the value of \(C_{6}\) and \(C_{7}\) as follows:

$$ \begin{aligned} C_{6}&=\frac{\Gamma (L+2)}{\Gamma (\alpha +L+1)} \biggl[1+ \frac{\beta \Gamma (\alpha )}{L+2} \biggr] \Biggl(B- \frac{1}{\Gamma (\alpha )}\sum _{\ell =0}^{L}(\alpha +L- \ell -1)^{(\alpha -1)}\Psi ( \ell +\alpha -1) \Biggr)\\ &\quad {}- \frac{A\beta }{L+2} \end{aligned} $$

and

$$ \begin{aligned} C_{7}=\beta A- \frac{\beta \Gamma (\alpha )\Gamma (L+2)}{\Gamma (\alpha +L+1)} \Biggl[B- \frac{1}{\Gamma (\alpha )}\sum_{\ell =0}^{L}(\alpha +L- \ell -1)^{(\alpha -1)}\Psi (\ell +\alpha -1) \Biggr]. \end{aligned} $$

Substituting the values of \(C_{6}\), \(C_{7}\) and q into (4.7), we obtain w in the form

$$ \begin{aligned} w(x)&= q(x)+ \frac{1}{\Gamma (\alpha )}\sum _{\ell =0}^{x- \alpha }(x-\ell -1)^{(\alpha -1)}\Psi (\ell +\alpha -1) \\ &\quad{}+\frac{u(x)}{\Gamma (\alpha )}\sum_{\ell =0}^{L}( \alpha +L- \ell -1)^{(\alpha -1)} \Psi (\ell +\alpha -1). \end{aligned} $$
(4.8)

 □

The definitions of Ulam stability for fractional difference equation are introduced in the sequel on the basis of [32, 37].

Definition 4.2

If, for every function \(v\in \mathbb{E}\) of

$$\begin{aligned} \bigl\vert {}^{RL}_{0}\Delta _{x}^{\alpha } v(x)-\Psi \bigl(x+\alpha -1, v(x+ \alpha -1)\bigr) \bigr\vert \leq \epsilon , \quad x\in [0, L]_{\mathbb{N}_{0}}, \end{aligned}$$
(4.9)

where \(\epsilon >0\), there exists a solution \(w\in \mathbb{E}\) of (4.1) and positive constant \(\delta _{1}>0\) such that

$$\begin{aligned} \bigl\vert v(x)-w(x) \bigr\vert \leq \delta _{1} \epsilon , \quad x\in [\alpha -2, \alpha +L]_{\mathbb{N}_{\alpha -2}}, \end{aligned}$$
(4.10)

then the discrete FBVP (4.1) is said to be Hyers–Ulam stable.

Definition 4.3

If, for every function \(v\in \mathbb{E}\) of

$$\begin{aligned} \bigl\vert {}^{RL}_{0}\Delta _{x}^{\alpha } v(x)-\Psi \bigl(x+\alpha -1, v(x+ \alpha -1)\bigr) \bigr\vert \leq \epsilon \Phi (x+\alpha -1),\quad x\in [0, L]_{ \mathbb{N}_{0}}, \end{aligned}$$
(4.11)

where \(\epsilon >0\), there is a solution \(w\in \mathbb{E}\) of (4.1) and positive constant \(\delta _{2}>0\) such that

$$\begin{aligned} \bigl\vert v(x)-w(x) \bigr\vert \leq \delta _{2} \epsilon \Phi (x+\alpha -1), \quad x\in [\alpha -2, \alpha +L]_{\mathbb{N}_{\alpha -2}}, \end{aligned}$$
(4.12)

then the discrete FBVP (4.1) is said to be Hyers–Ulam–Rassias stable.

Remark 4.4

A function \(v\in \mathbb{E}\) is a solution to (4.9) if and only if there exists \(g:[\alpha -2, \alpha +L]\rightarrow \mathbb{R}\) satisfying

\((H_{2})\):

\(\vert g(x+\alpha -1) \vert \leq \epsilon \), \(x \in [0, b]_{ \mathbb{N}_{0}}\),

\((H_{3})\):

\({}^{RL}_{0}\Delta _{x}^{\alpha } v(x)=\Psi (x+\alpha -1, v(x+\alpha -1))+g(x+ \alpha -1)\), \(x \in [0, b]_{\mathbb{N}_{0}}\).

A similar remark can be formulated for inequality (4.11).

Lemma 4.5

If v solves (4.9), then

$$\begin{aligned}& \Biggl\vert v(x) - q(x) - \frac{1}{\Gamma (\alpha )}\sum_{\ell =0}^{x- \alpha }(x- \ell -1)^{(\alpha -1)}\Psi \bigl(\ell +\alpha -1, v(\ell + \alpha -1)\bigr) \\& \qquad{} - \frac{u(x)}{\Gamma (\alpha )} \sum_{\ell =0}^{L}( \alpha +L-\ell -1)^{(\alpha -1)} \Psi \bigl(\ell +\alpha -1, v(\ell + \alpha -1) \bigr) \Biggr\vert \\& \quad \leq \frac{\epsilon \Gamma (\alpha +L+1)}{\Gamma (\alpha +1)\Gamma (L+1)}. \end{aligned}$$

Proof

If v solves the inequality (4.9), then from Remark 4.4 and Lemma 2.4, the solution to \((H_{3})\) satisfies

$$\begin{aligned} v(x) =& q(x) + \frac{1}{\Gamma (\alpha )}\sum_{\ell =0}^{x-\alpha }(x- \ell -1)^{(\alpha -1)}\Psi \bigl(\ell +\alpha -1, v(\ell +\alpha -1)\bigr) \\ &{} + \frac{u(x)}{\Gamma (\alpha )} \sum_{\ell =0}^{L}( \alpha +L- \ell -1)^{(\alpha -1)} \Psi \bigl(\ell +\alpha -1, v(\ell +\alpha -1) \bigr) \\ &{} +\frac{1}{\Gamma (\alpha )}\sum_{\ell =0}^{x-\alpha }(x- \ell -1)^{( \alpha -1)}g(\ell +\alpha -1). \end{aligned}$$

Hence,

$$\begin{aligned}& \Biggl\vert v(x) - q(x) - \frac{1}{\Gamma (\alpha )}\sum_{\ell =0}^{x- \alpha }(x- \ell -1)^{(\alpha -1)}\Psi \bigl(\ell +\alpha -1, v(\ell + \alpha -1)\bigr) \\& \qquad{} - \frac{u(x)}{\Gamma (\alpha )} \sum_{\ell =0}^{L}( \alpha +L-\ell -1)^{(\alpha -1)} \Psi \bigl(\ell +\alpha -1, v(\ell + \alpha -1) \bigr) \Biggr\vert \\& \quad = \Biggl\vert \frac{1}{\Gamma (\alpha )}\sum_{\ell =0}^{x-\alpha }(x- \ell -1)^{( \alpha -1)}g(\ell +\alpha -1) \Biggr\vert \\& \quad \leq \frac{1}{\Gamma (\alpha )}\sum_{\ell =0}^{x-\alpha }(x- \ell -1)^{( \alpha -1)} \bigl\vert g(\ell +\alpha -1) \bigr\vert \\& \quad \leq \frac{\epsilon }{\Gamma (\alpha )} \frac{\Gamma (x+1)}{\alpha \Gamma (x+1-\alpha )} \\& \quad \leq \frac{\epsilon \Gamma (\alpha +L+1)}{\Gamma (\alpha +1)\Gamma (L+1)}. \end{aligned}$$

 □

Theorem 4.6

Suppose that the hypothesis (\(H_{1}\)) together with the inequality (4.9) is satisfied. Then the discrete FBVP (4.1) is Hyers–Ulam stable provided that

$$ \begin{aligned} K< \frac{\Gamma (\alpha +1)\Gamma (L+1)}{2\Gamma (\alpha +L+1)}. \end{aligned} $$
(4.13)

Proof

With the help of solution (4.8) and Lemma 4.5, for \(x\in [\alpha -2, \alpha +L]_{\mathbb{N}_{\alpha -2}}\), we have

$$ \begin{aligned} & \bigl\vert v(x) - w(x) \bigr\vert \\ &\quad \leq \Biggl\vert v(x) - q(x) - \frac{1}{\Gamma (\alpha )}\sum _{ \ell =0}^{x-\alpha }(x-\ell -1)^{(\alpha -1)}\Psi \bigl(\ell +\alpha -1, w( \ell +\alpha -1)\bigr) \\ & \qquad {}- \frac{u(x)}{\Gamma (\alpha )} \sum_{\ell =0}^{L}( \alpha +L-\ell -1)^{(\alpha -1)} \Psi \bigl(\ell +\alpha -1, w(\ell + \alpha -1) \bigr) \Biggr\vert \\ &\quad \leq \Biggl\vert v(x) - q(x) - \frac{1}{\Gamma (\alpha )}\sum _{ \ell =0}^{x-\alpha }(x-\ell -1)^{(\alpha -1)}\Psi \bigl(\ell +\alpha -1, v( \ell +\alpha -1)\bigr) \\ & \qquad {}- \frac{u(x)}{\Gamma (\alpha )} \sum_{\ell =0}^{L}( \alpha +L-\ell -1)^{(\alpha -1)} \Psi \bigl(\ell +\alpha -1, v(\ell + \alpha -1) \bigr) \Biggr\vert \\ &\qquad {}+ \frac{1}{\Gamma (\alpha )}\sum_{\ell =0}^{x-\alpha }(x- \ell -1)^{(\alpha -1)} \\ &\qquad {}\times \bigl\vert \Psi \bigl(\ell +\alpha -1, v(\ell + \alpha -1)\bigr) - \Psi \bigl( \ell +\alpha -1, w(\ell +\alpha -1)\bigr) \bigr\vert \\ &\qquad {}+\frac{ \vert u(x) \vert }{\Gamma (\alpha )} \sum_{\ell =0}^{L}( \alpha +L-\ell -1)^{(\alpha -1)} \\ &\qquad{} \times \bigl\vert \Psi \bigl(\ell +\alpha -1, v(\ell + \alpha -1)\bigr) - \Psi \bigl( \ell +\alpha -1, w(\ell +\alpha -1)\bigr) \bigr\vert . \end{aligned} $$

It follows that

$$ \begin{aligned} \bigl\vert v(x) - w(x) \bigr\vert &\leq \frac{\epsilon \Gamma (\alpha +L+1)}{\Gamma (\alpha +1)\Gamma (L+1)} \\ &\quad{}+\frac{K}{\Gamma (\alpha )}\sum_{\ell =0}^{x-\alpha }(x- \ell -1)^{( \alpha -1)} \bigl\vert v(\ell +\alpha -1) - w(\ell +\alpha -1) \bigr\vert \\ &\quad{}+\frac{K \vert u(x) \vert }{\Gamma (\alpha )} \sum_{ \ell =0}^{L}( \alpha +L-\ell -1)^{(\alpha -1)} \bigl\vert v(\ell +\alpha -1) - w(\ell +\alpha -1) \bigr\vert \\ &\leq \frac{\epsilon \Gamma (\alpha +L+1)}{\Gamma (\alpha +1)\Gamma (L+1)}+ \frac{K \Vert v-w \Vert }{\Gamma (\alpha )} \biggl[ \frac{\Gamma (x+1)}{\alpha \Gamma (x+1-\alpha )} \biggr] \\ &\quad{}+\frac{K \Vert v-w \Vert \vert u(x) \vert }{\Gamma (\alpha )} \biggl[\frac{\Gamma (\alpha +L+1)}{\alpha \Gamma (L+1)} \biggr]. \end{aligned} $$

Therefore, we are left with

$$\begin{aligned} \Vert v - w \Vert \leq \frac{\epsilon \Gamma (\alpha +L+1)}{\Gamma (\alpha +1)\Gamma (L+1)}+ \frac{2K\Gamma (\alpha +L+1)}{\Gamma (\alpha +1)(L+1)} \Vert v-w \Vert . \end{aligned}$$

From the above inequalities, we have \(\Vert v- w \Vert \leq \delta _{1}\epsilon \), where \(\delta _{1}= \frac{\Gamma (\alpha +L+1)}{\Gamma (\alpha +1)\Gamma (L+1)-2K \Gamma (\alpha +L+1)}>0\). Thus, Eq. (4.1) is Hyers–Ulam stable. □

We assume the following.

(\(H_{4}\)):

Let \(\Phi \in [\alpha -2, \alpha +L]_{\mathbb{N}_{\alpha -2}}\rightarrow \mathbb{R}^{+}\) be an increasing function, and there exists a constant \(\lambda >0\) such that

$$\begin{aligned} \frac{\epsilon }{\Gamma (\alpha )}\sum_{\ell =0}^{x-\alpha }(x-\ell -1)^{( \alpha -1)}\Phi (\ell +\alpha -1)\leq \lambda \epsilon \Phi (x+ \alpha -1), \quad x\in [0, L]_{\mathbb{N}_{0}}. \end{aligned}$$

Lemma 4.7

If v solves (4.11), then

$$\begin{aligned}& \Biggl\vert v(x) - q(x)-\frac{1}{\Gamma (\alpha )}\sum_{\ell =0}^{x- \alpha }(x- \ell -1)^{(\alpha -1)}\Psi \bigl(\ell +\alpha -1, v(\ell + \alpha -1)\bigr) \\& \qquad{} - \frac{u(x)}{\Gamma (\alpha )} \sum_{\ell =0}^{L}( \alpha +L-\ell -1)^{(\alpha -1)} \Psi \bigl(\ell +\alpha -1, v(\ell + \alpha -1) \bigr) \Biggr\vert \\& \quad \leq \lambda \epsilon \Phi (x+\alpha -1). \end{aligned}$$

Proof

From inequality (4.11), for \(x\in [\alpha -2, \alpha +L]_{\mathbb{N}_{\alpha -2}}\), we can find a function \({}^{RL}_{0}\Delta _{x}^{\alpha } v(x)=\Psi (x+\alpha -1, v(x+\alpha -1))+g(x+ \alpha -1)\) and \(\vert g(x+\alpha -1) \vert \leq \epsilon \Phi (x+\alpha -1)\). It follows that

$$\begin{aligned}& \Biggl\vert v(x) - q(x) - \frac{1}{\Gamma (\alpha )}\sum_{\ell =0}^{x- \alpha }(x- \ell -1)^{(\alpha -1)}\Psi \bigl(\ell +\alpha -1, v(\ell + \alpha -1)\bigr) \\& \qquad{} - \frac{u(x)}{\Gamma (\alpha )} \sum_{\ell =0}^{L}( \alpha +L-\ell -1)^{(\alpha -1)} \Psi \bigl(\ell +\alpha -1, v(\ell + \alpha -1) \bigr) \Biggr\vert \\& \quad = \Biggl\vert \frac{1}{\Gamma (\alpha )}\sum_{\ell =0}^{x-\alpha }(x- \ell -1)^{( \alpha -1)}g(\ell +\alpha -1) \Biggr\vert \\& \quad \leq \frac{1}{\Gamma (\alpha )}\sum_{\ell =0}^{x-\alpha }(x- \ell -1)^{( \alpha -1)} \bigl\vert g(\ell +\alpha -1) \bigr\vert \\& \quad \leq \frac{\epsilon }{\Gamma (\alpha )}\sum_{\ell =0}^{x-\alpha }(x- \ell -1)^{(\alpha -1)}\Phi (\ell +\alpha -1) \\& \quad \leq \lambda \epsilon \Phi (x+\alpha -1). \end{aligned}$$

 □

Theorem 4.8

If the hypotheses (\(H_{1}\)), (\(H_{4}\)) and the inequality (4.13) are satisfied, then a discrete FBVP (4.1) is Hyers–Ulam–Rassias stable.

Proof

With the help of Lemma 2.6, Lemma 4.7 and solution (4.8), we obtain

$$ \begin{aligned} \Vert v- w \Vert \leq \delta _{2} \epsilon \Phi (x+ \alpha -1); \end{aligned} $$

where \(\delta _{2} = \frac{\lambda \Gamma (\alpha +1)\Gamma (L+1)}{\Gamma (\alpha +1)\Gamma (L+1)-2K\Gamma (\alpha +L+1)}>0\). Thus (4.1) is Hyers–Ulam–Rassias stable. □

5 Applications

Some suitable examples are presented to validate the theoretical results and numerical solutions to the discrete FBVP (3.1) and (4.1) with different applications by using Caputo and Riemann–Liouville fractional difference operators. Computational aspects regarding numerical values and diagrams are executed with MATLAB.

5.1 Steady-state heat equation

Consider the mathematical model of heat flow in a rod made out of a heat–conducting material, subject to an external heat source along its length and some boundary conditions at each end. Let \(w(x, t)\) denote the temperature distribution at the real position x varying only with any real time t, where \(a < x < b\) along some finite length of the rod. The general form of the one dimensional non-homogeneous heat equation with heat generation [4] is given as

$$ w_{t}(x,t)=\bigl(k(x)w_{x}(x,t) \bigr)_{x}+\psi (x,t), $$
(5.1)

where \(k(x)\) is the coefficient of heat conduction, which may vary with x, and \(\psi (x,t)\) is the heat source (or sink). Equation (5.1) is often called the diffusion equation. The heat equation (5.1) reduces to

$$ w_{t}(x,t)=k w_{xx}(x,t)+\psi (x,t), $$
(5.2)

subject to the initial condition

$$ w(x,0)=w(x), $$
(5.3)

and the Neumann–Dirichlet boundary conditions

$$ w_{x}(a,t)=A(t), \quad \text{and}\quad w(b,t)=B(t), $$
(5.4)

where the rod is assumed to be insulated at one end and the temperature is specified at the other end. In general, we expect the temperature distribution to change with time. However, if \(\psi (x,t)\), \(A(t)\) and \(B(t)\) are all time independent, Eqs. (5.2), (5.3) and (5.4) are called initial boundary value problems.

Now we are interested in computing the steady-state solution to the above problem and by setting \(w_{t} = 0\) in (5.2), we obtain the steady-state heat equation in x as the following ordinary differential equation:

$$ \textstyle\begin{cases} w^{\prime \prime }(x)=\Psi (x), \quad a< x< b, \\ w^{\prime }(a)=A , \qquad w(b)=B, \end{cases} $$
(5.5)

where \(\Psi (x)=-\frac{\psi (x)}{k}\), a, b, A and B are real valued constants. The above steady-state heat equation with Neumann–Dirichlet boundary conditions (5.5) is transformed to the discrete FBVPs (3.2) and (4.2).

Example 5.1

Suppose that \(\Psi (x)=x^{(8)}\), \(L=1\), \(A=0\) and \(B=0\) where we have different fractional orders α like \(\alpha =1.1\), \(\alpha =1.4\), \(\alpha =1.7\) and \(\alpha =2\). Then the discrete fractional order steady-state boundary value problems (3.2) and (4.2) become

$$ \textstyle\begin{cases} {}^{*}_{0}\Delta _{x}^{\alpha }w(x)= (x+\alpha -1 )^{(8)}, \quad x\in [0, 1]_{\mathbb{N}_{0}}, \\ \Delta w (\alpha -2 )=0, \qquad w (\alpha +1 )=0, \end{cases} $$
(5.6)

where \({}^{*}_{0}\Delta _{x}^{\alpha }={{}^{C}_{0}}\Delta _{x}^{\alpha }\) or \({{}^{*}_{0}}\Delta _{x}^{\alpha }={{}^{RL}_{0}}\Delta _{x}^{\alpha }\). The solutions of (5.6) can be formulated for different values of α by using the procedure discussed in the previous sections. Indeed, by using Definition 2.1 and Lemma 2.7 in (3.3) and (4.3) we obtain

$$\begin{aligned}& \begin{aligned} \frac{1}{\Gamma (\alpha )}\sum_{\ell =0}^{x-\alpha }(x- \ell -1)^{( \alpha -1)}(\ell +\alpha -1)^{(8)}&=\Delta ^{-\alpha }(x+ \alpha -1)^{(8)} \\ &=\frac{\Gamma (9)}{\Gamma (9+\alpha )}(x+\alpha -1)^{(\alpha +8)}, \end{aligned} \\& \frac{1}{\Gamma (\alpha )}\sum_{\ell =0}^{x-\alpha }(x-\ell -1)^{( \alpha -1)}(\ell +\alpha -1)^{(8)}=\frac{\Gamma (9)}{\Gamma (9+\alpha )} \cdot \frac{\Gamma (x+\alpha )}{\Gamma (x-8)}. \end{aligned}$$

Similarly, we find

$$\begin{aligned} \frac{1}{\Gamma (\alpha)}\sum_{\ell =0}^{L}(\alpha +L-\ell-1)^{( \alpha -1)}(\ell +\alpha -1)^{(8)}= \frac{\Gamma (9)}{\Gamma (9+\alpha )}\cdot \frac{\Gamma (2\alpha +1)}{\Gamma (\alpha -7)}. \end{aligned}$$

If \({{}^{*}_{0}}\Delta _{x}^{\alpha }= {{}^{C}_{0}}\Delta _{x}^{\alpha }\) then the analytical solution to (5.6) has the form

$$\begin{aligned} w(x)=p(x)+\frac{\Gamma (9)}{\Gamma (9+\alpha )} \biggl[ \frac{\Gamma (x+\alpha )}{\Gamma (x-8)} \biggr]- \frac{\Gamma (9)}{\Gamma (9+\alpha )} \biggl[ \frac{\Gamma (2\alpha +1)}{\Gamma (\alpha -7)} \biggr], \end{aligned}$$

where \(p(x)\) is defined in Theorem 3.1. If \({{}^{*}_{0}}\Delta _{x}^{\alpha }= {{}^{RL}_{0}}\Delta _{x}^{\alpha }\) then the analytical solution to (5.6) is

$$\begin{aligned} w(x)=q(x)+\frac{\Gamma (9)}{\Gamma (9+\alpha )} \biggl[ \frac{\Gamma (x+\alpha )}{\Gamma (x-8)} \biggr]+u(x) \frac{\Gamma (9)}{\Gamma (9+\alpha )} \biggl[ \frac{\Gamma (2\alpha +1)}{\Gamma (\alpha -7)} \biggr], \end{aligned}$$

where q and u are defined as in Theorem 4.1. For the values of the fractional orders \(\alpha \in (1,2]\) and \(x\in [0,1]\), the solutions with respect to the Riemann–Liouville operator exhibit better results than the Caputo difference operator, as seen in Fig. 1(a) and Fig. 1(b) and as shown in Table 1 and Table 2. When \(\alpha =2\), the graphs of both solutions are provided in Fig. 1(c). Note that both curves are overlapping. In three dimensions the solution surfacing over different values of α and x are shown in Fig. 2.

Figure 1
figure 1

Steady-state solutions curves for x using (a) CFDO, (b) RLFDO and (c) comparison of steady-state solutions at \(\alpha =2\) of Example 5.1

Figure 2
figure 2

Surface plots of corresponding steady-state solutions to different values of α and x using (a) CFDO and (b) RLFDO of Example 5.1

Table 1 Steady-state solutions to fractional order \(\alpha \in (1,2]\) varying \(x\in [0,1]\) using CFDO
Table 2 Steady-state solutions to fractional order \(\alpha \in (1,2]\) varying \(x\in [0,1]\) using RLFDO

Example 5.2

Suppose that \(\Psi (x)=x^{(6)}\), \(L=4\), \(A=0\) and \(B=1\) where we have different fractional orders α like \(\alpha =1.1\), \(\alpha =1.4\), \(\alpha =1.7\) and \(\alpha =2\). Then the discrete fractional order steady-state boundary value problems (3.2) and (4.2) become

$$ \textstyle\begin{cases} {}^{*}_{0}\Delta _{x}^{\alpha }w(x)= (x+\alpha -1 )^{(6)}, \quad x\in [0, 4]_{\mathbb{N}_{0}}, \\ \Delta x (\alpha -2 )=0, \qquad x (\alpha +4 )=1, \end{cases} $$
(5.7)

where \({}^{*}_{0}\Delta _{x}^{\alpha }\) is either \({{}^{C}_{0}}\Delta _{x}^{\alpha }\) or \({{}^{RL}_{0}}\Delta _{x}^{\alpha }\). By using a similar method to Example 5.1, we obtain the steady-state solutions to this problem. If \({{}^{*}_{0}}\Delta _{x}^{\alpha }= {{}^{C}_{0}}\Delta _{x}^{\alpha }\) then the analytical solution to (5.7) is

$$\begin{aligned} w(x)=p(x)+\frac{\Gamma (7)}{\Gamma (7+\alpha )} \biggl[ \frac{\Gamma (x+\alpha )}{\Gamma (x-6)} \biggr]- \frac{\Gamma (7)}{\Gamma (7+\alpha )} \biggl[ \frac{\Gamma (2\alpha +4)}{\Gamma (\alpha -2)} \biggr], \end{aligned}$$

where p is defined in Theorem 3.1. If \({{}^{*}_{0}}\Delta _{x}^{\alpha }= {{}^{RL}_{0}}\Delta _{x}^{\alpha }\) then the analytical solution to (5.7) is

$$\begin{aligned} w(x)=q(x)+\frac{\Gamma (7)}{\Gamma (7+\alpha )} \biggl[ \frac{\Gamma (x+\alpha )}{\Gamma (x-6)} \biggr]+u(x) \frac{\Gamma (7)}{\Gamma (7+\alpha )} \biggl[ \frac{\Gamma (2\alpha +4)}{\Gamma (\alpha -2)} \biggr], \end{aligned}$$

where q and u are defined as in Theorem 4.1. For the values of the fractional orders \(\alpha \in (1,2]\) and \(x\in [0,4]\), it is realized that the solutions with respect to the Riemann–Liouville operator exhibit better results than the Caputo difference operator, as seen in Fig. 3(a) and Fig. 3(b) and as shown in Table 3 and Table 4. As we see in Fig. 3, the trajectories of both solutions do not overlap for fractional orders \(\alpha \in (1,2]\) with x increasing. In three dimensions, the solution surfaces over different values of α and x are shown in Fig. 4.

Figure 3
figure 3

Steady-state solutions curve varies x using (a) CFDO and (b) RLFDO of Example 5.2

Figure 4
figure 4

Surface plots of corresponding steady-state solutions to different values of α and x using (a) CFDO and (b) RLFDO of Example 5.2

Table 3 Steady-state solutions to fractional order \(\alpha \in (1,2]\) varying \(x\in [0,4]\) using CFDO
Table 4 Steady-state solutions to fractional order \(\alpha \in (1,2]\) varying \(x\in [0,4]\) using RLFDO

5.2 Simple gravity pendulum

In [4], the following differential equation, which represents the motion of a simple pendulum, is considered:

$$ \frac{d^{2}w}{dx^{2}}+\frac{g}{\gamma }\sin \bigl(w(x)\bigr)=0, $$
(5.8)

where g is the acceleration due to the gravitational constant, γ is the length of the pendulum and w is the angular displacement.

Example 5.3

Let us consider the parameters \(A=0\), \(B=1\), \(L=1\) with fractional order \(\alpha =1.3\). Then the discrete fractional order boundary value problems (3.1) and (4.1) for the simple pendulum become

$$ \textstyle\begin{cases} {}^{*}_{0}\Delta _{x}^{1.3}w(x)= \Psi (x+0.3, w(x+0.3)), \quad x\in [0,1]_{ \mathbb{N}_{0}}, \\ \Delta w (-0.7 )=0, \qquad w (2.3 )=1, \end{cases} $$
(5.9)

where \({}^{*}_{0}\Delta _{x}^{1.3}={}^{C}_{0}\Delta _{x}^{1.3}\) or \({}^{*}_{0}\Delta _{x}^{1.3}= {}^{RL}_{0}\Delta _{x}^{1.3}\) and \(\Psi (x+0.3, w(x+0.3))=-\frac{g}{\gamma }\sin (w(x+0.3))\). For the values of \(g=9.8\mbox{ m/s}^{2}\), \(\gamma =60\mbox{ m}\), we choose \(K=\frac{g}{\gamma }=0.1633\). If \({}^{*}_{0}\Delta _{x}^{1.3}={}^{C}_{0}\Delta _{x}^{1.3}\), in this case, inequality (3.12) takes the form

$$ K \biggl[\frac{2\Gamma (\alpha +L+1)}{\Gamma (\alpha +1)\Gamma (L+1)} \biggr]\approx 0.7513< 1. $$

Therefore, from Theorem 3.2, we conclude that the boundary value problem (5.9) has a unique solution. Furthermore, if \({}^{*}_{0}\Delta _{x}^{1.3}= {}^{RL}_{0}\Delta _{x}^{1.3}\), we obtain

$$ \frac{\Gamma (\alpha +1)\Gamma (L+1)}{2\Gamma (\alpha +L+1)}\approx 0.2174. $$

If \(K=0.1633<0.2174\) and the inequality

$$\begin{aligned} \bigl\vert {}^{RL}_{0}\Delta _{x}^{1.3} v(x)-\Psi \bigl(x+0.3, v(x+0.3)\bigr) \bigr\vert \leq \epsilon , \quad x\in [0, 1]_{\mathbb{N}_{0}}, \end{aligned}$$

holds, then the boundary value problem (5.9) is Hyers–Ulam stable by Theorem 4.6.

5.3 Temperature distribution equation

In [8, 10, 16], the authors considered the following mathematical model, which describes the temperature distribution in lumped system of combined convection–radiation in a slab made of materials with variable thermal conductivity:

$$ \textstyle\begin{cases} D^{\alpha }w(x)-\eta w^{4}(x)=0,\quad 1< \alpha \leq 2, x\in [0,1], \\ w^{\prime }(0)=0 ,\qquad w(1)=1, \end{cases} $$
(5.10)

where \(w=\frac{\mathbb{T}-\mathbb{T}_{a}}{\mathbb{T}_{i} - \mathbb{T}_{a}}\) and \(x= \frac{t }{V\rho c_{a} /\mathbb{S}h}\) are dimensionless temperature and time, respectively. The parameter \(\eta = (\mathbb{T}-\mathbb{T}_{a} )\xi \), where V, \(\mathbb{S}\), ρ, c, \(\mathbb{T}_{i}\), \(\mathbb{T}_{a}\), \(c_{a}\) and h are the volume, surface area, density, specific heat, the initial temperature, temperature of the convection environment, specific heat at temperature \(\mathbb{T}_{a}\) and heat transfer coefficient of the lumped system, respectively.

Example 5.4

Suppose that \(\alpha =1.4\), \(L=1\), \(\eta = 4\times 10^{-10}\) and \(M=250\) with \(\Psi (x,w)=\eta w^{4}(x)\). Then we obtain the discrete fractional order heat transfer boundary value problem (5.10) and it takes the form

$$ \textstyle\begin{cases} {}^{C}_{0}\Delta _{x}^{1.4}w(x)=\eta w^{4} (x+0.4 ), \quad x\in [0,1], \\ \Delta w (-0.6 )=0, \qquad w (2.4 )=1. \end{cases} $$
(5.11)

The main result of (5.11) is discussed in Sect. 2. The Banach space is

$$ \mathbb{E}:= \bigl\{ w(x) | [-0.6, 2.4 ]_{\mathbb{N}_{-0.6}} \rightarrow \mathbb{R}, \Vert w \Vert \leq 500 \bigr\} . $$

We note that

$$ \frac{M\Gamma (\alpha +1)\Gamma (L+1)}{2\Gamma (\alpha +L+1)}=250(0.2083) \approx 52.0833. $$

It is clear that \(\vert \Psi (x,w) \vert \leq 25<52.0833\), whenever \(w\in [-500, 500]\). Therefore by Theorem 3.3, we conclude that the boundary value problem (5.11) has at least one solution.

6 Conclusion

This work made a study of a discrete fractional boundary value problem with the Caputo and Riemann–Liouville difference operators. The existence and uniqueness of solutions and various types of Hyers–Ulam stability are discussed for the addressed problem based on the properties of fractional operators, the contraction mapping principle and the Brouwer fixed point theorem. Theoretical results are complemented with suitable examples accompanied by numerical solutions for different values of α and x. The discrete FBVP appearing in mathematical models of engineering applications is solved via the Riemann–Liouville and Caputo difference operators. Subsequently, the dynamics exhibited by RLFDO and CFDO for the proposed models is illustrated through diagrams. Furthermore, the occurrence of overlapping or non-overlapping cases is presented in the graphical representations.

Results obtained in the present paper can be considered as a contribution to the developing field of discrete fractional boundary value problems describing mathematical and physical applications.

References

  1. Wu, J.H.: Theory and Applications of Partial Functional Differential Equations. Springer, New York (1996)

    MATH  Google Scholar 

  2. Kanwal, R.P.: Applications to partial differential equations. In: Linear Integral Equations. Birkhäuser, Boston (1997)

    MATH  Google Scholar 

  3. Evans, G.A., Blackledge, J.M., Yardley, P.D.: Analytic Methods for Partial Differential Equations. Springer, London (1999)

    MATH  Google Scholar 

  4. LeVeque, R.J.: Finite Difference Methods for Ordinary and Partial Differential Equations. SIAM, Philadelphia (2007)

    MATH  Google Scholar 

  5. Podlubny, I.: Fractional Differential Equations. Academic Press, New York (1999)

    MATH  Google Scholar 

  6. Atangana, A., Jafari, H., Belhaouari, S.B., Bayram, M.: Partial fractional equations and their applications. Math. Probl. Eng. 2015, 387205 (2015)

    Google Scholar 

  7. Mebrate, B.: Numerical solution of a one dimensional heat equation with Dirichlet boundary conditions. Am. J. Appl. Math. 3, 305–311 (2015)

    Google Scholar 

  8. Danish, M., Kumar, S., Kumar, S.: Exact solutions of three nonlinear heat transfer problems. Eng. Lett. 19, 1–6 (2011)

    Google Scholar 

  9. Khan, R.A.: Generalized approximation method for heat radiation equations. Appl. Math. Comput. 212, 287–295 (2009)

    MathSciNet  MATH  Google Scholar 

  10. Saeed, U., Rehman, M.: Assessment of Haar wavelet-quasilinearization technique in heat convection-radiation equations. Appl. Comput. Intell. Soft Comput. 2014, Article ID 454231 (2014)

    Google Scholar 

  11. Biazar, J., Ghazvini, H.: An analytical approximation to the solution of a wave equation by a variational iteration method. Appl. Math. Lett. 21, 780–785 (2008)

    MathSciNet  MATH  Google Scholar 

  12. Pardoux, E., Veretennikov, A.Y.: On the Poisson equation and diffusion approximation. Ann. Probab. 29, 1061–1085 (2001)

    MathSciNet  MATH  Google Scholar 

  13. Ahmad, M., Zada, A., Wang, X.: Existence, uniqueness and stability of implicit switched coupled fractional differential equations of ψ-Hilfer type. Int. J. Nonlinear Sci. Numer. Simul. 21, 1–11 (2020)

    MathSciNet  Google Scholar 

  14. Ahmad, M., Zada, A., Alzabut, J.: Stability analysis of a nonlinear coupled implicit switched singular fractional differential system with p-Laplacian. Adv. Differ. Equ. 2019, 436 (2019)

    MathSciNet  Google Scholar 

  15. Ahmad, M., Jiang, J., Zada, A., Shah, S.O., Xu, J.: Analysis of implicit coupled system of fractional differential equations involving Katugampola–Caputo fractional derivative. Complexity 2020, Article ID 9285686 (2020)

    MATH  Google Scholar 

  16. Ismail, M., Saeed, U., Alzabut, J., Rehman, M.: Approximate solutions for fractional boundary value problems via Green-CAS wavelet method. Mathematics 7, 1164 (2019)

    Google Scholar 

  17. Kilbas, A.A., Srivastava, H.M., Trujillo, J.J.: Theory and Applications of Fractional Differential Equations. Elsevier, New York (2006)

    MATH  Google Scholar 

  18. Zada, A., Alzabut, J., Waheed, H., Popa, I.-L.: Ulam–Hyers stability of impulsive integrodifferential equations with Riemann–Liouville boundary conditions. Adv. Differ. Equ. 2020, 64 (2020)

    MathSciNet  Google Scholar 

  19. Ray, S.S., Sahoo, S.: Generalized Fractional Order Differential Equations Arising in Physical Models. Chapman & Hall, London (2018)

    MATH  Google Scholar 

  20. Iswarya, M., Raja, R., Rajchakit, G., Cao, J., Alzabut, J., Huang, C.: Existence, uniqueness and exponential stability of periodic solution for discrete-time delayed bam neural networks based on coincidence degree theory and graph theoretic method. Mathematics 2019, 1055 (2019)

    Google Scholar 

  21. Yan, R.A., Sun, S.R., Han, Z.L.: Existence of solutions of boundary value problems for Caputo fractional differential equations on time scales. Bull. Iran. Math. Soc. 42, 247–262 (2016)

    MathSciNet  MATH  Google Scholar 

  22. Zhao, Y., Sun, S., Han, Z., Li, Q.: The existence of multiple positive solutions for boundary value problems of nonlinear fractional differential equations. Commun. Nonlinear Sci. Numer. Simul. 16, 2086–2097 (2011)

    MathSciNet  MATH  Google Scholar 

  23. Guo, Y., Shu, X., Li, Y., Xu, F.: The existence and Hyers–Ulam stability of solution for an impulsive Riemann–Liouville fractional neutral functional stochastic differential equation with infinite delay of order \(1 < \beta < 2\). Bound. Value Probl. 2019, 59 (2019)

    MathSciNet  Google Scholar 

  24. Ahmad, M., Zada, A., Alzabut, J.: Hyes–Ulam stability of coupled system of fractional differential equations of Hilfer–Hadamard type. Demonstr. Math. 52, 283–295 (2019)

    MathSciNet  MATH  Google Scholar 

  25. Salem, A., Alzahrani, F., Almaghamsi, L.: Fractional Langevin equations with nonlocal integral boundary conditions. Mathematics 402, 402 (2019)

    Google Scholar 

  26. Miller, K.S., Ross, B.: An Introduction to the Fractional Calculus and Fractional Differential Equations. Wiley, New York (1993)

    MATH  Google Scholar 

  27. Abdeljawad, T.: On Riemann and Caputo fractional differences. Comput. Math. Appl. 62, 1602–1611 (2011)

    MathSciNet  MATH  Google Scholar 

  28. Atici, F.M., Eloe, P.M.: Two-point boundary value problems for finite fractional difference equations. J. Differ. Equ. Appl. 17, 445–456 (2011)

    MathSciNet  MATH  Google Scholar 

  29. Atici, F.M., Eloe, P.M.: A transform method in discrete fractional calculus. Int. J. Difference Equ. 2, 165–176 (2007)

    MathSciNet  Google Scholar 

  30. Anastassiou, G.A.: Nabla discrete fractional calculus and nabla inequalities. Math. Comput. Model. 51, 562–571 (2010)

    MathSciNet  MATH  Google Scholar 

  31. Chen, H., Jin, Z., Kang, S.: Existence of positive solution for Caputo fractional difference equation. Adv. Differ. Equ. 2015, 44 (2015)

    MathSciNet  MATH  Google Scholar 

  32. Chen, F., Zhou, Y.: Existence and Ulam stability of solutions for discrete fractional boundary value problem. Discrete Dyn. Nat. Soc. 2013, Article ID 459161 (2013)

    MathSciNet  MATH  Google Scholar 

  33. Chen, C., Bohner, M., Jia, B.: Ulam–Hyers stability of Caputo fractional difference equations. Math. Methods Appl. Sci. 42, 7461–7470 (2019)

    MathSciNet  MATH  Google Scholar 

  34. Chen, C., Bohner, M., Jia, B.: Method of upper and lower solutions for nonlinear Caputo fractional difference equations and its applications. Fract. Calc. Appl. Anal. 22, 1307–1320 (2019)

    MathSciNet  MATH  Google Scholar 

  35. Chen, C., Bohner, M., Jia, B.: Existence and uniqueness of solutions for nonlinear Caputo fractional difference equations. Turk. J. Math. 44, 857–869 (2020)

    MathSciNet  Google Scholar 

  36. Selvam, A.G.M., Dhineshbabu, R.: Existence and uniqueness of solutions for a discrete fractional boundary value problem. Int. J. Appl. Math. 2, 283–295 (2020)

    Google Scholar 

  37. Selvam, A.G.M., Dhineshbabu, R.: Hyers–Ulam stability results for discrete antiperiodic boundary value problem with fractional order \(2<\delta\leq3\). Int. J. Eng. Adv. Technol. 9, 4997–5003 (2019)

    Google Scholar 

  38. Selvam, A.G.M., Dhineshbabu, R.: Ulam stability results for boundary value problem of fractional difference equations. Adv. Math. 9, 219–230 (2020)

    Google Scholar 

  39. Goodrich, C.S.: Existence and uniqueness of solutions to a fractional difference equation with nonlocal conditions. Comput. Math. Appl. 61, 191–202 (2011)

    MathSciNet  MATH  Google Scholar 

  40. Pan, Y., Han, Z., Sun, S., Hou, C.: The existence of solutions to a class of boundary value problems with fractional difference equations. Adv. Differ. Equ. 2013, 275 (2013)

    MathSciNet  MATH  Google Scholar 

  41. Alzabut, J., Abdeljawad, T., Baleanu, D.: Nonlinear delay fractional difference equations with applications on discrete fractional Lotka–Volterra competition model. J. Comput. Anal. Appl. 25, 889–898 (2018)

    MathSciNet  Google Scholar 

Download references

Acknowledgements

J. Alzabut would like to thank Prince Sultan University for funding this work through research group Nonlinear Analysis Methods in Applied Mathematics (NAMAM) group number RG-DES-2017-01-17.

Availability of data and materials

Not applicable.

Funding

Not applicable.

Author information

Authors and Affiliations

Authors

Contributions

All authors contributed equally and significantly to this paper. All authors have read and approved the final version of the manuscript.

Corresponding author

Correspondence to Jehad Alzabut.

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

Selvam, A.G.M., Alzabut, J., Dhineshbabu, R. et al. Discrete fractional order two-point boundary value problem with some relevant physical applications. J Inequal Appl 2020, 221 (2020). https://doi.org/10.1186/s13660-020-02485-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13660-020-02485-8

MSC

Keywords