Alexandria Engineering Journal (2016) xxx, xxx-xxx

HOSTED BY

Alexandria University Alexandria Engineering Journal

www.elsevier.com/locate/aej www.sciencedirect.com

ORIGINAL ARTICLE

Effect of partial slip on an unsteady MHD mixed convection stagnation-point flow of a micropolar fluid towards a permeable shrinking sheet

Aurangzaib a,% K. Bhattacharyyab, S. Shafiec

a Department of Mathematical Sciences, Federal Urdu University of Arts, Science & Technology, Gulshan-e-Iqbal, Karachi, Pakistan b Department of Mathematics, Institute of Science, Banaras Hindu University, Varanasi 221005, Uttar Pradesh, India c Department of Mathematical Sciences, Faculty of Science, Universiti Teknologi Malaysia JB, 81310 Skudai, Johor, Malaysia

Received 22 December 2014; revised 20 March 2016; accepted 7 April 2016

KEYWORDS

Unsteady; Partial slip effect; MHD;

Micropolar fluid; Permeable shrinking sheet

Abstract The objective of the present study was to investigate the partial slip effect on an unsteady two-dimensional mixed convection stagnation point flow towards a permeable shrinking sheet. The governing equations are reduced to a system of non-dimensional partial differential equations using a semi-similarity transformation, before being solved numerically by using Keller-box method. The features of the flow characteristics for different values of the governing parameters are analysed and discussed. The results indicate that the momentum, thermal and concentration boundary layer thicknesses increase with increasing mixed convection parameter for opposing flow, whereas the opposite effect is observed for assisting flow. The results also show that the surface velocity is higher when there is slip at a sheet compared to its absent. Further, the study indicates that the boundary layer thicknesses become thicker and thicker with increasing shrinking parameter, while the opposite effect is observed with increasing Hartmann number. Comparison with previously published work for special cases is performed and found to be in excellent agreement.

© 2016 Faculty of Engineering, Alexandria University. Production and hosting by Elsevier B.V. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/).

1. Introduction

The study of hydromagnetic stagnation point flow of an incompressible viscous fluid has attracted many researchers due to wide applications in industries and engineering. Some of the applications are the aerodynamic extrusion of plastic sheets, geothermal energy extractions, cooling of nuclear reac-

* Corresponding author.

E-mail address: zaib20042002@yahoo.com ( Aurangzaib).

Peer review under responsibility of Faculty of Engineering, Alexandria

University.

tor, plasma studies, cooling of underground electric cables, geothermal energy extractions, artificial fibres and blood flow problems. Hiemenz [1] was first to investigate the two-dimensional stagnation flow towards a stationary plate. On the other hand, the flow over a linearly stretching plate was first considered by Crane [2]. Chiam [3] combined Hiemenz and Crane's problems by considering the steady two-dimensional stagnation-point flow of an incompressible viscous fluid towards a stretching sheet. Kumari and Nath [4] studied the effect of the magnetic field on the stagnation point flow with heat transfer over a stretching sheet. Mahapatra and Gupta [5] investigated the steady two-dimensional

http://dx.doi.org/10.1016/j.aej.2016.04.018

1110-0168 © 2016 Faculty of Engineering, Alexandria University. Production and hosting by Elsevier B.V. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/).

stagnation-point flow of an incompressible viscous electrically conducting fluid towards a stretching sheet. Sharma and Singh [6] studied the steady MHD two-dimensional stagnation point flow towards a stretching sheet with variable thermal conductivity. Ashraf and Rashid [7] investigated the radiation effect on MHD two-dimensional stagnation point flow towards a heated shrinking sheet. Mostafa and Shimaa [8] investigated the combined effects of magnetic field and thermal radiation on a steady micropolar fluid near a stagnation point towards a moving surface. Bhattacharyya et al. [9] studied the boundary layer flow near a stagnation point with heat and mass transfer over a shrinking sheet with Soret and Dufour effects. Zeeshan et al. [10] investigated MHD natural convection flow of a water/ethylene glycol nanofluid with inverted cone through a porous medium. Ellahi et al. [11] examined MHD flow of a Jeffrey fluid on a peristaltic flow in a rectangular duct embedded in a porous medium. Khan et al. [12] investigated the peristaltic motion and heat transfer of Oldroyd fluid in a channel with inclined magnetic field. The effect of magnetic field on blood flow of Prandtl fluid in a porous walls through a tapered stenosed artery was investigated by Ellahi et al. [13]. Akbar et al. [14] studied the combined effects of magnetic field and heat flux for two different water based nanoparticles for the peristaltic flow in a symmetric vertical permeable channel. Recently, Kandelousi and Ellahi [15] examined the combined effects of MHD and ferrohydrodynamic on Fe3O4-plasma nanofluid flow in a vessel and obtained the results using lattice Boltzmann method.

