Journal of Applied Mathematics

Volume 2012 (2012), Article ID 287271, 12 pages

http://dx.doi.org/10.1155/2012/287271

## Airflow and Heat Transfer in the Slot-Vented Room with Radiant Floor Heating Unit

^{1}School of Civil Engineering, Hunan University, Hunan, Changsha 410082, China^{2}Department of Building Engineering, Hunan Institute of Engineering, Hunan, Xiangtan 411104, China^{3}Soochow University, Jiangsu, Suzhou 215137, China^{4}Institute of Built Environment and Control, Zhongkai University of Agriculture and Engineering, Guangdong, Guangzhou 510225, China

Received 13 January 2012; Accepted 20 February 2012

Academic Editor: M. F. El-Amin

Copyright © 2012 Xiang-Long Liu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

Radiant floor heating has received increasing attention due to its diverse advantages, especially the energy saving as compared to the conventional dwelling heating system. This paper presents a numerical investigation of airflow and heat transfer in the slot-vented room with the radiant floor heating unit. Combination of fluid convection and thermal radiation has been implemented through the thermal boundary conditions. Spatial distributions of indoor air temperature and velocity, as well as the heat transfer rates along the radiant floor and the outer wall, have been presented and analyzed covering the domains from complete natural convection to forced convection dominated flows. The numerical results demonstrate that the levels of average temperature in the room with lateral slot-ventilation are higher than those without slot-ventilation, but lower than those in the room with ceiling slot-ventilation. Overall, the slot-ventilation room with radiant floor heating unit could offer better indoor air quality through increasing the indoor air temperature and fresh air exchanging rate simultaneously. Concerning the airborne pollutant transports and moisture condensations, the performance of radiant floor heating unit will be further optimized in our future researches.

#### 1. Introduction

Radiant floor heating is a novel indoor heating technique that has received increasing attention in recent years. In the conventional heating methods, such as the air convection by radiators and the hot air supply by air conditioners, the hot air is always located in the upside or middle upside of the room instead of the occupied zone, leading to the obvious temperature differences between the upside and underside of the room. Under such condition, occupants usually feel uncomfortable. In the radiant floor heating, however, the effective temperature of the floor is higher than the temperature in the upside of the room, which can supply more heat to the underside of the room for warming the occupants’ feet instead of their heads that follows the rule of human’s physiologic adjustment. Therefore, the radiant floor heating improves the thermal comfort condition for the occupants as compared to the conventional heating methods. In addition, the radiant system has the advantages of energy saving, long life, wide selection of the heat sources, low cost, and environmental friendliness. The theoretical and experimental research [1, 2] showed that radiant floor system can decrease more than 30% of the energy consumption. The radiant floor system is also used for indoor cooling in summer only by changing the heat source to the cool source. In recent years, the indoor radiant heating or cooling system has been receiving increasing utilization in residential and commercial buildings. For example, among the newly built European buildings, many of them have been equipped with the radiant heating/cooling systems.

Radiant heating system is such a complex system that involves the various heat transfer mechanisms, including the heat conduction in the floor, the radiative heat transfer between radiant surface and other surfaces, convective heat transfer between the radiant surface and its neighboring air, and the buoyancy’s effects. The numerical method provides a convenient way to solve such problem. Ghaly and Elbarbary [3] and Shateyi [4] numerically studied the convection coupled with radiative heat transfer and buoyancy’s effects. Ma et al. [5] numerically studied the characteristics of the temperature and velocity distribution in the room, and their results showed that radiant heat flux accounted for 50–60% of the total heat flux. Sattari and Farhanieh [6] studied the effects of design parameters on performance of a typical radiant floor heating system using finite element method. Bozkir and Canbazoǧlu [7] carried out the experimental and numerical research on the radiant heating system with an attempt to find the characteristics. In practice, the adequate fresh air supply is usually required to ensure indoor air quality. As a result, the floor radiant heating system is usually combined with the slot-ventilation system. However, few of the research works focuses on the air flow and heat transfer in the slot-vented room with radiant floor heating unit.

This paper presents the numerical simulation of floor radiant heating system with three types of slot-ventilation, that is, lateral slot-ventilation (LSV), ceiling slot-ventilation (CSV) and no slot-ventilation (NSV). The indoor temperature, velocity, and heat flux distributions are calculated. This work is of great significance for the design and wider application of the floor radiant heating system.

#### 2. Mathematical Model

