 Open Access
 Total Downloads : 950
 Authors : P. R. Sharma, I. K. Dadheech
 Paper ID : IJERTV1IS3150
 Volume & Issue : Volume 01, Issue 03 (May 2012)
 Published (First Online): 30052012
 ISSN (Online) : 22780181
 Publisher Name : IJERT
 License: This work is licensed under a Creative Commons Attribution 4.0 International License
Effect of Joule Heating on Steady MHD Flow of Low Prandtl Fluid on a Porous Stretching Sheet
P. R. Sharma
Professor, Department of Mathematics, University of Rajasthan, Jaipur302004, India.
Abstract
Aim of this paper is to investigate the effect of Joule heating on steady MHD flow of low Prandtl viscous and conducting fluid on a porous stretching sheet. In this analysis, the governing equations (viz. continuity, momentum and energy equations) are transformed into a system of ordinary differential equations. The momentum equation is solved analytically and energy equation is solved numerically using fourth order RungeKutta method along with shooting method. The fluid velocity and temperature characteristics are discussed and presented through graphs. The effects of skinfriction coefficient (viscous drag) and Nusselt number (rate of heat transfer) are derived, discussed and presented through tables for various physical parameters.

Introduction
The fluid dynamics over a stretching sheet is important in many practical applications such as extrusion of plastic sheet, paper production, glass blowing, metal spinning and drawing plastic films. Crane [1] was the first to consider the boundary layer flow caused by a stretching sheet, which moves with a velocity varying linearly with the distance from a fixed point. The heat transfer aspect of this problem was investigated by Carragher and Crane [2], under the conditions when the ambient fluid is proportional to a power of the distance from a fixed point. Fluid flows and heat transfer characteristics on stretching sheet with variable temperature condition have been investigated by Grubka and Bobba [3]. Elbashbeshy [4] discussed the heat transfer over a stretching sheet with variable heat flux. Elbashbeshy and Bajid [5] analyzed the stretching problem with internal heat generation and suction or injection in porous medium.

K. Dadheech
Assistant Professor, Birla Institute of Technology, Mesra (Ranchi)
Ext. Center Jaipur, India.
Effect of variable thermal conductivity and heat source/sink on MHD flow near a stagnation point on a linear stretching sheet was investigated by Sharma and Singh [6]. In this similar context Tamanna Sultana et al. [7] studied heat transfer in a porous medium over a stretching surface with internal heat generation and suction or injection in the presence of radiation.
The low Prandtl fluid flows have become important due to industrial applications, for instance it is used to deal with the problem of cooling of nuclear reactor by fluid having very low Prandtl number [Michiyoshi et al. (8), Fumizawa (9)]. Liquid metals have small Prandtl number of order 0.001~0.1 [e.g. Pr = 0.01 is for Bismuth, Pr =
0.023 is for Mercury etc.] and generally used as coolant because of very high thermal conductivity. They have ability to transport heat even if small temperature difference exists between the surface and the fluid. Due to this reason, liquid metals are used as coolant in nuclear reactors for the disposal of waste heat [Sharma and Singh (10)]. Kay (11), Arunachalam and Rajappa (12), Chaim (13) etc. presented low Prandtl fluid flow taking various geometries.
An analysis of thermal boundary layer in an electrically conducting fluid over a linearly stretching sheet in the presence of the constant transverse magnetic field with suction or blowing at the sheet was carried out by Chaim [14]. In this paper the viscous and Joule dissipation and internal heat generation was taken into account in the energy equation. The problem of viscous dissipation, Joule heating and heat source/sink on nonDarcy natural convection flow over an isoflux permeable sphere in a porous medium is numerically analyzed by Yih [15].
In this present study, the effect of Joule heating is investigated on steady MHD flow on a porous stretching sheet with a fluid of low Prandtl number.


Formulation of the Problem
Consider steady twodimensional laminar viscous incompressible fluid over a porous stretching sheet in the presence of transverse magnetic field. The
stretching sheet has a uniform temperature T and
The momentum and energy equations can be transformed into ordinary differential equations by taking the similarity variable and dimensionless temperature parameter respectively as
w
linear velocityuw x . Stretching sheet is placed in
y
uw
and
T T
(6)
the plane
y 0
and xaxis is taken along the
Tw T
sheet. Two equal and opposite forces are applied along x axis hence the sheet is stretched linearly keeping the origin fixed as demonstrated in figure 1.
The governing equations of continuity, momentum and energy of laminar boundary layer are given as
Using eq. (6), the equations (2) and (3) along with the boundary conditions are reduced to
f /// ff // ( f / )2 M 2 f / 0 , (7)
// Pr f / M 2 Ec Pr( f / )2 0 , (8)
u u 0
, (1)
where f is nondimensional stream function. Prime
x y
u u v
x
u
y
2u
y 2
o B 2
0 u
, (2)
denotes the differentiation with respect to . Corresponding boundary conditions are reduced to
T T 2T 2 2
f f w , f 1, 1
/
at 0 ,
C p u x v y y 2

