Recent Developments in Integral Transforms, Special Functions, and Their Extensions to Distributions Theory
View this Special IssueResearch Article  Open Access
Veyis Turut, Nuran Güzel, "Multivariate Padé Approximation for Solving Nonlinear Partial Differential Equations of Fractional Order", Abstract and Applied Analysis, vol. 2013, Article ID 746401, 12 pages, 2013. https://doi.org/10.1155/2013/746401
Multivariate Padé Approximation for Solving Nonlinear Partial Differential Equations of Fractional Order
Abstract
Two tecHniques were implemented, the Adomian decomposition method (ADM) and multivariate Padé approximation (MPA), for solving nonlinear partial differential equations of fractional order. The fractional derivatives are described in Caputo sense. First, the fractional differential equation has been solved and converted to power series by Adomian decomposition method (ADM), then power series solution of fractional differential equation was put into multivariate Padé series. Finally, numerical results were compared and presented in tables and figures.
1. Introduction
Recently, differential equations of fractional order have gained much interest in engineering, physics, chemistry, and other sciences. It can be said that the fractional derivative has drawn much attention due to its wide application in engineering physics [1–9]. Some approximations and numerical techniques have been used to provide an analytical approximation to linear and nonlinear differential equations and fractional differential equations. Among them, the variational iteration method, homotopy perturbation method [10, 11], and the Adomian decomposition method are relatively new approaches [5–9, 12, 13].
The decomposition method has been used to obtain approximate solutions of a large class of linear or nonlinear differential equations [12, 13]. Recently, the application of the method is extended for fractional differential equations [6–9, 14].
Many definitions and theorems have been developed for multivariate Padé approximations MPA (see [15] for a survey on multivariate Padé approximation). The multivariate Padé Approximation has been used to obtain approximate solutions of linear or nonlinear differential equations [16–19]. Recently, the application of the unvariate Padé approximation is extended for fractional differential equations [20, 21].
The objective of the present paper is to provide approximate solutions for initial value problems of nonlinear partial differential equations of fractional order by using multivariate Padé approximation.
2. Definitions
For the concept of fractional derivative, we will adopt Caputo’s definition, which is a modification of the RiemannLiouville definition and has the advantage of dealing properly with initial value problems in which the initial conditions are given in terms of the field variables and their integer order, which is the case in most physical processes. The definitions can be seen in [3, 4, 22, 23].
3. Decomposition Method [24]
Consider The decomposition method requires that a nonlinear fractional differential equation (1) is expressed in terms of operator form as where is a linear operator which might include other fractional derivatives of order less than , is a nonlinear operator which also might include other fractional derivatives of order less than is the Caputo fractional derivative of order , and is the source function [24].
Applying the operator [3, 4, 22, 23], the inverse of the operator , to both sides of (5) Odibat and Momani [24] obtained From this, the iterates are determined in [24] by the following recursive way:
4. Multivariate Padé Aproximation [25]
Consider the bivariate function with Taylor series development around the origin. We know that a solution of unvariate Padé approximation problem for is given by Let us now multiply th row in and by and afterwards divide th column in and by . This results in a multiplication of numerator and denominator by . Having done so, we get if .
This quotent of determinants can also immediately be written down for a bivariate function . The sum shall be replaced by th partial sum of the Taylor series development of and the expression by an expression that contains all the terms of degree in . Here a bivariate term is said to be of degree . If we define Then it is easy to see that and are of the form We know that and are called Padé equations [25]. So the multivariate Padé approximant of order for is defined as
5. Numerical Experiments
In this section, two methods, ADM and MPA, shall be illustrated by two examples. All the results are calculated by using the software Maple12. The full ADM solutions of examples can be seen from Odibat and Momani [24].
Example 1. Consider the nonlinear timefractional advection partial differential equation [24]
subject to the initial condition
Odibat and Momani [24] solved the problem using the decomposition method, and they obtained the following recurrence relation [24]:
where are the Adomian polynomials for the nonlinear function . In view of (15), the first few components of the decomposition series are derived in [24] as follows:
and so on; in this manner, the rest of components of the decomposition series can be obtained [24].
The first three terms of the decomposition series are given by [24]
For (16) is
Now, let us calculate the approximate solution of (18) for and by using Multivariate Padé approximation. To obtain multivariate Padé equations of (18) for and , we use (10). By using (10), we obtain
So, the multivariate Padé approximation of order for (17), that is,
For (17) is
For simplicity, let ; then
Using (10) to calculate the multivariate Padé equations for (22) we get
recalling that , we get multivariate Padé approximation of order for (21), that is, For (17) is For simplicity, let ; then Using (10) to calculate the multivariate Padé equations and then recalling that , we get multivariate Padé approximation of order for (25), that is, Table 1, Figures 1(a), 1(b), 1(c), 2(a), 2(b), 2(c), and 2(d) shows the approximate solutions for (13) obtained for different values of using the decomposition method (ADM) and the multivariate Padé approximation (MPA). The value of is for the exact solution [24].

