Abstract and Applied Analysis

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

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

Research Article

Solution of Fractional Partial Differential Equations in Fluid Mechanics by Extension of Some Iterative Method

A. A. Hemeda

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

Received 23 September 2013; Accepted 7 November 2013

Academic Editor: Tie-cheng Xia

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.

An extension of the so-called new iterative method (NIM) has been used to handle linear and nonlinear fractional partial differential equations. The main property of the method lies in its flexibility and ability to solve nonlinear equations accurately and conveniently. Therefore, a general framework of the NIM is presented for analytical treatment of fractional partial differential equations in fluid mechanics. The fractional derivatives are described in the Caputo sense. Numerical illustrations that include the fractional wave equation, fractional Burgers equation, fractional KdV equation, fractional Klein-Gordon equation, and fractional Boussinesq-like equation are investigated to show the pertinent features of the technique. Comparison of the results obtained by the NIM with those obtained by both Adomian decomposition method (ADM) and the variational iteration method (VIM) reveals that the NIM is very effective and convenient. The basic idea described in this paper is expected to be further employed to solve other similar linear and nonlinear problems in fractional calculus.

1. Introduction

Recent advances of fractional differential equations are stimulated by new examples of applications in fluid mechanics, viscoelasticity, mathematical biology, electrochemistry, and physics. For example, the nonlinear oscillation of earthquake can be modeled with fractional derivatives [1], and the fluid-dynamic traffic model with fractional derivatives [2] can eliminate the deficiency arising from the assumption of continuum traffic flow. Based on experimental data fractional partial differential equations for seepage flow in porous media are suggested in [3], and differential equations with fractional order have recently proved to be valuable tools to the modeling of many physical phenomena [4]. Different fractional partial differential equations have been studied and solved including the space-time-fractional diffusion-wave equation [5–7], the fractional advection-dispersion equation [8, 9], the fractional telegraph equation [10], the fractional KdV equation [11], and the linear inhomogeneous fractional partial differential equations [12].

The NIM [13–15] is a suitable approach to provide analytical approximation to linear and nonlinear problems and it is particularly valuable as tool for scientists and applied mathematicians, because it provides immediate and visible symbolic terms of analytical solutions, as well as numerical approximate solutions to both linear and nonlinear differential equations without linearization or discretization. The NIM, proposed by Daftardar-Gejji and Jafari in 2006 [13] and improved by Hemeda [14], was successfully applied to a variety of linear and nonlinear equations such as algebraic equations, integral equations, integrodifferential equations, ordinary and partial differential equations of integer and fractional order, and systems of equations as well. NIM is simple to understand and easy to implement using computer packages and yields better results [15] than the existing ADM [16], homotopy perturbation method (HPM) [17], or VIM [18].

The objective of this work is to extend the application of the NIM to obtain analytical solutions to some fractional partial differential equations in fluid mechanics. These equations include wave equation, Burgers equation, KdV equation, Klein-Gordon equation, and Boussinesq-like equation. The NIM is a computational method that yields analytical solutions and has certain advantages over standard numerical methods. It is free from rounding-off errors as it does not involve discretization and does not require large computer obtained memory or power. The method introduces the solution in the form of a convergent fractional series with elegantly computable terms. The corresponding solutions of the integer order equations are found to follow as special cases of those of fractional order equations. The obtained results of the NIM are compared with those obtained by both ADM [19, 20] and the VIM [21–25] which confirm that this method is very effective and convenient to these equations and to other similar equations where it has the advantage that there is no need to calculate Adomian’s polynomials for the nonlinear problems as in the ADM. For more details see [26–31].

Throughout this work, fractional partial differential equations are obtained from the corresponding integer order equations by replacing the first-order or the second-order time derivative by a fractional in the Caputo sense [32] of order with or .

2. Preliminaries and Notations

In this section we give some basic definitions and properties of the fractional calculus theory which are used further in this work.

Definition 1. A real function , , is said to be in the space , if there exists a real number (>), such that , where , and it is said to be in the space if and only if , .

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

Properties of the operators can be found in [32–35]; we mention only the following: for , , , and ,

(1), (2), (3).

The Riemann-Liouville derivative has certain disadvantages when trying to model real-world phenomena with fractional differential equations. Therefore, we will introduce a modified fractional differential operator proposed by Caputo in this work on the theory of viscoelasticity [32].

