International Scholarly Research Notices
Volume 2014 (2014), Article ID 747098, 10 pages
http://dx.doi.org/10.1155/2014/747098

Research Article

Approximate Solutions for Flow with a Stretching Boundary due to Partial Slip

U. Filobello-Nino,¹ H. Vazquez-Leal,¹ A. Sarmiento-Reyes,² B. Benhammouda,³ V. M. Jimenez-Fernandez,¹ D. Pereyra-Diaz,¹ A. Perez-Sesma,¹ J. Cervantes-Perez,¹ J. Huerta-Chua,⁴ J. Sanchez-Orea,¹ and A. D. Contreras-Hernandez¹

¹Electronic Instrumentation and Atmospheric Sciences School, University of Veracruz, Circuito Gonzalo Aguirre Beltrán S/N, 91000 Xalapa, VER, Mexico
²National Institute for Astrophysics, Optics and Electronics, Luis Enrique Erro No. 1, Santa Maria Tonantzintla, 72840 Puebla, PUE, Mexico
³Abu Dhabi Men’s College, Higher Colleges of Technology, P.O. Box 25035, Abu Dhabi, UAE
⁴Civil Engineering School, University of Veracruz, Venustiano Carranza S/N, Colonia Revolucion, 93390 Poza Rica, VER, Mexico

Received 1 August 2014; Revised 26 October 2014; Accepted 26 October 2014; Published 24 November 2014

Academic Editor: Weimin Han

Copyright © 2014 U. Filobello-Nino 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.

The homotopy perturbation method (HPM) is coupled with versions of Laplace-Padé and Padé methods to provide an approximate solution to the nonlinear differential equation that describes the behaviour of a flow with a stretching flat boundary due to partial slip. Comparing results between approximate and numerical solutions, we concluded that our results are capable of providing an accurate solution and are extremely efficient.

1. Introduction

According to the classification of Prandtl, the fluid motion is divided into two regions. The first region is near the object where the effect of friction is important and is known as the boundary layer, while, for the second type, these effects can be neglected [1–3]. It is common to define the boundary layer as the region where the fluid velocity parallel to the surface is less than of the free stream velocity [1].

The boundary layer thickness increases from the edge along the surface on which the fluid moves. Even for the case of a laminar flow, the exact solution of equations describing the laminar boundary layer is very difficult to calculate and only few simple problems can straightforward be analysed [1, 3].

An interesting case is the one where a flow is induced into a viscoelastic fluid by a linearly stretched sheet [4–6] (see Figure 1). Extrusion of molten polymers through a slit die for the production of plastic sheets is an important process in polymer industry [4]. The process is normally complicated from the physical point of view, because it requires significant heat transfer between the sheet and a surrounding fluid that plays the role of a cooling medium. An important aspect of the flow is the extensibility of the sheet which can be employed to improve its mechanical properties along the sheet. To obtain better results is necessary to improve the cooling rate, whereby it is common to add some polymeric additives into water (which is one of the most employed fluids as cooling medium) in order to have a better control on the cooling rate. The flow due to a stretching boundary is also important in other engineering processes of interest, such as the glass fibre drawing and crystal growing among many others. Detail discussion about this topic can be found in [4]. This work assumes that the boundary conditions for the problem under study are adequately described by Navier’s condition, which states that the amount of relative slip is proportional to local shear stress. Unlike what happens with fluids like water, mercury, and glycerine, which do not require slip boundary conditions [7], there are cases where partial slip between the fluid and the moving surface may occur. Some known cases include emulsions, as mustard and paints and polymer solutions and clay [7].

Figure 1: Schematic showing a stretching boundary.

He [8, 9] proposed the standard HPM; it was introduced as a powerful tool to approach several kinds of nonlinear problems. The HPM can be considered as a combination between the classical perturbation technique and the homotopy (whose origin is in the topology), but not restricted to the limitations found in traditional perturbation methods. For instance, HPM method does need neither small parameter nor linearisation, just few iterations to obtain accurate results [5, 8–35]. The fundamentals of HPM convergence can be found in [21, 24, 25].

