ABSTRACT

In the present paper, the boundary layer flow of Walters Liquid B Model over a stretching plate has been considered to solve heat flow problem with variable conductivity. First, using similarity transformation, the velocity components have been obtained. Then, the heat flow problem has been considered in two ways: 1) prescribed surface temperature (PST), and 2) prescribed stretching plate heat flux (PHF) in case of variable conductivity. Due to variable conductivity, temperature profile has its two partone mean temperature and other temperature profile induced due to variable conductivity. The related results have been discussed with the help of graphs.

1. 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, glassfiber 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 it 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]. The following Table 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 [22] considered heat transfer problem with variable thermal conductivity in stagnation-point flow towards stretching sheet. N. Ahmad and K. Marwah [16,17] also studied boundary layer flow of Walters Liquid B Model with heat transfer for linear stretching plate with variable thermal conductivity numerically.

In almost all the flow problems over stretching sheet with heat transfer referred so far where closed form solution is obtained, the thermal conductivity of liquid has been taken constant. In this paper, we try to solve the boundary layer flow of Walters Liquid B Model over a stretching plate and heat transfer with variable thermal conductivity in the following two cases:

1) prescribed surface temperature (PST), and 2) prescribed stretching plate heat flux (PHF). The temperature field has been obtained in each case and the effect on it, of Visco-elasticity and Prandtl number has been discussed graphically.

2. Mathematical Formation and Solution

The problem considered here is the steady boundary layer flow due to a moving flat plate in a quiescent Walters Liquid B Model. 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 in the positive x-direction. Under these assumptions, the boundary layer along moving plate is governed by the equations:

(1)

(2)

where u, the horizontal velocity component; v, the vertical velocity component;, the kinematic viscosity; k₀, the coefficient of visco-elasticity The relevant boundary conditions are:

Introducing the dimensionless variables

the Equations (1) and (2) reduce to

(3)

(4)

where k₁=with boundary conditions

(5)

where bar has been dropped for convenience.

Setting the similarity solution of the form, we have

(6)

Putting u and v in the Equation (2), we have

(7)

which is non-linear differential equation of order three.

The boundary conditions (5) reduce to

(8)

(9)

Boundary conditions suggest that the velocity function may be of the form where r is complex number with positive real part. Thus,

Now, from the Equation (7), we get.

Therefore, the velocity components become as follows:

(8)

3. Heat Transfer Problem

In absence of viscous dissipation and heat generation, the energy equation for two dimensional heat flow is given by

(9)

subject to boundary conditions

(10)

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.

3.1. Case A: Prescribed Power Law Surface Temperature (PST)

Let the surface temperature be of the form

while the temperature out side the dynamic region be. Now, we define the dimensionless temperature by

where.

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 where.

Now, substituting u and v and changing the independent variable y to, we have

(11)

with boundary conditions

(12)

(13)

From Equation (11), we note that the heat transfer takes place in two parts, that is, one part of heat transfer is due to temperature difference and the other part is due to variable thermal conductivity. We denote the first part by q_m and second by q_v. Thus, equating the terms independent of e and the terms involving e, we have

(14)

(15)

and

(16)

, as (17)

The solution of the Equation (14) is

(18)

Equation (17) is a non-linear differential equation of order two. Let the solution of this equation be of the form

Putting this solution in the Equation (18), we have

(19)

The roots of this equation are 0 and 1/2. Therefore,

(20)

The general solution (20) of the Equation (16) is real only when. Therefore, the heat transfer due to variable thermal conductively takes place within the dynamic region. Hence, the boundary conditions (17) may be presented as

(21)

The solution (20) finally reduces to

(22)

3.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

at (23)

, as (24)

where D is a constant, k is the thermal conductivity. Now we define dimensionless temperature by

(25)

where

when

and.

Writing the Equation (11) in terms of, we get the following differential equation

(26)

together with boundary conditions:

(27)

Equating the terms independent of and the terms involving from Equation (26), we get the following two boundary value problems:

(28a)

, as (28b)

and,

(29a)

, as (29b)

The solution of the Equation (28a) together with boundary conditions (28b) is given by

(30)

The general solution of the Equation (29a) is

(31)

where C and D are arbitrary constants to be determined. Here, we observe that the dynamic region for this temperature field is as it was in PST case. Therefore, the boundary conditions may be taken as

(32)

Hence finally the solution becomes

(32)

4. Discussion and Results

A boundary layer flow of visco-elastic fluid (Walters Liquid B Model) over a stretching plate and heat transfer has been analyzed through out this paper. The boundary layer equations of momentum and heat transfer have been solved analytically. Two different analytical expressions have been obtained for dimensionless temperature fields for two general cases of boundary conditions namely 1) PST case, and 2) PHF respectively. The equation governing the induced temperature field due to variable thermal conductivity has been solved with modified boundary conditions in both the cases PST and PHF The related results have been discussed in the subsequent paragraphs in detail.

Figure 1 is a graph drawn for mean temperature field versus in PST case taking different values of visco-elastic parameter. It is observed that the temperature at the surface of stretching plate is invariant with respect to physical parameters. Further, we observe that as we move away the stretching plate within dynamic region, the temperature field increases as viscoelastic parameter increases. Physically, when increases, the fluid absorbs more heat which causes the increase of temperature.

Figure 2 exhibits variation pattern of mean temperature field q_m(h) with regard to Prandtl number Pr in PST case. We notice that as Prandtl number Pr increases, the mean temperature decreases.

Figure 3 exhibits the variation pattern of induced temperature field q_v in PST case. This temperature field is independent of k₀ and Pr. We see that the induced temperature decreases as we move away the plate. The variation pattern is represented by a parabolic curve having focus at h = 1. q_v decreases as we move away the stretching plate and it becomes almost zero at h = 1 which end upper edge of boundary layer thickness.

Figure 4 is the graph of mean temperature field versus in PHF case. Here we observe that this temperature field attains the maximum value at the surface of stretching plate which is more than PST case. The reason of this maximum value is that the direction of heat flow is from fluid to stretching plate. Temperature field is approaching to zero asymptotically. The variation pattern is almost same as in PST case.

Figure 5 shows the graph of versus for the different values of Prandtl number Pr in PHF case. It has been observed that as Pr increases, temperature field decreases absolutely. Temperature field is tending to zero asymptotically.

Equation (32) represents the temperature field induced by variable thermal conductivity in PHF case. The geometry of this equation is parabola having focus at h = 1. This temperature field is independent of k₁ and Pr. We see that the induced temperature decreases as we move away the plate. It varies from maximum at the plate to zero at the upper end of boundary layer.

Putting k₀ = 0, we get the results due to Naseem Ahmad, Z. U. Siddiqui and M. K. Mishra [23].

5. Conclusion

The boundary layer flow of Walters Liquid B over linear stretching plate has been studied together with heat transfer with variable thermal conductivity. The heat transfer has been studied by considering two cases: PST and PHS. The effect of visco-elasticity k₀ has been looked upon in the Figures 1 and 2. This paper is an extension of the work due to Naseem Ahmad, Z. U. Siddiqui and M. K. Mishra [23].

REFERENCES