Mathematical Problems in Engineering

Volume 2013, Article ID 186934, 7 pages

http://dx.doi.org/10.1155/2013/186934

## Approximate Analytical Solutions of Fractional Perturbed Diffusion Equation by Reduced Differential Transform Method and the Homotopy Perturbation Method

^{1}Jiangsu Maritime Institute, Nanjing 211100, China^{2}Institute of Science, PLA University of Science and Technology, Nanjing 211101, China

Received 23 March 2013; Accepted 6 May 2013

Academic Editor: Guanghui Wen

Copyright © 2013 Zhoujin Cui 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

The approximate analytical solutions of differential equations with fractional time derivative are obtained with the help of a general framework of the reduced differential transform method (RDTM) and the homotopy perturbation method (HPM). RDTM technique does not require any discretization, linearization, or small perturbations and therefore it reduces significantly the numerical computation. Comparing the methodology (RDTM) with some known technique (HPM) shows that the present approach is effective and powerful. The numerical calculations are carried out when the initial conditions in the form of periodic functions and the results are depicted through graphs. The two different cases have studied and proved that the method is extremely effective due to its simplistic approach and performance.

#### 1. Introduction

Various fields of science and engineering deal with dynamical systems (see [1, 2]), some of which can be described by fractional-order equations (see [3–5] and the references therein). The last two decades have witnessed a great progress in fractional calculus and fractional-order dynamical systems. It has been found that fractional calculus is a mathematical tool that works adequately for anomalous social and physical systems with nonlocal, frequency- and history-dependent properties, and for intermediate states such as soft materials, which are neither idea solid nor idea fluid (see [3, 4, 6]). Differential equations with fractional-order derivatives/integrals are called fractional differential equations, and they have found many successful applications in viscoelasticity, heat conduction, electromagnetic wave, diffusion wave, control theory, and so on (see [7–9] and the references therein).

Bistable systems play an important role in the study of spatial patterns. A typical example, which appears in population dynamics, leads to referred to as the classical Fisher-Kolmogorov (FK) equation. In investigating such spatial patterns, a key role is played by a model equation, which is simpler than the full equation describing the process. Recently, interest has turned to higher-order model equations involving bistable dynamics, such as the extended Fisher-Kolmogorov (EFK) equation proposed by Coullet et al. [10] in 1987 as well as by Dee and Van Saarloos [11] in 1988. The EFK equation has appeared in studies of phase transitions, for instance, near a Lifshitz point (cf. Zimmerman [12]). Another well-known equation of this type is the Swift-Hohenberg equation proposed in 1977.

Whereas the EFK equation and the SH equation involve fourth-order spatial derivatives, certain phase-field models lead to even-higher-order spatial gradients. The following equation studied by Gardner and Jones [13] and Caginalp and Fife [14], in which and are constants. In this paper, we consider the initial value problem (IVP) of perturbed diffusion equation with fractional time derivative where the fractional derivative in (5) is in the sense of Caputo and .

Historically, comparatively little was known about the extraordinary range of behavior presented by the solutions of nonlinear partial differential equations. Except for a few special cases, it is impossible to find a closed form solution for a fractional differential equation. Effective methods for solving such equations are needed. So approximate and numerical techniques must be used.

The homotopy perturbation method (HPM) is relatively new approach to provide an analytical approximation to linear or nonlinear problem. This method was first presented by He [15, 16] and applied to various nonlinear problems [17, 18]. The basic difference of this method from the other perturbation techniques is that it does not require small parameters in the equation, which overcomes the limitations of the traditional perturbation techniques. Recently, the application of this method is extended for fractional differential equations [5, 19–21]. Zhou [22] was the first one to use differential transform method (DTM) in engineering applications. He employed DTM in solution of initial boundary value problems in electric circuit analysis. This method constructs an analytical solution in the form of a polynomial. It is different from the traditional high order Taylor series method, which requires symbolic computation of necessary derivatives of the data functions. The differential transform is an iterative procedure for obtaining analytical differential equations. The concept of DTM has broadened to problems involving partial differential equations and systems of differential equations [23, 24]. Some researchers have lately applied DTM for analysis of uniform and non-uniform beams [25–27]. Recently, the application of reduced differential transform method is successfully extended to obtain analytical approximate solutions to linear and nonlinear ordinary differential equations of fractional order [28–30].

