Advances in Mathematical Physics

Volume 2016 (2016), Article ID 7304659, 7 pages

http://dx.doi.org/10.1155/2016/7304659

## Similarity Solutions for Multiterm Time-Fractional Diffusion Equation

^{1}Mathematics & Engineering Physics Department, Faculty of Engineering, Mansoura University, Mansoura 35516, Egypt^{2}Mathematics & Engineering Physics Department, Faculty of Engineering, Modern University for Technology and Information, Cairo 11585, Egypt

Received 12 January 2016; Revised 19 March 2016; Accepted 29 March 2016

Academic Editor: Xiao-Jun Yang

Copyright © 2016 A. Elsaid et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

Similarity method is employed to solve multiterm time-fractional diffusion equation. The orders of the fractional derivatives belong to the interval and are defined in the Caputo sense. We illustrate how the problem is reduced from a multiterm two-variable fractional partial differential equation to a multiterm ordinary fractional differential equation. Power series solution is obtained for the resulting ordinary problem and the convergence of the series solution is discussed. Based on the obtained results, we propose a definition for a multiterm error function with generalized coefficients.

#### 1. Introduction

Fractional differential equations (FDEs) appear in modeling many problems in the fields of science and engineering [1–3]. These problems are modeled using different types of fractional derivative operators which include Riemann-Liouville definition [1], Caputo definition [4], Riesz definition [5], Riesz-Feller definition [6], and the modified Riemann-Liouville definition recently proposed by Jumarie [7].

To obtain analytic solutions to fractional partial differential equations (FPDEs), two methods have been basically used: the first method is the application of both Laplace and Fourier transforms and the second method is the separation of variables technique [1]. But recently several semianalytic methods have been also utilized to present series solution to FPDEs such as Adomian decomposition method [8, 9], homotopy analysis method [10, 11], homotopy perturbation method [12, 13], variational iteration method [14, 15], and fractional differential transformation method [16, 17].

Fractional diffusion models are formulated using fractional derivative operators to replace regular derivatives. Different forms of fractional diffusion equations have been widely researched. For example, a time-fractional diffusion equation has been explicitly introduced in physics by Nigmatullin [18] to describe diffusion in special types of porous media which exhibit a fractal geometry. Giona et al. [19] also presented a time-fractional diffusion equation that describes relaxation phenomena in complex viscoelastic materials. Wyss [20] considered the time-fractional diffusion and wave equations and obtained the solution in closed form in terms of Fox functions. Gorenflo et al. [21] used the similarity method and the Laplace transform method to obtain the scale invariant solution of the time-fractional diffusion-wave equation in terms of the Wright function. Recently, fractional diffusion equations have been modeled and studied via local fractional derivatives as in the work of Liu et al. [22] and Zhao et al. [23].

Fractional diffusion equations have been also utilized to model anomalous diffusion where a particle plume spreads faster than predicted by classical models and may exhibit significant asymmetry. The anomalous diffusion processes differ from regular diffusion in that the dispersion of particles proceeds faster (superdiffusion) or slower (subdiffision) than for the regular case. Anomalous diffusion is one of the most popular phenomena in the theory of random walks and transport processes [24–26]. Also anomalous diffusion is a vast and exciting topic, in particular for “crowded” biological systems, where the transport of molecules plays a central functional role [27]. Lenzi et al. investigate the solutions for a fractional diffusion equation subjected to boundary conditions which can be connected to adsorption-desorption processes; the analytical solutions were obtained using the Green function approach [28]. Prehl et al. in [29] deal with diffusion equation with a focus on space-fractional diffusion equation.

The large number of applications that are modeled via fractional diffusion equations motivated mathematicians to develop and apply several techniques to obtain analytic and approximate solutions to this class of FPDEs. Some examples of this work include studying the existence and uniqueness of solution [30], presenting the fundamental solution to some types of fractional diffusion equations [31, 32], applying Green’s function approach [33], obtaining approximate solution via pseudospectral scheme [34], solution via radial basis functions [35], stochastic solution via Monte Carlo simulation [36], approximate solution based on the shifted Legendre-tau technique [37], solution by integral transform method [38], numerical solution by finite difference scheme [39], approximate analytical solution by homotopy analysis method [40], approximate analytical solution by optimal homotopy analysis method [41], approximate analytical solution via Adomian method [42], approximate analytical solution by variational iteration method [43], and approximate analytical solution via generalized differential transform method [44].