There are other modern alternatives to find approximate solutions to the differential equations that describe some nonlinear problems such as those based on variational approaches [36–39], tanh method [40], exp-function [41, 42], Adomian decomposition method [43–49], parameter expansion [50], homotopy analysis method [4, 51, 52], and perturbation method [53] among many others.

To figure out how HPM method works, consider a general nonlinear equation in the form

with the following boundary conditions:

where is a general differential operator, is a boundary operator, is a known analytical function, and is the domain boundary for .

Also, can be divided into two parts and , where is linear and nonlinear; from this last statement, (1) can be rewritten as

In a broad sense, a homotopy can be constructed in the following form [8, 9]:

or

where is a homotopy parameter, whose values are within range of and , and is the first approximation to the solution of (3) that satisfies the boundary conditions. Assuming that solution for (4) or (5) can be written as a power series of

Substituting (6) into (5) and equating identical powers for terms, it is possible to obtain the values for the sequence .

When , it yields in the approximate solution for (1) in the form

Another way to build a homotopy, which is relevant for this paper, is by considering the following general equation:

where and are the linear and nonlinear operators, respectively. It is desired that solution for describes, accurately, the original nonlinear system.

By the homotopy technique, a homotopy is constructed as follows [18]:

Again, it is assumed that solution for (9) can be written in the form (6); thus, taking the limit when results in the approximate solution for (8).

The variation of homotopic parameter within the range amounts to a deformation that begins from an initial equation with known solution until it becomes the equation to be solved. From a practical point of view, taking the limit is just setting .

2. Padé Approximant

Let be an analytical function with the Maclaurin’s expansion

Then the Padé approximant to of order which we denote by is defined by [54–57]

where we considered , and the numerator and denominator have no common factors.

The numerator and the denominator in (11) are constructed so that and and their derivatives agree at up to . That is,

From (12), we have

From (13), we get the following algebraic linear systems:

From (14), we calculate first all the coefficients , . Then, we determine the coefficients , from (15).

Note that for a fixed value of , the error (12) is the smallest when the numerator and denominator of (11) have the same degree or when the numerator has one degree higher than the denominator.

3. Laplace-Padé Resummation Method

Several approximate methods provide power series solutions (polynomial). Nevertheless, sometimes, this type of solutions lacks large domains of convergence. Therefore, Laplace-Padé resummation method [54] is used in literature to enlarge the domain of convergence of solutions or inclusive to find exact solutions.

The Laplace-Padé method can be explained as follows.

(1)First, Laplace transform is applied to power series.(2)Next, is substituted by in the resulting equation.(3)After that, we convert the transformed series into a meromorphic function by forming its Padé approximant of order . and are arbitrarily chosen, but they should be of smaller values than the order of the power series. In this step, the Padé approximant extends the domain of the truncated series solution to obtain better accuracy and convergence.(4)Then, is substituted by .(5)Finally, by using the inverse Laplace transform, we obtain the exact or approximate solution.

4. Formulation

Consider a two-dimensional stretching boundary (see Figure 1). Experiments show that the velocity of the boundary is approximately proportional to the distance from the orifice [7, 58], so that

where is a proportionality constant.

Let be the fluid velocities for the directions, respectively. In this case the boundary condition is adequately described by Navier’s condition, which states that the amount of relative slip is proportional to local shear stress:

where is a constant of proportionality and is the kinematic viscosity of the bulk fluid. The relevant expressions for this case are the Navier-Stokes equations:

and continuity

where and are density and pressure, respectively.

To solve equations (18)–(20), we have to consider boundary conditions (16) and (17), besides the fact that there is no lateral velocity or pressure gradient away from the stretching surface.

Next, we will show that it is possible to get an ordinary nonlinear differential equation from (18).

With this end, we note that component of velocity is negative () and by symmetry arguments it only depends on (see Figure 1). Therefore, it is possible to define a function , , such that [7]:

where

From (21) and continuity equation (20) we obtain

