# MHD Two-Fluid Flow and Heat Transfer between Two Inclined Parallel Plates in a Rotating System

Authors: P. Sri Ramachandra Murty; G. Balaji Prakash
Journal: International Scholarly Research Notices
 (2014)

---

## Abstract

Two-phase magnetohydrodynamic convective flow of electrically conducting fluid through an inclined channel is studied under the action of a constant transverse magnetic field in a rotating system. The fluids in the two phases are steady, incompressible, laminar, immiscible, and electrically conducting, having different densities, viscosities, and thermal and electrical conductivities. The transport properties of both the fluids are assumed constant. The bounding infinite inclined parallel plates are maintained at different constant temperatures, making an angleϕ with the horizontal. Approximate solutions for velocity and temperature distributions are obtained by using a straightforward regular perturbation technique. An in-depth study has been done on the effects of rotation parameter, Hartmann number, inclination angle, the ratio of electrical conductivities, and viscosities of two fluids on the flow. It is observed that the effect of increasing rotation is to decrease the primary velocity. Further it is noticed that as the rotation increases, the secondary velocity increases for smaller rotation, while for larger rotation it decreases. It is also found that the temperature distribution decreases as the rotation increases.


---

## Body

## 1. Introduction

The study of convective Hartmann flow and heat transfer between two parallel plates is receiving considerable interest in the current literature. Shail [1] studied the problem of Hartmann flow of a conducting fluid in a horizontal channel of insulated plates with a layer of nonconducting fluid overlying a conducting fluid in a two-fluid flow. Rudraiah et al. [2] and Umavathi [3] have done a detailed analysis on free and forced convective heat transfer to an electrically conducting fluid in a channel. Lohrasbi and Sahai [4], Malashetty and Leela [5, 6], and Raju and Murty [7] have studied the Hartmann flow characteristics of two fluids in a horizontal channel. Seth et al. [8] have presented the detailed analysis on Hartmann flow in a rotating system in the presence of inclined magnetic field with Hall effects. Chauhan and Rastogi [9, 10] have considered Hall current and heat transfer effects on MHD flow and MHD Couette flow in a channel partially filled with a porous medium in a rotating system.Numerous publications dealing with both the experimental and the theoretical aspects of the two-phase flow systems with or without considering the heat transfer problems associated with MHD power generators, MHD devices, and thermonuclear power generations have appeared in the literature. Abdul Mateen [11, 12] has studied the magnetohydrodynamic flow and transient magnetohydrodynamic flow of two immiscible fluids through a horizontal channel. Raju and Nagavalli [13] have studied the MHD two-layered unsteady flow and heat transfer through a horizontal channel in the presence of applied magnetic and electric fields in a rotating system. Recently, Murty and Linga Raju [14] have investigated magnetohydrodynamic two-phase flow and heat transfer between two parallel porous walls in a rotating system.Basically, the inclined geometry has enormous applications in the heat transfer technology like solar collector. Though such problems with inclined geometry are close to realistic practical situations, much attention has not been given to handling them except the studies by Malashetty and Umavathi [15], Malashetty et al. [16], S. Daniel and Y. S. Daniel [17], and Murty and Prakash [18]. Therefore, in this paper, we have studied two-phase magnetohydrodynamic convective flow of electrically conducting fluid through an inclined channel under the action of a constant transverse magnetic field when rotated by an angular velocity about an axis perpendicular to the plates.


## 2. Mathematical Formulation

