# Unsteady Stagnation Point Flow and Heat Transfer over a Stretching/Shrinking Sheet with Suction or Injection

Authors: M. Suali; N. M. A. Nik Long; N. M. Ariffin
Journal: Journal of Applied Mathematics
 (2012)

---

## Abstract

The unsteady stagnation point flow and heat transfer over a stretching/shrinking sheet with suction/injection is studied. The governing partial differential equations are converted into nonlinear ordinary differential equations using a similarity transformation and solved numerically. Both stretching and shrinking cases are considered. Results for the skin friction coefficient, local Nusselt number, velocity, and temperature profiles are presented for different values of the governing parameters. It is found that the dual solutions exist for the shrinking case, whereas the solution is unique for the stretching case. Numerical results show that the range of dual solutions increases with mass suction and decreases with mass injection.


---

## Body

## 1. Introduction

The study of flow and heat transfer over a stretching/shrinking sheet receives considerable attention from many researchers due to its variety of application in industries such as extrusion of a polymer in a melt spinning process, manufacturing plastic films, wire drawing, hot rolling, and glass fiber production. Sakiadis [1, 2] reported the flow field analysis where the stretched surface was assumed to move with uniform velocity, and similarity solutions were obtained for the governing equations. Crane [3] gave a closed-form solution for steady two-dimensional flow of an incompressible viscous fluid caused by the stretching of an elastic sheet, which moves in its own plane with a velocity which varies linearly with distance from a fixed point. Chiam [4] investigated the steady two-dimensional stagnation point flow towards a stretching surface in the case when the stretching velocity is identical to the free stream velocity. P.S. Gupta and A.S. Gupta [5] extended the work of Crane [3] by investigating the effect of mass transfer on a stretching sheet with suction or blowing for linear surface velocity subject to uniform temperature. Mahapatra and Gupta [6] investigated the steady two-dimensional magnetohydrodynamics (MHDs) stagnation point flow of an incompressible viscous electrically conducting fluid toward a stretching surface, the flow being permeated by a uniform transverse magnetic field. In studying similar problem to Chiam [4], Mahapatra and Gupta [7] observed that the structure of the boundary layer depends on the stretching sheet parameters and the angle of incidence. Mahapatra and Gupta [8] obtained the exact similarity solution of the Navier-Stokes equations by considering the combination of both stagnation flow and stretching surface and observed that the flow displays a boundary-layer structure when the stretching velocity of the surface is less than the free stream velocity. Lok et al. [9] extended the Mahapatra and Gupta’s work [8] to the oblique stagnation flow and found that the free stream obliqueness is the shift of the stagnation point toward the incoming flow and depends on the inclination angle. Miklavcic and Wang [10] obtained the solution for steady two-dimensional, as well as axisymmetric viscous, flow over a shrinking sheet. Fang and Zhang [11] obtained the closed-form solution for steady MHD flow over a shrinking surface subject to applied suction. Wang [12] investigated the steady two-dimensional and axisymmetric stagnation point flow with heat transfer over a shrinking sheet and found that solutions do not exist for the larger shrinking rates. Ishak et al. [13] extended the work of Mahapatra and Gupta [8] to the oblique stagnation flow and a shrinking sheet in a micropolar fluid. It was found in these problems that the solutions for a shrinking sheet are non-unique. Bachok et al. [14] considered the similarity solutions of stagnation point flow and heat transfer over a warm, laminar liquid flow to a melting stretching or shrinking sheet. Recently, Patel and Timol [15] analyzed the steady two-dimensional stagnation point flow and heat transfer of a non-Newtonian fluid. In the unsteady flow problem, Wang [16] studied the concept of the flow and investigated the behavior of liquid film on an unsteady stretching surface. Andersson et al. [17] presented a new similarity solution for the temperature field on an unsteady stretching surface by transforming the time-dependent thermal energy equation into an ordinary differential equation. Elbashbeshy and Bazid [18] presented an exact similarity solution for unsteady flow and heat transfer over a horizontal stretching surface and studied the effect of various governing parameters such as the Prandtl number and unsteadiness parameter which determine the velocity and temperature profiles and heat transfer coefficient. Later, Fang [19] extended the flow over a shrinking sheet to power law surface velocity and obtained the multiple solutions for certain mass transfer with controlling parameters. Tsai et al. [20] studied numerically the nonuniform heat source or sink effect on the flow and heat transfer from an unsteady stretching sheet through a quiescent fluid. The results showed that the heat transfer rate and the skin friction increase as the unsteadiness parameter increases, but, as the space- and temperature-dependent parameters for heat source/sink increase, the heat transfer rate and the skin friction decrease. Sajid and Hayat [21] investigated the effect of the MHD for two-dimensional and axisymmetric shrinking sheet. Later, Fang et al. [22] found that multiple solutions exist for a certain range of mass suction and unsteadiness parameters for the unsteady viscous flow over a continuously shrinking sheet. Recently, Adhikary and Misra [23] investigated the unsteady two-dimensional hydromagnetic flow and heat transfer of an incompressible viscous fluid. Ali et al. [24] presented the influence of rotation, unsteadiness, and mass suction parameters on the reduced skin friction coefficients as well as the lateral velocity and velocity profiles for the unsteady viscous flow over a shrinking sheet with mass transfer in a rotating fluid. Recently, Fang et al. [25] investigated the boundary layers of an unsteady incompressible stagnation point flow with mass transfer whereas Bhattacharyya [26] found that the dual solutions for velocity distribution exist for certain value of velocity ratio parameter for the dual solutions in unsteady stagnation point flow over a shrinking sheet. Very recently, Nik Long et al. [27] found that multiple solutions exist for a certain range of ratio of the shrinking velocity to the free stream velocity which again depends on the unsteadiness parameter for the unsteady stagnation point flow and heat transfer over a stretching/shrinking sheet. The aim of this work is to extend the paper by Nik Long et al. [27] to the case of unsteady stagnation point flow of an incompressible viscous fluid by considering both stretching and shrinking sheets with suction or injection. The momentum and heat equation are solved numerically, and the characteristics of the flow are obtained.


## 2. Mathematical Formulation

Consider the unsteady stagnation point flow over a stretching or shrinking sheet immersed in an incompressible viscous fluid of ambient temperatureT∞. It is assumed that the free stream velocity is in the form U∞(x,t)=ax(1-λt)-1, the sheet is stretched with velocity Uw(x,t)=bx(1-λt)-1, and the surface temperature is Tw(x,t)=T∞+cx(1-λt)-1. The x-axis runs along the sheet, and the y-axis is measured normal to it. Having these assumptions, along with the boundary-layer approximations and neglecting the viscous dissipation, the governing equations are given by
(2.1)∂u∂x+∂v∂y=0,(2.2)∂u∂t+u∂u∂x+v∂u∂y=∂U∞∂t+U∞∂U∞∂x+ν∂2u∂y2,∂T∂t+u∂T∂x+v∂T∂y=α∂2T∂y2,
with the boundary conditions(2.3)u=Uw,v=Vw,T=Tw,aty=0,u⟶U∞,T⟶T∞,asy⟶∞,
where u and v are velocity components in x and y directions; respectively, ν is kinematic viscosity; α is thermal diffusivity; T is fluid temperature. Introducing the following similarity transformations(2.4)ψ(x,y,t)=(aν1-λt)1/2xf(η),η(y,t)=(aν(1-λt))1/2y,θ(η)=T-T∞Tw-T∞,
where η is similarity variable and ψ is stream function defined as u=∂ψ/∂y and v=-∂ψ/∂x, thus we have(2.5)u=ax1-λtf′(η),v=-(aν1-λt)1/2f(η);
therefore, the mass transfer velocity Vw can take the form
(2.6)Vw(t)=-(aν1-λt)1/2f0,
where prime denotes differentiation with respect to η. With these values of u and v, (2.1) is identically satisfied whereas (2.2) reduce to the following nonlinear ordinary differential equations(2.7)f′′′+ff′′+1-f′2+A(1-f′-12ηf′′)=0,1Prθ′′+fθ′-f′θ-A(θ+η2θ′)=0.
The boundary conditions (2.3) become
(2.8)f(0)=f0,f′(0)=ba=ε,θ(0)=1,θ(∞)⟶0,f′(∞)⟶1,
where ε(=b/a) is the ratio of stretching/shrinking velocity parameter and free stream velocity parameter, f0>0 and f0<0 are the suction and injection parameters, respectively, Pr=v/α is the Prandtl number, and A=λ/a is unsteadiness parameter. The quantities of physical interest to be obtained are the effect of suction or injection on the skin friction coefficient f′′(0) and the local Nusselt number –θ′(0) which are defined as
(2.9)Cf=τwρU∞2/2,Nux=xqwk(Tw-T∞),
where the surface shear stress τw and the surface heat flux qw are defined as(2.10)τw=μ(∂u∂y)y=0,qw=-k(∂T∂y)y=0,
with μ and k being the dynamic viscosity and thermal conductivity, respectively. Using the similarity variables equation (2.4), we obtain
(2.11)12CfRex1/2=f′′(0),NuxRex1/2=-θ′(0),
where Rex=U∞x/ν is local Reynolds number. It should be noticed that (2.7) reduce to those of Wang [12] and Nik Long et al. [27] when A=0 (steady-state flow) and A=0.01, respectively.


## 3. Results and Discussion

The nonlinear ordinary differential equations (2.7) subject to the boundary conditions (2.8) are solved numerically using the shooting method where we convert the boundary value problem into an initial value problem. The numerical results are given to carry out a parametric study showing the influence of the nondimensional parameters: the unsteadiness parameters A, ε, Prandtl number Pr, and f0. For the validation of numerical results, the case A=0 and Pr=0.7 with no effect of suction or injection (f0=0) are considered and compared to those of Wang [12]. The quantitative comparison is shown in Table 1 for the variation of ε with parameter f0 and found to be in favorable agreement. Based on our computations, the critical values of ε, εc with parameter f0 are presented in Table 2 and compare well with the previous results reported by Nik Long et al. [27]. It is observed that the values of |εc| increase as f0 increases. Hence, suction delays the boundary layer separation, whereas injection accelerates it.Table 1
Variation off′′(0) with ε when f0=0.
AεWang [12]PresentFirst solutionSecond solutionFirst solutionSecond solution03−4.2765450.21.051131.0511300.11.146561.146561−0.11.308601.308602−0.251.402241.402240−0.51.495671.495672−1.01.3288201.3288170−1.151.082230.1167021.0822320.116702Table 2
The values ofεc for different values of A and f0.
Af0Nik Long et al. [27]Present−0.3−0.265432−20−0.4393980.3−0.638668−0.3−0.816678−0.50−0.9565660.3−1.1396570.010−1.2536−1.253600−0.3−1.4641100.50−1.6011110.3−1.725011−0.3−1.96992220−2.1389740.3−2.383360The solution domains forf′′(0) are shown in Figures 1 and 2 for different mass transfer parameters f0 and ε, respectively. It is seen that there exist dual solutions for Ac<A<0 where Ac is a critical number dependent on f0. The solution is unique for (A>0)∪(A=Ac), and there is no solution for (A<Ac). The existence of the dual solution in our case is in agreement with [25, 28, 29]. The effect of ε and the mass suction (f0>0) enlarge the solution domain, and the mass injection (f0<0) reduces it. The velocity gradient decreases with decreasing values of A for the first and second solutions from positive to negative values of f′′(0), which indicates that there exist multiple regions for the momentum boundary layer. Figures 3 and 4 show the relationship between the skin friction coefficient and the local Nusselt number with ε and for the various values of f0. It is observed that there is no solutions for ε<εc, the dual solutions exist for εc<ε≤-1.0, and the solutions are unique for ε>-1.0, which means that the boundary layer separates from the surface. Figures 5 and 6 show the skin friction coefficient and the local Nusselt number with f0, respectively, for ε=-1.18 and A=0.01,0.04. For these values of ε, critical value (f0,c) of f0 being negative at which, there is a saddle node bifurcation. The dual solutions exist beyond the critical number which is f0>f0,c, but there is no solution for f0<f0,c. This behavior of solution indicates that the solutions exist for a certain interval of ε the case of injection, whereas no such interval appears for the suction, with both branches of solutions continuing to large value of f0. The skin friction coefficient f′′(0) and the local Nusselt number –θ′(0) increase as the parameter A increases and decrease as f0 increases. The samples of velocity and temperature profiles are presented in Figures 7 and 8, respectively, for different values of f0, ε=-1.18, and A=0.01. These figures show that the boundary conditions (2.8) for (2.7) are satisfied and approached infinity asymptotically, which support the numerical results presented in Figures 3, 4, 5, and 6 and exhibit the existence of the dual solution. We observed that the temperature profiles increase as the values of f0 increase. The same behavior of solution is observed for the velocity profiles.Figure 1
Solution domain for momentum boundary layers under different mass transfer parameters whenε=0.01.Figure 2
Solution domain for momentum boundary layers for different values ofε when f0=0.Figure 3
Skin friction coefficientf′′(0) as a function of ε for different values of f0 when A=0.01.Figure 4
Local Nusselt number–θ′(0) as a function of ε for different values of f0 when A=0.01.Figure 5
Skin friction coefficientf′′(0) as a function of f0 for different values of A.Figure 6
Local Nusselt number–θ′(0) as a function of f0 for different values of A when Pr=0.7.Figure 7
The velocity profilesf′(η) for different values of f0 when A=0.01, ε=-1.18.Figure 8
The temperature profilesθ(η) for different values of f0 when A=0.01, ε=-1.18, and Pr=0.7.


## 4. Conclusion

A numerical study is performed for the problem of unsteady two-dimensional stagnation point flow and heat transfer over a stretching and shrinking sheets with suction or injection. The similarity transformation is used to reduce the partial differential equations into nonlinear ordinary differential equations. The effects of the mass suction or mass injection parametersf0, governing parameters A, and ε on the fluid flow and heat transfer characteristics have been discussed, and the numerical results obtained are comparable well with the previously reported cases. It was found that the suction (f0>0) widens the range of ε for which solution exists, whereas the injection (f0<0) acts in the opposite manner. Moreover, for the shrinking sheet, the dual solutions exist whereas for the stretching sheet the solution is found to be unique for all ε.



---

*Source: 781845-2012-06-05.xml*