B0 u
(3)
f / 0 , 0
at , (9)
where u, v are the velocity components along x and y axes, respectively, is the coefficient of
where
f v
u
is injection and
viscosity, is the density, is electrical conductivity, B0 is magnetic field intensity, Cp is
specific heat at constant pressure and is the thermal conductivity.
The corresponding boundary conditions are given
w w
0
suction velocity at the plate for
f w 0 , respectively.
f w 0
and
by
u uw x , v vw ,T Tw
, at y = 0
The governing boundary layer equation (7) along with the boundary conditions (9) admits a solution [Ahmad and Mubeen (16)] of the form given by
u 0 , T T ,
at y
(4)
1 1
Positive and negative value of vw indicates the
f ( )
f w r
e r , (10)
r
blowing and suction, respectively and obviously
vw 0
is corresponding to nonporous sheet.
where
Since, the low Prandtl fluid is taken into account the order of viscous dissipation term is considered less than the order of Joule heating i.e.
2
r f w
f 2 4(M 2 1)
w
2
. (11)
0
u
C y
B 2u 2 .
Further the energy equation (8) with the boundary conditions (9) is solved using RungeKutta fourth
p


Method of Solution
The continuity equation (1) is satisfied by introducing the stream function
order method along with shooting technique [Conte and Boor (17) and Sharma and Singh (6)].

Skinfriction Coefficient
The shear stress at the sheet is given by
x, y
uw xf ()
such that
u
w y
. (12)
y0
u ,
y
v
x
. (5)
The skinfriction coefficient at the sheet is defined as
C w 2x
u
f 1 2
2 w
3 1/p>
2 Re 2
f // (0)
(13)
velocity of fluid and therefore temperature of the fluid near the sheet is higher. It is observed from figure 7 that with the increase in Prandtl number temperature profiles decreases.
The skinfriction coefficient and Nusselt number
C f
f // (0) .
are presented by equations (13) and (15), which are directly proportional to f // (0) and / 0,

Nusselt Number
The local heat flux is given as
w y
q (x) T . (14)
y 0
respectively. The effect of physical parameters on these two are presented in the table 2 and table 3.
It is observed from table2, for the both cases (suction and injection) the value of f // (0) decreases as the value of Hartmann number increases. It is noted from table 3 the numerical
values of / 0 decrease as fw increases. Further
Now the Nusselt Number at the sheet is defined as
in case of suction, the value of / 0 decreases
x
Nu
qw (x)
(Re
) 12 / (0)
(15)
as the value of Prandtl number Pr increases and
(Tw T )
reverse result is found in case of injection. Also when Hartmann number increases the values of
x
Nu x (0) .
/
/ 0 increase.

Particular Case
In the absence magnetic field (i.e. M = 0), the results of present paper are reduced to those obtained by Chaim [13]. It is seen from Table 1 that the numerical results / (0) of present paper are in good agreement.

Results and Discussion
The energy equation (8) with the boundary conditions (9) is solved using RungeKutta fourth order method along with shooting technique for different values of the parameters M, Pr and fw taking step size 0.001.
The numerical calculations are presented in the form of tables and graphs for various values of physical parameters.
Figure 2 shows the velocity profiles for different values of suction and injection parameter fw with M
= 0.5. It can be seen easily that the velocity profiles decrease with the increase of fw. The Hartmann number (M) shows an importance of magnetic field sets in Lorentz force, which results in retarding force on the velocity field and therefore as M
Figure 1. Schematic diagram of the problem.
1.2
0.8
fw = 0.5
f fw = 0.1
fw = 0.0
increases the velocity profiles decrease. As observed from figures 3 and 4 it is also noted that the velocity profiles decrease, while taking increment in M for fw = 0.5 (figure 3, in case of suction) and for fw =0.5 (figure 4, in case of
0.4
0
fw = 0.1
fw = 0.5
injection). Figures 5 and 6 show the effect of the Hartmann number on temperature profiles at Pr=0.01 for fw = 0.5 (figure 5, in case of suction) and for fw =0.5 (figure 6, in case of injection), respectively and it is observed that temperature profiles increase with the increase in Hartmann number. Because the magnetic field retards the
0 0.5 1
Figure 2. Velocity distribution when M = 0.5
1.2 1.2
0.8
f
0.4
0
M = 0.5
M = 1.0
M = 1.5
M = 2.5
0 0.5 1
0.8
0.4
0
M = 0.5
M = 1.0
M = 1.5
0 0.5 1
Figure 3. Velocity distribution when
f w = 0.5.
Figure 6. Temperature distribution when
f w = 0.5.
1.2 1.2
0.8
f M = 0.5
M = 1.0
0.8
Pr = 0.01
Pr = 0.05
0.4
0
M = 1.5
M = 2.5
0 0.5 1
0.4
0
Pr = 0.07
0 0.5 1
Figure 4. Velocity distribution when
f w = 0.5.
Figure 7. Temperature distribution when M = 0.5.
Table 1. Numerical values of / (0) for various values of Pr are compared with the results obtained by Chaim (1998)
– / (0)
Pr
f w
Chaim [13]
Present work
0.023
0.0
0.0224886
0.02302313
0.023
1.5
0.0457422
0.04576912
0.1
0.0
0.0912924
0.0913340
0.1
1.5
0.1955034
0.19554759
1.0
0.0
0.5819767
0.5836673
1.0
1.5
1.760407
1.76268901
Figure 5. Temperature distribution when f w = 0.5.
Table 2. Numerical values of
f // (0) for various
2. Carragher, P. and Crane, L. J., Heat transfer on a continuous stretching sheet. Z. Angew. Math. Mech., Vol. 62, 1982, pp. 564573.
values of M and fw.
fw =0.5
M
f // (0)
0.5
1.395644
1.0
1.686141
1.5
2.070027
2.5
2.954163
fw = 1.0
0.5
1.724745
1.0
2.000
1.5
2.320829
2.5
3.238613
fw = 1.5
0.5
2.096291
1.0
2.350781
1.5
2.702562
2.5
3.545085
fw = 0.5
0.5
0.895644
1.0
1.186141
1.5
1.570027
2.5
2.454163
fw = 1
0.5
0.724745
1.0
1.000
1.5
1.370829
2.5
2.238613
fw = 1.5
0.5
0.596291
1.0
0.850781
1.5
1.202562
2.5
2.045085
Table 3. Numerical values of / (0) for various values of physical parameters
fw =0.5
Pr
M
/ 0
.01
0.5
0.4050120
.01
1.0
0.4038703
.01
1.5
0.4024835
fw = 1.0
.01
0.5
0.4073560
.01
1.0
0.4063797
.01
1.5
0.4051136
.07
0.5
0.4481621
fw = 0.5
.01
0.5
0.4003308
.01
1.0
0.3988151
.01
1.5
0.3971187
fw = 1
.01
0.5
0.3979683
.01
1.0
0.3962657
.01
1.5
0.3944089
.07
0.5
0.3812571

