• Views 1,641
• Citations 12
• ePub 45
• PDF 1,181
`Mathematical Problems in EngineeringVolume 2013, Article ID 543848, 10 pageshttp://dx.doi.org/10.1155/2013/543848`
Research Article

## Fractional Variational Iteration Method and Its Application to Fractional Partial Differential Equation

1Faculty of Science and Technology, Universiti Sains Islam Malaysia, 71800 Nilai, Malaysia
2Department of Mathematics, Faculty of Science, University Putra Malaysia, 43400 UPM, Serdang, Selangor Darul Ehsan, Malaysia

Received 13 March 2013; Accepted 5 June 2013

Academic Editor: Mufid Abudiab

Copyright © 2013 Asma Ali Elbeleze 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

We use the fractional variational iteration method (FVIM) with modified Riemann-Liouville derivative to solve some equations in fluid mechanics and in financial models. The fractional derivatives are described in Riemann-Liouville sense. To show the efficiency of the considered method, some examples that include the fractional Klein-Gordon equation, fractional Burgers equation, and fractional Black-Scholes equation are investigated.

#### 1. Introduction

The topic of fractional calculus (theory of integration and differentiation of an arbitrary order) was started over 300 years ago. Recently, fractional differential equations have attracted many scientists and researchers due to the tremendous use in fluid mechanics, mathematical biology, electrochemistry, and physics. For example, differential equations with fractional order have recently proved to be suitable tools to modeling of many physical phenomena [1] and the fluid-dynamic traffic model with fractional derivative [2], and nonlinear oscillation of earthquake can be modeled with fractional derivatives [3].

There are several types of time fractional differential equations. (1) Fractional Klein-Gordon equations This model is obtained by replacing the order time derivative with the fractional derivative of order . The linear and nonlinear Klein-Gordon equations are used to modeling many problems in classical and quantum mechanics and condensed matter physics. For example, nonlinear sine Klein-Gordon equation models a Josephson junction [4, 5]. (2) Fractional Burger's equation In general, fractional Burger's model is derived from well-known Burger's equation model by replacing the ordinary time derivatives to fractional order time derivatives. Reference [6] has investigated unsteady flows of viscoelastic fluids with fractional Burger's model and fractional generalized Burger's model through channel (annulus) tube and solutions for velocity field.(3)Fractional Black-Scholes European option pricing equationsIn financial model the fractional Black-Scholes equation is obtained by replacing the order of derivative with a fractional derivative order [10]. where is the European call option price at asset price and at time , is the maturity, is the risk-free interest rate, and represents the volatility function of underlying asset.The payoff functions are where and are the value of the European call and put options, respectively, denotes the expiration price for the option, and the function gives the large value between and . The Black-Scholes equation is one of the most significant mathematical models for a financial market. This equation is used to submit a reasonable price for call or put options based on factors such as underlying stock volatility and days to expiration.

Formerly, [7] investigated approximate analytical solution of fractional nonlinear Klein-Gordon equation (1) when by using HPM, while [8] solved this equation by using HAM also when . Reference [9] solved the coupled Klein-Gordon equation with time fractional derivative by ADM. References [10, 11] solved fractional Black-Scholes equations by using HPM using Sumudu and Laplace transforms, respectively. Reference [12] gave the exact solution of fractional Burgers equation, while [13] used DTM to find the approximate and exact solution of space- and time fractional Burgers equations. Reference [14] solved this equation by using VIM.

The variational iteration method [1529] is one of approaches to provide an analytical approximation solutions to linear and nonlinear problems. The fractional variational iteration method with Riemann-Liouville derivative was proposed by Wu and Lee [30] and applied to solve time fractional and space fractional diffusion equations. Furthermore Wu [31] explained a possible use of the fractional variational iteration method as a fractal multiscale method. Recently fractional variational iteration method has been used to obtain approximate solutions of fractional Riccati differential equation [32].

The objective of this paper is to extend the application of the fractional variational iteration method to obtain analytical approximate solution for some fractional partial differential equations. These equations include fractional Klein-Gordon equation (1), Burgers equation (2), and fractional Black-Scholes equations (3).

Motivated and inspired by the ongoing research in this field, we will consider the following time fractional differential equation: with initial condition where is modified Riemann-Liouville derivative [3335] of order defined in Section 2, and are continuous functions, are linear and nonlinear operators, and is unknown function.

