Unsteady Boundary Layer Flow past a Stretching Plate and Heat Transfer with Variable Thermal Conductivity

An unsteady boundary layer flow of viscous incompressible fluid over a stretching plate has been considered to solve heat flow problem with variable thermal conductivity. First, using similarity transformation, the velocity components have been obtained, and then the heat flow problem has been attempted in the following two ways: 1) prescribed stretching surface temperature (PST), and 2) prescribed stretching surface heat flux (PHF) Flow and temperature fields have been analyzed through graphs. The expressions for skin friction and coefficient of convective heat transfer Nusselt number in PST and PHF cases have been derived.


Introduction
Due to number of applications in industrial manufacturing process, the problem of boundary layer flow past a stretching plate has attracted considerable attention of researchers during the past few decades.Examples of such technological process are hot rolling, wire drawing, glass-fiber and paper production.In the process of drawing artificial fibers the polymer solution emerges from orifice with a speed which increases from almost zero at the orifice up to a plateau value at which remains constant.The moving fiber, which is of great technical importance, is governed by the rate at which the fiber is cooled and this, in turn affects the final properties of the yarn.A number of works are presently available that follow the pioneering classical work of Sakiadis [1], F. K. Tsou, E. M. Sparrow, R. J. Goldstein [2] and Crane [3].Table 1 lists some relevant works that pertain to cooling liquids, i.e., heat transfer for stretching surface.
There are liquid metals whose thermal conductivity varies with temperature in an approximately linear manner in the range from 0˚F to 400˚F.In 1996, T. C. Chiam [21] considered heat transfer problem with variable thermal conductivity in stagnation-point flow towards stretching sheet.Naseem Ahmad and Kavita Marwah [22] also studied boundary layer flow of Walters Liquid B Model with heat transfer for linear stretching plate with variable conductivity numerically.In 2010, Ahmad and Mishra investigated unsteady boundary layer flow and heat transfer over a stretching sheet [23].
In almost all the problems of stretching sheet with heat transfer where closed form solution is obtained, the thermal conductivity of liquid has been taken constant.The present paper is the extension of the work done by N. Ahmad and M. Mishra [23] assuming that the thermal conductivity varies in linear manner with temperature.The temperature field has two parts: one mean temperature, other is due to variable thermal conductivity.Both the parts mean temperature and temperature due to variable thermal conductivity have been analyzed thoroughly for some new recommendations.

Mathematical Formation and Solution
The problem considered here is the unsteady boundary layer flow due to a stretching flat plate in a quiescent viscous incompressible fluid.The flow is two dimensional where x-axis is along the plane of moving plate and y-axis is normal to it, respectively.We assume that the surface is moving continuously with the velocity ; , 1 and t < 1/a in the positive x-direction.Under these assumptions, the boundary layer flow along moving plate is governed by the equations where u, the horizontal velocity component; v, the vertical velocity component;  , the kinematic viscosity The relevant boundary conditions are: the Equations (2.1) and (2.2) reduce to with boundary conditions 0, , 0 , 0 where bar has been dropped for convenience.
Setting the similarity solution of the form Putting u and v in the Equation (2.4), we have and the relevant boundary conditions become Boundary conditions suggest that the velocity function f may be of the form   ry f y e    where r is unknown to be determined.Thus .Therefore, we have the velocity components as follows: and using continuity equation, we have Copyright © 2012 SciRes.WJM

Skin Friction
The wall shear stress at the stretching plate is given by 0 Thus, the skin friction is

Heat Transfer Problem
In absence of viscous dissipation and heat generation, the energy equation for two dimensional heat flow is given subject to boundary conditions 0, where T P is plate temperature, T  is temperature of surrounding fluid, C P is specific heat at constant pressure and k is thermal conductivity.

Case A: Prescribed Power Law Surface Temperature (PST)
Let the surface temperature be of the form 2 0, while the temperature outside the dynamic region be .Now, we define the dimensionless temperature by For liquid metals, it has been found that the thermal conductivity varies with temperature in an approximately linear manner in the range from 0˚F to 400˚F.Therefore, we assume k as . Now, substituting u and v in the Equation (4.1) and changing the independent variable y to ry   , we have The Equation (4.3) can be rewritten as From Equation (4.5), we note that the heat transfer takes place in two parts according to and 0  0  .If 0  , then we have the main heat transfer due to constant thermal conductivity i.e.
  and if 0  , then we get the first correction equation to main heat transfer as

