Mathematical Problems in Engineering

Volume 2013 (2013), Article ID 195360, 10 pages

http://dx.doi.org/10.1155/2013/195360

Research Article

Mixed Convection Boundary Layers with Prescribed Temperature in the Unsteady Stagnation Point Flow toward a Stretching Vertical Sheet

Nik Mohd Asri Nik Long,1,2 Lee Feng Koo,2,3 Tze Jin Wong,2,3 and Melini Suali¹

¹Department of Mathematics, Universiti Putra Malaysia, 43400 Serdang, Selangor, Malaysia

²Institute for Mathematical Research, Universiti Putra Malaysia, 43400 Serdang, Selangor, Malaysia

³Department of Basic Science and Engineering, Faculty of Agriculture and Food Sciences, Universiti Putra Malaysia, Sarawak Campus, 97008 Bintulu, Sarawak, Malaysia

Received 26 February 2013; Accepted 26 April 2013

Academic Editor: Mohamed Abd El Aziz

Copyright © 2013 Nik Mohd Asri Nik Long 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.

Abstract.

Mixed convection boundary layer caused by time-dependent velocity and the surface temperature in the two-dimensional unsteady stagnation point flow of an in-compressible viscous fluid over a stretching vertical sheet is studied. The transformed nonlinear boundary layer equations are solved numerically using the shooting technique in cooperation with Runge-Kutta-Fehlberg (RKF) method. Different step sizes are used ranging from 0.0001 to 1. Numerical results for the skin friction coefficient and local Nusselt number are presented for both assisting and opposing flows. It is found that the dual solutions exist for the opposing flow, whereas the solution is unique for the assisting flow. Important features of the flow characteristics are displayed graphically. Comparison with the existing results for the steady case show an excellent agreement.

1. Introduction

Stagnation point flow, describing the fluid motion near the stagnation region of a solid surface, exists in both cases of fixed or moving body in a fluid. Flow in incompressible viscous fluid on stretching sheet plays an important role in engineering applications such as extrusion of polymers, continuous casting, cooling of metallic plates, glass fiber production, hot rolling, paper production, wire drawing, and aerodynamic extrusion of plastic sheets [1–6]. Sakiadis [7] studied the boundary layer flow problem on a continuously moving surface with constant velocity. Crane [8] obtained the exact analytical solution for the two-dimensional stretching surface in a quiescent fluid. The effects of a variable surface temperature and linear surface stretching were examined by Grubka and Bobba [9]. Similar problem to that of Grubka and Bobba [9] for the case of power-law surface velocity and three different thermal boundary conditions were considered by Ali [10] and extended to the stretching surface subject to suction or injection [11]. Meanwhile, Ingham [12] investigated the existence of the dual solutions of the boundary layer equations of a continuously moving vertical plate with temperature inversely proportional to the distance up to the plate. Wang [13, 14] studied the behavior of liquid film on an unsteady stretching surface where the similarity transformation was employed to transform the governing partial differential equations to a nonlinear ordinary differential equation with an unsteadiness parameter. The boundary layer flow due to a stretching surface in vertical direction in steady, viscous, and incompressible fluid whereby the buoyancy forces are taken into consideration was discussed by Daskalakis [15], Lin and Chen [16], Ali and Al-Yousef [17], Chamkha [18], and Ishak et al. [19]. Vajravelu [20] studied the flow and heat transfer characteristics in a viscous fluid over a nonlinearly stretching sheet. Elbashbeshy and Bazid [21, 22] analyzed the heat transfer over an unsteady stretching surface. Tsai et al. [23] studied the nonuniform heat source or sink effect on the flow and heat transfer from an unsteady stretching sheet through a quiescent fluid numerically. Ishak and Nazar [24] extended Grubka and Bobba [9] and Ali’s [10] work to the axisymmetric laminar boundary layer flow along a continuous stretching cylinder immersed in an incompressible viscous fluid. Patil et al. [25] obtained the numerical results for two-dimensional mixed convection flow along a vertical semi-infinite power-law stretching sheet. Bachok et al. [26] adopted the RKF method to solve the free convection boundary layer flow near a continuously moving vertical plate subject to suction and injection, whereas Bachok and Ishak [27] discussed the variations of skin friction coefficient and heat transfer rate at the surface subject to the nonlinear stretching/shrinking sheet immersed in a viscous fluid. Aziz and Hashim [28] investigated the effects of viscous dissipation on flow and heat transfer in a thin liquid film on an unsteady stretching sheet using homotopy analysis method. Chamkha and Ahmed [29] presented the effects of heat generation/absorption and chemical reaction on unsteady magnetohydrodynamics flow heat and mass transfer near a stagnation point of a three-dimensional porous body in the presence of a uniform magnetic field. Suali et al. [30] studied the behavior of unsteady stagnation point flow of an incompressible viscous fluid by considering both stretching and shrinking sheets with suction and injection. Very recently, Liu and Megahed [31] analyzed the effects of internal heat generation on the flow and heat transfer on an unsteady stretching sheet in the presence of thermal radiation and variable heat flux by means of homotopy perturbation method.