To solve the problem (1)-(2), we consider the FVIM in this work. This method is based on variational iteration method [19, 36] and modified Riemann-Liouville derivatives proposed by Jumarie.

This paper is organized as follows. in Section 2 some basic definitions of fractional calculus theory are given. In Section 3, the solution procedure of the fractional iteration method is given; we present the application of the FVIM for some fractional partial differential equations in Section 4. The conclusions are drawn in Section 5.

#### 2. Fractional Calculus

##### 2.1. Fractional Derivative via Fractional Difference

Definition 1. The left-sides Riemann-Liouville fractional integral operator of order , of a function , , is a defined as

Definition 2. The modified Riemann-Liouville derivative [34, 35] is defined as where , , and .

Definition 3. Let , denote a continuous (but not necessarily differentiable) function, and let denote a constant discretization span. Define the forward operator by the equality Then the fractional difference of order , , of is defined by the expression and its fractional derivative of order is defined by the limit

Equation (11) is defined as Jumarie fractional derivative of order which is equivalent to (8). For more details we refer the reader to [35].

For , some properties of the fractional modified Riemann-Liouville derivative.

Fractional Leibnitz product law: fractional Leibnitz formulation: The fractional integration by parts formula:

Definition 4. Fractional derivative of compounded function [34, 35] is defined as

Definition 5 (see [34, 35]). The integral with respect to is defined as the solution of the fractional differential equation

Lemma 6 (see [34, 35]). Let denote a continuous function; then the solution of (2) is defined as that is,

For example, with in (7), one obtains

Definition 7. The Mittag-Leffler function with is defined by the following series representation, valid in the whole complex plane [37]:

#### 3. Fractional Variational Iteration Method

To describe the solution procedure of fractional variational iteration method, we consider the time-fractional differential equations (1)–(3).

According to variational iteration method we construct the following correction function: where is the general Lagrange multiplier which can be defined optimally via variational theory [22] and is the restricted variation, that is, .

By using (7), we obtain a new correction functional Making the above functional stationary the following conditions can be obtained: Now, we can get the coefficients of to zero: So, the generalized Lagrange multiplier can be identified as Then we obtain the following iteration formula by substituting (25) in (23): where and is an initial approximation which can be freely chosen if it satisfies the initial and boundary conditions of the problem.

#### 4. Applications

In this section, we have applied fractional variational iteration method (FVIM) to fractional partial differential equations.

Example 8. In this example we consider the following fractional nonlinear Klein-Gordon differential equation: subject to initial condition Substituting , and in (1). Construction the iteration formula as follows: Taking the initial value we can derive the first approximate as follows: Thus, the approximate solution is
In Figures 1 and 2 we have shown the surface of corresponding to the values for FVIM and HPM; the two figures indicate that the differences among VIM and HPM, and the exact solution in Example 8 are negligible when while when the results of VIM and HPM somewhat diverge from the exact solution.

Figure 1: The surface shows the solution for (27) with initial condition (28): FVIM results are, respectively, (a) and (c) ; HPM [7] results are, respectively, (b) and (d) .
Figure 2: The surface shows the solution for (27) with initial condition (28): (a) FVIM when , (b) HPM [7] when , and (c) exact solution.

Example 9. We consider the one-dimensional linear inhomogeneous fractional Burger equation subject to initial condition By construction the iteration formula as follows: Taking the initial value we can derive the first approximate as follows: So, the exact solution follows immediately. The exact solution is obtained by using two iterations and this is dependent on proper selection of initial guess .

Example 10. We consider the following fractional Black-Scholes option pricing equation [38] as follows: where is the risk-free interest rate subject to initial condition The exact solution for special case is given by By construction the iteration formula as follows: Taking the initial value we can derive the first approximate as follows: so that the solution of the problem is given by where (z) is Mittag-Leffler function in one parameter. Equation (41) represents the closed form solution of the fractional Black-Scholes equation (36). Now for the standard case , this series has the closed form of the solution , which is an exact solution of the given Black-Scholes equation (36) for .

In Figure 3 we have shown the surface of corresponding to the value ( for FVIM&HPM and for FVIM ).

Figure 3: The surface shows the solution for (36) with initial condition (37): (a) FVIM (), (b) HPM [10] (), and (c) FVIM ().

#### 5. Conclusion