Clearly (20) is satisfied under transformations (21), (22), and (23), while (18) adopts the simpler form;

To deduce the boundary conditions of (24), we see from Figure 1 that and and therefore

Finally, substituting (16) and (23) into (17) we obtain

5. Approximate Solution for a Two-Dimensional Viscous Flow Equation

Next, we present some solution methods to the study problem.

5.1. HPM Method

In this section, HPM is used to find approximate solutions for (24). Identifying the linear part as

and the nonlinear as

we initiate the HPM method by constructing a homotopy based on (9), in the form

Assuming that the solution has the form [8, 9]:

Then, substituting (30) into (29) and equating terms having identical powers of we obtain

In order to fulfil the boundary conditions from (24) given by (25)-(26), we find that , , , , , , , , , , , , and so on. We have assumed the value for as some adequate constant , which adopts the following values for the corresponding as follows: , , , , , and [6]. Thus, the results obtained from above equations are

and so on.

To exemplify, we will consider the seventh order approximation by substituting solutions (32) into (30) and calculating the limit when :

The above approximation is valid for all values of and corresponds to the results reported in [5] for the cases when and . Clearly, if higher order approximations are considered, better accuracy is obtained but the resulting expressions could be too long.

5.2. HPM Laplace-Padé Scheme (LPHPM)

Next, we study the case ; it means that (see (26)). This is an interesting case because the following exact solution was reported in [59]:

Thus, substituting in (33) we obtain the following series solution:

As we will see (35) is accurate only for small values of . To guarantee the validity of the approximate solution (35) for large values of this variable, the series solution is transformed using the Padé approximation and Laplace transform (see Section 3). As first step, Laplace transform [54] is applied to (35):

and then, is substituted by in the equation to obtain

Following Laplace-Padé scheme, Padé approximant is applied to obtain

here is substituted by ; the result is

Finally, by means of the inverse Laplace transform applied to (39), we obtain the exact solution (34) for (24).

5.3. HPM Padé Scheme (PHPM)

Another way to recover lost information from the truncated series (33) is by means of applying the Padé approximant . For the particular case when (see (35)), we obtain the following rational approximation for (24) (using dummy variables and )

and therefore, the final result is

6. Discussion

Figures 3 and 5 show a comparison between Runge Kutta 4 (RK4) numerical solution for different values of and approximations given by Laplace-Padé and Padé methods applied to (33). From Figures 4 and 6 it can be noticed that the relative error obtained by our approximations was low. As a matter of fact, the largest error for LPHPM is −0.006 when and for PHPM is −0.0250 when . In a broad sense, if higher order approximations are considered, higher accuracy is obtained from LP-HPM and PHPM methods. In particular, it is noteworthy that Laplace-Padé method allows recovering the exact solution (34) when from the truncated series (33). This contrasts with Figure 2, which shows the comparison between RK4 and HPM approximations given by (33) for and . It is evident that HPM series is accurate only for a restricted domain of values for the independent variable . Also it is worth comparing the cumbersome approximation (35) for , with the handy expressions (34) and (41) obtained by LHPM and PHPM, respectively. In [5] HPM was employed to solve (24) with a good approximation for a restricted domain of values for and small values for , while in [6] the same equation was solved using perturbation method (PM) for small values of . Indeed it is known that PM provides, in general, better results for small values of the perturbation parameter. Unlike the LPHPM and PHPM schemes were employed to obtain accurate solutions for different values of , having low relative errors for values within the range , as shown in Figures 4 and 6.

Figure 2: Fourth order Runge Kutta numerical solution for (24) (solid line) and HPM solution (33) (symbols).

Figure 3: Fourth order Runge Kutta solution for (24) (symbols) and LPHPM (solid line).

Figure 4: Relative error for different cases of LPHPM.

Figure 5: Fourth order Runge Kutta numerical solution for (24) (symbols) and PHPM (solid line).

Figure 6: Relative error for different cases of PHPM.