The steady 2D model of the turbulent airflow and heat transfer in the room is developed in this study. The schematic of radiant heating system and coordinate system is shown in Figure 1. The radiant floor is assumed as the constant temperature ; the outside wall is cooled at the constant heat flux ; other walls are adiabatic. The air velocity at the inlet is , and the inflow temperature is . The standard two equations are used to describe the turbulent flow and the governing equations are given as follows.

Continuity equation: Momentum equations: where , are the fluid velocities in - and -directions respectively; is the pressure of the fluid; represents the acceleration of gravity. is effective viscosity coefficient, which is equal to the sum of molecular viscosity and turbulent eddy viscosity . Turbulent eddy viscosity is determined by the local turbulent kinetic energy and the dissipation rate : Energy equation in terms of temperature is given below Transport equations for turbulent kinetic energy and dissipation rate equations are in (2.5) is the production term given by The closure coefficients , and in the above turbulent model are set as constants 0.09, 1.44, 1.92, 1.0, 1.0, and 1.3, respectively [8].

The DO radiation model is employed to handle the radiative heat transfer, and the absorbing and scattering of air are omitted: where is the distance between the two points of different wall surfaces and is the absorption coefficient. Each of the wall surfaces is assumed as grey-body radiation surfaces, and the radiation emissivity at the floor is , and are taken as 0.8 at the other surfaces.

#### 3. Numerical Method

Based on the control volume and SIMPLEC method [9], the governing equations (2.1)–(2.7) subject to the boundary conditions are solved in the iterative manner. The grid is , with the depressed mesh at near the wall to match the high velocity and temperature gradients. In the calculation, the Boussinesq assumption is used to consider the buoyancy’s effect, and the thermal expansion coefficient is 0.0033. The thermophysical properties of air are constant except the density. The converged solution is obtained when the following convergence criteria are satisfied for the dependent variables:

#### 4. Results and Discussion

The calculation conditions are as follows: the constant wall temperature = 30°C; the inlet velocity = 0.3, 0.6 and 0.9 m/s; the cooling heat flux of the outside wall = −100 W/m^{2}. The air temperature at the inlet is equal to the mass-averaged temperature without slot-ventilation so that no additional heat load is induced into the room. The indoor air temperature and velocity distributions, as well as the radiative and convective heat flux between radiant floor and outside wall, are calculated under the different slot-ventilation modes (i.e., LSV, CSV, and NSV).

##### 4.1. Velocity and Temperature Distributions

Figure 2 shows the temperature and velocity distributions for the three slot ventilation modes at = 30°C, = −100 W/m^{2}, and = 0.6 m/s. The air temperature in the room is quite uniform (the temperature difference is less than 1°C), except in the area near the floor. The air temperature is lower in the area near outside wall. As no fresh air is supplied for NSV (Figure 2(a)), the air near outside wall is firstly cooled down and then moves downward due to the density difference in the vertical direction. The cooled air flows along the floor due to the stagnation effect as it touches the floor and is warmed up at the floor. Then the heated air moves upwards along the inside wall due to the buoyancy’ effect. As a result, the clockwise natural convection of indoor air is formed. It is clearly seen that only in the region near the wall the obvious air flow is observed, and the velocity is quite low elsewhere. As shown in Figure 2(b), the indoor air temperature is higher under LSV mode as compared to the case with NSV. Because both the air inlet and air outlet are installed at the same inside wall (one is near the ceiling and the other is near the floor), the attaching jet flow is formed along the ceiling as the air is supplied. The air is then cooled near the outside wall and goes downwards so that the flow pattern similar to the natural convection is formed, but the air velocity is higher for LSV. Due to the higher air velocity along the walls, the convective heat transfer between the air and the wall is enhanced for LSV. Under CSV, (Figure 2(c)), the indoor air temperature is further increased, and the air flow pattern is quite different from those under above-mentioned two modes. As the air is supplied from the middle of ceiling, the air from the inlet is continuously expanded and decelerated along the flow direction. The supplied air is divided into two parts; one flows towards the outlet and is exhausted, and the other flows towards to the outside wall which is cooled and goes down, forming the same flow pattern as the natural convection. As the magnitude of the air velocity is close to that of natural convection velocity, the flow pattern is the combined result of the forced convection and the natural convection.

