Mathematical modelling and analysis of three dimensional darcy

International Journal of Mechanical Engineering and Technology (IJMET), ISSN 0976 –
6340(Print), ISSN 0976 – 6359(Online) Volume 4, Issue 2, March - April (2013) © IAEME
562
MATHEMATICAL MODELLING AND ANALYSIS OF THREE
DIMENSIONAL DARCY –BRINKMAN (D-B) MODEL IN AN
INCLINED RECTANGULAR POROUS BOX
Dr. R. P. Sharma
Dept. of Mechanical Engineering, Birla Institute of Technology, Mesra, Ranchi, 835215
India
ABSTRACT
In this paper, numerical studies on three- dimensional natural convection in an
inclined differentially heated porous box employing Darcy-Brinkman- flow model are
presented. The relative effects of inertia and viscous forces on natural convection in porous
media are examined. When Brinkman viscous terms are included in the buoyancy term of the
momentum equation, no-slip conditions for velocity at the walls are satisfied. The flow and
temperature fields become three-dimensional due to Brinkman viscous term. The governing
equations for the present studies are obtained by setting Da≠0 and Fc/Pr = 0 in the general
governing equations for D-B flow description. The system is characterized by the non-
dimensional parameters: Rayleigh number (Ra), vertical aspect ratio (ARY), horizontal aspect
ratio (ARZ), Darcy number (Da) and the angle of inclination (φ). Numerical solutions have
been obtained employing the SAR scheme for 200 < Ra < 2000, 0.2 < ARY< 5, 0.2 ARZ < 5,
10-5
< Da < 10-2
and -60o
<φ< 60o
.
Keywords: Brinkman, buoyancyterm, no-slip conditions, viscous term etc
1.0 INTRODUCTION
Studies on natural convection heat transfer in porous media employing non-Darcy
extensions to describe the fluid flow which include Forchheimer non-linear inertial terms,
convective terms and Brinkman viscous terms, have been reported in the literature. Ghanet,
al. [1] included the viscous terms due to Brinkman as an extension to the Darcy equation to
describe the fluid flow. They concluded that the average Nusselt number shows a maximum
when the aspect ratio is around unity.Tong and Subramaniam [2] developed boundary layer
solutions to Brinkman extended Darcy flow model based on the modified Oseen technique
and the flow field is found to be governed by a new parameter. Tong and Subramaniam found
INTERNATIONAL JOURNAL OF MECHANICAL ENGINEERING
AND TECHNOLOGY (IJMET)
ISSN 0976 – 6340 (Print)
ISSN 0976 – 6359 (Online)
Volume 4, Issue 2, March - April (2013), pp. 562-567
© IAEME: www.iaeme.com/ijmet.asp
Journal Impact Factor (2013): 5.7731 (Calculated by GISI)
www.jifactor.com
IJMET
© I A E M E

563
that a pure Darcy analysis is applicable when the parameter defined in their study is less than
10-4
.Sen [3] considered the influence of Brinkman viscous terms and convective terms in the
limit of vanishingly small aspect ratio. Lauriat and Prasad [4] considered the Brinkman
extended Darcy flow model along with convective terms for AR > 1. The numerical study [4]
concluded that the influence of the Brinkman terms is significant for high Ra and Darcy
number and the influence of convective terms on the Nusselt number is insignificant.Kwok
and Chen [5] observed from their experimental investigation that the flow becomes unstable
at a critical Rayleigh number of 66.2. They performed a linear stability analysis to study the
effects of including the Brinkman terms and variable viscosity separately. The study
concluded that the effect of variable viscosity has a profound influence on the stability of the
thermal convection, whereas the effect of no-slip condition at the wall has minimal influence
on the stability.Chang and Hsiao [6] investigated the influence of anisotropic permeability
and thermal conductivity on natural convection in a vertical cylinder filled with anisotropic
porous material.Holst [7] numerically solved transient, three-dimensional natural convection
in a porous box. They concluded that, compared to the two-dimensional values, the three-
dimensional heat flow under certain instances is higher at high Rayleigh numbers. At low
Rayleigh numbers, the situation may be reversed depending on the initial conditions. Horne
[8] investigated the tendency of the flow to be two or three-dimensional. Horne concluded
that the flow pattern is more complicated than thought previously.In a numerical study,
Singh, et, al., [9] examined the influence of the Brinkman-extended Darcy model on the fluid
flow and heat transfer in a confined fluid overlying a porous layer.Vasseur, et al., [10] studied
analytically and numerically the thermally driven flow in a thin, inclined, rectangular cavity
filled with a fluid saturated porous layer.Sharma R P & Sharma R V has worked on
modelling&simulation of three –dimensional natural convection in a porous box and
concluded that three-dimensional average Nusselt values are lower than two-dimensional
values. [11]
2.0 MATHEMATICAL MODELLING
2.1 Governing Equation
Fig. 1 Physical model and co-ordinate system