A steady laminar and fully developed two-phase magnetohydrodynamic convective flow driven by a common constant pressure gradient-
∂
p
/
∂
x and temperature gradient Δ
T
=
T
w
1
-
T
w
2 has been considered in the presence of a constant magnetic field applied transversely to the direction of the flow. The physical configuration is shown in Figure 1, which consists of two infinite inclined parallel plates maintained at different constant temperatures, extending in the x and z directions, making an angle ϕ with the horizontal. The regions 0
≤
y
≤
h
1 and -
h
2
≤
y
≤
0 are occupied by two different electrically conducting incompressible fluids having density ρ
i, viscosity μ
i, electrical conductivity σ
i, and thermal conductivity K
i. The whole system is rotated with an angular velocity Ω in a counterclockwise direction about y-axis perpendicular to the plates. The transport properties of both the fluids are assumed constant. The suffix i  (i
=
1,2) represents the values for phases I and II, respectively.Figure 1
Physical configuration.With these assumptions, the governing equations of motion and energy for Boussinesq fluids as in Malashetty et al. [16] for both phases are
(1)
μ
i
d
2
u
i
d
y
2
+
ρ
i
g
β
i
T
i
-
T
w
2
sin
⁡
ϕ
-
σ
i
B
0
2
u
i
=
∂
p
∂
x
+
2
Ω
ρ
i
w
i
,
μ
i
d
2
w
i
d
y
2
-
σ
i
B
0
2
w
i
=
-
2
Ω
ρ
i
u
i
,
d
2
T
i
d
y
2
+
μ
i
K
i
d
u
i
d
y
2
+
d
w
i
d
y
2
+
σ
i
B
0
2
u
i
2
+
w
i
2
K
i
=
0
,
where u
i and w
i are the primary and secondary velocity components along x and z directions, respectively, T
i is the temperature, β
i is the coefficient of thermal expansion, and g is the acceleration due to gravity. The fluid and the thermometric boundary conditions are unchanged by the addition of electromagnetic field. The no-slip condition requires that the velocity must be vanishing at the wall. In addition, the fluid velocity, shear stress, temperature, and heat flux must be continuous across the interface.The boundary and interface conditions are(2)
u
1
h
1
=
0
,
w
1
h
1
=
0
;
u
1
0
=
u
2
0
,
w
1
0
=
w
2
0
;
u
2
-
h
2
=
0
,
w
2
-
h
2
=
0
,
μ
1
d
u
1
d
y
=
μ
2
d
u
2
d
y
,
μ
1
d
w
1
d
y
=
μ
2
d
w
2
d
y
,
at
y
=
0
,
T
1
h
1
=
T
w
1
,
T
1
0
=
T
2
0
,
T
2
-
h
2
=
T
w
2
,
K
1
d
T
1
d
y
=
K
2
d
T
2
d
y
at
y
=
0
.
In making these equations dimensionless, the following transformations are used:
(3)
u
i
*
=
u
i
u
-
1
,
y
i
*
=
y
i
h
i
,
θ
=
T
-
T
w
2
Δ
T
,
m
=
μ
1
μ
2
,
K
=
K
1
K
2
,
h
=
h
2
h
1
,
n
=
ρ
2
ρ
1
,
b
=
β
2
β
1
,
s
=
σ
2
σ
1
,
Gr
=
g
β
1
h
1
3
Δ
T
ν
1
2
,
M
=
B
0
h
1
σ
1
μ
1
,
Pr
⁡
=
μ
1
C
p
K
1
,
Ec
=
u
-
1
2
C
p
Δ
T
,
Re
=
u
-
1
h
1
ν
1
,
P
=
h
1
2
μ
1
u
-
1
∂
p
∂
x
,
R
2
=
Ω
h
1
2
ν
.
With the above nondimensional quantities, the governing equations (1) become
(4)
d
2
u
i
d
y
2
+
Gr
Re
A
sin
⁡
ϕ
θ
i
-
B
M
2
u
i
=
CP
+
2
R
2
w
i
,
d
2
w
i
d
y
2
-
B
M
2
w
i
=
-
2
R
2
u
i
,
d
2
θ
i
d
y
2
+
PrEc
D
d
u
i
d
y
2
+
d
w
i
d
y
2
+
PrEc
F
M
2
u
i
2
+
w
i
2
=
0
,
where A
=
b
m
n
h
2,  B
=
m
s
h
2,  C
=
m
h
2,  D
=
(
K
/
m
),  F
=
K
h
2
s , and A, B, C, D, and F are all equal to 1 for phase I.The nondimensional forms of the boundary and interface conditions (2) become
(5)
u
1
1
=
0
,
w
1
1
=
0
;
u
1
0
=
u
2
0
,
w
1
0
=
w
2
0
;
u
2
-
1
=
0
,
w
2
-
1
=
0
,
d
u
1
d
y
=
1
m
h
d
u
2
d
y
,
d
w
1
d
y
=
1
m
h
d
w
2
d
y
at
y
=
0
,
θ
1
1
=
1
,
θ
1
0
=
θ
2
0
,
θ
2
-
1
=
0
,
d
θ
1
d
y
=
1
K
h
d
θ
2
d
y
at
y
=
0
.The asterisks have been dropped for simplicity. Further writingq
1
=
u
1
+
i
w
1 and q
2
=
u
2
+
i
w
2, (4) can be written in complex form as
(6)
d
2
q
i
d
y
2
+
Gr
Re
A
sin
⁡
ϕ
θ
i
-
B
M
2
q
i
=
CP
-
2
i
R
2
q
i
,
(7)
d
2
θ
i
d
y
2
+
PrEc
D
d
q
i
d
y
d
q
-
i
d
y
+
PrEc
F
M
2
q
i
q
-
i
=
0
which are to be solved subject to the boundary and interface conditions:
(8)
q
1
1
=
0
,
q
1
0
=
q
2
0
,
q
2
-
1
=
0
,
d
q
1
d
y
=
1
m
h
d
q
2
d
y
at
y
=
0
,
θ
1
1
=
1
,
θ
1
0
=
θ
2
0
,
θ
2
-
1
=
0
,
d
θ
1
d
y
=
1
K
h
d
θ
2
d
y
at
y
=
0
.


