Special Issue

## Mathematical Modeling of Heat and Mass Transfer in Energy Science and Engineering 2014

View this Special Issue

Research Article | Open Access

Volume 2014 |Article ID 435697 | https://doi.org/10.1155/2014/435697

Hui Chen, "Mixed Convection Unsteady Stagnation-Point Flow towards a Stretching Sheet with Slip Effects", Mathematical Problems in Engineering, vol. 2014, Article ID 435697, 7 pages, 2014. https://doi.org/10.1155/2014/435697

# Mixed Convection Unsteady Stagnation-Point Flow towards a Stretching Sheet with Slip Effects

Accepted30 Mar 2014
Published17 Apr 2014

#### Abstract

The paper studies the unsteady mixed convection flow of an incompressible viscous fluid about a stagnation point on a stretching sheet in presence of velocity and thermal slips. The governing equations are transformed into the ordinary differential equations by using similarity transformations. The transformed equations are solved numerically by an efficient shooting method. The characteristics of the flow and heat transfer features for governing parameters are analyzed and discussed for both the assisting and opposing flows. It is found that dual solutions exist for certain range of buoyancy parameter which again depend on the unsteadiness parameter and the slip parameters (i.e., and ). The numerical results show that the increase of unsteadiness parameter and the slip effects cause increment in the existence range of similarity solution. The effects of unsteadiness parameter, the velocity ratio parameter, and the velocity and thermal slip parameters on the velocity and temperature distributions are analyzed and discussed.

#### 1. Introduction

The stagnation-point flow due to a stretching sheet has received much attention because of its important practical applications in industry and practical applications, such as extrusion of polymers, glass fiber, the cooling of metallic plate, and the aerodynamics. As a hot topic in fluid mechanics, the two-dimensional stagnation flow was first studied by Hiemenz [1]. The result has been later extended to axisymmetric case by Homann [2] and improved by Howarth [3]. Following these works, various aspects of stagnation-point flow and heat transfer have been studied and many literatures have been generated on this problem [411]. The flow becomes time dependent in certain aspects, but the physical situation described in the above studies is under the condition of a steady state. Consequently, it is necessary to consider the unsteadiness of the flow. Nazar et al. [12] considered an unsteady boundary layer flow in the region of the stagnation point on a stretching sheet, while Bhattacharyya [13, 14] investigated the unsteady stagnation-point flow over a shrinking sheet. Sharma and Singh [15] also investigated an unsteady flow near a stagnation point on a stretching sheet in the presence of a time-dependent free stream. Some important properties of unsteady flows on a stretching sheet were described by Bachok et al. [16], Ishak et al. [17], and Hayat et al. [18].

The mixed convection in stagnation flow is a topic of significance in fluid mechanics when the buoyancy forces due to the temperature difference between the surface and the free stream become large, in the sense that both the flow and thermal fields are greatly affected by the buoyancy forces. Much interest has been given to the free and forced convection stagnation-point flows of a viscous fluid. Devi et al. [19] studied the unsteady laminar mixed convection in two-dimensional stagnation-point flows around heated surfaces by taking both cases of an arbitrary wall temperature and arbitrary surface heat flux variations. The unsteady mixed convection flow of a micropolar fluid was studied by Lok et al. [20], where they found the smooth transition from the initial unsteady-state flow to the final steady-state flow. Recently, Ishak et al. [21] reported the existence of dual solutions for both assisting and opposing flows of an electrically conducting fluid past a vertical permeable flat plate. Hayat et al. [22] investigated the effects of mixed convection unsteady stagnation-point flow of viscous fluid with variable free stream velocity.