In the present paper, we investigate the behavior of the mixed convection boundary layer flow with prescribed temperature in the unsteady stagnation point flow over a stretching vertical sheet. The momentum and energy equations are solved numerically, and the characteristics of the flow are analyzed.

2. Basic Equation

Consider two-dimensional incompressible viscous fluid near the unsteady stagnation point flow on a stretching vertical surface. Assume that the fluid occupies the half plane and the vertical surface located at the Cartesian coordinate with origin , and the -axis is along the direction of the surface and the -axis is normal to it. Two equal forces are impulsively applied along the -axis of the vertical stretching sheet where , , , and are denoted as the velocity of the stretching sheet, prescribed temperature, uniform temperature of the ambient fluid, and the free stream velocity to the boundary layer, respectively. Both assisting and opposing flows are taken into consideration, and their configuration is shown in Figure 1. Under this assumption together with the Boussinesq and boundary layer approximations, the governing equations can be written as follows:

where and are velocity components in and direction, , , , , and are gravity acceleration, fluid temperature, thermal diffusivity, kinematic viscosity, and thermal expansion coefficient, and the sign ± corresponds to the assisting and opposing buoyant flows, respectively, subjected to the boundary conditions Nik Long et al. [32]:

and , , and are denoted as

where , , and are positive constants. Equation (1) describes the unsteady flow over a stretching surface.

(a) Assisting flow

(b) Opposite flow

Figure 1: Physical model and coordinate system.

Introduce the following similarity transformations:

where is similarity variable and is stream function defined as

Substituting (4) into (5) yields

where the prime denotes the differentiation with respect to . Substituting (4) into (1) leads to the following nonlinear ordinary differential equations:

and the boundary condition (2) becomes

where is the ratio of stretching velocity parameter, is the unsteadiness parameter, is the Prandtl number and the constant , () is the buoyancy parameter defined as , is local Grashof number, and is the local Reynolds number. The buoyancy parameter, , can be written as

The solution of (7) subject to the boundary conditions (9) is given by

where and .

The quantities of physical interest are the skin friction coefficient, , and the Local Nusselt number, , which are defined by

where is the fluid density. The wall shear stress, , and the heat flux, , are defined as

where and are the dynamic viscosity and the thermal conductivity, respectively. The relation between and is defined as . Using the similarity variables (4), we obtain

which represent the skin friction coefficient and heat transfer rate at the surface, respectively.

3. Result and Discussion

The transformed nonlinear ordinary differential equations (7) subject to the boundary condition (9) are solved numerically using RKF method with shooting technique. The numerical results are given to carry out a parametric study showing the influence of the nondimensional parameters; that is, unsteadiness parameter, , Prandtl number, , and buoyancy parameter, . For the validation of the numerical results, the case of (steady flow) is considered and compared with those of Ishak et al. [33]. The quantitative comparison is tabulated in Table 1, and it is found to be in excellent agreement.

Table 1: Value of and for and at various .

(a)  

Buoyancy assisting flow

Pr

