Research Article | Open Access

Volume 2015 |Article ID 671527 | https://doi.org/10.1155/2015/671527

M. Aghakhani, M. Suhatril, M. Mohammadhassani, M. Daie, A. Toghroli, "A Simple Modification of Homotopy Perturbation Method for the Solution of Blasius Equation in Semi-Infinite 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 Semi-Infinite Domains

Revised28 Aug 2015
Accepted09 Sep 2015
Published08 Nov 2015

#### 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 Runge-Kutta 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) , Adomian Decomposition Method (ADM) , Variational Iteration Method (VIM) , Differential Transformation Method (DTM) , and Homotopy Perturbation Method (HPM)  can be mentioned. The homotopy perturbation method, first proposed by He , 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) , Generalized Homotopy Method (GHM) , and -expansion method .

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  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  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. . 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. . In their study, Sumudu transform was utilized for transformation of partial differential equations into a new form. Dong et al.  obtained the solution of strongly nonlinear mixed Volterra-Fredholm 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 semi-infinite 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 well-known classical Blasius first derived by Blasius  and dates back about a century, which describes the velocity profile of two-dimensional 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 , Laminar boundary layers at the interface of cocurrent parallel streams , or the flow near the leading edge of a very long, steadily operating conveyor belt . 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 . 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 .

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 , 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 .

#### 3. Application of HPM-Padé 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

 Order Error 0.33548 0.33267 0.33213 0.332052 0.3320549 Numerical  0.33205734
 Exact,Howarth  VIM,He  HPM-Padé(present) Relative error (%)(VIM) Relative error (%)(HPM-Padé) 0 0 0 0 0 0 1 0.16577 0.19319 0.16557 16.68 0.12 2 0.65003 0.6794 0.65001 4.51 0.003 3 1.39682 1.39106 1.39679 0.41 0.002 4 2.30576 2.24573 2.30572 2.6 0.001 5 3.28329 3.17748 3.28309 3.22 0.006 6 4.27964 4.14688 4.27767 3.10 0.04 7 5.27926 5.13359 5.26736 2.76 0.22
 Exact,Howarth  VIM,He  HPM-Padé(present) Relative error (%)(VIM) Relative error (%)(HPM-Padé) 0 0 0 0 0 0 1 0.32979 0.35064 0.32977 6.32 0.006 2 0.62977 0.61218 0.62976 2.79 0.001 3 0.84605 0.7964 0.84603 5.87 0.002 4 0.95552 0.90185 0.95551 5.62 0.001 5 0.99150 0.95523 0.99152 3.62 0.002 6 0.99868 0.98032 0.99883 1.84 0.015 7 0.99992 0.99158 0.99943 0.83 0.48
 Exact,Howarth  VIM,He  HPM-Padé(present) Relative error (%)(VIM) Relative error (%)(HPM-Padé) 0 0.33206 0.5436 0.33205 63.7 0.003 1 0.32301 0.27141 0.32300 15.97 0.003 2 0.26675 0.22748 0.26675 14.72 0 3 0.16136 0.14117 0.16135 12.52 0.006 4 0.06424 0.07469 0.06422 16.27 0.031 5 0.01591 0.036 0.01586 126.27 0.314 6 0.0024 0.01645 0.00110 — — 7 0.00022 0.00723 0.00060 — —

#### 4. Application of HPM-Padé 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 Runge-Kutta 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 68, the result of 8th order HPM-Padé 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.

 Order Error 5th order Runge-Kutta Fehlberg (present)
 Numerical(5th order Runge-Kutta Fehlberg) HPM-Padé(present) Relative error (%)(HPM-Padé) 0 0 0 0 1 0.786198 0.78620 0.00047 2 1.218546 1.21855 0.00053 3 1.432728 1.43273 0.0002 4 1.533086 1.53308 0.00022 5 1.578851 1.57884 0.00029 6 1.599437 1.59945 0.0011 7 1.612470 1.61280 0.002
 Numerical(5th order Runge-Kutta Fehlberg) HPM-Padé(present) Relative error (%)(HPM-Padé) 0 1 1 0 1 0.587153 0.587153 0.000021 2 0.301784 0.301783 0.000056 3 0.144016 0.144016 0.00012 4 0.066244 0.066243 0.0013 5 0.029956 0.029949 0.00025 6 0.013469 0.013434 0.025 7 0.006119 0.006005 1.8
 Numerical(5th order Runge-Kutta Fehlberg) HPM-Padé(present) Relative error (%)(HPM-Padé) 0 0.00016 1 0.00064 2 0.00015 3 0.00017 4 0.0043 5 0.067 6 0.59 7 3.7