In all the above studies, the assumption of the flow field obeys the conventional no-slip condition at the boundary. However, the assumption of the conventional no-slip condition at the boundary is not true and should be replaced by partial slip boundary condition in certain situations [25]. With a slip at the wall boundary, the flow behavior and the shear stress in the fluid are quite different from those in the no-slip cases. Wang [26] gave an exact solution of the Navier-Stokes equations for the flow due to a stretching boundary with slip. Later, he [27] considered the effect of stagnation slip flow on the heat transfer from a moving plate. Ariel et al. [28, 29] studied the effects of slip on the flow of an elastic-viscous fluid with some other physical features. Many researchers have also investigated the slip flows in different configurations recently [3035]. It can be pointed out here that less work has been done on the mixed convection flow with slip effect at the boundary. Cao and Baker [36] considered the slip effects on the mixed convective flow and heat transfer from a vertical plate and reported the local nonsimilarity solutions. Mukhopadhyay [37, 38] investigated effects of slip on unsteady mixed convective flow and heat transfer past a stretching surface and a porous stretching surface. Bhattacharyya et al. [39] studied the mixed convective flow adjacent to a vertical permeable stretching sheet in porous medium with slip effects. The similarity solution of the mixed convection boundary layer flow near the stagnation-point on a vertical surface with the slip effect was studied by Aman et al. [40]. Very recently, Nik Long et al. [41] studied mixed convection boundary layer caused by time-dependent velocity and the surface temperature in the two-dimensional unsteady stagnation-point flow over a stretching vertical sheet with the no-slip boundary condition.

Motivated by the above studies, in this paper we investigate the behaviour of the mixed convection unsteady stagnation-point flow towards a stretching sheet with slip effect on the boundary. The momentum and energy equations are solved numerically by a shooting method. The effects of the key parameters on the flow and heat transfer characteristics are analyzed and discussed. To the best of our knowledge, this problem has not been studied before.

#### 2. Flow Analysis

Consider an unsteady two-dimensional flow of a viscous and incompressible fluid in the vicinity of a stagnation-point towards a vertical stretching sheet. The sheet stretching velocity is where is a parameter showing the unsteadiness of the problem and is a constant with for a stretching sheet. The free stream velocity is , where is the strength of the stagnation flow. The surface temperature of the stretching sheet varies with the distance as , where is the constant free stream temperature with . The particular forms of the above expressions for , , and have been chosen in order to transform the governing partial differential equations into a set of ordinary differential equations, thereby facilitating the exploration of the effects of the controlling parameters. It should be noticed that the expressions for , , and are valid for time , and , , have dimension time−1.

Using the boundary layer approximations, the governing equations for this problem are where and are the velocities in the and directions, is the kinematic viscosity, is the acceleration due to gravity, is the thermal expansion coefficient, is the fluid temperature, and is the thermal diffusivity. The boundary conditions with partial slip are given by Here and are the velocity slip factor and the thermal slip factor and and are the initial values of velocity and thermal slip factors, respectively.

To obtain the similarity solution, we now introduce the following similarity transformations: where is the steam function which automatically assures mass conservation (1). The velocity components are obtained as Substituting (5) into (2) and (3), we get where is the unsteadiness parameter, is the Prandtl number, is the buoyancy or mixed convection parameter defined as with being the local Grashof number, and is the local Reynolds number. It should be noticed that is a constant with and corresponding to assisting flow and opposing flow, respectively. The corresponding boundary conditions (4) become where and are the dimensionless velocity and thermal slip parameters, respectively, and is the ratio of stretching parameter and free steam velocity parameter.

The physical quantities of interest are the skin friction coefficient and the local Nusselt number , which are defined as where the surface shear stress and the surface heat flux are given by Using the similarity transformations (5) and the relation (13) into (12), we get with being the local Reynolds number.

#### 3. Results and Discussion

