Mathematical Problems in Engineering

Volume 2014 (2014), Article ID 435697, 7 pages

http://dx.doi.org/10.1155/2014/435697

Research Article

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

Hui Chen1,2

¹School of Mechanical Engineering and Automation, Fuzhou University, Fuzhou, Fujian 350108, China

²The State Key Laboratory of Fluid Power Transmission and Control, Zhejiang University, Hangzhou, Zhejiang 310027, China

Received 29 January 2014; Accepted 30 March 2014; Published 17 April 2014

Academic Editor: Zhijun Zhang

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.

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 [4–11]. 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 [30–35]. 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⁻¹.

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.

Table 1: Values of (0) and (0) for different values of Pr when , , , and for assisting flow.

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

Table 2: Values of (0) and (0) for different values of Pr when , , , and for opposing flow.

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.

Figure 1: Variation of with for different values of under no-slip boundary condition.

Figure 2: Variation of with for different values of under no-slip boundary condition.

Figure 3: Variation of for different values of under slip boundary condition.

Figure 4: Variation of with for different values of under slip boundary condition.

Figure 5: Influence of unsteady parameter on velocity profiles .

Figure 6: Influence of unsteady parameter on temperature profiles .

Figure 7: Influence of velocity ratio parameter on velocity profiles .

Figure 8: Influence of velocity ratio parameter on temperature profiles .

Figure 9: Influence of velocity slip parameters on velocity profiles .

Figure 10: Influence of velocity slip parameters on temperature profiles .

Figure 11: Influence of thermal slip parameters on velocity profiles .

Figure 12: Influence of thermal slip parameters on temperature profiles .

It is seen from Figures 1–4 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 5–8, 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 9–12. 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 5–12 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 5–12 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 1–4 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 5–12. 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

K. Hiemenz, “Die Grenzschicht in Einem in Dem Gleichformingen Flussigkeitsstrom Eingetauchten Gerade Kreiszlinder,” Dingler's Polytechnic Journal, vol. 326, pp. 321–410, 1911.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

K. Bhattacharyya, “Dual solutions in unsteady stagnation-point flow over a shrinking sheet,” Chinese Physics Letters, vol. 28, no. 8, Article ID 084702, 2011.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

H. I. Andersson, “Slip flow past a stretching surface,” Acta Mechanica, vol. 158, no. 1-2, pp. 121–125, 2002.

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.

C. Y. Wang, “Stagnation slip flow and heat transfer on a moving plate,” Chemical Engineering Science, vol. 61, no. 23, pp. 7668–7672, 2006.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.