Definition 3. The fractional derivative of in the Caputo sense is defined as

Also, we need here two of its basic properties.

Lemma 4. If , , and , , then

The Caputo fractional derivative is considered here because it allows traditional initial and boundary conditions to be included in the formulation of the problem [36]. In this work, we consider the one-dimensional linear inhomogeneous fractional partial differential equations in fluid mechanics, where the unknown function is assumed to be a causal function of time, that is, vanishing for . The fractional derivative is taken in Caputo sense as follows.

Definition 5. For to be the smallest integer that exceeds , the Caputo time-fractional derivative operator of order is defined as

3. New Iterative Method (NIM)

To illustrate the basic idea of the NIM, consider the following general functional equation [13, 14, 26, 27]:

where is a nonlinear operator from a Banach space and is a known function (element) of the Banach space . We are looking for a solution of (5) having the series form

The nonlinear operator can be decomposed as

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

The required solution for (5) can be obtained recurrently from the recurrence relation:

Then

The n-term approximate solution of (5) and (6) is given by

Remark 6. If is a contraction, that is, , , then

Proof. From (9), we have

and the series absolutely and uniformly converges to a solution of (5) [37], which is unique in view of the Banach fixed point theorem [38]. For more details about NIM see [39].

3.1. Reliable Algorithm

To illustrate the basic idea of the reliable algorithm, we consider the general fractional partial differential equation of arbitrary fractional order:

with the initial conditions

where is a nonlinear function of , (partial derivatives of with respect to and ) and is the source function. In view of the fractional calculus and the properties of the fractional integral operators, the initial value problem (14a) and (14b) is equivalent to the following fractional integral equation:

where and . We get the solution of (15) by employing the recurrence relation (9).

Remark 7. When the general functional equation (5) is linear, the recurrence relation (9) can be simplified in the form

Proof. From the properties of integration, in case is an integral operator, we have

In case is a differential operator, we obtain the same result.

3.2. Convergence of NIM

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

and the recurrence relation (9) becomes

If , , then

Thus as , which proves the convergence of the NIM for solving the general functional equation (5).

Remark 8. For linear problems, we get the solution of (5) by employing the recurrence relation (16) in place of the recurrence relation (9).

4. Numerical Examples

4.1. Linear Problems

To incorporate the above discussion, three linear fractional partial differential equations will be studied. The NIM is used to obtain the exact solution of these problems.

Example 1. Consider the following one-dimensional linear inhomogeneous fractional wave equation:

subject to the initial condition

Problem (21a) and (21b) is solved in [40] by using the ADM; the first few components of solution are as follows:

The solution for (21a) and (21b) in series form is given by

Canceling the noise terms and keeping the nonnoise terms in (23) yield the exact solution of (21a) and (21b) given by which is easily verified.

Also, problem (21a) and (21b) is solved in [40] by using the VIM. By beginning with , the following approximations can be obtained:

As in (23), canceling the noise terms and keeping the nonnoise terms yield the exact solution of (21a) and (21b).

According to the NIM and by (15) and (16), we obtain.

Therefore, the initial value problem (21a) and (21b) is equivalent to the integral equation:

Let ; we can obtain the following first few components of the new iterative solution for (21a) and (21b):

The solution in series form is given by

Canceling the noise terms and keeping the non-noise terms in (28) yield the exact solution of (21a) and (21b).

From (23), (24), and (28), it is clear that the three methods are the same in solving (21a) and (21b).

Example 2. Consider the following one-dimensional linear inhomogeneous fractional Burgers equation

subject to the initial condition

Problem (29a) and (29b) is solved in [40] by using the ADM; the first few components of solution are as follows:

Therefore, the exact solution is given by .

Also, problem (29a) and (29b) is solved in [40] by using the VIM with ; the following approximations can be obtained:

The exact solution follows immediately.

According to the NIM, by (15) and (16), we obtain

Therefore, the initial value problem (29a) and (29b) is equivalent to the integral equation:

Let ; we obtain the following first few components of the new iterative solution for (29a) and (29b):

Therefore, the exact solution follows immediately.

Example 3. Consider the following one-dimensional linear inhomogeneous fractional Klein-Gordon equation:

subject to the initial conditions