The nonlinear ordinary differential (7) subject to the boundary conditions (8)–(11) have been numerically solved using an improved shooting method described by Zheng et al. [42]. The results show the influences of some important nondimensional parameters on the feature of the flow and the heat transfer characteristics. In order to validate the numerical method used in this study and to judge the accuracy of the present analysis, the present values of the skin friction coefficient and the local Nusselt number for different values of are compared with those obtained by Ishak et al. [23] and Pal [24] when and (absence of the slip effect), (the steady-state flow), , and for assisting and opposing flow. The quantitative comparisons are shown in Tables 1 and 2 and found to be in very good agreement.

 Pr Ishak et al. [23] Pal [24] Present work (0) (0) (0) (0) (0) (0) 0.72 0.3645 1.0931 0.36449 1.09331 0.36449 1.09311 6.8 0.1804 3.2902 0.18042 3.28957 0.18042 3.28957 20.0 0.1175 5.6230 0.11750 5.62014 0.11750 5.62013 40.0 0.0873 7.9463 0.08724 7.93831 0.08724 7.93830 60.0 0.0729 9.7327 0.07284 9.71801 0.07284 9.71800 80.0 0.0640 11.2413 0.06394 11.21875 0.06394 11.21873 100.0 0.0578 12.5726 0.05773 12.54113 0.05773 12.54109
 Pr Ishak et al. [23] Pal [24] Present work (0) (0) (0) (0) (0) (0) 0.72 −0.3852 1.0293 −0.38519 1.02925 −0.38519 1.02925 6.8 −0.1832 3.2466 −0.18323 3.24609 −0.18323 3.24608 20.0 −0.1183 5.5923 −0.11831 5.58960 −0.11831 5.58959 40.0 −0.0876 7.9227 −0.08758 7.91491 −0.08758 7.91489 60.0 −0.0731 9.7126 −0.07304 9.69818 −0.07304 9.69817 80.0 −0.0642 11.2235 −0.06408 11.20118 −0.06408 11.20117 100.0 −0.0579 12.5564 0.05783 12.52519 −0.05783 12.52514

The skin friction coefficient, local Nusselt number, velocity, and temperature profiles are shown in Figures 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, and 12. Figures 1 and 2 show the variation of skin friction coefficient and local Nusselt number with for the unsteadiness parameter in the absence of slip when the velocity ratio parameter and the Prandtl number . Moreover, the variations of and in the presence of slip (i.e., , ) are, respectively, shown in Figures 3 and 4.

It is seen from Figures 14 that the dual solutions exist for the buoyancy opposing flow (), whereas the solution is unique for the assisting flow (). It is worth mentioning that for all values of the parameters as shown in Figures 3 and 4, which means that the heat is transferred from the hot surface to the cool fluid. In the dual solutions range, we identify the first solution and second solution in the following discussion by the principle that the first solution has the higher values of and for a given than the second solution. It is observed that no solutions exist when , where is the critical value for which (7) have no solutions. Further, these figures show that the solutions could be obtained for all values of (assisting flow), while for (opposing flow), solutions exist only when . Based on our computations, we find that , respectively, for under nonslip condition shown in Figures 1 and 2, and , respectively, for under slip condition shown in Figures 3 and 4. Thus, it can be concluded that the solution domain expands as the unsteadiness parameter increases. Moreover, the velocity and thermal slip cause more increment in the existence range of similarity solution. It is also noteworthy that the impact of on in Figures 2 and 4 is more pronounced than on the in Figures 1 and 3.

In Figures 58, velocity profiles and temperature profiles are shown for different values of the unsteadiness parameter and velocity ratio parameter , respectively, in presence of slip (, ) for both the first and second solution branches. It is seen from Figures 5 and 7 that, for the first solution branch, the velocity increases with unsteadiness parameter or the velocity ratio parameter and this implies an accompanying reduction of the thickness of the momentum boundary layer. The opposite trend can be observed for the second solution branch. From Figures 6 and 8, it is observed that the temperature decreases with the increase of or for the first solution branch while a different trend is observed for the second solution branch.

The samples of the and for the selected values of the slip parameters are investigated in Figures 912. The effect of velocity slip parameter is shown in Figures 9 and 10, while the influence of thermal slip parameter is shown in Figures 11 and 12. It is clear that effects of the two type slips on velocity and temperature profiles are opposite. For the first solution branch, an increase in velocity slip parameter would decrease and increase , while an increase in thermal slip parameter would increase and decrease . For the second solution branch, the velocity and temperature profiles ( and ) show the opposite trend.

It is evident from Figures 512 that the first solution displays the thinner boundary layer thickness compared with the second solution. From Figures 5, 7, 9, and 10, it is interesting to note that for the second solution branch the value of initially decreases with to a negative value and for small it starts to increase and ultimately it becomes the positive value 1. Thus the velocity profiles exhibit reverse flow () near the wall (). Moreover, the samples of velocity and temperature profiles presented in Figures 512 show that the boundary conditions (11) are asymptotically satisfied, which supports the validity of the obtained numerical results.