In this paper, we try to find an approximate analytical solution of (5) for fractional time derivatives with the help of powerful analytical method. We use the reduced differential transform method (RDTM) and homotopy perturbation method (HPM) to obtain the solutions and compare them with each other. We know that the HPM method is based on the use of homotopy parameter for classification of most favorable values of parameters in between . While, RDTM technique does not require any parameter, discretization, linearization, or small perturbations and therefore it reduces significantly the numerical computation. For the standard cases, comparing the methodology with some known techniques shows that the present approach is effective and powerful.

#### 2. Preliminaries

There are several approaches to define the fractional calculus, for example, Riemann-Liouville, Gruunwald-Letnikow, Caputo, and Generalized Functions approach. For the readers’ convenience, definitions of fractional integral/derivative and some preliminary results are given in this section.

*Definition 1 (see [6]). *The fractional integral of order of a function is given by
provided that the right side is point-wise defined on , where is the well-known gamma function.

*Definition 2 (see [6]). *The Caputo derivative of order of a continuous function is defined to be
where , provided that the right side is point-wise defined on .

*Definition 3 (see [6]). *The Riemann-Liouville derivative of order of a continuous function is defined to be
where , provided that the right side is point-wise defined on .

Riemann-Liouville fractional derivative is mostly used by mathematicians but this approach is not suitable for real world physical problems since it requires the definition of fractional order initial conditions, which have no physically meaningful explanation yet. Caputo introduced an alternative definition, which has the advantage of defining integer order initial conditions for fractional order differential equations. We have chosen the Caputo fractional derivative because it allows traditional initial and boundary conditions to be included in the formulation of the problem. And some other properties of fractional derivative can be found in [4, 6].

The relation between the Riemann-Liouville operator and Caputo operator is given by By using (9), from (8), we obtain fractional derivative in the Caputo sense as follows:

Let us expand the analytical and continuous function in terms of a fractional power series as follows: where is the order of fraction and is the fractional differential transform of .

Since the initial conditions are implemented to the integer order derivatives, the transformations of the initial conditions are defined as follows:

The following theorems that can be deduced from (7) and (8) are given below.

Theorem 4. *If , then .*

Theorem 5. *If , then
*

Theorem 6. *If , then .*

Theorem 7. *If , then .*

#### 3. Solution of the Problem by the RDTM and HPM

To illustrate the basic ideas of HPM [15], we consider the following general nonlinear differential equation: with boundary conditions where is a general differential operator, is a boundary operator, is a known analytic function, and is the boundary of the domain .

The operator can be generally divided into two parts and , where is linear, while is nonlinear. Therefore (14) can be written as follows: We construct a homotopy of (14) which satisfies which is equivalent to where is an embedding parameter, and is an initial guess approximation of (14) which satisfies the boundary conditions. It follows from (17) and (18) that Thus, the changing process of from to is just that of from to . In topology this is called deformation and and are called homotopic. Here the embedding parameter is introduced much more naturally, unaffected by artificial factors; further it can be considered as a small parameter for . So it is very natural to assume that the solution of (18) and (19) can be expressed as The approximate solutions of the original equations can be obtained by setting ; that is, The convergence of series (21) has been proved by He in his paper [15].

Because of the knowledge of various perturbation methods that low-order approximate solution leads to high accuracy, there requires no infinite series. Then after a series of recurrent calculation by using Mathematica software, we will get approximate solutions of fractional model (5).

According to HPM, we construct the following homotopy of (5) where the homotopy parameter is considered to be small, .

Assuming the solution of (22) has the form when , (23) becomes the approximate solution of (5). Substituting (23) into (22) and equating the terms with identical powers of , we obtain the following set of linear differential equations: and so on. The method is based on applying the operators (the inverse operators of the Caputo derivative ) on both sides of the above linear differential equations. Using this selection we obtain the successive approximations and the solution may be obtained.

Finally, we approximate the analytical solutions of by the truncated series

For the illustration of the methodology of the reduced differential transform method, we write (5) in the standard operator form with initial condition where is the fractional time derivative operator and is linear operator and is nonlinear operator.