564
The physical model is shown in fig.1 is a parallelepiped box of length L, width B and height
H filled with fluid saturated porous medium.
Governing equations for natural convection in the porous box comprising of conservation of
mass, momentum and energy are as follows:
0=
∂
∂
+
∂
∂
+
∂
∂
z
w
y
v
x
u
(1)






∂
∂
+
∂
∂
+
∂
∂
+





+
∂
∂
−=
′
+ 2
2
2
2
2
2
sin
z
u
y
u
x
u
g
x
p
uv
K
K
u
K
µφρρ
µ
(2)






∂
∂
+
∂
∂
+
∂
∂
+





+
∂
∂
−=
′
+ 2
2
2
2
2
2
cos
z
v
y
v
x
v
g
y
p
vv
K
K
v
K
µφρρ
µ
(3)






∂
∂
+
∂
∂
+
∂
∂
+





∂
∂
−=
′
+ 2
2
2
2
2
2
z
w
y
w
x
w
z
p
wv
K
K
w
K
µρ
µ
(4)
u 





∂
∂
+
∂
∂
+
∂
∂
=
∂
∂
+
∂
∂
+
∂
∂
2
2
2
2
2
2
z
T
y
T
x
T
z
T
w
y
T
v
x
T
α (5)
Governing equations are rendered dimensionless introducing the following non-dimensional
variables:
L
x
X =
L
y
Y =
L
z
Z =
L
u
U
/α
=
L
v
V
/α
=
L
w
W
/α
=
ch
c
TT
TT
−
−
=θ
K
p
P
/µα
=
mρ
ρ
ρ = (6)
The dimensionless conservation of mass, momentum and energy are derived and non-
dimensional parameters Ra, the Rayleigh number, Fc, the Forchheimer number, Pr, the
Prandtl number and Da, the Darcy number are defined by -
ρ = 1 - β∆T (θ-0.5);
να
β TKgL
Ra
∆
= ;
L
K
Fc
′
= ;
α
v
=Pr ; 2
L
K
Da = (7)
Hydrodynamic Boundary Conditions
(i) With Brinkman terms
U = V = W = 0 at X = 0,1 for 0 ≤ Y ≤ ARY and 0 ≤ Z ≤ ARZ
U = V = W = 0 at Y = 0,ARY for 0 ≤ X ≤ 1 and 0 ≤ Z ≤ ARZ (8)
U = V = W = 0 at Z = 0,ARY for 0 ≤ X ≤ 1 and 0 ≤ Y ≤ ARZ

