Research Article  Open Access
R. Yulita Molliq, M. S. M. Noorani, "Solving the Fractional RosenauHyman Equation via Variational Iteration Method and Homotopy Perturbation Method", International Journal of Differential Equations, vol. 2012, Article ID 472030, 14 pages, 2012. https://doi.org/10.1155/2012/472030
Solving the Fractional RosenauHyman Equation via Variational Iteration Method and Homotopy Perturbation Method
Abstract
In this study, fractional RosenauHynam equations is considered. We implement relatively new analytical techniques, the variational iteration method and the homotopy perturbation method, for solving this equation. The fractional derivatives are described in the Caputo sense. The two methods in applied mathematics can be used as alternative methods for obtaining analytic and approximate solutions for fractional RosenauHynam equations. In these schemes, the solution takes the form of a convergent series with easily computable components. The present methods perform extremely well in terms of efficiency and simplicity.
1. Introduction
Recent advances of fractional differential equations are stimulated by new examples of applications in fluid mechanics, viscoelasticity, mathematical biology, electrochemistry, and physics. For example, the nonlinear oscillation of earthquake can be modeled with fractional derivatives [1], and the fluiddynamic traffic model with fractional derivatives [2] can eliminate the deficiency arising from the assumption of continuum traffic flow. Based on experimental data fractional partial differential equations for seepage flow in porous media are suggested in [3], and differential equations with fractional order have recently proved to be valuable tools to the modeling of many physical phenomena [4]. Fractional partial differential equations also have studied and successfully solved such as the spacetime fractional diffusionwave equation [5–7], the fractional advectiondispersion equation [8, 9], the fractional KdV equation [10], and the linear inhomogeneous fractional partial differential equations [11].
Most nonlinear differential equations are usually arising from mathematical modeling of many physical systems. In most cases, it is very difficult to achieve analytic solutions of these equations. Perturbation techniques are widely used in science and engineering to handle nonlinear problems and do great contribution to help us understand many nonlinear phenomena. However, perturbation techniques are based on the existence of small/large parameter. Therefore, these techniques are not valid for strongly nonlinear problems.
The homotopy perturbation method (HPM) is the new approach for finding the approximate analytical solution of linear and nonlinear problems. The method was first proposed by He [12, 13] and was successfully applied to solve nonlinear wave equation by He [14–16]. The convergence of Homotopy perturbation series to the exact solution is considered in [17]. Similarly, applying the variational iteration method, created by He [3, 18, 19], consists in constructing the appropriate correction functional connected with the considered equation. The correction functional contains a Lagrange multiplier, the determination of which leads to a recurrence formula. Convergence of the VIM method is discussed by Tatari and Dehghan in [20]. Both of the methods examined have found application in determining the approximate solutions of different technical problems [21]. Adaptation of the VIM method for solving fractional heatwavelike equation and fractional ZakharovKuznetsov equation were discussed by Yulita and colleagues in [22, 23]. Whereas, Chun [24] obtained the numerical solution of heat conduction problem by VIM. Recently, the application of the VIM method for solving kuramoto and Sivashinsky equations was presented by Porshokouhi and Ghanbari in [25]. For the application of HPM, this method used for solving fractional vibration equation [26] and partial differential equations of fractional order in finite domains [27].
In the present paper, VIM and HPM will be applied for solving fractional RosenauHay nam equation which written as subject to the initial condition where , is a parameter describing the order of the fractional derivative (), is the time, and is the spatial coordinate. Fractional RH equation when has appeared in the study of the formation of patterns in liquid drops [28].
2. Basic Definitions
Fractional calculus unifies and generalizes the notions of integerorder differentiation and fold integration [4, 29]. We give some basic definitions and properties of fractional calculus theory which will be used in this paper:
Definition 2.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 and only if , .
The RiemannLiouville fractional integral operator is defined as follows.
Definition 2.2. The RiemannLiouville fractional integral operator of order , of a function , , is defined as
In this paper only real and positive values of will be considered.
Properties of the operator can be found in [29] and we mention only the following: for , , , and : (1), (2), (3). The ReimannLiouville derivative has certain disadvantages when trying to model realworld phenomena with FDEs. Therefore, we will introduce a modified fractional differential operator proposed by Caputo in his work on the theory of viscoelasticity [30].
Definition 2.3. The fractional derivative of in Caputo sense is defined as In addition, we also need the following property.
Lemma 2.4. If , and , , then The Caputo differential derivative is considered here because the initial and boundary conditions can be included in the formulation of the problems [4]. The fractional derivative is taken in the Caputo sense as follows.
Definition 2.5. For to be the smallest integer that exceeds , the Caputo fractional derivative operator of order is defined as For more information on the mathematical properties of fractional derivatives and integrals, one can consult [4, 29].
3. Basic Idea of He’s Variational Iteration Method
To clarify the basic ideas of VIM, we consider the following differential equation: where , is the Caputo fractional derivative of order , , is a nonlinear function, and is the source function. According to VIM, we can write down a correction functional as follows: where is a general Lagrangian multiplier which can be optimally identified via the variational theory [31]. is considered as a restricted variation [32], that is, and the subscript indicates the th approximation. We have where is considered as restricted variations, that is, . For , we have Thus, we obtain the following stationary conditions: Solving this system of equations yields Furthermore, substituting (3.6) to (3.2), the iteration formula of VIM can be written as follows: In this case, we begin with the initial approximation: The correction functional (3.8) will give several approximations, and therefore the exact solution is obtained as
4. Basic Idea of Homotopy Perturbation Method
Consider the following 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, generally speaking, be divided into two parts, and , where is linear and is nonlinear, therefore (4.1) can be written as By using homotopy technique, one can construct a homotopy which satisfies or where is an embedding parameter, and is the initial approximation of (4.1) which satisfies the boundary conditions. Clearly, we have or the changing process of from zero to unity is just that of changing from to . This is called deformation, and also, and are called homotopic in topology. If the embedding parameter is considered as a small parameter, applying the classical perturbation technique, we can assume that the solution of (4.3) and (4.4) can be given as a power series in , that is, and setting results in the approximate solution of (4.1) as
5. Application of HPM and VIM Methods
We consider the application of VIM to fractional RosenauHynam equations which is rewritten as follows: where with the initial conditions of
5.1. VIM Implement for Fractional RosenauHynam Equation
According to the formula (3.8), the iteration formula for (4.9) is given by The iteration starts with an initial approximation which is initial condition in (5.2). Furthermore, using the iteration formula in (5.3), we can directly obtain other components as and so on.
The exact solution of this equation is given by [33] where is arbitrary constant [28].
6. Analysis of the Homotopy Perturbation Method (HPM)
Now applying the classical perturbation technique for solving (5.1) with initial condition in (5.2). To solve (5.1) by the homotopy perturbation method, we construct the following homotopy: or where is an embedding parameter. If , then (6.2) becomes a linear equation, And when , then (6.2) turns out to be (5.1). Assume the solution of (6.2) to be in the form Substituting (6.4) into (6.2) and equating the terms with identical powers of , we obtain the following set of linear differential equations: and so on. Equations (6.6) and (6.7) can be solved by applying the operator , which is the inverse of the operator and then by simple computation, Thus, the solution reads as follows and so on.
In this manner, the rest of the components of the homotopy perturbation series can be obtained. Finally, we approximate the analytical solutions of by the truncated series where .
7. Numerical Results and Discussion
Tables 1 and 2 show the approximate solutions for (5.1) obtained for different values of using the decomposition method and the variational iteration method in different values of , that is, and , respectively. Tables 3 and 4 show the absolute error of (5.1) when in different value of , that is, and , respectively. Figures 1 and 2 show the approximate solutions for (5.1) in different values of using the fifth iterates of VIM, the fifth terms of HPM when , and exact solution, respectively. From Tables 3 and 4 show that the approximate solution using the VIM is more accurate than the approximate solution obtained using the HPM. It is to be noted that only the fifth iterates of the variational iteration solution and only fifth terms of the homotopy perturbation series were used in evaluating.




