Abstract and Applied Analysis

Volume 2013 (2013), Article ID 617010, 9 pages

http://dx.doi.org/10.1155/2013/617010

Research Article

New Iterative Method: An Application for Solving Fractional Physical Differential Equations

A. A. Hemeda

Department of Mathematics, Faculty of Science, Tanta University, Tanta 31527, Egypt

Received 12 September 2012; Accepted 7 March 2013

Academic Editor: Soon Y. Chung

Copyright © 2013 A. A. Hemeda. 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.

The new iterative method with a powerful algorithm is developed for the solution of linear and nonlinear ordinary and partial differential equations of fractional order as well. The analysis is accompanied by numerical examples where this method, in solving them, is used without linearization or small perturbation which con…firm the power, accuracy, and simplicity of the given method compared with some of the other methods.

1. Introduction

Considerable attention has been devoted to the study of the fractional calculus during the past three decades and its numerous applications in the area of physics and engineering. The applications of fractional calculus used in many fields such as electrical networks, control theory of dynamical systems, probability and statistics, electrochemistry, chemical physics, optics, and signal processing can be successfully modelled by linear or nonlinear fractional differential equations.

So far there have been several fundamental works on the fractional derivative and fractional differential equations [1–3]. These works are to be considered as an introduction to the theory of fractional derivative and fractional differential equations and provide a systematic understanding of the fractional calculus such as the existence and uniqueness [4, 5]. Recently, many other researchers have paid attention to existence result of solution of the initial value problem and boundary problem for fractional differential equations [4–6].

Finding approximate or exact solutions of fractional differential equations is an important task. Except for a limited number of these equations, we have difficulty in finding their analytical solutions. Therefore, there have been attempts to develop new methods for obtaining analytical solutions which reasonably approximate the exact solutions. Several such techniques have drawn special attention, such as Adomain’s decomposition method [7], homotopy perturbation method [8–10], homotopy analysis method [11, 12], variational iteration method [13–17], Chebyshev spectral method [18, 19], and new iterative method [20–22]. Among them, the new iterative method provides an effective procedure for explicit and numerical solutions of a wide and general class of differential systems representing real physical problems. The new iterative method is more superior than the other nonlinear methods, such as the perturbation methods where this method does not depend on small parameters, such that it can find wide application in nonlinear problems without linearization or small perturbation.

The motivation of this paper is to extend the application of the new iterative method proposed by Daftardar-Gejji and Jafari [20–22] to solve linear and nonlinear ordinary and partial differential equations of fractional order. This motivation is based on the importance of these equations and their applications in various subjects in physical branches [10, 11, 14, 23–25].

There are several definitions of a fractional derivative of order [3, 26]. The two most commonly used definitions are Riemann-Liouville and Caputo. Each definition uses Riemann-Liouville fractional integration and derivative of whole order. The difference between the two definitions is in the order of evaluation. Riemann-Liouville fractional integration of order is defined as

The next two equations define Riemann-Liouville and Caputo fractional derivatives of order ,  respectively, as

where ,  .

Caputo fractional derivative first computes an ordinary derivative followed by a fractional integral to achieve the desired order of fractional derivative. Riemann-Liouville fractional derivative is computed in the reverse order. Therefore, Caputo fractional derivative allows traditional initial and boundary conditions to be included in the formulation of the problem.

From properties of   and ,  it is important to note that

where   is Caputo derivative operator of order ,

2. Basic Idea of New Iterative Method

For the basic idea of the new iterative method, we consider the following general functional equation [20–22]:

where is a nonlinear operator from a Banach space and   is a known function. We have been looking for a solution of (4) having the series form

The nonlinear operator can be decomposed as

From (5) and (6), (4) is equivalent to

We define the following recurrence relation:

Then,

If ,  , then

and the series absolutely and uniformly converges to a solution of (4) [27], which is unique, in view of the Banach fixed point theorem [28]. The n-term approximate solution of (4) and (5) is given by