References
1. Crane, L. J., Flow past a stretching plate. ZAMP, Vol. 21, 1970, pp. 645647.
 Grubka, L. J. and Bobba K. M., Heat transfer characteristics of a continuously stretching surface with variable temperature. Transactions of ASME Journal of Heat and Mass Transfer, Vol. 107, 1985, pp. 248250.

Elbashbeshy, E. M. A., Heat transfer over a stretching surface with variable heat flux. J. Physics D: Appl. Physics, Vol. 31, 1998, pp. 19511955.

Elbashbeshy, E. M. A. and Bajid, M. A. A., Heat transfer in a porous medium over a stretching surface with internal heat generation and suction or injection. Appl. Math. Comput., Vol. 158, 2004, pp. 699703.

Sharma P. R. and Singh G., Effect of variable thermal conductivity and heat source/sink on MHD flow near a stagnation point on a linearly stretching sheet. JAFM, Vol. 2, 2009, pp. 13 21.

Sultana T., Saha S., Rahman M. M. and Saha Gautam, Heat transfer in porous medium over a stretching surface with internal heat generation and suction or injection in the presence of radiation. Journal of Mechanical Engineering, Vol. 40, 2009, pp. 2228.

Michiyoshi, I., takahashi, I. And Serizawa, A., Natural convection heat transfer from a horizontal cylinder to mercury under magnetic field. Int. J. Heat Mass Transfer, Vol. 19, 1976, pp. 10211029.

Fumizawa, M., Natural convection experiment with Nak under transverse magnetic field. Journal of Nuclear Science and Technology, Vol. 17, 1980, pp. 98105.

Sharma P. R. and Singh G., Steady MHD natural convection flow with variable electrical conductivity and heat generation along an isothermal vertical plate. Tamkang J. Science Engineering, Vol. 13 , 2010, pp. 235242.

Kay, W. M., Convective Heat and Mass Transfer. McGraw Hill Book Co., New York, 1966.

Arunachalam, M. and Rajappa, N. R., Forced convection in liquid metals with variable thermal conductivity and capacity. Acta Mechanica, Vol. 31, 1978, pp. 2531.

Chaim, T. C, Heat transfer in a fluid with variable thermal conductivity over stretching sheet. Acta Mechanica, Vol. 129, 1998, pp. 6372.

Chaim, T. C., MHD heat transfer over a non isothermal stretching sheet. Acta Mechanica, Vol. 122, 1977, pp. 169179.

Yih, K. A., Viscous and Joule heating effects on nonDarcy MHD natural convection flow over a permeable sphere in porous media with internal heat generation. Int. Commun. Heat Mass, Vol. 27, 2000, pp. 591600.

Ahmad N. and Mubeen A., Boundary layer flow and heat transfer for the stretching plate with suction. Int. Comm. Heat Mass Transfer, Vol. 22, 1995, pp. 895906.

Conte, S. D. and Boor, C., Elementary Numerical Analysis. McGrawHill Book Co., New York, 1981.