In recent times, the boundary layer flow due to a shrinking sheet is received considerable attention due to its wide engineering applications. Miklaviciic and Wang [16] studied the steady viscous flow over a shrinking sheet for both two-dimensional and axisymmetric cases. The effect of suction on MHD boundary layer flow with heat and mass transfer past a shrinking sheet was examined by Muhaimin et al. [17]. Fang and Zhang [18] obtained an exact solution for steady MHD flow with suction over a shrinking sheet. Since the vorticity (rotation or non-potential) flow over a shrinking sheet is not confined within the boundary layer and the flow is likely to exist, whether an adequate suction on the boundary is imposed or a stagnation point (which contains the vorticity) is considered. The stagnation point flow towards a shrinking sheet for both two-dimensional and axisymmetric cases was first investigated by Wang [19]. He found that solutions do not exist for large shrinking rates and may be non-unique in two-dimensional case. Fang et al. [20] numerically investigated the viscous flow over an unsteady shrinking sheet with mass transfer. Fan et al. [21] extended Wang's problem to unsteady flow. Suali et al. [22] studied heat transfer on an unsteady stagnation point flow with suction over a stretching/shrinking sheet. Recently, Bhattacharyya [23] discussed the unsteady boundary layer flow and heat transfer near a stagnationpoint towards a shrinking/stretching sheet.

All of the above mentioned studies are restricted to the flows of Newtonian fluids. These Newtonian fluids cannot describe some engineering and industrial processes which are made up of materials having an internal structure. The micropolar fluids are described as a non-Newtonian fluid consisting of dumb-bell molecule, fluid suspension, polymer fluids, animal blood. The presence of smoke or dust in a gas may also model using micropolar fluid dynamics. The theory of micropolar fluid model introduced by Eringen [24] exhibits the local effects

arising from the microstructure and micro motion of the fluid elements. Comprehensive reviews of the subject and its applications can be found in the review articles of Peddieson and McNitt [25], Ariman et al. [26,27] and books by tukaszewicz [28] and Eringen [29]. Guram and Smith [30] and Gorla [31] investigated the steady micropolar fluid near a stagnation point. Nazar et al. [32] examined the steady two-dimensional stagnation-point flow over a stretching sheet. Ishak et al. [33] studied the two-dimensional mixed convection boundary layer flow of a micropolar fluid near the stagnation point over a stretching sheet. Ishak et al. [34] numerically investigated the steady two-dimensional stagnation-point flow of a micropolar fluid over a shrinking sheet. Ashraf and Bashir [35] numerically investigated the effect of MHD on the steady two-dimensional stagnation-point flow of a micropolar fluid towards a heated shrinking sheet. Aman et al. [36] studied the mixed convection stagnation point flow on a vertical surface in the presence of slip effects. Reddy et al. [37] considered the boundary layer flow near a stagnation point of electrically conducting micropolar fluid over a stretching sheet with viscous dissipation, chemical reaction and heat source/sink. Hayat et al. [38] examined the unsteady stagnation point flow of second grade fluid over a porous stretching sheet. Ellahi et al. [39] investigated the unsteady flow of a micropolar fluid through a composite stenosis in arterial blood flow with slip velocity. In other paper, Ellahi et al. [40] studied micropolar fluid with heat and mass transfer of blood flow in a porous walls through a tapered stenosis arteries. Recently, Aurangzaib and Shafie [41] investigated the unsteady MHD flow near a stagnation-point of a micropolar fluid past a shrinking sheet with thermophoresis and slip effects.

In some practical situations, it is important to replace the no-slip condition by the partial slip condition because when fluid flows in micro electro mechanical system (MEMS), the no-slip boundary condition at the solid-fluid interface is no longer applicable. The non-equilibrium region near the interface is more accurately described by the slip flow model. Wang [42] found closed form solution of two-dimensional viscous flow due to stretching sheet with the effects of slip and suction. Aman and Ishak [43] investigated the effect of partial slip on the steady incompressible flow over shrinking permeable sheet. Das [44] studied the effect of partial slip and temperature dependent fluid properties on steady hydro-magnetic micropolar fluid flow over an inclined permeable plate in the presence of chemical reaction, thermal radiation and heat source/sink. Das [45] investigated the steady MHD mixed convection boundary layer stagnation point flow of a micropolar fluid towards a shrinking sheet in the presence of partial slip. Ellahi and Hameed [46] investigated MHD flow of a second grade fluid with heat transfer in a channel with nonlinear slip effects. Zeeshan and Ellahi [47] obtained the series solution of MHD slip flow and heat transfer of a non-Newtonian third grade fluid in a porous pipe. The mixed convection flow and heat transfer past a porous vertical cylinder were investigated by Ellahi et al. [48]. Recently, Chen [49] considered the unsteady mixed convection stagnation point flow on a stretching sheet with slip effect.

