Mass Transpiration in Nonlinear MHD Flow Due to Porous Stretching Sheet.

Motivated from numerous practical applications, the present theoretical and numerical work investigates the nonlinear magnetohydrodynamic (MHD) laminar boundary layer flow of an incompressible, viscous fluid over a porous stretching sheet in the presence of suction/injection (mass transpiration). The flow characteristics are obtained by solving the underlying highly nonlinear ordinary differential equation using homotopy analysis method. The effect of parameters corresponding to suction/injection (mass transpiration), applied magnetic field, and porous stretching sheet parameters on the nonlinear flow is investigated. The asymptotic limits of the parameters regarding the flow characteristics are obtained mathematically, which compare very well with those obtained using the homotopy analysis technique. A detailed numerical study of the laminar boundary layer flow in the vicinity of the porous stretching sheet in MHD and offers a particular choice of the parametric values to be taken in order to practically model a particular type of the event among suction and injection at the sheet surface.

Introduction

The steady, laminar MHD boundary layer flows driven by moving boundaries are among the classical problems of theoretical fluid mechanics (see Schlichting[1]). The usual stretching sheet problems arise in polymer extrusion processes that involve the cooling of continuous strips extruded from a dye by drawing them horizontally through a stagnant cooling fluid[2].

The phenomena of momentum transfer in steady, laminar boundary layer flows have received much attention due to their wide applications which include drawing of plastics and elastic sheets, metal and polymer expulsion forms, paper creation, and cooling of metallic sheets etc.[2,3].

Viscous fluid flow past a linear stretching sheet is a classical problem of laminar boundary layer flow. Blasius[4] first discovered the boundary layer flow on a flat plate using similarity transformations. Sakiadis[5,6] investigated the steady laminar boundary layer flow on a moving plate in a quiescent liquid and obtained both closed form as well as approximate solutions. Crane[7] considered the flow due to stretching of plastic sheets in the polymer industry and obtained an analytical solution of the laminar boundary layer equations. Recently, Al-Housseiny and Stone[8] have investigated the laminar boundary layer flow due to motion of stretching sheets by taking into account both the fluid motion as well as the motion of the sheet.

The control of the boundary layer flow due to stretching sheet can be enhanced by introducing magnetohydrodynamic (MHD) effects. This can be done by taking an electrically conducting fluid above the sheet and applying a magnetic field perpendicular to the plane of the sheet. In this connection, Pavlov[9] was the first to investigate the MHD flow past a stretching sheet using the Hartman formulation. He found that the applied magnetic field and permeability cause depletion of the boundary layer thickness near the sheet surface. Chakrabarti and Gupta[10] extended the classical work of Crane[7] to include the effect of a transverse magnetic field and obtained the analytical solutions to the MHD flow over a stretching sheet.

Siddheshwar and Mahabaleshwar[3] investigated MHD viscoelastic fluid flow and heat transfer over a stretching sheet in the presence of radiations using Chandrasekhar formulation[11]. Thereafter, many authors have investigated the MHD boundary layer flows past a stretching/shrinking sheet with different control parameters and conditions[12-23].

Boundary layer flows through saturated porous media and MHD has gained significant attention in the recent times because of their promising engineering applications, such as in moisture transport in thermal insulation, ceramic processing, extraction of geothermal energy, nuclear reactor cooling systems, underground nuclear waste disposal, ground water pollution control, and filtration processes.

The objective of the present study is to investigate the boundary layer fluid flow over a permeable stretching sheet. The fluid under consideration is taken as electrically conducting and subjected to an externally imposed magnetic field normal to the stretching surface. The flow, mass transfer, and MHD effects are examined by applying the well known homotopy analysis method11 [24] (hereafter referred to as HAM)[13,25,26] to the highly nonlinear ordinary differential equations governing the flow. The HAM was originally proposed and developed by Liao[24,25,27-31]. This method is suitable for investigations of the present type as is evident from its use in the earlier studies[13] regarding ease of its applicability and accuracy for obtaining the solution of nonlinear ordinary differential equations. There is another powerful variant of HAM known as optimal homotopy perturbation method (OHPM)[32,33] available in literature to obtain analytic solutions of nonlinear differential equations. However, we will not use OHPM in the present work since HAM is sufficient for our purpose. The present study focuses on obtaining nonlinear numerical calculations relevant for visualizing the changes happening in the boundary layer flow generated due to stretching of the sheet by application of the external magnetic field, the porosity of the medium around the sheet, and suction/injection of the fluid at the sheet surface.

