Mathematical Problems in Engineering

Volume 2012 (2012), Article ID 454023, 26 pages

doi:10.1155/2012/454023

Research Article

Homotopy Analysis Solution of Hydromagnetic Mixed Convection Flow Past an Exponentially Stretching Sheet with Hall Current

Mohamed Abd El-Aziz1,2 and Tamer Nabil1,3

¹Mathematics Department, Faculty of Science, King Khalid University, Abha, P.O. Box 9004, Saudi Arabia

²Mathematics Department, Faculty of Science, Helwan University, P.O. Box 11795, Cairo, Egypt

³Basic Science Department, Faculty of Computers and Informatics, Suez Canal University, Ismailia, Egypt

Received 9 May 2012; Revised 18 July 2012; Accepted 11 September 2012

Academic Editor: Claude Lamarque

Copyright © 2012 Mohamed Abd El-Aziz and Tamer Nabil. 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 effect of thermal radiation on steady hydromagnetic heat transfer by mixed convection flow of a viscous incompressible and electrically conducting fluid past an exponentially stretching continuous sheet is examined. Wall temperature and stretching velocity are assumed to vary according to specific exponential forms. An external strong uniform magnetic field is applied perpendicular to the sheet and the Hall effect is taken into consideration. The resulting governing equations are transformed into a system of nonlinear ordinary differential equations using appropriate transformations and then solved analytically by the homotopy analysis method (HAM). The solution is found to be dependent on six governing parameters including the magnetic field parameter M, Hall parameter m, the buoyancy parameter , the radiation parameter R, the parameter of temperature distribution a, and Prandtl number Pr. A systematic study is carried out to illustrate the effects of these major parameters on the velocity and temperature distributions in the boundary layer, the skin-friction coefficients, and the local Nusselt number.

1. Introduction

