Solving a Higher-Dimensional Time-Fractional Diffusion Equation via the Fractional Reduced Differential Transform Method

Abstract

In this study, exact and approximate solutions of higher-dimensional time-fractional diffusion equations were obtained using a relatively new method, the fractional reduced differential transform method (FRDTM). The exact solutions can be found with the benefit of a special function, and we applied Caputo fractional derivatives in this method. The numerical results and graphical representations specified that the proposed method is very effective for solving fractional diffusion equations in higher dimensions.

1. Introduction

Fractional calculus is a generalization of integration and differentiation to noninteger-order fundamental operator ${}_a{\mathbf{D}}_t^{\alpha }$ where a and t are the bounds of the operation and $\alpha \in \mathbf{R}$; this notation was designed by Harold T. Davis. Diverse definitions for fractional derivatives have been proposed such as Riemann–Liouville, Caputo, Hadamard, Erdélyi–Kober, Grünwald–Letnikov, Marchaud, and Riesz, to name a few. The three greatest regular definitions for the universal fractional differintegral are the Caputo, the Riemann–Liouville, and the Grünwald–Letnikov definition [1,2,3].

In this study, we used the Caputo fractional derivative; the binary significant explanations for that are the initial conditions for fractional-order differential equations in a form connecting only the limit values of integer-order derivatives at the lower terminal initial time [3]. Similarly, the fractional derivative of a constant function is zero.

Up to now, there have been diverse methods for solving fractional differential equations using different definitions of fractional derivatives. Let us indicate some of these applications: Almeida et al. [4] applied fractional differential equations for modeling particular real phenomena, though Bulut et al. [5] considered the nonlinear time-fractional Burgers equation via the improved Bernoulli subequation function method. Atangana et al. [6] considered an advection–dispersion model with a fractional order and fractal dimension. Alshammari et al. [7] proposed residual power series (RPS) to find the numerical solution of a class of fractional Bagley–Torvik problems (FBTP) arising in a Newtonian fluid. Similarly, Yépez-Martínez et al. [8] solved the nonlinear coupled spacetime-fractional mKdV partial differential equation using Feng’s first integral method. The definition of the beta fractional derivative to find exact and approximate solutions of time-fractional diffusion equations in different dimensions was modified in [9,10]. Youssri [11] adopted the spectral Tau method for solving the nonlinear Riccati initial-value problem with a new generalized Caputo FF derivative. Youssri et al. [12] presented the numerical solutions of the fractional pantograph differential equations (FPDEs) using generalized Lucas polynomials (GLPs). Abd-Elhameed et al. [13] presented an explicit formula that approximates the fractional derivatives of Chebyshev polynomials of the first kind in the Caputo sense. Abd-Elhameed and Youssri [14] derived novel formulae for the high-order derivatives of Chebyshev polynomials of the fifth kind.

Keskin and Oturanc [15] suggested the fractional reduced differential transform method (FRDTM). The FRDTM is one of the best common methods for solving fractional partial differential equations as the FRDTM is a generalization of the reduced differential transform method (RDTM), which in turn is a generalization of the differential transform method (DTM) for solving different types of differential equations. Scholars frequently try to find different methods to simplify the resulting solutions and decrease the solution steps, making the progress of mathematical techniques required to complete the greatest consequences. The applicability of the FRDTM to some diverse categories of fractional differential equations has been obtainable as follows: Mukhtar et al. [16] applied the FRDTM to solve nonlinear fractional Burgers equations in different dimensions. Gupta [17] offered the approximate analytical solutions of the Benney–Lin equation with a fractional time derivative. Srivastava et al. [18] used the FRDTM to obtain the exact solution of a mathematical model for the generalized time-fractional-order biological population model, for multiterm time-fractional diffusion equations. Abuasad et al. [19] suggested a modified method of the FRDTM, and the FRDTM approximate solution of the time-fractional Korteweg–de Vries equation was offered by Ebenezer et al. [20].

An application of the FRDTM to a system of linear and nonlinear fractional partial differential equations was organized by Singh [21]. Rawashdeh [22] applied the FRDTM to solve nonlinear fractional partial differential equations such as the spacetime-fractional Burgers’ equations and the time-fractional Cahn–Allen equations. Singh and Kumar [23] used the FRDTM to find approximate solutions of time-fractional-order multidimensional Navier–Stokes equations. The fractional Helmholtz equations were considered to find exact and approximate solutions via the FRDTM by Abuasad et al. [24]. Furthermore, Abuasad et al. [25] found approximate solutions of the fractional SSIS epidemic model using the fractional multistep differential transformed method.