In view of the ADM [40], the first few components of solution for (35a) and (35b) are derived as follows:

The solution in series form is given by

In view of the VIM [40], with , the following approximations for (35a) and (35b) are obtained:

According to the NIM, by (15) and (16), we can obtain

Therefore, the initial value problem (35a) and (35b) is equivalent to the integral equation:

Let ; we can obtain the following first few components of the new iterative solution for (35a) and (35b):

The solution for (35a) and (35b) in series form is given by

From (37), (38), and (42), the ADM, the VIM, and the NIM give the same solution for the classical Klein-Gordon equation (35a) and (35b) in the case which is given by

Canceling the noise terms and keeping the non-noise terms in (43) yield the exact solution of (35a) and (35b) for the special case , which is given by , which is easily verified.

4.2. Nonlinear Problems

For nonlinear equations in general, there exists no method that yields the exact solution and therefore only approximate solutions can be derived. In this subsection, we use the NIM to provide approximate solutions for two kinds of nonlinear time-fractional partial differential equations.

Example 4. Consider the following one-dimensional nonlinear homogeneous time-fractional KdV equation:

subject to the initial condition

The time-fractional KdV equation (44a) and (44b) is solved in [11, 40] by using the ADM. The solution in series form is found as

where

Also, problem (44a) and (44b) is solved in [40] by using the VIM with ; the first few approximate solutions are

According to the NIM, by (15) and (9), we can obtain

Therefore, the initial value problem (44a) and (44b) is equivalent to the integral equation:

Let ; we obtain the following first few components of the new iterative solution for (44a) and (44b):

The solution in series form is given by

where

Now, the approximate solution for (44a) and (44b) obtained by the NIM in (51) is the same solution obtained by both the ADM in (45) and the VIM in (47). Therefore, all the three methods provide the same approximate solution for the time-fractional KdV equation.

Example 5. In this example, we consider the one-dimensional nonlinear homogeneous time-fractional Boussinesq-like equation:

subject to the initial conditions

The time-fractional Boussinesq-like equation (53a) and (53b) is solved by the ADM in [40]. The fourth-order term approximate solution in series form is given by

Also, problem (53a) and (53b) is solved by the VIM in [40]. The fourth-order term approximate solution, with , is given by

According to the NIM, by (15) and (9), we can obtain the following first approximations:

The fourth-order term approximate solution in series form is given by

where

It is interesting to point out that, for the case of , the approximate solution:

follows immediately upon replacing by in the decomposition solution (54), the variational iteration solution (55), or the new iteration solution (57), which converges to the exact solution of the Boussinesq-like equation (53a) and (53b), when , .

The advantage of the NIM in the nonlinear problems, Examples 4 and 5, is that there is no need to calculate Adomian’s polynomials as done in the ADM, which means that the first method is simple and easy of the procedure of calculations over the second. Also, for more details about the solutions of the fractional differential equations by the VIM and some new asymptotic methods, you can see [41–43].

5. Conclusion

NIM has been known as a powerful tool for solving many functional equations such as ordinary, partial differential equations, integral equations, integrodifferential equations, and so many other equations. In this work, we have presented a general framework of the NIM for the analytical treatment of fractional partial differential equations in fluid mechanics. The present work shows the validity and great potential of the NIM for solving linear and nonlinear fractional partial differential equations. All of the examples show that the results of the NIM are in excellent agreement with those obtained by both the ADM and the VIM without calculating Adomian’s polynomials for the nonlinear problems as in the ADM. The basic idea described in this work is expected to be further employed to solve other similar linear and nonlinear problems in fractional calculus.

Conflict of Interests

The author declared that there is no conflict of interests regarding the publication of this paper.

References

J. H. He, “Nonlinear oscillation with fractional derivative and its applications,” in Proceedings of the International Conference on Vibrating Engineering 98, Dalian, China, 1988.

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.

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.

I. Podlubny, Fractional Differential Equations, Academic Press, New York, NY, USA, 1999.

K. Al-Khaled and S. Momani, “An approximate solution for a fractional diffusion-wave equation using the decomposition method,” Applied Mathematics and Computation, vol. 165, no. 2, pp. 473–483, 2005.

F. Mainardi, Y. Luchko, and G. Pagnini, “The fundamental solution of the space-time fractional diffusion equation,” Fractional Calculus & Applied Analysis, vol. 4, no. 2, pp. 153–192, 2001.