The importance of providing analytical solutions, although approximate, with good accuracy is that numerical solutions only provide a qualitative idea of the problem to be solved. Besides, just as it was shown in one of our case studies, a solution obtained by a numerical method like Runge Kutta could hide the case of an exact solution. This possibility was successfully explored by LPHPM for case study . Finally, another reason why we are interested in obtaining analytical approximation solutions is that numerical algorithms could give some problems, such as numerical instabilities and oscillations, among others. This means that the numerical solutions may not correspond to the real solution of the original differential equation [60].

Just as it was seen with our HPM solutions, one disadvantage of approximating them with polynomials is its tendency to oscillate; this gives rise to the fact that the obtained solutions diverge, especially for the case of problems defined on open intervals. This can be attributed to the possibility that the radius of convergence may not be sufficiently large to contain the boundaries of the domain of study [56, 61].

In order to improve the aforementioned, the use of rational functions was proposed, through the Padé approximation method.

The Padé approximation is an extension of the Taylor polynomial approximation but for rational functions. When the denominator of the approximation is a zero degree polynomial function, the Padé approximation is reduced to a Maclaurin polynomial.

Padé-approximate extrapolation technique consists in approximating a truncated series (such as those resulting of HPM) by a rational function; the latter extends the range of validity of the initial polynomial. The above results are particularly relevant if the truncated series represent the solution of a differential equation.

We observe that even though the series has a finite region of convergence, Padé approximant allows obtaining the limit of the function under study as takes large values, if (see Figures 3 and 5) [55]. In fact, the rational functions whose denominator and numerator have the same values of and , or a degree almost identical, give rise to results better than those obtained by methods based on polynomial functions [61].

Finally, the convergence is uniform in any compact region in the case of Padé approximation, whereas truncated Maclaurin series is valid only in a near neighbourhood of zero point. The above explains why HPM-Padé has larger intervals of convergence in comparison with HPM series [55].

As it was already mentioned, LPHPM applies Padé approximant to the resulting expression, derived of applying Laplace transform to the HPM truncated series, and the rational function obtained in this way acquires the benefits of the abovementioned Padé method.

Nevertheless, one possible advantage of LPHPM is that the application of inverse Laplace transform at the last step of the method could result in the exact solution of the problem.

7. Conclusions

This work showed that some nonlinear problems may be adequately approximated using the coupling of the HPM method with Laplace-Padé and Padé methods to deal with HPM truncated power series. For instance, the flow induced by a stretching sheet is adequately described by our approximations given by LPHPM and PHPM (see Figure 3 through Figure 6). Figures 4 and 6 illustrate the relative error and show that the proposed solutions are highly accurate. Since that this procedure is, in principle, applicable to other similar problems, we conclude that LPHPM and PHPM are methods with high potential in the search for analytical approximate solutions for nonlinear problems.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

The authors gratefully acknowledge the financial support from the National Council for Science and Technology of Mexico (CONACyT) through Grant CB-2010-01 no. 157024.

References

W. F. Hughes and J. A. Brighton, Dinámica De Los Fluidos, McGraw-Hill, New York, NY, USA, 1967.

R. Resnick and D. Halliday, Física, vol. 1, John Wiley & Sons, 1977.

L. D. Landau and E. M. Lifschitz, Fluid Mechanics, vol. 6, Pergamon Press, 2nd edition, 1987.

V. Aliakbar, A. Alizadeh-Pahlavan, and K. Sadeghy, “The influence of thermal radiation on MHD flow of Maxwellian fluids above stretching sheets,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 3, pp. 779–794, 2009.

M. M. Rashidi and D. D. Gangi, “Homotopy perturbation method for solving flow in the extrusion processes,” IJE Transactions A: Basics, vol. 23, pp. 267–272, 2010.

C. Y. Wang, “Flow due to a stretching boundary with partial slip—an exact solution of the Navier-Stokes equations,” Chemical Engineering Science, vol. 57, no. 17, pp. 3745–3747, 2002.

U. Filobello-Nino, H. Vazquez-Leal, Y. Khan et al., “A handy exact solution for flow due to a stretching boundary with partial slip,” Revista Mexicana de Física E, vol. 59, no. 1, pp. 51–55, 2013.

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.

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