2.1. Convergence of the Method

Now we analyze the convergence of the new iterative method for solving any general functional equation (4). Let , where is the exact solution, is the approximate solution, and is the error in the solution of (4); obviously satisfies (4), that is,

and the recurrence relation (8) becomes

If ,  , then

Thus as , which proves the convergence of the new iterative method for solving the general functional equation (4). For more details, you can see [29].

3. Suitable Algorithm

In this section, we introduce a suitable algorithm for solving nonlinear partial differential equations using the new iterative method. Consider the following nonlinear partial differential equation of arbitrary order:

where is a nonlinear function of and (partial derivatives of with respect to and ) and is the source function. In view of the new iterative method, the initial value problem (14a) and (14b) is equivalent to the integral equation

where

Remark 1. When the general functional equation (4) is linear, the recurrence relation (8) can be simplified in the form

Proof. From the properties of integration and by using (8) and (16b), we have

Therefore, we get the solution of (15) by employing the recurrence relation (8) or (17).

4. Applications

To illustrate the effectiveness of the proposed method, several test examples are carried out in this section.

Example 2. In this example, we consider the following initial value problem in the case of the inhomogeneous Bagely-Torvik equation [23, 24]:

where . The exact solution of this problem is .

By applying the technique described in Sections 2 and 3, the initial value problem (19) is equivalent to the integral equation

Let . In view of recurrence relation (17), we have the following first approximations:

and so on. In the same manner the rest of components can be obtained. The 6-term approximate solution for (19) is

Remark 3. In Example 2. we have used the recurrence relation (17). If we used the recurrence relation (8) in place of (17), we obtain the same result.

In Figure 1, we have plotted the 6-term approximate solution with the corresponding exact solution for (19). It is remarkable to note that the two solutions are almost equal.

Figure 1: Plots of the approximate solution and the exact solution for (19).

Comparing these obtained results with those obtained by new Jacobi operational matrix in [23, 24], we can confirm the simplicity and accuracy of the given method.

Example 4. Consider the following fractional Riccati equation [10]:

The exact solution when is .

By applying the technique described in Sections 2 and 3, the initial value problem (23) is equivalent to the integral equation

Let . In view of recurrence relation (8), we have the following first approximations:

and so on. The 4-term approximate solution for (23) is

In Figure 2, we have plotted the 4-term approximate solution for (23) for different values of with the corresponding exact solution. It is remarkable to note that the approximate solution,  in case ,  and the exact solution are almost equal (continuous curve) whenever the approximate solution, in cases , is of high agreement with the exact solution (dashed and dotted curves, resp.).

Figure 2: Plots of the approximate solution for different values of and the exact solution for (23).

Comparing the obtained results with those obtained by homotopy analysis method, in case , in [10], we can confirm the simplicity and accuracy of the given method.

Example 5. Consider the following initial value problem with fractional order [23, 24]:

The exact solution for this problem is .

As in Example 4, the initial value problem (27) is equivalent to the integral equation

Let . In view of recurrence relation (8), we have the following first approximations:

and so on. The 4-term approximate solution and the corresponding exact solution for (27) are plotted in Figure 3. It is remarkable to note that the two solutions are almost equal.

Figure 3: Plots of the approximate solution and the exact solution for (27).

Comparing these obtained results with those obtained by new Jacobi operational matrix in [23, 24], we can confirm the simplicity and accuracy of the given method.

Example 6. Consider the following fractional order wave equation in 2-dimensional space [14]:

The exact solution for this problem when is

The initial value problem (30) is equivalent to the integral equation

Let .  In view of recurrence relation (17), we have the following first approximations:

and so on. The n-term approximate solution for (30) is

In closed form this gives:

which is the exact solution for the given problem. When , the above n-term approximate solution for (30) becomes

In closed form, this gives

which is the same result obtained by variational iteration method in [14].