Mathematical Model

We consider the steady two dimensional flow of an electrically conducting, incompressible Newtonian fluid through a porous medium over a stretching sheet issuing from a slit at the origin of the rectangular coordinates as shown in Fig. 1. The sheet is assumed to be horizontal and coincides with the plane y = 0. Two equal and opposite forces are applied to the sheet to stretch it along the x-axis. The velocity of the stretching is (U(x),0), where

Figure 1

Diagram showing flow of an electrically conducting, incompressible Newtonian fluid through a porous medium over a stretching sheet issuing from the slit at the origin of the rectangular coordinates (x, y) with suction/injection. The system is subjected to a vertical magnetic field H0.

such that n > −1 is the stretching parameter; U0 and L are the characteristic scales for measuring the velocity of the stretching and distance, respectively. The porous medium is assumed to have permeabilityand is subjected to an external vertical magnetic fieldwhere κ0 and H0 denote the characteristic scales for measuring the permeability of the porous medium and magnetic field, respectively. The electrical conductivity σ of the fluid is assumed to be small so that the induced magnetic field in the fluid is weak as compared to the applied magnetic field. The sheet is permeable and subjected to suction velocity (0,V(x)) (see [[1], Ch. 11, pp. 302]), where

As a convention, V(x) > 0 implies suction while V(x) < 0 implies injection of the fluid at y = 0.

The stretching of the sheet induces a fluid velocity field (u(x, y), v(x, y)) which satisfies the equation of continuityand the boundary layer approximation of the momentum equationwhere ν is the kinematic viscosity and ρ is the fluid density. The boundary conditions for the flow are given bywhere U(x) and V(x) are as in (1) and (4).

We use similarity transformations to convert the system (5) and (6) in simpler form. Let ψ be the stream function for the flow, such that and . Introducing the Reynolds number and the similarity variable η, such thatand

where f(η) denotes the dimensionless form of the stream function, we havewhereis the dimensionless measure of suction/injection known as the mass transpiration parameter[22].

Using (9) and (10) in (5) and (6), we have the following nonlinear third order systemwhere 0 ≤ η < ∞. The stretching parameter n, the magnetic parameter appear in (12), where

The symbol is the dimensionless form of the porosity of the medium so that is the dimensionless measure of permeability. Also, the parameter is related to the Chandrasekhar number Q by , where Q = σL2H02/(ρν). Observe that Q ≥ 0 and . If no magnetic field is applied, so that . Increasing the permeability of the porous medium lowers the value of and hence of . If the me0dium is not porous, the term containing porosity is absent which correspond to and . On the other hand, if the medium is not porous and no magnetic field is applied then . Also note that V < 0 for suction and V > 0 for injection. Finally, recall that n > −1.

Solution

To solve the nonlinear system (12), we use the well known HAM[25,31] which is described as follows. For p ∈ [0,1] as the homotopy embedding parameter, we consider the following boundary value problemwhere α ≠ 0 is the unknown scalar to be determined. For p = 0, the system (14) gives the linear system f ′′′−α2f ′ = 0 and for p = 1, it is the nonlinear system (12). Now assume a solution of (14) in the form

and compare the like powers of p to obtain the following sequence of boundary value problems

The system (15) and (16) can be solved recursively starting with k = 0, where at each later stage, one needs to solve a linear boundary value problem using the solutions obtained from all of the previous stages. We have obtained the solutions of first four of these problems in closed form as follows.where

The expressions for the functions f for k > 3 are lengthy; so, we omit them. The order of approximation of f using HAM is the integer N such that . In the present work, we have obtained up to seventh order approximation of f. All numerical computations have been done using computer programming in MATLAB.

The HAM allows us to choose the scalar α appropriately. We take

where we have f'0(0) = −α from (16a). From (32), we have , which can be solved numerically to obtain an approximate value of α, depending upon N and the dimensionless parameters. For N = 1, (32) gives the following closed form expression

which is valid for all n ≥ 0 as holds since and 0 ≤ β < 2 for all n ≥ 0. For −1 < n < 0, β ∈ (−∞,0). In this case, (33) holds for . For , (33) holds either for all β ∈ (−∞,0) and or for β ≥ β0 and , where

The present formula (33) recovers the value of f ′′(0) as in Crane[7] on taking and β = 1, Pavlov[9] on taking V→0 and β = 1, Gupta and Gupta[34] for β = 1, and Hayat et al.[14] and Rashidi[15] for V→0.