According to the RDTM and Theorems 6 and 7, we can construct the following iteration formula for (26): For the convenience of readers, we can give the first few nonlinear term are From initial condition, we write Substituting (30) into (28) and by a straight forward iterative calculations, we get the following values: and so on. Then the inverse transformation of the set of values gives approximation solution as where is order of approximation solution.

Therefore, the exact solution of problem is given by

#### 4. Numerical Results and Discussion

Here, we take three different values of and comparing the results of RDTM and HPM in the form of three- and two dimensional figures for each case, we would see that RDTM and HPM solutions are in excellent agreement. In this section we assume that .

*Case Study 1. *. See Figures 1 and 2.

*Case Study 2*. . See Figures 3 and 4.

#### 5. Conclusions

In this paper, RDTM and HPM are successfully applied to find the solution of the fractional differential equation and standard motion with different initial conditions. Unlike the traditional methods, the solutions here are given in series 9 form. The approximate solution to the equation was computed without any need for special transformations, linearization, or discretization. In addition, we compare these two methods and show that the results of the RDTM method are in excellent agreement with results of the HPM method and the obtained numerical solutions are shown graphically. We use the mathematical simulation to calculate the functions obtained from the RDTM and HPM. It was shown that RDTM and HPM methods are powerful tools for solving analytically of nonlinear equations.

#### Acknowledgments

This project was supported by the National Natural Science Foundation of China (no. 41206163), Open Research Program of State Key Laboratory of Tropical Oceanography (South China Sea Institute of Oceanology Chinese Academy of Sciences) (no. LTO1205), and the Pre-research Foundation of PLA University of Science and Technology.

#### References