## 3. Solutions

The governing equations of momentum (6) along with the energy equations (7) are to be solved subject to the boundary and interface conditions (8) for the velocity and temperature distributions. Here, we consider the Eckert number very small. Hence, the product PrEc (=
ε) is very small and can be exploited to use the regular perturbation method. The solutions are assumed in the form
(9)
q
i
,
θ
i
=
q
i
0
,
θ
i
0
+
ε
q
i
1
,
θ
i
1
+
⋯
,
where q
i
0
,
θ
i
0 are solutions for the case ε equal to zero. q
i
1
,
θ
i
1 are perturbed quantities relating to q
i
0
,
θ
i
0, respectively. Substituting the above solution in (6) and (7) and equating the coefficients of similar powers of ε to zero, we get the zeroth and the first order equations as follows.Zeroth order equations:(10)
d
2
q
i
0
d
y
2
+
Gr
Re
A
sin
⁡
ϕ
θ
i
0
-
B
M
2
q
i
0
=
CP
-
2
i
R
2
q
i
0
,
d
2
θ
i
0
d
y
2
=
0
.
First order equations:
(11)
d
2
q
i
1
d
y
2
+
Gr
Re
A
sin
⁡
ϕ
θ
i
1
-
B
M
2
q
i
1
=
-
2
i
R
2
q
i
1
,
d
2
θ
i
1
d
y
2
+
D
d
q
i
0
d
y
d
q
-
i
0
d
y
+
F
M
2
q
i
0
q
-
i
0
=
0
.
The corresponding Boundary conditions (8) reduce to
(12)
q
10
1
=
0
,
q
10
0
=
q
20
0
,
q
20
-
1
=
0
,
d
q
10
d
y
=
1
m
h
d
q
20
d
y
at
y
=
0
,
θ
10
1
=
1
,
θ
10
0
=
θ
20
0
,
θ
20
-
1
=
0
,
d
θ
10
d
y
=
1
K
h
d
θ
20
d
y
at
y
=
0
,
(13)
q
11
1
=
0
,
q
11
0
=
q
21
0
,
q
21
-
1
=
0
,
d
q
11
d
y
=
1
m
h
d
q
21
d
y
at
y
=
0
,
θ
11
1
=
0
,
θ
11
0
=
θ
21
0
,
θ
21
-
1
=
0
,
d
θ
11
d
y
=
1
K
h
d
θ
21
d
y
at
y
=
0
.
We observe that (10) and (11) are linear and coupled and therefore can be solved exactly. Here, we consider q
10
=
u
10
+
i
w
10,  q
20
=
u
20
+
i
w
20,  q
11
=
u
11
+
i
w
11, and q
21
=
u
21
+
i
w
21.Solutions of the zeroth order equations (10) using boundary conditions (12) are
(14)
θ
10
=
y
+
K
h
1
+
K
h
,
θ
20
=
1
+
y
K
h
1
+
K
h
,
(15)
u
10
=
c
1
e
A
5
y
+
c
2
e
-
A
5
·
y
cos
⁡
A
6
y
+
A
14
+
A
12
y
,
w
10
=
-
c
1
e
A
5
y
-
c
2
e
-
A
5
·
y
sin
A
6
y
-
A
15
-
A
13
y
,
u
20
=
c
3
e
A
23
y
+
c
4
e
-
A
23
·
y
cos
⁡
A
24
y
+
A
30
+
A
32
y
,
w
20
=
-
c
3
e
A
23
y
-
c
4
e
-
A
23
·
y
sin
⁡
A
24
y
-
A
31
-
A
33
y
.
Solutions of the first order equations (11) using boundary conditions (13) are
(16)
θ
11
=
c
5
y
+
c
6
+
B
56
e
2
A
5
y
+
B
57
cos
⁡
⁡
2
A
6
y
hhh
+
B
58
e
-
2
A
5
y
+
B
132
e
A
5
y
sin
⁡
A
6
y
+
B
98
y
2
hhh
+
B
133
e
A
5
y
cos
⁡
⁡
A
6
y
+
B
134
e
-
A
5
y
sin
⁡
A
6
y
hhh
+
B
135
e
-
A
5
y
cos
⁡
⁡
A
6
y
+
B
136
e
A
5
y
y
cos
⁡
⁡
A
6
y
hhh
+
B
137
e
A
5
y
y
sin
⁡
A
6
y
+
B
138
e
-
A
5
y
y
cos
⁡
⁡
A
6
y
hhh
+
B
139
e
-
A
5
y
y
sin
⁡
A
6
y
+
B
96
y
3
+
B
97
y
4
+
i
B
99
sin
⁡
2
A
6
y
+
B
140
e
A
5
y
sin
⁡
A
6
y
hhhhh
+
B
141
e
A
5
y
cos
⁡
⁡
A
6
y
+
B
142
e
-
A
5
y
sin
⁡
A
6
y
hhhhh
+
B
143
e
-
A
5
y
cos
⁡
⁡
A
6
y
+
B
144
e
A
5
y
y
cos
⁡
A
6
y
hhhhh
+
B
145
e
A
5
y
y
sin
⁡
A
6
y
+
B
146
e
-
A
5
y
y
cos
⁡
A
6
y
h
h
h
h
h
+
B
147
e
-
A
5
y
y
sin
⁡
A
6
y
,
θ
21
=
c
7
y
+
c
8
+
E
67
e
2
A
23
y
+
E
68
cos
⁡
⁡
2
A
24
y
hhh
+
E
69
e
-
2
A
23
y
+
E
135
e
A
23
y
sin
⁡
A
24
y
hhh
+
E
136
e
A
23
y
cos
⁡
⁡
A
24
y
+
E
137
e
-
A
23
y
sin
⁡
A
24
y
hhh
+
E
138
e
-
A
23
y
cos
⁡
⁡
A
24
y
-
E
107
y
4
-
E
108
y
3
hhh
+
E
139
e
A
23
y
y
cos
⁡
⁡
A
24
y
+
E
140
e
A
23
y
y
sin
⁡
A
24
y
hhh
+
E
141
e
-
A
23
y
y
cos
⁡
⁡
A
24
y
-
E
109
y
2
h
h
h
+
E
142
e
-
A
23
y
y
sin
⁡
A
24
y
+
i
E
110
sin
⁡
2
A
24
y
+
E
143
e
A
23
y
sin
⁡
A
24
y
hhhhh
+
E
144
e
A
23
y
cos
⁡
⁡
A
24
y
+
E
145
e
-
A
23
y
sin
⁡
A
24
y
hhhhh
+
E
146
e
-
A
23
y
cos
⁡
A
24
y
+
E
147
e
A
23
y
y
cos
⁡
⁡
A
24
y
hhhhh
+
E
148
e
A
23
y
y
sin
⁡
A
24
y
+
E
149
e
-
A
23
y
y
cos
⁡
⁡
A
24
y
hhhhh
+
E
150
e
-
A
23
y
y
sin
⁡
A
24
y
,
u
11
=
c
9
e
A
5
y
+
c
10
e
-
A
5
y
cos
⁡
⁡
A
6
y
+
F
95
e
2
A
5
y
+
F
96
e
-
2
A
5
y
-
F
97
cos
⁡
⁡
2
A
6
y
+
F
98
sin
⁡
2
A
6
y
+
F
99
e
A
5
y
y
sin
⁡
A
6
y
+
F
100
e
A
5
y
y
cos
⁡
⁡
A
6
y
-
F
101
e
-
A
5
y
y
cos
⁡
⁡
A
6
y
-
F
107
y
4
-
F
108
y
3
+
F
102
e
-
A
5
y
y
sin
⁡
A
6
y
+
F
103
e
A
5
y
y
2
sin
⁡
A
6
y
+
F
104
e
A
5
y
y
2
cos
⁡
⁡
A
6
y
-
F
109
y
2
-
F
110
y
+
F
105
e
-
A
5
y
y
2
sin
⁡
A
6
y
-
F
106
e
-
A
5
y
y
2
cos
⁡
⁡
A
6
y
-
F
111
,
w
11
=
c
10
e
-
A
5
y
-
c
9
e
A
5
y
sin
⁡
A
6
y
-
F
112
e
2
A
5
y
-
F
113
e
-
2
A
5
y
-
F
114
cos
⁡
⁡
2
A
6
y
-
F
115
sin
⁡
2
A
6
y
+
F
116
e
A
5
y
y
cos
⁡
⁡
A
6
y
+
F
117
e
A
5
y
y
sin
⁡
A
6
y
-
F
125
y
3
-
F
126
y
2
-
F
118
e
-
A
5
y
y
cos
⁡
⁡
A
6
y
+
F
119
e
-
A
5
y
y
sin
⁡
A
6
y
+
F
120
e
A
5
y
y
2
sin
⁡
A
6
y
-
F
127
y
-
F
128
+
F
121
e
A
5
y
y
2
cos
⁡
⁡
A
6
y
+
F
122
e
-
A
5
y
y
2
sin
⁡
A
6
y
-
F
123
e
-
A
5
y
y
2
cos
⁡
⁡
A
6
y
-
F
124
y
4
,
u
21
=
c
11
e
A
23
y
+
c
12
e
-
A
23
y
cos
⁡
⁡
A
24
y
+
G
96
e
2
A
23
y
+
G
97
e
-
2
A
23
y
+
G
98
cos
⁡
⁡
2
A
24
y
+
G
112
-
G
99
sin
⁡
2
A
24
y
+
G
100
e
A
23
y
y
sin
⁡
A
24
y
+
G
101
e
A
23
y
y
cos
⁡
⁡
A
24
y
+
G
108
y
4
+
G
109
y
3
+
G
102
e
-
A
23
y
y
cos
⁡
⁡
A
24
y
+
G
103
e
-
A
23
y
y
sin
⁡
A
24
y
+
G
104
e
A
23
y
y
2
sin
⁡
A
24
y
+
G
110
y
2
+
G
105
e
A
23
y
y
2
cos
⁡
⁡
A
24
y
+
G
106
e
-
A
23
y
y
2
cos
⁡
⁡
A
24
y
+
G
107
e
-
A
23
y
y
2
sin
⁡
A
24
y
+
G
111
y
,
w
21
=
c
12
e
-
A
23
y
-
c
11
e
A
23
y
sin
⁡
A
24
y
-
G
113
e
2
A
23
y
-
G
114
e
-
2
A
23
y
+
G
115
cos
⁡
⁡
2
A
24
y
+
G
129
+
G
116
sin
⁡
2
A
24
y
+
G
117
e
A
23
y
y
sin
⁡
A
24
y
+
G
118
e
A
23
y
y
cos
⁡
⁡
A
24
y
+
G
128
y
+
G
127
y
2
+
G
119
e
-
A
23
y
y
cos
⁡
⁡
A
24
y
+
G
120
e
-
A
23
y
y
sin
⁡
A
24
y
+
G
121
e
A
23
y
y
2
sin
⁡
A
24
y
+
G
126
y
3
+
G
122
e
A
23
y
y
2
cos
⁡
⁡
A
24
y
+
G
123
e
-
A
23
y
y
2
sin
⁡
A
24
y
+
G
124
e
-
A
23
y
y
2
cos
⁡
⁡
A
24
+
G
125
y
4
.
The constants appearing in (15) and (16) are not given for the sake of brevity. Since the problem contains too many nondimensional parameters, for the sake of conciseness, we fix P
=
-
5.0, b
=
1.0, Re
=
5.0, n
=
1.5, and K
=
1.0. In the figures, all the other parameters except the varying one are chosen from the set (M, Gr, ϕ, m, h, s, R) = (2.0, 5.0, 30°, 0.5, 1.0, 2.0, 2.0).