(a)
(b)
(c)
(a)
(b)
(c)
(d)
Example 2. Consider the nonlinear timefractional hyperbolic equation [24]
subject to the initial condition
Odibat and Momani [24] solved the problem using the decomposition method, and they obtained the following recurrence relation in [24]:
where are the Adomian polynomials for the nonlinear function . In view of (30), the first few components of the decomposition series are derived in [24] as follows:
and so on; in this manner the rest of components of the decomposition series can be obtained.
The first three terms of the decomposition series (7) are given [24] by
For (43) is
Now, let us calculate the approximate solution of (33) for and by using multivariate Padé approximation. To obtain multivariate Padé equations of (33) for and , we use (10). By using (10), we obtain
So, the multivariate Padé approximation of order for (33), that is,
For (32) is
For simplicity, let ; then
Using (10) to calculate multivariate Padé equations of (37) for and , we use (10). By using (10), we obtain
recalling that , we get multivariate Padé approximation of order for (36), that is, For (32) is For simplicity, let ; then Using (10) to calculate multivariate Padé equations of (41) for and , we use (10). By using (10), we obtain
recalling that , we get multivariate Padé approximation of order for (40), that is, Table 2, Figures 3(a), 3(b), 3(c), 4(a), 4(b), 4(c), and 4(d) show the approximate solutions for (28) obtained for different values of using the decomposition method (ADM) and the multivariate Padé approximation (MPA). The value of is for the exact solution [24].