J. H. He, “Homotopy perturbation method for solving boundary value problems,” Physics Letters A, vol. 350, no. 1-2, pp. 87–88, 2006.

J.-H. He, “Recent development of the homotopy perturbation method,” Topological Methods in Nonlinear Analysis, vol. 31, no. 2, pp. 205–209, 2008.

A. Belendez, C. Pascual, M. L. Alvarez, D. I. Méndez, M. S. Yebra, and A. Hernández, “Higher order analytical approximate solutions to the nonlinear pendulum by He's homotopy method,” Physica Scripta, vol. 79, no. 1, Article ID 015009, 24 pages, 2009.

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.

M. El-Shahed, “Application of He's homotopy perturbation method to Volterra's integro-differential equation,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 6, no. 2, pp. 163–168, 2005.

J.-H. He, “Some asymptotic methods for strongly nonlinear equations,” International Journal of Modern Physics B, vol. 20, no. 10, pp. 1141–1199, 2006.

D. D. Ganji, H. Babazadeh, F. Noori, M. M. Pirouz, and M. Janipour, “An application of homotopy perturbation method for non linear blasius equation to boundary layer flow over a flat plate,” International Journal of Nonlinear Science, vol. 7, pp. 309–404, 2009.

D. D. Ganji, H. Mirgolbabaei, M. Miansari, and M. Miansari, “Application of homotopy perturbation method to solve linear and non-linear systems of ordinary differential equations and differential equation of order three,” Journal of Applied Sciences, vol. 8, no. 7, pp. 1256–1261, 2008.

A. Fereidoon, Y. Rostamiyan, M. Akbarzade, and D. D. Ganji, “Application of He's homotopy perturbation method to nonlinear shock damper dynamics,” Archive of Applied Mechanics, vol. 80, no. 6, pp. 641–649, 2010.

P. R. Sharma and G. Methi, “Applications of homotopy perturbation method to partial differential equations,” Asian Journal of Mathematics & Statistics, vol. 4, no. 3, pp. 140–150, 2011.

H. Aminikhah, “Analytical approximation to the solution of nonlinear Blasius' viscous flow equation by LTNHPM,” ISRN Mathematical Analysis, vol. 2012, Article ID 957473, 10 pages, 2012.

H. Vázquez-Leal, U. Filobello-Niño, R. Castañeda-Sheissa, L. Hernández-Martínez, and A. Sarmiento-Reyes, “Modified HPMs inspired by homotopy continuation methods,” Mathematical Problems in Engineering, vol. 2012, Article ID 309123, 19 pages, 2012.

H. Vazquez-Leal, R. Castaneda-Sheissa, U. Filobello-Nino, A. Sarmiento-Reyes, and J. Sanchez Orea, “High accurate simple approximation of normal distribution integral,” Mathematical Problems in Engineering, vol. 2012, Article ID 124029, 22 pages, 2012.

U. Filobello-Nino, H. Vazquez-Leal, R. Castaneda-Sheissa et al., “An approximate solution of Blasius equation by using HPM method,” Asian Journal of Mathematics and Statistics, vol. 2012, 10 pages, 2012.

J. Biazar and H. Aminikhah, “Study of convergence of homotopy perturbation method for systems of partial differential equations,” Computers & Mathematics with Applications, vol. 58, no. 11-12, pp. 2221–2230, 2009.

J. Biazar and H. Ghazvini, “Convergence of the homotopy perturbation method for partial differential equations,” Nonlinear Analysis: Real World Applications, vol. 10, no. 5, pp. 2633–2640, 2009.

U. Filobello-Nino, H. Vazquez-Leal, Y. Khan et al., “HPM applied to solve nonlinear circuits: a study case,” Applied Mathematical Sciences, vol. 6, no. 87, pp. 4331–4344, 2012.