In many industrial manufacturing processes, the problem of flow and heat transfer in two-dimensional boundary layer on a continuous stretching surface, moving in an otherwise quiescent fluid medium, has attracted considerable attention during the last few decades. Examples may be found in continuous casting, glass-fiber production, hot rolling, wire drawing, paper production, drawing of plastic films, metal and polymer extrusion, and metal spinning. In recent years MHD flow problems have become more important in industry, since many metallurgical processes involve the cooling of continuous strips or filaments. By drawing them in an electrically conducting fluid in the presence of a magnetic field, the rate of cooling can be controlled. Another important application of hydromagnetics to metallurgy is the purification of molten metals from nonmetallic inclusions by the application of a magnetic field. The flow past a moving or stretching surface in an ambient fluid differs from that of the classical Blasius problem of flow past a stationary surface. The moving surface sucks the fluid and pumps it back in the downstream direction. Consequently, both the surface shear stress and the heat transfer are significantly enhanced. Sakiadis [1] firstly studied the boundary layer flow over a stretched surface moving with a constant velocity. Crane [2] extended this concept to a stretching sheet with linearly varying surface speed and presented an exact analytical solution for the steady two-dimensional stretching of a surface in a quiescent fluid. This problem has later been extensively studied in various directions including porous surface, non-Newtonian fluids, magnetohydrodynamic fluid, heat transfer, mass transfer, porous medium, and slip effects. Some interesting investigations are mentioned in the references [3–9]. It is known that the buoyancy force can produce significant changes in the velocity and temperature distribution and, hence, in the heat transfer rate from the surface. The effect of buoyancy force over continuous moving surfaces through an otherwise quiescent fluid was investigated by Chen and Strobel [10] and Fan et al. [11] for horizontal surfaces, by Chen [12], Ali and Al-Yousef [13, 14] and Abd El-Aziz and Salem [15] for vertical surfaces and by Moutsoglou and Chen [16], Strobel and Chen [17], and Chen [18] for vertical and inclined surfaces. A new dimension is added to the study of flow and heat transfer in a viscous fluid over a stretching surface by considering the effect of thermal radiation. We know that the radiation effect is important under many nonisothermal situations. If the entire system involving the polymer extrusion process is placed in a thermally controlled environment, then radiation could become important. Radiative heat transfer flow is very important in manufacturing industries for the design of reliable equipment, nuclear plants, gas turbines, and various propulsion devices for aircraft, missiles, satellites, and space vehicles. Also, the effect of thermal radiation on the convective flows is important in the content of space technology and processes involving high temperature. The knowledge of radiation heat transfer in the system can perhaps lead to a desired product with a sought characteristic. However, in all the stretching sheet problems (both hydrodynamic and hydromagnetic) mentioned earlier, radiation effect has not been considered. Extensive literature that deals with flows in the presence of radiation effects is now available (see e.g, England and Emery [19], Gorla and Pop [20], Raptis [21], and Abd El-Aziz [22–24]). All the above-mentioned investigations were limited to a continuous surface moving with a constant, linear, or nonlinear velocity. However very little attention is given to the flow over an exponentially stretching sheet. Magyari and Keller [25] analyzed the steady free convection flow and heat transfer from an exponentially stretching vertical surface. Elbashbeshy [26] examined the flow and heat transfer characteristics by considering exponentially stretching continuous surface. Viscoelastic boundary layer flow over an exponential stretching continuous sheet has been examined by [27, 28]. Recently, Partha et al. [29] studied the effect of viscous dissipation on the mixed convection heat transfer from an exponentially stretching surface. Very recently, Abd El-Aziz [30] examined the problem of viscous dissipation effect on mixed convection flow of a micropolar fluid over an exponentially stretching sheet. When the conducting fluid is an ionized gas, and the strength of the applied magnetic field is large, the conductivity normal to the magnetic field is reduced due to the free spiraling of electrons and ions about the magnetic lines of force before suffering collisions and a current is induced in a direction normal to both electric and magnetic fields. This phenomenon is called Hall effect. In all of the previous investigations, the Hall term was ignored in applying Ohm's law as it has no marked effect for small and moderate values of the magnetic field. When the medium is rarefied or if a strong magnetic field is present, the conductivity of the fluid is anisotropic and the effect of Hall current cannot be neglected. The study of MHD viscous flows with Hall current has important applications in problems of power generators and Hall accelerators as well as flight magnetohydrodynamics. The current trend for the application of magnetohydrodynamics is toward a strong magnetic field (so that the influence of electromagnetic force is noticeable) and towards a low density of the gas (such as in space flight and in nuclear fusion research) [31]. Under these conditions the Hall current becomes important. With the above understanding, Pop and Watanabe [32] studied the problem of free convection flow of an electrically conducting viscous fluid without neglecting the Hall effect. Abo-Eldahab and Abd El-Aziz [33] analyzed the problem of MHD free convection flow of an electrically conducting and heat generating/absorbing fluid past a semi-infinite vertical plate taking into consideration the effects of Hall and ion-slip currents. Megahed et al. [34] investigated the effects of Hall currents on a steady free convection flow and mass transfer past a semi-infinite plate past a viscous incompressible electrically conducting fluid using similarity analysis. Abo-Eldahab and Abd El-Aziz [35] studied the effect of Hall current and Ohmic heating on mixed convection boundary layer flow of a micropolar fluid from a rotating cone. Abo-Eldahab and Abd El-Aziz [36] presented an analysis for the effects of viscous dissipation and Joule heating on the flow of an electrically conducting and viscous incompressible fluid past a semi-infinite plate in the presence of a strong transverse magnetic field and heat generation/absorption. Saha et al. [37] analyzed the effect of Hall current on the steady, laminar, natural convection boundary layer flow of MHD viscous and incompressible fluid from a semi-infinite heated permeable vertical flat plate. However, relatively little work has been done on the effect of Hall current on the boundary layer flow of an electrically conducting viscous fluid past a stretching surface. The effect of Hall current on a steady, laminar, hydromagnetic boundary layer flow of an electrically conducting and heat generating/absorbing fluid along a stretching sheet is considered by Abo-Eldahab et al. [38] and Salem and Abd El-Aziz [39]. Recently, Abd El Aziz [40] studied the unsteady flow and heat transfer over a stretching surface with Hall effect. To the best of our knowledge, no analytical or numerical results have been reported for the effects of Hall current on the convective heat transfer past an exponentially stretching sheet in the presence of thermal radiation. Recently, Sajid and Hayat [41] gave analytical solution for the problem of radiation effect on heat transfer of a viscous fluid over an exponentially stretching sheet using the homotopy analysis method (HAM) (Liao [42, 43]). This method is based on a fundamental concept in topology, that is, homotopy [44] which is widely used in numerical techniques (Chan and Keller [45], Grigolyuk and Shalashilin [46]). In this paper, HAM is employed to find an analytical solution of the problem of Hall effect on MHD mixed convection flow past a vertical stretching surface in the presence of radiation.

2. Analysis