Motivated by the above discussions, the effect of partial slip on an unsteady MHD mixed convection stagnation point flow with heat and mass transfer in a micropolar fluid towards a permeable shrinking sheet is investigated in this paper. The unsteadiness is caused by the impulsive motion of the external

stream or of the shrinking surface. The governing equations are solved numerically by using Keller-box method. The obtained results are then compared with previous published results available in the literature.

2. Mathematical formulations

Consider an unsteady two-dimensional mixed convection stagnation point flow of an incompressible electrically conducting micropolar fluid towards a permeable shrinking sheet. A transverse magnetic field B0 is applied in the y-direction and induced magnetic field is negligibly assumed to be very small magnetic Reynolds number. The flow configuration is shown schematically in Fig. 1. It is assumed that the velocity uw(x), the temperature Tw(x) and the concentration Cw(x) of the sheet are proportional to the distance x from the stagnation point, where Tw(x) > Tx and Cw(x) > Cx with Tx and Cx being the uniform temperature and concentration of the ambient fluid. The velocity of the flow external to the boundary layer is ue(x). Under the usual boundary layer approximation, along with Boussinesq approximations the governing equations are governed by the following equations: du dv 0 dx dy '

du du du due f i + k\ d2u k dN d t dx dy e dx \ p J dy2 p dy + gb(T— Ti)+gb *(C — Ci) rB,

H--(Ue - u),

dN dN dN y d2N k

— + u — + v — = — --: 2N + —

dx dy pj dy

dT dT dT_ k d2T dt dx dy pcp dy2 '

— +u —+v @£-D —

dt dx dy dy2

Assisting flow

Opposing flow

where (u, v) are the velocity components along the (x, y) axes, respectively, cp is the specific heat at constant pressure, k is the vortex viscosity, r is the electrical conductivity of the fluid, B0 is externally imposed magnetic field strength in the y-direction, p is the density of the fluid, D is the molecular diffusivity, k is the coefficient of thermal expansion, N is the components of microrotation or angular velocity whose rotation is in the direction of the x, y - planes, y is the spin gradient viscosity given by y —(i + K=2)j — i(1 + K/2)j, j is the micro-inertia density, and i is the coefficient of viscosity of the fluid. According to Xu [50] we assume that j is constant therefore j — t/a.

The appropriate physical boundary conditions for the problem under study are given by

t < 0 : u(t, x,y) — 0, v(t, x,y) — 0, N(t, x,y) — 0, T(t, x, y)

— Ti, C(t, x,y) — Ci for any x,y,

t P 0 : u — bx + L —, v — v0, N ——n —, T — Tw, C dy dy

— Cw at y — 0,

u ue x

(x), N ! 0, T ! Ti, C ! Ci as y-

The sheet velocity uw(x), temperature Tw(x) and concentration Cw(x) are assumed to vary proportional to the distance x from the stagnation point and the velocity of the flow external to the boundary layer is ue(x) are taken of the form:

uw(x) — bx, Tw(x) — Ti + cx, Cw(x)

— Ctx + dx and ue(x) — ax. (7)

where a; b; c and d are positive constants, and n is a constant such that 0 6 n 6 1. It is noted the case when n — 0 is called a strong concentration. In this case, N — 0 near the wall implying that the concentrated particle flows in which the microelements close to the wall surface are unable to rotate [51]. The case n — 1 /2 presents the vanishing of the antisymmetric part of the stress tensor and denotes weak concentrations whereas n — 1 is used for modelling of turbulent boundary layer flows. In this study, we only consider strong concentration (n — 0).

Introducing the following non-dimensional variables as

w — paffit nxf(n, g), g — p y, N x h(n, g) n — VI - s — at, 0(£, g)

_ (T-Ti)

( Tw T1 )

>/(n g)—ë-cù ■

where 0 < n 6 1 and the continuity Eq. (1) is satisfied by introducing a stream function W such that

u — IT, v — .

The transformed nonlinear partial differential equations are as follows:

(1 + f + 2g(1 — nX + Kh, + n2 (1 + ffm — f + M2(1 — f) + 1(9 + S/j)

— 2n(1 - n2) f,

Figure 1 Physical model and coordinate system.

(1 + K/2)hg„ + 2(1 - n2)(h + ghg) + n2(Jhg - fgh - K(2h + fj)

- 2^(1 - n2) g,

Pdgg + 2g(1 - n2)hg + n2(fhg -fgh) - 2n(1 - n2) ;

/„ + 2g(1 - n2)/g + n2(//g -fg/) - 2n(1 - n2) ,

subject to the boundary conditions (6) which become

M, 0)-f0 ,fg(n, 0)-e + /„(n, 0),

h(n, 0) - -nfgg(n, 0), h(n, 0) - 1,/(n, 0) - 1,

fg(n, g) ! 1, h(n, g) ! 0, h(n, g) ! 0,/(n, g) ! 0 as g

(10) (11) (12)

where f0 --Vw/^/av is the suction (f0 > 0) or injection (f0 < 0) parameter, M2 - B0r/pa is the Hartmann number, k - Grx/Re\ is the mixed convection parameter, Grx - [gb(Tw - Tœ)x3/D2] is the local Grashof number, Rex - uwx/t is the local Reynolds number, d - [fi'(C„ - C1)/b(Tw - Tœ)] is the constant concentration buoyancy parameter, Pr - i-p/k is the Prandtl number, a - L^/a/t is the slip parameter, Sc - t/D is Schmidt number, and e - b/a with e > 0 and e < 0 correspond to a stretching and shrinking sheet, respectively. It should be noted that dual solutions is only possible when (n - 1) which can be seen in Figs. 17-20.

The quantities of physical interest are the skin-friction coefficient in the x-direction C/x, and the Nusselt number Nux and the Sherwood number Shx are defined respectively by,

Cfx T ; sw

(l + k) — + kN dy

then CfxRe

- n [1 + (1 - n)Kfgg(n, 0),

1w --k(@^) , Nux

VdWy=0 (Tw -

(k) then Nux/Re\

-- hg(n, 0),

-pD(Sh - mw

VdWy=0' x (Cw - Ci) \.PD

then Shx/Re\

-j/g (n, 0).

Figure 2 Variation of k on the velocity f

3. Numerical method

The nonlinear partial differential Eqs. (9)-( 1 2) subject to boundary conditions ( 3) were solved numerically using an implicit finite difference scheme, known as the Keller-box method. This method has been found to be very suitable in dealing with nonlinear parabolic problems. By introducing the new dependent variables, these equations are first written as a system of first-order equations which are then expressed in finite difference forms using central differences. Since the system of equations is nonlinear, it is linearized using Newton's method before putting them in matrix-vector form. The resulting linear system is solved along with the boundary conditions by the block-tridiagonal-elimination method. Details of the method can be found in the book by Cebeci and Bradshaw [52]. In this study, we used MATLAB software. Here, the value of step size Dg = 0.01 has been used and the convergence criterion was set to 5 x 10—5, which give accuracy to four decimal places. The satisfaction of the outer boundary condition is achieved by considering the boundary layer thickness = 10. However for dual solutions, we considered different boundary layer thicknesses . The correctness of the present numerical scheme, a comparison of the obtained numerical results corresponding to the skin friction coefficient /"(0) is checked with the previous published results of Wang [19] and Ishak et al. [34] in Table 1. Thus, it is seen from Table 1 that the numerical results are in close agreement.

4. Results and discussion

Figs. 2-5 illustrate the effects for the mixed convection parameter k on the velocity, microrotation, temperature and concentration of the fluid for both assisting and opposing flows. From Figs. 2 and 3, it is observed that the velocity and microrotation profiles decrease for opposing flow and it makes the momentum and microrotation boundary layer thicknesses thicker and thicker, while the reverse trend is observed for the assisting flow. Negative values of velocity close to the sheet indicate the region of reverse flow. This is due to the fact that an increase in the value of k leads to an increase in the temperature difference (Tw — T1). This leads to an increase in the convection currents as a result of which fluid velocity increases and thus, boundary layer thickness decreases in case of assisting flow. It is also observed from Fig. 3 that the microrotation profile is negative near the sheet and it gradually increases, i.e. accelerates the fluid far away from the sheet. The temperature of the fluid decreases with an increase in the value of mixed convection parameter until it achieves a constant value, namely zero for the assisting flow as shown in Fig. 4, whereas the opposite trend is observed for the opposing flow. This is not surprising since the fluid receives the heat from the surface and then the heat energy changed into other forms such as kinetic energy. The concentration boundary layer thickness decreases with increasing mixed convection parameter for the assisting flow as shown in Fig. 5, while the reverse trend is observed for

Table 1 Values of f(0) for stretching sheet in the absence of buoyancy force Hartmann number M2 and slip parameter a.

e Wang [19] Ishak et al. [34] Present (n = 1)

K = 0 K = 1, n = 0.5 K = 0 K = 1, n = 0.5

0 1.23258 1.232588 1.006404 1.232588 1.006405

0.1 1.14656 1.146561 0.936163 1.146562 0.936162

0.2 1.05113 1.051130 0.858244 1.051131 0.858246

0.5 0.71330 0.713295 0.582403 0.713300 0.582403

1 0 0 0 0 0

2 -1.88731 -1.887307 -1.540979 -1.887307 -1.540980

3 -4.276541 -3.491781 -4.276540 -3.491782

4 -7.086378 -5.786003 -7.086380 -5.786005

5 -10.26475 -10.264749 -8.381133 -10.264750 -8.381132

Figure 3 Variation of k on the microrotation h.

0.8 0.6 0.4 0.2 0,

K = M =1, Pr= 0.71, Sc = 0.94, 6=1,J =1, f=-fl.7J, a = 0.2, f =

V \ Vs " V^ Vs \ «A N N \ mV \ \ v ».=1,3,5 -Assisting flow ---Opposing flow -

* = 0,1,3,5

Figure 5 Variation of k on the concentration /.

Figure 4 Variation of k on the temperature h. Figure 6 Variation of e on the velocity f

the case of opposing flow. It is evident from these figures that all curves approach the far field boundary conditions (13) asymptotically.

Figs. 6-9 illustrate the effect of shrinking parameter e on the velocity, microrotation, temperature and concentration profiles. It is observed from Figs. 6 and 7 that the velocity and microrotation profiles decrease with an increase in the value of magnitude of shrinking parameter e and consequently the thicknesses of momentum and microrotation boundary layer increases. Reverse flow is observed near the surface as the magnitude of shrinking parameter increases. On the other hand, the temperature and concentration of fluid increase with an increase in the value of magnitude of shrinking parameter e

as shown in Figs. 8 and 9. Thus, the thermal and concentration boundary layer thicknesses become thicker and thicker.

The influence of Hartmann number M2 on the velocity, microrotation, temperature and concentration profiles is displayed in Figs. 10-13 in the presence of slip and in the absence of slip at the boundary. From Figs. 10 and 11, it is observed that for both slip and no-slip cases, the velocity and microrotation profiles along the sheet increase with increasing Hartmann number and it makes the momentum and microrotation boundary layer thicknesses thinner and thinner. This is due to the fact that the magnetic force enhances the fluid motion in the boundary layer because of the last term of momentum equation (a^ — u) remains positive in the boundary layer region. Here an electrically conducting fluid

Figure 11 Variation of M2 and a on the microrotation h.

Figure 12 Variation of M2 and a on the temperature h.

Figure 7 Variation of e on the microrotation h.

Figure 8 Variation of e on the temperature h.

K = 1, Pr= 0.71, Sc = 0

I, 1 = 1, 8=2, e=-0.75J =1, (=0

Figure 9 Variation of e on the concentration /.

-06 = 0.0

m % m m. m ---a =0.5 "

, M2 = 0,1,5, 7

Figure 10 Variation of M2 and a on the velocity f

Figure 13 Variation of M2 and a on the concentration /.

produces a drag-like force called the Lorentz force associated with the magnetic field makes the boundary layer thinner. As expected, the velocity and microrotation profiles are lower for the case of no-slip than those for the presence of slip. Figs. 12 and 13 show that the temperature and concentration of fluid decrease with increasing Hartmann number for slip as well as for no-slip conditions and cause to decrease the thermal and concentration boundary layer thicknesses. On the other hand, in the presence of slip, the temperature and concentration of fluid are lower compared to the case of no-slip. This means that the heat and mass are transferred from hot surface to the cool fluid.

Figs. 14-16 demonstrate the effect of material parameter K on the reduced skin friction, the reduced Nusselt number and the reduced Sherwood number for suction or injection. It is observed from these results that the reduced skin friction,

Figure 16 Variation of K and f on the reduced Sherwood number versus k.

reduced Nusselt number and reduced Sherwood number decrease with an increase in the value of micropolar parameter K. This discloses the fact that the micropolar fluids offer great resistance to the fluid motion as compared with Newtonian fluid. This phenomenon also reflects the fact that the increase in the value of micropolar parameter results in an enhancement of the total viscosity in the fluid flow. Therefore, micropolar fluid is a very efficient fluid medium in boundary layer for controlling the laminar flow. It is also observed from these figures that the suction parameter enhances the reduced skin friction, the reduced Nusselt number and the reduced Sherwood number as shown in Table 2. Opposite effect is observed for

Table 2 Computed values of the reduced skin friction,

reduced Nusselt number and reduced I Sherwood number for

various values f0, k and K when M2 = 1, d = 1, Sc = 0.94,

Pr = 0.71, e = -0.75 a = 0.2, n = 0.8.

fo k K f„(i, 0) -№, 0) -/g(n, 0)

-1 0 0 1.469512 0.532935 0.562585

1 1.133350 0.489613 0.514304

3 0.816596 0.446069 0.466404

2 0 1.799685 0.559284 0.592588

1 1.362586 0.509448 0.536666

3 0.968109 0.460173 0.482152

1 0 0 2.264289 1.172826 1.397221

1 1.616277 1.105574 1.318972

3 1.097223 1.040262 1.244493

2 0 2.555103 1.195132 1.423502

1 1.810903 1.122571 1.338763

3 1.221325 1.052368 1.258434

fn -0.5

-1 -1.5 -2

M1 = l,Pr=l, Sc = 1, X = -0.1, S= -0.1, f = 2.4, * = -1.1S, a = 0.2

II/ K= 0, 0.5, 1

^ / ^ / -First solution

— Second solution

1 2 3 4 5

Figure 17 Variation of K on the velocity f.

Figure 18 Variation of K on the microrotation h.

injection. Further, it is seen from these results that the value of/(£,0), 0) and -/'(n,0) is always positive, which implies that the fluid exerts a drag force on the sheet and heat and mass are transferred from hot sheet to cold fluid, respectively.

Figs. 17-20 show the dual velocity, microrotation, temperature and concentrations for different values of K. Examining the dual velocity profiles (Fig. 17) for different values of K, it

Figure 19 Variation of K on the temperature h.

Figure 20 Variation of K on the concentration /.

reveals that due to the material parameter, the momentum boundary layer thickness increases for the first solution and reduces for the second solution. On the other hand, the microrotation profiles decrease with increasing K for both first and second solutions as shown in Fig. 18. From Figs. 19 and 20, it is found that the temperature and concentration of fluid increase with K for both solutions. Further, the second solution exhibits a larger boundary layer thickness compared to the first solution.

5. Conclusions

The combined effects of heat and mass transfer on an unsteady MHD mixed stagnation point flow of a micropolar fluid towards a permeable shrinking sheet in the presence of slip effect are investigated. The numerical results were compared with those from previous studied to verify their validity. In the light of the present investigation, it is found that the momentum, microrotation, thermal and concentration boundary layer thicknesses increase with increasing mixed convection parameter for opposing flow, whereas reverse effect is observed for assisting flow. In addition, the momentum, microrotation, thermal and concentration boundary layer thicknesses decrease with an increase in the value of Hartmann number for both cases of slip and no-slip conditions. Further, the shrinking parameter enhanced the momentum, thermal and concentration boundary layer thicknesses. Furthermore, the dual solutions is only possible for when n = 1.

Acknowledgement

We would like to express our very sincerely thanks to the

reviewers for very useful suggestions to improve the version

of the paper. The first author would like to acknowledge the

Dean of Science for the financial support for this research.

References

[1] K. Hiemenz, Die Grenzschicht an einem in den gleichförmigen Flussigkeitsstrom eingetauchten geraden Kreiszylinder, Dinglers Polym. J. 326 (1911) 321-330.

[2] L.J. Crane, Flow past a stretching plate, Z. Angew. Math. Phys. 21 (1970) 645-647.

[3] T.C. Chiam, Stagnation-point flow towards a stretching plate, J. Phys. Soc. Jpn. 63 (1994) 639-646.

[4] M. Kumari, N. Nath, Flow and heat transfer in a stagnationpoint flow over a stretching sheet with a magnetic field, Mech. Res. Commun. 26 (4) (1999) 469-478.

[5] T.R. Mahapatra, A.S. Gupta, Magnetohydrodynamic stagnation-point flow towards a stretching sheet, Acta Mech. 152 (2001) 191-196.

[6] P.R. Sharma, G. Singh, Effects of variable thermal conductivity and heat source/sink on MHD flow near a stagnation point on a linearly stretching sheet, J. Appl. Fluid Mech. 2 (2008) 13-21.

[7] M. Ashraf, M. Rashid, Radiation effect on MHD boundary layer stagnation point flow towards a heated shrinking sheet, World Appl. Sci. J. 13 (7) (2011) 1748-1756.

[8] A.A.M. Mostafa, E.W. Shimaa, MHD stagnation point flow of a micropolar fluid towards a moving surface with radiation, Meccanica 47 (2012) 1119-1130.

[9] K. Bhattacharyya, G.C. Layek, G.C. Seth, Soret and Dufour effects on convective heat and mass transfer in stagnation-point flow towards a shrinking surface, Phys. Scr. 89 (9) (2014) 1-10.

[10] A. Zeeshan, R. Ellahi, M. Hassan, Magnetohydrodynamic flow of water/ethylene glycol based nanofluids with natural convection through a porous medium, Eur. Phys. J. Plus 129 (261) (2014) 1-10.

[11] R. Ellahi, M.M. Bhatti, A. Riaz, M. Sheikholeslami, Effects of magnetohydrodynamics on peristatic flow of Jeffrey fluid in a rectangular duct through a porous medium, J. Porous Media 17 (2014) 143-157.

[12] A.A. Khan, R. Ellahi, M.M. Gulzar, M. Sheikholeslami, Effects of heat transfer on peristaltic motion of Oldroyd fluid in the presence of inclined magnetic field, J. Magnet. Mag. Mat. 372 (2014) 97-106.

[13] R. Ellahi, S.U. Rahman, S. Nadeem, K. Vafai, The blood flow of Prandtl fluid through a tapered stenosed arteries in permeable walls with magnetic field, Commun. Thero. Phys. 63 (2015) 353358.

[14] N.S. Akbar, M. Raza, R. Ellahi, Influence of induced magnetic field and heat flux with the suspension of carbon nanotubes for the peristaltic flow in a permeable channel, J. Magnet. Mag. Mat. 381 (2015) 405-415.

[15] M.S. Kandelousi, R. Ellahi, Simulation of ferrofluid flow for magnetic drug targeting using the lattice Boltzmann method, Zeitsch. für Naturfor. A 70 (2015) 115-124.

[16] M. Miklavcic, C.Y. Wang, Viscous flow due to a shrinking sheet, Q. Appl. Math. 64 (2006) 283-290.

[17] Muhaimin, R. Kandasamy, A.B. Khamis, Effects of heat and mass transfer on nonlinear MHD boundary layer flow over a shrinking sheet in the presence of suction, Appli. Math. Mech. 29 (10) (2008) 1309-1317.

[18] T. Fang, J. Zhang, Closed-form exact solutions of MHD viscous flow over a shrinking sheet, Commun. Nonlinear Sci. Numer. Sim. 14 (7) (2009) 2853-2857.

[19] C.Y. Wang, Stagnation flow towards a shrinking sheet, Int. J. Non-Linear Mech. 43 (2008) 377-382.

[20] T. Fang, J. Zhang, S.S. Yao, Viscous flow over an unsteady shrinking sheet with mass transfer, Chin. Phys. Lett. 26 (2009) 014703-1-014703-4.

[21] T. Fan, H. Xu, I. Pop, Unsteady stagnation flow and heat transfer towards a shrinking sheet, Int. Commun. Heat Mass Trans. 37 (2010) 1440-1446.

[22] M. Suali, N.M.A. Nik Long, N.M. Ariffin, Unsteady stagnation point flow and heat transfer over a stretching/shrinking sheet with suction or injection, J. Appl. Math. 2012 (2012) 1-12.

[23] K. Bhattacharyya, Heat transfer analysis in unsteady boundary layer stagnation-point flow towards a shrinking/stretching sheet, Ain Shams Eng. J. 4 (2013) 259-264.

[24] A.C. Eringen, Theory of micropolar fluids, J. Math. Mech. 16 (1966) 1-18.

[25] J. Peddieson, R.P. McNitt, Boundary layer theory for a micropolar fluid, Recent Adv. Eng. Sci. 5 (1970) 405-416.

[26] T. Ariman, M.A. Turk, N.D. Sylvester, Microcontinuum fluid mechanics—review, Int. J. Eng. Sci. 11 (1973) 905-930.

[27] T. Ariman, M.A. Turk, N.D. Sylvester, Application of microcontinuum fluid mechanics, Int. J. Eng. Sci. 12 (1974) 273-283.

[28] G. tukaszewicz, Micropolar Fluids: Theory and Application, Birkhauser, Basel, 1999.

[29] A.C. Eringen, Microcontinuum Field Theories II: Fluent Media, Springer, New York, 2001.

[30] G.C. Guram, C. Smith, Comput Stagnation flows of micropolar fluids with strong and weak interactions, Math. Appl. 6 (1980) 213-223.

[31] R.S.R. Gorla, Micropolar boundary layer at a stagnation point, Int. J. Eng. Sci. 21 (1983) 25-34.

[32] R. Nazar, N. Amin, D. Filip, I. Pop, Stagnation point flow of a micropolar fluid towards a stretching sheet, Int. J. Non-Linear Mech. 39 (2004) 1227-1235.

[33] A. Ishak, R. Nazar, I. Pop, Mixed convection stagnation point flow of a micropolar fluid towards a stretching sheet, Meccanica 43 (2008) 411-428.

[34] A. Ishak, Y.Y. Lok, I. Pop, Stagnation-point flow over a shrinking sheet in a micropolar fluid, Chem. Eng. Commun. 197 (2010) 1417-1427.

[35] M. Ashraf, S. Bashir, Numerical simulation of MHD stagnation point flow and heat transfer of a micropolar fluid towards a heated shrinking sheet, Int. J. Num. Meth. Fluids 69 (2) (2012) 384-398.

[36] F. Aman, A. Ishak, I. Pop, Slip effects on mixed convective stagnation-point flow and heat transfer over a vertical surface, in: Prosiding Seminar Kebangsaan Aplikasi Sains dan Matematik, 2013.

[37] K.S.N. Reddy, M.S. Babu, S.V.K. Varma, N.B. Reddy, MHD stagnation point flow of a micropolar fluid over a stretching surface with heat source/sink, chemical reaction and viscous dissipation, Int. J. Eng. Inven. 3 (10) (2014) 15-26.

[38] T. Hayat, M. Qasim, S.A. Shehzad, A. Alsaedi, Unsteady stagnation point flow of second grade fluid with variable free stream, Alexandria Eng. J. 53 (2014) 455-461.

[39] R. Ellahi, S.U. Rahman, M.M. Gulzar, K. Vafai, A mathematical study of non-Newtonian micropolar fluid in arterial blood flow through composite stenosis, Appl. Math. Inform. Sci. 8 (2014) 1567-1573.

[40] R. Ellahi, S.U. Rahman, S. Nadeem, N.S. Akbar, Influence of heat and mass transfer on micropolar fluid of blood flow through a tapered stenosed arteries with permeable walls, J. Comp. Theor. Nanosci. 11 (2014) 1156-1163.

[41] Aurangzaib, S. Shafie, Slip effect on an unsteady MHD stagnation-point flow of a micropolar fluid towards a shrinking sheet with thermophoresis effect, Int. J. Comp. Meth. Eng. Sci. Mech. 16 (2015) 285-291.

[42] C.Y. Wang, Analysis of viscous flow due to a stretching sheet with surface slip and suction, Nonlinear Anal.: Real World Appl. 10 (2009) 375-380.

[43] F. Aman, A. Ishak, Boundary layer flow and heat transfer over a permeable shrinking sheet with partial slip, J. Appl. Sci. Res. 6 (8) (2010) 111-115.

[44] K. Das, Slip effects on heat and mass transfer in MHD micropolar fluid flow over an inclined plate with thermal radiation and chemical reaction, Int. J. Numer. Meth. Fluids 70 (1) (2012) 96-113.

[45] K. Das, Slip effects on MHD mixed convection stagnation point flow of a micropolar fluid towards a shrinking vertical sheet, Comput. Math. Appl. 63 (2012) 255-267.

[46] R. Ellahi, M. Hameed, Numerical analysis of steady non-Newtonian flows with heat transfer analysis, MHD and nonlinear slip effects, Int. J. Num. Meth. Heat Fluid Flow 22

(2012) 24-38.

[47] A. Zeeshan, R. Ellahi, Series solutions for nonlinear partial differential equations with slip boundary conditions for non-Newtonian MHD fluid in porous space, Appl. Math. Inf. Sci. 7

(2013) 257-265.

[48] R. Ellahi, A. Riaz, S. Abbasbandy, T. Hayat, K. Vafai, A study on the mixed convection boundary layer flow and heat transfer over a vertical slender cylinder, Therm. Sci. 18 (2014) 12471258.

[49] H. Chen, Mixed convection unsteady stagnation-point flow towards a stretching sheet with slip effects, Math. Probl. Eng. 2014 (2014) 1-7.

[50] H. Xu, Series solutions of unsteady boundary layer flow of a micropolar fluid near the forward stagnation point of a plane surface, Acta Mech. 184 (2006) 87-91.

[51] S.K. Jena, M.N. Mathur, Similarity solution for laminar free convection flow of a thermomicropolar fluid past a nonisothermal flat plate, Int. J. Eng. Sci. 19 (1981) 1431-1439.

[52] T. Cebeci, P. Bradshaw, Physical and Computational Aspects of Convective Heat Transfer, Springer, New York, 1988.