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.
Unsteady; Partial slip effect; MHD; Micropolar fluid; Permeable shrinking sheet
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 reactor, 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 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. Miklavčič 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 stagnation-point 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] and [27] and books by Łukaszewicz [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.
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 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 , the temperature and the concentration of the sheet are proportional to the distance from the stagnation point, where and with and being the uniform temperature and concentration of the ambient fluid. The velocity of the flow external to the boundary layer is . Under the usual boundary layer approximation, along with Boussinesq approximations the governing equations are governed by the following equations:
|
(1) |
|
(2) |
|
(3) |
|
(4) |
|
(5) |
where are the velocity components along the axes, respectively, is the specific heat at constant pressure, is the vortex viscosity, is the electrical conductivity of the fluid, is externally imposed magnetic field strength in the y -direction, is the density of the fluid, is the molecular diffusivity, is the coefficient of thermal expansion, is the components of microrotation or angular velocity whose rotation is in the direction of the – planes, is the spin gradient viscosity given by , is the micro-inertia density, and is the coefficient of viscosity of the fluid. According to Xu [50] we assume that is constant therefore .
|
Figure 1. Physical model and coordinate system.
|
The appropriate physical boundary conditions for the problem under study are given by
|
|
(6) |
|
The sheet velocity , temperature and concentration are assumed to vary proportional to the distance from the stagnation point and the velocity of the flow external to the boundary layer is are taken of the form:
|
(7) |
where and are positive constants, and is a constant such that . It is noted the case when is called a strong concentration. In this case, near the wall implying that the concentrated particle flows in which the micro-elements close to the wall surface are unable to rotate [51]. The case presents the vanishing of the anti-symmetric part of the stress tensor and denotes weak concentrations whereas is used for modelling of turbulent boundary layer flows. In this study, we only consider strong concentration .
Introducing the following non-dimensional variables as
|
(8) |
where and the continuity Eq. (1) is satisfied by introducing a stream function such that
|
The transformed nonlinear partial differential equations are as follows:
|
(9) |
|
(10) |
|
(11) |
|
(12) |
subject to the boundary conditions (6) which become
|
(13) |
where is the suction or injection parameter, is the Hartmann number, is the mixed convection parameter, is the local Grashof number, is the local Reynolds number, is the constant concentration buoyancy parameter, is the Prandtl number, is the slip parameter, is Schmidt number, and with and correspond to a stretching and shrinking sheet, respectively. It should be noted that dual solutions is only possible when which can be seen in Figure 17, Figure 18, Figure 19, Figure 2 and Figure 20.
|
Figure 2. Variation of on the velocity .
|
The quantities of physical interest are the skin-friction coefficient in the x -direction and the Nusselt number and the Sherwood number are defined respectively by,
|
|
|
The nonlinear partial differential Eqs. (9), (10), (11) and (12) subject to boundary conditions (13) 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 has been used and the convergence criterion was set to , which give accuracy to four decimal places. The satisfaction of the outer boundary condition is achieved by considering the boundary layer thickness . 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 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.
Wang [19] | Ishak et al. [34] | Present | |||
---|---|---|---|---|---|
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 2, Figure 20, Figure 3, Figure 4 and Figure 5 illustrate the effects for the mixed convection parameter on the velocity, microrotation, temperature and concentration of the fluid for both assisting and opposing flows. From Figure 2 and Figure 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 leads to an increase in the temperature difference . 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 the case of opposing flow. It is evident from these figures that all curves approach the far field boundary conditions (13) asymptotically.
|
Figure 3. Variation of on the microrotation .
|
|
Figure 4. Variation of on the temperature .
|
|
Figure 5. Variation of on the concentration .
|
Figure 6, Figure 7, Figure 8 and Figure 9 illustrate the effect of shrinking parameter on the velocity, microrotation, temperature and concentration profiles. It is observed from Figure 6 and Figure 7 that the velocity and microrotation profiles decrease with an increase in the value of magnitude of shrinking parameter 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 as shown in Figure 8 and Figure 9. Thus, the thermal and concentration boundary layer thicknesses become thicker and thicker.
|
Figure 6. Variation of on the velocity .
|
|
Figure 7. Variation of on the microrotation .
|
|
Figure 8. Variation of on the temperature .
|
|
Figure 9. Variation of on the concentration .
|
The influence of Hartmann number on the velocity, microrotation, temperature and concentration profiles is displayed in Figure 10, Figure 11, Figure 12 and Figure 13 in the presence of slip and in the absence of slip at the boundary. From Figure 10 and Figure 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 remains positive in the boundary layer region. Here an electrically conducting fluid 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. Figure 12 and Figure 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.
|
Figure 10. Variation of and on the velocity .
|
|
Figure 11. Variation of and on the microrotation .
|
|
Figure 12. Variation of and on the temperature .
|
|
Figure 13. Variation of and on the concentration .
|
Figure 14, Figure 15 and Figure 16 demonstrate the effect of material parameter 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, reduced Nusselt number and reduced Sherwood number decrease with an increase in the value of micropolar parameter . 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 injection. Further, it is seen from these results that the value of and 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.
|
Figure 14. Variation of and on the reduced skin friction versus .
|
|
Figure 15. Variation of and on the reduced Nusselt number versus .
|
|
Figure 16. Variation of and on the reduced Sherwood number versus .
|
−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 |
Figure 17, Figure 18, Figure 19, Figure 2 and Figure 20 show the dual velocity, microrotation, temperature and concentrations for different values of . Examining the dual velocity profiles (Fig. 17) for different values of , it 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 for both first and second solutions as shown in Fig. 18. From Figure 19 and Figure 20, it is found that the temperature and concentration of fluid increase with for both solutions. Further, the second solution exhibits a larger boundary layer thickness compared to the first solution.
|
Figure 17. Variation of on the velocity .
|
|
Figure 18. Variation of on the microrotation .
|
|
Figure 19. Variation of on the temperature .
|
|
Figure 20. Variation of on the concentration .
|
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 .
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.
Published on 12/04/17
Licence: Other
Are you one of the authors of this document?