Consider the steady mixed convection boundary layer flow past a heated semi-infinite vertical wall stretching with velocity and a given temperature distribution moving through a quiescent viscous, incompressible, and electrically conducting fluid with constant temperature . The positive coordinate is measured along the stretching sheet in the direction of motion and the positive coordinate is measured normal to the sheet in the outward direction toward the fluid. The leading edge of the stretching sheet is taken as coincident with -axis (see Figure 1). The fluid is considered to be a gray, absorbing-emitting radiation but nonscattering medium and the Rosseland approximation is used to describe the radiative heat flux in the energy equations. Also, the flow is subjected to a strong transverse magnetic field with a constant intensity along the positive -direction. The magnetic Reynolds is assumed to be small enough () so that the induced magnetic field can be neglected. This assumption is justified since the magnetic Reynolds number is generally very small for weakly ionized gases [47]. In general, for an electrically conducting fluid, Hall current affects the flow in the presence of a strong magnetic field. The effect of Hall current gives rise to a force in the -direction, which induces a crossflow in that direction and hence the flow becomes three-dimensional. To simplify the problem, we assume that there is no variation of flow quantities in -direction. This assumption is considered to be valid if the surface is of infinite extent in the -direction. If the Hall term is retained in generalized Ohm's law, then the following expression holds [48]:

where is the current density vector, is the velocity vector, is the magnetic induction vector, is the cyclotron frequency of electrons, is the electron collision time, and is the electrical conductivity. The ion-slip and thermoelectric effects are not included in (2.1). Further it is assumed that and where and are cyclotron frequency and collision time for ions, respectively. In addition, we consider the case of a short circuit problem in which the applied electric field and for partially ionized gas, the electron pressure gradient may be neglected. Assuming the plate to be electrically nonconducting, the generalized Ohm's law under the above conditions gives everywhere in the flow. Hence under these assumptions, equating the and components in (2.1) and solving for the current density components and , we have

Here , , and are the -, -, and -components of the velocity vector and is the Hall parameter.

Figure 1: Schematic representation of the physical model and coordinates system.

Finally, we assume the fluid is isotropic, homogeneous, and has the scalar constant viscosity and electric conductivity. Under the above assumptions and invoking the Boussinesq approximation, the boundary layer equations governing the flow and heat transfer of a viscous and incompressible fluid can be written as

where is the fluid temperature, is the kinematic coefficient of viscosity with being the fluid viscosity and is the fluid density, is the thermal diffusivity with being the fluid thermal conductivity and is the heat capacity at constant pressure, and is the radiative heat flux.

The radiative heat flux under Rosseland approximation [41] has the form

where is the Stefan-Boltzmann constant and is the mean absorption coefficient.

We assume that the temperature difference within the flow is sufficiently small such that may be expressed as a linear function of temperature. This is accomplished by expanding in a Taylor series about and neglecting higher-order terms, thus

In view of (2.7) and (2.8), (2.6) reduces to

The associated boundary conditions are

The stretching surface is assumed to have an exponential velocity distribution of the form [30]:

Here is a constant and is the reference length. The exponential velocity (2.11) is valid only when which occurs very near to the slot [28].

Also the surface temperature of the stretching sheet is assumed to be in the form:

where and are the parameters of temperature distribution on the stretching surface and is the ambient temperature. The special case corresponds to the well-known but important particular case of the isothermal plate.

Introduce the dimensionless variables , , , and as follows:

In (2.14) the stream function is defined by and , such that the continuity equation (2.3) is satisfied automatically and is given by (2.12).

In terms of these new variables, the velocity components can be expressed as

where prime denotes ordinary differentiation with respect to . Now substituting (2.7)–(2.10) in (2.4), (2.5), and (2.9) we obtain the following locally similar ordinary differential equations:

The boundary conditions (2.10) then turn into

where is the mixed convection or buoyancy parameter with being the local Grashof number and is the local Reynolds number, is the local magnetic field parameter, is the radiation parameter, and is the Prandtl number.

We notice that when , (2.18), (2.19), and (2.20) are uncoupled and a purely forced convection situation results. In this case, the flow field is not affected by the thermal field. Also, in the absence of Hall (), buoyancy (), and radiation () effects, we notice that (2.18)–(2.20) reduce to those of Magyari and Keller [25].

For practical applications, the major physical quantities of interest are the local skin-friction coefficient in the -direction

local skin-friction coefficient in the -direction,

and the local Nusselt number,

3. HAM Solutions

In order to solve the governing nonlinear system (2.18)–(2.20) subject to the boundary conditions (2.21) we employ the homotopy analysis method [43]. According to the boundary conditions (2.21), it is reasonable to assume that , , and can be expressed by the following set of base functions:

such that

where , and are constant coefficients. The rule of solution expression provides us with a starting point. It is under the rule of solution expression that initial approximations, auxiliary linear operators, and the auxiliary functions are determined. So, according to the rule of solution expression, we choose the initial guess and auxiliary linear operator [49–51] in the following forms:

in which the auxiliary linear operators have the following properties:

where are constants. Let denote the embedding parameter and let , , and indicate nonzero auxiliary parameters. We then construct the following equations.

3.1. Zeroth-Order Deformation Equations

Consider

where the nonlinear operators , and , respectively, are

Obviously when and , the above zeroth-order deformation equations (3.5) have the solutions:

Expanding , , and in Taylor's series with respect to , we have

where

3.2. th-Order Deformation Equations

Differentiating the zeroth-order deformation equations (3.5) -times with respect to , then setting , and finally dividing them by , we obtain the th-order deformation equations:

subject to the boundary conditions:

where

If we let , , and as the special solutions of (3.11), the general solution is given by

where the integral constants are determined by the boundary conditions (3.12). In this way it is easy to solve the linear nonhomogeneous equations (3.11) by using Maple one after the other in the order .

3.3. Convergence of the Analytic Solution

Liao [43] showed that whenever a solution series converges, it will be one of the solutions of considered problem. Therefore, it is important to ensure that the solutions series are convergent. The solutions series (3.9) contain the nonzero auxiliary parameters , , and , which can be chosen properly by plotting the so-called -curves to ensure the convergence of the solutions series and rate of approximation of the HAM solution. To plot the convergence curve of one of the auxiliary parameters , , and we must choose the values of other two parameters. In the present work, the optimal HAM [52] is used to obtain the optimal values of the auxiliary parameters by means of the minimum of the residual squares of the governing equations. The interval on -axis for which the -curve becomes parallel to the -axis is recognized as the set of admissible values of for which the solution series converges. To see the range for admissible values of , , and for the present problem, -curves of , , and are shown in Figures 2, 3, and 4 for 14th-order of approximation when , , , , and . According to these figures, the convergence ranges for , , and are , , and . To assure the convergence of the HAM solution, the values of , , and should be chosen from these regions. The region for the values of , , and is strongly dependent on the values of involving parameters. Obviously our calculations show that the series (3.9) converge in the whole region of when and .

Figure 2: curve for 14th-order of approximation.

Figure 3: curve for 14th-order of approximation.

Figure 4: curve for 14th-order of approximation.

4. Results and Discussion

This section describes the graphical results of some interesting parameters for velocity and temperature profiles. Figures 5, 6, and 7 present typical tangential velocity , transverse velocity , and temperature profiles for , , , , and and various values of the magnetic parameter    and . Figures 5 and 7 demonstrate that the velocity of the fluid diminishes, whilst the temperature distribution enhances within the boundary layer the magnetic parameter rising from to . This is due to the fact that the application of a transverse magnetic field results in a drag-like force called the Lorentz force. This force tends to slow down the movement of the fluid along surface and increases in the temperature. On the other hand, as increases a crossflow in the transverse direction is greatly induced due to the Hall effect. Accordingly the transverse velocity increases as increases as shown in Figure 6.

Figure 5: Tangential velocity profiles for various values of with , , , , and Pr = 0.72.

Figure 6: Transverse velocity profiles for various values of with , , , , and Pr = 0.72.

Figure 7: Temperature profiles for various values of with , , , , and Pr = 0.72.

Figures 8, 9, and 10 illustrate the influence of the Hall parameter on the tangential velocity , transverse velocity , and temperature profiles in the boundary layer. Figure 8 shows that the tangential velocity increases while the temperature decreases with increasing . This is due to the fact that the effective conductivity () decreases with increasing which in turn reduces the magnetic damping force on . Also it is shown from these figures that the velocity and temperature profiles approach their classical hydrodynamic values when the Hall parameter increases to since the magnetic force terms approach zero value for very large values of . On the other hand Figure 9 shows that the transverse flow in the -direction first increases gradually with , reaching a maximum profile for and then decreases for being equal to zero when becomes very large. This is due to the fact that for large values of , the term () is very small and hence the resistive effect of the magnetic field is diminished.

Figure 8: Tangential velocity profiles for various values of with , , , , and Pr = 0.72.

Figure 9: Transverse velocity profiles for various values of with , , , , and Pr = 0.72.

Figure 10: Temperature profiles for various values of with , , , , and Pr = 0.72.

In Figures 11, 12, and 13, the influence of the radiation parameter on the profiles of the tangential velocity, transverse velocity, and temperature is presented, respectively, for , , , , and Pr = 0.72. From Figure 13 it is obvious that the temperature is greatly increased as is decreased. This is due to the fact that a decrease in the values of for given and leads to a decrease in the Rosseland radiation absorptivity . According to (2.6) and (2.7), the divergence of the radiative heat flux increases as decreases which in turn increases the rate of radiative heat transferred to the fluid and hence the fluid temperature increases. In view of this explanation, the effect of radiation becomes more pronounced as    and can be ignored when . The increase in the fluid temperature has a direct effect on the buoyancy force which in turn induces more flow in the boundary layer causing the velocity (tangential and transverse) of the fluid there to increase as obvious from Figures 11 and 12. Also, it is seen from Figures 11–13 that the larger the , the thinner the momentum and thermal boundary layer thickness.

