Applied Mathematics and Mechanics: (English Edition)
Applied Mathematics and Mechanics: (English Edition)
Applied Mathematics and Mechanics: (English Edition)
Applied Mathematics
and Mechanics
(English Edition)
Abstract
The combined eect of conduction-convection-radiation on natural convection ow of an optically thick Newtonian uid with gray radiant properties, conned in
a porous media square cavity with Darcy-Brinkman-Forchheimer drag is studied numerically. For a gray uid, Rosseland diusion approximation is considered. It is assumed that
(i) the temperature of the left vertical wall varies linearly with height, (ii) the right vertical
and top walls are at a lower temperature, and (iii) the bottom wall is uniformly-heated.
The governing equations are solved using the alternate direct implicit method together
with the successive over relaxation technique. The investigation of the eect of governing
parameters, namely, the Forschheimer resistance (), the temperature dierence (), and
the Plank number (Rd), on the ow pattern and heat transfer characteristics is carried
out. It can be seen that the reduction of ow and heat transfer occur as the Forschheimer
resistance is increased. On the other hand, both the ow strength and heat transfer increase as the temperature ratio is increased.
Key words conduction-convection-radiation, Darcy-Brinkman-Forchheimer drag, porous
medium, enclosure
Chinese Library Classification O361.3
2010 Mathematics Subject Classification
76W05
Nomenclature
aR ,
V,
Da,
x, y,
F,
X, Y ,
K,
k,
kr ,
g,
P r,
qw ,
qr ,
Ra,
Rd,
T,
TH ,
688
T0 ,
U, V ,
t,
cp ,
u, v,
,
,
,
,
,
,
,
dimensionless time;
stream function;
dimensionless stream function;
vorticity function;
dimensionless temperature;
porosity parameter;
dimensionless vorticity function.
Greek symbols
,
,
,
,
,
,
,
thermal diusivity;
coecient of thermal expansion;
Forschheimer resistance;
temperature ratio;
density of the uid;
kinematic viscosity;
Stefan-Boltzman constant;
Introduction
During the last three decades, a wide range of studies are conducted on conjugate heat
transfer phenomena in porous media. It is of practical importance in a number of applications,
for example, high performance insulation for buildings and cold storage installations. Nield
and Bejan[1] presented an extensive review of existing studies on such topics. Reviews can
also be found in Ingham and Pop[2] , Pop and Ingham[3] , Vafai[4] , and Al-Amiri[5] . Numerical
investigations of non-Darcian eects on transient conjugate natural convection-conduction heat
transfer from a two-dimensional vertical plate n embedded in a high-porosity medium were
carried out by Hung et al.[6] . The results showed that the inertial eects on heat transfer
characteristics are negligible at the beginning. However, these eects become increasingly
important over longer periods of time. Furthermore, the two-dimensional transient conjugate
free convection was investigated analytically and numerically by Vynnycky and Kimura[7] and
Kimura et al.[8] using vertical plate in a porous medium.
Sen[9] considered Darcy-Brinkman convective ow in a shallow porous rectangular cavity
with adiabatic upper and lower plate boundaries and dierentially heated sidewalls. An important non-Darcian study has also been presented by Lauriat and Prasad[10] for a heated
vertical porous cavity. The inertia and viscous forces on natural convection were examined
via the Darcy-Brinkman-Forchheimer model. A numerical investigation of transient free convection in a two-dimensional square cavity lled with a porous medium was conducted by
Saeid and Pop[11] . These results are in good agreement with the results obtained by Walker
and Homsy[12] , Bejan[13] , Weber[14] , Gross et al.[15] , and Manole and Lage[16] . Also, the results show that the time required to reach steady state is shorter for high Rayleigh number
and longer for low Rayleigh number. Finally, the classical Darcy formulation as well as the
Darcy-Brinkman model, Darcy-Forchheimer model, and Darcy-Brinkman-Forchheimer model,
has thoroughly been reviewed by Nield and Bejan[1] .
It is worth mentioned that heat transfer in uids which absorb and emit radiant energy is
a task complicated by the coupled, non1inear physical phenomena of internal radiation (see
Viskant and Anderson[17] ) and natural convection (see Bejan[13] ). There are several specic
problems of scientic and technological interest in which heat transfer via radiation and convection occurs in optically thick liquids, such as the ow of the earths mantle (see Matyska et
al.[18] ), the ow of oxide melts during crystal growth (see Rosenberger[19]), and the processing
of molten glass (see Gardon[20]). In general, the interaction of radiant heat transfer and natural convection in nite enclosures has received considerable attention. Lauriat[21] analyzed the
eect of combined radiation-convection in gray uids ow enclosed in vertical cavities. Studies
on natural convection of a participating uid in a rectangular cavity heated by incident radiant energy through a transparent side wall were performed by Webb and Viskanta[22] . Fusegi
and Farouk[23] studied the laminar and turbulent natural convection-radiation interactions in
a square enclosure lled with a non-gray gas for a wide range of Grashof number and Prandtl
number. Y
ucel et al.[24] and Tan and Howell[25] investigated the ow of an absorbing, emitting, and scattering uid in a dierentially heated square enclosure. Salinger et al.[26] studied
the eects of optical thickness and thermal gradients on the stability and structure of ows
in a cylindrical container heated from below. Debry et al.[27] examined the performance of
the approximations for modeling a representative problem of heat transfer and buoyant ow
in optically thick uids considering Rosseland diusion approximation (also known as P1).
This approximation is not applicable when the uid is optically thin medium (see Siegel and
Howell[28] ). Recently, Gerardin et al.[29] modied this diusion approximation and applied for
simulation of radiative transfer in a 3D geometry for an absorbing scattering medium.
Here, we propose to investigate the natural convection ow of an optically thick radiative
uid saturated porous media conned in a square cavity. To see the eect of radiation here
assume the Rosseland diusion approximation. In order to account for the viscous and inertial
eects in the momentum equation, the medium is considered to be porous due to DarcyBrinkman-Forchheimer drag. Finally, we assume that the temperature of left vertical wall
varies linearly with height, right vertical and top walls are cooled, and the bottom wall is
uniformly-heated. This has importance for applications in developing technology and industry,
such as prevention of sub-oil water pollution, storage of nuclear waste, and geothermal energy
systems (for which Cheng[30] has provided an extensive review).
Mathematical formulation
690
where F is the empirical constant in the second order resistance, is thermal expansion coefcient, is the density, is the dynamic viscosity, = ck p is thermal diusivity in which cp
is Molar specic heat at constant pressure, k is the coecient of thermal conductivity, and t
is the time. In equation (4), qr is the divergence of the radiant ux vector represents the
contribution of internal radiant transfer to the energy balance.
In equations (2) and (3), K is the measure of permeability of the porous medium (a packed
bed of spheres), dened by
K=
3
d2 ,
180(1 )2
where d is the diameter of the solid sphere, and is known as the porosity of the media and
dened by
=
Vf
.
Vc
Fig. 1
In the energy equation (4), qr is the diusion heat ux for radiative heat transfer, which has
the same form as that of the Fourier conduction law. Using the Rosseland mean absorption
coecient dened by Siegel and Howell[28] , radiative ux vector can be written as
qr = kr T,
(5)
16T 3
,
3aR
(6)
where and aR are the Stefan-Boltzmann constant and the mean absorption coecient, respectively. The energy equation together with equation (5) can be written as
cp
T
t
+u
T
T
+v
= k2 T + (kr T ).
x
y
(7)
Here, it is worth mentioned that Hossain and Wilson[31] investigated the problem of an
enclosure lled with a uid-saturated porous medium, along with generation of heat depending
on the uid temperature, but without eect of Forchheimer drag and internal radiation. On the
other hand, Derby et al.[27] studied the buoyant ow of an optically thick uid representative of
molten glass or an oxide crystal melt is solved the governing equations using the nite element
method. They considered four dierent radiation models, namely, the Rosseland diusion
approximation with and without the radiation of slip, the P1, the approximation, and a rigorous
numerical treatment, considering that the container top and bottom are cool and isothermal
and that the side-wall is heated with a parabolic temperature distribution.
The solution procedure discussed here is based on equations involving
u=
,
y
v=
,
x
(8)
v
u
.
x y
(9)
The vorticity equation obtained by eliminating the pressure between the momentum equations
is given as below:
+u
+v
t
x
y
F 2 (|V |v) (|V |u)
T
= 2
.
+ g
K
x
y
x
K
(10)
Now taking H as the reference length, and T0 the temperature of the undisturbed uid in the
cavity as the reference temperature, we consider the following dimensionless dependent and
independent variables
x
y
X = , Y = , = 2 t,
H
H
H
H
U = u, V = v, = ,
(11)
= H , = T T0 .
TH T0
Using the above dimensionless variables in vorticity equation (10) and the energy equation (7),
one may get the following dimensionless equations:
+U
+V
= 2 ( + |V |)
X
Y
|V |
|V | Ra
U
,
+
V
X
Y
P r X
1
4
+U
+V
=
1 + Rd(1 + )3 2
X
Y
Pr
3
2 2
4
.
Rd(1 + )2
+
+
Pr
X
Y
(12)
(13)
(14)
In a similar manner, from equation (8), we nd the non-dimensional form of the velocity components
U=
,
Y
V =
.
X
(15)
692
g(TH T0 )H 3
Ra
=
, Pr = ,
H2
1
F 2 H
=
=
, = ,
K
Da
K
Rd = 4T0 , = TH T0 .
kaR
T0
(16)
U = V = 0 for 0 Y 1 at X = 0,
U = V = = 0 for 0 X 1 at Y = 1,
U = V = 0, = 1 for 0 X 1 at Y = 0,
= 1 Y for 0 Y 1 at X = 0,
= 1 Y /A for 0 Y A at X = 1,
= 0 for A Y 1 at X = 1,
(17)
where
A=
L
1
=
.
H
10
Now, it is clear that for the present model, the non-dimensional parameters of interest
are the Rayleigh number Ra, the Prandtl number P r, the porosity parameter (1/Da), the
Forschheimer drag coecient , the Plank constant Rd, and the temperature ratio . In the
present investigation, pertaining to argon gas, the value of the Prandtl number is chosen as 0.7,
and the Rayleigh number is taken to be 2 105 .
Method of solution
For numerical simulation, here, successive over relaxation method with residual tolerance of
order 105 is applied on the stream function equation (14). Considering H to be the reference
height of the cavity, we take the uniform mesh size h = H/jmax , where jmax is the maximum
number of equi-spaced intervals along coordinate axes. The relaxation parameter r is obtained
from the relation (see [33]), i.e.,
1 (1 )
,
(18)
r=2
where
=
cos
jmax 1
+ cos
jmax 1
2
(19)
From these calculated values of stream function, the velocity components are updated using
equation (12) at each time step. For transient vorticity transport and energy equations (15)
and (16), we use alternate direct implicit (ADI) method. Given the values of variables and
in the ow eld at any time step, we compute the values of and at the next time step from
equations (15) and (16). Forward time central space discretization is used for the transient
and diusion terms in the ADI method, whereas for the convective terms, the ADI method is
modied using second upwind dierencing technique. However, like other implicit methods, the
stability requirement for vorticity at implicit wall boundaries imposes a restriction on the
time step, in a form similar to that of explicit schemes. For uniform mesh size h, is given
by (see [33])
1
2( h22 )
|u|+|v|
h
(20)
where is the parameter that stands for the coecient of diusion term of the transport
equations, whose solutions are required (see [33]). However, at solid wall boundaries, the time
step restriction reduces to
h2
.
4
(21)
Since the Prandtl number is xed to be 0.7, thus for grid size of 71 71, equation (21) requires
that a time step of = 106 would be appropriate for the entire computation in the proposed
method. The computations, to reach to the steady state, are carried out till
m+1
m
(i,j) (i,j)
< 108
m
(i,j)
is satised. Here, the superscript m refers to the number of time step and (i, j) are spatial
grid locations along X- and Y -axis, respectively. It is numerically testied that the steady state
regime for all ow variables lies well within this range.
As before, like in Hossain and Wilson[31] , a grid dependence study has been carried out for
the present thermally-driven cavity ow taking Rd =1.0, 2.0, and 4.0 while =0.5, Da = 1 000,
= 1.0, P r = 0.7, and Ra = 2 105 with meshes of 41 41, 55 55, 61 61, and 71 71
points. Thus, the numerical values obtained for min and max are entered in Table 1. Finally,
for computational economy, a 61 61 mesh is used throughout the present simulations.
Table 1
Rd
1.0
2.0
4.0
Numerical values of min and max for Rd = 1.0, 2.0, and 4.0 while = 0.5, Da = 1 000,
= 1.0, P r = 0.7, and Ra = 2 105
41 41
55 55
61 61
71 71
min
max
min
max
min
max
min
max
15.382
14.203
12.378
1.303 5
0.969 1
0.284 4
15.431
14.195
12.362
1.356 7
0.961 2
0.264 1
15.430
14.196
12.358
1.347 6
0.960 2
0.262 5
15.402
14.188
12.345
1.336 6
0.951 9
0.250 8
694
drag ( = 0) was investigated by Hossain and Wilson[31] for P r = 0.7 and = 10. However, knowing the temperature distribution, we can evaluate the net radiative heat ux in the
medium. However, the limitation to the accuracy of the optically thick limit approximation
should be recognized. According to Ozisik[34] , this approximation breaks down in the intermediate vicinity of the boundaries because it does not take into account radiation from the boundary
surface. This is a serious restriction, since heat transfer at the wall is important. Thus, with
black body condition for optically thick limit the net wall heat ux forms the following relations:
2
3
q
Rd(1
+
|
=
1
+
)
,
L
X=0
3
X X=0
qR = 1 + Rd(1 + |X=1 )3
,
3
X X=1
(22)
q
Rd(1
+
|
=
1
+
)
,
Y =0
B
3
Y Y =0
q = 1 + 2 Rd(1 + |
3
.
T
Y =1 )
3
Y Y =1
Here, qL , qR , qB , and qT are, respectively, the net heat ux from the left, right, bottom, and top
surface of the cavity. It follows that the net average heat transfer corresponding to the above
net heat transfer coecients are calculated as
1
1
qL dY, QR =
qR dY,
QL =
0
0
(23)
1
1
qB dX, QT =
qT dX,
QB =
0
where QL , QR , QB , and QT are, respectively, the net average heat ux from the left, right,
bottom, and top surface of the cavity.
Fig. 2
Fig. 3
Streamlines for dierent while = 103 , Rd=0.0, = 0.5, P r = 0.7, and Ra = 2 105
Isotherms for dierent while = 103 , Rd=0.0, = 0.5, P r = 0.7, and Ra = 2 105
696
The labeled contours of isotherms also depict how the temperature is distributed along the
vertical wall. It can be seen that the pattern of temperature distribution is not much aected
due to the increasing eect of the Forchheimer resistance .
Figures 4(a)4(d) show the heat transfer rate of the solid walls at Ra = 2 105 , P r = 0.7,
Rd = 0.0, and = 1 000 for dierent values of . Figure 4(a) shows the heat transfer rate of
the left wall. It can now be observed that for = 0, a wave like structure of heat transfer curve
appears. However, the distribution transforms to a single peak curve for = 50. Moreover,
this peak value decreases with the increase in the value of . Thus, heat transfer along the
left wall numerically decreases with the increase in the value of . Figure 4(b) shows the heat
ux distribution for the right wall. It can be seen that for all values of , the heat transfer
curve shoots down from its maximum value in the region 0 Y A, which is the region of
steep temperature variation along the right wall. However, the curves drop to the minimum
close to Y = 0.2. The decrease in the local maximum of the curves in the region 0.2 Y 1
show that the heat transfer rate decreases for increasing values of . Figure 4(c) shows the
heat transfer curves for the bottom wall. Comparing Figs. 4(a) and 4(c), we see that the heat
ux distribution is dierent for left and bottom walls which may well be attributed to the
imposed surface heating conditions. As the bottom wall is at a constant temperature and the
heat transfer shoots to maximum close to X = 1, which may be due to the reason that X = 1
is the region of steep temperature dierence along the right wall. This behaviour is due to the
temperature conditions imposed along these walls. Figure 4(d) now shows the heat transfer
rate of the top wall. For all values of Forschheimer resistance , the heat ux is unity along the
left and zero along the right end of the top wall due to the boundary conditions. In between, a
peak value appears in each curve. This corresponds to the region where the uid coming from
Fig. 4
Rate of heat transfer qw from dierent walls for =0.0, 25.0, and 50.0 while = 103 , Rd=0.0,
= 0.5, P r = 0.7, and Ra = 2 105
the left wall meets the top surface. However, like the other three walls, the decrease in this
peak value shows that the heat transfer decreases with increasing values of . Thus, from the
above discussion, we conclude that the heat transfer rate for all the four walls decreases with
the increase in the value of Forchheimer resistance, and it is due to decreased uid motion as
an eect of .
4.2 Eect of surface temperature
Figure 5 shows the pattern of streamlines at Ra = 2 105 , P r = 0.7, and = Rd = 1.0
for (a) = 0.5, (b) = 1.0, and (c) = 1.5, respectively. It can now be seen that with the
increase in the value of , the volume ow rate and the size of the primary cell increases. This is
now expected since the increase in the value of corresponds to an increase in the temperature
distribution close to the walls at a higher temperature and/or alternatively a decrease in the
temperature distribution of the uid close to the walls at a lower temperature distribution.
This makes sense why the strength of the secondary cell decreases for increasing values of .
Further, not only the size, but also the strength of ow in the primary cells increase from 0.56
Fig. 5
Streamlines for dierent while = 103 , =1, Rd=1.0, P r = 0.7, and Ra = 2 105
to 0.90 as is increased from 0.5 to 1.5. Figure 6 shows isotherms at Ra = 2105, P r = 0.7, and
= Rd = 1.0 for (a) = 0.5, (b) = 1.0, and (c) = 1.5, respectively. One can see that the
the temperature distribution remains smoother, and does not vary signicantly as increases.
This may be due to the reason that the ow pattern also does not change very signicantly.
Since corresponds to the temperature dierence, we discern that an increase in would
increase the temperature dierence between the walls. Figures 7(a)7(d) now show the heat
ux distribution of the left, right, bottom and top walls, respectively, at Ra = 2105 , P r = 0.7,
= Rd = 1.0, = 1 000 for dierent values of . Figure 7(a) shows the heat ux distribution
of the left wall. From this gure, we now discern that the wave like structure that appeared in
Fig. 4(a) for = Rd = 0, does not appear for non-zero values of and Rd. However, one peak
698
Fig. 6
Fig. 7
Isotherms for dierent while = 103 , =1, Rd=1.0, P r = 0.7, and Ra = 2 105
Rate of heat transfer from dierent walls for =0.5, 1.0, and 1.5 while Rd = 1.0, = 103 ,
=1.0. P r = 0.7, and Ra = 2 105
value appears which numerically increases for increasing values of . Figure 7(b) shows the
heat ux distribution of the right wall. It can now again be seen that the heat transfer drops
to its minimum close to Y = 0.2. However, heat transfer slightly increases as increases. This
is now expected since an increase in causes the temperature gradient along the lower part
of the right wall also to increase. Nevertheless, unlike Fig. 4(b), this peak value in Fig. 7(b)
asymptotically drops to zero in the range 0.2 Y 1, showing that the heat transfer increases
for increasing values of . Figure 7(c) depicts the heat ux distribution of the bottom wall. It is
now easy to understand that the heat ux distributions of the left wall reaches to its maximum
close to X = 1 due to imposed boundary temperature around the right bottom. This maximum
value also increases as increases. Figure 7(d) nally represents the heat ux of the top wall,
which shows that like other three walls, the heat transfer rate of the right wall also increases
as increases. However, comparing it with Fig. 4(d), we see that its value at X = 0 is not
unity, rather it is 1.7 which is due to non-zero value of Rd, that makes contribution in heat
transfer given by equations (22). Also, the heat ux in this case is greater than the one in the
case of Rd = 0 given in Fig. 4(d).
4.3 Eect of Plank number Rd
Figures 8 and 9 respectively show the streamlines and isotherms at Ra = 2 105, P r = 0.7,
= 1.0, = 103 , = 0.5 for (a) Rd = 0.0, (b) Rd = 2.0, and (c) Rd = 4.0. Figures 8(a) 8(c)
reveal that there is an overall decrease in the strength of ow both in primary and secondary
cells, with the increase in the value of Rd. Moreover, the transition of ow from two cell to
one cell pattern also suggests that the contribution of Rd causes a reduction in the convective
ow. The increase in the radiation eects may enhance the uid temperature, which ultimately
reduces the temperature dierence between the walls and the ambient uid, thereby reducing
the overall ow rate. Figures 9(a) 9(c) indicate that the temperature distribution in the ow
region is not much aected by the variation of Rd.
Fig. 8
Streamlines for dierent Rd while =0.5, = 103 , =1.0, P r = 0.7, and Ra = 2 105
700
Fig. 9
Isotherms for dierent Rd while = 0.5, = 103 , = 1.0, P r = 0.7, and Ra = 2 105
The consideration of black body conditions in optically thick limits of solid walls leads
to a heat transfer characteristics, which is also complemented with the eect of Rd. This
net/enhanced heat transfer rate with the added eect of Rd is dened in equations (22). Figure
10 represents the heat ux distributions of the (a) left, (b) right, (c) bottom, and (d) top solid
walls at Ra = 2 105 , P r = 0.7, = 103 , = 1.0, = 0.5, for Rd = 0.0, 2.0, and 4.0,
respectively. It is now clear from these gures that the heat transfer rate of the four walls
increases with the increase in Rd. However, Fig. 10(d) needs some attention. The heat ux
distribution at the left end top is unity for Rd = 0, but it increases for non zero values of Rd.
This is due to the enhanced eect of Rd on the heat ux.
Conclusions
Fig. 10
Rate of heat transfer from dierent walls for Rd =0.0, 2.0, and 4.0 while = 0.5, = 103 ,
=1.0, P r = 0.7, and Ra = 2 105
References
[1] Nield, D. A. and Bejan, A. Convection in Porous Media, 3rd ed., Springer Science + Business
Media, New York (2006)
[2] Ingham, D. B. and Pop, I. Transport Phenomenon in Porous Media, Pergamon, Oxford (1998)
[3] Pop, I. and Ingham, D. B. Convective Heat Transfer : Mathematical and Computational Modeling
of Viscous Fluids and Porous Media, Pergamon, Oxford (2001)
[4] Vafai, K. Handbook of Porous Media, 2nd ed., Taylor and Francis Group, New York (2005)
[5] Al-Amiri, A. Natural convection in porous enclosures: the application of the two-energy equation
model. Numerical Heat Transfer : Part A, 41, 817834 (2002)
[6] Hung, C. I., Chen, C. O. K., and Cheng, P. Transient conjugate natural convection heat transfer
along a vertical plate n in a high-porosity medium. Numerical Heat Transfer : Part A, 15, 133
148 (1989)
[7] Vynnycky, M. and Kimura, S. Transient conjugate free convection due to a vertical plate in a
porous medium. Int. J. Heat Mass Tran., 38, 219231 (1995)
[8] Kimura, S., Kiwata, T., and Okajima, A. Conjugate natural convection in porous media. Advanced
in Water Resources, 20, 111126 (1997)
[9] Sen, A. K. Natural convection in a shallow porous cavitythe Brinkman model. Int. J. Heat Mass
Tran., 30, 855868 (1987)
[10] Lauriat, G. and Prasad, V. Non-Darcian eects on natural convection in a vertical porous enclosure. Int. J. Heat Mass Tran., 32, 21352148 (1989)
[11] Saeid, N. H. and Pop, I. Transient free convection in a square cavity lled with a porous medium.
Int. J. Heat Mass Tran., 47, 19171924 (2004)
702
[12] Walker, K. L. and Homsy, G. M. Convection in a porous cavity. J. Fluid Mech., 87, 474499
(1978)
[13] Bejan, A. Convection Heat Transfer, Wiley-Interscience, New York (1984)
[14] Weber, J. E. The boundary-layer regime for convection in a vertical porous layer. Int. J. Heat
Mass Tran., 18, 569573 (1975)
[15] Gross, R. J., Baer, M. R., and Hickox, C. E. Application of ux-corrected transport (FCT) to
high Rayleigh number natural convection in a porous medium. Proceedings of the International
Heat Transfer Conference on Heat Transfer, Lalifornia, San Francissco (1986)
[16] Manole, D. M. and Lage, J. L. Numerical benchmark results for natural convection in a porous
medium cavity. American Society of Mechanical Engineers, 216, 5560 (1992)
[17] Viskanta, R. and Anderson, E. E. Heat transfer in semitransparent solids. Advances in Heat
Transfer, Academic Press, New York (1975)
[18] Matyska, C., Moser, J., and Yuen, D. A. The potential inuence of radiative ransfer on the
formation of megaplumes in the lower mantle. Earth and Planetary Science Letters, 125, 255266
(1994)
[19] Rosenberger, F. E. Fundamentals of Crystal Growth I: Macroscopic Equilibrium and Transport
Concepts, Springer-Verlag, Berlin (1979)
[20] Gardon, R. A review of radiant heat transfer in glass. J. Am. Ceram. Soc., 44, 305312 (1961)
[21] Lauriat, G. Combined radiation-convection in gray uids enclosed n vertical cavities. Journal of
Heat Transfer, 104, 609615 (1982)
[22] Webb, B. W. and Viskanta, R. Analysis of radiation-induced natural convection in rectangluar
enclosures. Journal of Thermophysics, 1, 146153 (1987)
[23] Fusegi, T. and Farouk, B. Laminar and turbulent natural convection diation interactions in a
square enclosure lled with a nongray gas. Numerical Heat Transfer: Part A, 15, 303322 (1989)
[24] Y
ucel, A., Acharya, S., and Williams, M. L. Natural convection and radiation in square enclosures.
Numerical Heat Transfer: Part A, 15, 261278 (1989)
[25] Tan, Z. and Howell, J. R. Combined radiation and natural convection in a two-dimensional participating square medium. Int. J. Heat Mass Tran., 34, 785793 (1991)
[26] Salinger, A. G., Brandon, S., Aris, R., and Derby, J. J. Buoyancy-driven ows of a radiatively
participating uid in a vertical cylinder heated from below. Proc. R. Soc. Lond. A, 442, 313341
(1993)
[27] Derby, J. J., Brandon, S., and Salinger, A. G. The diusion and P1 approximations for modeling
buoyant ow of an optically thick uid. Int. J. Heat Mass Tran., 41, 14051415 (1998)
[28] Siegel, R. and Howell, J. R. Thermal Radiation Heat Transfer, 2nd ed., McGraw-Hill, New York
(1981)
[29] Gerardin, J., Seiler, N., Ruyer, P., Trovalet, L., and Boulet, P. P1 approximation, MDA and IDA
for the simulation of radiative transfer in a 3D geometry for an absorbing scattering medium.
Journal of Quantitative Spectroscopy and Radiative Transfer, 113, 140149 (2012)
[30] Cheng, P. Heat transfer in geothermal systems. Advanced in Heat Transfer, 4, 1105 (1978)
[31] Hossain, M. A. and Wilson, M. Natural convection ow in a uid-saturated porous medium
enclosed by non-isothermal walls with heat generation. Int. J. Therm. Sci., 41, 447454 (2002)
[32] Saha, S. C. and Gu, Y. T. Free convection in a triangular enclosure with uid-saturated porous
medium and internal heat generation. ANZIAM Journal, 53, C127C141 (2011)
[33] Roache, P. J. Computational Fluid Dynamics, 2nd ed., Hermosa, Albuquerque, New Mexico (1998)
[34] Ozisik, M. N. Radiative Transfer and Interactions with Conduction and Convection, John-Wiley
& Sons, New York (1973)