M. A. Noor and S. T. Mohyud-Din, “Homotopy perturbation method for solving Thomas-Fermi equation using Pade approximants,” International Journal of Nonlinear Science, vol. 8, no. 1, pp. 27–31, 2009.

X. Feng, L. Mei, and G. He, “An efficient algorithm for solving Troesch's problem,” Applied Mathematics and Computation, vol. 189, no. 1, pp. 500–507, 2007.

S. H. Mirmoradi, I. Hosseinpour, S. Ghanbarpour, and A. Barari, “Application of an approximate analytical method to nonlinear Troesch's problem,” Applied Mathematical Sciences, vol. 3, no. 32, pp. 1579–1585, 2009.

M. Madani, M. Fathizadeh, Y. Khan, and A. Yildirim, “On the coupling of the homotopy perturbation method and Laplace transformation,” Mathematical and Computer Modelling, vol. 53, no. 9-10, pp. 1937–1945, 2011.

M. Fathizadeh, M. Madani, Y. Khan, N. Faraz, A. Yildirim, and S. Tutkun, “An effective modification of the homotopy perturbation method for MHD viscous flow over a stretching sheet,” Journal of King Saud University—Science, vol. 25, no. 2, pp. 107–113, 2013.

Y. Khan, Q. Wu, N. Faraz, and A. Yildirim, “The effects of variable viscosity and thermal conductivity on a thin film flow over a shrinking/stretching sheet,” Computers & Mathematics with Applications, vol. 61, no. 11, pp. 3391–3399, 2011.

N. Faraz and Y. Khan, “Analytical solution of electrically conducted rotating flow of a second grade fluid over a shrinking surface,” Ain Shams Engineering Journal, vol. 2, no. 3-4, pp. 221–226, 2011.

Y. Khan, H. Vazquez-Leal, and Q. Wu, “An efficient iterated method for mathematical biology model,” Neural Computing and Applications, vol. 23, no. 3-4, pp. 677–682, 2013.

S. Thiagarajan, A. Meena, S. Anitha, and L. Rajendran, “Analytical expression of the steady-state catalytic current of mediated bioelectrocatalysis and the application of He's homotopy perturbation method,” Journal of Mathematical Chemistry, vol. 49, no. 8, pp. 1727–1740, 2011.

L. M. B. Assas, “Approximate solutions for the generalized KdV–Burgers' equation by He's variational iteration method,” Physica Scripta, vol. 76, no. 2, pp. 161–164, 2007.

J.-H. He, “Variational approach for nonlinear oscillators,” Chaos, Solitons & Fractals, vol. 34, no. 5, pp. 1430–1439, 2007.

M. Kazemnia, S. A. Zahedi, M. Vaezi, and N. Tolou, “Assessment of modified variational iteration method in BVPs high-order differential equations,” Journal of Applied Sciences, vol. 8, no. 22, pp. 4192–4197, 2008.

R. Noorzad, A. T. Poor, and M. Omidvar, “Variational iteration method and homotopy-perturbation method for solving Burgers equation in fluid dynamics,” Journal of Applied Sciences, vol. 8, no. 2, pp. 369–373, 2008.

D. J. Evans and K. R. Raslan, “The tanh function method for solving some important non-linear partial differential equations,” International Journal of Computer Mathematics, vol. 82, no. 7, pp. 897–905, 2005.

F. Xu, “A generalized soliton solution of the Konopelchenko-Dubrovsky equation using He's exp-function method,” Zeitschrift fur Naturforschung A Journal of Physical Sciences, vol. 62, no. 12, pp. 685–688, 2007.

J. Mahmoudi, N. Tolou, I. Khatami, A. Barari, and D. D. Ganji, “Explicit solution of nonlinear ZK-BBM wave equation using Exp-function method,” Journal of Applied Sciences, vol. 8, no. 2, pp. 358–363, 2008.

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.

E. Babolian and J. Biazar, “On the order of convergence of Adomian method,” Applied Mathematics and Computation, vol. 130, no. 2-3, pp. 383–387, 2002.