(a)
(b)
(c)
(a)
(b)
(c)
8. Conclusions
The fundamental goal of this work has been to construct an approximate solution of nonlinear partial differential equations of fractional order. For computations and plots, the Mathematica and Maple packages were used. The goal has been achieved by using the variational iteration method (VIM) and the homotopy perturbation method (HPM). The methods were used in a direct way without using linearization or restrictive assumptions. There are four important points that were gotten. First, the VIM and the HPM provide the solutions in terms of convergent series with easily computable components. Second, the approximate solution in (5.1) using the VIM converges faster than the approximate solution using the HPM. Third, the variational iteration method handles nonlinear equations without any need for the socalled He’s polynomials. Finally, the recent appearance of fractional differential equations as models in some fields of applied mathematics makes it necessary to investigate methods of solution for such equations (analytical and numerical) and we hope that this work is a step in this direction.
Acknowledgment
The financial support received from UKM Grant ERGS/1/2011/STG/UKM/01/13 is gratefully acknowledged.
References
 J. H. He, “Semiinverse method of establishing generalized variational principles for fluid mechanics with emphasis on turbomachinery aerodynamics,” International Journal of Turbo and Jet Engines, vol. 14, no. 1, pp. 23–28, 1997. View at: Google Scholar
 J. H. He, “Some applications of nonlinear fractional differential equations and their approximations,” Bulletin of Science, Technology & Society, vol. 15, no. 2, pp. 86–90, 1999. View at: Google Scholar
 J. H. He, “Approximate analytical solution for seepage flow with fractional derivatives in porous media,” Computer Methods in Applied Mechanics and Engineering, vol. 167, no. 12, pp. 57–68, 1998. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 I. Podlubny, Fractional Differential Equations, vol. 198 of Mathematics in Science and Engineering, Academic Press, San Diego, Calif, USA, 1999.
 K. AlKhaled and S. Momani, “An approximate solution for a fractional diffusionwave equation using the decomposition method,” Applied Mathematics and Computation, vol. 165, no. 2, pp. 473–483, 2005. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 F. Mainardi, Y. Luchko, and G. Pagnini, “The fundamental solution of the spacetime fractional diffusion equation,” Fractional Calculus & Applied Analysis, vol. 4, no. 2, pp. 153–192, 2001. View at: Google Scholar  Zentralblatt MATH
 A. Hanyga, “Multidimensional solutions of timefractional diffusionwave equations,” Proceedings of the Royal Society of London, Series A, vol. 458, no. 2020, pp. 933–957, 2002. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 F. Huang and F. Liu, “The time fractional diffusion equation and the advectiondispersion equation,” The Australian & New Zealand Industrial and Applied Mathematics Journal, vol. 46, no. 3, pp. 317–330, 2005. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 F. Huang and F. Liu, “The fundamental solution of the spacetime fractional advectiondispersion equation,” Journal of Applied Mathematics & Computing, vol. 18, no. 12, pp. 339–350, 2005. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 S. Momani, “An explicit and numerical solutions of the fractional KdV equation,” Mathematics and Computers in Simulation, vol. 70, no. 2, pp. 110–118, 2005. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 L. Debnath and D. D. Bhatta, “Solutions to few linear fractional inhomogeneous partial differential equations in fluid mechanics,” Fractional Calculus & Applied Analysis, vol. 7, no. 1, pp. 21–36, 2004. View at: Google Scholar  Zentralblatt MATH
 J. H. He, “Homotopy perturbation technique,” Computer Methods in Applied Mechanics and Engineering, vol. 178, no. 34, pp. 257–262, 1999. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 J. H. He, “A coupling method of a homotopy technique and a perturbation technique for nonlinear problems,” International Journal of NonLinear Mechanics, vol. 35, no. 1, pp. 37–43, 2000. View at: Publisher Site  Google Scholar
 J. H. He, “Periodic solutions and bifurcations of delaydifferential equations,” Physics Letters A, vol. 347, no. 4–6, pp. 228–230, 2005. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 J. H. He, “Application of homotopy perturbation method to nonlinear wave equations,” Chaos, Solitons and Fractals, vol. 26, no. 3, pp. 695–700, 2005. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 J. H. He, “Limit cycle and bifurcation of nonlinear problems,” Chaos, Solitons and Fractals, vol. 26, no. 3, pp. 827–833, 2005. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 J. H. He, Nonperturbative methods for strongly nonlinear problems [Dissertation], deVerlag im Internet GmbH, Berlin, Germany, 2006.
 J. H. He, “A new approach to nonlinear partial differential equations,” Communications in Nonlinear Science and Numerical Simulation, vol. 2, no. 4, pp. 230–235, 1997. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 J. H. He, “Variational iteration method for delay differential equations,” Communications in Nonlinear Science and Numerical Simulation, vol. 2, no. 4, pp. 235–236, 1997. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 M. Tatari and M. Dehghan, “On the convergence of He's variational iteration method,” Journal of Computational and Applied Mathematics, vol. 207, no. 1, pp. 121–128, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 Z. Odibat and S. Momani, “Numerical methods for nonlinear partial differential equations of fractional order,” Applied Mathematical Modelling, vol. 32, no. 1, pp. 28–39, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 R. Yulita Molliq, M. S. M. Noorani, and I. Hashim, “Variational iteration method for fractional heat and wavelike equations,” Nonlinear Analysis: Real World Applications, vol. 10, no. 3, pp. 1854–1869, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 R. Yulita Molliq, M. S. M. Noorani, I. Hashim, and R. R. Ahmad, “Approximate solutions of fractional ZakharovKuznetsov equations by VIM,” Journal of Computational and Applied Mathematics, vol. 233, no. 2, pp. 103–108, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 C. Chun, “Variational iteration method for a reliable treatment of heat equations with illdefined initial data,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 9, no. 4, pp. 435–440, 2008. View at: Google Scholar
 M. G. Porshokouhi and B. Ghanbari, “Application of He's variational iteration method for solution of the family of KuramotoSivashinsky equations,” Journal of King Saud University  Science, vol. 23, no. 4, pp. 407–411, 2011. View at: Publisher Site  Google Scholar
 S. T. MohyudDin and A. Yildirim, “An algorithm for solving the fractional vibration equation,” Computational Mathematics and Modeling, vol. 23, no. 2, pp. 228–237, 2012. View at: Publisher Site  Google Scholar
 A. M. A. ElSayed, A. Elsaid, I. L. ElKalla, and D. Hammad, “A homotopy perturbation technique for solving partial differential equations of fractional order in finite domains,” Applied Mathematics and Computation, vol. 218, no. 17, pp. 8329–8340, 2012. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 P. Rosenau and J. M. Hyman, “Compactons: solitons with finite wavelength,” Physical Review Letters, vol. 70, no. 5, pp. 564–567, 1993. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 R. Gorenflo and F. Mainardi, “Fractional calculus: integral and differential equations of fractional order,” in Fractals and Fractional Calculus, A. Carpinteri and F. Mainardi, Eds., Springer, New York, NY, USA, 1997. View at: Google Scholar
 M. Caputo, “Linear models of dissipation whose Q is almost frequency independent II,” Geophysical Journal of the Royal Astronomical Society, vol. 13, no. 5, pp. 529–539, 1967. View at: Publisher Site  Google Scholar
 M. Inokuti, H. Sekine, and T. Mura, “General use of the Lagrange multiplier in nonlinear mathematical physics,” in Variational Method in the Mechanics of Solids, S. Nemat Nasser, Ed., pp. 156–162, Pergamon Press, New York, NY, USA, 1978. View at: Google Scholar
 B. A. Finlayson, The Method of Weighted Residuals and Variational Principles, Academic Press, New York, NY, USA, 1972.
 P. A. Clarkson, E. L. Mansfield, and T. J. Priestley, “Symmetries of a class of nonlinear thirdorder partial differential equations,” Mathematical and Computer Modelling, vol. 25, no. 89, pp. 195–212, 1997. View at: Publisher Site  Google Scholar  Zentralblatt MATH
Copyright
Copyright © 2012 R. Yulita Molliq and M. S. M. Noorani. 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.