Mathematical Problems in Engineering

Volume 2015 (2015), Article ID 671527, 7 pages

http://dx.doi.org/10.1155/2015/671527

## A Simple Modification of Homotopy Perturbation Method for the Solution of Blasius Equation in Semi-Infinite Domains

^{1}Young Researchers and Elite Club, Islamic Azad University, Ilkhchi Branch, Ilkhchi, Iran^{2}Department of Civil Engineering, University of Malaya, 50603 Kuala Lumpur, Malaysia^{3}Department of Civil Engineering, University of Tabriz, Tabriz, Iran

Received 8 July 2015; Revised 28 August 2015; Accepted 9 September 2015

Academic Editor: Gerhard-Wilhelm Weber

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.

#### 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) [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 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 [14] 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 [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 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

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 HPM-Padé approximation are compared with Howarth’s numerical solution [28]. Furthermore, as it can be seen from Tables 2–4, purposed modification of HPM-Padé method is more accurate than the variational iteration method by He [29] and is valid for a wider range of the solution domain.