Advances in Mathematical Physics
Volume 2016 (2016), Article ID 4253649, 8 pages
http://dx.doi.org/10.1155/2016/4253649

Research Article

MHD Flow due to the Nonlinear Stretching of a Porous Sheet

Tarek M. A. El-Mistikawy

Department of Engineering Mathematics & Physics, Faculty of Engineering, Cairo University, Giza 11221, Egypt

Received 13 September 2015; Accepted 20 December 2015

Academic Editor: Luigi C. Berselli

Copyright © 2016 Tarek M. A. El-Mistikawy. 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 MHD flow due to the nonlinear stretching of a porous sheet is investigated. A closed form solution is obtained when the stretching rate is inversely proportional to the distance from the origin. Otherwise a uniformly valid asymptotic expansion, for large magnetic interaction number , is developed. It coincides with a homotopy perturbation expansion for the problem. The asymptotic/homotopy perturbation expansion gives results in excellent agreement with accurate numerical results, for large as well as small values of . For large , the expansion, being asymptotic, needs a small number of terms, regardless of the mass transfer rate or the degree of nonlinearity. For small , the expansion is a homotopy perturbation one. It needs considerably increasing number of terms with higher injection rates and/or with stretching rates approaching the inverse proportionality. It may even fail.

1. Introduction

The problem of the MHD flow due to the nonlinear stretching of a nonporous sheet was first formulated and analyzed by Chiam [1]. He chose a special magnetic field that allowed the problem to be cast in a self-similar form. He also explored the asymptotic solution for large magnetic interaction number and obtained a three-term expansion for the stream function that was marred by the presence of secular terms.

The problem has recently come under extensive study. The aim has been to construct approximate expansion solutions through newly developed analytical methods, with the help of computer symbolic manipulation software facilities. Hayat et al. [2] has used Adomian’s decomposition method to construct an expansion in powers of the similarity coordinate . Yildirim and Sezer [3], Fathizadeh et al. [4], and Sushila et al. [5] have arrived at the same power expansion through variants of the homotopy perturbation method. The power expansion has given accurate results only at a narrow range of . This range has been extended through the use of Padé approximants in [2, 3].

Raftari et al. [6] have constructed a different homotopy problem, on which they have imposed an unnecessary restriction, and have obtained a two-term homotopy perturbation expansion involving secular terms.

In this paper, the problem of the MHD flow due to the nonlinear stretching, according to a power law of index , of a porous sheet is studied. When , a closed form solution is obtained. When , a uniformly valid asymptotic expansion as is obtained and is shown to be identical to a homotopy perturbation expansion of the problem. The present asymptotic/homotopy perturbation expansion alleviates problems facing previous works on Chiam’s problem [1–6]. Consequences of the expansion being asymptotic or a homotopy perturbation one are assessed with regard to the number of terms retained, and the accuracy achieved.

2. Formulation of the Problem

An electrically conducting fluid is driven by a symmetrically stretching insulated sheet, lying along the -axis. The stretching rate is proportional to . The sheet is porous allowing fluid flow of rate proportional to in the normal -direction, along which a magnetic field of strength proportional to is applied. The power law index may be positive or negative.

For , the boundary layer equations write [1]

with the boundary conditions

() are the velocity components in the () directions, () are (stretching, suction, and magnetic field) coefficients, () are the fluid (density, kinematic viscosity, and electric conductivity), respectively.

The following self-similar problems can be formulated, with denoting the stream function.

3. Case of

We introduce the following similarity variables and parameters:

The problem takes the form

The last condition corresponds to at the sheet.

The problem admits the following solution:

Thus, no such self-similar solution is possible, when .

Differentiation gives

The surface shear is represented by

The flow is of a source type. The stream lines are lines of constant , that is, rays emanating from the origin. The source strength measured by the rate of fluid discharge (to the first quadrant) is represented by

This discharge must come from a slit at the origin.

