Research Article  Open Access
Numerical Methods for Fractional Order Singular Partial Differential Equations with Variable Coefficients
Abstract
We implement relatively analytical methods, the homotopy perturbation method and the variational iteration method, for solving singular fractional partial differential equations of fractional order. The process of the methods which produce solutions in terms of convergent series is explained. The fractional derivatives are described in Caputo sense. Some examples are given to show the accurate and easily implemented of these methods even with the presence of singularities.
1. Introduction
In the last decades, fractional calculus found many applications in various fields of engineering and physical sciences such as physics, chemistry, biology, economy, viscoelasticity, electrochemistry, electromagnetic, relaxation processes, diffusion, control, porous media, and many more; see, for example, [1–6].
The numerical solution of singular differential equations of integer order has been a hot topic in numerical and computational mathematics for a long time [7, 8]. Singular partial differential equations of fractional order, as generalizations of classical singular partial differential equations of integer order, are increasingly used to model problems in physics and engineering. Consequently, considerable attention has been given to the solution of singular partial differential equations of fractional order.
In most of these equations analytical solutions are either quite difficult or impossible to achieve, so approximations and numerical techniques must be used. Several methods have been used to solve singular differential equations such as variational iteration method [8], homotopy perturbation method [7], and homotopy analysis method [9].
The homotopy perturbation method [10–12] and variational iteration method [13–16] are relatively new approach to provide an analytical approximation to linear and nonlinear problems, and they are particularly valuable as tools for scientists and applied mathematicians, because it provides immediate and visible symbolic terms of analytic solutions, as well as numerical approximate solutions to both linear and nonlinear differential equations.
Recently, the application of these methods extended for fractional differential equations. He [15] was the first to apply the variational iteration method to fractional differential equations.
The objective of the present paper is to extend the application of homotopy perturbation method to provide approximate solutions and to make comparison with that obtained by the variational iteration method for singular partial differential equations of fractional order: where , , subject to initial conditions and boundary conditions where is the fractional derivative in the Caputo sense, , and and , are continuous. The paper is organized as follows. In Section 2, some basic definitions and properties of fractional calculus theory are given. In Section 3, the basic idea of HPM is given. In Section 4, the basic idea of variational iteration method is given. In Sections 5 and 6 the analysis of HPM and VIM for singular fractional partial differential equations, respectively, is given. Some examples are given in Section 7. Concluding remarks are listed in Section 8.
2. Preliminaries
In this section, we give some basic definitions and properties of fractional calculus theory which is used in this paper.
Definition 1. A real function , is said to be in space , if there exists a real number , such that where , and it is said to be in the space if , .
Definition 2. The RiemannLiouville fractional integral operator of order of a function , is defined as in particular .
For and , some properties of the operator (i),(ii),(iii).
Definition 3. The Caputo fractional derivative of , is defined as
Lemma 4. if , , , then the following two properties hold: (i), (ii).
3. Homotopy Perturbation Method
To illustrate the basic idea of this method, we 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 .
In general, the operator can be divided into two parts and , where is linear, while is nonlinear. Equation (6) therefore can be rewritten as follows: By the homotopy technique [10, 11], we construct a homotopy which satisfies or where is an embedding parameter and is an initial approximation of (6) which satisfies the boundary conditions.
From (9) and (10) we have The change in the process of from zero to unity is just that of from to . In topology this is called deformation and and are called homotopic.
Now, assume that the solution of (9) and (10) can be expressed as Setting results in the approximate solution of (6). Consider
4. The Variational Iteration Method
To illustrate the basic concepts of VIM, we consider the following differential equation: where is a linear operator, is a nonlinear operator, and is a nonhomogeneous term.
According to VIM, one constructs a correction functional as follows: where is a general Lagrange multiplier, and denotes restricted variation; that is, .
5. Analysis of Homotopy Perturbation Method
To illustrate the basic concepts of HPM for (1) with (2) and (3), we use the view of He [10, 11], where the following homotopy was constructed for (1): or where is an embedding parameter. If , (17) becomes linear fractional differential equation and when , (17) turns out to be the original equation. In view of basic assumption of HPM, solution of (1) can be expressed as a power series in When , we get the approximate solution of (1) The convergence of series (20) has been proved in [12].
Substituting (19) into (17) and then equating the terms with identical power of , we obtain the following series of linear equations: with the initial conditions (2) and boundary conditions (3).
The initial approximation can be chosen in the following manner. The equations in the system (21) can be solved by applying the operator and by simple computation; we approximate the series solution of HPM by the following Nterm truncated series: which is the approximate solution of (1) with (2) and (3).
6. Analysis of Variational Iteration Method
To solve the fractional singular partial differential equations by using the variational iteration method, with initial and boundary conditions (2) and (3), we construct the following correction functional: or where and is the RiemannLiouville fractional integral operator of order , with respect to variable , is a general Lagrange multiplier which can be identified optimally by variational theory [17], and are considered as restricted variation; that is, .
Making the above correction functional stationary, the following condition can be obtained: which yields to Lagrange multiplier We obtain the following iteration formula by substitution of (28) in (25): That is, This yields the following iteration formula: The initial approximation can be chosen by the following manner which satisfies initial conditions: We can obtain the following firstorder approximation by substitution of (32) in (31): Finally, by substituting the constant values of and in (33), we have the results as the approximate solutions of (1) with (2) and (3).
7. Applications
In this section we have applied homotopy perturbation method and variational iteration method to fractional singular partial differential equations with known exact solution.
Example 5. Consider the following fourthorder fractional singular partial differential equation:
With initial conditions
and with boundary conditions
The exact solution in special case is
First, we construct
Substituting (19) in (38) and then equating the terms with same powers of , we get the series
Now applying the operator to (39) and using initial and boundary conditions yield
so that the solution of the problem (34) is given by
Now, we solve by variational iteration method. According to variational iteration method, the formula (31) for (34) can be expressed in the following form:
Suppose that an initial approximation has the following form which satisfies the initial conditions:
Now by iteration formula (42), we obtain the first approximation:
and second approximation:
Table 1 shows the approximate solution for (34) obtained for different values of using the homotopy perturbation method and the variational iteration method. The value is the only case for which we know the exact solution and our approximate solution using the variational iteration method is more accurate than approximate solution obtained using homotopy perturbation method. It is to be noted that only the thirdorder term series was used in evaluating the approximate solutions for Table 1.