To compare the first order approximation of f ′′(0) by (33) with the higher order approximations, we have obtained Fig. 2(a) which shows the variation of with V for and n = 1, 5, 10. The solid thick curves and the thin dashed curves have been obtained for N = 6 and N = 1, respectively. Clearly for each n, the formula (33) gives good approximation to f ′′(0). To quantify this approximation, we define

Figure 2

The solid thick curves and the thin dashed curves have been obtained for N = 6 and N = 1, respectively.

which is the relative difference in the values of f ′′(0) = −α obtained using first order approximation with respect to the N-th order approximation. Figure 2(b) shows the variation of E6 with V for the same parametric values as in Fig. 2(a). It is clear from the Fig. 2(a) that the relative difference E6 = 0 for n = 1, and it remains less than 5% for the other two values of n.

Table 1 shows the numerical value of α obtained with the higher order approximations for V → 0, , and n = 10. Clearly the method converges for N = 3 within the tolerance of 10−2. For N = 7 and V → 0, we have f ′′(0) ≈ −1.2349677, which is close to the corresponding numerical value −1.234875 obtained by Vajravelu and Cannon[35] and Cortell[36].

Table 1

Dependence of α on N and V for and n = 10.

	N→	1	2	3	4	5	6	7
Vc
2	α	2.5954481	2.5786412	2.5699411	2.5656288	2.5635278	2.5624998	2.5619887
1	α	1.8399457	1.8231932	1.8158511	1.8127370	1.8113701	1.8107333	1.8104193
0	α	1.2431631	1.2371618	1.2358824	1.2353717	1.2351441	1.2350297	1.2349677
−1	α	0.8399457	0.8947983	0.8499919	0.8673482	0.8628021	0.8624300	0.8639661
−2	α	0.5954481	0.5205083	0.5457379	0.5283839	0.5349155	0.5315083	0.5330783

To compare the third order approximation (N = 3) of HAM to the solution with the higher order approximations, we have plotted f(η) with respect to η in Fig. 3 for n = 10 and . We have taken since the convergence of the method is comparatively rapid for   > 0. The other parametric values are chosen to test the extreme case where the error can possibly be maximum. The three curves in each subfigure correspond to N = 3, 4, and 5, respectively as shown in the legend. Clearly, the three curves are indistinguishable for each value of V which shows that the method converges for N = 3 and justifies that the third order analytic approximation to f using HAM are sufficient to describe the solution correctly. For rest of the numerical calculations, we have taken 4 ≤ N ≤ 7.

Figure 3

f(η) vs η for n = 10 and . The different subfigures are obtained for different values of V.

Asymptotic Analysis

To understand the full parametric dependence of the present boundary layer flow, we obtain approximate analytic solution f for the following extreme cases.

Case I:

Let for some nonzero real number r when is large, where . Assume

where g < 0. The boundary condition f(0) = V gives F(0) = 0, and f ′(0) = 1 implies , so that g = r and F′(0) = 1. Now using (12), we have at the leading order

where r < 0, which for sufficiently large requires −1− 2r = 0 or r = −1/2. So, for , we get

The solution of (36) is given by . Thus, we have the following asymptotic solution for large .

The function f − V as obtained from (37) and (12) for large is shown in Fig. 4(a) for V = 1.5 and n = 1.5. In each case, the points marked * correspond to the asymptotic solution (37), while the solid curves correspond to the numerical solution of (12). The asymptotic solution is in a maximum relative error of 0.6% for , which decreases rapidly on increasing . Thus the two solutions are in good agreement for .

Figure 4

The points marked * correspond to (a) the asymptotic solution (37) and (b) the asymptotic solution (39). The other curves correspond to the numerical solution of (12).

Case II: V >> 1

For V >> 1, we take η = Vη*, where and h is a nonzero scalar. Assume

for large V, where q < 0. Using the condition f(0) = V implies G(0) = 0 and f′(0) = 1 gives G′(0) = V, so that G′(0) = 1 and h = q. Now from (12), we have

where q < 0, which for V >> 1 requires for a balance, 1 + q = 0 or q = −1. Thus, for V → ∞, we havewhich solves to G = 1 − e− = 1 − e−. Thus, in this case the asymptotic solution is given by

Figure 4b demonstrates a comparison between the asymptotic solution f − V as obtained from (39) marked * and the numerical solution of (12) (solid and dashed curves). The fixed parametric values are n = 1.5 and . Clearly, the two solutions are in very good agreement for V ≥ 4 within a relative error not exceeding 0.6%.

Results and Discussion

The numerical results have been obtained for a wide range of parameters. For most of the numerical calculations, we have taken −2 ≤ V ≤ 2, 0 ≤ n ≤ 5, and . Even higher values of are permissible since the method converges faster as becomes larger, which can be seen from (33) through (16a)–(16d), since α = |f ′′(0)| which in turn rises with for large . Also, we have taken at least the fourth order approximation of f, that is N ≥ 4 for the rest of the numerical calculations to meet the convergence issues using HAM. For obtaining the streamline patterns, we have taken 7th order approximation (N = 7) to f.

Skin friction coefficient at sheet wall

If τ denotes the shear stress near the stretching sheet due to the fluid flow, we have

Then the coefficient C of the frictional drag which is also known as the skin friction parameter is defined for x > 0 as

Note that the skin friction coefficient C depends upon f ′′(0), n, , V, and the horizontal distance x from the slit, where f ′′(0) further depends upon , n, , and V.

In the literature, f ′′(0) is generally taken as a measure of the skin friction for a fixed x, which in the present situation holds only if there is no suction/injection, that is, V ≈ 0. If this is the case, we have

We first discuss the variation of C with x in the following three cases.

Case I: n ≥ 1. Figure 5 depicts the variation of |C| with x/L for and . Each subfigure corresponds to one value of V among , , , , and . The different curves (labeled with different color and style as shown in the legend) in each subfigure correspond to different values of n ≥ 1 among 1, 1.5, 2.0, 2.5 and 3.0.

Figure 5

|C| vs x/L for and . Each subfigure corresponds to one value of V among , , and , , and . The different curves (labeled with different color and style as shown in the legend) in each subfigure correspond to different values of n ≥ 1 among 1, 1.5, 2.0, 2.5 and 3.0. Each black mark • corresponds to x0/L and each asterisk shows the local maximum corresponding to x = x.

We first explain the curves for V = −2.0 which correspond to suction. We observe from (41) that C = 0 at x = x0 (see the black marks • in Fig. 5), where