## 4. Results and Discussion

Two-phase magnetohydrodynamic convective flow between two infinite inclined parallel plates in a rotating system is studied analytically. The resulting differential equations are solved using perturbation method to obtain approximate solutions for temperature distribution and primary and secondary velocity distributions. Here, we note that whenR
=
0, that is, in the absence of rotation, these results are in agreement with that of Malashetty et al. [16].The effect of rotation parameterR on primary velocity u and secondary velocity w is shown in Figures 2 and 3, respectively. From Figure 2, it is observed that the primary velocity u decreases with the increase in the rotation parameter. The rotation parameter R defines the relative magnitude of the Coriolis force and the viscous force in the regime. As the high magnitude Coriolis forces oppose the buoyancy force, the velocity will be decreased. From Figure 3, it is concluded that as the rotation parameter R increases in (0, 2), the secondary velocity w also increases but outside this range as R increases, it decreases. Therefore, by increasing the rotation parameter R, the secondary flow becomes oscillatory.Figure 2
Primary velocity profiles for different values of rotation parameterR.Figure 3
Secondary velocity profiles for different values of rotation parameterR.Figures4 and 5 show the effect of the angle of inclination ϕ on primary velocity u and secondary velocity w, respectively. As the angle of inclination ϕ increases, both the primary and secondary velocities increase because the magnitude of the buoyancy force increases with increase in the inclination angle. The effect of the ratio of the viscosities m on primary and secondary velocities is shown in Figures 6 and 7, respectively. The smaller the value of the viscosity of the fluid in the lower phase compared to the fluid in the upper phase, the larger the primary and secondary flow fields.Figure 4
Primary velocity profiles for different values of inclination angleϕ.Figure 5
Secondary velocity profiles for different values of inclination angleϕ.Figure 6
Primary velocity profiles for different values of ratio of viscositiesm.Figure 7
Secondary velocity profiles for different values of ratio of viscositiesm.Figures8 and 9 represent the effect of the ratio of heights h on primary and secondary velocities, respectively. The smaller the height of the upper phase compared to the lower phase, the larger the primary as well as secondary flow fields. The effect of the ratio of the electrical conductivity s on primary velocity u and secondary velocity w is shown in Figures 10 and 11, respectively. We have observed that as the ratio of the electrical conductivity s increases, the primary velocity u increases but the secondary velocity w decreases.Figure 8
Primary velocity profiles for different values of ratio of heightsh.Figure 9
Secondary velocity profiles for different values of ratio of heightsh.Figure 10
Primary velocity profiles for different values of ratio of electrical conductivitiess.Figure 11
Secondary velocity profiles for different values of ratio of electrical conductivitiess.Figures12 and 13 show the effect of Hartmann number M on primary and secondary velocities, respectively. The effect of increasing Hartmann number M is to decrease both the primary and secondary velocities. This is because an increase in applied magnetic field strength causes greater interaction between the fluid motion and the magnetic field, therefore, an increase in the Lorentz force. Since this force opposes the buoyancy force, both the velocities will be decreased. The effect of Grashof number Gr on primary and secondary velocities is shown in Figures 14 and 15, respectively. We find that an increase in the value of Grashof number Gr increases both the primary and secondary velocities.Figure 12
Primary velocity profiles for different values of Hartmann numberM.Figure 13
Secondary velocity profiles for different values of Hartmann numberM.Figure 14
Primary velocity profiles for different values of Grashof number Gr.Figure 15
Secondary velocity profiles for different values of Grashof number Gr.The effect of rotation parameter on temperatureθ can be seen in Figure 16. From the figure, it is evident that the temperature decreases with the increase in the rotation parameter R. The rotation parameter R  R
2
=
Ω
h
i
2
/
ν defines the relative magnitude of the Coriolis force and the viscous force in the regime. As the high magnitude Coriolis forces oppose the buoyancy force, the velocity will be decreased leading to a reduction in the viscous and Joule dissipation and so to a reduction in the temperature.Figure 16
Temperature profiles for different values of rotation parameterR.Figure17 shows the effect of the angle of inclination ϕ on temperature θ. As the angle of inclination ϕ increases, the temperature θ also increases because the magnitude of the buoyancy force increases with the increase in the inclination angle. The effect of the ratio of viscosities m on the temperature θ is shown in Figure 18. From the figure, it is observed that the less viscous fluid in the lower phase adds the heat transfer. Figure 19 exhibits the effect of the ratio of heights h on the temperature θ. From the figure, it is noticed that the smaller the height of the upper phase compared to the lower phase, the larger the magnitude of the temperature. Figure 20 shows the effect of the ratio of electrical conductivities s on the temperature θ. From the figure, it is concluded that as the ratio of electrical conductivities s increases, the temperature also increases.Figure 17
Temperature profiles for different values of inclination angleϕ.Figure 18
Temperature profiles for different values of ratio of viscositiesm.Figure 19
Temperature profiles for different values of ratio of heightsh.Figure 20
Temperature profiles for different values of ratio of electrical conductivitiess.Figure21 represents the effect of Hartmann number M on the temperature θ. From the figure, it is clear that the effect of increasing M is to decrease the temperature. This is because an increase in the applied magnetic field strength causes greater interaction between the fluid motion and the magnetic field, therefore, an increase in the Lorentz force. Since this force opposes the buoyancy force, the temperature will be decreased. The effect of Grashof number Gr on temperature θ is shown in Figure 22; we observe from this figure that an increase in the value of Grashof number Gr increases the temperature θ.Figure 21
Temperature profiles for different values of Hartmann numberM.Figure 22
Temperature profiles for different values of Grashof number Gr.



---

*Source: 256898-2014-12-17.xml*