#### 5. Conclusion

In the present paper, a simple modification of the homotopy perturbation method is proposed for the solution of the Blasius equation in semi-infinite domains. The equation in semi-infinite 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 RP004D-11AET. The authors would like to acknowledge the support.

1. S. Liao, “On the homotopy analysis method for nonlinear problems,” Applied Mathematics and Computation, vol. 147, no. 2, pp. 499–513, 2004.
2. G. Adomian, Solving Frontier Problems of Physics: The Decomposition Methoc [i.e Method], Kluwer Academic Publishers, Dordrecht, The Netherlands, 2013.
3. 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 Site | Google Scholar
4. M.-J. Jang, C.-L. Chen, and Y.-C. Liy, “On solving the initial-value problems using the differential transformation method,” Applied Mathematics and Computation, vol. 115, no. 2-3, pp. 145–160, 2000.
5. J.-H. He, “Homotopy perturbation method: a new nonlinear analytical technique,” Applied Mathematics and Computation, vol. 135, no. 1, pp. 73–79, 2003.
6. N. Herisanu, V. Marinca, and Gh. Madescu, “An analytical approach to non-linear dynamical model of a permanent magnet synchronous generator,” Wind Energy, vol. 18, pp. 1657–1670, 2015. View at: Publisher Site | Google Scholar
7. H. Vazquez-Leal, “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
8. M. N. Alam, M. A. Akbar, and S. T. Mohyud-Din, “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
9. 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
10. E. Yusufoğlu, “An improvement to homotopy perturbation method for solving system of linear equations,” Computers & Mathematics with Applications, vol. 58, no. 11-12, pp. 2231–2235, 2009. View at: Publisher Site | Google Scholar
11. 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
12. D. Kumar, J. Singh, and S. Kumar, “Numerical computation of fractional multi-dimensional 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
13. C. Dong, Z. Chen, and W. Jiang, “A modified homotopy perturbation method for solving the nonlinear mixed Volterra-Fredholm integral equation,” Journal of Computational and Applied Mathematics, vol. 239, pp. 359–366, 2013. View at: Publisher Site | Google Scholar | MathSciNet
14. 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
15. 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
16. O. E. Potter, “Laminar boundary layers at the interface of co-current parallel streams,” The Quarterly Journal of Mechanics and Applied Mathematics, vol. 10, pp. 302–311, 1957.
17. 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
18. J.-H. He, “A simple perturbation approach to Blasius equation,” Applied Mathematics and Computation, vol. 140, no. 2-3, pp. 217–222, 2003.
19. 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
20. 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
21. 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
22. 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
23. 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
24. 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
25. A. Ebaid and N. Al-Armani, “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
26. O. Costin, T. E. Kim, and S. Tanveer, “A quasi-solution 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
27. J.-H. He, “Homotopy perturbation technique,” Computer Methods in Applied Mechanics and Engineering, vol. 178, no. 3-4, pp. 257–262, 1999. View at: Publisher Site | Google Scholar | MathSciNet
28. 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
29. 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
30. G. A. Baker and P. R. Graves-Morris, Essentials of Padé Approximants, Academic Press, New York, NY, USA, 1975.
31. 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

#### More related articles

We are committed to sharing findings related to COVID-19 as quickly as possible. We will be providing unlimited waivers of publication charges for accepted research articles as well as case reports and case series related to COVID-19. Review articles are excluded from this waiver policy. Sign up here as a reviewer to help fast-track new submissions.