Ishak et al. [33]
Present
Ishak et al. [33]
Present

0.72
0.3645
0.3645
1.0931
1.0931

6.8
0.1804
0.1804
3.2902
−3.2896

20
0.1175
0.1175
5.6230
5.6201

40
0.0873
0.0872
7.9463
7.9383

60
0.0729
0.0728
9.7327
9.7180

80
0.0640
0.0639
11.2413
11.2187

100
0.0578
0.0577
12.5726
12.5411

(b)  

Buoyancy opposing flow

Pr

Ishak et al. [33]
Present
Ishak et al. [33]
Present

0.72
−0.3852
−0.3852
1.0293
1.0293

6.8
−0.1832
−0.1832
3.2466
3.2461

20
−0.1183
−0.1183
5.5923
5.5896

40
−0.0876
−0.0876
7.9227
−7.9149

60
−0.0731
−0.0730
9.7126
−9.6982

80
−0.0642
−0.0641
11.2235
11.2012

100
−0.0579
−0.0578
12.5564
12.5252

Figures 2 and 3 present the variations of skin friction coefficient, , and local Nusselt number, , with respect to buoyancy parameter, , for different values of , and . For , the variations are displayed in Figures 4 and 5. It is observed that for the assisting flow, an increase in would increase, and for all cases (Figures 2 to 5), while an increase in would decrease (Figures 2 and 4) and increase (Figures 3 and 5). For the opposing flow, as and increase the decreases and increases for the first solution, but opposite behavior can be seen for the second solution. It is also noted that for the opposing flow, the dual solutions exist for and no solutions for . For , there exist a unique solution; the boundary layer separation occurs.

(a) Assisting flow

(b) Opposing flow

Figure 2: Comparison of skin friction, for different when .

(a) Assisting flow

(b) Opposing flow

Figure 3: Comparison of local Nusselt number, for different when .

(a) Assisting flow

(b) Opposing flow

Figure 4: Comparison of skin friction, for different when .

(a) Assisting flow

(b) Opposing flow

Figure 5: Comparison of local Nusselt number, for different when .

The effects of unsteadiness parameter, and , on and for fix and are shown in Figures 6 and 7. For the assisting flow, as and increase, the decreases and increases. While for the opposing flow, as increases the and decrease for the first solution, and for the second solution, and decrease until they archive their minimum values, and and then they increase until . For , there exist no solution. An increase in would increase but only give little effect on for the first solution, and and decrease for the second solution.

(a) Assisting Flow

(b) Opposing Flow

Figure 6: Comparison of skin friction, for different when and .

(a) Assisting flow

(b) Opposing flow

Figure 7: Comparison of local Nusselt number, for different when and .

Figures 8 and 9 exhibit the variations of the skin friction coefficient, and local Nusselt number, for some value of . As and increase, the decreases and increases for the assisting flow. While for the opposing flow as increases and decreases for the first solution and increases and decreases for the second solution. For the increasing values of , decreases and increases for both first and second solutions.

(a) Assisting flow

(b) Opposing flow

Figure 8: Comparison of Skin friction, for different when and .

(a) Assisting flow

(b) Opposing flow

Figure 9: Comparison of local Nusselt number, for different when and .

The velocity and temperature profile for different values of and when and are displayed in Figures 10 and 11. The existance of the dual solutions for the opposing flow are clearly displayed. It is seen that for the assisting flow, the velocity profile increases to its maximum value, , then decreases monotonically to a constant value. For the opposing flow, the velocity profile shows the opposite trend. The temperature profile decreases monotonically with and becomes zero outside the boundary layer which satisfies the condition at infinity, . These properties support the validation of the present results. It is also noticed that an increase in values of and will increase the heat transfer rate at the surface. For the second solution, the temperature profile increases to attain its maximum value , then decreases gradually to a constant.

(a) Assisting flow

(b) Opposing flow

Figure 10: Velocity profile, for different when and .

(a) Assisting flow

(b) Opposing flow

Figure 11: Temperature profile, for different when and .

4. Conclusion