#### 4. Conclusions

In this paper, we have studied the mixed convection unsteady boundary layer flow and heat transfer about a stagnation-point towards a stretching sheet in the presence of both velocity and thermal slip conditions at the boundary. The governing equations are reduced to the ordinary differential equations by the similarity transformation and then numerically solved to obtain the influence of key parameters on the skin friction coefficient, the local Nusselt number, the velocity, and temperature profiles. The existence and duality of solutions are displayed in Figures 14 with the buoyancy parameter , and the range of for which the similarity solution exists increases with the unsteadiness parameter . Moreover, the velocity and thermal slip parameters cause more increment in the existence range of similarity solution. The effects of the unsteadiness parameter , velocity ratio parameter , velocity slip parameter , and thermal slip parameter on the velocity and temperature profiles are shown in Figures 512. It is noticed that the reverse flow occurs near the sheet. Different flow behavior is observed with the first and second solution branches.

#### Conflict of Interests

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

#### Acknowledgments

This work is supported by the National Natural Science Foundation of China (Grant no. 51305080), Fujian Provincial Natural Science Foundation of China (Grant no. 2012J05089), and Visiting Scholar Foundation of Key Lab in University (Grant no. GZKF-201217).

#### References

1. K. Hiemenz, “Die Grenzschicht in Einem in Dem Gleichformingen Flussigkeitsstrom Eingetauchten Gerade Kreiszlinder,” Dingler's Polytechnic Journal, vol. 326, pp. 321–410, 1911. View at: Google Scholar
2. F. Homann, “Die Einfluss grosse Zähigkeit bei der Strömung um der Zylinder und um die Kugel,” Zeitschrift für Angewandte Mathematik und Mechanik, vol. 16, no. 3, pp. 153–164, 1936. View at: Publisher Site | Google Scholar
3. L. Howarth, “On the solution of the laminar boundary layer equations,” Proceedings of the Royal Society of London, vol. 164, no. 919, pp. 547–579, 1938. View at: Publisher Site | Google Scholar
4. T. C. Chiam, “Stagnation-point flow towards a stretching plate,” Journal of the Physical Society of Japan, vol. 63, no. 6, pp. 2443–2444, 1994. View at: Publisher Site | Google Scholar
5. T. R. Mahapatra and A. S. Gupta, “Heat transfer in stagnation-point flow towards a stretching sheet,” Heat and Mass Transfer, vol. 38, no. 6, pp. 517–521, 2002. View at: Publisher Site | Google Scholar
6. G. C. Layek, S. Mukhopadhyay, and S. A. Samad, “Heat and mass transfer analysis for boundary layer stagnation-point flow towards a heated porous stretching sheet with heat absorption/generation and suction/blowing,” International Communications in Heat and Mass Transfer, vol. 34, no. 3, pp. 347–356, 2007. View at: Publisher Site | Google Scholar
7. A. Ishak, R. Nazar, and I. Pop, “Mixed convection stagnation point flow of a micropolar fluid towards a stretching sheet,” Meccanica, vol. 43, no. 4, pp. 411–418, 2008.
8. K. Bhattacharyya, S. Mukhopadhyay, and G. C. Layek, “Reactive solute transfer in magnetohydrodynamic boundary layer stagnation-point flow over a stretching sheet with suction/blowing,” Chemical Engineering Communications, vol. 199, no. 3, pp. 368–383, 2012. View at: Publisher Site | Google Scholar
9. K. Bhattacharyya, Md. G. Arif, and W. Ali Pramanik, “MHD boundary layer stagnation-point flow and mass transfer over a permeable shrinking sheet with suction/blowing and chemical reaction,” Acta Technica, vol. 57, no. 1, pp. 1–15, 2012. View at: Google Scholar | MathSciNet
10. K. Bhattacharyya and K. Vajravelu, “Stagnation-point flow and heat transfer over an exponentially shrinking sheet,” Communications in Nonlinear Science and Numerical Simulation, vol. 17, no. 7, pp. 2728–2734, 2012. View at: Publisher Site | Google Scholar
11. K. Bhattacharyya, “Heat transfer in boundary layer stagnation-point flow towards a shrinking sheet with non-uniform heat flux,” Chinese Physics B, vol. 22, no. 7, Article ID 074705, 6 pages, 2013. View at: Google Scholar
12. R. Nazar, N. Amin, D. Filip, and I. Pop, “Unsteady boundary layer flow in the region of the stagnation point on a stretching sheet,” International Journal of Engineering Science, vol. 42, no. 11-12, pp. 1241–1253, 2004.
13. K. Bhattacharyya, “Dual solutions in unsteady stagnation-point flow over a shrinking sheet,” Chinese Physics Letters, vol. 28, no. 8, Article ID 084702, 2011. View at: Publisher Site | Google Scholar
14. K. Bhattacharyya, “Heat transfer in unsteady boundary layer stagnation-point flow towards a shrinking sheet,” Ain Shams Engineering Journal, vol. 44, no. 2, pp. 259–264, 2013. View at: Google Scholar
15. P. R. Sharma and G. Singh, “Unsteady flow about a stagnation point on a stretching sheet in the presence of variable freestream,” Thammasat International Journal of Science and Technology, vol. 13, no. 1, pp. 11–16, 2008. View at: Google Scholar
16. N. Bachok, A. Ishak, and R. Nazar, “Flow and heat transfer over an unsteady stretching sheet in a micropolar fluid,” Meccanica, vol. 46, no. 5, pp. 935–942, 2011.
17. A. Ishak, R. Nazar, and I. Pop, “Heat transfer over an unsteady stretching permeable surface with prescribed wall temperature,” Nonlinear Analysis: Real World Applications, vol. 10, no. 5, pp. 2909–2913, 2009.
18. T. Hayat, M. Qasim, and Z. Abbas, “Homotopy solution for the unsteady three-dimensional MHD flow and mass transfer in a porous space,” Communications in Nonlinear Science and Numerical Simulation, vol. 15, no. 9, pp. 2375–2387, 2010.
19. C. D. S. Devi, H. S. Takhar, and G. Nath, “Unsteady mixed convection flow in stagnation region adjacent to a vertical surface,” Wärme- und Stoffübertragung, vol. 26, no. 2, pp. 71–79, 1991. View at: Publisher Site | Google Scholar
20. Y. Y. Lok, N. Amin, and I. Pop, “Unsteady mixed convection flow of a micropolar fluid near the stagnation point on a vertical surface,” International Journal of Thermal Sciences, vol. 45, no. 12, pp. 1149–1157, 2006. View at: Publisher Site | Google Scholar
21. A. Ishak, R. Nazar, N. Bachok, and I. Pop, “MHD mixed convection flow near the stagnation-point on a vertical permeable surface,” Physica A: Statistical Mechanics and its Applications, vol. 389, no. 1, pp. 40–46, 2010. View at: Publisher Site | Google Scholar
22. T. Hayat, S. A. Shehzad, A. Rafique, and M. Y. Malik, “Mixed convection unsteady stagnation point flow over a stretching sheet with heat transfer in the presence of variable free stream,” International Journal for Numerical Methods in Fluids, vol. 68, no. 4, pp. 483–493, 2012.
23. A. Ishak, R. Nazar, and I. Pop, “Mixed convection boundary layers in the stagnation-point flow toward a stretching vertical sheet,” Meccanica, vol. 41, no. 5, pp. 509–518, 2006.
24. D. Pal, “Heat and mass transfer in stagnation-point flow towards a stretching surface in the presence of buoyancy force and thermal radiation,” Meccanica, vol. 44, no. 2, pp. 145–158, 2009.
25. H. I. Andersson, “Slip flow past a stretching surface,” Acta Mechanica, vol. 158, no. 1-2, pp. 121–125, 2002.
26. C. Y. Wang, “Flow due to a stretching boundary with partial slip—an exact solution of the Navier-Stokes equations,” Chemical Engineering Science, vol. 57, no. 17, pp. 3745–3747, 2002. View at: Publisher Site | Google Scholar
27. C. Y. Wang, “Stagnation slip flow and heat transfer on a moving plate,” Chemical Engineering Science, vol. 61, no. 23, pp. 7668–7672, 2006. View at: Publisher Site | Google Scholar
28. P. D. Ariel, T. Hayat, and S. Asghar, “The flow of an elastico-viscous fluid past a stretching sheet with partial slip,” Acta Mechanica, vol. 187, no. 1–4, pp. 29–35, 2006.
29. P. D. Ariel, “Two dimensional stagnation point flow of an elastico-viscous fluid with partial slip,” Zeitschrift für Angewandte Mathematik und Mechanik, vol. 88, no. 4, pp. 320–324, 2008.
30. P. D. Ariel, “Axisymmetric flow due to a stretching sheet with partial slip,” Computers & Mathematics with Applications, vol. 54, no. 7-8, pp. 1169–1183, 2007.
31. T. Fang, J. Zhang, and S. Yao, “Slip MHD viscous flow over a stretching sheet—an exact solution,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 11, pp. 3731–3737, 2009. View at: Publisher Site | Google Scholar
32. K. Bhattacharyya, S. Mukhopadhyay, and G. C. Layek, “Slip effects on an unsteady boundary layer stagnation-point flow and heat transfer towards a stretching sheet,” Chinese Physics Letters, vol. 28, no. 9, Article ID 094702, 5 pages, 2011. View at: Publisher Site | Google Scholar
33. S. Mukhopadhyay and R. S. R. Gorla, “Effects of partial slip on boundary layer flow past a permeable exponential stretching sheet in presence of thermal radiation,” Heat and Mass Transfer, vol. 48, no. 10, pp. 1773–1781, 2012. View at: Publisher Site | Google Scholar
34. S. Mukhopadhyay, P. R. De, K. Bhattacharyya, and G. C. Layek, “Slip effects on mixed convection flow along a stretching cylinder,” International Journal of Heat and Technology, vol. 30, no. 2, pp. 19–24, 2012. View at: Google Scholar
35. K. Bhattacharyya, K. Vajravelu, and T. Hayat, “Slip effects on the parametric space and the solution for boundary layer flow of casson fluid over a porous stretching/shrinking sheet,” International Journal of Fluid Mechanics Research, vol. 40, no. 6, pp. 482–493, 2013. View at: Publisher Site | Google Scholar
36. K. Cao and J. Baker, “Slip effects on mixed convective flow and heat transfer from a vertical plate,” International Journal of Heat and Mass Transfer, vol. 52, no. 15-16, pp. 3829–3841, 2009.
37. S. Mukhopadhyay, “Effects of slip on unsteady mixed convective flow and heat transfer past a stretching surface,” Chinese Physics Letters, vol. 27, no. 12, Article ID 124401, 4 pages, 2010. View at: Publisher Site | Google Scholar
38. S. Mukhopadhyay, “Effects of slip on unsteady mixed convective flow and heat transfer past a porous stretching surface,” Nuclear Engineering and Design, vol. 241, no. 8, pp. 2660–2665, 2011. View at: Publisher Site | Google Scholar
39. K. Bhattacharyya, M. S. Uddin, and G. C. Layek, “Effect of partial slip on boundary layer mixed convective flow adjacent to a vertical permeable stretching sheet in porous medium,” Acta Technica, vol. 58, no. 1, pp. 27–39, 2013. View at: Google Scholar
40. F. Aman, A. Ishak, and I. Pop, “Mixed convection boundary layer flow near stagnation-point on vertical surface with slip,” Applied Mathematics and Mechanics, vol. 32, no. 12, pp. 1599–1606, 2011.
41. N. M. A. Nik Long, L. F. Koo, T. J. Wong, and M. Suali, “Mixed convection boundary layers with prescribed temperature in the unsteady stagnation point flow toward a stretching vertical sheet,” Mathematical Problems in Engineering, vol. 2013, Article ID 195360, 10 pages, 2013. View at: Publisher Site | Google Scholar | MathSciNet
42. L. C. Zheng, C. Liang, and X. X. Zhang, “A numerical method for solving the boundary layer equations of laminar natural convention about a vertical plate,” Journal of University of Science and Technology Beijing, vol. 14, no. 1, pp. 33–35, 2007. View at: Publisher Site | Google Scholar

Copyright © 2014 Hui Chen. 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.