565
Thermal Boundary Conditions
θ = 0 at X = 0 for 0 ≤ Y ≤ ARY and 0 ≤ Z ≤ ARZ
θ = 1 at X = 1 for 0 ≤ Y ≤ ARY and 0 ≤ Z ≤ ARZ (9)
0
Y
=
∂
θ∂
at Y = 0, ARY for 0 ≤ X ≤ 1 and 0 ≤ Z ≤ ARZ
0
Z
=
∂
θ∂
at Z = 0, ARZ for 0 ≤ X ≤ 1 and 0 ≤ Y ≤ ARY
Where ARY, the vertical aspect ratio and ARZ the horizontal aspect ratio are defined as
ARY = H/L (10)
ARZ = B/L (11)
Boundary conditions on temperature (θ) are the same as given by Eq (9). The average
Nusselt number based on the characteristic length, L of the box is defined as,
k
Lh
Nu = (12)
The average Nusselt number at X = 0 and X = 1 is obtained by using a suitable numerical
method.
RESULTS & DISCUSSION
Variation of average Nusselt number with ARZ for Ra = 1000, ARY =1.0 and Da =
10-2
is shown in Fig. 2. The three-dimensional effects are pronounced for ARZ< 1 even for
inclined porous enclosure.Variation of average Nusselt number (Nu) with angle of inclination
(φ) for Ra = 500, 1000, 2000 and Da upto 10-2
are shown in Fig. 3. From this figure, it is
evident that in all cases the average Nusselt number (Nu) is increasing up to an angle (φ) = -
30o
and then it is decreasing. This is due to the fact that with increasing inclination, the
multicellular fluid motion in the porous matrix gets changed gradually to unicellular one. The
multicellular circular cells (Fig. 4) and buoyancy forces are equally significant upto the
critical angle of inclination (φc). As these increase slowly uptoφcritical, the average Nusselt
number (Nu) shows marginal increase. After critical angle of inclination (≥φc), as the motion
turns to unicelluar one (Fig. 5) and the average Nusselt number (Nu) values start declining
with increasing inclination angle. Multicellular convection augments the heat transferred
through the porous material. The main mechanism of heat exchange is due to the motion of
the fluid in direction perpendicular to the isothermal walls whereas in the unicellular mode
this occurs only in regions close to the isothermal wall (Fig. 5).
In unicellular circulation situation the temperature difference between fluid and hot
surface (i.e. the driving force for heat transfer) is smaller compared to multicellular
circulation situation. Further, in multicellular system the fluid mixing is more efficient. As a
result, average heat transfer rate during multicellular circulation is more than that during
unicellular situation.

566
0 1 2 3 4 5
1.0
1.5
2.0
2.5
3.0
3.5
4.0
4.5
5.0
Nu
ARZ
φ=-60
o
.
φ=-30
o
.
φ=0
o
.
φ=30
o
.
φ=60
o
.
-60 -40 -20 0 20 40 60
1
2
3
4
5
6
Nu
φ
Ra=200
Ra=500
Ra=1000
Ra=2000
Fig. 2: Variation of Nu with ARz for Fig.3: Variation of Nu with ϕ for
Ra=1000,ARy=1.0, Da=10-2 and ϕ=0° Da=10-2, ARy=1.0, and ARz=1.0
Fig. 4: Iso-vector-potential ( ) plot for Fig. 5: Iso-vector-potential ( ) lines for
Ra = 1000, ARY =ARZ =1.0, Ra = 1000, ARY =ARZ =1.0,
Da = 10-2 and φ = -30o. Da = 10-2 and φ = 0o.
0.00 0.10 0.20 0.30 0.40 0.50 0.60 0.70 0.80 0.90 1.00
0.00
0.10
0.20
0.30
0.40
0.50
0.60
0.70
0.80
0.90
1.00
0.00 0.10 0.20 0.30 0.40 0.50 0.60 0.70 0.80 0.90 1.00
0.00
0.10
0.20
0.30
0.40
0.50
0.60
0.70
0.80
0.90
1.00