Figure 3 shows the temperature distribution at the different locations for the cases shown in Figure 2. In this figure, , are dimensionless distance , . It is found that the average temperature for CSV is highest, while the average temperature for NSV is lowest. For all the cases, the temperature gradients along the height direction are quite small except in the area near the floor and ceiling, which indicates the uniform temperature distribution in most of the area. The temperature gradients near the floor are quite large, especially in the area closer to outside wall (e.g., ). This is because the air cooled at the outside wall produces the greater temperature difference as it is passing along the floor. Near the ceiling, because of the influences of the radiation, convection, cooling from the outside wall, heat transfer mechanism is complex. As a result, the obvious temperature gradient is also observed near the ceiling.

Figure 4 shows the velocity distribution at different locations for the cases shown in Figure 2. It is clearly found that the influence of the slot-ventilation on the flow pattern is obvious. Under these three conditions, air velocity in the occupied zone is all less than 0.2 m/s, while the air velocity of the area near the ceiling and floor is as high as 0.5 m/s. Compared to NSV, the velocity near the floor and the ceiling for LSV and CSV is much higher, and CSV leads to the more complex flow pattern.

Figure 5 shows the temperature and velocity distributions for LSV at = 30°C, = −100 W/m^{2}, and = 0.3, and 0.9 m/s. It is clear that the airflow pattern and temperature distributions are similar to the case shown in Figure 2(b). As the inlet velocity increases ( = 0.9 m/s), (Figure 5(a)), more air flows towards the outside wall, leading to the more intensive convective heat transfer between the indoor air and the outside wall. As a result, the air is cooled to the lower temperature. For this reason, the bulk temperature of the indoor air is lower. In reverse, as the inlet velocity is lower ( = 0.3 m/s) (Figure 5(b)), less air reaches the outside wall, and the bulk temperature of the indoor air is higher because of the weaker convective heat transfer along the outside wall.

Figure 6 shows the temperature and velocity distributions for CSV at = 30°C, = −100 W/m^{2}, and = 0.3 and 0.9 m/s. The airflow pattern and temperature distributions are similar to the case shown in Figure 2(c). As the inlet velocity increases ( = 0.9 m/s) (Figure 6(a)), more air flows are supplied, leading to the more intensive convective heat transfer between the air and the outside wall and the floor. However, because the direction of the inlet velocity is towards the floor, the convective and radiative heat transfer along the floor are dominant. As a result, the bulk temperature of the air in the room is higher. In reverse, as the inlet velocity is lower ( = 0.3 m/s) (Figure 6(b)), the convective heat transfer along both the outside wall and the floor decreases, which decreases the bulk air temperature in the room as compared to the case of = 0.9 m/s. However, the air temperature at = 0.3 m/s is higher than that for = 0.6 m/s. Therefore, the air temperature is related to the intensities of the convective heat transfer along the outside wall and the floor. There exists an optimal inlet velocity for CSV to increase the bulk temperature of the indoor air.

##### 4.2. Characteristics of Radiative and Convective Heat Transfer

The different behaviors of the indoor air temperature under the different conditions are the combined results of the different heat transfer mechanisms. The heat fluxes along the floor and outside wall are studied and the Nusselt numbers for the convective heat transfer *Nu _{c}* and the radiative heat transfer

*Nu*along the wall surface are given. Radiative heat flux and convective heat flux are defined as where is net radiative heat flux between the cells at the outside wall or floor and at other walls; is the distance of two points; and are the angles to the surface normal, and is heat conductivity. The radiative and convective heat transfer coefficients are defined as and : where is the bulk temperature of indoor air,

_{r}*Nu*can be calculated through , in which is characteristic length, which is taken as 3.6 m for the floor and 2.9 m for the outside wall.

Figures 7 and 8, respectively, present the distributions of *Nu _{c}* and

*Nu*distribution along with floor and outside wall for the cases shown in Figure 2. It is found that

_{r}*Nu*is greater than

_{r }*Nu*in most of the areas under the different conditions, which indicates that the radiative heat transfer is stronger than the convective heat transfer. For NSV, the radiative heat flux accounts for about 56% of the total heat flux. For LSV and CSV, the percentages of the radiative heat transfer are about 54% and 61%, respectively. In the area near the outside wall, the intense convection is observed due to the local vortex and higher temperature gradients (Figure 3) and thus

_{c}*Nu*is greater

_{c}*Nu*, indicating the stronger convective heat transfer in this area. The magnitude of the radiative heat flux increases with the increased ; therefore, the radiation is more significant in the area closer to the low-temperature outside wall. Under LSV and NSV, the air velocity and hence