A. Hanyga, “Multidimensional solutions of time-fractional diffusion-wave equations,” The Royal Society of London A, vol. 458, no. 2020, pp. 933–957, 2002.

F. Huang and F. Liu, “The time fractional diffusion equation and the advection-dispersion equation,” The ANZIAM Journal, vol. 46, no. 3, pp. 317–330, 2005.

F. Huang and F. Liu, “The fundamental solution of the space-time fractional advection-dispersion equation,” Journal of Applied Mathematics & Computing, vol. 18, no. 1-2, pp. 339–350, 2005.

S. Momani, “Analytic and approximate solutions of the space-and time-fractional telegraph equations,” Applied Mathematics and Computation, vol. 170, no. 2, pp. 1126–1134, 2005.

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.

L. Debnath and D. D. Bhatta, “Solutions to few linear fractional inhomogeneous partial differential equations in fluid mechanics,” Fractional Calculus & Applied Analysis, vol. 7, no. 1, pp. 21–36, 2004.

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 the nth-order integro-differential equations,” Information B, vol. 16, no. 6, pp. 3841–3852, 2013.

S. Bhalekar and V. Daftardar-Gejji, “New iterative method: application to partial differential equations,” Applied Mathematics and Computation, vol. 203, no. 2, pp. 778–783, 2008.

G. Adomian, Solving Frontier Problems of Physics: The Decomposition Method, Kluwer, 1994.

J.-H. He, “Homotopy perturbation technique,” Computer Methods in Applied Mechanics and Engineering, vol. 178, no. 3-4, pp. 257–262, 1999.

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, pp. 1839–1852, 2012.

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.

T. Öziş and A. Yıldırım, “Comparison between Adomian's method and He's homotopy perturbation method,” Computers & Mathematics with Applications, vol. 56, no. 5, pp. 1216–1224, 2008.

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.

J.-H. He, “Variational principle for nano thin film lubrication,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 4, no. 3, pp. 313–314, 2003.

J.-H. He, “Variational principles for some nonlinear partial differential equations with variable coefficients,” Chaos, Solitons & Fractals, vol. 19, no. 4, pp. 847–851, 2004.

J.-H. He, “Variational iteration method—some recent results and new interpretations,” Journal of Computational and Applied Mathematics, vol. 207, no. 1, pp. 3–17, 2007.

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.

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

S. Bhalekar and V. Daftardar-Gejji, “Solving evolution equations using a new iterative method,” Numerical Methods for Partial Differential Equations, vol. 26, no. 4, pp. 906–916, 2010.

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 nonlinear partial differential equations,” Chaos, Solitons and Fractals, vol. 39, no. 3, pp. 1297–1303, 2009.

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.

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.

M. Caputo, “Linear models of dissipation whose Q is almost frequency indepedent—part II,” Geophysical Journal International, vol. 13, no. 5, pp. 529–539, 1967.

A. Y. Luchko and R. Groreo, The Initial Value Problem for Some Fractional Differential Equations with Caputo Derivative, Preprint series A80-98, Fachbreich Mathematik und Informatik, Freic Universitat Berlin, 1998.

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

K. B. Oldham and J. Spanier, The Fractional Calculus, Academic Press, New York, NY, USA, 1974.

I. Podlubny, “Geometric and physical interpretation of fractional integration and fractional differentiation,” Fractional Calculus & Applied Analysis, vol. 5, no. 4, pp. 367–386, 2002.

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, 2nd edition, 1999.

A. A. Hemeda, “New iterative method: an application for solving fractional physical differential equations,” Abstract and Applied Analysis, vol. 2013, Article ID 617010, 9 pages, 2013.

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.

G.-C. Wu and D. Baleanu, “Variational iteration method for the Burgers' flow with fractional derivatives—new Lagrange multipliers,” Applied Mathematical Modelling, vol. 37, no. 9, pp. 6183–6190, 2013.

G. C. Wu, “New trends in the variational iteration method,” Communications in Fractional Calculus, vol. 2, pp. 59–75, 2011.

J.-H. He, “Asymptotic methods for solitary solutions and compactons,” Abstract and Applied Analysis, vol. 2012, Article ID 916793, 130 pages, 2012.