4. Case of

We introduce the following similarity variables and parameters:

The problem takes the form

Chiam’s nonporous sheet problem [1] corresponds to and .

The case of ( and ) admits the exact solution

where is the positive root of the quadratic equation

This solution was obtained by Chakrabarti and Gupta [7]. It degenerates to that of Andersson [8] for a nonporous sheet , to that of P. S. Gupta and A. S. Gupta [9] in the absence of the magnetic field , and to that of Crane [10] when .

4.1. Asymptotic Limit as

We are interested in the limiting behavior of the flow as the magnetic interaction number grows indefinitely. In the limit as , having the fixed value at . The electromagnetic force term is, thus, of leading order. It can be balanced by the viscous term in a contracting region in which .

Introducing the stretched variable , the problem becomes

where primes denote differentiation with respect to .

The straightforward expansions of would involve secular terms. To overcome this difficulty we follow the lead of He [11] and strain the coefficient 1 of . To put this in effect, we rewrite (8a) as

Introducing the asymptotic expansions

and equating like powers of , we arrive at problems for . The problem for is

with the solution

The problem for is, then,

with the solution

The problem for is

The right hand side of (15a) produces a secular term, the removal of which requires the regularity condition

leading to the solution

Likewise, the problem for gives the regularity condition

and the solution

The problem for , is

In general, the solution for , is of the form

Expressions to determine ’s are given in Appendix. It is noted that, of all ’s, only appears in the expressions for ’s. Moreover, the following properties of ’s can be proved by induction:

(a)for , has as a factor,(b)for even , has as a factor,(c)for , is inversely proportional to ; that is,

where ’s are given by the same expressions in Appendix for ’s but with set equal to unity.

For , the regularity condition is

which combines with (9b) and (10b) to give the following equation for :

In the case of (), for due to property (a) of . Therefore,

Equations (23a) and (23b) are identical to the exact solution expressed by (7a) and (7b), as we set .

4.2. Homotopy Perturbation Solution

We now seek a homotopy perturbation solution to the problem described by (8a), (8b), (8c), and (8d). To that end, we formulate an associated homotopy problem

where is the embedding parameter, varying from 0 to 1.

Setting , the homotopy problem (24a), (24b), (24c), and (24d) acquires the same form as the physical problem (8a), (8b), (8c), and (8d) and coincides with it at the final step of the homotopy perturbation procedure, when we set .

As both problems are treated as perturbation problems, the resulting expansions will be identical.

The present homotopy perturbation expansion is free from secular terms and may involve as many terms as one chooses. Allowing a secular term, Raftari et al. [6] obtained a two-term homotopy perturbation expansion, for the case of , that imposed a restriction on and , namely, , and that could not be extended to include further terms.

Through variants of the homotopy perturbation method, Yildirim and Sezer [3], Fathizadeh et al. [4], and Sushila et al. [5] obtained the same polynomial solution to the case of that had been obtained earlier by Hayat et al. [2] via Adomian’s decomposition method. Results presented in [4, 5] showed deviation from numerically calculated results, beyond a narrow range of . This range was extended through the use of Padé approximants in [2, 3]. In contrast, the present homotopy perturbation solution, involving negative exponentials, is valid for all values of .

Keeping enough terms of the present homotopy perturbation expansion, it is expected to get accurate results for all values of . However, the expansion, being an asymptotic one as as well, should require few terms to give accurate results for large .

5. Numerical Procedures

To put the above expansions into use, we cut short the summations in (10a) and (21) to the same finite upper limit . The following iterative procedure is then adopted.

(1)For the given and , determine all ’s using the formulas in Appendix.(2)Choose a starting value for .(3)Calculate ’s using (20).(4)Introduce ’s into (22) and solve it for a new .(5)If convergence of is not reached, that is, , go to step  .(6)Otherwise, evaluate and its derivatives for the given .