- G. Wen and Z. Duan, “Dynamics behaviors of weighted local-world evolving networks with extended links,”
*International Journal of Modern Physics C*, vol. 20, no. 11, pp. 1719–1735, 2009. View at Publisher · View at Google Scholar · View at Scopus - G. Wen, Z. Duan, W. Yu, and G. Chen, “Consensus of second-order multi-agent systems with delayed nonlinear dynamics and intermittent communications,”
*International Journal of Control*, vol. 86, no. 2, pp. 322–331, 2013. View at Google Scholar - S. G. Samko, A. A. Kilbas, and O. I. Marichev,
*Fractional Integrals and Derivatives. Theory and Applications*, Gordon and Breach, London, UK, 1993. View at Zentralblatt MATH · View at MathSciNet - V. Kiryakova,
*Generalized Fractional Calculus and Applications*, vol. 301 of*Pitman Research Notes in Mathematics Series*, Longman Scientific & Technical, Harlow, UK, 1994. View at Zentralblatt MATH · View at MathSciNet - Z. Cui, Z. Yang, and Z. Rui, “Application of homotopy perturbation method to nonlinear fractional population dynamics models,”
*International Journal of Applied Mathematics and Computation*, vol. 4, pp. 403–412, 2012. View at Google Scholar - I. Podlubny,
*Fractional Differential Equations*, vol. 198 of*Mathematics in Science and Engineering*, Academic Press, San Diego, Calif, USA, 1999. View at Zentralblatt MATH · View at MathSciNet - R. Hilfer,
*Applications of Fractional Calculus in Physics*, World Scientific, Singapore, 2000. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. Xu and W. Tan, “Intermediate processes and critical phenomena: theory, method and progress of fractional operators and their applications to modern mechanics,”
*Science in China G*, vol. 49, no. 3, pp. 257–272, 2006. View at Publisher · View at Google Scholar · View at Scopus - Z. Cui, P. Yu, and W. Zong, “Dynamical behaviors and chaos in a new fractional-order financial system,” in
*Proceedings of the 5th International Workshop on Chaos-Fractals Theories and Applications (IWCFTA '12)*, pp. 109–113, 2012. - P. Coullet, C. Elphick, and D. Repaux, “Nature of spatial chaos,”
*Physical Review Letters*, vol. 58, no. 5, pp. 431–434, 1987. View at Publisher · View at Google Scholar · View at MathSciNet - G. T. Dee and W. Van Saarloos, “Bistable systems with propagating fronts leading to pattern formation,”
*Physical Review Letters*, vol. 60, no. 25, pp. 2641–2644, 1988. View at Publisher · View at Google Scholar · View at Scopus - W. Zimmermann, “Propagating fronts near a Lifshitz point,”
*Physical Review Letters*, vol. 66, no. 11, p. 1546, 1991. View at Publisher · View at Google Scholar · View at Scopus - R. A. Gardner and C. K. R. T. Jones, “Traveling waves of a perturbed diffusion equation arising in a phase field model,”
*Indiana University Mathematics Journal*, vol. 39, no. 4, pp. 1197–1222, 1990. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - G. Caginalp and P. Fife, “Higher-order phase field models and detailed anisotropy,”
*Physical Review B*, vol. 34, no. 7, pp. 4940–4943, 1986. View at Publisher · View at Google Scholar · View at MathSciNet - J.-H. He, “Homotopy perturbation technique,”
*Computer Methods in Applied Mechanics and Engineering*, vol. 178, no. 3-4, pp. 257–262, 1999. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - J. H. He, “A new perturbation technique which is also valid for large parameters,”
*Chaos, Solitons and Fractals*, vol. 26, pp. 827–833, 2005. View at Google Scholar - J. H. He, “Homotopy perturbation method for bifurcation of nonlinear problems,”
*International Journal of Nonlinear Sciences and Numerical Simulation*, vol. 6, no. 2, pp. 207–208, 2005. View at Google Scholar · View at Scopus - X. Li, M. Xu, and X. Jiang, “Homotopy perturbation method to time-fractional diffusion equation with a moving boundary condition,”
*Applied Mathematics and Computation*, vol. 208, no. 2, pp. 434–439, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - S. Momani and Z. Odibat, “Homotopy perturbation method for nonlinear partial differential equations of fractional order,”
*Physics Letters A*, vol. 365, no. 5-6, pp. 345–350, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - S. Das, P. K. Gupta, and Rajeev, “A fractional predator-prey model and its solution,”
*International Journal of Nonlinear Sciences and Numerical Simulation*, vol. 10, no. 7, pp. 873–876, 2009. View at Google Scholar · View at Scopus - S. Das and P. K. Gupta, “A mathematical model on fractional Lotka-Volterra equations,”
*Journal of Theoretical Biology*, vol. 277, pp. 1–6, 2011. View at Publisher · View at Google Scholar · View at MathSciNet - J. K. Zhou,
*Differential Transformation and Its Applications for Electrical CircuIts*, Huazhong University Press, Wuhan, China, 1986. - M.-J. Jang, C.-L. Chen, and Y.-C. Liu, “Two-dimensional differential transform for partial differential equations,”
*Applied Mathematics and Computation*, vol. 121, no. 2-3, pp. 261–270, 2001. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - F. Ayaz, “Solutions of the system of differential equations by differential transform method,”
*Applied Mathematics and Computation*, vol. 147, no. 2, pp. 547–567, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - S. H. Ho and C. K. Chen, “Analysis of general elastically end restrained non-uniform beams using differential transform,”
*Applied Mathematical Modelling*, vol. 22, no. 4-5, pp. 219–234, 1998. View at Google Scholar · View at Scopus - O. O. Ozdemir and M. O. Kaya, “Flapwise bending vibration analysis of a rotating tapered cantilever Bernoulli-Euler beam by differential transform method,”
*Journal of Sound and Vibration*, vol. 289, no. 1-2, pp. 413–420, 2006. View at Publisher · View at Google Scholar · View at Scopus - S. Çatal, “Solution of free vibration equations of beam on elastic soil by using differential transform method,”
*Applied Mathematical Modelling*, vol. 32, no. 9, pp. 1744–1757, 2008. View at Publisher · View at Google Scholar · View at Scopus - A. Arikoglu and I. Ozkol, “Solution of fractional differential equations by using differential transform method,”
*Chaos, Solitons and Fractals*, vol. 34, no. 5, pp. 1473–1481, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Y. Keskin and G. Oturanç, “Reduced differential transform method for partial differential equations,”
*International Journal of Nonlinear Sciences and Numerical Simulation*, vol. 10, no. 6, pp. 741–749, 2009. View at Google Scholar · View at Scopus - Y. Keskin and G. Oturanç, “Reduced differential transform method for generalized KdV equations,”
*Mathematical & Computational Applications*, vol. 15, no. 3, pp. 382–393, 2010. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet