Mathematical Problems in Engineering

Volume 2013 (2013), Article ID 186934, 7 pages

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

Research Article

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

Zhoujin Cui,¹ Zisen Mao,² Sujuan Yang,² and Pinneng Yu²

¹Jiangsu Maritime Institute, Nanjing 211100, China

²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.

(a)

(b)

Figure 1: The comparison of the results of the RDTM and HPM for and at and .

(a)

(b)

(c)

(d)

Figure 2:   Results of the HPM for , , and and .

Case Study 2. . See Figures 3 and 4.

(a)

(b)

Figure 3: The comparison of the results of the RDTM and HPM for and at and .  

(a)

(b)

(c)

(d)

Figure 4: Results of the HPM for , , and and .

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.

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.

S. G. Samko, A. A. Kilbas, and O. I. Marichev, Fractional Integrals and Derivatives. Theory and Applications, Gordon and Breach, London, UK, 1993.

V. Kiryakova, Generalized Fractional Calculus and Applications, vol. 301 of Pitman Research Notes in Mathematics Series, Longman Scientific & Technical, Harlow, UK, 1994.

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.

I. Podlubny, Fractional Differential Equations, vol. 198 of Mathematics in Science and Engineering, Academic Press, San Diego, Calif, USA, 1999.

R. Hilfer, Applications of Fractional Calculus in Physics, World Scientific, Singapore, 2000.

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.

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.

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.

W. Zimmermann, “Propagating fronts near a Lifshitz point,” Physical Review Letters, vol. 66, no. 11, p. 1546, 1991.

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.

G. Caginalp and P. Fife, “Higher-order phase field models and detailed anisotropy,” Physical Review B, vol. 34, no. 7, pp. 4940–4943, 1986.

J.-H. He, “Homotopy perturbation technique,” Computer Methods in Applied Mechanics and Engineering, vol. 178, no. 3-4, pp. 257–262, 1999.

J. H. He, “A new perturbation technique which is also valid for large parameters,” Chaos, Solitons and Fractals, vol. 26, pp. 827–833, 2005.

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.

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.

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.

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.

S. Das and P. K. Gupta, “A mathematical model on fractional Lotka-Volterra equations,” Journal of Theoretical Biology, vol. 277, pp. 1–6, 2011.

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.

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.

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.

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.

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.

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.

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.

Y. Keskin and G. Oturanç, “Reduced differential transform method for generalized KdV equations,” Mathematical & Computational Applications, vol. 15, no. 3, pp. 382–393, 2010.