The main benefits of the FRDTM are that it can be applied to diverse types of linear and nonlinear PFDEs in different dimensions. The multistep differential transform method (MsDTM) can overcome the central complications of the differential transform method (DTM) and reduced differential transform method (RDTM), which are that the achieved series solution often converges in the real irrelevant space and the range of convergence is a precise slow procedure or completely divergent given a wider space [25]. Abdou [26] used the FRDTM to develop a scheme to study the numerical solution of time-fractional nonlinear evolution equations under initial conditions. Comparing the DTM, RDTM, FRDTM, and MsDTM, we found that the DTM is an upgraded method of the Taylor series method, which requests extra computational work for large orders, and it decreases the size of the computational domain [27]. Meanwhile, the RDTM is simpler than the DTM, and the whole number of calculations needed in the RDTM is less than that in the DTM [28]. The FRDTM is an improved method of the RDTM for fractional-order derivatives. The MsDTM can overcome the difficulties of the DTM and RDTM, where the series solutions often converge in a very insignificant space and the range of convergence is a long procedure or completely divergent in a long time span.

The importance of this research lies in finding the exact and approximate solutions for higher-dimensional time-fractional diffusion equations using a relatively new method and comparing the exact solutions of nonfractional diffusion equations with the approximate solutions for different values of the fractional derivatives. The novel aspect of this research is the explanation of the FRDTM with simple and sequential steps so that every researcher can apply and understand the method directly without referring to other literature. The other distinctive aspect of this research is to find a general formula for the solutions, which reach the exact solution in the first and second examples, and to compare the nonfractional exact solutions with the approximate fractional solutions using figures in two and three dimensions. The central physical purpose of accepting and studying diffusion equations of fractional orders is to define the phenomena of anomalous diffusion, usually met in transport processes through complex and/or disordered media with fractal supports [29]. The time-fractional diffusion Equations (12)–(14) can be interpreted as a model of the diffusion of a particle under the action of the external force.

This paper is organized as follows: After giving simple properties and definitions of the fractional derivative in Section 2, we introduce the proposed method in Section 3. Section 4 presents the exact and approximate solutions of three examples of time-fractional diffusion equations. Section 5 is the conclusion.

2. Preliminaries and Fractional Derivative Order

The unique functions of mathematical physics are found to be very useful for finding solutions of initial- and boundary-value problems governed by partial differential equations and fractional differential equations, and they play a significant and exciting role as solutions of fractional-order differential equations [30]. Many special functions have attracted the attention of researchers, such as the Wright function, the error function, and the Millin–Ross function. In this paper, our attention is focused on only two types of these special functions: the Mittag–Leffler function and the Gamma function. We used the Mittag–Leffler function since after finding the solution in a compact form, we can write the exact solution by using the definition of the Mittag–Leffler function, while the Gamma function is an essential part of the definition of fractional derivatives.

2.1. Mittag–Leffler Function

The Mittag–Leffler (M-L) function is named after a Swedish mathematician who defined and studied it in 1903. The M-L function is a straight generalization of the exponential function ${\mathrm{e}}^x$. The one-parameter M-L function in powers series is given by the formula [3]:

$$E_{\gamma }\left(x\right)=\sum _{k=0}^{\infty }\frac{x^k}{\Gamma (\gamma k+1)},(\gamma >0).$$

(1)

For selected integer values of $\gamma$, we obtain:

$$\begin{array}{ccc}& & E_0\left(x\right)=\frac{1}{1-z},E_1\left(x\right)={\mathrm{e}}^x,\hfill \\ & & E_2\left(x\right)=cosh\left(\sqrt{x}\right).\hfill \end{array}$$

In powers series, the two-parameter M-L function is defined by:

$$E_{\gamma ,\beta }\left(x\right)=\sum _{k=0}^{\infty }\frac{x^k}{\Gamma (\gamma k+\beta )},(\gamma >0,\beta >0).$$

(2)

For special choices of the parameters $\gamma$ and $\beta$, we obtain the famous traditional functions:

$$\begin{array}{ccc}& & E_{1,1}\left(x\right)=E_1\left(x\right)={\mathrm{e}}^x,E_{1,2}\left(x\right)=\frac{e^x-1}{x},\hfill \\ & & E_{2,1}\left(x^2\right)=cosh\left(x\right),E_{2,2}\left(x^2\right)=\frac{sinh\left(x\right)}{x}.\hfill \end{array}$$

2.2. Caputo Fractional Derivative

Let $a\in \mathbb{R}$, then the (left-sided) Caputo fractional derivative $\left({}^c{D}_{a+}^{\alpha }y\right)\left(x\right)$ (the small c represents the Caputo derivative) of order $\alpha \in {\mathbb{R}}^{+}$ is well defined as [31]:

$$\left({}^c{D}_{a+}^{\alpha }y\right)\left(x\right)=\frac{1}{\Gamma (n-\alpha )}{\int }_a^x\frac{y^{\left(n\right)}\left(t\right)}{{(x-t)}^{\alpha -n+1}}\mathrm{d}t,$$

(3)

for $(n-1<\alpha <n;x\ge a),n\in \mathbb{N}$ and $\Gamma \left(x\right)$ is the Gamma function. For the ease of presentation, we symbolize the Caputo fractional derivative as ${\mathbf{D}}_x^{\alpha }f\left(x\right)$.

3. Fractional Reduced Differential Transform Method for n+1 Variables

This section gives the basic definitions and properties of the FRDTM [16,18,32,33]. Consider a function $f(t,x_1,x_2,\dots ,x_n)$ to be analytical and continuously differentiable with respect to $(n+1)$ variables in the domain of interest, such that:

$$f(t,x_1,x_2,\dots ,x_n)=m_1\left(x_1\right)m_2\left(x_2\right)\cdots m_n\left(x_n\right)h\left(t\right).$$

(4)

Then, from the properties of the DTM and motivated by the components of the form $x_1^{i_1}{x}_2^{i_2}\cdots x_n^{i_n}{t}^{\alpha j}$, we write the general solution function $f(t,x_1,x_2,\dots ,x_n)$ as an infinite linear combination of such components:

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill f(t,x_1,x_2,\dots ,x_n)=& \sum _{i_1=0}^{\infty }{m}_1\left(i_1\right)x_1^{i_1}\sum _{i_2=0}^{\infty }{m}_2\left(i_2\right)x_2^{i_2}\cdots \sum _{i_n=0}^{\infty }{m}_n\left(i_n\right)x_n^{i_n}\sum _{j=0}^{\infty }h\left(j\right)t^{\alpha j}\hfill \\ \hfill =& \sum _{i_1=0}^{\infty }\sum _{i_2=0}^{\infty }\cdots \sum _{i_n=0}^{\infty }\sum _{j=0}^{\infty }F(i_1,i_2,\dots ,i_n,j)x_1^{i_1}{x}_2^{i_2}\cdots x_n^{i_n}{t}^{\alpha j},\hfill \end{array}\end{array}$$

(5)

where $F(i_1,i_2,\dots ,i_n,j)=m_1\left(i_1\right)m_2\left(i_2\right)\cdots m_n\left(i_n\right)h\left(j\right)$ is referred to as the spectrum of $f(t,x_1,x_2,\dots ,x_n)$. Furthermore, the lowercase $f(t,x_1,x_2,\dots ,x_n)$ is used for the original function, while its fractional reduced transformed function is represented by the uppercase $F_k(x_1,x_2,\dots ,x_n)$, which is called the T-function.

3.1. Step 1: Finding the Fractional Reduced Transformed Function

Let $f(t,x_1,x_2,\dots ,x_n)$ be analytical and continuously differentiable with respect to $n+1$ variables t, $x_1,x_2,\dots ,x_n$ in the domain of interest, then the FRDTM in n dimensions of $f(t,x_1,x_2,\dots ,x_n)$ is given by:

$$F_k(x_1,x_2,\dots ,x_n)=\frac{1}{\Gamma (k\alpha +1)}{\left[{\mathbf{D}}_t^{\alpha k}\left(f(t,x_1,x_2,\dots ,x_n)\right)\right]}_{t=t_0},$$

(6)