In this study, the numerical solutions for both assisting and opposing flows of the mixed convection boundary layer caused by time-dependent velocity and surface temperature in two-dimensional unsteady stagnation point flow over a stretching vertical sheet are obtained. The similar transformation is advocated to reduce the governing partial differential equations into nonlinear ordinary differential equations which are solved numerically using RKF method via the shooting technique. Comparison with the existing results showed the excellent agreement. Dual solutions are found to exist for the opposing flow whereas the solution is unique for the assisting flow. The boundary layer separation occurs at for the opposing flow. The effect of buoyancy parameter, , the unsteadiness parameter, , and Prandtl number, , of the fluid has been displayed graphically and discussed in details.

Acknowledgment

This project is supported by the University Putra Malaysia for the Research University Grant scheme, Project no. 05-02-12-1834RU.

References

E. M. A. Elbashbeshy, “Heat transfer over a stretching surface with variable surface heat flux,” Journal of Physics D, vol. 31, no. 16, pp. 1951–1954, 1998.

T. Fang, C. F. F. Lee, and J. Zhang, “The boundary layers of an unsteady incompressible stagnation-point flow with mass transfer,” International Journal of Non-Linear Mechanics, vol. 46, no. 7, pp. 942–948, 2011.

T. G. Fang, J. Zhang, and S. S. Yao, “Viscous flow over an unsteady shrinking sheet with mass transfer,” Chinese Physics Letters, vol. 26, no. 1, Article ID 014703, 4 pages, 2009.

A. Ishak, R. Nazar, and I. Pop, “Steady and unsteady boundary layers due to a stretching vertical sheet in a porous medium using darcy-brinkman equation model,” International Journal of Applied Mechanics and Engineering, vol. 11, no. 3, pp. 623–637, 2006.

A. Ishak, R. Nazar, and L. Pop, “Heat transfer over an unsteady stretching surface with prescribed heat flux,” Canadian Journal of Physics, vol. 86, no. 6, pp. 853–855, 2008.

Y. Y. Lok, A. Ishak, and I. Pop, “MHD stagnation-point flow towards a shrinking sheet,” International Journal of Numerical Methods for Heat and Fluid Flow, vol. 21, no. 1, pp. 61–72, 2011.

B. C. Sakiadis, “Boundary layer behavior on continuous solid surfaces,” The American Institute of Chemical Engineers Jounal, vol. 7, pp. 221–225, 1961.

L. J. Crane, “Flow past a stretching plate,” Zeitschrift für Angewandte Mathematik und Physik ZAMP, vol. 21, no. 4, pp. 645–647, 1970.

L. J. Grubka and K. M. Bobba, “Heat transfer characteristics of a continuous, stretching surface with variable temperature,” Journal of Heat Transfer, vol. 107, no. 1, pp. 248–250, 1985.

M. E. Ali, “Heat transfer characteristics of a continuous stretching surface,” Warme-Und Stoffubertragung, vol. 29, no. 4, pp. 227–234, 1994.

M. E. Ali, “On thermal boundary layer on a power-law stretched surface with suction or injection,” International Journal of Heat and Fluid Flow, vol. 16, no. 4, pp. 280–290, 1995.

D. B. Ingham, “Singular and non-unique solutions of the boundary-layer equations for the flow due to free convection near a continuously moving vertical plate,” ZAMP Zeitschrift für Angewandte Mathematik und Physik ZAMP, vol. 37, no. 4, pp. 559–572, 1986.

C. Wang, “Liquid film on an unsteady stretching sheet,” Quarterly of Applied Mathematics, vol. 48, pp. 601–610, 1990.

C. Wang, “Analytic solutions for a liquid film on an unsteady stretching surface,” Heat and Mass Transfer, vol. 42, no. 8, pp. 759–766, 2006.

J. E. Daskalakis, “Free convection effects in the boundary layer along a vertically stretching at surface,” Canadian Journal of Physics, vol. 70, no. 12, pp. 1587–1597, 1993.

C. R. Lin and C. K. Chen, “Exact solution of heat transfer from a stretching surface with variable heat flux,” Heat and Mass Transfer, vol. 33, no. 5-6, pp. 477–480, 1998.

