Effects of conduction on magneto hydrodynamics mixed convection flow in triangular enclosures

Mathematical Theory and Modeling www.iiste.org
ISSN 2224-5804 (Paper) ISSN 2225-0522 (Online)
Vol.3, No.5, 2013
49
Effects of Conduction on Magneto-Hydrodynamics Mixed
Convection Flow in Triangular Enclosures
1
Muhammad Sajjad Hossain*, 2
M. A. Hoque and 3
Md. Kamruzzaman
1
Assistant Professor of Mathematics, 2
Lecturer in Mathematics, 3
Lecturer in Physics
Department of Arts and Sciences , Ahsanullah University of Science & Technology (AUST)
Dhaka-1208, Bangladesh.
*Corresponding author’s email: msh80_edu@yahoo.com
ABSTRACT
We study the significance of conduction with magneto-hydrodynamics (MHD) mixed convection flow in
triangular enclosures. The enclosures consist of one temperature of the left moving wall, which has constant flow
speed. A uniform magnetic field is applied in the horizontal direction normal to the moving wall. The results
indicate that both the flow and the thermal fields strongly depend on the Hartmann number Ha, Prandlt number
Pr for fixed Reynolds number Re, at the convective regimes. It is also observed that, the parameters Prandtl
number Pr influence on the flow fields and have significant effect on the thermal fields at the convective regimes.
According to the results it is found that both the flow and thermal fields inside the enclosures strongly depend on
the relevant dimensionless groups. For the purpose of comparison of the effect of MHD on heat transfer, the
results in terms of average Nusselt number and average temperature are shown in tabular form.
Keywords: Streamfunction, Isotherms, Triangular Enclosures, MHD
1. INTRODUCTION
Mixed convection is a combination of forced and free convections which is the general case of convection when
a flow is determined simultaneously by both an outer forcing system (i.e., outer energy supply to the fluid-
streamlined body system) and inner volumetric (mass) forces, viz., by the nonuniform density distribution of a
fluid medium in a gravity field. The most vivid manifestation of mixed convection is the motion of the
temperature stratified mass of air and water areas of the Earth that the traditionally studied in geophysics.
However, mixed convection is found in the systems of much smaller scales, i.e., in many engineering devices.
On heating or cooling of channel walls, and at the small velocities of a fluid flow that are characteristic of a
laminar flow, mixed convection is almost always realized. Many mathematicians, versed engineers and
researchers studied the problems of MHD mixed convection for different types of cavity. Amongst them, Oreper
and Szekely (1983) studied the effect of an externally imposed magnetic field on buoyancy driven flow in a
rectangular cavity. They found that the presence of a magnetic field can suppress natural convection currents and
that the strength of the magnetic field is one of the important factors in determining the quality of the crystal.
Ozoe and Maruo (1987) investigated magnetic and gravitational natural convection of melted silicon two-
dimensional numerical computations for the rate of heat transfer. Rudraiah et al. (1995a) investigated the effect
of surface tension on buoyancy driven flow of an electrically conducting fluid in a rectangular cavity in the
presence of a vertical transverse magnetic field to see how this force damps hydrodynamic movements. At the
same time, Rudraiah et al. (1995b) also studied the effect of a magnetic field on free convection in a rectangular
enclosure. The problem of unsteady laminar combined forced and free convection flow and heat transfer of an
electrically conducting and heat generating or absorbing fluid in a vertical liddriven cavity in the presence of a
magnetic field was formulated by Chamkha (2002).Recently, Rahman et al. (2009a) studied the effect of a heat
conducting horizontal circular cylinder on MHD mixed convection in a lid-driven cavity along with Joule
Heating. The author considered the cavity consists of adiabatic horizontal walls and differentially heated vertical
walls, but it also contains a heat conducting horizontal circular cylinder located somewhere inside the cavity. The
results indicated that both the flow and thermal fields strongly depend on the size and locations of the inner
cylinder, but the thermal conductivity of the cylinder has significant effect only on the thermal field. Rahman et
al.(2009b) investigated the effect on magnetohydrodynamic (MHD) mixed convection around a heat conducting
Horizontal circular cylinder placed at the center of a rectangular cavity along with joule heating. They observed
that the streamlines, isotherms, average Nusselt number, average fluid temperature and dimensionless
temperature at the cylinder center strongly depend on the Richardson number, Hartmann number and the cavity
aspect ratio. Very recently, Rahman et al. (2010) made numerical investigation on the effect of magneto
hydrodynamic (MHD) mixed convection flow in a vertical lid driven square enclosure including a heat
conducting horizontal circular cylinder with Joule heating. The author found that the Hartmann number,
Reynolds number and Richardson number have strong influence on the streamlines and isotherms.
To the best of our knowledge in the light of the above literatures, it has been pointed out that no work has been
paid on MHD mixed convection inside a triangular enclosure yet.The aim proposed research is to investigate the

Vol.3, No.5, 2013
50
effects of Magneto hydrodynamic (MHD) mixed convection flow in a triangular enclosures. Results will be
presented for different non-dimensional governing and physical parameters in terms of streamlines, isotherms,
heat transfer rate as well as the average temperature of the fluid in the enclosures.
2. PHYSICAL MODEL
The physical model considered here is shown in Fig. 1, along with the important geometric parameters. The heat
transfer and the fluid flow in a two-dimensional triangular enclosure with MHD flow whose left wall and bottom
wall are subjected to cold Tc and hot Th temperatures respectively while the inclined walls are kept adiabatic.
The fluid was assumed with Prandtl number (Pr = 0.71, 2.0, 3.0, 6.0), Reynolds number (Re = 50), Hartmann
number (Ha = 5-50) and Newtonian, and the fluid flow is considered to be laminar. The properties of the fluid
were assumed to be constant.
Figure 1: Schematic diagram of the physical system
3. MATHEMATICAL FORMULATION
The fundamental laws used to solve the fluid flow and heat transfer problems are the conservation of mass
(continuity equations), conservation of momentums (momentum equations), and conservation of energy (energy
equations), which constitute a set of coupled, nonlinear, partial differential equations. For laminar incompressible
thermal flow, the buoyancy force is included here as a body force in the v-momentum equation. The governing
equations for the two-dimensional steady flow after invoking the Boussinesq approximation and neglecting
radiation and viscous dissipation can be expressed as:
0=
∂
∂
+
∂
∂
Y
V
X
U
(1)
2 2
2 2
1U U P U U
U V
X Y X Re X Y
 ∂ ∂ ∂ ∂ ∂
+ = − + + 
∂ ∂ ∂ ∂ ∂ 
(2)
V
Ha
Ra
Y
V
X
V
Y
P
Y
U
V
X
U
U
Re
Pr)(
Re
1 2
2
2
2
2
−+
∂
∂
+
∂
∂
+
∂
∂
−=
∂
∂
+
∂
∂
θ (3)
2 2
2 2
1
U V
X Y Re Pr X Y
θ θ θ θ ∂ ∂ ∂ ∂
+ = + 
∂ ∂ ∂ ∂ 
(4)
The dimensionless parameters appearing in the equations (4.5) through (4.8) are Reynolds number, ,Re
ν
Lui=
Grashof number,
2
3
ν
β TLg
Gr
∆
= , Prandtl number, ,Pr
α
ν
= Hartmann number,
µ
σ 22
02 LB
Ha = ,
Rayleigh number, Pr×= GrRa , temperature difference, ch TTT −=∆ and thermal diffusivity,
B0
g
Th
Adiabatic wall
X
Y
v
Tc

Vol.3, No.5, 2013
51
pCρ
κ
α = , where non-dimensional variables are as follows:
)(
)(
,,,,, 2
000 ch
c
TT
TT
v
p
P
v
v
V
v
u
U
L
y
Y
L
x
X
−
−
====== θ
ρ
Where X and Y are the coordinates varying along horizontal and vertical directions, respectively, U and V are the
velocity components in the X and Y directions, respectively, θ is the dimensionless temperature and P is the
dimensionless pressure.
The appropriate dimensionless boundary conditions for the above equations () can be written as:
At the left vertical wall: U=0, V=1, 0=θ
At the bottom wall: U=0, V=0, 1=θ
At the inclined wall: U=0, V=0, 0=
∂
∂
N
θ
Where N is the non-dimensional distances either along X or Y direction acting normal to the surface and K is the
dimensionless thermal conductivity.
According to Singh and Sharif (2003), the average Nusselt number at the heated wall of the enclosures based on
the non-dimensional variables may be expressed as dX
Y
Nu
y 0
1
0 =
∫ 





∂
∂
=
θ
and the bulk average temperature
defined as /av dV Vθ θ= ∫ , where is the enclosures volume.
4. FINITE ELEMENT TECHNIQUE
The Galerkin finite element method (1973, 1999) is used to solve the non-dimensionalgoverning equations along
with boundary conditions for the considered problem. The equation of continuity has been used as a constraint
due to mass conservation and thisrestriction may be used to find the pressure distribution. The basic unknowns
for the above differential equations are the velocity components U, V, the temperature θ and the pressure P. The
velocity components, the temperature distributions and linear interpolation for the pressure distribution
according to their highest derivative orders in the differential eqs. (1) - (4) are as follows:
( ),U X Y N Uα α=
,
( ),V X Y N Vα α=
,
( , )X Y Nα αθ θ= ,
( ),P X Y H Pλ λ=
where α = 1, 2, … …, 6; λ= 1, 2, 3; Nα are the element interpolation functions for the velocity components and
the temperature, and Hλ are the element interpolation functions for the pressure.
To derive the finite element equations, the method of weighted residuals finite-element formulation is applied to
the equations (1) – (4) as
0=





∂
∂
+
∂
∂
∫ dA
Y
V
X
U
NA α (5)
2 2
2 2
1
A A A
U U P U U
N U V dA H dA N dA
X Y X Re X Y
α λ α
 ∂ ∂ ∂ ∂ ∂   
+ = − + +    
∂ ∂ ∂ ∂ ∂     
∫ ∫ ∫ (6)
V

Vol.3, No.5, 2013
52
VdAN
Ha
dANRadA
Y
V
X
V
N
dA
Y
P
HdA
Y
V
V
X
V
UN
AAA
AA
∫∫∫
∫∫
−+
∂
∂
+
∂
∂
+
∂
∂
−=
∂
∂
+
∂
∂
ααα
λα
θ
Re
Pr)(
Re
1
)()(
2
2
2
2
2
(7)
2 2
2 2
1
A A
N U V dA N dA
X Y Re Pr X Y
α α
θ θ θ θ ∂ ∂ ∂ ∂ 
+ = +  
∂ ∂ ∂ ∂   
∫ ∫ (8)
So, the coefficients in element matrices are in the form of the integrals over the element area as,
,x xA
K N N dAα βαβ
= ∫
,
,y yA
K N N dAα βαβ
= ∫
,
,x xA
K N N N dAα β γαβγ
= ∫
,y yA
K N N N dAα β γαβγ
= ∫
,
A
K N N dAαβ α β= ∫ ,
, ,xx x xA
S N N dAα βαβ
= ∫
, ,yy y yA
S N N dAα βαβ
= ∫
,
,x xA
M H H dAα µαµ
= ∫
,
,y yA
M H H dAα µαµ
= ∫
0
0
u xS
Q N S dSαα
= ∫ ,
0
0
v yS
Q N S dSαα
= ∫ ,
1 wwSw
Q N q dSθ αα
= ∫
2 wwSws
Q N q dSθ αα
= ∫
Using reduced integration technique of Reddy (1985) and Newton-Raphson iteration technique by first writing
the unbalanced values from the set of the finite element equations, the nonlinear algebraic Eqs. (5) - (8) so
obtained are modified by imposition of boundary conditions. These modified nonlinear equations are transferred
into linear algebraic equations by using Newton’s method. Finally, these linear equations are solved by using
Triangular Factorization method.
5. MESH GENERATION
The present numerical technique will discretize the computational domain into unstructured triangles by
Delaunay Triangular method. The Delaunay triangulation is a geometric structure that has enjoyed great
popularity in mesh generation since the mesh generation was in its infancy. In two dimensions, the Delaunay
triangulation of a vertex set maximizes the minimum angle among all possible triangulations of that vertex set.
Figure 2: Current mesh structure of elements for triangular enclosures.
6. GRID INDEPENDENCE TEST
A grid refinement study is performed for a square open enclosures to obtain grid refinement solution with Pr =
0.71, Re = 50, Ha = 20 and Ra = 104
. Figure 3 shows the convergence of the average Nusselt number, Nuav at the
heated surface with grid refinement. It is observed that grid independence is achieved with 13398 elements
where there is insignificant change in Nu with further increase of mesh elements.

Mathematical Theory and Modeling
Vol.3, No.5, 2013
Figure 3: Convergence of average Nusselt number with grid refinement for
104
.
7. CODE VALIDATION
The computed average Nusselt numbers for the constant surface temperature and different Rayleigh Numbers are
tabulated in Table 1 for numerical code validation against the numerical results of Basak et al. (2007), for pure
natural convection with Pr = 0.71 in a triangular enclosures. It is seen the present computation and those of
Basak et al.(2007), are good agreement.
Ra
103
104
105
Table 1: Comparison between Basak et.al.(2007) and present work
8. RESULTS AND DISCUSSION
The MHD mixed convection flow and temperature fields as well as heat transfer rates
inside the triangular enclosures has been numerically investigated. The MHD mixed convection phenomenon
inside the enclosure is strongly influenced by governing as well as physical parameters, namely and Hartman
number Ha, Reynolds number Re, Prandtl number
physical properties of fluid strongly influenced on the heat transfer rate, and the numerical computation were
performed for a range of Ha, Re and
mentioned parameters, computations are performed at
Moreover, the results of this study are presented in terms of streamlines and isotherms. Furthermore, the heat
transfer effectiveness of the enclosure is displayed in terms of average Nusselt number
average bulk temperature θav.
8.1 EFFECTS OF HARTMAN NUMBER
The special effects of Hartmann number
Figure 4 - 7, while Pr = 0.71 and Re
also decreases sharply with increasing of
propensity to slow down the movement of the buoyancy induced flow in the enclosures. The effects of Hartmann
number Ha on isotherms are shown in the Figure 4
that the higher values isotherms are more tighte
6.32
6.325
6.33
6.335
6.34
6.345
6.35
0
Nu
Mathematical Theory and Modeling
0522 (Online)
53
Figure 3: Convergence of average Nusselt number with grid refinement for Pr = 0.71, Re= 50,
ble 1 for numerical code validation against the numerical results of Basak et al. (2007), for pure
= 0.71 in a triangular enclosures. It is seen the present computation and those of
Basak et al.(2007), are good agreement.
Nuav
Present work Basak et al. (2007)
5.49 5.40
5.77 5.56
7.08 7.54
RESULTS AND DISCUSSION
The MHD mixed convection flow and temperature fields as well as heat transfer rates
, Prandtl number Pr, Rayleigh number Ra. As the nature of flow and thermo
and Pr (Ha =5-50, Re = 50, Pr = 0.71-6). In addition, for each value of
mentioned parameters, computations are performed at Ra =5 103
that focuses on pure mixed convection.
fer effectiveness of the enclosure is displayed in terms of average Nusselt number Nu
EFFECTS OF HARTMAN NUMBER
The special effects of Hartmann number Ha on the flow and thermal fields at Ra= 1e
Re = 50 are kept fixed. In this folder, the size of the counter clockwise vortex
also decreases sharply with increasing of Ha. This is because application of a slanting magnetic field has the
slow down the movement of the buoyancy induced flow in the enclosures. The effects of Hartmann
on isotherms are shown in the Figure 4 – 7 while Re = 50 and Pr = 0.71 are kept fixed. It is observed
that the higher values isotherms are more tightened at the vicinity of the heated wall of the enclosures for the
5000 10000 15000 20000
Elements Number
×
www.iiste.org
= 50, Ha=20 and Ra =
ble 1 for numerical code validation against the numerical results of Basak et al. (2007), for pure
= 0.71 in a triangular enclosures. It is seen the present computation and those of
The MHD mixed convection flow and temperature fields as well as heat transfer rates and bulk temperature
. As the nature of flow and thermo
6). In addition, for each value of
that focuses on pure mixed convection.
Nu and the dimensionless
= 1e3
to 1e4
is illustrate in
= 50 are kept fixed. In this folder, the size of the counter clockwise vortex
. This is because application of a slanting magnetic field has the
slow down the movement of the buoyancy induced flow in the enclosures. The effects of Hartmann
= 0.71 are kept fixed. It is observed
ned at the vicinity of the heated wall of the enclosures for the
20000

Vol.3, No.5, 2013
54
considered value of Ra. A similar trend is found for the higher values of Ha. From this figure it can easily be
seen that isotherms are almost parallel at the vicinity of the bottom horizontal hot wall but it changes to parabolic
shape with the increasing distance from inclined wall for the highest value of Ha (=20.0) at Pr = 0.71, indicating
that most of the heat transfer process is carried out by conduction. In addition, it is noticed that the isothermal
layer near the heated surface becomes slightly thin with the decreasing Ha. However in the remaining area near
the left wall of the enclosures, the temperature gradients are small in order to mechanically driven circulations.
The effects of Ha on average Nusselt number Nu at the heated surface and average bulk temperature θ av in the
enclosures is illustrated in Fig. 8 & 9 with Re =50 and Pr =0.71. From this figure, it is found that the average
Nusselt number Nu decreases with increasing Pr. It is to be highlighted that the highest heat transfer rate occurs
for the lowest values of Ha (= 5). The average bulk temperature of the fluid in the enclosures is high for higher
values of Ha (10 - 50). It is to be mentioned here that the highest average temperature is documented for the
higher value of Ha in the free convection dominated region.
Table 2: Average Nusselt numbers for different Hartmann number while Ha = 5, 10, 20 and 50, Re=50 and
Pr=0.71.
Nuav
Ha Ra=103
Ra=5×103
Ra=104
5 4.38836694 6.434903 5.764501
10 4.36317068 6.413537 5.805674
20 4.26278628 6.331658 5.913316
50 3.79764503 5.890227 5.878074
Table 3: Average bulk temperature for different Hartmann number while Ha= 5, 10, 20 and 50, Re=50 and
Pr=0.71.
θav
Ha Ra=103
Ra=5×103
Ra=104
5 0.593528 0.566768 0.57386
10 0.593913 0.567008 0.570757
20 0.598198 0.568812 0.561785
50 0.633216 0.585029 0.533999

Vol.3, No.5, 2013
55
Ra=103
Ra=5×103
Ra=104
Isotherm Streamline
Fig 4: Isotherms and streamlines patterns for Re = 50, Pr = 0.71 and Ha = 5
298.54
300.00
302.12
305.00
305.76
310.00
315.00
295.00
290.00
1.00 3.00
5.00
7.00
2.00
1.007.00
1.00
8.006.00
2.00
295.00
297.28
299.09
300.00
302.59
303.54
305.00
310.00
315.00
290.00
280.00
2.00
10.00
16.00
18.00
2.00
10.00
6.00
16.00
18.00
14.00
0.00
2.00
297.61
298.66
300.00
303.45
305.00
310.00
315.00
295.00
290.00
295.00
296.68
5.00
10.00
15.00
20.00
25.00
15.0015.00
25.00
25.00
5.00
10.00
25.00

Vol.3, No.5, 2013
56
Ra=103
Ra=5×103
Ra=104
Isotherm Streamline
298.57
300.00301.69
305.00
305.85
310.00
315.00
295.00
280.00
1.00
3.00 5.00
7.00
2.00
5.00
7.00
7.00
1.00
1.00
3.00
1.00
295.00
297.32
299.31
300.00
301.27303.49
305.00
310.00
315.00
280.00290.00
2.00
4.00
12.00
18.00
4.00
8.00
10.00
18.00
16.00
8.00
4.00
2.00
296.19
297.49
298.76
300.00
303.04
305.00
310.00
315.00
295.00
290.00
285.00290.00
5.00
5.00
15.00
20.00
25.00
15.00
20.00
5.00
5.00
20.00
10.00
20.00

Vol.3, No.5, 2013
57
Ra=103
Ra=5×103
Ra=104
Isotherm Streamline
318.10
310.00
306.14
305.00
302.28
300.00
297.69
295.00
280.00
1.00
2.00
4.00
6.00
1.006.00
4.00
6.00
3.00
1.00
295.00
298.01
300.00
300.73
302.81
303.80
310.00
315.00
290.00
280.00
302.81
2.00
8.00
12.00
16.00
14.00 2.00
16.00
2.00
2.00
295.00
296.66
298.23
299.43
300.00
302.01
305.00
310.00
315.00
290.00
290.00
290.00
5.00
10.00
15.00
20.00
10.00
15.00
20.0015.00
5.00
15.00

Vol.3, No.5, 2013
58
Ra=103
Ra=5×103
Ra=104
Isotherm Streamline
295.00
300.00
305.00
306.89
308.77
310.00
315.00
319.73
285.00
0.23
0.64
1.34
1.34
3.00
2.00
3.00
1.00
4.00
3.00
1.00
1.34
295.00
298.21
300.00
302.92
305.00
305.58
310.00
315.00
290.00
285.00280.00
0.48
2.00
4.00
6.00
10.00
10.00
0.90
6.00
8.00
10.00
8.00
12.006.00
10.004.00
4.00
0.90
0.90
295.00
297.11
298.94
300.00
302.21
302.51
305.00
310.00
315.00
290.00
285.00
290.00
2.24
2.96
5.00
10.00
15.00
10.00
5.00
2.96
10.00
15.00

Vol.3, No.5, 2013
59
Fig. 8: Effect of average Nusselt number and Rayleigh number while Pr = 0.71 and Re = 50.
Fig. 9: Effect of average bulk temperature and Rayleigh number while Pr = 0.71and Re = 50.
Ra
Nu
2000 4000 6000 8000 10000
3.5
4
4.5
5
5.5
6
6.5
Ha= 5
Ha= 10
Ha= 20
Ha= 50
Ra
2000 4000 6000 8000 10000
0.5
0.55
0.6
0.65
Ha=5
Ha=10
Ha=20
Ha=50
θav

Vol.3, No.5, 2013
60
8.2 EFFECTS OF PRANDLT NUMBER
The property of Prandtl number on flow fields for the considered Ra is depicted in Fig.10 - 13, while Re = 50
and Ha = 20 are kept fixed. Now, from the figure, it is seen that a small vortex is formed at the center of the
enclosures. It is evident that the size of the vortex is increased with the increasing value of Pr. The influence of
Prandtl number Pr on thermal fields for the selected value of Ra is displayed in Fig.12 for Re = 50 and Ha = 20.
As expected, the thermal boundary layer decreases in thickness as Pr increases. This is reflected by the denser
deserting of isotherms close to the heated wall and cooled wall as Pr increases. The spread of isotherms at low
values of Pr is due to a strong stream wise conduction that decreases the stream wise temperature gradient in the
field. At the increase of Pr the isotherms become linear in the enclosures and impenetrable to respective wall.
From the Fig. 13, it can easily be seen that the isotherms for different values of Pr at Ra = 104
are clustered near
the heated horizontal wall and cooled vertical heated surface of the enclosure, indicating steep temperature
gradient along the vertical direction in this region. Besides, the thermal layer near the heated wall becomes
thinner moderately with increasing Pr.
The effects of Prandtl number on average Nusselt number Nu at the heated surface and average bulk temperature
θ av in the enclosures are illustrated in Fig. 14 & 15 with Re = 50 and Ha = 20. From this figure, it is found that
the average Nusselt number Nu goes up penetratingly with increasing Pr. It is to be highlighted that the highest
heat transfer rate occurs for the highest values of Pr (= 6). The average bulk temperature of the fluid in the
enclosures declines very mildly for higher values of Pr (0.71- 6.0). But for the lower values of Pr, the average
bulk temperature is linear. It is to be mentioned here that the highest average temperature is documented for the
lower value of Pr (= 0.71) in the free convection dominated region.
Table 4: Average Nusselt numbers for different Prandtl number while Pr = 0.71, 2.0, 3.0 and 6.0, Re=50 and
Ha=20.
Nuav
Pr Ra=103
Ra=5×103
Ra=104
0.71 4.262786 6.331658 7.510701
2.0 4.998317 7.454621 8.857028
3.0 6.807629 10.118253 12.023071
6.0 8.111548 11.920293 13.940204
Table 5: Average bulk temperature for different Prandtl number while Pr = 0.71, 2.0, 3.0 and 6.0, Re=50 and
Ha=20.
θav
Pr Ra=103
Ra=5×103
Ra=104
0.71 0.598198 0.568812 0.558697
2.0 0.584269 0.550002 0.538424
3.0 0.547777 0.504842 0.490760
6.0 0.524512 0.475957 0.459601

Vol.3, No.5, 2013
61
Ra=10
3
Ra=5×10
3
Ra=10
4
Isotherm Streamline
298.09
300.00
302.28
305.00
306.51
310.00
315.00
295.00
290.00
285.00
285.00
1.0
2.0
4.0
6.0
2.06.0
4.0
6.0
1.0
1.0
295.00
297.77
299.17
300.00
303.18
305.00
310.00
315.00
290.00
285.00280.00
2.0
4.0
6.0
8.0
12.0
14.0
16.0
8.0
10.0
14.016.0
14.08.0
16.0
6.0
8.0
4.0
2.0
2.0
295.00
296.75
297.77
299.18
300.00
300.67
303.04
305.00
310.00
315.00290.00
285.00
2.49
5.00 15.00
25.00
15.00
25.00
20.00
2.49
20.00

Vol.3, No.5, 2013
62
Ra=10
3
Ra=5×10
3
Ra=10
4
Isotherm Streamline
295.00
297.51
298.86
300.00
302.24
305.00
310.00
315.00
290.00
285.00280.00285.00
1.0
2.0
4.0
6.0
2.0
2.0
7.0
7.0
5.0
1.0
7.0
292.70
295.00
297.26
298.10
300.00
305.00
310.00315.00
290.00
290.00
285.00
2.00
4.00
6.00
10.00
12.00
14.00 16.00
4.00
6.00
10.00
8.00
12.00
18.008.00
291.89
293.32
295.00
296.82
297.87
300.00
300.00
305.00
315.00
290.00
285.00
2.68
10.00
15.00
20.00
15.00
25.00 2.68
15.00
15.00
2.68

Vol.3, No.5, 2013
63
Ra=10
3
Ra=5×10
3
Ra=10
4
Isotherm Streamline
295.00
297.04
298.09
300.00
301.70
305.00
310.00
315.00
290.00
280.00
0.4
2.0
3.66.0
6.0
2.0
4.0
8.0
3.61.5
292.01
294.09
295.00
296.67
297.19
300.00
305.00
310.00
290.00
294.09
280.00
2.02.0
4.06.08.0
4.0
8.0
10.0
12.0
14.0
4.0
16.0
8.0
6.0
6.0
16.08.014.0
4.0
290.00
292.60
294.06
295.00
295.90
296.51
300.00
310.00
280.00
1.72
5.00
10.00
15.00
20.00
5.00
10.00
20.00
20.0010.0025.00
25.00

Vol.3, No.5, 2013
64
Ra=10
3
Ra=5×10
3
Ra=10
4
Isotherm Streamline
293.33
295.74
297.06
298.33
298.97
300.00
300.52
300.80
310.00
290.00
0.50
1.00
2.00 3.36
2.00
3.36
4.50
5.50
6.50
2.002.00
8.003.72
5.50
2.00
290.00
292.00
293.30
294.43
295.00
296.12
296.98
297.78
315.00285.00
292.00
2.00
4.00
4.00
8.00
12.00
14.00
18.00
4.00
4.00
10.00
14.00
290.00
291.81
293.32
294.12
294.74
295.00
296.05
315.00
285.00
5.00
10.00
10.00
15.00
20.00
5.00
20.00
25.0020.00

Vol.3, No.5, 2013
65
Fig.14: Effect of average Nusselt number and Rayleigh number while Re = 50 and Ha = 20
Fig. 15: Effect of average bulk temperature and Rayleigh number while Re = 50 and Ha=20.
Ra
Nu
2000 4000 6000 8000 10000
2
4
6
8
10
12
14
16
18
Pr= 0.71
Pr=2
Pr= 3
Pr= 6
Ra
2000 4000 6000 8000 10000
0.5
0.55
0.6
0.65
Pr= 0.71
Pr= 2
Pr= 3
Pr=6
θav

Vol.3, No.5, 2013
66
9. CONCLUSIONS
Two-dimensional laminar steady state MHD mixed convection flow in triangular enclosures with the bottom
wall at iso-heat flux has been studied numerically. Numerical results have been shown by considering Prandtl
numbers of 0.71, 2.00, 3.0 and 6.00, Hartman numbers of 5, 10, 20 and 50, also Rayleigh numbers of 103
to 104
keeping fixed Reynolds number 50, where the stream line of flow and thermal fields as well as characteristics of
heat transfer process particularly its expansion has been evaluated. The results of the numerical analysis lead to
the following conclusions:
The MHD has significant effect on the flow and thermal distributions in the triangular enclosures. The
average Nusselt number (Nu) at the hot wall is the highest for the Pr = 6.0 when Rayleigh number 104
,
whereas the lowest heat transfer rate for the Ha = 50 when Rayleigh number 103
. Moreover, the average
Nusselt number with MHD effect is considered higher than those obtained without MHD effect for
different value of Pr number. However, the average bulk temperature at the hot wall is the lowest for
the Pr = 6.0 when Rayleigh number 104
, whereas the highest rate for the Ha = 50 when Rayleigh
number 103
The heat transfer rate increases steadily with the increase of Ra due to increasing of Ra number refers
that the fluid flow switch from conductive to convective heat transfer. The bulk Temperature decreases
smoothly with the increase of Ra because at constant conductivity bulk temperature is inversely
proportional to the Nu (Average Heat Transfer rate).
The heat transfer rate augments with the increase of Ra because when Pr number increases, the
convective heat transfer is dominated inside the flow field. The bulk Temperature remains linear at Pr =
0.71 because low Pr number conduction heat transfer dominates in the enclosures. But higher Pr bulk
temperature decreases as conduction turns into convection.
With the increase of Ha number electromagnetic conduction increases that slow down the heat transfer
rate though it increases at forced convection dominated region for a while but fall steadily in mixed
convection dominated region.
ACKNOWLEDGEMENTS
We would like to express our gratitude to the Department of Arts and Sciences, Ahsanullah University of
Science and Technology (AUST), Dhaka, Bangladesh for providing computing facility during this work.
NOMENCLATURE
θav Average temperature
B0 Magnetic Induction
Cp Specific heat at constant pressure (J/kg K)
G Gravitational acceleration (m/s2
)
Gr Grashof number
H Convective heat transfer coefficient (W/m2
K)
k Thermal conductivity of fluid(W/m K)
K Tthermal conductivity ratio fluid
N non-dimensional distance
Nu Average Nusselt number
Nulocal Local Nusselt number
P Non-dimensional pressure
p pressure
Pr Prandtl number
Ra Rayleigh number
Re Reynolds number
T Non-dimensional temperature
U Dimensionless horizontal velocity
u velocity in x-direction (m/s)
V dimensionless vertical velocity
v velocity in y-direction (m/s)
V0 Lid velocity
x, y Cartesian coordinates
X, Y dimensionless Cartesian coordinates
Greek symbols
β Coefficient of thermal expansion (K-1
)
ρ Density of the fluid (kg/m3
)
α Thermal diffusivity (m2
/s)
∆θ Temperature difference
θ Fluid temperature
µ Dynamic viscosity of the fluid (Pa s)
ν Kinematic viscosity of the fluid (m2
/s)
σ Fluid electrical conductivity(Ω-1
m-1
)

Vol.3, No.5, 2013
67
REFERENCES
Chamkha, A. J. (2002). Hydromagnetic combined convection flow in a vertical lid-driven enclosures with
internal heat generation or absorption, Numer. Heat Transfer, Part A, Vol. 41, pp. 529-546.
Dechaumphai P. (1999). Finite Element Method in Engineering, 2nd ed. (Chulalongkorn Univ. Press, Bangkok,).
Oreper, G. M., Szekely, J. (1983). The effect of an externally imposed magnetic field on buoyancy driven flow in
a rectangular enclosures, J. of Crystal Growth, Vol. 64, pp. 505-515.
Ozoe, H., and Maruo, M. (1987), Magnetic and gravitational natural convection of melted silicon-two
dimensional numerical computations for the rate of heat transfer, JSME, Vol. 30, pp. 774-784.
Rahman, M. M., Alim, M. A. (2010). MHD mixed convection flow in a vertical lid-driven square enclosure
including a heat conducting horizontal circular cylinder with Joule heating, Nonlinear Analysis: Modeling and
Control, Vol. 15(2), pp. 199-211.
Rahman, M. M., Alim, M. A., Chowdhury, M. K. (2009b), Megnetohydrodynamic mixed convection around a
heat conducting horizontal circular cylinder in a rectangular liddriven enclosures with Joule heating, J. Sci Res.
1 (3), pp. 461-472.
Rahman, M.M., Mamun M.A.H., Saidur, R., and Nagata Shuichi (2009a). Effect Of a Heat Conducting
Horizontal Circular Cylinder on MHD Mixed Convection in A Lid- Driven Enclosures Along With Joule Heating,
Int. J. of Mechanical and Materials Eng. (IJMME), Vol. 4, No. 3, pp.256-265.
Reddy J.N. (1985). An Introduction to the Finite Element Method, McGraw-Hill, New York, 1985.
Rudraiah, N., Barron, R. M. (1995b). Venkatachalappa, M., and Subbaraya, C. K., “Effect of magnetic field on
free convection in a rectangular enclosure, Int. J. Engng. Sci., Vol. 33, pp. 1075-1084.
Rudraiah, N., Venkatachalappa, M., and Subbaraya, C. K. (1995a), Combined surface tension and buoyancy-
driven convection in a rectangular open enclosures in the presence of magnetic field, Int. J. Non-linear Mech.,
Vol. 30(5), pp. 759-770.
Singh, S., and Sharif, M. A. R. (2003). Mixed convective cooling of a rectangular cavity with inlet and exit
openings on differentially heated side walls, Numer. Heat Transfer, Part A Vol. 44, pp. 233–253.
Taylor C. and Hood P. (1973). Computer and Fluids, 1 (1), 73 .http://dx.doi.org/10.1016/0045-7930(73)90027-3

Effects of conduction on magneto hydrodynamics mixed convection flow in triangular enclosures

More Related Content

What's hot (16)

Viewers also liked (13)

Similar to Effects of conduction on magneto hydrodynamics mixed convection flow in triangular enclosures (20)

More from Alexander Decker (20)

Recently uploaded (20)

Effects of conduction on magneto hydrodynamics mixed convection flow in triangular enclosures