Using (20), we can transfer (22) to an algebraic equation of degree in and can derive a Newton iterative formula for to be used in step  ; namely,

The two quantities of physical significance are the surface shear represented by

and the rate of flow entrainment represented by

These two quantities, together with , will be the focus of our calculations.

The results of the abovementioned procedure are tested against corresponding numerical results. The problem described by (6a), (6b), (6c), and (6d) is solved numerically using Keller’s two-point, second-order-accurate, finite-difference scheme [12]. A uniform step size is used on a finite domain . The value of is chosen sufficiently large in order to insure the asymptotic satisfaction of the far-field condition (6d). The nonlinear terms are quasi-linearized, and an iterative procedure is implemented, terminating when the maximum error in and becomes less than .

The accuracy of the numerical procedure was established by comparison with the exact solution for with different combinations of and , both producing identical results to, at least, the tenth significant figure.

The numerical solution calculates as the positive solution of the quadratic equation

which is the result of substitution of (26b) into (22). This guides the choice of a suitable starting value for , in step   of the iterative procedure described above.

6. Results

In the results presented below for the expansion solutions, two factors are of concern: the choice of the cutoff number and the accuracy in comparison to the numerical solutions.

First, it should be noted that has to be ≥3 for the effect of to be retained. As increases, progressive “saturation” in the significant figures of is achieved. Where to stop depends on the flow parameters , , and . To illustrate the dependence on , consider the case of (), . Table 1 presents results at consecutive ’s for , 1, and 0.01. (We proceed in odd values of , to avoid programming complexities.) It is noted that, for the larger , saturation in 6 significant figures is attained at . The intermediate and smaller ’s achieve the same level of saturation at and 21, respectively.

Table 1: Dependence of on [ (), ].

100
5
1.00498 3457
−10.0664 3912
0.0994219 5793

7
1.00498 3454
−10.0664 3909
0.0994219 5832

9
1.00498 3454
−10.0664 3909
0.0994219 5832

Numerical
1.00498 3447
−10.0664 2517
0.0994218 2036

1
9
1.39844 8416
−1.52527 1655
0.683382 0804

11
1.39845 4500
−1.52527 4499
0.683379 5495

13
1.39845 4580
−1.52527 4777
0.683379 4078

Numerical
1.39845 3909
−1.52527 3038
0.683378 5009

0.01
19
9.54667 9245
−1.15300 4023
0.944193 2777

21
9.54668 0372
−1.15300 4063
0.944193 1944

23
9.54668 0386
−1.15300 4075
0.944193 1820

Numerical
9.54666 8853
−1.15300 3010
0.944192 0274

Also presented in Table 1 are the corresponding numerical results. As expected, for the larger , few terms of the expansion, which is asymptotic, give excellent agreement with the numerical solution. For the smaller ’s, the expansion is no longer asymptotic but should be considered a homotopy perturbation one. Its agreement with the numerical solution is as excellent.

Does this agreement extend to other values of and ? The answer is given in Tables 2, 3, and 4, in which, for each indicated value of or , the first line gives the expansion solution and, if possible, the value of for saturation in 6 significant figures, while the second line gives the numerical solution.

Table 2: Validity of the expansion solution at different values of and , ().

100
1.0
0.01

−2
−9.11655
5
−1.89021
—

—
—

—

−9.11654

−1.89021
−0.844839

−0.714243
−0.558696

0.249190

−1
−9.57913
5
−0.895514
−1.11742
29
−0.0471147
—

—

−9.57911

−0.895515
−1.11742

−0.0471154
−0.778042

0.504566

0
−10.0664
5
0.0994220
−1.52527
13
0.683379
−1.15300
17
0.944193

−10.0664

0.0994218
−1.52527

0.683379
−1.15300

0.944192

1
−10.5786
5
1.09460
−2.09551
9
1.48920
−1.74219
11
1.59722

−10.5786

1.09460
−2.09551