Figure 11: Tangential velocity profiles for various values of with , , , , and Pr = 0.72.

Figure 12: Transverse velocity profiles for various values of with , , , , and Pr = 0.72.

Figure 13: Temperature profiles for various values of with , , , , and Pr = 0.72.

The effect of the buoyancy parameter on the tangential velocity , transverse velocity , and temperature profiles is displayed in Figures 14, 15, and 16, respectively, for , , , , and Pr = 0.72. From Figure 14 it is seen that the tangential velocity increases with the positive values (assisting) of the buoyancy parameter and decreases with the negative values (opposing flow) of as compared to the case of pure forced convection . This is due to the fact that a positive induces a favorable pressure gradient that enhances the fluid flow in the boundary layer, while a negative produces an adverse pressure gradient that slows down the fluid motion. The effect of buoyancy parameter on the transverse velocity is the same as that on the tangential velocity as shown in Figure 15. Further, it is clear from Figure 16 that the effect of buoyancy parameter is to increase the temperature in the case of opposing flow and decrease it in the case of assisting flow. The reason for this trend is that the positive (negative) buoyancy force accelerates (decelerates) the fluid in the boundary layer (as mentioned earlier) which results in thinner (thicker) thermal boundary layer.

Figure 14: Tangential velocity profiles for various values of with , , , , and Pr = 0.72.

Figure 15: Transverse velocity profiles for various values of with , , , , and Pr = 0.72.

Figure 16: Temperature profiles for various values of with , , , , and Pr = 0.72.

The effect of the wall temperature distribution characterized by the parameter on the profiles of , , and is depicted in Figures 17, 18, and 19 for , , , , and Pr = 0.72. It is seen from Figure 19 that the temperature decreases as increases. The decreased temperature has a direct effect in decreasing the thermal buoyancy forces, which in turn decrease the tangential and lateral velocities and , respectively, as obvious from Figures 17 and 18. Also, Figure 19 shows a change in the direction of heat flow as a result of the peak which occurs in the temperature profile when . The presence of the peak indicates that the temperature attains its maximum value in the body of the fluid close to the surface and not at the surface and hence heat is expected to transfer to the surface from the ambient fluid.

Figure 17: Tangential velocity profiles for various values of with , , , , and Pr = 0.72.

Figure 18: Transverse velocity profiles for various values of with , , , , and Pr = 0.72.

Figure 19: Temperature profiles for various values of with , , , , and Pr = 0.72.

