JAM Journal of Applied Mathematics 1687-0042 1110-757X Hindawi Publishing Corporation 528717 10.1155/2013/528717 528717 Research Article Stagnation-Point Flow toward a Vertical, Nonlinearly Stretching Sheet with Prescribed Surface Heat Flux Wong Sin Wei 1 Awang M. A. Omar 1 Ishak Anuar 2 Liu Chein-Shan 1 Institute of Mathematical Sciences University of Malaya 50603 Kuala Lumpur Malaysia um.edu.my 2 School of Mathematical Sciences Faculty of Science and Technology Universiti Kebangsaan Malaysia 43600 Bangi, Selangor Malaysia ukm.my 2013 11 3 2013 2013 27 11 2012 31 01 2013 31 01 2013 2013 Copyright © 2013 Sin Wei Wong et al. 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.

An analysis is carried out to study the steady two-dimensional stagnation-point flow of an incompressible viscous fluid towards a stretching vertical sheet. It is assumed that the sheet is stretched nonlinearly, with prescribed surface heat flux. This problem is governed by three parameters: buoyancy, velocity exponent, and velocity ratio. Both assisting and opposing buoyant flows are considered. The governing partial differential equations are transformed into a system of ordinary differential equations and solved numerically by finite difference Keller-box method. The flow and heat transfer characteristics for different values of the governing parameters are analyzed and discussed. Dual solutions are found in the opposing buoyant flows, while the solution is unique for the assisting buoyant flows.

1. Introduction

The study of fluid flow and heat transfer due to a stretching surface has significant application in the industrial processes, for example, in polymer sheet extrusion from a die, drawing of plastic films, and manufacturing of glass fiber. The quality of the final product greatly depends on the heat transfer rate at the stretching surface as explained by Karwe and Jaluria [1, 2].

Different from the flow induced by a stretching horizontal plate (see Crane , Weidman and Magyari , and Weidman and Ali ), the effect of the buoyancy force could not be neglected for the vertical plate. There are several works that reported the flow and heat transfer characteristics that are brought about by the buoyancy force . Ramachandran et al.  studied the effect of buoyancy force on the stagnation point flows past a vertically heated surface at rest and found that dual solutions exist in the buoyancy opposing flow region. In the present paper, in addition to the flow under the influence of buoyancy force as discussed by Ramachandran et al. , we discuss the consequent flow and heat transfer characteristics that are also brought about by the stretching sheet with power-law velocity variation. It is worth mentioning that the problems of the stagnation-point flow toward a stretching sheet have been considered by many authors , by considering various flow configurations as well as surface heating conditions.

2. Problem Formulation

Consider a mixed convection stagnation-point flow towards a vertical nonlinearly stretching sheet immersed in an incompressible viscous fluid, as shown in Figure 1. The Cartesian coordinates (x,y) are taken such that the x-axis is measured along the sheet oriented in the upwards or downwards direction and the y-axis is normal to it. It is assumed that the wall stretching velocity is given by Uw=axm and the far field inviscid velocity distribution in the neighborhood of the stagnation point (0,0) is given by U(x)=bxm, V(y)=-bym. The surface heat flux is in the form of qw(x)=cx(5m-3)/2 (see Merkin and Mahmood ), where a, b, c, and m are constants. This qw(x) ensured that the buoyancy parameter is independent of x. For the assisting flow, as shown in Figure 1(a), the x-axis points upwards in the same direction of the stretching surface such that the external flow and the stretching surface induce flow and heat transfer in the velocity and thermal boundary layers, respectively. On the other hand, for the opposing flow, as shown in Figure 1(b), the x-axis points vertically downwards in the same direction of the stretching surface such that the external flow and the stretching surface also induce flow and heat transfer, respectively, in the velocity and thermal boundary layers. The steady boundary layer equations, with Boussinesq approximation, are (1)  ux+vy=0,(2)u  ux+v  uy=UdUdx+ν  2uy2  +  gβ(T-T),(3)uTx+vTy=α2Ty2, subject to the boundary conditions (4)u=Uw(x),v=0,Ty=-qwkaty=0,uU(x),  TTasy, where u    and    v are the velocity components along the x- and   y-axes, respectively, g is the acceleration due to gravity, α is the thermal diffusivity of the fluid, ν is the kinematic viscosity, β is the coefficient of thermal expansion, ρ is the fluid density, and T is the far field ambient constant temperature.