1.48920
−1.74219

1.59721

2
−11.1155
5
2.09003
−2.81024
11
2.36110
−2.50877
11
2.40646

−11.1154

2.09003
−2.81023

2.36110
−2.50877

2.40646

Table 3: Validity of the expansion solution at different values of and , .

100
1
0.01

−1.5
−9.96665
5
0.0999166
−0.652972
269
0.912433
—

—

−9.96672

0.0999173
−0.652978

0.912440
0.696751

1.70859

−1
−9.98335
5
0.0998335
−0.851110
113
0.851110
—

—

−9.98340

0.0998340
−0.851114

0.851114
−0.0140447

1.40447

−0.5
−10.0000
5
0.0997506
−1.01756
53
0.804077
—
—
—

−10.0001

0.0997510
−1.01756

0.804080
−0.379674

1.24358

0
−10.0167
5
0.0996680
−1.16333
33
0.766078
−0.634631
147
1.13644

−10.0167

0.0996683
−1.16333

0.766080
−0.634630

1.13645

0.5
−10.0333
5
0.0995857
−1.29435
19
0.734312
−0.835754
55
1.05725

−10.0333

0.0995859
−1.29435

0.734313
−0.835755

1.05725

1
−10.0499
3
0.0995037
−1.41421
3
0.707107
−1.004988
3
0.995037

−10.0499

0.0995037
−1.41421

0.707107
−1.004988

0.995037

1.5
−10.0664
5
0.0994220
−1.52527
13
0.683379
−1.15300
17
0.944193

−10.0664

0.0994218
−1.52527

0.683379
−1.15300

0.944192

Table 4: Least possible (in whole numbers) for , .

2.5
3
−0.984756
17
0.850785

−0.984762

0.850788

5
4
−0.439138
69
0.990757

−0.439173

0.990754

7.5
6
−0.742926
81
0.694277

−0.742974

0.694285

10
8
−0.944398
93
0.569562

−0.944470

0.569572

Table 2 includes cases of suction and injection . Agreement with the numerical solutions to at least five significant figures is observed. At the larger , the asymptotic nature of the expansions allows a fixed small value of , regardless of the value of . At the smaller ’s, the homotopy perturbation expansions are aided by suction and progressively hurt by injection, to the extent that no solution is possible as indicated by the void entries. Either saturation in enough significant figures could not be attained, or convergence of the iterative process for could not be reached. As an illustration, consider the case of and . The following are some results that show no tendency toward saturation: , , , , , …, , , , , , , …, , , , , …, , , , , …, . The ellipses indicate -ranges of no convergence to a positive . When , and 62, converges to the same negative value . Otherwise, the left hand side of (22) drawn versus , in the range indicated by the numerical solution, does not cross the zero line.

Table 3 includes cases of stretching sheets with progressively increasing rate and progressively decreasing rate , which are associated with . At the larger , the asymptotic expansions allow a fixed small value of , regardless of the value of . At the smaller ’s, the homotopy perturbation expansions saturate readily when . However, as decreases negatively they tend toward saturation with painful slowness, hence the void entries. As for accuracy, the expansion and numerical solutions are in excellent agreement, differing, at worst, by 1 in the fifth significant figure.

Other cases of progressively decreasing stretching rate , which are associated with , are given in Table 4. The indicated values of are the least whole numbers below which no solution was possible either numerically or through the expansions. As gets closer to ( increases), the least ’s get higher, the homotopy perturbation solutions require more terms for saturation, and their agreement with the numerical solutions deteriorates.

(It should be noted that the results presented in [1–4] for the case of are incorrect, as they were associated wrongly with .)

7. Conclusion