567
CONCLUSIONS
The flow description is within the framework of Darcy-Brinkman flow model. The
non-dimensional parameters required to describe the system are Ra, ARY, ARZ, Da andφ. The
numerical solutions have been obtained for 100 < a < 2000, 10-5
< Da < 10-2
, 0.2 < ARY,
ARZ< 5 and – 60o
<φ< 60o
.Due to inclusion of Brinkman viscous term, which satisfies the no-
slip boundary condition, the flow becomes three-dimensional particularly for low horizontal
aspect ratio, ARZ< 1. The three-dimensional average Nusselt number values are lower than
that of two-dimensional average Nusselt number values for ARZ< 1. For Da < 10-5
, the
average Nusselt number values are close that of Darcy’s values. The critical angle of
inclination for Darcy – Brinkman flow description is – 30o
irrespective of Darcy number and
Rayleigh number.
REFERENCES
[1] C.M. lvy B.K.C Ghan and J.M. Barry, Natural convection in enclosed porous media with
rectangular boundaries, ASME trans. Journal of Heat Transfer, Vol. 92. pp.21 - 27,
1970.
[2] T.W. Tong and E. Subramanian , A boundary layer analysis for naural convection in
porous enclosure use of Brinkman extended Darcy flow model, International Journal of
Heat and Mass Transfer, Vol. 28 pp.563-571, 1985.
[3] A.K. Sen, Natural convection in a shallow porous cavity- the Brinkman model,
International Journal of Heat and Mass Transfer, Vol. 30 pp.855-869, 1987.
[4] G. Lauriat and V. Prasad, Natural convection in a vertical porous cavity: a numerical
study for Brinkman-extended Darcy formulation, ASME Trans. Journal of Heat Transfer,
Vol.109 pp.688-696, 1987.
[5] L.P.Kwok and C.F. Chen, Stability of thermal convection in a verical porous layer,
ASME Trans, Journal of Heat Transfer, Vol.109 pp.889-892, 1987.
[6] Wen-Jeng Chang and Chi-Feng Hsiao, Natural convection in a vertical cylinder filled
with anistropic porous media, International Journal of Heat Mass Transfer, Vol. 36, pp.
3361-3367, 1993.
[7] P.H. Holst, Transient three-dimensional natural convection in confined porous media,
International Journal of Heat and Mass Transfer, pp. 73-90, 1972.
[8] R.N. Horne, 3-D natural convection in a confined porous medium heated from below,
Journal of Fluid Mechanics, Vol. 92, pp. 751-766, 1979.
[9] A.K.Singh, E. Leonardy and G.R.Thope, Three-dimensional natural convection in a
confined fluid overlaying a porous layer, Trans. ASME Journal of Heat Transfer , Vol.
115, pp. 631-638 , 1993.
[10] P.Vasseur, M.G. Satish and L. Robillard, Natural Convection in a thin, Inclined Porous
Layer exposed to a constant Heat Flux , International Journal of Heat and Mass Transfer ,
Vol. 3 pp.537-549, 1987.
[11] R.P. Sharma, R.V. Sharma, “Modelling & simulation of three-dimensional natural
convection in a porous media”, International Journal of Mechanical Engineering and
Technology (IJMET), Volume 3, Issue 2, 2012, pp. 712-721, ISSN Print:
0976 – 6340, ISSN Online: 0976 – 6359.
[12] Dr. R. P. Sharma and Dr. R. V. Sharma, “A Numerical Study of Three-Dimensional
Darcybrinkman-Forchheimer (Dbf) Model in a Inclined Rectangular Porous Box”,
International Journal of Mechanical Engineering & Technology (IJMET), Volume 3,
Issue 2, 2012, pp. 702 - 711, ISSN Print: 0976 – 6340, ISSN Online: 0976 – 6359.

Mathematical modelling and analysis of three dimensional darcy

More Related Content

What's hot (17)

Viewers also liked (9)

Similar to Mathematical modelling and analysis of three dimensional darcy (20)

More from IAEME Publication (20)

Recently uploaded (20)

Mathematical modelling and analysis of three dimensional darcy