APPROXIMATE SOLUTIONS FOR FLOW WITH A STRETCHING BOUNDARY DUE TO PARTIAL SLIP U. Filobello-Ninoa , H. Vazquez-Leala,∗, A. Sarmiento-Reyesb , B. Benhammoudac , V.M. Jimenez-Fernandeza , D. Pereyra-Diaza , A. Perez-Sesmaa , J. Cervantes-Pereza , J. Huerta-Chuad , J. Sanchez-Oreaa , A.D. Contreras-Hernandeza a Electronic Instrumentation and Atmospheric Sciences School, Universidad Veracruzana, Cto. Gonzalo Aguirre Beltr´ an S/N, Xalapa, Veracruz, 91000, M´exico. b National Institute for Astrophysics, Optics and Electronics, Luis Enrique Erro No. 1, Sta. Maria Tonantzintla, 72840, Puebla, Mexico. c Abu Dhabi Mens College, Higher Colleges of Technology, P.O. Box 25035, Abu Dhabi, United Arab Emirates. d Civil Engineering School, Universidad Veracruzana, Venustiano Carranza S/N, Col. Revolucion, Poza Rica, Veracruz, 93390, Mexico. Abstract In this article, the homotopy perturbation method (HPM) is coupled with versions of Laplace-Pad´e and Pad´e 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 to provide an accurate solution and are extremely efficient. Keywords: Boundary layer, fluid mechanics, homotopy perturbation method, Pad´e method, Laplace Pad´e method, partial slip, extrusion processes, viscous. 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,2,3]. It is common to define the boundary layer as the region where the fluid velocity parallel to the surface is less than 99% 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,52,53] (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 ∗ Corresponding author: [email protected] 1 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 others 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 of slip boundary conditions [58], there are cases where partial slip between the fluid and the moving surface may occur. Some known cases include emulsions, as mustard and paints; polymer solutions and clay [58]. v u Y X U=bX Figure 1: Schematic showing a stretching boundary. Ji-Huan He [5,6] 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,6,20-30,33-38,39,41,42,44-49,52]. The fundamentals of HPM convergence can be found in [33,36,37]. 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 [7-9, 31], tanh method [10], exp-function [11,12], Adomians decomposition method [13-18,40], parameter expansion [19], homotopy analysis method [4,32,43], and perturbation method [50] among many others. To figure out how HPM method works, consider a general nonlinear equation in the form A(u) − f (r) = 0, r ∈ Ω, with the following boundary conditions 2 (1) B(u, ∂u/∂n) = 0, r ∈ Γ, (2) where A is a general differential operator, B is a boundary operator, f (r) a known analytical function, and Γ is the domain boundary for Ω. Also, A can be divided into two parts L and N, where L is linear and N nonlinear; from this last statement, (1) can be rewritten as L(u) + N(u) − f (r) = 0. (3) In a broad sense, a homotopy can be constructed in the following form[5,6] H(v, p) = (1 − p)[L(v) − L(u0 )] + p[L(v) + N(v) − f (r)] = 0, p ∈ [0, 1], r ∈ Ω, (4) H(v, p) = L(v) − L(u0 ) + p[L(u0 ) + N(v) − f (r)] = 0, p ∈ [0, 1], r ∈ Ω. (5) or where p is a homotopy parameter, whose values are within range of 0 and 1, u0 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 p v = v0 + v1 p1 + v2 p22 + · · · . (6) Substituting (6) into (5) and equating identical powers for p terms, it is possible to obtain the values for the sequence u0 , u1 , u2 , . . . When p → 1, it yields in the approximate solution for (1) in the form v = v0 + v1 + v2 + v3 + · · · . (7) Another way to build a homotopy, which is relevant for this paper, is by considering the following general equation L(v) + N(v) = 0, (8) where L(v) and N(v) are the linear and no linear operators, respectively. It is desired that solution for L(v) = 0 describes, accurately, the original nonlinear system. By the homotopy technique, a homotopy is constructed as follows [28] (1 − p)L(v) + p[L(v) + N(v)] = 0. (9) Again, it is assumed that solution for (9) can be written in the form (6); thus, taking the limit when p → 1, results in the approximate solution for (8). The variation of homotopic parameter within the range [0, 1] amounts to a deformation that begins from an initial equation with known solution until it becomes to the equation to be solved. From a practical point of view, taking the limit p → 1, is just setting p = 1. 3 2. Pad´ e Approximant Let u(t) be an analytical function with the Maclaurin’s expansion u (t) = ∞ X n=0 un tn , 0 ≤ t ≤ T . (10) Then the Pad´e approximant to u (t) of order [L, M] which we denote by [L/M]u (t) is defined by [55,56,57,59] p0 + p1 t + . . . + pL tL [L/M]u (t) = , (11) 1 + q1 t + . . . + qM tM where we considered q0 = 1, and the numerator and denominator have no common factors. The numerator and the denominator in (11) are constructed so that u (t) and [L/M]u (t) and their derivatives agree at t = 0 up to L + M. That is u(t) − [L/M]u (t) = O tL+M +1 . (12) From (12), we have u (t) M X n=0 n qn t − L X n=0 pn tn = O tL+M +1 . From (13), we get the following algebraic linear systems uL q1 + . . . + uL−M +1qM = −uL+1 uL+1 q1 + . . . + uL−M +2qM = −uL+2 .. . u L+M −1 q1 + . . . + uL qM = −uL+M , and p0 = u 0 p1 = u1 + u0 q1 .. . p = u + u q + ...+ u q . L L L−1 1 0 L (13) (14) (15) From (14), we calculate first all the coefficients qn , 1 ≤ n ≤ M. Then, we determine the coefficients pn , 0 ≤ n ≤ L from (15). Note that for a fixed value of L+ M + 1, 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´ e resummation method Several approximate methods provide power series solutions (polynomial). Nevertheless, sometimes, this type of solutions lacks of large domains of convergence. Therefore, LaplacePad´e resummation method [55] is used in literature to enlarge the domain of convergence of solutions or inclusive to find exact solutions. The Laplace-Pad´e method can be explained as follows: 4 1. First, Laplace transforma is applied to power series. 2. Next, s is substituted by 1/t in the resulting equation. 3. After that, we convert the transformed series into a meromorphic function by forming its Pad´e approximant of order [N/M]. N and M are arbitrarily chosen, but they should be of smaller values than the order of the power series. In this step, the Pad´e approximant extends the domain of the truncated series solution to obtain better accuracy and convergence. 4. Then, t is substituted by 1/s. 5. Finally, by using the inverse Laplace s 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 U is approximately proportional to the distance from the orifice X [51,58], so that U = bX. (16) where b is a proportionality constant. Let (u, v) be the fluid velocities for the (X, Y ) 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. ∂u (X, 0), (17) ∂Y where k 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. u(X, 0) − U = kν uuX + vuY + pX − ν(uXX + uY Y ) = 0, ρ (18) uvX + vvY + pY − ν(vXX + vY Y ) = 0, ρ (19) and continuity uX + vY = 0, (20) where ρ and p are density and pressure, respectively. To solve equations (18) –(20), we have to consider boundary conditions (16) and (17), besides 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). 5 With this end, we note that Y component of velocity is negative (v < 0) and by symmetry arguments it only depends on Y (see Figure 1). Therefore, it is possible to define a function y(x) ≥ 0, 0 ≤ x < ∞, such that [58]. √ v = − bνy(x), (21) where x = Y r a ν (22) From (21) and continuity equation (20) we obtain u = bXy ′ (x), (23) Clearly (20) is satisfied under transformations (21), (22) and (23), while (18), adopts the simpler form. y ′′′ (x) − (y ′(x))2 + y(x)y ′′(x) = 0. (24) To deduce the boundary conditions of (24), we see from Figure 1, that v(Y = 0) = 0 and lim u(X, Y ) = 0 and therefore Y →∞ y(0) = 0, y (∞) = 0, ′ (25) (26) Finally, substituting (16) and (23) into (17) we obtain √ y ′(0) = k bνy ′′(0) + 1. (27) 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 L = y ′′′ , (28) N = −(y ′ )2 + yy ′′, (29) and the nonlinear as we initiate the HPM method by constructing a homotopy based on (9), in the form 6 (1 − p)y ′′′ + p(y ′′′ − (y ′)2 + yy ′′) = 0. (30) Assuming that the solution has the form [5,6]. y = y0 + y1 p1 + y2 p22 + · · · . (31) Then, substituting (31) into (30), and equating terms having identical powers of p we obtain p0 p1 p2 p3 p4 p5 p6 p7 y0′′′ = 0, y1′′′ − (y0′ )2 + y0 y0′′ = 0, y1 y0′′ − 2y0′ y1′ + y2′′′ + y0 y1′′ = 0, 2 y2′′ + y1 y1′′ − y1′ + y2 y0′′ + y3′′′ − 2y0′ y2′ = 0, y1′′ + y0 y3′′ + y1 y2′′ − 2y0′ y3′ + y3 y0′′ + y4′′′ − 2y1′ y2′ = 0, 2 y3′′ − y2′ + y5′′′ − 2y1′ y3′ + y4 y0′′ + y2 y2′′ − 2y0′ y4′ + y0 y4′′ + y3 y1′′ = 0, y5 y0′′ + y0 y5′′ + y6′′′ − 2y0′ y5′ + y4 y1′′ + y1 y4′′ + y2 y3′′ − 2y1′ y4′ +y3y2′′ − 2y2′ y3′ = 0, 2 : y4′′ − y3′ + y4 y2′′ + y1 y5′′ + y0 y6′′ + y6 y0′′ − 2y2′ y4′ + y5 y1′′ + y3 y3′′ −2y0′ y6′ + y7′′′ − 2y1′ y5′ = 0. : : : : : : : In order to fulfil the√boundary conditions from (24) given by (25)–(27); we find that y0 (0) = 0, y0′ (0) = 1 + k bνa, y0′′ (0) = a, y1 (0) = 0, y1′ (0) = 0, y1′′(0) = 0, y2(0) = 0, y2′ (0) = 0, y2′′(0) = 0, y3(0) = 0, y3′ (0) = 0, y3′′(0) = 0, and so on. We have assumed the value for y ′′ (0) as some adequate constant a, which adopts the following values for the corresponding k as follows {k = 0, a = −1}, {k = 0.3, a = −0.701}, {k = 1, a = −0.430}, {k = 2, a = −0.284}, {k = 5, a = −0.145}, {k = 20, a = −0.0438} [53]. Thus, the results obtained from above equations are 7 y0 = y1 = y2 = y3 = y4 = 1 x 1+ x+k a , 2 1 1 1 2 1 2 2 1+ xk + x + k a + 2k + x a x3 , 6 4 20 4 6 ax 1 1 1 2 2 2 1 + − xk − x + k a + 2k − x a , 720 7 56 7 7 1 3 3 3 1 2 2 3 4 x 4 1+ xk + x k+ x k + x + k a4 2520 8 160 16 1760 1 2 3 2 3 3 3 3 x + 4k + x k + xk a + 160 8 8 1 2 3 1 2 2 + 6k + x + xk a + 4k + x a , 16 8 8 9 13 4 29 2 3 317 5 x 1+ xk + xk + x − 45360 40 440 640640 317 4 101 3 2 5 xk + x k + k a5 + 3520 45760 87 2 2 13 3 317 4 101 3 4 + x k + 5k + xk + x + x k a4 440 10 45760 1760 87 2 39 2 101 3 3 + x k + xk + x + 10k a3 440 20 3520 9 x 29 2 13 13 2 2 − 10k + x + xk a + 5k + x a . 45360 440 10 40 8 (32) (33) (34) (35) (36) y5 y6 1 589 3 3 37 2 4 419 4 2 19 5 = 1+ xk + xk + x k + xk 1247400 17472 208 26880 24 3431 3431 x6 + k 6 + kx5 a6 + 19009536 1118208 37 2 3 589 3 2 3431 5 95 4 + x k + 6k 5 + xk + x + xk 52 5824 1118208 24 419 4 x k a5 x11 + 13440 1 111 2 2 419 4 589 3 95 3 4 4 + x k + 15k + x + x k + xk a 1247400 104 26880 5824 12 95 37 589 3 3 + 20k 3 + xk 2 + x2 k + x a 12 52 17472 19 37 2 95 2 2 xk + x + 15k a + x + 6k a x11 . + 24 208 24 1 599 4 3 351 6 1693 3 4 3257 2 5 = 1+ − xk − xk − xk − xk 97297200 60928 224 26880 3360 124291 6 124291 7 7 35137 5 2 7 xk − x x− x a +k − 2437120 46305280 926105600 3257 2 4 1693 3 3 35137 5 1053 5 124291 6 + − xk − xk − x k− xk − x 672 6720 1218560 112 46305280 3257 2 3 1 1797 4 2 6 13 6 a + − x + xk x k + 21k 5 +7k − 60928 97297200 336 1693 3 2 35137 5 1797 4 5265 4 5 − xk − x − x k− xk a 4480 2437120 60928 224 599 4 1693 3 1755 3 4 13 3257 2 2 4 x k + 35k − x − x k− xk a x + − 336 60928 6720 56 1 5265 2 3257 2 1693 3 3 3 + 35k − xk − x k− x a 97297200 224 672 26880 1053 3257 2 4 13 2 + − xk + 21k − x a x . 112 3360 9 (37) (38) y7 23 1875011 395 7 = − 1 + k8 − x8 + xk 9081072000 808351272960 552 392303 270805 5 3 1875011 x7 k − x6 k 2 + xk − 35145707520 4792596480 68465664 1037917 4 4 5075 3 5 8393 2 6 + a8 x15 xk + xk + xk 51349248 1351296 225216 1037917 8393 2 5 2765 6 23 x4 k 3 + x k + 8k 7 + xk − 9081072000 808351272960 37536 552 1875011 270805 5 2 25375 3 4 xk − x7 + xk + 1351296 35145707520 22821888 392303 25375 3 3 1037917 4 2 23 − x6 k a7 x15 − xk + xk 2396298240 9081072000 675648 8558208 270805 5 2765 5 41965 2 4 392303 + x k+ xk + xk − x6 22821888 184 75072 4792596480 6 15 23 41965 2 3 25375 3 2 6 +28k a x − x k + 56k 5 + xk 9081072000 56304 675648 270805 5 13825 4 1037917 4 41965 2 2 5 15 + x + xk + x k+ xk ax 68465664 552 12837312 75072a 1037917 4 25375 3 13825 3 23 70k 4 + x + x k+ xk − 9081072000 51349248 1351296 552 5075x3 2765xk 56k 3 2765xk 2 8393x2 k + + + + a4 x15 + a 184a 37536a 1351296a 552a2 23 8393 6 2 − x + 28k a2 x15 9081072000 225216 23 395x − + 8k ax15 , 9081072000 552 (39) and so on. With the end to exemplify, we will consider the seventh order approximation by substituting solutions (32)–(39) into (31) and calculating the limit when p → 1. ! 7 X (40) y = lim y i pi . p→1 i=0 The above approximation is valid for all values of k ≥ 0, and corresponds to the results reported in [52] for the cases when k = 0 and k = 20. Clearly, if higher order approximations are considered, better accuracy is obtained but the resulting expressions could be too long. 5.2. HPM Laplce Pad´e scheme (LPHPM) Next, we study the case k = 0, it means that y ′ (0) = 1 (see (27)). This is an interesting case because the following exact solution was reported in [54]. 10 y(x) = 1 − exp(−x). (41) Thus, substituting k = 0 in (40) we obtain the following series solution y(x) = −0.5x2 + x − 7.878129877 × 10−11 x19 + 1.13972148927 × 10−11 x20 +1.57691357363 × 10−10 x18 + 2.07320134446 × 10−13 x21 −1.35120811977 × 10−13 x22 + 5.87481791204 × 10−15 x23 +0.16666667x3 − 0.041666666x4 + 0.00833333333x5 −0.00138888888x6 + 0.000198412698x7 − 0.0000248015873x8 −2.755731 × 10−7 x10 + 0.000002755x9 + 2.5052108 × 10−8 x11 −2.087675 × 10−9 x12 + 1.605904 × 10−10 x13 −1.147074559 × 10−11 x14 + 7.647163731 × 10−13 x15 −4.779477332 × 10−14 x16 − 5.073836906 × 10−11 x17 (42) As we will see (42) is accurate only for small values of x. To guarantee the validity of the approximate solution (42) for large values of this variable, the series solution is transformed using the Pad´e approximation and Laplace transform (see section 3). As first step, Laplace transform [55] is applied to (42). −1 1 9.583359 × 106 2.77283070001 × 107 1.009599 × 106 + − + + s3 s2 s20 s21 s19 8 8 7 1.51875891 × 10 1.51875891 × 10 1 1 1.0592181 × 10 − + + 4− 5 + 22 23 24 s s s s s 1 1 1 1 1 0.9999999998 1 0.9999999999 1 − 13 + + 6 − 7 + 8 − 9 − 11 + 10 + s s s s s s s12 s s14 0.9999999997 1 1 18047 − + 16 − 17 − 18 , 15 s s s s (43) then, s is substituted by 1/x in the equation to obtain −x3 + x2 − 9.583359 × 106 x20 + 2.77283070001 × 107 x21 +1.009599 × 106 x19 + 1.0592181 × 107 x22 − 1.51875891 × 108 x23 +x4 − x5 + x6 − x7 + x8 − x9 − x11 + x10 + 0.9999999998x12 − x13 +0.9999999999x14 − 0.9999999997x15 + x16 − x17 − 18047x18 +1.51875891 × 108 x24 . (44) Following Laplace-Pad´e scheme, Pad´e approximant [6/6] is applied to obtain x2 , 1+x here x is substituted by 1/s, the result is 11 (45) 1 , s(s + 1) (46) Finally, by means of the inverse Laplace transform applied to (46), we obtain the exact solution (41) for (24). 5.3. HPM Pad´e scheme (PHPM) Another way to recover lost information from the truncated series (40), is by means of applying the Pad´e approximant [7/7]. For the particular case when k = 0 (see (42)), we obtain the following rational approximation for (24) (using dummy variables a1 and b1) a1 = 1.15625115625 × 10−7 x7 + 0.000174825174825x5 +0.0320512820513x3 + x, b1 = 1 + 0.1153846153x2 + 0.001456876456x4 + 0.5x + 0.01602564102x3 +0.00008741258741x5 + 0.00000323750323x6 +5.781255781 × 10−8 x7 , therefore, the final result is y(x) = a1 b1 (47) 6. Discussion Figure 3 and Figure 5 show a comparison between Runge Kutta 4 (RK4) numerical solution for different values of k, and approximations given by Laplace Pad´e and Pad´e methods applied to (40). From Figure 4 and Figure 6 can be noticed that the relative error obtained by our approximations were low. As a matter of fact, the largest error for LPHPM is 0.006 when k = 5, and for PHPM is -0.0250 when k = 0. In a broad sense, if higher order approximations are considered, better the accuracy of LP-HPM and PHPM methods. In particular, is noteworthy that Laplace Pad´e method allows recovering the exact solution (41) when k = 0 from the truncated series (40). This contrasts with the Figure 2, which shows the comparison between RK4 and HPM approximations given by (40) for k = 0, 0.3, 1, 2, 5 and k = 20. It is evident that HPM series is accurate only for a restricted domain of values for the independent variable x. Also it is worth to compare the cumbersome approximation (42) for k = 0, with the handy expressions (41) and (47) obtained by LHPM and PHPM respectively. In [52] HPM was employed to solve (24) with a good approximation for a restricted domain of values for x and small values for k; while in [53], the same equation was solved using perturbation method (PM) for small values of k. Indeed 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 k; having low relative errors for x values within the range 0 ≤ x ≤ 30, as shown in Figure 4 and Figure 6. 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. 12 Besides, just as it was shown in one of our cases study, 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 k = 0. Finally, other reason why we are interested in obtaining analytical approximation solutions is that numerical algorithms could give some problems, such as numerical instabilities, oscillations, among others. This means that the numerical solutions may not correspond to the real solution of the original differential equation [61]. Just as it was seen with our HPM solutions, one disadvantage of approximate with polynomials is its tendency to oscillate, this gives rise 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 [57,60]. In order to improve the aforementioned, the use of rational functions was proposed, through the Pad´e approximation method. The Pad´e 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´e approximation is reduced to a Maclaurin Polynomial. Pad´e-approximate extrapolation technique consists in approximating a truncated series (such as those resulting of HPM) by a rational function, this 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´e approximant let obtain the limit of the function under study as x takes large values, if L = M (see Figures 3 and 5) [56]. In Fact, the rational functions whose denominator and numerator have the same values of L and M, or a degree almost identical, give rise to a results better than those obtained by methods based on polynomial functions [60]. Finally, the convergence is uniform in any compact region in the case of Pad´e approximation, whereas truncated Mclaurin series is valid only in a near neighbourhood of zero point. The above explains why HPM-Pad´e has larger intervals of convergence in comparison with HPM series [56]. Such as it was explained LPHPM employ Pad´e approximant to the results of apply Laplace transform to the truncated series of HPM, and the rational function obtained in this way acquire the benefits of Pad´e method above mentioned. 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´e, and Pad´e 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). Figure 4 and Figure 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 13 k k k k k k Figure 2: Fourth order Runge Kutta numerical solution for (24) (solid line) and HPM solution (40) (symbols). k k k k k k Figure 3: Fourth order Runge Kutta solution for (24) (symbols) and LPHPM (solid line). 14 k k k k k k Figure 4: Relative error for different cases of LPHPM. k k k k k k Figure 5: Fourth order Runge Kutta numerical solution for (24) (symbols) and PHPM (solid line). 15 k k k k k k Figure 6: Relative error for different cases of PHPM. problems, we conclude that LPHPM and PHPM are methods with high potential in the search for analytical approximate solutions for nonlinear problems. Acknowledgements We gratefully acknowledge the financial support from the National Council for Science and Technology of Mexico (CONACyT) through grant CB-2010-01 #157024. Disclosure Policy The authors declare that there is no conflict of interests regarding the publication of this paper. References [1] Hughes W F and Brighton J A, Din´amica De Los Fluidos, Mc Graw Hill. 1967. [2] Resnick R and Halliday D, F´ısica Vol. 1, John Wiley and Sons. 1977. [3] Landau L D and Lifschitz E M, Fluid Mechanics Vol. 6 2nd Edition, Pergamon Press. 1987. [4] V. Aliakbar, A. Alizadeh- Pahlavan and K. Sadaghy, The influence of thermal radiation on MHD flow of Maxwellian fluids above stretching sheets, Nonlinear Sciences and Numerical Simulation, Elsevier, 779-794. 2007. [5] He, J.H. A coupling method of a homotopy technique and a perturbation technique for nonlinear problems. Int. J. Non-Linear Mech., 351, 37-43, 1998, DOI: 10.1016/S00207462(98)00085-7. [6] He, J.H. Homotopy perturbation technique. Comput. Methods Applied Mech. Eng., 178, 257-262, 1999, DOI: 10.1016/S0045-7825(99)00018-3. [7] Assas, L.M.B. Approximate solutions for the generalized K-dV-Burgers equation by He’s Variational iteration method. Phys. Scr., 76, 161-164. 2007, DOI: 10.1088/00318949/76/2/008 16 [8] He, J.H. Variational approach for nonlinear oscillators. Chaos, Solitons and Fractals, 34, 1430-1439, 2007, DOI: 10.1016/j.chaos.2006.10.026. [9] Kazemnia, M., S.A. Zahedi, M. Vaezi and N. Tolou, Assessment of modified variational iteration method in BVPs high-order differential equations. Journal of Applied Sciences, 8, 4192-4197, 2008, DOI:10.3923/jas.2008.4192.4197. [10] Evans, D.J. and K.R. Raslan. The Tanh function method for solving some important nonlinear partial differential. Int. J. Computat. Math., 82, 897-905, 2005, DOI: 10.1080/00207160412331336026. [11] Xu, F. A generalized soliton solution of the Konopelchenko-Dubrovsky equation using Exp-function method. Zeitschrift Naturforschung - Section A Journal of Physical Sciences, 62(12), 685-688, 2007. [12] Mahmoudi, J., 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, 8, 358-363, 2008, DOI:10.3923/jas.2008.358.363. [13] Adomian, G. A review of decomposition method in applied mathematics. Mathematical Analysis and Applications. 135, 501-544, 1998. [14] Babolian, E. and J. Biazar. On the order of convergence of Adomian method. Applied Mathematics and Computation, 130(2), 383-387. 2002, DOI: 10.1016/S0096-3003(01)001035. [15] Kooch, A. 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, 7, 57-67. 2012, DOI:10.3923 / tasr.2012.57.67 [16] Kooch, A. 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, 11, 3421-3428, 2011, DOI:10.3923/jas.2011.3421.3428 [17] Vanani, S. K., S. Heidari and M. Avaji. A low-cost numerical algorithm for the solution of nonlinear delay boundary integral equations. Journal of Applied Sciences, 11, 3504-3509 2011, DOI:10.3923jas.2011.3504.3509. [18] Chowdhury, S. H. A comparison between the modified homotopy perturbation method and Adomian decomposition method for solving nonlinear heat transfer equations. Journal of Applied Sciences, 11, 1416-1420, 2011, DOI:10.3923/jas.2011.1416.1420. [19] Zhang, L.-N. and L. Xu. Determination of the limit cycle by He’s parameter expansion for oscillators in a potential. Zeitschrift f¨ ur Naturforschung - Section A Journal of Physical Sciences, 62(7-8), 396-398, 2007. [20] He, J.H. Homotopy perturbation method for solving boundary value problems. Physics Letters A, 350(1-2), 87-88, 2006. [21] He, J.H. Recent Development of the Homotopy Perturbation Method. Topological Methods in Nonlinear Analysis, 31(2), 205-209, 2008. [22] Belendez, A., C. Pascual, M.L. Alvarez, D.I. M´endez, M.S. Yebra and A. Hern´andez. High order analytical approximate solutions to the nonlinear pendulum by He’s homotopy method. Physica Scripta, 79(1), 1-24 2009, DOI: 10.1088/0031-8949/79/01/015009. [23] He, J.H. A coupling method of a homotopy and a perturbation technique for nonlinear problems. International Journal of Nonlinear Mechanics, 35(1), 37-43, 2000. 17 [24] El-Shaed, M. Application of He’s homotopy perturbation method to Volterra’s integro differential equation. International Journal of Nonlinear Sciences and Numerical Simulation, 6, 163-168, 2005. [25] He, J.H. Some Asymptotic Methods for Strongly Nonlinear Equations. International Journal of Modern Physics B, 20(10), 1141-1199, 2006 DOI: 10.1142/S0217979206033796. [26] Ganji, D.D, H. Babazadeh, F Noori, M.M. Pirouz, M Janipour. An Application of Homotopy Perturbation Method for Non linear Blasius Equation to Boundary Layer Flow Over a Flat Plate, ACADEMIC World Academic Union, ISNN 1749-3889(print), 1749-3897 (online). International Journal of Nonlinear Science, 7(4), 309-404, 2009. [27] Ganji, D.D., H. Mirgolbabaei , Me. Miansari and Mo. 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, 8, 1256-1261 2008, DOI:10.3923/jas.2008.1256.1261. [28] Fereidon, A., Y. Rostamiyan, M. Akbarzade and D.D. Ganji. Application of Hes homotopy perturbation method to nonlinear shock damper dynamics. Archive of Applied Mechanics, 80(6), 641-649, 2010, DOI: 10.1007/s00419-009-0334-x. [29] Sharma, P.R. and G. Methi. Applications of homotopy perturbation method to partial differential equations. Asian Journal of Mathematics & Statistics, 4, 140-150, 2011, DOI:10.3923/ajms.2011.140.150. [30] Aminikhah Hossein. Analytical Approximation to the Solution of Nonlinear Blasius Viscous Flow Equation by LTNHPM. International Scholarly Research Network ISRN Mathematical Analysis, Volume 2012, Article ID 957473, 10 pages, 2011, doi: 10.5402/2012/957473. [31] Noorzad, R., A. Tahmasebi Poor and M. Omidvar, Variational iteration method and homotopy-perturbation method for solving Burgers equation in fluid dynamics. Journal of Applied Sciences, 8, 369-373, 2008. [32] Patel, T., M.N. Mehta and V.H. Pradhan. The numerical solution of Burgers equation arising into the irradiation of tumour tissue in biological diffusing system by homotopy analysis method. Asian Journal of Applied Sciences, 5, 60-66, 2012, DOI:10.3923/ajaps.2012.60.66. [33] Vazquez-Leal H, U. Filobello-Nio, R. Castaeda-Sheissa, L. Hernandez Martinez and A. Sarmiento-Reyes. Modified HPMs inspired by homotopy continuation methods. Mathematical Problems in Engineering, Vol. 2012, 20 pages, Article ID 309123, 2012, DOI: 10.155/2012/309123. [34] Vazquez-Leal H., R. Castaneda-Sheissa, U. Filobello-Nino, A. Sarmiento-Reyes, and J. Sanchez-Orea. High accurate simple approximation of normal distribution related integrals. Mathematical Problems in Engineering, Vol. 2012, 22 pages Article ID 124029, 2012, DOI: 10.1155/2012/124029. [35] U. Filobello-Nino, H. Vazquez-Leal, R. Castaneda-Sheissa, A. Yildirim, L. Hernandez Martinez, D. Pereyra Diaz, A. Perez Sesma and C. Hoyos Reyes. An approximate solution of Blasius equation by using HPM method. Asian Journal of Mathematics and Statistics, Vol. 2012, 10 pages, 2012, DOI: 10.3923 /ajms.2012, ISSN 1994-5418. [36] Biazar, J. and H. Aminikhan. Study of convergence of homotopy perturbation method for systems of partial differential equations. Computers and Mathematics with Applications, 58(11-12), 2221-2230, 2009. [37] Biazar, J. and H. Ghazvini. Convergence of the homotopy perturbation method for partial differential equations. Nonlinear Analysis: Real World Applications, 10(5), 263318 2640, 2009. [38] U. Filobello-Nino, H. Vazquez-Leal, Y. Khan, A. Yildirim, D. Pereyra Diaz, A. Perez Sesma, L. Hernandez Martinez, J. Sanchez Orea, R. Castaneda-Sheissa, and F. Rabago Bernal. HPM Applied to Solve Nonlinear Circuits: A Study Case. Applied Mathematics Sciences, 6(87), 4331 - 4344, 2012. [39] Noor Muhammad Aslam, Syed Tauseef Mohyud- Din. Homotopy approach for perturbation method for solving Thomas- Fermi equation using Pade Approximants. International of nonlinear science, 8, 5 pages, 2009. [40] Y.C. Jiao, Y. Yamamoto, et al. An Aftertreatment Technique for Improving the Accuracy of Adomian,s Decomposition Method. An international Journal computers and mathematics with applications, 43, 783-798, 2001. [41] Xinlong Feng, Liquan Mei, and Guoliang He. An efficient algorithm for solving Troesch’s problem. Applied Mathematics and Computation, 189(1), 500-507, 2007. [42] S. H Mirmoradia, I. Hosseinpoura, S. Ghanbarpour, and A. Barari. Application of an approximate analytical method to nonlinear Troesch’s problem. Applied Mathematical Sciences, 3(32), 1579-1585, 2009. [43] Hany N. Hassana and Magdy 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, 34, 977-989, 2011. [44] Mohammad Madani, Mahdi Fathizadeh, Yasir Khan, and Ahmet Yildirim. On the coupling of the homotopy perturbation method and Laplace transformation. Mathematical and Computer Modelling, 53 (910), 1937-1945, 2011. [45] M Fathizadeh, M Madani, Yasir Khan, Naeem Faraz, Ahmet Yildirim and Serap Tutkun. An effective modification of the homotopy perturbation method for mhd viscous flow over a stretching sheet. Journal of King Saud University-Science, 25(2), 107-113, 2011. [46] Yasir Khan, Qingbiao Wu, Naeem Faraz, and Ahmet Yildirim. The effects of variable viscosity and thermal conductivity on a thin film flow over a shrinking/stretching sheet. Computers and Mathematics with Applications, 61(11), 3391-3399, 2011. [47] Naeem Faraz and Yasir Khan. Analytical solution of electrically conducted rotating flow of a second grade fluid over a shrinking surface. Ain Shams Engineering Journal, 2(34), 221 -226, 2011. [48] Yasir Khan, Hector Vazquez Leal and Q.Wu. An efficient iterated method for mathematical biology model, Neural Computing and Applications, 23(3-4), 677-682, 2012. [49] S. Thiagarajan, A. Meena, S. Anitha, and L. Rajendran. Analytical expression of the steadystate catalytic current of mediated bioelectrocatalysis and the application of He’s Homotopy perturbation method, J Math Chem Springer, 49(8), 1727-1740, 2011, DOI 10.1007/ s10910-011-9854-z. [50] U. Filobello-Nino, H. Vazquez-Leal, Y. Khan, A. Yildirim, V.M. Jimenez-Fernandez, A.L. Herrera May, R. Castaneda-Sheissa, and J.Cervantes-Perez. Using Perturbation method and Laplace-Pad´e approximation to solve nonlinear problems, Miskolc Mathematical Notes, 14(1), 89101, 2013. [51] J. Vleggar, Laminar boundary-layer behaviour on continuous accelerated surfaces. Chemical Engineering Science, 32, 1517-1525, 1977. [52] M.M. Rashidi and D.D. Gangi J. Homotopy perturbation method for solving flow in the extrusion processes. IJE Transactions A: Basics, 23, 267-272, 2010. 19 [53] C.Y.Wang. Flow due to a stretching boundary with partial slip-an exact solution of the Navier-Stokes equations. Chemical Engineering Science, 23, 267-272, 2002. [54] Crane, L.J Flow past a stretching plate. Zeitschrift f¨ ur Angewandte Mathematik und Physik, 21, 645-647, 1970. [55] Brahim Benhammouda, Hector Vazquez-Leal, L. Benhamouche, H. Zahed , Y. Abo Einaga, A. Sarmiento-Reyes, A. Marin-Hernandez and J. Huerta-Chua. Procedure for Exact Solutions of Stiff Ordinary Differential Equations Systems, British Journal of Mathematics & Computer Science, 4(23), 3252-3270, 2014. [56] Alev Kelleci and Ahmet Yildirim, An efficient numerical method for solving coupled Burgers’ equation by combining homotopy perturbation and pade techniques, Numerical Methods for Partial Differential Equations, 27(4), 982995, 2011. [57] Behrouz Raftari, Hojatollah Adibi, Ahmet Yildirim, Solution of the MHD Falkner-Skan flow by Adomian decomposition method and Pad´e approximants, International Journal of Numerical Methods for Heat & Fluid Flow, 22 (8), 1010 - 1020, 2012. [58] Filobello-Nino, U, Vazquez-Leal, H, Khan, Y, Perez-Sesma, A, Diaz-Sanchez, A, HerreraMay, A, Pereyra-Diaz, D, Castaneda-Sheissa, R, Jimenez-Fernandez, V.M., & CervantesPerez, J. A handy exact solution for flow due to a stretching boundary with partial slip. Revista mexicana de f´ısica E, 59(1), 51-55, 2013. [59] Hector Vazquez-Leal, Brahim Benhammouda, Uriel Filobello-Nino, Arturo SarmientoReyes, Victor Manuel Jimenez-Fernandez, Jose Luis Garcia-Gervacio, Jesus Huerta-Chua, Luis Javier Morales-Mendoza and Mario Gonzalez-Lee, Direct application of Pad´e approximant for solving nonlinear differential equations, SpringerPlus, 3(1), 563, 2014, doi:10.1186 / 2193-1801-3-563. [60] J. Faires, Richard Burden, Numerical Methods, 4th Edition, Brooks/Cole, Boston. 2012 [61] Hector Vazquez-Leal and Francisco Guerrero, Application of series method with Pad´e and Laplace-Pad´e resummation methods to solve a model for the evolution of smoking habit in Spain, Computational and Applied Mathematics, 33(1), 181-192, 2014. 20
© Copyright 2024