Available online at www.sciencedirect.com
ScienceDirect Procedia Engineering 121 (2015) 1193 – 1199
9th International Symposium on Heating, Ventilation and Air Conditioning (ISHVAC) and the 3rd International Conference on Building Energy and Environment (COBEE)
Effect of surface thermal radiation on natural convection and heat transfer in a cavity containing a horizontal porous layer Yuancheng Wanga, *, Jun Yanga, Xiaojing Zhanga ,Yu Pana a
College of Thermal Energy Engineering, Shandong Jianzhu University, Jinan 250101, China
Abstract In this paper, the effects of surface radiation on heat transfer and natural convection in a cavity containing a horizontal porous layer have been studied numerically. The governing equations for the momentum and heat transfer in both free fluid and porous medium were solved by the finite element method. The radiative heat transfer is calculated by making use of the radiosity of the surfaces that assumed to be grey. Comparisons with experimental and numerical results in the literature have been carried out. Effects of thermal radiation on natural convection and heat transfer in both free fluid and porous medium were analyzed. It was found that surface thermal radiation can significantly change the temperature fields in both the regions of free flow and porous medium. The mean temperature at the interface decreases and the temperature gradients are created on the upper two corners of the porous medium region as Ra increases. © 2015 2015The TheAuthors. Authors. Published by Elsevier Ltd.is an open access article under the CC BY-NC-ND license © Published by Elsevier Ltd. This Peer-review under responsibility of the organizing committee of ISHVACCOBEE 2015. (http://creativecommons.org/licenses/by-nc-nd/4.0/). Peer-review under responsibility of the organizing committee of ISHVAC-COBEE 2015
Keywords: natural convection; heat transfer; surface thermal radiation; cavity; porous medium
1. Introduction Natural convection and heat transfer in composite porous/fluid domain exist in many natural phenomena and engineering applications, such as grain storage, air conditioning systems and insulation used in buildings. During the past several decades, several experiments and numerical simulations have been presented to describe the phenomena of natural convection and heat transfer with and without thermal radiation in a cavity [1-8]. The published literature also contains several studies of natural convection in cavities that are partially filled with porous medium, however, none appears to include the interaction of thermal radiation and natural convection [9-12]. Thermal radiation E-mail address:
[email protected]
1877-7058 © 2015 The Authors. Published by Elsevier Ltd. This is an open access article under the CC BY-NC-ND license
(http://creativecommons.org/licenses/by-nc-nd/4.0/). Peer-review under responsibility of the organizing committee of ISHVAC-COBEE 2015
doi:10.1016/j.proeng.2015.09.137
1194
Yuancheng Wang et al. / Procedia Engineering 121 (2015) 1193 – 1199
between the walls and roof of the buildings and grain storage silos affects the natural convection flow and heat transfer, and influence the moisture transportation, so it is important to investigate the interaction of surface thermal radiation with natural convection in a cavity. In this paper, we focus on natural convection and heat transfer with thermal radiation in a cavity partially filled with a porous medium in which a differentially heated cavity in which the two side walls are held at constant temperatures, and the upper and lower walls are deemed to be adiabatic, and a classical investigation of the thermal physical parameters that affect the natural convection and temperature fields is conducted. 2. Mathematical formulation The system studied in this work consists of a differentially heated cavity of width W and height H, and the porous medium is half of the whole cavity. 2.1 Governing equations In this analysis we shall consider the fluid to be air, which are incompressible and laminar. The thermophysical properties of the fluid are assumed constant, except for the density in the buoyancy term in the momentum equations. The porous medium is considered as bulk wheat, which are homogeneous, isotropic and in local thermodynamic equilibrium with air. The governing conservation equation for the free fluid and the porous medium will be written separately. For the fluid region, we have
Continuity: v
Momentum:
0
(1)
U vv P 2 v p U o gE T To
Thermal energy:
U o ca v T
(2)
k 2T
(3)
in which v is the velocity vector, p is the pressure, T is the temperature of the fluid and porous medium, the density of the fluid,
Uo
is the density of the fluid at the temperature
U is
To , g is the gravity vector, E is the
coefficient of volumetric expansion of the fluid, c, k and P are the specific heat, thermal conductivity and viscosity of the fluid, respectively. The third term on the right hand side of the equation (2) arises as a result of Boussinesq’s approximation. For the porous region, the governing equations are
Continuity: Momentum:
vD
0
U v D v D P 2 v D p U o gE T To PK 1 v D
(4) (5)
1195
Yuancheng Wang et al. / Procedia Engineering 121 (2015) 1193 – 1199
Thermal energy:
U p c p v D T
k eff 2T
(6)
in which v D is the volume averaged Darcian or superficial velocity of the fluid through the porous medium, K is the permeability tensor, which in this work is taken to be isotropic, U p is the density of porous medium, c p is specific heat of porous medium, keff is the effective thermal conductivity of the porous medium. The boundary conditions for the equations (1) to (6) are
v
n y T
0, W and y
0
0 ; Q kn y T
Q k eff n T
Where
0, H ; T
0 when x
vD
when
porous
y
k n T
fluid
when
0
Y
Th when x 0 ; T
when
y
Yint erface T ;
H ;v porous
Tc when x W ;
v D when Y
T
fluid
when
Yint erface Y
;
Yint erface
.
n y is unit normal vector, the argument of which is positive in the direction leaving the surface. Q is the
net radiation flux on the surface. There is an exchange of radiant energy between the two side walls and the roof of the cavity and the interface of the saturated porous medium and fluid. The fluid above the porous region is deemed to be transparent to thermal radiation and not participating in radiative heat transfer. The porous medium is opaque to thermal radiation at its upper surface which can absorb and emit radiation. 2.2. Validation of the model The model was validated by comparing simulation results with two sets of published data. One set, that of Beckermann et al.[10] deals with heat transfer in a cavity partially filled with a porous medium, but that excludes the effects of thermal radiation. Beckermann et al.[10] presented both experimental and numerical studies on natural convection in a rectangular cavity filled to half of its width with a vertical layer of saturated porous medium, namely glass beads. The aspect ratio, H/W, of the cavity is 1.0, H is the surface emissivity, Ra is the Rayleigh number, Ra
ca U 2 gE (Th Tc ) H 3 /( kP ) , Pr is Prandtl number, Da is the Darcy number, Da
K / H 2 and
K is the permeability of the porous medium. Rk is the ratio of thermal conductivities between the porous medium k eff / k , k is the thermal conductivity of fluid. The surface emissivity, Rayleigh, Darcy and and fluid and Rk Prandtl numbers are summarized in Table 1. In this paper, the results of experiment 2 and experiment 4 of the experiments were compared with our numerical results.
Table 1. Summary of experimental conditions from Beckermann et al. [10] Test number Experiment 2 Experiment 4
Test fluid water glycerin
Ra 3.208×107 3.471×106
Pr 6.97 12630
Da 1.296×10-5 7.112×10-7
Rk 1.383 2.259
H 0.38 0.38
1196
Yuancheng Wang et al. / Procedia Engineering 121 (2015) 1193 – 1199
A comparison between measured and predicted dimensionless temperatures by the present numerical model is shown in Figs. 1. It can be seen that agreement between the two sets of values is good, especially in the Fig.1a. The discrepancies occur in Fig.1b, which is the comparison between measured and predicted results for glycerin. This is possibly due to the inaccuracies in determining the exact position of the movable thermocouple probe, and nonuniformities in the porosity at the walls, which is expected to produce a considerable difference between the numerical models. In addition, the viscosity of glycerin varies by almost an order of magnitude over the temperature range in the experiments [10]. 1.0
0.9
0.9
0.8
0.8
0.7
0.7
0.6
0.6
0.5
0.5
T
T
1.0
0.4
0.4 0.3 0.2 0.1 0.0
Analysis,y/L=0.847 Analysis,y/L=0.496 Analysis,y/L=0.055 Experiment,y/L=0.847 Experiment,y/L=0.496 Experiment,y/L=0.055
0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0
X
(a) experiment 2
0.3 0.2 0.1
Analysis,y/L=0.875 Analysis,y/L=0.500 Analysis,y/L=0.141 Experiment,y/L=0.875 Experiment,y/L=0.500 Experiment,y/L=0.141
0.0 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 &
(b) experiment 4
Fig.1 The dimensionless temperature comparison between the experimental and predicted results for experiment 2 and experiment 4
3. Results and discussion A cavity with an air fluid overlying a porous medium, which occupies half of the whole cavity, was chosen as prototype system for investigation. The system comprises the classical differentially heated cavity, in which the isothermal vertical side-walls are maintained at constant but different temperatures, the upper and lower surfaces of the cavity are adiabatic. The three cases are investigated with the following parameters: Tc = 288.5K, Th = 298.5K, T0=(Th+Tc)/2, Rk = 5.31, Da= 5.78e-7 and H =0.0, 0.1,0.5,1.0. In order to demonstrate the effect of surface radiation on the flow fields, the isotherm plots from the simulations for cases with and without radiation at several emissivities are presented in figure 2. Comparing with results without radiation, the results with radiation show that temperature rises at the interface and decreases at the top wall in the free flow region. Because of the natural convection and thermal radiation at the interface between the free flow and porous medium, the mean temperature in porous medium is higher than that without thermal radiation, and the mean temperature in the porous medium increases slightly with the increase of emissivity on the surfaces. The prediction in figure 2 also indicates that the appearance of surface thermal radiation can significantly shifts the line of the mean temperature (T=293.5K) towards the right in the porous medium region and change the temperature fields even if the given emissivity is quite small. However, the emissivity from 0.5 to 1.0, the changes are not as significant as that as emissivity increases from 0.0 to 0.5.
Yuancheng Wang et al. / Procedia Engineering 121 (2015) 1193 – 1199
(a) H =0.0
(c)
H =0.5
(b)
H =0.1
(d)
H =1.0
1197
Fig. 2. the isotherms for the nonradiation and radiation with several emissivities (Ra =106ˈRk=5.31, Da= 5.78e-7)
The effects of emissivity on mean convective Nusselt number are shown in Figs 3(a)-(c) for Ra=106,105and104, respectively. The general effect of surface emissivity on the mean convective Nusselt number, Nu c , on the both hot and cold side wall in the cavity is different between the free flow region and porous medium region. The
Nu c on
the side wall in the porous medium region is much lower than that in the free flow region, and changes very little because of the weaker natural convection and lack of thermal radiation in the porous medium. The thermal radiation on the surfaces makes the temperature gradient near the hot side wall increase slightly in the free flow region, and thus rises the Nu c on the hot side slightly with the increase of emissivity. However, as radiative energy is proportional to the fourth power of absolute temperature, the radiative heat exchange between the hot wall and the insulated (top) wall becomes stronger than that between the cold and insulated walls, the temperature gradients are
1198
Yuancheng Wang et al. / Procedia Engineering 121 (2015) 1193 – 1199
weakened near the hot wall but strengthened near the cold wall. Thus
Nu c at the hot wall is lower than that at the
2.0 1.9 1.8 1.7 1.6 1.5 1.4 1.3 1.2 1.1 1.0 0.9 0.8 0.7 0.6 0.5 0.4
hot side wall in the free flow region hot side wall in the porous medium region cold side wall in the free flow cold side wall in the porous medium region
Nuc
Nuc
cold wall, and the temperature gradient near the cold side increase slightly in the free flow region with further increase of emissivity. Effect on radiative Nusselt number due to variation of emissivity is shown in Fig. 3(d). As expected, Nu r increases rapidly with the increase of emissivity. This indicates that surface radiation plays an important part in heat transfer in the cavity, especially at higher emissivity.
0.0
0.2
0.4
H
0.6
0.8
4.4 4.2 4.0 3.8 3.6 3.4 3.2 3.0 2.8 2.6 2.4 2.2 2.0 1.8 1.6 1.4 1.2 1.0 0.8 0.6 0.0
1.0
(a) Ra =104ˈRk=5.31, Da= 5.78e-7
20
4.0
18
3.5
16
H
0.6
0.8
1.0
14 hot side wall in the free flow region hot side wall in the porous medium region cold side wall in the free flow cold side wall in the porous medium region
2.5 2.0
12
Nur
3.0
Nuc
0.4
hot side wall in free flow region,Ra=1e6 hot side wall in free flow region,Ra=1e5 hot side wall in free flow region,Ra=1e4 cold side wall in free flow region,Ra=1e6 cold side wall in free flow region,Ra=1e5 cold side wall in free flow region,Ra=1e4
22
4.5
0.2
(b) Ra =105ˈRk=5.31, Da= 5.78e-7 24
5.0
hot side wall in the free flow region hot side wall in the porous medium region cold side wall in the free flow cold side wall in the porous medium region
10 8
1.5
6
1.0
4 2
0.5
0 0.0
0.2
0.4
H
0.6
0.8
(c) Ra =106ˈRk=5.31, Da= 5.78e-7
1.0
0.0
0.2
0.4
H
0.6
0.8
1.0
(d) Ra=104, 105,106, Rk=5.31, Da= 5.78e-7
Fig. 3 The effect of surface emissivity on the mean convective and radiative Nusselt numbers with several Rayleigh numbers
4. Conclusions The combined radiation-natural convection heat transfer in a differentially heated cavity partially filled with a porous medium has been investigated numerically. Comparisons with experimental and numerical simulation results from the literature have been carried out to check the accuracy of the present numerical method. The investigation shows the effects of thermal radiation on natural convection and heat transfer in both the free fluid and porous medium.
Yuancheng Wang et al. / Procedia Engineering 121 (2015) 1193 – 1199
1199
The surface thermal radiation alters the distribution of the temperature in the free flow and porous medium regions even if the given emissivity is quite small, but there are no remarkable changes with the increase of emissivity from 0.5 to 1.0. Comparing with nonradiation, the temperature rises at the interface and decreases at the top insulated wall, resulting in the temperature in porous medium being higher than that without thermal radiation, and the temperature in the porous medium increases slightly with the increase of emissivity on the surfaces. Acknowledgements This study is supported by the National Natural Science Foundation of China (No. 51276102), Shandong Provincial Natural Science Foundation of China (ZR2011EEM011). The authors also wish to express their appreciation to Shandong Provincial Key Laboratory of Building Energy-Saving Technique and Key Laboratory of Renewable Energy Utilization Technologies in Buildings of the National Education Ministry for providing the support for this study. References [1] G. de Vahl Davis, Natural convection of air in a square cavity: a benchmark numerical solution, Int. J. Numer. Meth. Fl. 3(1983) 249-264. [2] N.C. Markatos, K.A. Pericleous, Laminar and turbulent natural convection in an enclosed cavity, Int. J. Heat. Mass. Tran. 27 (1984) 755-772. [3] R.A. W.M. Henkes, C.J. Hoogendoorn, Scaling of the laminar natural convection flow in a heated square cavity, Int. J. Heat. Mass. Tran. 36 (1992) 2913-2923. [4] A.D. Orazio, M. Coreione, G.P. Celata, Application to natural convection enclosed flows of a Lattice Boltzmann BGK model coupled with a general purpose thermal boundary condition, Int. J. Therm. Sci. 43 (2004) 575-586. [5] C. Balaji, S.P. Venkateshan, Interaction of surface radiation with free convection in a square cavity. Int. J. Heat. Fluid. Fl. 14 (1993) 260-267. [6] N. Ramesh, P. Venkateshan, Effect of surface radiation on natural convection in a square enclosure, Int. J. Therm. Sci. 13 (1999) 299-301. [7] E. H. Ridouane, M. Hasaoui, A. Amahmid, A. Raji, Interaction between natural convection and radiation in a square cavity heated from below, Numerical Heat Transfer, Part A. 45 (2004) 289-311. [8] H. Bouali, A. Mezrhab, H. Amaoui, M. Bouzidi, Radiation-natural convection heat transfer in an inclined rectangular enclosure, Int. J. Therm. Sci. 45 (2006) 553-566. [9] T. Nishimura, T. Takumi, M. Shiraisi, Y. Kawamura, H. Ozoe, Numerical analysis of natural convection in a rectangular enclosure horizontally divided into fluid and porous regions, Int. J. Heat. Mass. Tran. 29 (1986) 889-898. [10] C. Beckermann, S. Ramadhyani, R. Viskkanta, Natural convection flow and heat transfer between a fluid layer and a porous layer inside a rectangular enclosure, Journal of Heat Transfer. 109 (1987) 363-370. [11] J. A. Ochoa-Tapia, S. Whitaker, Momentum transfer at the boundary between a porous medium and a homogeneous fluid – I. Theoretical development, Int. J. Heat. Mass. Tran. 38 (1995) 2635-2646. [12] J. A. Ochoa-Tapia, S. Whitaker, Momentum transfer at the boundary between a porous medium and a homogeneous fluid – II. Comparison with experiment, Int. J. Heat. Mass. Tran. 38 (1995) 2647-2655.