Nusselt Number
The coefficient of convective heat transfer is given by 0 Therefore the Nusselt number is (see Table 2):

Case B: Prescribed Power Law Surface Heat
Flux (PHF Case) The power law heat flux on the surface of stretching plate is considered to be a quadratic power of x in the form 2 where D is a constant, k is the thermal conductivity.Now together with boundary conditions: Equating the terms independent of and the terms involving from Equation (4.16), we get the following two boundary value problems: and where the nomenclature m g is main heat transfer when thermal conductivity is constant and c g is correction to the heat flow due to variation in thermal conductivity in PHF case.The solution of the Equation (4.18) together with the boundary conditions (4.18) is The general solution of the Equation (4.19) is   where we again observe that the dynamic region for this temperature field is 0 1    .Hence, the boundary conditions (4.17  3) Recalling (*), we have a b e at    

Discussion and Results
The problem of unsteady boundary layer flow of viscous incompressible fluid overstretching plate has been analyzed.The velocity field has been obtained by similarity transformation method.Later, the heat flow problem has been studied by considering PST and PHF cases.We summarize the results as in Figure 1 to Figure 6. Figure 1 shows that horizontal component of velocity u.It increases as time progresses within the dynamical region [0, 1].We also see that u is maximum in the immediate neighborhood of stretching plate and it starts decreasing as .In fully developed flow, as time goes on progressing, the velocity progresses too.       is independent of Prandtl number P r .Figure 5 is the graph of ReC f versus time.We see that as time progresses the skin friction increases .We mean that as time progresses the velocity increases, in turn skin friction increases.
Figure 6 is temperature field in PHF case.We see that it approaches to zero asymptotically.
We see the g m for P r = 5.10 of liquid oxygen at 56˚K, P r = 2.15 of para-hydrogen at 14˚K and P r = 1.54 of liquid ammonia at 10˚C.It has been observed that g m increases absolutely as P r increases, but for P r = 9.42 of water at 10˚C, g m behaves in a different way due to its density.
Nusselt number in PST case is given by the Equation (4.12) where we observe referring Table 2 that it depends on Prandtl number and time t as well.From Table 1 we find that as Prandtl number increases the Nusselt number also increases but Nusselt number decreases as time increases.
Referring Table 3 for PHF case, we see the same pattern as in PST case from Table 2.
We have not got Prandtl number from the temperature field due to variation in thermal conductivity.

1 
gamma function Equation (4.8) is a non linear differential equation of order two.Let the solution of this equation be of the form:  .Therefore, the heat transfer in case 0  takes place within the dynamic region 0 1    .
Writing the Equation (4.2) in terms of   g  , we get the following differential equation

Figure 2 Figure 3
Figure 2 is a graph of v versus y for different instant of time.The vertical component v is almost constant within [0, 1] and later it starts increasing.v progresses as we march away from the slit and it also increases as time progresses.Figure 3 is to study the variation of mean temperature field m  with respect to Prandtl number P r .We see that as P r increases, m  decreases within dynamical region [0,1] in PST case.As r ul P v  decreases, i.e. kinematic viscosity  increases in turn viscosity of fluid increases.In case of more viscosity, generally flow of heat becomes slow.It is supported by our study.Figure 4 is expressing the trend of temperature field

Figure 4 Figure 1 .
Figure 1.Velocity field at different instant of time.

Figure 2 .
Figure 2. Velocity component v at different instant of time.

Figure 3 .
Figure 3. Mean temperature field for different values of Prandtl number P r in PST case.

Figure 4 .
Figure 4. Temperature distribution due to variation in thermal conductivity.

Figure 5 .
Figure 5. Skin friction with respect to time t.

Figure 6 .
Figure 6.Temperature field for different values of Prandtl number in PHF case.