The problem of the MHD flow due to a nonlinearly stretching porous sheet has been investigated. A closed form solution has been obtained when the stretching rate is inversely proportional to distance (the power law index ). Otherwise, uniformly valid asymptotic expansion, for large magnetic interaction numbers , has been developed. This expansion has been shown to coincide with a homotopy perturbation expansion for the problem. The expansion gives results in excellent agreement with accurate numerical results, for large as well as small values of . For large , the expansion being asymptotic, a small number of terms are needed regardless of the value of or the suction rate . For small , the expansion is a homotopy perturbation one. It needs considerably increasing number of terms as or grows negatively and may even fail.

The problem belongs to the class of problems identified below, whose asymptotic expansions and homotopy perturbation expansions are identical. The study helps realize the potential and limitations of asymptotic and homotopy perturbation methods.

Let be a sought after vector function which is defined on a domain with boundary and governed by the following physical problem

where is a small parameter and and are left hand side operators, while and are right hand side operators.

We form the associated homotopy problem

where is the embedding parameter, varying from 0 to 1.

Setting , the homotopy problem acquires the same form as the physical problem and coincides with it at the final step of the homotopy perturbation method, when we set .

As both problems are treated as perturbation problems, the resulting expansions will be identical.

Appendix

The following expressions determine the coefficients in , . They involve summations, the limits of which may contain the forms , , , , or where , and are integers. When the upper limit is less than the lower limit, the summation is to be zeroed.

Given

then for even and

and for odd and

with

The boundary conditions (18c) and then (18b) determine and then .

The expressions for ’s are identical to those for ’s but with .

Conflict of Interests

The author declares that there is no conflict of interests regarding the publication of this paper.

References

T. C. Chiam, “Hydromagnetic flow over a surface stretching with a power-law velocity,” International Journal of Engineering Science, vol. 33, no. 3, pp. 429–435, 1995.

T. Hayat, Q. Hussain, and T. Javed, “The modified decomposition method and Padé approximants for the MHD flow over a non-linear stretching sheet,” Nonlinear Analysis: Real World Applications, vol. 10, no. 2, pp. 966–973, 2009.

A. Yildirim and S. A. Sezer, “Non-perturbative solution of the magnetohydrodynamic flow over a nonlinear stretching sheet by homotopy perturbation method-padé technique,” Zeitschrift für Naturforschung, vol. 65, no. 12, pp. 1106–1110, 2010.

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.

Sushila, J. Singh, and Y. S. Shishodia, “An efficient analytical approach for MHD viscous flow over a stretching sheet via homotopy perturbation sumudu transform method,” Ain Shams Engineering Journal, vol. 4, no. 3, pp. 549–555, 2013.

B. Raftari, S. T. Mohyud-Din, and A. Yildirim, “Solution to the MHD flow over a non-linear stretching sheet by homotopy perturbation method,” Science China: Physics, Mechanics and Astronomy, vol. 54, no. 2, pp. 342–345, 2011.

A. Chakrabarti and A. S. Gupta, “Hydromagnetic flow and heat transfer over a stretching sheet,” Quarterly of Applied Mathematics, vol. 37, no. 1, pp. 73–78, 1979.

H. I. Andersson, “An exact solution of the Navier-Stokes equations for magnetohydrodynamic flow,” Acta Mechanica, vol. 113, no. 1–4, pp. 241–244, 1995.

P. S. Gupta and A. S. Gupta, “Heat and mass transfer on a stretching sheet with suction and blowing,” The Canadian Journal of Chemical Engineering, vol. 55, no. 6, pp. 744–746, 1977.

L. J. Crane, “Flow past a stretching plate,” Journal of Applied Mathematics and Physics, vol. 21, no. 4, pp. 645–647, 1970.

J.-H. He, “Addendum: new interpretation of homotopy perturbation method,” International Journal of Modern Physics B, vol. 20, no. 18, pp. 2561–2568, 2006.

H. B. Keller, “Accurate difference methods for linear ordinary differential systems subject to linear constraints,” SIAM Journal on Numerical Analysis, vol. 6, pp. 8–30, 1969.