The variations with positive values of the buoyancy parameter of the local skin-friction coefficient in the -direction in terms of , the local skin-friction coefficient in the -direction in terms of , and the local Nusselt number in terms of for various values of the radiation parameter and the wall temperature distribution characterized by the parameter are shown in Figures 20, 21, and 22, respectively. Figure 20 shows that for a fixed and both values of , the local skin-friction coefficient increases for negative values of ( and decreases for positive values of ( as compared to the constant wall temperature case (). Also, the effect of on is more pronounced for higher values of . In addition, for given and the local skin-friction coefficient is greatly increased as is increased. As mentioned earlier, the positive buoyancy force acts like a favorable pressure gradient which accelerates the fluid in the boundary layer. This results in thinner boundary layer and hence in higher velocity gradient at the surface. Therefore, the skin friction coefficient increases with . Further, it is noted from Figure 20 that a single value of is obtained for all and values when (the forced convection flow). This is because (2.18) and (2.20) are uncoupled when ; that is, the solutions of the flow and thermal fields are independent when buoyancy force is absent. Accordingly the parameters of the thermal field have no effect on the flow field. On the other hand, Figure 20 shows for given that decreases with increasing for all values of but the effect of is more pronounced for positive values of and larger values of .

Figure 20: Local skin-friction coefficient for various values of and for , , , and Pr = 0.72.

Figure 21: Local skin-friction coefficient for various values of and for , , , and Pr = 0.72.

Figure 22: Local Nusselt number for various values of and for , , , and Pr = 0.72.

The effects of , , and on the local skin-friction coefficient in the -direction are the same as those on as depicted in Figure 21. It is seen further that decreases with for all -values but this trend is found to be more noticeable at a positive value of and larger values of . On the other hand, for forced convection flow , the local skin-friction coefficient has a unique value of despite the values of and for the same reason mentioned previously.

Figure 22 reveals that the local Nusselt number is increased with increasing the value of for all values of and ; in other words, heat transfer rate can be enhanced by enlarging the surface temperature variation. Also, the local Nusselt number increases with for all and a negative value of . On the other hand, the effect of radiation parameter on for a constant wall temperature () and a positive value of is opposite to that of negative value of . Namely, the heat transfer rate is reduced as is increased for a positive value of and all values of . It is noted that negative heat transfer rates are obtained for which indicate that heat is transferred from the fluid to the stretching surface in spite of the excess of surface temperature over that of the free stream fluid. On the other hand, increases with for values of while the opposite trend is true for .

The variations of the local skin-friction coefficient in the -direction , the local skin-friction coefficient in the -direction , and the local Nusselt number as a function of the buoyancy parameter for two values of and and different values of Hall parameter are presented in Figures 23, 24, and 25. It is clear from these figures that for fixed and the skin-friction coefficient in the -direction and the local Nusselt number decrease while the local skin-friction coefficient in the -direction increases with increasing the magnetic parameter . Further, it is remarkable to note that the effect of the magnetic parameter on and for higher values of ( and ) is much less than that for lower values of () while the opposite trend is noticed for the local skin-friction coefficient in the -direction . On the other hand, for all and values Figures 23 and 25 show that the local skin-friction in the -direction and the local Nusselt number increase as increases. Further, from Figure 24 it is interesting to note that for all values of and the local skin friction coefficient in the -direction increases as increases from to and then decreases as increases from to and this result agrees well with the profiles of introduced in Figure 9.

Figure 23: Local skin-friction coefficient for various values of and for , , , and Pr = 0.72.

Figure 24: Local skin-friction coefficient for various values of and for , , and Pr = 0.72.

Figure 25: Local Nusselt number for various values of and for , , , and Pr = 0.72.

Finally, Table 1 shows a good agreement of our results for heat transfer rate obtained by homotopy analysis method with the numerical results reported by Magyari and Keller [25] for forced convective flow () of a viscous incompressible fluid over an exponentially stretching sheet in the absence of magnetic field , Hall current , and radiation at selected values of and .

Table 1: Comparison of the values of obtained by HAM for and and various values of Pr and with those of Magyari and Keller [25].

Pr

Reference [25]
HAM
Reference [25]
HAM
Reference [25]
HAM
Reference [25]
HAM

0.5
0.175815
0.174455
0.330493
0.330824
0.594338
0.594112
1.008405
1.007963

1
0.299876
0.299795
0.549643
0.549589
0.954782
0.955325
1.560294
1.560759

3
0.634113
0.635700
1.122188
1.123344
1.869075
1.869397
2.938535
2.934606

5
0.870431
0.874551
1.521243
1.524815
2.500135
2.500382
3.886555
3.889363

10
1.308613
1.309138
2.257429
2.256776
3.660379
3.673949
5.628198
5.625053

Acknowledgments

This research is supported by King Khalid University (KKU), Saudi Arabia, under the project no. KKU-SCI-11-10. This support is highly appreciated and acknowledged. The authors are very thankful to the reviewers for their encouraging comments and constructive suggestions to improve the presentation of this work.

References

B. C. Sakiadis, “Boundary layer behavior on continuous solid surface, the boundary layer on a continuous flat surface,” AIChE Journal, vol. 7, pp. 221–225, 1961.

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

F. K. Tsou, E. M. Sparrow, and R. J. Goldstein, “Flow and heat transfer in the boundary layer on a continuous moving surface,” International Journal of Heat and Mass Transfer, vol. 10, no. 2, pp. 219–235, 1967.

L. J. Grubka and K. M. Bobba, “Heat transfer characteristics of a continuous, stretching surface with variable temperature,” Journal of Heat Transfer, vol. 107, no. 1, pp. 248–250, 1985.

L. E. Erickson, L. T. Fan, and V. G. Fox, “Heat and mass transfer on a moving continuous flat plate with suction or injection,” Industrial and Engineering Chemistry Fundamentals, vol. 5, no. 1, pp. 19–25, 1966.

E. M. Abo-Eldahab and M. A. El Aziz, “Blowing/suction effect on hydromagnetic heat transfer by mixed convection from an inclined continuously stretching surface with internal heat generation/absorption,” International Journal of Thermal Sciences, vol. 43, no. 7, pp. 709–719, 2004.

R. Cortell, “Flow and heat transfer of a fluid through a porous medium over a stretching surface with internal heat generation/absorption and suction/blowing,” Fluid Dynamics Research, vol. 37, no. 4, pp. 231–245, 2005.

I. C. Liu, “Flow and heat transfer of an electrically conducting fluid of second grade in a porous medium over a stretching sheet subject to a transverse magnetic field,” International Journal of Non-Linear Mechanics, vol. 40, no. 4, pp. 465–474, 2005.

A. Ali and A. Mehmood, “Homotopy analysis of unsteady boundary layer flow adjacent to permeable stretching surface in a porous medium,” Communications in Nonlinear Science and Numerical Simulation, vol. 13, no. 2, pp. 340–349, 2008.

T. S. Chen and F. A. Strobel, “Buoyancy effects in boundary layer adjacent to a continuous, moving horizontal flat plate,” Journal of Heat Transfer, vol. 102, no. 1, pp. 170–172, 1980.

J. R. Fan, J. M. Shi, and X. Z. Xu, “Similarity solution of mixed convection over a horizontal moving plate,” Heat and Mass Transfer, vol. 32, no. 3, pp. 199–206, 1997.

C. H. Chen, “Laminar mixed convection adjacent to vertical, continuously stretching sheets,” Heat and Mass Transfer, vol. 33, no. 5-6, pp. 471–476, 1998.

M. Ali and F. Al-Yousef, “Laminar mixed convection from a continuously moving vertical surface with suction or injection,” Heat and Mass Transfer, vol. 33, no. 4, pp. 301–306, 1998.

M. Ali and F. Al-Yousef, “Laminar mixed convection boundary layers induced by a linearly stretching permeable surface,” International Journal of Heat and Mass Transfer, vol. 45, no. 21, pp. 4241–4250, 2002.

M. A. El-Aziz and A. M. Salem, “MHD-mixed convection and mass transfer from a vertical stretching sheet with diffusion of chemically reactive species and space or temperature dependent heat source,” Canadian Journal of Physics, vol. 85, no. 4, pp. 359–373, 2007.

A. Moutsoglou and T. S. Chen, “Buoyancy effects in boundary layers on inclined, continuous, moving sheets,” Journal of Heat Transfer, vol. 102, no. 2, pp. 371–373, 1980.

F. A. Strobel and T. S. Chen, “Buoyancy effects on heat and mass transfer in boundary layers adjacent to inclined, continuous, moving sheets,” Numerical Heat Transfer, vol. 3, no. 4, pp. 461–481, 1980.

C. H. Chen, “Mixed convection cooling of a heated, continuously stretching surface,” Heat Mass Transfer, vol. 36, no. 1, pp. 79–86, 2000.

W. G. England and A. F. Emery, “Thermal radiation effects on laminar free convection boundary layer of an absorbing gas,” Journal of Heat Transfer, vol. 31, pp. 37–44, 1969.

R. S. R. Gorla and I. Pop, “Conjugate heat transfer with radiation from a vertical circular pin in a non-newtonian ambient medium,” Wärme- und Stoffübertragung, vol. 28, no. 1-2, pp. 11–15, 1993.

A. Raptis, “Radiation and free convection flow through a porous medium,” International Communications in Heat and Mass Transfer, vol. 25, no. 2, pp. 289–295, 1998.

M. Abd-El Aziz, “Thermal radiation effects on magnetohydrodynamic mixed convection flow of a micropolar fluid past a continuously moving semi-infinite plate for high temperature differences,” Acta Mechanica, vol. 187, no. 1–4, pp. 113–127, 2006.

M. Abd El-Aziz, “Thermal-diffusion and diffusion-thermo effects on combined heat and mass transfer by hydromagnetic three-dimensional free convection over a permeable stretching surface with radiation,” Physics Letters A, vol. 372, no. 3, pp. 263–272, 2008.

M. Abd El-Aziz, “Radiation effect on the flow and heat transfer over an unsteady stretching sheet,” International Communications in Heat and Mass Transfer, vol. 36, no. 5, pp. 521–524, 2009.

E. Magyari and B. Keller, “Heat and mass transfer in the boundary layers on an exponentially stretching continuous surface,” Journal of Physics D, vol. 32, no. 5, pp. 577–585, 1999.

E. M. A. Elbashbeshy, “Heat transfer over an exponentially stretching continuous surface with suction,” Archives of Mechanics, vol. 53, no. 6, pp. 643–651, 2001.

S. K. Khan and E. Sanjayanand, “Viscoelastic boundary layer flow and heat transfer over an exponential stretching sheet,” International Journal of Heat and Mass Transfer, vol. 48, no. 8, pp. 1534–1542, 2005.

E. Sanjayanand and S. K. Khan, “On heat and mass transfer in a viscoelastic boundary layer flow over an exponentially stretching sheet,” International Journal of Thermal Sciences, vol. 45, no. 8, pp. 819–828, 2006.

M. K. Partha, P. V. S. N. Murthy, and G. P. Rajasekhar, “Effect of viscous dissipation on the mixed convection heat transfer from an exponentially stretching surface,” Heat and Mass Transfer, vol. 41, no. 4, pp. 360–366, 2005.

M. A. El-Aziz, “Viscous dissipation effect on mixed convection flow of a micropolar fluid over an exponentially stretching sheet,” Canadian Journal of Physics, vol. 87, no. 4, pp. 359–368, 2009.

G. W. Sutton and A. Sherman, Engineering Magnetohydrodynamics, McGraw-Hill, New York, NY, USA, 1965.

I. Pop and T. Watanabe, “Hall effects on magnetohydrodynamic free convection about a semi-infinite vertical flat plate,” International Journal of Engineering Science, vol. 32, no. 12, pp. 1903–1911, 1994.

E. M. Abo-Eldahab and M. A. El Aziz, “Hall and ion-slip effects on MHD free convective heat generating flow past a semi-infinite vertical flat plate,” Physica Scripta, vol. 61, no. 3, pp. 344–348, 2000.

A. A. Megahed, S. R. Komy, and A. A. Afify, “Similarity analysis in magnetohydrodynamics: hall effects on free convection flow and mass transfer past a semi-infinite vertical flat plate,” International Journal of Non-Linear Mechanics, vol. 38, no. 4, pp. 513–520, 2003.

E. M. Abo-Eldahab and M. A. El Aziz, “Hall current and ohmic heating effects on mixed convection boundary layer flow of a micropolar fluid from a rotating cone with power-law variation in surface temperature,” International Communications in Heat and Mass Transfer, vol. 31, no. 5, pp. 751–762, 2004.

E. M. Abo-Eldahab and M. A. El Aziz, “Viscous dissipation and Joule heating effects on MHD-free convection from a vertical plate with power-law variation in surface temperature in the presence of Hall and ion-slip currents,” Applied Mathematical Modelling, vol. 29, no. 6, pp. 579–595, 2005.

L. K. Saha, M. A. Hossain, and R. S. R. Gorla, “Effect of Hall current on the MHD laminar natural convection flow from a vertical permeable flat plate with uniform surface temperature,” International Journal of Thermal Sciences, vol. 46, no. 8, pp. 790–801, 2007.

E. M. Abo-Eldahab, M. A. El-Aziz, A. M. Salem, and K. K. Jaber, “Hall current effect on MHD mixed convection flow from an inclined continuously stretching surface with blowing/suction and internal heat generation/absorption,” Applied Mathematical Modelling, vol. 31, no. 9, pp. 1829–1846, 2007.

A. M. Salem and M. Abd El-Aziz, “Effect of Hall currents and chemical reaction on hydromagnetic flow of a stretching vertical surface with internal heat generation/absorption,” Applied Mathematical Modelling, vol. 32, no. 7, pp. 1236–1254, 2008.

M. Abd El-Aziz, “Flow and heat transfer over an unsteady stretching surface with Hall effect,” Meccanica, vol. 45, no. 1, pp. 97–109, 2010.

M. Sajid and T. Hayat, “Influence of thermal radiation on the boundary layer flow due to an exponentially stretching sheet,” International Communications in Heat and Mass Transfer, vol. 35, no. 3, pp. 347–356, 2008.

S. J. Liao, The Proposed Homotopy Analysis Technique for the Solution of Nonlinear Problems [Ph.D. thesis], Shanghai Jiao Tong University, 1992.

S. Liao, Beyond Perturbation: Introduction to the Homotopy Analysis Method, Chapman & Hall, Boca Raton, Fla, USA, 2003.

P. J. Hilton, An Introduction to Homotopy Theory, Cambridge University Press, 1953.

T. F. C. Chan and H. B. Keller, “Arc-length continuation and multigrid techniques for nonlinear elliptic eigenvalue problems,” SIAM Journal on Scientific and Statistical Computing, vol. 3, no. 2, pp. 173–194, 1982.

E. E. Grigolyuk and V. I. Shalashilin, The Continuation Method Applied to Nonlinear Problems in Solid Mechanics, Kluwer, 1991.

G. W. Sutton and A. Sherman, Engineering Magnetohydrodynamics, McGraw-Hill, New York, NY, USA, 1965.

T. G. Cowling, Magnetohydrodynamics, Interscience, New York, NY, USA, 1957.

M. Kumari and G. Nath, “Unsteady MHD mixed convection flow over an impulsively stretched permeable vertical surface in a quiescent fluid,” International Journal of Non-Linear Mechanics, vol. 45, no. 3, pp. 310–319, 2010.

A. R. Ghotbi, H. Bararnia, G. Domairry, and A. Barari, “Investigation of a powerful analytical method into natural convection boundary layer flow,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 5, pp. 2222–2228, 2009.

Z. Ziabakhsh, G. Domairry, M. Mozaffari, and M. Mahbobifar, “Analytical solution of heat transfer over an unsteady stretching permeable surface with prescribed wall temperature,” Journal of the Taiwan Institute of Chemical Engineers, vol. 41, no. 2, pp. 169–177, 2010.

S. Liao, “An optimal homotopy-analysis approach for strongly nonlinear differential equations,” Communications in Nonlinear Science and Numerical Simulation, vol. 15, no. 8, pp. 2003–2016, 2010.