M. Ali and F. Al-Yousef, “Laminar mixed convection from a continuously moving vertical surface with suction or injection,” Heat and Mass Transfer, vol. 33, no. 4, pp. 301–306, 1998.

A. J. Chamkha, “Hydromagnetic three-dimensional free convection on a vertical stretching surface with heat generation or absorption,” International Journal of Heat and Fluid Flow, vol. 20, no. 1, pp. 84–92, 1999.

A. Ishak, R. Nazar, and I. Pop, “Mixed convection on the stagnation point flow toward a vertical, continously stretching sheet,” Journal of Heat Transfer, vol. 129, no. 8, pp. 1087–1090, 2007.

K. Vajravelu, “Viscous flow over a nonlinearly stretching sheet,” Applied Mathematics and Computation, vol. 124, no. 3, pp. 281–288, 2001.

E. M. A. Elbashbeshy and M. A. A. Bazid, “Heat transfer over an unsteady stretching surface with internal heat generation,” Applied Mathematics and Computation, vol. 138, no. 2-3, pp. 239–245, 2003.

E. M. A. Elbashbeshy and M. A. A. Bazid, “Heat transfer over an unsteady stretching surface,” Heat and Mass Transfer, vol. 41, no. 1, pp. 1–4, 2004.

R. Tsai, K. H. Huang, and J. S. Huang, “Flow and heat transfer over an unsteady stretching surface with non-uniform heat source,” International Communications in Heat and Mass Transfer, vol. 35, no. 10, pp. 1340–1343, 2008.

A. Ishak and R. Nazar, “Laminar boundary layer flow along a stretching cylinder,” European Journal of Scientific Research, vol. 36, no. 1, pp. 22–29, 2009.

P. M. Patil, S. Roy, and A. J. Chamkha, “Mixed convection flow over a vertical power-law stretching sheet,” International Journal of Numerical Methods for Heat and Fluid Flow, vol. 20, no. 4, pp. 445–458, 2010.

N. Bachok, A. Ishak, R. Nazar, and I. Pop, “Ingham problem for free convection near a continuously moving vertical permeable plate,” IMA Journal of Applied Mathematics, vol. 77, no. 4, pp. 578–589, 2011.

N. Bachok and A. Ishak, “Similarity solutions for the stagnation-point flow and heat transfer over a nonlinearly stretching/shrinking sheet,” Sains Malaysiana, vol. 40, no. 11, pp. 1297–1300, 2011.

R. C. Aziz and I. Hashim, “Liquid film on unsteady stretching sheet with general surface temperature and viscous dissipation,” Chinese Physics Letters, vol. 27, no. 11, Article ID 110202, 4 pages, 2010.

A. J. Chamkha and S. E. Ahmed, “Similarity solution for unsteady MHD flow near a stagnation point of a three-dimensional porous body with heat and mass transfer, heat generation/absorption and chemical reaction,” Journal of Applied Fluid Mechanics, vol. 4, no. 2, pp. 87–94, 2011.

M. Suali, N. M. A. Nik Long, and N. M. Ariffin, “Unsteady stagnation point flow and heat transfer over a stretching/shrinking sheet with suction or injection,” Journal of Applied Mathematics, vol. 2012, Article ID 781845, 12 pages, 2012.

I. C. Liu and A. M. Megahed, “Homotopy perturbation method for thin film flow and heat transfer over an unsteady stretching sheet with internal heating and variable heat flux,” Journal of Applied Mathematics, vol. 2012, Article ID 418527, 12 pages, 2012.

N. M. A. Nik Long, M. Suali, A. Ishak, N. Bachok, and N. M. Arifin, “Unsteady stagnation point flow and heat transfer over a stretching/shrinking sheet,” Journal of Applied Sciences, vol. 11, no. 20, pp. 3520–3524, 2011.

A. Ishak, R. Nazar, and I. Pop, “Mixed convection boundary layers in the stagnation-point flow toward a stretching vertical sheet,” Mechanica, vol. 41, no. 5, pp. 509–518, 2006.