where $k=0,1,2,\cdots$, with the time-fractional derivative.

3.2. Step 2: Finding the Inverse of the Fractional Reduced Transformed Function

The inverse FRDTM of $F_k(x_1,x_2,\dots ,x_n)$ is defined by:

$$f(t,x_1,x_2,\dots ,x_n):=\sum _{k=0}^{\infty }{F}_k(x_1,x_2,\dots ,x_n){(t-t_0)}^{k\alpha }.$$

(7)

From (6) and (7), we have:

$$f(t,x_1,x_2,\dots ,x_n)=\sum _{k=0}^{\infty }\frac{1}{\Gamma (k\alpha +1)}{\left[{\mathbf{D}}_t^{\alpha k}\left(f(t,x_1,x_2,\dots ,x_n)\right)\right]}_{t=t_0}{(t-t_0)}^{k\alpha }.$$

In particular, for $t_0=0$, the above equation becomes:

$$f(t,x_1,x_2,\dots ,x_n)=\sum _{k=0}^{\infty }\frac{1}{\Gamma (k\alpha +1)}{\left[{\mathbf{D}}_t^{\alpha k}\left(f(t,x_1,x_2,\dots ,x_n)\right)\right]}_{t=0}{t}^{k\alpha }.$$

(8)

3.3. Step 3: Finding the Approximate Solution

The inverse transformation of the set of values ${\left\{F_k(x_1,x_2,\dots ,x_n)\right\}}_{k=0}^m$ gives an approximate solution as:

$${\tilde{f}}_m(t,x_1,x_2,\dots ,x_n)=\sum _{k=0}^m{F}_k(x_1,x_2,\dots ,x_n)t^{\alpha k},$$

(9)

where m is the order of the approximate solution.

3.4. Step 4: Finding the Exact Solution

The exact solution using the FRDTM is given by:

$$f(t,x_1,x_2,\dots ,x_n)=\underset{m\to \infty }{lim}{\tilde{f}}_m(t,x_1,x_2,\dots ,x_n).$$

(10)

In

Table 1. Fundamental operations of the FRDTM for $n+1$ variables.

Original Function	Transformed Function
1. $f=c_{1}u\pm c_{2}v$	$F_{\alpha k}=c_1{U}_{\alpha k}\pm c_2{V}_{\alpha k}$
2. $f=uv$	$F_{\alpha k}={\sum }_{i=0}^k{U}_{\alpha i}{V}_{\alpha (k-i)}$
3. $f={\mathbf{D}}_t^{m\alpha }u$	${\displaystyle F_{\alpha k}=\frac{\Gamma \left(\alpha \right(k+m)+1)}{\Gamma (k\alpha +1)}{U}_{\alpha (k+m)}}$
4. ${\displaystyle f=\frac{{\partial }^{h}u}{\partial x_i^h}}$	${\displaystyle F_{\alpha k}=\frac{{\partial }^h{U}_{\alpha k}}{\partial x_i^h}}$, $i=1,2,\dots ,n$
5. $f=x_i^m{t}^r$	$F_{k\alpha }=x_i^m\delta (\alpha k-r)$, $i=1,2,\dots ,n$
6. $f=x_i^m{t}^{r}u$	${\displaystyle F_{\alpha k}=x_i^m\sum _{i=0}^k\delta (\alpha i-r)U_{\alpha (k-r)}}$, $i=1,2,\dots ,n$

, we provide certain properties of the FRDTM, where $\delta (k-m)$ is defined by:

$$\delta (k-m)=\left\{\begin{array}{cc}1,\hfill & k=m\hfill \\ 0,\hfill & k\ne m,\hfill \end{array}\right.$$

(11)

where $f=f(t,x_1,x_2\dots ,x_n),u=u(t,x_1,x_2\dots ,x_n),F_k=F_k(x_1,x_2\dots ,x_n),$ and $U_k=U_k(x_1,x_2\dots ,x_n).$

We prove Property 3 from Table 1 in two dimensions; other proofs of the properties can be found in [33,34,35,36].

If $w(x,y)={\mathbf{D}}_x^{n\alpha }u(x,y)$, then ${\displaystyle W_k\left(y\right)=\frac{\Gamma \left(\alpha \right(k+n)+1)}{\Gamma (k\alpha +1)}U(k+n)\left(y\right)}$.

From Equation (6), we have:

$$\begin{array}{ccc}\hfill W_k\left(y\right)& =& \frac{1}{\Gamma (k\alpha +1)}{\left[{\mathbf{D}}_x^{\alpha k}\left({\mathbf{D}}_x^{n\alpha }u(x,y)\right)\right]}_{x=x_0},\hfill \\ & =& \frac{1}{\Gamma (k\alpha +1)}\left[{\mathbf{D}}_x^{\alpha (k+n)}u(x,y)\right){]}_{x=x_0},\hfill \\ & =& \frac{\Gamma \left(\right(k+n)\alpha +1)}{\Gamma \left(\right(k+n)\alpha +1)\Gamma (k\alpha +1)}\left[{\mathbf{D}}_x^{\alpha (k+n)}u(x,y)\right){]}_{x=x_0},\hfill \\ & =& \frac{\Gamma \left(\right(k+n)\alpha +1)}{\Gamma (k\alpha +1)}U(k+n)\left(y\right).\hfill \end{array}$$

4. Numerical Examples

The purpose of this paper is to apply the FRDTM to find exact and approximate solutions for time-fractional diffusion equations in two, three, and four dimensions. The time-fractional diffusion equation is gained from the standard diffusion equation by consistently changing the first-order time derivative with a specified fractional derivative.

To show the effectiveness of the proposed method for finding exact and approximate solutions, we apply the FRDTM in two-, three-, and four-dimensional time-fractional diffusion equations.

$$\frac{{\partial }^{\beta }f(X,t)}{\partial t^{\beta }}=D\Delta f(X,t)-\nabla ·\mathbf{W}\left(X\right)f(X,t),0<\beta \le 1,D>0,$$

(12)

subject to initial and boundary conditions:

$$f(X,0)=\varphi \left(X\right),X\in \Omega ,$$

(13)

$$f(X,t)=\phi (X,t),X\in \partial \Omega ,t\ge 0.$$

(14)

Here, ${\partial }^{\beta }/\partial t^{\beta }(·)$ is the m-$\beta$-derivative of order $\beta$. $\Delta$ is the Laplace operator. ∇ is the Hamilton operator. $\Omega =[0,L_1]\times [0,L_2]\times \cdots \times [0,L_d]$ is the spatial domain of the problem. d is the dimension of the space, $X=(x_1,x_2,\cdots ,x_d)$. $\partial \Omega$ is the boundary of $\Omega$. $f(X,t)$ denotes the probability density function of finding a particle at X in time t. The positive constant D depends on the temperature, the friction coefficient, the universal gas constant, and lastly, the Avogadro constant. $\mathbf{W}\left(X\right)$ is the external force. Equation (12) can be interpreted as a model of the diffusion of a particle under the action of the external force $\mathbf{W}\left(X\right)$ [37]. In this section, we establish the applicability of the proposed method through test examples.

4.1. Example 1: Two-Dimensional Time-Fractional Diffusion Equations

Let $D=1,\mathbf{W}\left(X\right)=-1,\Omega =[0,1]$, then Equation (12) can be written as:

$$\frac{{\partial }^{\beta }f(x,t)}{\partial t^{\beta }}=\frac{{\partial }^{2}f(x,t)}{\partial x^2}+\frac{\partial f(x,t)}{\partial x},0<\beta \le 1.$$

(15)

The initial condition is given by:

$$f(x,0)=e^x,x\in [0,1].$$

(16)

Applying the appropriate properties given in Table 1 to Equation (15), we obtain the following recurrence relation:

$$F_{k+1}\left(x\right)=\frac{\Gamma (k\beta +1)}{\Gamma \left(\beta \right(k+1)+1)}\left(\frac{{\partial }^{2}f(k-1)}{\partial x\partial x}+\frac{\partial f(k-1)}{\partial x}\right),$$

(17)

for $k=0,1,2,\cdots$. From (17), we find the inverse transform coefficients of $x^{k\alpha }$ as:

$$\begin{array}{ccc}\hfill & & F_0=e^x,\hfill \\ & & F_1=\frac{2e^x}{\Gamma (\beta +1)},\hfill \\ & & F_2=\frac{4e^x}{\Gamma (2\beta +1)},\hfill \\ & & F_3=\frac{8e^x}{\Gamma (3\beta +1)},\cdots ,\hfill \end{array}$$