(a)
(b)
(c)
(a)
(b)
(c)
(d)
6. Concluding Remarks
The fundamental goal of this paper has been to construct an approximate solution of nonlinear partial differential equations of fractional order by using multivariate Padé approximation. The goal has been achieved by using the multivariate Padé approximation and comparing with the Adomian decomposition method. The present work shows the validity and great potential of the multivariate Padé approximation for solving nonlinear partial differential equations of fractional order from the numerical results. Numerical results obtained using the multivariate Padé approximation and the Adomian decomposition method are in agreement with exact solutions.
References
 J. H. He, “Nonlinear oscillation with fractional derivative and its applications,” in Proceedings of International Conference on Vibrating Engineering, pp. 288–291, Dalian, China, 1998. View at: Google Scholar
 J. H. He, “Some applications of nonlinear fractional differential equations and their approximations,” Bulletin of Science and Technology, vol. 15, no. 2, pp. 86–90, 1999. View at: Google Scholar
 Y. Luchko and R. Gorenflo, The Initial Value Problem for Some Fractional Differential Equations with the Caputo Derivative, Series A0898, Fachbreich Mathematik und Informatik, Freic Universitat Berlin, Berlin, Germany, 1998.
 I. Podlubny, Fractional Differential Equations, vol. 198 of Mathematics in Science and Engineering, Academic Press, San Diego, Calif, USA, 1999. View at: Zentralblatt MATH  MathSciNet
 S. Momani, “Nonperturbative analytical solutions of the space and timefractional Burgers equations,” Chaos, Solitons & Fractals, vol. 28, no. 4, pp. 930–937, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 Z. M. Odibat and S. Momani, “Application of variational iteration method to nonlinear differential equations of fractional order,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 7, no. 1, pp. 27–34, 2006. View at: Google Scholar
 S. Momani and Z. Odibat, “Analytical solution of a timefractional NavierStokes equation by Adomian decomposition method,” Applied Mathematics and Computation, vol. 177, no. 2, pp. 488–494, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 S. Momani and Z. Odibat, “Numerical comparison of methods for solving linear differential equations of fractional order,” Chaos, Solitons & Fractals, vol. 31, no. 5, pp. 1248–1255, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 Z. M. Odibat and S. Momani, “Approximate solutions for boundary value problems of timefractional wave equation,” Applied Mathematics and Computation, vol. 181, no. 1, pp. 767–774, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 G. Domairry and N. Nadim, “Assessment of homotopy analysis method and homotopy perturbation method in nonlinear heat transfer equation,” International Communications in Heat and Mass Transfer, vol. 35, no. 1, pp. 93–102, 2008. View at: Publisher Site  Google Scholar
 G. Domairry, M. Ahangari, and M. Jamshidi, “Exact and analytical solution for nonlinear dispersive $K(m,p)$ equations using homotopy perturbation method,” Physics Letters A, vol. 368, no. 34, pp. 266–270, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 G. Adomian, “A review of the decomposition method in applied mathematics,” Journal of Mathematical Analysis and Applications, vol. 135, no. 2, pp. 501–544, 1988. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 G. Adomian, Solving Frontier Problems of Physics: The Decomposition Method, vol. 60 of Fundamental Theories of Physics, Kluwer Academic Publishers, Dordrecht, The Netherlands, 1994. View at: MathSciNet
 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  MathSciNet
 Ph. Guillaume and A. Huard, “Multivariate Padé approximation,” Journal of Computational and Applied Mathematics, vol. 121, no. 12, pp. 197–219, 2000. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 V. Turut, E. Çelik, and M. Yiğider, “Multivariate Padé approximation for solving partial differential equations (PDE),” International Journal for Numerical Methods in Fluids, vol. 66, no. 9, pp. 1159–1173, 2011. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 V. Turut and N. Güzel, “Comparing numerical methods for solving timefractional reactiondiffusion equations,” ISRN Mathematical Analysis, vol. 2012, Article ID 737206, 28 pages, 2012. View at: Google Scholar  Zentralblatt MATH  MathSciNet
 V. Turut, “Application of Multivariate padé approximation for partial differential equations,” Batman University Journal of Life Sciences, vol. 2, no. 1, pp. 17–28, 2012. View at: Google Scholar
 V. Turut and N. Güzel, “On solving partial differential eqauations of fractional order by using the variational iteration method and multivariate padé approximation,” European Journal of Pure and Applied Mathematics. Accepted. View at: Google Scholar
 S. Momani and R. Qaralleh, “Numerical approximations and Padé approximants for a fractional population growth model,” Applied Mathematical Modelling, vol. 31, no. 9, pp. 1907–1914, 2007. View at: Publisher Site  Google Scholar
 S. Momani and N. Shawagfeh, “Decomposition method for solving fractional Riccati differential equations,” Applied Mathematics and Computation, vol. 182, no. 2, pp. 1083–1092, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 K. B. Oldham and J. Spanier, The Fractional Calculus, Academic Press, New York, NY, USA, 1974. View at: MathSciNet
 M. Caputo, “Linear models of dissipation whose Q is almost frequency independent. Part II,” Journal of the Royal Astronomical Society, vol. 13, no. 5, pp. 529–539, 1967. View at: Publisher Site  Google Scholar
 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
 A. Cuyt and L. Wuytack, Nonlinear Methods in Numerical Analysis, vol. 136 of NorthHolland Mathematics Studies, NorthHolland Publishing, Amsterdam, The Netherlands, 1987. View at: MathSciNet
Copyright
Copyright © 2013 Veyis Turut and Nuran Güzel. 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.