A. Kooch and M. Abadyan, “Efficiency of modified Adomian decomposition for simulating the instability of nano-electromechanical switches: comparison with the conventional decomposition method,” Trends in Applied Sciences Research, vol. 7, no. 1, pp. 57–67, 2012.

A. Koochi and M. Abadyan, “Evaluating the ability of modified adomian decomposition method to simulate the instability of freestanding carbon nanotube: comparison with conventional decomposition method,” Journal of Applied Sciences, vol. 11, no. 19, pp. 3421–3428, 2011.

S. K. Vanani, S. Heidari, and M. Avaji, “A low-cost numerical algorithm for the solution of nonlinear delay boundary integral equations,” Journal of Applied Sciences, vol. 11, no. 20, pp. 3504–3509, 2011.

S. H. Chowdhury, “A comparison between the modified homotopy perturbation method and adomian decomposition method for solving nonlinear heat transfer equations,” Journal of Applied Sciences, vol. 11, no. 7, pp. 1416–1420, 2011.

Y. C. Jiao, Y. Yamamoto, C. Dang, and Y. Hao, “An aftertreatment technique for improving the accuracy of Adomian's decomposition method,” Computers & Mathematics with Applications, vol. 43, no. 6-7, pp. 783–798, 2002.

L.-N. Zhang and L. Xu, “Determination of the limit cycle by He's parameter expansion for oscillators in a potential,” Zeitschrift für Naturforschung—Section A Journal of Physical Sciences, vol. 62, no. 7, pp. 396–398, 2007.

T. Patel, M. N. Mehta, and V. H. Pradhan, “The numerical solution of Burger's equation arising into the irradiation of tumour tissue in biological diffusing system by Homotopy Analysis Method,” Asian Journal of Applied Sciences, vol. 5, no. 1, pp. 60–66, 2012.

H. N. Hassan and M. A. El-Tawil, “An efficient analytic approach for solving two-point nonlinear boundary value problems by homotopy analysis method,” Mathematical Methods in the Applied Sciences, vol. 34, no. 8, pp. 977–989, 2011.

U. Filobello-Nino, H. Vazquez-Leal, Y. Khan et al., “Using perturbation methods and Laplace-Padé approximation to solve nonlinear problems,” Miskolc Mathematical Notesc, vol. 14, no. 1, pp. 89–101, 2013.

B. Benhammouda, H. Vazquez-Leal, L. Benhamouche et al., “Procedure for exact solutions of stiff ordinary differential equations systems,” The British Journal of Mathematics & Computer Science, vol. 4, no. 23, pp. 3252–3270, 2014.

A. Kelleci and A. Yıldırım, “An efficient numerical method for solving coupled Burgers' equation by combining homotopy perturbation and Pade techniques,” Numerical Methods for Partial Differential Equations, vol. 27, no. 4, pp. 982–995, 2011.

B. Raftari, H. Adibi, and A. Yildirim, “Solution of the MHD Falkner-Skan flow by Adomian decomposition method and Padé approximants,” International Journal of Numerical Methods for Heat and Fluid Flow, vol. 22, no. 8, pp. 1010–1020, 2012.

H. Vazquez-Leal, B. Benhammouda, U. Filobello-Nino, et al., “Direct application of Padé approximant for solving nonlinear differential equations,” SpringerPlus, vol. 3, article 563, 2014.

J. Vleggaar, “Laminar boundary-layer behaviour on continuous, accelerating surfaces,” Chemical Engineering Science, vol. 32, no. 12, pp. 1517–1525, 1977.

L. J. Crane, “Flow past a stretching plate,” Zeitschrift für Angewandte Mathematik und Physik, vol. 21, no. 4, pp. 645–647, 1970.

H. Vazquez-Leal and F. Guerrero, “Application of series method with Padé and Laplace-Padé resummation methods to solve a model for the evolution of smoking habit in Spain,” Computational & Applied Mathematics, vol. 33, no. 1, pp. 181–192, 2014.

J. Faires and R. Burden, Numerical Methods, Brooks/Cole, Boston, Mass, USA, 4th edition, 2012.