Mathematical Problems in Engineering

Volume 2010 (2010), Article ID 318298, 11 pages

http://dx.doi.org/10.1155/2010/318298

## Modified Variational Iteration Method for Free-Convective Boundary-Layer Equation Using Padé Approximation

^{1}HITEC University, Taxila Cantt, Pakistan^{2}Department of Mathematics, Ege University, Bornova 35100, Izmir, Turkey^{3}Department of Mathematics, University of Dayton, 300 College Park, Dayton, OH 45469-2316, USA

Received 16 October 2009; Accepted 17 January 2010

Academic Editor: Jihuan He

Copyright © 2010 Syed Tauseef Mohyud-Din 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.

#### Abstract

This paper is devoted to the study of a free-convective boundary-layer flow modeled by a system of nonlinear ordinary differential equations. We apply a modified variational iteration method (MVIM) coupled with He's polynomials and Padé approximation to solve free-convective boundary-layer equation. It is observed that the combination of MVIM and the Padé approximation improves the accuracy and enlarges the convergence domain.

#### 1. Introduction

The boundary-layer flows of viscous fluids are of utmost importance for industry and applied sciences. These flows can be modeled by systems of nonlinear ordinary differential equations on an unbounded domain, see [1–4] and the references therein. Keeping in view the physical importance of such problems, there is a dire need of extension of some reliable and efficient technique for the solution of such problems. He [1, 2, 5–15] developed the variational iteration (VIM) and homotopy perturbation (HPM) methods which are very efficient and accurate and are [1, 2, 4–42] being used very frequently for finding the appropriate solutions of nonlinear problems of physical nature. In a later work, Ghorbani and Nadjfi [24] introduced He’s polynomials which are calculated for He’s homotopy perturbation method. It is also established [24] that He’s polynomials are compatible with Adomian’s polynomials but are easier to implement and are more user friendly. Recently, Mohyud-Din, Noor and Noor [4, 33–36] made the elegant coupling of He’s polynomials and the correction functional of variational iteration method (VIM) and found the solutions of number of nonlinear singular and nonsingular problems. It is observed that [4, 33–36] the modified version of VIM is very efficient in solving nonlinear problems. The basic motivation of this paper is the extension of the modified variational iteration method (MVIM) coupled with Padé approximation to solve a free-convective boundary-layer flow modeled by a system of nonlinear ordinary differential equations. Numerical and figurative illustrations show that it is a promising tool to solve nonlinear problems. It needs to be highlighted that Herisanu and Marinca [41] suggested an optimal variational iteration algorithm. It needs to be highlighted that He in his latest article “The variational iteration method which should be followed” [42] presented a very comprehensive and detailed study on various aspects of variational iteration method in connection with partial differential equations, ordinary differential equations, fractional differential equations, fractal-differential equations, and difference-differential equations.

#### 2. Modified Variational Iteration Method (MVIM)

To illustrate the basic concept of the modified variational iteration method (MVIM), we consider the following general differential equation:

where is a linear operator, is a nonlinear operator, and is the forcing term. According to variational iteration method [1, 2, 4, 10–23, 28, 33–39, 41, 42], we can construct a correction functional as follows:

where is a Lagrange multiplier [1, 2, 10–15, 42], which can be identified optimally via variational iteration method. The subscripts denote the th approximation; is considered as a restricted variation. That is, (2.2) is called a correction functional. Now, we apply He’s polynomials [24]

which is the coupling of variational iteration method and He’s polynomials and is called the modified variational iteration method (MVIM) [4, 33–36]. The comparison of like powers of gives solutions of various orders.

#### 3. Mathematical Model

Let us consider the problem of cooling of a low-heat-resistance sheet that moves downwards in a viscous fluid [3]:

subject to

where and are the velocity components in the *-* and *-*directions, respectively.is the temperature, is the temperature of the surrounding fluid, is the kinematic viscosity, is the thermal diffusivity, is the acceleration due to gravity, and is the coefficient of thermal expansion. Using the similarity variables

where is the stream function defined by and , and are the similarity functions dependent on and , (3.1) is transformed to

subject to the boundary conditions

where the primes denote differentiation with respect to , and is the Prandtl number.

#### 4. The Padé Approximation

We denote Padé approximants to by

where is polynomial of degree at most and is a polynomial of degree at most *. *The former power series is

And we write the and as

so

and the coefficients of and are determined by the equation. From (4.4), we have

which system of homogeneous equations with unknown quantities. We impose the normalization condition

We can write out (4.5) as

From (4.7), we can obtain . Once the values of are all known (4.8) gives an explicit formula for the unknown quantities . For the diagonal approximants like have the most accurate approximants by built-in utilities of Maple.

#### 5. Solution Procedure

Consider problems (3.4)–(3.5) formulated in Section 3 and is related to the free-convective boundary-layer flow.

The correction functional is given by

Making the correction functional stationary, the Lagrange multipliers can easily be identified

Consequently,