Example 7. Consider the following fractional order heat equation in 2-dimensional space [11]:

The exact solution for this problem when is

The initial value problem (38) is equivalent to the integral equation

Let . In view of recurrence relation (17), we have the following first approximations:

and so on. The n-term approximate solution for (38) is

When , The n-term approximate solution for (38) becomes

In closed form, this gives

which is the exact solution for the given problem.

The obtained results in this example are the same as these obtained in [11] by the homotopy perturbation method, in case , but with the simplicity of the given method.

Example 8. In this last example, we consider the following fractional order nonlinear wave equation [25]:

The exact solution for this problem when is where ,  .

The initial value problem (45) is equivalent to the integral equation

Let . In view of recurrence relation (8), we have

and so on. The 3-term approximate solution and the corresponding exact solution for (45) are plotted in Figure 4(a), in case , for  ., in Figure 4(b), in case , for ., and in Figure 4(c), in case . It is remarkable to note that in the first two figures all the solutions are almost equal.

(a)

(b)

(c)

Figure 4: (a) Plots of the approximate solution for different values of and the exact solution, in case ; for (45). (b) Plots of the approximate solution for different values of and the exact solution, in case ; for (45). (c) Plots of the approximate solution, in case for (45).

Comparing these results with those obtained by the modification homotopy perturbation method in [25], we can confirm the accuracy and simplicity of the given method.

5. Conclusion

In this paper, the new iterative method with suitable algorithm is successfully used to solve linear and nonlinear ordinary and partial differential equations with fractional order. It is clear that the computations are easy and the solutions agree well with the corresponding exact solutions and more accurate than the solutions obtained by other methods. Moreover, the accuracy is high with little computed terms of the solution which confirm that this method with the given algorithm is a powerful method for handling fractional differential equations.

References

K. Miller and B. Ross, An Introduction to the Fractional Calculus and Fractional Diffrential Equations, John Wiley & Sons, New York, NY, USA, 1993.

K. B. Oldham and J. Spanier, The Fractional Calculus, Academic Press, London, UK, 1974.

I. Podlubny, Fractional Differential Equations, vol. 198, Academic Press, San Diego, Calif, USA, 1999.

M. Amairi, M. Aoun, S. Najar, and M. N. Abdelkrim, “A constant enclosure method for validating existence and uniqueness of the solution of an initial value problem for a fractional differential equation,” Applied Mathematics and Computation, vol. 217, no. 5, pp. 2162–2168, 2010.

J. Deng and L. Ma, “Existence and uniqueness of solutions of initial value problems for nonlinear fractional differential equations,” Applied Mathematics Letters, vol. 23, no. 6, pp. 676–680, 2010.

E. Girejko, D. Mozyrska, and M. Wyrwas, “A sufficient condition of viability for fractional differential equations with the Caputo derivative,” Journal of Mathematical Analysis and Applications, vol. 381, no. 1, pp. 146–154, 2011.

S. S. Ray and R. K. Bera, “Solution of an extraordinary differential equation by Adomian decomposition method,” Journal of Applied Mathematics, no. 4, pp. 331–338, 2004.

M. Dehghan, J. Manafian, and A. Saadatmandi, “Solving nonlinear fractional partial differential equations using the homotopy analysis method,” Numerical Methods for Partial Differential Equations, vol. 26, no. 2, pp. 448–479, 2010.

I. Hashim, O. Abdulaziz, and S. Momani, “Homotopy analysis method for fractional IVPs,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 3, pp. 674–684, 2009.

Z. Odibat, S. Momani, and H. Xu, “A reliable algorithm of homotopy analysis method for solving nonlinear fractional differential equations,” Applied Mathematical Modelling, vol. 34, no. 3, pp. 593–600, 2010.

A. A. Hemeda, “Homotopy perturbation method for solving partial differential equations of fractional order,” International Journal of Mathematical Analysis, vol. 6, no. 49–52, pp. 2431–2448, 2012.

