Research Article  Open Access
M. Aghakhani, M. Suhatril, M. Mohammadhassani, M. Daie, A. Toghroli, "A Simple Modification of Homotopy Perturbation Method for the Solution of Blasius Equation in SemiInfinite Domains", Mathematical Problems in Engineering, vol. 2015, Article ID 671527, 7 pages, 2015. https://doi.org/10.1155/2015/671527
A Simple Modification of Homotopy Perturbation Method for the Solution of Blasius Equation in SemiInfinite Domains
Abstract
A simple modification of the homotopy perturbation method is proposed for the solution of the Blasius equation with two different boundary conditions. Padé approximate is used to deal with the boundary condition at infinity. The results obtained from the analytical method are compared to Howarth’s numerical solution and fifth order RungeKutta Fehlberg method indicating a very good agreement. The proposed method is a simple and reliable modification of homotopy perturbation method, which does not require the existence of a small parameter, linearization of the equation, or computation of Adomian’s polynomials.
1. Introduction
Nonlinear differential equations are ubiquitous in science and technology. However, finding analytical solutions for this class of equations always has been a challenging task. In the recent years, several approximate methods were proposed for the analytical solution of nonlinear differential equations that do not depend on the existence of a small or large parameter in the equation. Among them, Homotopy Analysis Method (HAM) [1], Adomian Decomposition Method (ADM) [2], Variational Iteration Method (VIM) [3], Differential Transformation Method (DTM) [4], and Homotopy Perturbation Method (HPM) [5] can be mentioned. The homotopy perturbation method, first proposed by He [5], combines ease of implementation of perturbation methods with the flexibility of the homotopy analysis method. In addition, in the recent years, some novel methods for approximate solution of nonlinear differential equations emerged, such as Optimal Homotopy Asymptotic Method (OHAM) [6], Generalized Homotopy Method (GHM) [7], and expansion method [8].
In the past, many scientists attempted to suggest an improvement to the homotopy perturbation method. Their studies were mainly focused on enlarging convergence radius and accelerating convergence of the solution and new suggestions for homotopy construction as well as alternation of linear and nonlinear part of homotopy based on the applied problem. Jafari and Aminataei [9] introduced a new treatment for homotopy perturbation method, which improves results from HPM. They discussed convergency of the proposed method. In order to demonstrate efficiency, accuracy, and superiority of the suggested method, the new modification was applied to some experiments in their paper. Yusufoğlu [10] purposed an alternation of HPM for exact solution of system of linear equations. He introduced an accelerating parameter to solve the system linear equation. Through embedding the accelerating parameters, convergence improved and required iterations reduced to one iteration. An effective method for convergence improvement of HPM for the solution of fractional differential equations was suggested by Hosseinnia et al. [11]. They proposed a method to select the linear part in the HPM to keep the inherent stability of fractional equations. To illustrate the improvement, Riccati fractional differential equations were solved. Results indicated the accuracy and effectiveness of the method compared to traditional HPM. A modification of HPM to solve fractional multidimensional diffusion equations is introduced by Kumar et al. [12]. In their study, Sumudu transform was utilized for transformation of partial differential equations into a new form. Dong et al. [13] obtained the solution of strongly nonlinear mixed VolterraFredholm integral equation by using an improvement to the homotopy perturbation method. Although the traditional HPM is divergent in these kinds of problems, the proposed method is convergent and leads to the exact solution.
In the present paper, a simple and reliable modification of HPM is employed for the solution of two different forms of nonlinear Blasius equation in a semiinfinite domain. For that, we consider two forms of the Blasius equation arising in fluid flow inside the velocity boundary layer as follows.
The first form of the Blasius equation is as follows: And the second form is as follows:
It can be seen that two equations are the same except for boundary conditions. The first form of the equation is the wellknown classical Blasius first derived by Blasius [14] and dates back about a century, which describes the velocity profile of twodimensional viscous laminar flow over a finite flat plate. This form of the Blasius equation is the simplest form and the origin of all boundary layer equations in fluid mechanics. The second form of the equation, proposed more recently, arises in the steady free convection about a vertical flat plate embedded in a saturated porous medium [15], Laminar boundary layers at the interface of cocurrent parallel streams [16], or the flow near the leading edge of a very long, steadily operating conveyor belt [17]. Despite the fact that the Blasius equation is a simple third order nonlinear differential equation with smooth and monotonic solution, the hope that solution techniques employed to handle this equation can be extended to the difficult equations motivated several researchers to conduct abundant studies around it by means of different analytical and numerical approaches [18–21]. It is worth mentioning that despite that this problem is a century old, it is still being investigated by several authors, and some recent paper was published about it [22–26].
In the recent years, the homotopy perturbation method is successfully applied to a wide variety of problems in science and technology. Analytical solution of the Blasius equation plays a very important role in the design and optimization of fluid devices. In the present paper, a simple modification of this analytical method was employed for the solution of the Blasius equation with two different forms of boundary conditions. For that, first, the equations are transferred to corresponding initial value problems with appropriate boundary conditions. Then, the initial value problems are solved by the homotopy perturbation method. Finally, the diagonal Padé transformation is used to handle the boundary condition at infinity and enlarge the convergence radius of the resulting series. For the convenience of the reader, we first introduce the homotopy perturbation method in the following section.
2. Homotopy Perturbation Method
The combination of the perturbation method and the homotopy method is called the homotopy perturbation method (HPM). This method lacks the limitations of the traditional perturbation methods and can take the full advantage of the traditional perturbation techniques. To illustrate the basic idea of this method [5], consider the following nonlinear differential equation:subjected to the boundary conditions where is a general differential operator, is a boundary operator, is a known analytical function, and is the boundary of the domain. denotes differentiation along the normal drawn outwards from . The operator , generally speaking, can be divided into two parts of and , where is the linear part, while is a nonlinear one. Hence, (3) can be rewritten as follows:We can construct the following homotopy which does not require that (3) has a small parameter:in which is called homotopy parameter and is an initial approximation of (3) satisfying the given conditions. From (3), we haveSolution of (3) can be written as a power series in as follows:Setting results in the approximate solution in the form ofThe convergence of (9) has been proven by He [27].
3. Application of HPMPadé to the First Form of the Blasius Equation
In this section, we consider the following form of the Blasius equation:
Due to the boundary condition at infinity, direct application homotopy perturbation method is not appropriate for this equation and will not yield the desired results. Consequently, a modification in the equation or boundary condition is necessary. Equation (10) is reformulated as follows: where is an unknown to be determined later. This modification of boundary conditions converts (10) from a boundary value problem to the corresponding initial value problem, which can be easily handled by HPM. This procedure is similar to the shooting method for numerical solution of boundary value problems in which the given boundary value problem is converted into an initial value problem.
Following the outline given in Section 2, a homotopy is constructed for (11):It is assumed that (12) has a solution in the following form:Substituting from (13) into (12) and equating terms with identical powers of , we haveand ….
Equations (14) can be rapidly solved, resulting in the following:and ….
Therefore, according to (9), the solution of (12) reads
The unknown initial curvature in (16), , can be determined by imposing the boundary condition of (10) at infinity, that is, For that, Padé approximants of (16) were formed. Padé approximants have the advantage of converting a function into a rational function in order to obtain more information about . A [] Padé approximate to function is the quotient of two polynomials, and of degrees and , respectively. It is well known that if a function is free of singularities on the real axis, the Padé approximant will usually converge on the entire real axis [30]. Following the procedure suggested by Boyd [31], in order to determine the unknown parameter in (16), the boundary condition at infinity was imposed to the diagonal Padé approximant. Then, the roots of Padé approximant are used for the calculation of the unknown value. The result of this calculation is presented in Table 1. It can be seen from this table that the value obtained for agrees very well with the numerical results. In Tables 2–4, , , and obtained from 12th order HPMPadé approximation are compared with Howarth’s numerical solution [28]. Furthermore, as it can be seen from Tables 2–4, purposed modification of HPMPadé method is more accurate than the variational iteration method by He [29] and is valid for a wider range of the solution domain.