Applying the modified variational iteration method (MVIM), we get

Comparing the coefficient of like powers of , we get

The series solution is given by

It is observed in Figures 1 and 2 that the flow has a boundary-layer structure and the thickness of this boundary-layer decreases with increase in the Prandtl number, as expected. This is due to the inhibiting influence of the viscous forces.

Figure 3 shows the increase of the Prandtl number, * σ,* that results in the decrease, as expected, of temperature distribution at a particular point of the flow region, that is, there would be a decrease of the thermal boundary-layer thickness with the increase of values of

*implying a slow rate of thermal diffusion. Thus higher Prandtl number*

*σ**leads to faster cooling of the plane sheet.*

*σ*#### 6. Conclusions

In this study, we employed modified variational iteration method (MVIM) coupled with Padé approximation to solve a system of two nonlinear ordinary differential equations that describes a free-convective boundary-layer in glass-fiber production process. The results show strong effects of the Prandtl number on the velocity and temperature profiles since the two model equations are coupled.

#### Acknowledgment

The authors are highly grateful to the referee for his/her very constructive comments.

#### References

- J. H. He, “An elementary introduction of recently developed asymptotic methods and nanomechanics in textile engineering,”
*International Journal of Modern Physics B*, vol. 22, no. 21, pp. 3487–4578, 2008. View at Publisher · View at Google Scholar - J.-H. He, “Some asymptotic methods for strongly nonlinear equations,”
*International Journal of Modern Physics B*, vol. 20, no. 10, pp. 1141–1199, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - S. A. Kechil and I. Hashim, “Non-perturbative solution of free-convective boundary-layer equation by adomian decomposition method,”
*Physics Letters A*, vol. 363, pp. 110–114, 2007. View at Publisher · View at Google Scholar - S. T. Mohyud-Din, M. A. Noor, and K. I. Noor, “Some relatively new techniques for nonlinear problems,”
*Mathematical Problems in Engineering*, vol. 2009, Article ID 234849, 25 pages, 2009. View at Google Scholar · View at MathSciNet - J.-H. He, “Recent development of the homotopy perturbation method,”
*Topological Methods in Nonlinear Analysis*, vol. 31, no. 2, pp. 205–209, 2008. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - J.-H. He, “Comparison of homotopy perturbation method and homotopy analysis method,”
*Applied Mathematics and Computation*, vol. 156, no. 2, pp. 527–539, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - J. H. He, “Homotopy perturbation method for bifurcation of nonlinear problems,”
*International Journal of Nonlinear Science and Numerical Simulation*, vol. 6, no. 2, pp. 207–208, 2005. View at Google Scholar - J.-H. He, “The homotopy perturbation method nonlinear oscillators with discontinuities,”
*Applied Mathematics and Computation*, vol. 151, no. 1, pp. 287–292, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - J.-H. He, “A 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 · View at Google Scholar · View at MathSciNet - 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. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - 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. View at Publisher · View at Google Scholar - J.-H. He, “Variational iteration method for autonomous ordinary differential systems,”
*Applied Mathematics and Computation*, vol. 114, no. 2-3, pp. 115–123, 2000. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - J.-H. He, “The variational iteration method for eighth-order initial-boundary value problems,”
*Physica Scripta*, vol. 76, no. 6, pp. 680–682, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - 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. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - J.-H. He and X.-H. Wu, “Construction of solitary solution and compacton-like solution by variational iteration method,”
*Chaos, Solitons & Fractals*, vol. 29, no. 1, pp. 108–113, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - S. Abbasbandy, “A new application of He's variational iteration method for quadratic Riccati differential equation by using Adomian's polynomials,”
*Journal of Computational and Applied Mathematics*, vol. 207, no. 1, pp. 59–63, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - S. Abbasbandy, “Numerical solution of non-linear Klein-Gordon equations by variational iteration method,”
*International Journal for Numerical Methods in Engineering*, vol. 70, no. 7, pp. 876–881, 2007. View at Publisher · View at Google Scholar · View at MathSciNet - M. A. Abdou and A. A. Soliman, “Variational iteration method for solving Burger's and coupled Burger's equations,”
*Journal of Computational and Applied Mathematics*, vol. 181, no. 2, pp. 245–251, 2005. View at Publisher · View at Google Scholar · View at MathSciNet - M. A. Abdou and A. A. Soliman, “New applications of variational iteration method,”
*Physica D*, vol. 211, no. 1-2, pp. 1–8, 2005. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - T. A. Abassy, M. A. El-Tawil, and H. El Zoheiry, “Solving nonlinear partial differential equations using the modified variational iteration Padé technique,”
*Journal of Computational and Applied Mathematics*, vol. 207, no. 1, pp. 73–91, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - B. Batiha, M. S. M. Noorani, and I. Hashim, “Variational iteration method for solving multispecies Lotka-Volterra equations,”
*Computers & Mathematics with Applications*, vol. 54, no. 7-8, pp. 903–909, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - J. Biazar and H. Ghazvini, “He's variational iteration method for fourth-order parabolic equations,”
*Computers & Mathematics with Applications*, vol. 54, no. 7-8, pp. 1047–1054, 2007. View at Publisher · View at Google Scholar · View at MathSciNet - D. D. Ganji, H. Tari, and M. B. Jooybari, “Variational iteration method and homotopy perturbation method for nonlinear evolution equations,”
*Computers & Mathematics with Applications*, vol. 54, no. 7-8, pp. 1018–1027, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - A. Ghorbani and J. S. Nadjfi, “He's homotopy perturbation method for calculating Adomian's polynomials,”
*International Journal of Nonlinear Science and Numerical Simulation*, vol. 8, no. 2, pp. 229–332, 2007. View at Google Scholar - M. Javidi and A. Golbabai, “A variational iteration method for solving parabolic partial differential equations,”
*Computers & Mathematics with Applications*, vol. 54, no. 7-8, pp. 1003–1009, 2007. View at Google Scholar - J. Lu, “Variational iteration method for solving two-point boundary value problems,”
*Journal of Computational and Applied Mathematics*, vol. 207, no. 1, pp. 92–95, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - W.-X. Ma, H. Y. Wu, and J. S. He, “Partial differential equations possessing Frobenius integrable decompositions,”
*Physics Letters A*, vol. 364, no. 1, pp. 29–32, 2007. View at Publisher · View at Google Scholar · View at MathSciNet - S. Momani and Z. Odibat, “Application of He's variational iteration method to Helmholtz equation,”
*Chaos, Solitons & Fractals*, vol. 27, no. 5, pp. 1119–1123, 2006. View at Google Scholar - T. Mohyud-Din, M. A. Noor, and K. I. Noor, “Travelling wave solutions of seventh-order generalized KdV equations using He's polynomials,”
*International Journal of Nonlinear Sciences and Numerical Simulation*, vol. 10, no. 2, pp. 223–229, 2009. View at Google Scholar - S. T. Mohyud-Din and M. A. Noor, “Homotopy perturbation method for solving fourth-order boundary value problems,”
*Mathematical Problems in Engineering*, vol. 2007, Article ID 98602, 15 pages, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - S. T. Mohyud-Din and M. A. Noor, “Homotopy perturbation method for solving partial differential equations,”
*Zeitschrift für Naturforschung A*, vol. 64, pp. 157–170, 2009. View at Google Scholar - S. T. Mohyud-Din, M. A. Noor, and K. I. Noor, “Parameter-expansion techniques for strongly nonlinear oscillators,”
*International Journal of Nonlinear Sciences and Numerical Simulation*, vol. 10, no. 5, pp. 581–583, 2009. View at Google Scholar - S. T. Mohyud-Din, M. A. Noor, and K. I. Noor, “On the coupling of polynomials with correction functional,”
*International Journal of Modern Physics B*. In press. - S. T. Mohyud-Din, “Solving heat and wave-like equations using He's polynomials,”
*Mathematical Problems in Engineering*, vol. 2009, Article ID 427516, 12 pages, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. A. Noor and S. T. Mohyud-Din, “Variational iteration method for solving higher-order nonlinear boundary value problems using He's polynomials,”
*International Journal of Nonlinear Sciences and Numerical Simulation*, vol. 9, no. 2, pp. 141–157, 2008. View at Google Scholar - M. A. Noor and S. T. Mohyud-Din, “Variational homotopy perturbation method for solving higher dimensional initial boundary value problems,”
*Mathematical Problems in Engineering*, vol. 2008, Article ID 696734, 11 pages, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. Rafei and H. Daniali, “Application of the variational iteration method to the Whitham-Broer-Kaup equations,”
*Computers & Mathematics with Applications*, vol. 54, no. 7-8, pp. 1079–1085, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - N. H. Sweilam, “Harmonic wave generation in non linear thermoelasticity by variational iteration method and Adomian's method,”
*Journal of Computational and Applied Mathematics*, vol. 207, no. 1, pp. 64–72, 2007. View at Publisher · View at Google Scholar · View at MathSciNet - M. Tatari and M. Dehghan, “On the convergence of He's variational iteration method,”
*Journal of Computational and Applied Mathematics*, vol. 207, no. 1, pp. 121–128, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - L. Xu, “He's homotopy perturbation method for a boundary layer equation in unbounded domain,”
*Computers & Mathematics with Applications*, vol. 54, no. 7-8, pp. 1067–1070, 2007. View at Publisher · View at Google Scholar · View at MathSciNet - N. Herisanu and V. Marinca, “A modified variational iteration method for strongly nonlinear problems,”
*Nonlinear Science Letters A*, vol. 1, pp. 183–192, 2010. View at Google Scholar - J. H. He, G. C. Wu, and F. Austin, “The variational iteration method which should be followed,”
*Nonlinear Science Letters A*, vol. 1, pp. 1–30, 2010. View at Google Scholar