Variational iteration method has been known as a powerful method for solving many fractional equations such as partial differential equations, integrodifferential equations, and so many other equations. In this paper, based on the variational iteration method and modified Riemann-Liouville derivative, we have presented a general framework of fractional variational iteration method for analytical and numerical treatment of fractional partial differential equations in fluid mechanics and in financial models. All of the examples concluded that the fractional variational iteration method is powerful and efficient in finding analytical approximate solutions as well as numerical solutions. For example, the results of Examples 8 and 10 illustrate that the present method is in excellent agreement with those of HPM and exact solution, where the obtained solution is shown graphically. Further, in Example 9 we got the exact solution in two iterations. The basic idea described in this paper is expected to be further employed to solve other similar linear and nonlinear problems in fractional calculus. Maple has been used for presenting graph of solution in the present paper.

#### Acknowledgments

The second author gratefully acknowledges that this research partially supported by Ministry of Higher Education (MOHE), Malaysia under the ERGS Grant 5527068.

#### References

1. I. Podlubny, Fractional Differential Equations, vol. 198 of Mathematics in Science and Engineering, Academic Press, San Diego, Calif, USA, 1999.
2. 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.
3. J. H. He, “Nonlinear oscillation with fractional derivative and its applications,” in Proceedings of the International Conference on Vibrating Engineering, pp. 288–291, Dalian, Chaina, 1998.
4. A. Barone, F. Esposito, C. J. Magee, and A. C. Scott, “Theory and applications of the sine-gordon equation,” La Rivista del Nuovo Cimento, vol. 1, no. 2, pp. 227–267, 1971.
5. E. Yusufoglu, “The variational iteration method for studying the Klein-Gordon equation,” Applied Mathematics Letters, vol. 21, no. 7, pp. 669–674, 2008.
6. M. Khan, S. Hyder Ali, and H. Qi, “On accelerated flows of a viscoelastic fluid with the fractional Burgers' model,” Nonlinear Analysis. Real World Applications. An International Multidisciplinary Journal, vol. 10, no. 4, pp. 2286–2296, 2009.
7. A. K. Golmankhaneh, A. K. Golmankhaneh, and D. Baleanu, “On nonlinear fractional KleinGordon equation,” Signal Processing, vol. 91, no. 3, pp. 446–451, 2011.
8. M. Kurulay, “Solving the fractional nonlinear Klein-Gordon equation by means of the homotopy analysis method,” Advances in Difference Equations, p. 2012, 2012.
9. E. Hesameddini and F. Fotros, “Solution for time-fractional coupled Klein-Gordon Schrodinger equation using decomposition method,” International Mathematical Forum, vol. 7, no. 21–24, pp. 1047–1056, 2012.
10. A. A. Elbeleze, A. Kiliçman, and B. M. Taib, “Homotopy perturbation method for fractional black-scholes european option pricing equations using Sumudu transform,” Mathematical Problems in Engineering, vol. 2013, Article ID 524852, 7 pages, 2013.
11. S. Kumar, A. Yildirim, Y. Khan, H. Jafari, K. Sayevand, and L. Wei, “Analytical solution of fractional Black-Scholes European option pricing equation by using Laplace transform,” Journal of Fractional Calculus and Applications, vol. 2, no. 8, pp. 1–9, 2012.
12. C. Xue, J. Nie, and W. Tan, “An exact solution of start-up flow for the fractional generalized Burgers' fluid in a porous half-space,” Nonlinear Analysis. Theory, Methods & Applications A: Theory and Methods, vol. 69, no. 7, pp. 2086–2094, 2008.
13. M. Kurulay, “The approximate and exact solutions of the space and time-fractional Burggres equations,” International Journal of Research and Reviews in Applied Sciences, vol. 3, no. 3, pp. 257–263, 2010.
14. Z. Odibat and S. Momani, “The variational iteration method: an efficient scheme for handling fractional partial differential equations in fluid mechanics,” Computers & Mathematics with Applications, vol. 58, no. 11-12, pp. 2199–2208, 2009.
15. D. D. Ganji and A. Sadighi, “Application of homotopy-perturbation and variational iteration methods to nonlinear heat transfer and porous media equations,” Journal of Computational and Applied Mathematics, vol. 207, no. 1, pp. 24–34, 2007.
16. M. Rafei, D. D. Ganji, H. Daniali, and H. Pashaei, “The variational iteration method for nonlinear oscillators with discontinuities,” Journal of Sound and Vibration, vol. 305, no. 4-5, pp. 614–620, 2007.
17. D. D. Ganji, M. Jannatabadi, and E. Mohseni, “Application of He's variational iteration method to nonlinear Jaulent-Miodek equations and comparing it with ADM,” Journal of Computational and Applied Mathematics, vol. 207, no. 1, pp. 35–45, 2007.
18. S. Momani and S. Abuasad, “Application of He's variational iteration method to Helmholtz equation,” Chaos, Solitons & Fractals, vol. 27, no. 5, pp. 1119–1123, 2006.
19. 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.
20. J. He, “Semi-inverse 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.
21. J. H. He and X. H. Wu, “Variational iteration method: new development and applications,” Computers & Mathematics with Applications, vol. 54, no. 7-8, pp. 881–894, 2007.
22. M. Inokuti, H. Sekine, and T. Mura, “General use of the Lagrange multiplier in non-linear mathematical physics,” in Variational Method in the Mechanics of Solids, S. Nemat-Nasser, Ed., pp. 156–162, Pergamon Press, Oxford, UK, 1978.
23. A. M. Wazwaz, “The variational iteration method for solving linear and nonlinear systems of PDEs,” Computers & Mathematics with Applications, vol. 54, no. 7-8, pp. 895–902, 2007.
24. A. M. Wazwaz, “The variational iteration method: a reliable analytic tool for solving linear and nonlinear wave equations,” Computers & Mathematics with Applications, vol. 54, no. 7-8, pp. 926–932, 2007.
25. A. M. Wazwaz, “The variational iteration method: a powerful scheme for handling linear and nonlinear diffusion equations,” Computers & Mathematics with Applications, vol. 54, no. 7-8, pp. 933–939, 2007.
26. Z. Odibat, “Reliable approaches of variational iteration method for nonlinear operators,” Mathematical and Computer Modelling, vol. 48, no. 1-2, pp. 222–231, 2008.
27. E. Yusufoglu, “Variational iteration method for construction of some compact and noncompact structures of Klein-Gordon equations,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 8, no. 2, pp. 153–158, 2007.
28. J. Biazar and H. Ghazvini, “He's variational iteration method for solving hyperbolic differential equations,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 8, no. 3, pp. 311–314, 2007.
29. H. Ozer, “Application of the variational iteration method to the boundary value problems with jump discontinuities arising in solid mechanics,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 8, no. 4, pp. 513–518, 2007.
30. G. C. Wu and E. W. M. Lee, “Fractional variational iteration method and its application,” Physics Letters A, vol. 374, no. 25, pp. 2506–2509, 2010.
31. G. C. Wu, “New trends in the variational iteration method,” Communications in Fractional Calculus, vol. 2, pp. 59–75, 2011.
32. M. Merdan, “On the solutions fractional riccati differential equation with modified Riemann-Liouville derivative,” International Journal of Differential Equations, vol. 2012, Article ID 346089, 17 pages, 2012.
33. G. Jumarie, “Stochastic differential equations with fractional Brownian motion input,” International Journal of Systems Science, vol. 24, no. 6, pp. 1113–1131, 1993.
34. G. Jumarie, “Laplace's transform of fractional order via the Mittag-Leffler function and modified Riemann-Liouville derivative,” Applied Mathematics Letters, vol. 22, no. 11, pp. 1659–1664, 2009.
35. G. Jumarie, “Table of some basic fractional calculus formulae derived from a modified Riemann-Liouville derivative for non-differentiable functions,” Applied Mathematics Letters, vol. 22, no. 3, pp. 378–385, 2009.
36. A.-M. Wazwaz, “The variational iteration method for solving two forms of Blasius equation on a half-infinite domain,” Applied Mathematics and Computation, vol. 188, no. 1, pp. 485–491, 2007.
37. F. Mainardi, “On the initial value problem for the fractional diffusion-wave equation,” in Waves and Stability in Continuous Media, S. Rionero and T. Ruggeeri, Eds., pp. 246–251, World Scientific, Singapore, 1994.
38. V. Gülkac, “The homotopy perturbation method for the Black-Scholes equation,” Journal of Statistical Computation and Simulation, vol. 80, no. 12, pp. 1349–1354, 2010.