Research Article | Open Access

Volume 2014 |Article ID 398286 | 8 pages | https://doi.org/10.1155/2014/398286

# Numerical Methods for Fractional Order Singular Partial Differential Equations with Variable Coefficients

Accepted22 Jan 2014
Published10 Mar 2014

#### 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, [16].

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 [1012] and variational iteration method [1316] 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 Riemann-Liouville 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 N-term 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 Riemann-Liouville 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 first-order 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 fourth-order 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 third-order term series was used in evaluating the approximate solutions for Table 1.

 0.1 0.000144 0.000141 0.000139 0.000139 0.000136 0.000136 0.000136 0.2 0.001149 0.001127 0.034155 0.001108 0.001089 0.001089 0.001089 0.3 0.003867 0.003793 0.003740 0.003729 0.003668 0.003668 0.003667 0.4 0.009134 0.008960 0.007492 0.008809 0.008663 0.008664 0.008664 0.5 0.017759 0.017421 0.017175 0.017128 0.016845 0.016844 0.016845 0.6 0.030519 0.029939 0.029516 0.029435 0.028949 0.028948 0.028948 0.7 0.048149 0.047233 0.046566 0.046438 0.045671 0.045671 0.045671 0.8 0.071336 0.069978 0.068989 0.068801 0.067663 0.067663 0.067663 0.9 0.100708 0.098792 0.097396 0.097130 0.095524 0.095524 0.095524

Example 6. Consider the following singular two-dimensional 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 third-order term series was used in evaluating the approximate solutions for Table 2.

 0.1 0.389368 0.389267 0.394618 0.391901 0.397339 0.394672 0.397339 0.2 0.389368 0.389927 0.394618 0.391901 0.397339 0.394672 0.397339 0.3 0.389368 0.389268 0.394618 0.391902 0.397339 0.394672 0.397339 0.4 0.389370 0.389270 0.394620 0.391904 0.397341 0.394674 0.397339 0.5 0.389376 0.389276 0.394627 0.391910 0.397347 0.394681 0.397341 0.6 0.389393 0.389293 0.394644 0.391927 0.397364 0.394698 0.397347 0.7 0.389432 0.389332 0.394682 0.391965 0.397404 0.394736 0.397344 0.8 0.389511 0.389409 0.394762 0.392044 0.397483 0.394816 0.397483 0.9 0.389655 0.389555 0.394991 0.392192 0.397619 0.394963 0.397632

#### 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

1. 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 II-Computational 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
2. 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
3. R. Hilfer, Applications of Fractional Calculus in Physics, World scientific, Singapore, 2000.
4. 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
5. I. Podlubny, “Fractional differential equation,” in Mathematics in Science and Engineering, vol. 198, Academic Press, San Diego, Calif, USA, 1999. View at: Google Scholar
6. K. S. Miller and B. Ross, An Introduction to the Fractional Calculus and Fractional Differential Equations, John Wiley & Sons, New York, NY, USA, 1993.
7. S. T. Mohydud-Din, A. Yildirim, M. M. Hosseini, and Y. Khan, “A study on systems of variable-coefficient singular parabolic partial differential equations,” World Applied Sciences Journal, vol. 10, no. 11, pp. 1321–1327, 2010. View at: Google Scholar
8. M. Abokhald, “Varational iteration method for nonlinear singular two-point boundary value problems arising in human physiology,” Journal of Mathematics, vol. 2013, Article ID 720134, 4 pages, 2013. View at: Publisher Site | Google Scholar
9. 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
10. J. H. He, “Homotopy perturbation technique,” Computer Methods in Applied Mechanics and Engineering, vol. 178, no. 3-4, pp. 257–262, 1999.
11. J. H. He, “Coupling method of a homotopy technique and a perturbation technique for non-linear problems,” International Journal of Non-Linear Mechanics, vol. 35, no. 1, pp. 37–43, 2000. View at: Publisher Site | Google Scholar
12. J. H. He, Non-pertubation methods for strongly nonlinear problems [Dissertation], de-Verlag im Internet Gmbh, Berlin, Germany, 2006.
13. 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.
14. 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
15. 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. 1-2, pp. 57–68, 1998.
16. J. H. He, “Variational iteration method—a kind of non-linear analytical technique: some examples,” International Journal of Non-Linear Mechanics, vol. 34, no. 4, pp. 699–708, 1999.
17. 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. View at: Google Scholar

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.