or in general,

$$F_k=\frac{2^k{e}^x}{\Gamma (1+k\beta )},\mathrm{where}k\ge 0.$$

(18)

After a small number of iterations, the differential inverse transform of ${\left\{F_k\left(y\right)\right\}}_{k=0}^{\infty }$ will provide the resulting series solution:

$$\begin{array}{ccc}\hfill f(x,t)& =& \sum _{k=0}^{\infty }{F}_k\left(x\right)t^{k\beta }\hfill \\ & =& e^x+\frac{2e^x}{\Gamma (\beta +1)}{t}^{\beta }+\frac{4e^x}{\Gamma (2\beta +1)}{t}^{2\beta }\hfill \\ & & +\frac{8e^x}{\Gamma (3\beta +1)}{t}^{3\beta }+\cdots ,\hfill \end{array}$$

which can be written in compact form,

$$f(x,t)=e^x\sum _{k=0}^{\infty }\frac{{\left(2t^{\beta }\right)}^k}{\Gamma (1+k\beta )}.$$

(19)

By the M-L function, we find the exact solution of Equation (15) subject to (16):

$$f(x,t)=e^x{E}_{\beta }\left(2t^{\beta }\right),$$

(20)

where $0<\beta \le 1$ and $E_{\beta }\left(z\right)$ is the one-parameter M-L function (1), and this is exactly the same solution obtained using the FVHPIM with the modified Riemann–Liouville derivative [37]. In comparison with the approximate solution obtained by the HDM via the modified beta derivative equation, the FRDTM gives the direct exact solution with simple computations [9]. In the case of $\beta =1$, we have $E_1\left(2t\right)=e^{2t}$. Then, the exact solution of nonfractional Equation (15) when $\beta =1$ is:

$$f(x,t)=e^{x+2t}.$$

(21)

Figure 1 shows the exact solution of nonfractional order and the three-dimensional plot for the approximate solution by the FRDTM of the fractional order ($\beta =0.8$), while Figure 2 depicts the approximate solutions of the fractional orders ($\beta =0.6,0.4$). In these four figures, it is interesting to note that with a constant range for both variables t and x, the lower the value of the fractional order, the greater the value of the approximate solutions is at the highest value of the variable t. Figure 3 depicts solutions in two-dimensional plots for different values of $\beta$. Through this figure, we can notice that the lower the fractional order, the more the approximate solutions move away from the exact solution of nonfractional order, and their value increases with the value of the variable t being constant. Figure 4 shows solutions in two-dimensional plots for different values of x. In this figure, we notice, with a constant value of the nonfractional order $\beta$, that the greater the values of the variable x, the greater the values of the approximate solutions are with a constant range of the variable t.

4.2. Example 2: Three-Dimensional Time-Fractional Diffusion Equations

Let $D=1,\Omega =[0,1]\times [0,1],\mathbf{W}=-(x,y)$ in Equation (12), then we have the following TFDE:

$$\frac{{\partial }^{\beta }f(x,y,t)}{\partial t^{\beta }}=\frac{{\partial }^{2}f(x,y,t)}{\partial x^2}+\frac{{\partial }^{2}f(x,y,t)}{\partial y^2}+x\frac{\partial f(x,y,t)}{\partial x}+y\frac{\partial f(x,y,t)}{\partial y}+2f(x,y,t),$$

(22)

with initial condition:

$$f(x,y,0)=x+y.$$

(23)

Using the suitable properties from Table 1 for Equation (22), we achieve the next recurrence relation:

$$F_{k+1}(x,y)=\frac{\Gamma (k\beta +1)}{\Gamma \left(\beta \right(k+1)+1)}(\frac{{\partial }^{2}w\left(k\right)}{\partial x\partial x}+\frac{{\partial }^{2}f\left(k\right)}{\partial y\partial y}+x\frac{\partial f\left(k\right)}{\partial x}+y\frac{\partial f\left(k\right)}{\partial y}+2f\left(k\right)),$$

(24)

where $k=0,1,2,\cdots$. The inverse transform coefficients of $t^{k\beta }$ are as follows:

$$\begin{array}{ccc}\hfill & & F_0=x+y,\hfill \\ & & F_1=\frac{3(x+y)}{\Gamma (\beta +1)},\hfill \\ & & F_2=\frac{9(x+y)}{\Gamma (2\beta +1)},\hfill \\ & & F_3=\frac{27(x+y)}{\Gamma (3\beta +1)},\cdots .\hfill \end{array}$$

More generally,

$$U_k=(x+y)\frac{{\left(3\right)}^k}{\Gamma (1+k\beta )}.$$

(25)

Again, if we continue in the same manner, and after a few iterations, the differential inverse transform of ${\left\{F_k(x,y)\right\}}_{k=0}^{\infty }$ will give the following series solution:

$$\begin{array}{ccc}\hfill f(x,y,t)& =& \sum _{k=0}^{\infty }{F}_k(x,y)t^{k\beta }\hfill \\ & =& (x+y)+\frac{3(x+y)}{\Gamma (\beta +1)}{t}^{\beta }+\frac{9(x+y)}{\Gamma (2\beta +1)}{t}^{2\beta }\hfill \\ & & +\frac{27(x+y)}{\Gamma (3\beta +1)}{t}^{3\beta }+\cdots .\hfill \end{array}$$

In compact form,

$$f(x,y,t)=(x+y)\sum _{k=0}^{\infty }\frac{{\left(3t^{\beta }\right)}^k}{\Gamma (1+k\beta )},$$

(26)

and using the M-L function, we obtain the exact solution:

$$f(x,y,t)=(x+y)E_{\beta }\left(3t^{\beta }\right),$$

(27)

where $0<\beta \le 1$ and $E_{\beta }\left(z\right)$ is the one-parameter M-L function (1), which is exactly the same result obtained using the FVHPIM via the m-R-L derivative [37]. In the case of $\beta =1$, $E_1\left(3t\right)=e^{3t}$, the exact solution of the nonfractional Equation (22) is:

$$u(x,y)=(x+y)e^{3t}.$$

(28)

Figure 5 shows the exact solution of nonfractional order and the three-dimensional plot of the approximate solution of the FRDTM ($\beta =0.9$), while Figure 6 depicts the approximate solutions for ($\beta =0.7,0.5$). Figure 7 depicts solutions in two-dimensional plots for different values of $\beta$. Figure 8 shows solutions in two-dimensional plots for different values of x.

4.3. Example 3: Four-Dimensional Time-Fractional Diffusion Equations

Let $D=1$, $\Omega =[0,1]\times [0,1]\times [0,1]$, $\mathbf{F}(x,y,z)=-(x,y,z)$ in Equation (12), then we have the following TFDE:

$$\begin{array}{c}\frac{{\partial }^{\beta }u(x,y,z,t)}{\partial t^{\beta }}=\Delta u(x,y,z,t)+x\frac{\partial u(x,y,z,t)}{\partial x}+\hfill \\ \hfill y\frac{\partial u(x,y,z,t)}{\partial y}+z\frac{\partial u(x,y,z,t)}{\partial z}+3u(x,y,z,t),0<\beta \le 1,\end{array}$$

(29)

with the initial condition,

$$u(x,y,z,0)={(x+y+z)}^2.$$

(30)

Using the appropriate properties from Table 1 for Equation (29), we obtain the following recurrence relation:

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill F_{k+1}(x,y,z)=& \frac{\Gamma (k\beta +1)}{\Gamma \left(\beta \right(k+1)+1)}(\frac{{\partial }^{2}w\left(k\right)}{\partial x\partial x}+\frac{{\partial }^{2}f\left(k\right)}{\partial y\partial y}+\frac{{\partial }^{2}w\left(k\right)}{\partial z\partial z}\hfill \\ \hfill +& x\frac{\partial f\left(k\right)}{\partial x}+y\frac{\partial f\left(k\right)}{\partial y}+z\frac{\partial f\left(k\right)}{\partial z}+3f\left(k\right)),\hfill \end{array}\end{array}$$

(31)

where $k=0,1,2,\cdots$. The inverse transform coefficients of $t^{k\beta }$ are as follows:

$$\begin{array}{ccc}\hfill F_0& =& {(x+y+z)}^2,\hfill \\ \hfill F_1& =& \frac{5{(x+y+z)}^2+6}{\Gamma (\beta +1)},\hfill \\ \hfill F_2& =& \frac{25{(x+y+z)}^2+48}{\Gamma (2\beta +1)},\hfill \\ \hfill F_3& =& \frac{125{(x+y+z)}^2+294}{\Gamma (3\beta +1)},\hfill \\ \hfill F_4& =& \frac{625{(x+y+z)}^2+1632}{\Gamma (4\beta +1)},\hfill \\ \hfill F_5& =& \frac{3125{(x+y+z)}^2+8646}{\Gamma (5\beta +1)},\cdots .\hfill \end{array}$$

If we continue in the same manner, and after a few iterations, the differential inverse transform of ${\left\{F_k(x,y,z)\right\}}_{k=0}^{\infty }$ will give the following series solution:

$$\begin{array}{ccc}\hfill f(x,y,z,t)& =& \sum _{k=0}^{\infty }{F}_k(x,y,z)t^{k\beta }\hfill \\ & =& {(x+y+z)}^2+\frac{5{(x+y+z)}^2+6}{\Gamma (\beta +1)}{t}^{\beta }+\frac{25{(x+y+z)}^2+48}{\Gamma (2\beta +1)}{t}^{2\beta }\hfill \\ & & +\frac{125{(x+y+z)}^2+294}{\Gamma (3\beta +1)}{t}^{3\beta }+\cdots .\hfill \end{array}$$

In the case of $\beta =1$, the tenth-order approximate solution of nonfractional Equation (29) is given by:

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill f_{10}(x,y,Z,t)=& \frac{t^{10}\left(9765625{(x+y+z)}^2+29119728\right)}{3628800}\hfill \\ \hfill +& \frac{t^9\left(1953125{(x+y+z)}^2+5800326\right)}{362880}\hfill \\ \hfill +& \frac{t^8\left(390625{(x+y+z)}^2+1152192\right)}{40320}\hfill \\ \hfill +& \frac{t^7\left(78125{(x+y+z)}^2+227814\right)}{5040}\hfill \\ \hfill +& \frac{1}{720}{t}^6\left(15625{(x+y+z)}^2+44688\right)\hfill \\ \hfill +& \frac{1}{120}{t}^5\left(3125{(x+y+z)}^2+8646\right)\hfill \\ \hfill +& \frac{1}{24}{t}^4\left(625{(x+y+z)}^2+1632\right)\hfill \\ \hfill +& \frac{1}{6}{t}^3\left(125{(x+y+z)}^2+294\right)\hfill \\ \hfill +& \frac{1}{2}{t}^2\left(25{(x+y+z)}^2+48\right)\hfill \\ \hfill +& t\left(5{(x+y+z)}^2+6\right)+{(x+y+z)}^2.\hfill \end{array}\end{array}$$

(32)

Equation (29) has been solved using the FVHPIM via m-R-L derivative [37], and the exact solution is:

$$u(x,y,z,t)=\left((3+{(x+y+z)}^2)\right)E_{\beta }\left(5t^{\beta }\right)-3E_{\beta }\left(3t^{\beta }\right).$$

(33)

Figure 9 shows the exact solution of nonfractional order and the three-dimensional plot of the approximate solution of the FRDTM ($\beta =1$), while Figure 10 depicts the approximate solutions for ($\beta =0.9,0.7$). Figure 11 depicts solutions in two-dimensional plots for different values of $\beta$. Figure 12 shows solutions in two-dimensional plots for different values of x.

5. Conclusions

Finding an exact solution is often considered difficult in most cases. By applying the FRDTM in Section 4.1 and Section 4.2, we were able to find exact solutions in the case of the two- and three-dimensional time-fractional diffusion equations, then we plotted the approximate solutions for different values of the fractional-order $\beta$ in the three- and two-dimensional time-fractional diffusion equation, and we also depicted the approximate solutions for different values of x. An approximate solution in the four-dimensional time-fractional diffusion equation was found in Section 4.3, and we compared it with the exact solution of a nonfractional differential equation, then we plotted the approximate solutions for different values of the fractional-order $\beta$ in three- and two-dimensions. Furthermore, we depicted the approximate solutions for different values of x. The graphical representations of the exact and approximate solutions showed the power of the FRDTM for solving different dimensions of the time-fractional diffusion equation. The computations of this paper were carried out by using the computer package of Mathematica 9.