Physical model and coordinate system.

The continuity equation (1) can be satisfied automatically by introducing a stream function ψ such that u=ψ/y and v=-ψ/x. The momentum and energy equations are transformed by the similarity variables (5)η=(Uνx)1/2y,ψ=[νxU]1/2f(η),θ(η)=k(T-T)qw(Uνx)1/2 into the following nonlinear ordinary differential equations: (6)    f′′′+m+12ff′′+m(1-f2)+λθ=0,    1Prθ′′+m+12fθ-(2m-1)fθ=0. The transformed boundary conditions are (7)f(0)=0,f(0)=ε,θ(0)=-1,f(η)1,θ(η)0asη, where ε=a/b. Here primes denote differentiation with respect to η, λ=Grx/Rex5/2 is the buoyancy or mixed convection parameter, Pr=ν/α is the Prandtl number, Grx=gβqwx4/(kν2) is the local Grashof number, and Rex=Ux/ν is the local Reynolds number. We note that λ is a constant, with λ>0 corresponds to the assisting flow and λ<0 denotes the opposing flow whilst λ=0 is for forced convective flow. For the forced convection flow (λ=0), the corresponding temperature problem possesses a larger similarity-solution domain than the mixed convection (λ0) one. The reason is that in the forced convection case, the existence of similarity solutions does not require the restriction of the applied wall heat flux to the special form qw(x)=cx(5m-3)/2. Namely, in the forced convection case, the much weaker assumption qw(x)=cxn suffices for the similarity reduction of the problem, where the flux exponent n does not depend on the velocity exponent m in any way.

The main physical quantities of interest are the values of f′′(0), being a measure of the skin friction, and the non-dimensional wall temperature θ(0). Our main aim is to find how the values of f′′(0) and θ(0) vary in terms of the parameters λ and m.

3. Results and Discussion

Equations (6) subject to the boundary conditions (7) are integrated numerically using a finite difference scheme known as the Keller box method . Numerical results are presented for different physical parameters. To conserve space, we consider the Prandtl number as unity throughout this paper. The results presented here are comparable very well with those of Ramanchandran et al. . For no buoyancy effects λ=0 and m=1, comparison of the values of f′′(0) was made with those of Wang  as presented Table 1, which shows a favourable agreement.

Comparison of the values of f′′(0) with those of Wang  when m=1 and λ=0.

ε Wang  Present
5 −10.26475 −10.26475
2 −1.88731 −1.88730
1 0 0
0.5 0.71330 0.71329
0.2 1.05113 1.05113
0.1 1.14656 1.14656
0 1.232588 1.23259

Figures 2 and 3 show the skin friction coefficient f′′(0) against buoyancy parameter λ for some values of velocity exponent parameter m when velocity ratio parameter is ε=0.5 and ε=1. Two branches of solutions are found. The solid lines are the upper branch solutions and the dash lines are the lower branch solutions. With increasing m, the range of λ for which the solution exists increases. Also from both figures of the upper branch solutions, the skin friction is higher for the assisting flow (λ>0) compared to the opposing flow (λ<0). This implies that increasing the buoyancy parameter λ increases the skin friction coefficient f′′(0) whilst for ε=1, the values of f′′(0) as shown in Figure 3 are positive for λ>0 and negatives for λ<0. Physically, this means that positive f′′(0) implies the fluid exerts a drag force on the sheet and negative implies the reverse. Similarly this also happens for ε=0.5 but at different values of λ.

Variation of the skin friction coefficient f′′(0) with λ for various values of m when ε=0.5.

Variation of the skin friction coefficient f′′(0) with λ for various values of m when ε=1.

As seen in Figures 2 and 3, there exists a critical value of velocity ratio λc such that for λ<λc there will be no solutions, for λc<λ<0 there will be dual solutions, and when λ>0, the solution is unique. Our numerical computations presented in Figure 2 show that for the velocity ratio ε=0.5, λc=-8.331, −2.677, and −0.7411 for m = 2, 1, and 0.5, respectively. On the other hand, for the velocity ratio ε=1 shown in Figure 3, λc=-14.98, −4.764, and −1.301 for m = 2, 1, and 0.5, respectively. The dual solutions exhibit the normal forward flow behavior and also the reverse flow where f(η)<0. From these two results, it seems that an increase in velocity ratio parameter ε leads to an increase of the critical values of |λc|. This increases the dual solutions range of (6)-(7).

Figures 4 and 5 display the variations of the wall temperature θ(0) against the buoyancy parameter λ, for some values of m when ε=0.5 and ε=1, respectively. Both figures clearly show that the wall temperature increases as m increases for the upper branch solutions. For the lower branch solutions, the wall temperature becomes unbounded as λ0-.

Variation of the wall temperature θ(0) with λ for various values of m when ε=0.5.

Variation of the wall temperature θ(0) with λ for various values of m when ε=1.

The velocity and temperature profiles for ε=0.5 when λ=-0.5 are presented in Figures 6 and 7, respectively. Figure 6 shows that the velocity increases as m increases for the upper branch solutions, while opposite trend is observed for the lower branch solutions. In Figure 7, for the upper branch solutions, it is seen that an increase in m tends to decrease the temperature. Besides that, the temperature is higher for the lower branch solution than the upper branch solution at all points, near and away from the solid surface. It can be seen from Figures 6 and 7 that all profiles approach the far field boundary conditions (7) asymptotically, thus supporting the numerical results obtained. Finally, the typical streamlines for m=1,ε=1, and   λ=-2 for both solution branches are shown in Figures 8 and 9.

Velocity profile f(η) for various values of m when ε=0.5 and λ=-0.5.

Temperature profile θ(η) for various values of m when ε=0.5 and λ=-0.5.

Streamlines for the upper branch solutions when m=1,ε=1, and λ=-2.

Streamlines for the lower branch solutions when m=1,ε=1, and λ=-2.

4. Conclusions

The problem of mixed convection stagnation-point flow towards a nonlinearly stretching vertical sheet immersed in an incompressible viscous fluid was investigated numerically. The effects of the velocity exponent parameter m, buoyancy parameter λ, and velocity ratio parameter ε on the fluid flow and heat transfer characteristics were discussed. It was found that for the assisting flow, the solution is unique, while dual solutions were found to exist for the opposing flow up to a certain critical value λc. Moreover, increasing the velocity exponent parameter m is to increase the range of the buoyancy parameter λ for which the solution exists.

Nomenclature a , b , c :

Constants

f :

Dimensionless stream function

g :

Acceleration due to gravity

G r x :

Local Grashof number

k :

Thermal conductivity

m :

Velocity exponent parameter

Pr :

Prandtl number

q w :

Surface heat flux

R e x :

Local Reynolds number

T :

Fluid temperature

T :

Ambient temperature

u , v :

Velocity components along the x- and y-directions, respectively

U :

Free stream velocity

U w :

Stretching velocity

x , y :

Cartesian coordinates along the surface and normal to it, respectively.

Greek Letters α :

Thermal diffusivity

β :

Thermal expansion coefficient

ε :

Velocity ratio parameter

η :

Similarity variable

θ :

Dimensionless temperature

λ :

Buoyancy parameter

ν :

Kinematic viscosity

ρ :

Fluid density

Ψ :

Stream function.

Subscripts w :

Condition at the solid surface

:

Condition far away from the solid surface.

Superscript :

Differentiation  with respect to

Acknowledgments

The authors wish to thank the anonymous reviewers for their valuable comments and suggestions. The third author would like to acknowledge the financial supports received from the Universiti Kebangsaan Malaysia under the incentive grants UKM-GUP-2011-202 and DIP-2012-31.

Karwe M. V. Jaluria Y. Fluid flow and mixed convection transport from a moving plate in rolling and extrusion processes Journal of Heat Transfer 1988 110 3 655 661 2-s2.0-0024063986 Karwe M. V. Jaluria Y. Numerical simulation of thermal transport associated with a continuously moving flat sheet in materials processing ASME Journal of Heat Transfer 1991 113 3 612 619 2-s2.0-0026202548 Crane L. J. Flow past a stretching plate Journal of Applied Mathematics and Physics 1970 21 4 645 647 2-s2.0-33646099212 10.1007/BF01587695 Weidman P. D. Magyari Generalized Crane flow induced by continuous surfaces stretching with arbitrary velocities Acta Mechanica 2010 209 3-4 353 362 2-s2.0-72849120567 10.1007/s00707-009-0186-z Weidman P. D. Ali M. E. Aligned and nonaligned radial stagnation flow on a stretching cylinder European Journal of Mechanics B/Fluids 2011 30 1 120 128 2-s2.0-78649934512 10.1016/j.euromechflu.2010.08.001 Chen C. H. Laminar mixed convection adjacent to vertical, continuously stretching sheets Heat and Mass Transfe 1998 33 5-6 471 476 2-s2.0-0032047482 Ali M. Al-Yousef F. Laminar mixed convection boundary layers induced by a linearly stretching permeable surface International Journal of Heat and Mass Transfer 2002 45 21 4241 4250 2-s2.0-0036779220 10.1016/S0017-9310(02)00142-4 Ishak A. Nazar R. Pop I. Mixed convection boundary layers in the stagnation-point flow toward a stretching vertical sheet Meccanica 2006 41 5 509 518 2-s2.0-34848923849 10.1007/s11012-006-0009-4 Ishak A. Nazar R. Pop I. Hydromagnetic flow and heat transfer adjacent to a stretching vertical sheet Heat and Mass Transfer 2008 44 8 921 927 2-s2.0-43149103424 10.1007/s00231-007-0322-z Ali F. M. Nazar R. Arifin N. M. Pop I. Effect of Hall current on MHD mixed convection boundary layer flow over a stretched vertical flat plate Meccanica 2011 46 5 1103 1112 2-s2.0-79961097988 MR2822386 10.1007/s11012-010-9371-3 Pal D. Heat and mass transfer in stagnation-point flow towards a stretching surface in the presence of buoyancy force and thermal radiation Meccanica 2009 44 2 145 158 2-s2.0-63049090053 10.1007/s11012-008-9155-1 Ramachandran N. Chen T. S. Armaly B. F. Mixed convection in stagnation flows adjacent to vertical surfaces ASME Journal of Heat Transfer 1988 110 2 373 377 2-s2.0-0024016105 Ayub M. Zaman H. Sajid M. Hayat T. Analytical solution of stagnation-point flow of a viscoelastic fluid towards a stretching surface Communications in Nonlinear Science and Numerical Simulation 2008 13 9 1822 1835 2-s2.0-41149138549 10.1016/j.cnsns.2007.04.021 MR2406957 ZBL1221.76027 Bachok N. Ishak A. Pop I. On the stagnation-point flow towards a stretching sheet with homogeneous-heterogeneous reactions effects Communications in Nonlinear Science and Numerical Simulation 2011 16 11 4296 4302 2-s2.0-79957916315 10.1016/j.cnsns.2011.01.008 Bachok N. Ishak A. Pop I. Boundary layer stagnation-point flow and heat transfer over an exponentially stretching/shrinking sheet in a nanofluid International Journal of Heat and Mass Transfer 2012 55 8122 8128 10.1016/j.ijheatmasstransfer.2012.08.051 Bachok N. Ishak A. Pop I. Melting heat transfer in boundary layer stagnation-point flow towards a stretching/shrinking sheet Physics Letters A 2010 374 40 4075 4079 2-s2.0-77956266234 10.1016/j.physleta.2010.08.032 Bhattacharyya K. Dual solutions in boundary layer stagnation-point flow and mass transfer with chemical reaction past a stretching/shrinking sheet International Communications in Heat and Mass Transfer 2011 38 7 917 922 2-s2.0-79959934401 10.1016/j.icheatmasstransfer.2011.04.020 Hayat T. Abbas Z. Sajid M. MHD stagnation-point flow of an upper-convected Maxwell fluid over a stretching surface Chaos, Solitons and Fractals 2009 39 2 840 848 2-s2.0-62549147533 10.1016/j.chaos.2007.01.067 Ibrahim W. Shankar B. Nandeppanavar M. M. MHD stagnation point flow and heat transfer due to nanofluid towards a stretching sheet International Journal of Heat and Mass Transfer 2013 56 1 9 10.1016/j.ijheatmasstransfer.2012.08.034 Ishak A. Jafar K. Nazar R. Pop I. MHD stagnation point flow towards a stretching sheet Physica A 2009 388 17 3377 3383 2-s2.0-67349124437 10.1016/j.physa.2009.05.026 Labropulu F. Li D. Pop I. Non-orthogonal stagnation-point flow towards a stretching surface in a non-Newtonian fluid with heat transfer International Journal of Thermal Sciences 2010 49 6 1042 1050 2-s2.0-77950300211 10.1016/j.ijthermalsci.2009.12.005 Layek G. C. Mukhopadhyay S. Samad Sk. A. 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 2007 34 3 347 356 2-s2.0-33847623306 10.1016/j.icheatmasstransfer.2006.11.011 Yian L. Y. Amin N. Pop I. Mixed convection flow near a non-orthogonal stagnation point towards a stretching vertical plate International Journal of Heat and Mass Transfer 2007 50 23-24 4855 4863 2-s2.0-34848906837 10.1016/j.ijheatmasstransfer.2007.02.034 Mahapatra T. R. Dholey S. Gupta A. S. Oblique stagnation-point flow of an incompressible visco-elastic fluid towards a stretching surface International Journal of Non-Linear Mechanics 2007 42 3 484 499 2-s2.0-34249800352 10.1016/j.ijnonlinmec.2007.01.008 Paullet J. Weidman P. Analysis of stagnation point flow toward a stretching sheet International Journal of Non-Linear Mechanics 2007 42 9 1084 1091 2-s2.0-35648972784 10.1016/j.ijnonlinmec.2007.06.003 Rosali H. Ishak A. Pop I. Stagnation point flow and heat transfer over a stretching/shrinking sheet in a porous medium International Communications in Heat and Mass Transfer 2011 38 8 1029 1032 2-s2.0-80052168819 10.1016/j.icheatmasstransfer.2011.04.031 Yacob N. A. Ishak A. Pop I. Melting heat transfer in boundary layer stagnation-point flow towards a stretching/shrinking sheet in a micropolar fluid Computers and Fluids 2011 47 1 16 21 2-s2.0-79957619935 10.1016/j.compfluid.2011.01.040 MR2949391 Merkin J. H. Mahmood T. Mixed convection boundary layer similarity solutions: prescribed wall heat flux Journal of Applied Mathematics and Physics 1989 40 1 51 68 2-s2.0-34249969299 10.1007/BF00945309 MR980464 ZBL0667.76126 Wang C. Y. Stagnation flow towards a shrinking sheet International Journal of Non-Linear Mechanics 2008 43 5 377 382 2-s2.0-43249115375 10.1016/j.ijnonlinmec.2007.12.021 Keller H. B. Bramble J. A new difference scheme for parabolic problems Numerical Solution of Partial-Differential Equations 1970 2 New York, NY, USA Academic Press 327 350 MR0277129 ZBL0243.65060