A. A. Hemeda, “Homotopy perturbation method for solving systems of nonlinear coupled equations,” Applied Mathematical Sciences, vol. 6, no. 93–96, pp. 4787–4800, 2012.

S. Esmaeili, M. Shamsi, and Y. Luchko, “Numerical solution of fractional differential equations with a collocation method based on Müntz polynomials,” Computers & Mathematics with Applications, vol. 62, no. 3, pp. 918–929, 2011.

A. A. Hemeda, “Variational iteration method for solving wave equation,” Computers & Mathematics with Applications, vol. 56, no. 8, pp. 1948–1953, 2008.

A. A. Hemeda, “Variational iteration method for solving non-linear partial differential equations,” Chaos, Solitons and Fractals, vol. 39, no. 3, pp. 1297–1303, 2009.

A. A. Hemeda, “Variational iteration method for solving nonlinear coupled equations in 2-dimensional space in fluid mechanics,” International Journal of Contemporary Mathematical Sciences, vol. 7, no. 37–40, pp. 1839–1852, 2012.

M. G. Sakar, F. Erdogan, and A. Yıldırım, “Variational iteration method for the time-fractional Fornberg-Whitham equation,” Computers & Mathematics with Applications, vol. 63, no. 9, pp. 1382–1388, 2012.

E. H. Doha, A. H. Bhrawy, and S. S. Ezz-Eldien, “Efficient Chebyshev spectral methods for solving multi-term fractional orders differential equations,” Applied Mathematical Modelling, vol. 35, no. 12, pp. 5662–5672, 2011.

E. H. Doha, A. H. Bhrawy, and S. S. Ezz-Eldien, “A Chebyshev spectral method based on operational matrix for initial and boundary value problems of fractional order,” Computers & Mathematics with Applications, vol. 62, no. 5, pp. 2364–2373, 2011.

V. Daftardar-Gejji and H. Jafari, “An iterative method for solving nonlinear functional equations,” Journal of Mathematical Analysis and Applications, vol. 316, no. 2, pp. 753–763, 2006.

A. A. Hemeda, “New iterative method: application to nth-order integro-differential equations,” International Mathematical Forum, vol. 7, no. 47, pp. 2317–2332, 2012.

A. A. Hemeda, “Formulation and solution of th-order derivative fuzzy integrodifferential equation using new iterative method with a reliable algorithm,” Journal of Applied Mathematics, vol. 2012, Article ID 325473, 17 pages, 2012.

A. Saadatmandi and M. Dehghan, “A new operational matrix for solving fractional-order differential equations,” Computers & Mathematics with Applications, vol. 59, no. 3, pp. 1326–1336, 2010.

E. H. Doha, A. H. Bhrawy, and S. S. Ezz-Eldien, “A new Jacobi operational matrix: an application for solving fractional differential equations,” Applied Mathematical Modelling, vol. 36, no. 10, pp. 4931–4943, 2012.

B. Ghazanfari, A. G. Ghazanfari, and M. Fuladvand, “Modification of the homotopy perturbation method for numerical solution of nonlinear wave and system of nonlinear wave equations,” The Journal of Mathematics and Computer Science, vol. 3, no. 2, pp. 212–224, 2011.

M. Caputo, “Linear methods of dissipation whose Q is almost frequency independent, part II,” Journal of the Royal Society of Medicine, vol. 13, pp. 529–539, 1967.

Y. Cherruault, “Convergence of Adomian's method,” Kybernetes, vol. 18, no. 2, pp. 31–38, 1989.

A. J. Jerri, Introduction to Integral Equations with Applications, Wiley-Interscience, New York, NY, USA, 2nd edition, 1999.

S. Bhalekar and V. Daftardar-Gejji, “Convergence of the new iterative method,” International Journal of Differential Equations, vol. 2011, Article ID 989065, 10 pages, 2011.