_{r}*Nu*decrease with the decreased . The velocity under CSV increases firstly and then decreases along , and

_{c}*Nu*has the same tendency. Near the inside wall, the local vortex is formed that increases the convective heat transfer and thus

_{c}*Nu*increases to some extent. There is no significant difference between

_{c}*Nu*under the three conditions. For LSV,

_{r}*Nu*is of the largest magnitude as the result of the strongest convective heat transfer. As ,

_{c }*Nu*for CSV is greater than that for NSV. It is also found from Figure 8 that

_{c}*Nu*is greater than

_{r}*Nu*for all the three conditions except in the area near air inlet for the CSV (), which suggests that the radiation dominates the heat transfer mechanism. Under LSV,

_{c}*Nu*is larger as as compared to the two other cases, indicating the strongest convection heat transfer in this area. As , for CSV, the air velocity significantly increases near the inlet, which increases the convection heat transfer.

_{c}#### 5. Conclusions

The numerical stimulations of the airflow and heat transfer in the slot-vented room with radiant floor heating unit, including LSV, CSV, and NSV, are performed. The temperature, velocity, and the behaviors of the radiative and convective heat transfer are calculated. The results show that the air temperatures in the room are quite uniform under these three air-supply modes and the average temperature difference is less than 1°C. Under LSV, the forced ventilation has the same moving direction with the natural convection and thus increases the air velocity. Under CSV, the flow pattern is more complex as the combined result of the forced ventilation and natural convection. Radiation is main heat transfer mechanism at the floor and the outside wall under all the three conditions. Radiative heat flux accounts for 50–60% of the total heat flux, which is largest under CSV. LSV offers the strongest convection along the floor and outside wall, which increases the convective heat transfer, especially at cold outside wall. Therefore, the average temperature for LSV is lower than that for CSV.

Concerning the airborne pollutant transports and moisture condensations, the performance of radiant floor heating unit will be further optimized in our future researches.

#### Acknowledgments

The authors gratefully acknowledge the funding for this project provided by International Science and Technology Cooperative Project of China (2010DFB63830), Major Science and Technology Project of Hunan Province (2010FJ1013), and National Natural Science Foundation of China (51178477).

#### References

- I. B. Kilkis, S. S. Sager, and M. Uludag, “A simplified model for radiant heating and cooling panels,”
*Simulation Practice and Theory*, vol. 2, no. 2, pp. 61–76, 1994. View at Publisher · View at Google Scholar · View at Scopus - A. K. Athienitis and T. Y. Chen, “Experimental and theoretical investigation of floor heating with thermal storage,”
*ASHRAE Transactions*, vol. 99, no. 1, pp. 1049–1057, 1993. View at Google Scholar - A. Y. Ghaly and E. M. E. Elbarbary, “Radiation effect on MHD free-convection flow of a gas at a stretching surface with a uniform free stream,”
*Journal of Applied Mathematics*, vol. 2, no. 2, pp. 93–103, 2002. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - S. Shateyi, “Thermal radiation and buoyancy effects on heat and mass transfer over a semi-infinite stretching surface with suction and blowing,”
*Journal of Applied Mathematics*, vol. 2008, Article ID 414830, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - L. D. Ma, W. Q. Tao, Y. Dai, and D. X. Sun, “Numerical simulation of temperature distribution and turbulent flow in the radiant floor heating room,”
*Journal of Engineering Thermophysics*, vol. 26, no. 3, pp. 501–503, 2005 (Chinese). View at Google Scholar · View at Scopus - S. Sattari and B. Farhanieh, “A parametric study on radiant floor heating system performance,”
*Renewable Energy*, vol. 31, no. 10, pp. 1617–1626, 2006. View at Publisher · View at Google Scholar · View at Scopus - O. G. Bozkir and S. Canbazoǧlu, “Unsteady thermal performance analysis of a room with serial and parallel duct radiant floor heating system using hot airflow,”
*Energy and Buildings*, vol. 36, no. 6, pp. 579–586, 2004. View at Publisher · View at Google Scholar · View at Scopus - N. C. Markatos, “The mathematical modelling of turbulent flows,”
*Applied Mathematical Modelling*, vol. 10, no. 3, pp. 190–220, 1986. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - S. V. Partankar,
*Numerical Heat Transfer and Fluid Flow*, McGraw–Hill, 1980.