Multiterm FPDEs provide more generalizations to fractional order models. Thus many authors have dealt with them either analytically or numerically. For example, the authors of [45] developed a direct solution technique for solving the linear multiterm FDEs with constant coefficients using a spectral tau method. In [46], an effective finite element method is presented for the multiterm time-space Riesz fractional advection-diffusion equations. A numerical method is proposed in [47] to reduce the multiterm FDEs to a system of algebraic equations based on Bernstein polynomials basis. Also in [48], the authors present a numerical method for solving the multiterm time-fractional wave-diffusion equations. A study on the nonhomogeneous generalized multiterm fractional heat propagation and fractional diffusion-convection equation in three-dimensional space is presented in [49].

Though symmetry methods have been applied to different types of linear and nonlinear partial differential equations, the research on using these methods for solving FPDEs is still in the initial stage. The work reported on the Riemann-Liouville fractional derivative includes [50] where scaling transformations are derived for the time-fractional heat equation to reduce it to an ordinary FDE but with Erdelyi-Kober fractional differential operator. Also, similarity solutions for the time-fractional nonlinear conduction equations are presented in [51] to reduce the considered problems to ordinary FDEs that are solved by analytic and numerical techniques. The fractional derivative proposed by Jumarie is studied in [52] where the Lie group method is applied to a space-time-fractional diffusion equation. Choksi and Timol [53] derived similarity solution for fractional diffusion-wave equation using Lie group transformation and Erdelyi-Kober fractional differential operator. Through the invariants of the group of scaling transformations, Duan et al. [54] derived the integrodifferential equation for the similarity variable; then, by virtue of Mellin transform, the fundamental solution of the FDE was expressed in terms of Fox functions. Finally, El Kinani and Ouhadan [55] used Lie symmetry analysis to reduce the number of independent variables of time-fractional partial differential equations; then they employed symmetry properties to construct some exact solutions.

In this work, we use similarity methods to solve multiterm time-fractional diffusion equation with fractional derivatives defined in Caputo sense. The FPDE is reduced to an ordinary FDE which is solved by the power series method. This paper is arranged as follows. The definition and properties of Caputo fractional derivative are listed in Section 2. In Section 3, we illustrate how the similarity method technique is employed to transform a multiterm FPDE into a multiterm FDE. In Section 4, power series method is applied to solve the obtained FDE. The conclusion of this work is summarized in Section 5.

#### 2. Fractional Calculus

*Definition 1. *A real function , , is said to be in the space , , if there exists a real number , such that , where , and it is said to be in the space if , .

*Definition 2. *The Riemann-Liouville fractional integral operator of order of a function , is defined as

*Definition 3. *The fractional derivative in Caputo sense of , , is defined as [1]

One basic property of Caputo fractional derivative is

For more details on Caputo fractional derivative definition and properties see [1, 56, 57].

#### 3. Similarity Method Solution

In this section, we illustrate the technique for using similarity methods in solving multiterm FPDEs with Caputo fractional derivative.

*Problem 4. *Consider the following problem:

where are constants, with an initial condition . Here denotes the partial fractional derivative of order of with respect to the time variable in the Caputo sense. To solve (4), first we perform its scaling transformation using similarity methods; see [51, 58]. Consider the new independent and dependent variables denoted by , , and defined in the following way:where is called the scaling parameter and , , and are arbitrary constants to be determined such that (4) remains invariant under this transformation. From Caputo definition for , one may easily verify thatwhere . Also, for , we have

Hence, by substituting (6) and (7) into (4), we get

From (8), it is clear that, by setting and , then (4) is invariant under transformation (5). The characteristic equation associated with transformation (5) is given by

At , this shows that can be expressed aswhere

By using formula (10), we havewhere .

Again, by using formula (10), we have

Substituting (11) and (12) into (4), the resulting FDE is given bywith the initial condition .

#### 4. Power Series Solution

In this section, we develop the series solution of (13). Let

Using (3), we have

Substituting (15) and (16) into (13) and comparing the coefficients of identical powers on both sides, we get

For the initial condition , we get . Hence, takes the formBased on these results, we propose a definition for multiterm error function with generalized coefficients and discuss its convergence.

*Definition 5. *A multiterm error function with generalized coefficients is defined in the formwhere

By this definition, the solution of the ordinary FDE obtained in (13) is given bywhere and are arbitrary constants.

To check the radius of convergence of the multiterm error function with generalized coefficients (19) and (20), we evaluate the following limit:

Hence, the series solution converges for all

We consider the case in (18). Figures 1, 2, and 3 illustrate the effect of changing the orders of fractional derivatives , , and on the behavior of Taylor series representation of the solution function at different values of the constants , .