Example 6. Consider the following singular twodimensional partial differential equation of fractional order:
With initial conditions
and with boundary conditions
and the exact solution in special case is
According to HPM, we construct the following homotopy:
Substituting (19) in (50) and then equating the terms with same powers of , we get the series
Now applying the operator to (51) and using initial and boundary conditions yield
so that the solution of the problem (46) is given by
Now, we solve the problem (46) by variational iteration method. According to variational iteration method, the formula (31) for (46) can be expressed in the following form:
Suppose that an initial approximation has the following form which satisfies the initial conditions
Now by iteration formula (42), we obtain the following approximations:
The second approximation takes the following form:
Table 2 shows the approximate solution for (46) obtained for different values of using the homotopy perturbation method and the variational iteration method. The value is the only case for which we know the exact solution , and our approximate solution using homotopy perturbation method is more accurate than approximate solution obtained using the variational iteration method. It is to be noted that only the thirdorder term series was used in evaluating the approximate solutions for Table 2.

8. Conclusion
The essential goal of this work has been to construct an approximate solution of singular partial differential equations of fractional order. The goal has been done by using the homotopy perturbation method and the variational iteration method. The methods provide the solutions in terms of convergent series with easily computable component even with the presence of singularities. Two examples were presented in Section 7 to show the accuracy of these methods; in Example 5 it seems that the approximate solution using the variational iteration method converges faster than approximate solution using homotopy perturbation method while the approximate solution in Example 6 using homotopy perturbation method converges faster than the approximate solution obtained using the variational iteration method to the exact solution. The fact that the proposed methods solve nonlinear problems without using Adomian’s polynomial can be considered as advantage of these techniques over the decomposition method.
Conflict of Interests
The authors declare that they have no conflict of interests regarding the publication of this paper.
Acknowledgment
The authors express their sincere thanks to the referees for the careful and noteworthy reading of the paper and for their very helpful suggestions that improved the paper substantially. The authors also gratefully acknowledge that this work is partially supported by the University Putra Malaysia Grant no. 5527179.
References
 K. Diethelm and A. D. Freed, “On the solution of nonlinear fractional order differential equations used in the modeling of viscoelasticity,” in Scientific Computing in Chemical Engineering IIComputational Fluid Dynamics, Reaction Engineering and Molecular Properties, F. Keil, W. Mackens, H. Voss, and J. Werther, Eds., pp. 217–224, Springer, Heidelberg, Germany, 1999. View at: Google Scholar
 R. Metzler, W. Schick, H. G. Kilian, and T. F. Nonnenmacher, “Relaxation in filled polymers: a fractional calculus approach,” The Journal of Chemical Physics, vol. 103, no. 16, pp. 7180–7186, 1995. View at: Publisher Site  Google Scholar
 R. Hilfer, Applications of Fractional Calculus in Physics, World scientific, Singapore, 2000.
 K. B. Oldham and J. Spanier, “The fractional calculus,” in Mathematics in Science and Engineering, vol. 198, Academic Press, San Diego, Calif, USA, 1974. View at: Google Scholar
 I. Podlubny, “Fractional differential equation,” in Mathematics in Science and Engineering, vol. 198, Academic Press, San Diego, Calif, USA, 1999. View at: Google Scholar
 K. S. Miller and B. Ross, An Introduction to the Fractional Calculus and Fractional Differential Equations, John Wiley & Sons, New York, NY, USA, 1993.
 S. T. MohydudDin, A. Yildirim, M. M. Hosseini, and Y. Khan, “A study on systems of variablecoefficient singular parabolic partial differential equations,” World Applied Sciences Journal, vol. 10, no. 11, pp. 1321–1327, 2010. View at: Google Scholar
 M. Abokhald, “Varational iteration method for nonlinear singular twopoint boundary value problems arising in human physiology,” Journal of Mathematics, vol. 2013, Article ID 720134, 4 pages, 2013. View at: Publisher Site  Google Scholar
 S. Abbasbsndy and E. Shivanian, “Solution of singular linear vibrational BVPs by the homotopy analysis method,” Journal of Numerical Mathematics and Stochastics, vol. 1, no. 1, pp. 77–84, 2009. View at: Google Scholar
 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, “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, Nonpertubation methods for strongly nonlinear problems [Dissertation], deVerlag im Internet Gmbh, Berlin, Germany, 2006.
 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
 J. H. He, “Approximate solution of non linear differential equations with convolution product nonlinearities,” Computer Methods in Applied Mechanics and Engineering, vol. 167, pp. 69–73, 1998. 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
 J. H. He, “Variational iteration method—a kind of nonlinear analytical technique: some examples,” International Journal of NonLinear Mechanics, vol. 34, no. 4, pp. 699–708, 1999. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 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. NematNasser, Ed., pp. 156–162, Pergamon Press, Oxford, UK, 1978. View at: Google Scholar
Copyright
Copyright © 2014 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.