Clearly, C < 0 for x < x0 and C > 0 for x > x0. The skin friction coefficient C is negative for x in the interval (0,x0], where |C| decreases rapidly from ∞ to 0. For x > x0, C is positive and starts increasing with further increase of x till a maximum is reached at about x/L = x/L (see the points marked as in Fig. 5), where (x/L)( is a real positive root of the following polynomial equation in t

Beyond x/L = x/L, |C| falls down with further increase in x/L and approaches the value 0 at about x/L = 1. The entire pattern of these curves shifts towards right in (x/L,C)-plane on increasing n. Thus the magnitude of skin friction |C| may increase or decrease on increasing n, depending upon the horizontal distance from the slit. Similar is the variation of C with x/L for the other negative values of V. The dependence of |C| on suction can be easily understood from Fig. 6, which has been drawn for n = 1.5 and the fixed parametric values of and same as in Fig. 5. Here in Fig. 6, each curve corresponds to one value of V. It is clear from the first subfigure of Fig. 6 that on increasing suction, the point of maximum (x,|C(x)|) shifts downward in the (x/L,C) plane and for any fixed value of x/L sufficiently away from the slit, |C| decreases on increasing suction. This is due to the comparative values of the two terms −f ′′(0) and (n−1)V/(x/L)( in shear stress τ or in C as can be seen from (41).

Figure 6

|C| vs x/L for , , and n = 1.5.

We now get back to Fig. 5 to explain the case of injection (V > 0) for which the behavior of C is different from the case of suction. Here, C remains positive for all x > 0. For any fixed value of n  ≥  1, the skin friction coefficient C is a decreasing function of x/L in the considered range of x/L and C. By correlating Figs. 5 and 6 it can be observed that the skin friction decreases on increasing injection. This observation is also in accordance with the physical expectation that the injection at the sheet wall enhances the vertical component of the fluid velocity near it, which in turn results in lowering of the frictional drag near the wall.

Case II: 0 ≤ n < 1. Observe from (41) that for (n−1) < 0, the term (n−1)V is positive for suction while it is negative for injection. Here, the effect of change of n is more prominent in the case of injection than that in suction in view of the variation of |C| with x/L (see the second column of Fig. 7).

Figure 7

|C| vs x/L for , . For suction V = −1 and for injection V = 1.

Case III: −1 < n < 0. On referring back to (41) once more, one can check that here, for a fixed value of n between −1 and 0, a substantial variation of |C| with x/L can be seen. This variation is comparatively lesser in the case 0 ≤ n < 1 as can be observed by comparing the first and second columns of Fig. 7. For any fixed value of V and x/L, the magnitude of |C| decreases with x/L near the slit and increases for a distance sufficiently away from the slit as n is varied from 0 to n → −1. Note from (33) that

which is independent of suction and injection. Also note that for , |f ′′(0)| rises rapidly (hence x0 falls rapidly) as n is varied from 0 to −1. Using (44) in (41), we have the following

From (45), we have for and , as n → −1+ for all values of V and all x > 0. On comparing all the three columns of Fig. 7, it can be inferred that the unusual behavior of significant increase or decrease of |C| with x/L is apparent in the case of (i) suction for n > 1 (ii) injection for 0 < n < 1 and (iii) both suction and injection for n < 0. These observation are useful to model the nonlinearity of the underlying boundary layer flow for the two cases of suction and injection using an appropriate value of the stretching sheet parameter n.

Effect of Reynolds number

To see the effect of Reynolds number on the skin friction near the wall, we have obtained Fig. 8, which shows the variation of |C| with for various values of n and V. The fixed parametric values are and x/L = 0.5. For a fixed value of V, the different curves correspond to different values of n = 0, 0.5, 1, 1.2, 2, and 3. Observe from the expression for C in (41) that for a given value of x/L, C = 0 for , where

Figure 8

C vs for different values of n. The fixed parametric values are and x/L = 0.5.

For , C < 0, whereas C > 0 for . However, for V = −1, C remains positive for all permissible values of . For a fixed value of n in the range 0 ≤ n ≤ 1, |C| decreases continuously with in the case of suction (V = −1), while in the case of injection (V = 1), C increases with initially, attains rapidly a maximum and then it decreases slowly with . For n > 1 and V = −1, |C| decreases rapidly with to 0 and then starts increasing with and slowly attains a maximum value at larger Reynolds number, after which it decreases further on increasing . The variation for the case of injection is similar as it is for n < 1 except that on increasing n, |C| increases for a fixed value of .

Role of magnetic parameter

The effect of on C can be observed from (33). Clearly, |f′′(0)| is an increasing function of . Larger contributes towards greater resistance to the flow in x-direction. This indicates that the skin friction near the sheet wall is expected to rise for x/L away from 0 when is sufficiently large. To see this effect explicitly, we have obtained Fig. 9 which shows variation of C with x/L for , where we have taken V = −1.5 for suction and V = 1.5 for injection. In each of the case for suction/injection, the subfigures have been drawn for n = 1, 1.5, 2, and 5. In each of the subfigures, the different curves correspond to , 5, 10, 20, 50, 102, 103, and 104. Each arrow head denotes the direction of increase of the magnetic parameter .

Figure 9

C vs x/L for different values of and n. Here, V = −1.5 for suction and V = 1.5 for injection.

For the case of suction, for n = 1 and for each value of , |C| is a decreasing function of x/L. However, the variation of |C| with x/L is different for n > 1, where |C| decreases on increasing x/L. The variation of |C| for n = 1.5, 2, and 5 is similar, but the maximum of |C| decreases on increasing n.

For the case of injection, C remains positive for all the considered values of n and . In either of the cases of suction and injection, we see that |C| increases significantly on increasing , which shows enhancement in the skin friction caused by increasing the applied magnetic field and decreasing the permeability of the porous medium.

Streamline pattern

Figure 10 shows the streamline pattern in the (x/L,y/L)-plane for (a) suction (V = −1.5) and (b) injection (V = 1.5) with the fixed parametric values of and . Each subfigure corresponds to one value of n among 0, 0.5, 1, 1.5, 2.0, and 3.0. We define the boundary layer thickness to be h(x), such that as y → h(x), and for y > h(x), u = 0. So in the region y > h(x), the stream lines will be vertical. In the presence of suction at the boundary, the fluid particles in contact with the sheet start moving with the velocity of the sheet, advance towards left-upward before they are along the path normal to the sheet surface. The particles farther from the slit move through greater distance along the inclined curves. So the boundary layer thickness h(x) increases with the distance from the slit. This happens for all the considered values of n and is due to the strengthened fluid sheet velocity. These observations together show that the boundary layer flow in the presence of suction is significantly affected by changing the stretching parameter n.

Figure 10

Streamline pattern in (x/L,y/L)-plane for (a) V = −1.5 and (b) V = 1.5. The fixed parametric values are, and .

On the other hand, in the presence of injection at the sheet (Fig. 10(b)), the fluid particles in contact with the sheet start moving with the velocity of the sheet, advance a little towards right-upward and attain paths normal to the sheet. Here, the boundary layer thickness h(x) is smaller than in the case of suction, where at a fixed distance from the slit, h(x) diminishes with increase in n. For each value of n, h(x) increases with the distance from the slit.

To observe the effect of porosity and applied magnetic field on the boundary layer flow pattern, we have obtained Fig. 11 which shows the streamline pattern in the xy-plane, for different values of n, V = −1 and . The various patterns (a)–(d) correspond to , 1, 2, and 5 respectively. Clearly, for a fixed value of n, the boundary layer thickness decreases on incrementing from 0.5 to 5.

Figure 11

Streamline pattern in (x/L,y/L)-plane for V = −1 and .

It will be useful to see the explicit dependence of the velocity profiles on for the two cases of suction and injection. This has been shown in Fig. 12, which shows the variation of the two components of velocity, u/U(x) and v/V(x) with η for n = 1.5 and . Here, we have taken V = −1 for suction and V = 1 for injection. Each curve is drawn for one value of among 0.5, 1, 2, 5, 10, and 20, and the arrowhead in each subfigure shows the direction of increase of . Clearly, on increasing , the horizontal velocity profile shifts downwards in the (η,u/U(x))-plane under suction as well as injection. On the other hand, the profile for the velocity component v/V(x) goes up on increasing for the case of suction, while opposite to this variation occurs under injection. On correlating all the four subfigures, we see that under suction, the applied magnetic field and the permeability of the porous medium tend to enhance the fluid velocity in the direction of suction, while it hinders the fluid motion along the sheet. Under the case of injection, the applied magnetic field hinders both the vertical as well as horizontal movement of the fluid in the vicinity of the sheet wall. The value of y corresponding to least value η for which u/U(x) becomes zero, is the boundary layer thickness. So increase in decreases the boundary layer thickness.

Figure 12

Velocity profiles u/U(x) and v/V(x) vs η for suction (V = −1) and injection (V = 1) and .

Conclusions

The present work generalizes the classical work of Crane[7], Pavlov[9], Gupta and Gupta[34], Hayat et al.[14], and Rashidi[15] on the flow of Newtonian fluid driven by stretching sheet with external magnetic field through porous medium with suction/injection. The underlying dynamical system is described by the nonlinear boundary layer equations, which are transformed into a system of nonlinear ordinary differential equations via similarity transformations. The resulting system is solved analytically and numerically using highly efficient HAM.

The second order analytical solution suggests that the axial and transverse velocities of the boundary layer flow can be approximated by a finite linear combination of the basis functions {exp{−kη}}, where the coefficients are functions of the dimensionless parameters n, V, , and . The method allows a very good approximation of the second derivative f ′′(0) through the formula (33). The numerical results have been obtained for a wide range of parametric values.

An increase in the values of the results in pulling down of the horizontal velocity profiles of the flow. This may be due to an increase in the Lorentz and Darcy forces in the fluid, which oppose the fluid motion. An increase in the parameter leads to rise in the skin friction in all the cases of suction and injection.

An increase in the numerical values of V and results in enhancement in the horizontal component of the fluid velocity due to which the laminar MHD boundary layer gets depleted.

The present numerical results are useful in modeling the nonlinearity of the underlying boundary layer flow for the two cases of suction and injection using an appropriate value of the stretching sheet parameter n. More precisely, it may be useful to take n > 1 for the case of suction, 0 < n < 1 for injection, and −1 < n < 0 both for suction and injection in order to observe the modifications occurring in the boundary layer flow.

The asymptotic behavior of the solution corresponding to the boundary layer flow near the stretching sheet surface in the presence of suction/injection is discussed separately for (i) large and (ii) V large. The approximated analytic solution is found to be in excellent good agreement with the numerical solution obtained without any approximation for of the order of 102 and V of the order of 101.

Since the linear stretching can be acknowledged in practice only with great care and meticulous effort, the consideration of the nonlinear stretching sheet is more useful.