Paper 19
Paper 19
Paper 19
' ' c c
'
' ' ' = +
' ' c c
(1)
( )
2
0 2
r
p h
q T T
C k Q T T
t y y
' ' c c c
'
' =
' ' ' c c c
(2)
With the following boundary conditions
0: 0, 0
0: , 0
0,
h
h
t u T T for y h
T q
t u U at y
y k
u T T at y h
' ' ' ' ' s = = s s
' c
' ' ' > = = =
' c
'
' ' ' = = =
(3)
where g is the acceleration due to gravity, the volumetric co-
efficient of thermal expansion, v the kinematic viscosity,
the density, k the thermal conductivity,
p
C the specific heat at
constant pressure, q the constant heat flux,
r
q the radiative
heat flux in y' direction, T' the fluid temperature, and u'
is the fluid velocity.
The radiative heat flux term is simplified by making use of the
Rosseland approximation Siegel and Howell [15] as
4
4
3
r
T
q
k y
o
-
' c
=
' c
(4)
Where o is the Stefan - Boltzmann constant and k
-
is the
mean absorption coefficient. It should be noted that by using
the Rosseland approximation we limit our analysis to optically
thick fluids. If temperature differences within the flow are
sufficiently small such that
4
T' may be expressed as a linear
function of the temperature, Then the Taylor series for
4
T'
about
h
T' , after neglecting higher order terms, is given by
4 3 4
4 3
h h
T T T T ' ' ' ' ~
(5)
It is emphasized here that equation (5) is widely used in
computational fluid dynamics involving radiation absorption
problems Chung [5] in expressing the term
4
T' as a linear
function.
In view of Equations (4) and (5), Equation (2) reduces to
3 2 2
2 2
16
3
h
p
T T T T
C k
t y k y y
o
-
' ' ' ' c c c c
= +
' ' ' ' c c c c
(6)
In order to solve the governing equations in dimensionless
form, we introduce the following non-dimensional quantities:
International Journal of Engineering Research (ISSN : 2319-6890)
Volume No.2, Issue No.2, pp : 113-118 01 April 2013
IJER@2013 Page 115
( )
3
2
2 2 2 2
0 0
3
, , , ,
/
, , Pr ,
4
h
p
p h
T T y t u g h q
y t u Gr
h h U hq k Uk
C
B h Q h h kk
M K R
k C k T
v |
u
v
o
|
v v o
-
'
' ' ' '
= = = = =
= = = = =
'
(7)
where Gr is the thermal Grashof number, Pr the Prandtl
number, R the radiation parameter, | heat source parameter, t
the dimensionless time, u the dimensionless velocity, y the
dimensionless coordinate axis normal to the plate, the
coefficient of viscosity and is the dimensionless temperature.
Then in view of Equations (7), Equations (1), (6) and (3)
reduces to the following non-dimensional form of equations:
2
2
1 u u
Mu u Gr
y t K
u
c c
+
c c
(8)
( )
2
2
3 4 3 Pr 3 Pr R R R
y t
u u
|u
c c
+
c c
(9)
The initial and boundary conditions are
0: 0, 0 0 1
0: 1, 1 0
0, 0 1
t u for y
t u at y
y
u at y
u
u
u
s = = s s
c
> = = =
c
= = =
(10)
3 SOLUTION OF THE PROBLEM
Equation (8) (9) are coupled, non linear partial differential
equations and these cannot be solved in closed form using
the initial and boundary conditions (10). However, these
equations can be reduced to a set of ordinary differential
equations, which can be solved analytically. This can be done
by representing the velocity, temperature and concentration of
the fluid in the neighbourhood of the fluid in the
neighbourhood of the plate as
( ) ( ) ( )
2
0 1
0
nt
f f y e f y c c = + + (11)
Substituting (11) in Equation (8) (9) and equating the
harmonic and non harmonic terms, and neglecting the higher
order terms of
( )
2
0 c , we obtain
0 1 0 0
u u Gr | u '' = (12);
1 4 1 5 1
u u | | u '' =
(13)
0 9 0
0 u | u '' = (14);
1 10 1
0 u | u '' = (15)
The corresponding boundary conditions can be written as
0 1
0 1
0 1 0 1
1, 0, 1, 0, 0
0, 0, 0, 0, 1
u u at y
y y
u u as y
u u
u u
c c
= = = = =
c c
(16)
The solutions of Equations (12) - (15) under the initial and
boundary conditions (16) by perturbation technique is given
by
( )
6 2 1 2
1 2 3 4
,
m y m y m y m y
u y t Z e Z e Z e Z e = + + +
( )
2 1
1 2
,
m y m y
y t De De u = +
APPENDIX
The constant not given Brevity of the space
4 RESULTS AND DISCUSSION
An exact solution to the problem of natural convection in
unsteady Couette flow between two long vertical parallel
plates in the presence of constant heat flux and thermal
radiation have been presented in the preceding section. In
order to get the physical insight into the problem, the
numerical values of the temperature field, the velocity field,
the skin-friction, the Nusselt number, the volume flow rate and
the vertical heat flux are computed for different values of the
system parameters such as Radiation parameter (R), Grashof
number
( ) Gr , Prandtl number (Pr) and heat source
parameter
( ) | . Figure (1) presents the velocity profiles for
both air and water
( ) Pr 7.0 = in the case of pure convection
( ) 100 R for different values of Gr . It is seen that the
velocity of air and water increases with increasingGr . At a
smaller Gr the velocity distribution is monotonic, but at a
higher time it passes through a maximum near the moving
plate when the buoyancy effect partly suppresses the inertial
effects of the plate velocity. Moreover, the velocity of air is
greater than the velocity of water. Physically this is possible
because fluids with high Prandtl number have greater
viscosity, which makes the fluid thick and hence move slowly.
Figure (2) presents the velocity profiles of air for different
values of K . It is observed that the velocity increases with
increasing K . Physically this is possible because as the
Grashof number or time increases, the contribution from the
buoyancy force near the moving hot plate become more
significant and hence a small rise in the fluid velocity near the
plate is observed. Figure (3) illustrate the influences of M
International Journal of Engineering Research (ISSN : 2319-6890)
Volume No.2, Issue No.2, pp : 113-118 01 April 2013
IJER@2013 Page 116
(Magnetic parameter) on velocity profiles respectively. It is
found that the velocity decreases with increase of magnetic
parameter M. It is also found that the velocity decreases away
from the plate and becomes minimum and finally takes
asymptotic value. Finally here we also see that point of
separation takes place for different values of magnetic
parameter. Figure (4) present typical profiles for the velocity,
for various values of a heat source
( ) | respectively. As
shown, the velocity decreasing with increasing| . In the
event that the strength of the presence of a heat source
( ) |
effect causes a reduction in the thermal state of the fluid, thus
producing lower thermal boundary layers. The effect of the
Prandtl number on the velocity shown in figure (5). As the
Prandtl number increases, the velocity decreases. Further,
figure (6) it is observed that the fluid velocity decreases with
increasing value of R. This result may be explained by the fact
that an increase in the radiation parameter
4
h
kk
R
T o
-
| |
=
|
'
\ .
for
fixed k and
h
T
'
means an increase in the Rosseland mean
absorption coefficient k
-
. When radiation is present, the
momentum boundary layer was found to be thicken, which is
in agreement with the observation made earlier with regard to
the temperature variations of air. Figure (7) present typical
profiles for the temperature for various values of a heat source
( ) | respectively. As shown, the temperature decreasing with
increasing. Figure (8) shows that the temperature profiles for
different values of Prandtl number
( ) Pr . It is observed that an
increase in the Prandtl number results a decrease of the
thermal boundary layer thickness and in general lower average
temperature within the boundary layer. The reason is that
increasing values of Prandtl number equivalent to increase the
thermal conductivities and therefore heat is able to diffuse
away from the heated plate more rapidly. Hence in the case of
increasing Prandtl numbers, the boundary layer is thinner and
the heat transfer is reduced. Boundary layer suction is the
technique in which air pumps is used to extract the boundary
layer at the wing. Further, increment of suction parameter
decreases the fluid temperature. Figure (9) shows the
temperature profiles for different values of the Radiation
parameter R, it is noticed that an increase in the radiation
parameter results decrease in the temperature with in boundary
layer, as well as decreased the thickness of the temperature
boundary layers. Figure (10) presents the skin-friction
variation with M in the pure convection case for different
values of Gr at the moving plate. It is observed that the skin-
friction increases with increasing Gr .
5 CONCLUSIONS
The temperature of the fluid increases with increasing time
whereas it decreases due to an increase in the value of
radiation parameter. In the case of pure convection (i.e. in the
absence of radiation), the velocity of the fluid increases with
increasing Grashof number, but falls owing to an increase in
the Prandtl number. The velocity of the fluid increases with
increasing Grashof number and time but it decreases owing to
an increase in the value of the radiation parameter. The skin-
friction at the moving plate increases with increasing values of
Grashof number and time for air flows.
We may conclude therefore, that the interaction between
the radiation, buoyancy forces and the applied shear induced
by a uniform vertical motion of the hot wall can affect the
configuration of the flow field significantly.
REFERENCES
[1] Alam M S, Rahman M M and Sattar M A: MHD Free
convection heat and mass transfer flow past an inclined
surface with heat generation, Thamasat. Int. J. Sci. Tech.
11 (4) (2006), pp. 1 8.
[2] Attia H A: The Effect of Variable Properties on the
Unsteady Couette Flow with Heat Transfer Considering
the Hall Effect, Communications in Nonlinear Science
and Numerical Simulation,13, (2008), pp. 1596- 1604.
[3] Chamkha A J: Unsteady MHD convective heat and mass
transfer past a semi- infinite vertical permeable moving
plate with heat absorption, Int. J. Eng. Sci. 24 (2004), pp.
217 230.
[4] Choi C K, Chung T J and Kim M C: Buoyancy Effects in
Plane Couette Flow Heated Uniformly From Below with
Constant Heat Flux, International Journal of Heat and
Mass Transfer, Vol. 47, (2004), pp. 2629-2636.
[5] Chung T J, Computational Fluid Dynamics, Cambridge
University Press, 2002.
[6] Cookey C I, Ogulu A and Omubo-Pepple V M: Influence
of viscous dissipation and radiation on unsteady MHD
free-convection flow past an infinite heated vertical plate
in a porous medium with time-dependent suction, Int. J.
Heat Mass Transfer, 46, (2003), pp. 2305 2311.
[7] El-Hakiem M A: MHD oscillatory flow on free-convection
radiation through a porous medium with constant
velocity, J. Magn. And Magnetic Mater. 220 (2, 3),
(2000), pp. 271 276.
International Journal of Engineering Research (ISSN : 2319-6890)
Volume No.2, Issue No.2, pp : 113-118 01 April 2013
IJER@2013 Page 117
[8] Grief R, Habib I S and Lin J C: Laminar convection of
radiating gas in a vertical channel, J. Fluid Mech. 46,
(1971), pp. 513 520.
[9] Hossain M A, Molla M M and Yaa L S: Natural
convection flow along a vertical wavy surface
temperature in the presence of heat generation/absorption,
Int. J. Thermal Science, 43, (2004), pp. 157 163.
[10] Jain N C and Gupta P: Three Dimensional Free
Convection Couette Flow with Transpiration Cooling,
Journal of Zhejiang University SCIENCE A, 3, (2006),
pp. 340-346.
[11] Jha B K: Natural Convection in Unsteady MHD Couette
Flow, Heat and Mass Transfer, 37, (2001), pp. 329-331.
[12]Kumar J, Lakshmana Rao C and Massoudi M: Couette
Flow of Granular Materials, International Journal of Non-
Linear Mechanics, 38, (2003), pp. 11-20.
[13] Makinde O D: Free-convection flow with thermal
radiation and mass transfer past a moving vertical porous
plate, Int. Comm., Heat Mass Transfer, 32, (2005), pp.
1411 1419.
[14] Narahari M: Effects of Thermal Radiation and Free
Convection Currents on the Unsteady Couette Flow
Between Two Vertical Parallel Plates with Constant Heat
Flux at one Boundary, Wseas transactions on heat and
mass transfer, 1 (5), (2010)
[15] Siegel R and Howell J R, Thermal Radiation Heat
Transfer, 4th Edition, Taylor & Francis, (2002).
[16] Soundalgekar V M: Hall Effects in MHD Couette Flow
with Heat Transfer, IEEE Transactions on Plasma
Science, Vol. PS-14, (5), (1986), pp. 579-583.
0 0.2 0.4 0.6 0.8 1
-0.2
0
0.2
0.4
0.6
0.8
1
1.2
1.4
y
Figure 1: Velocity profiles for different values of Gr (Pure convection case)
u
Gr=5.0,10.0,15.0,20.0
R=100,|=1.0,K=1.0,M=1.0
Pr=7.0
Pr=0.71
0 0.2 0.4 0.6 0.8 1
-0.2
0
0.2
0.4
0.6
0.8
1
1.2
y
Figure 2: Velocity profiles for different values of K
u
K=1.0,2.0,3.0,4.0
R=100.0,|=1.0,Gr=10.0,M=1.0,Pr=0.71
0 0.2 0.4 0.6 0.8 1
0
0.2
0.4
0.6
0.8
1
1.2
1.4
y
Figure 3: Velocity profiles for different values of M
u
R=5.0,|=2.0,K=1.0,Gr=15.0,Pr=0.71
M=1.0,2.0,3.0,4.0
0 0.2 0.4 0.6 0.8 1
0
0.2
0.4
0.6
0.8
1
1.2
1.4
y
Figure 4: Velocity profiles for different values of |
u
|=1.0,2.0,3.0,4.0
R=2.0,M=1.0,K=1.0,Gr=15.0,Pr=0.71
0 0.2 0.4 0.6 0.8 1
-0.2
0
0.2
0.4
0.6
0.8
1
1.2
y
Figure 5: Velocity profiles for different values of Pr
u
Pr=1.0,2.0,3.0,4.0
R=100.0,|=1.0,K=1.0,M=1.0,Gr=10.0
0 0.2 0.4 0.6 0.8 1
0
0.2
0.4
0.6
0.8
1
1.2
1.4
y
Figure 6: Velocity profiles for different values of R
u
R=1.0,2.0,3.0,4.0
|=2.0,M=1.0,K=1.0,Gr=15.0,Pr=0.71
International Journal of Engineering Research (ISSN : 2319-6890)
Volume No.2, Issue No.2, pp : 113-118 01 April 2013
IJER@2013 Page 118
0 0.2 0.4 0.6 0.8 1
0
0.2
0.4
0.6
0.8
1
y
Figure 7.: Temperature profiles for different values of Pr
u
Pr=0.71,R=1.0
|=1.0,2.0,3.0,4.0
0 0.2 0.4 0.6 0.8 1
0
0.05
0.1
0.15
0.2
0.25
0.3
0.35
0.4
0.45
y
Figure 8: Temperature profiles for different values of Pr
u
|=20.0,R=1.0
Pr=0.7,0.8,0.9,1.0
0 0.2 0.4 0.6 0.8 1
0
0.05
0.1
0.15
0.2
0.25
0.3
0.35
y
Figure 9: Temeprature profiles for different values of R
u
Pr=0.71,|=30.0
R=1.0,2.0,3.0,4.0
1 1.5 2 2.5 3 3.5 4
-2.2
-2
-1.8
-1.6
-1.4
-1.2
-1
-0.8
-0.6
-0.4
M
Figure 10: Skin friction for different values of Gr versus M
t
0
R=1.0,Pr=0.71,|=1.0,K=1.0
Gr=1.0,2.0,3.0,4.0