4. Application of HPMPadé to the Second Form of the Blasius Equation
The second form of the Blasius equation is as follows:Similar to the previous section, in order to solve (17), this equation was modified as follows:
In this equation, is unknown initial curvature of (17) and will be determined later by the application of the boundary condition of (17) at infinity. According to the instructions in Section 2, by embedding , a homotopy was constructed for (17) as follows:assuming that this equation has a solution in the following form:Substituting the above equation in (19) and rearranging will result in the following:and ….
Ordinary differential equations (21) can be solved easily. The result is as follows:and ….
Consequently, the solution of (18) can be written as follows:
In order to determine the unknown initial curvature in (23), boundary conditions of (17) at infinity, , must be applied. For that, Padé approximants of (23) which enlarge convergence radius of the solution were used. Then, was determined from to the Padé approximants. Initial curvature of (23) obtained from this method is compared to fifth order RungeKutta Fehlberg numerical method in Table 5. It is worth mentioning that examining the behavior of (17) reveals that its initial curvature must be negative. Therefore, the negative root of Padé approximants is selected. In Tables 6–8, the result of 8th order HPMPadé approximation is presented against that of exact (numerical) method. It can be seen that there is a good agreement between the results of the proposed method and numerical solution.




5. Conclusion
In the present paper, a simple modification of the homotopy perturbation method is proposed for the solution of the Blasius equation in semiinfinite domains. The equation in semiinfinite domain is transferred into equivalent initial value problems, which results in appearance of an unknown coefficient. In order to determine the unknown coefficient, the boundary condition of the problem at infinity is imposed to Padé approximant of the solution. The results are in very good agreement with numerical and previous data available in the literature.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgment
The study presented herein was made possible by the University of Malaya Research Grant UMRG RP004D11AET. The authors would like to acknowledge the support.
References
 S. Liao, “On the homotopy analysis method for nonlinear problems,” Applied Mathematics and Computation, vol. 147, no. 2, pp. 499–513, 2004. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 G. Adomian, Solving Frontier Problems of Physics: The Decomposition Methoc [i.e Method], Kluwer Academic Publishers, Dordrecht, The Netherlands, 2013.
 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
 M.J. Jang, C.L. Chen, and Y.C. Liy, “On solving the initialvalue problems using the differential transformation method,” Applied Mathematics and Computation, vol. 115, no. 23, pp. 145–160, 2000. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 J.H. He, “Homotopy perturbation method: a new nonlinear analytical technique,” Applied Mathematics and Computation, vol. 135, no. 1, pp. 73–79, 2003. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 N. Herisanu, V. Marinca, and Gh. Madescu, “An analytical approach to nonlinear dynamical model of a permanent magnet synchronous generator,” Wind Energy, vol. 18, pp. 1657–1670, 2015. View at: Publisher Site  Google Scholar
 H. VazquezLeal, “Generalized homotopy method for solving nonlinear differential equations,” Computational and Applied Mathematics, vol. 33, no. 1, pp. 275–288, 2014. View at: Publisher Site  Google Scholar  MathSciNet
 M. N. Alam, M. A. Akbar, and S. T. MohyudDin, “A novel (G′/G)expansion method and its application to the Boussinesq equation,” Chinese Physics B, vol. 23, no. 2, Article ID 020203, 2014. View at: Publisher Site  Google Scholar
 M. A. Jafari and A. Aminataei, “Improved homotopy perturbation method,” International Mathematical Forum, vol. 5, no. 29–32, pp. 1567–1579, 2010. View at: Google Scholar  MathSciNet
 E. Yusufoğlu, “An improvement to homotopy perturbation method for solving system of linear equations,” Computers & Mathematics with Applications, vol. 58, no. 1112, pp. 2231–2235, 2009. View at: Publisher Site  Google Scholar
 S. H. Hosseinnia, A. Ranjbar, and S. Momani, “Using an enhanced homotopy perturbation method in fractional differential equations via deforming the linear part,” Computers and Mathematics with Applications, vol. 56, no. 12, pp. 3138–3149, 2008. View at: Publisher Site  Google Scholar
 D. Kumar, J. Singh, and S. Kumar, “Numerical computation of fractional multidimensional diffusion equations by using a modified homotopy perturbation method,” Journal of the Association of Arab Universities for Basic and Applied Sciences, vol. 17, pp. 20–26, 2015. View at: Publisher Site  Google Scholar
 C. Dong, Z. Chen, and W. Jiang, “A modified homotopy perturbation method for solving the nonlinear mixed VolterraFredholm integral equation,” Journal of Computational and Applied Mathematics, vol. 239, pp. 359–366, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 H. Blasius, “Grenzschichten in Flussigkeiten mit kleiner Reibung,” Zeitschrift für Angewandte Mathematik und Physik, vol. 56, pp. 1–37, 1908. View at: Google Scholar
 P. Cheng and W. J. Minkowycz, “Free convection about a vertical flat plate embedded in a porous medium with application to heat transfer from a dike,” Journal of Geophysical Research, vol. 82, no. 14, pp. 2040–2044, 1977. View at: Publisher Site  Google Scholar
 O. E. Potter, “Laminar boundary layers at the interface of cocurrent parallel streams,” The Quarterly Journal of Mechanics and Applied Mathematics, vol. 10, pp. 302–311, 1957. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 J. A. Ackroyd, “On the laminar compressible boundary layer with stationary origin on a moving flat wall,” Mathematical Proceedings of the Cambridge Philosophical Society, vol. 63, no. 3, pp. 871–888, 1967. View at: Publisher Site  Google Scholar
 J.H. He, “A simple perturbation approach to Blasius equation,” Applied Mathematics and Computation, vol. 140, no. 23, pp. 217–222, 2003. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 S. Abbasbandy, “A numerical solution of Blasius equation by Adomian's decomposition method and comparison with homotopy perturbation method,” Chaos, Solitons & Fractals, vol. 31, no. 1, pp. 257–260, 2007. View at: Publisher Site  Google Scholar
 L.T. Yu and C.K. Chen, “The solution of the blasius equation by the differential transformation method,” Mathematical and Computer Modelling, vol. 28, no. 1, pp. 101–111, 1998. View at: Publisher Site  Google Scholar
 J. H. He, “Approximate analytical solution of Blasius' equation,” Communications in Nonlinear Science and Numerical Simulation, vol. 4, no. 1, pp. 75–78, 1999. View at: Publisher Site  Google Scholar
 B. I. Yun, “Constructing uniform approximate analytical solutions for the Blasius problem,” Abstract and Applied Analysis, vol. 2014, Article ID 495734, 6 pages, 2014. View at: Publisher Site  Google Scholar  MathSciNet
 I. Ahmad and M. Bilal, “Numerical solution of blasius equation through neural networks algorithm,” American Journal of Computational Mathematics, vol. 4, no. 3, pp. 223–232, 2014. View at: Publisher Site  Google Scholar
 V. Marinca and N. Herişanu, “The optimal homotopy asymptotic method for solving Blasius equation,” Applied Mathematics and Computation, vol. 231, pp. 134–139, 2014. View at: Publisher Site  Google Scholar  MathSciNet
 A. Ebaid and N. AlArmani, “A new approach for a class of the blasius problem via a transformation and adomian's method,” Abstract and Applied Analysis, vol. 2013, Article ID 753049, 8 pages, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 O. Costin, T. E. Kim, and S. Tanveer, “A quasisolution approach to nonlinear problems—the case of the Blasius similarity solution,” Fluid Dynamics Research, vol. 46, no. 3, 2014. View at: Publisher Site  Google Scholar  MathSciNet
 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  MathSciNet
 L. Howarth, “On the solution of the laminar boundary layer equations,” Proceedings of the Royal Society A, vol. 164, no. 919, pp. 547–579, 1938. View at: Publisher Site  Google Scholar
 J. H. He, “Approximate analytical solution of Blasius' equation,” Communications in Nonlinear Science and Numerical Simulation, vol. 3, no. 4, pp. 260–263, 1998. View at: Publisher Site  Google Scholar
 G. A. Baker and P. R. GravesMorris, Essentials of Padé Approximants, Academic Press, New York, NY, USA, 1975.
 J. P. Boyd, “Padé approximant algorithm for solving nonlinear ordinary differential equation boundary value problems on an unbounded domain,” Computers in Physics, vol. 11, no. 3, article 299, 1997. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2015 M. Aghakhani 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.