American Society of Heating, Refrigerating and Air-Conditioning Engineers

Flow Pattern and Thermal Comfort in Office Environment with Active Chilled Beams Journal: American Society of Heating, Refrigerating and Air-Conditioning Engineers, Inc. Manuscript ID: Draft Publication: HVAC&R Research Keywords: Commercial Buildings, Air Distribution Equipment, Research & Development

Page of 0 American Society of Heating, Refrigerating and Air-Conditioning Engineers 0 0 Flow Pattern and Thermal Comfort in Office Environment with Active Chilled Beams Hannu Koskela Henna Häggblom Risto Kosonen, PhD Mika Ruponen, PhD ABSTRACT In modern offices, the heat load per floor area has increased. With high cooling loads, the possibility of draft problems increases. The purpose of this paper was to study the flow patterns and draft risk in office environment where cooling and air distribution is implemented with active chilled beams. The study is based on experiments in a laboratory mock-up room in three load conditions: summer, winter and midseason (spring/autumn). Thermal plumes from heat sources and warm or cold windows had a notable effect on the flow pattern and velocity distribution in the occupied zone. Areas with increased draft risk were found in locations where the supply jet turns down to the occupied zone. Draft risk can also be high at the floor level as a result of a circulating flow pattern in the room. This paper concentrates on measurement and modeling results in a single-person office room. Comparisons are made with corresponding results in open-plan office. INTRODUCTION In modern offices, the heat load per floor area has increased due to higher density of workstations, increased heat load from equipment and high solar load from large unshaded windows. With high cooling loads, the air distribution becomes more difficult and the possibility of draft problems increases. The flow pattern in the room becomes more unstable due to interactions between cool inlet jets and warm convection flows from heat sources and large eddies appear into the room spaces (Müller et al. 00). The outdoor condition can have a notable effect on the room air flows. Cool window and outer wall surfaces can create downward flows and heaters or warm surfaces upward flows that affect the room flow pattern.

Page of 0 0 0 The purpose of this study was to examine the flow patterns and draft risk in two typical office spaces: a single-person office room and an open-plan office. Cooling and air distribution were implemented with active chilled beams. The study was based on experiments in a laboratory mock-up room in three different outdoor load conditions: summer, winter and midseason (spring/autumn) in Scandinavian type of climate. This paper concentrates on measurement and modeling results in a single-person office room. Comparisons are made with corresponding results in open-plan office. The results of the open-plan office case have been reported in more detail by Koskela et al. (00). METHODS The experiments were carried out in the environmental chamber at the Finnish Institute of Occupational Health in Turku, Finland. The two test room layouts are shown in Figure. The test room was built to represent a typical building module of a modern flexible office building. The width of the module is defined by the distance between the construction beams of the building; typically. m. One wall of the test room had six windows of size. m x. m. Their surface temperature was controlled by blowing air into the chamber behind them. Convective heaters were placed under the windows and were only used in the winter test case. The division of the windows and installation of exposed chilled beams enables walls to be built in several locations. The single-person office room was separated with a wall from the right end of the module. Figure. Dimensions and layout of the test rooms: open-plan office (left) and single-person office room (right).

Page of 0 American Society of Heating, Refrigerating and Air-Conditioning Engineers 0 0 The single-person room had one workstation with a computer and display located at the centre of a side wall. A simple plate was installed as a flow obstacle on the opposite wall to represent a typical location of a cupboard or a shelf. The room had one light fitting above the workstation at the height of. m. The open plan office layout had eight workstations arranged symmetrically in two groups. The internal heat load to the room was produced with computers, displays, light fittings and dummies representing persons. The purpose was to have approximately the same heat load per floor area in singleperson room and open-plan office test cases. Therefore, somewhat higher load values were used in the single-person room case. The persons were simulated with painted cylinders (height. m, diameter 0. m) similar to those defined in DIN () containing light bulbs of total 0 W in the open-plan office and light bulbs of total 0 W in the single-person room. Each workstation had a PC and display with power consumption adjusted to 0 W in the open-plan office and to 0 W in the single-person room. The lamps had a power consumption of 0 W each. Direct solar load was simulated by placing heater panels (size.0 m x 0. m, height 0.0 m) on the floor as shown in Figure. In the single-person room, the panels close to the windows were used in the summer conditions. The other panels further off in the room were used in the spring/autumn conditions, when the elevation angle of the sun is lower. The heaters under the windows in the winter test case had dimensions of. m x 0. m x 0. m. They produced mainly convective heat load due to stainless steel covering on all vertical surfaces. The exposed chilled beams were installed asymmetrically in the test rooms in the centerline of every second window. The distance from the ceiling was 0. m. The dimensions of the beams were. m x 0. m, height 0. m. The active length of the beams was.0 m and it was located symmetrically in relation to the centerline of the room. The active chilled beam model Halton CCE was selected to represent a typical unit with exposed installation in the room. In the device, outdoor air supply is combined with cooling of re-circulated room air. Outdoor air supply is typically introduced through small nozzles inside of the beam. Outdoor air jets induce room air through a heat exchanger, where it is cooled and the mixture (referred as inlet jets) is blown into the room through supply air openings. The flow rate of induced room air is typically - times the outdoor air flow rate. The heat exchanger consists of a cooling coil, which is cooled by water flow. In

Page of 0 0 0 this case, the inlet jets are blown from inlet jet openings in the upper surface of the beam. The supply of inlet jets creates plane jets to both sides of the beam that normally attach to the ceiling utilizing the Coanda effect. Figure shows the operation principle of an active chilled beam. Figure. Operating principle of the active chilled beam. The measurements of air velocity were carried out by using ultrasonic anemometers (Kaijo Denki WA- 0, accuracy ± 0.0 m/s), moved by an automated traversing system. The averaging time in each measurement point was s. The measurement grid density was 0. m x 0. m. Additional measurements were done using Dantec N0 flow analyzer with hot sphere sensors. The flow pattern was visualized using smoke and video recorded during experiments. The cooling power of the chilled beams was determined by measuring the cooling water flow rate and the rise of water temperature in the heat exchanger. In the single-person office room, distributions of air velocity and temperature were measured in horizontal and vertical planes shown in Figure. The horizontal planes were at the heights of. m and. m. The planes were. m wide and extended from the outer wall to the inner wall. Measurement height. m was used with the automated traversing system instead of the standard height. m because of obstacles. In the open-plan office case, only horizontal planes at heights of 0. m,. m and. m were measured. The locations of the planes were selected to cover the four workstations in the central part of the room.

Page of 0 American Society of Heating, Refrigerating and Air-Conditioning Engineers 0 0 Figure. Vertical and horizontal measurement planes in the single-person office room. The CFD simulations were carried out by using Ansys CFX software with SST turbulence model following the guidelines given by Nielsen et al. (00). The grid was unstructured with inflation layers on the room and heat source surfaces and had 000 nodes in the single-person office room case and 0 000 nodes in the open plan office case. The grid density was cm in the supply opening, cm in the supply jet area and 0 cm in other parts of the room. No radiation model was used. The convective parts of the heat loads were given to the surfaces as heat fluxes. The radiation part was distributed to room surfaces based on their approximate view factors. The mean air speed was calculated from mean air velocity and turbulent kinetic energy using a correction formula reported by Koskela et al. (00). ISO standard gives design criteria for maximum mean air speed in the office environment by defining three categories for different draft rate (DR) levels (ISO 00): Category A (DR 0 %): summer 0. m/s, winter 0.0 m/s Category B (DR 0 %): summer 0. m/s, winter 0. m/s Category C (DR %): summer 0. m/s, winter 0. m/s These values are based on the assumption that the room temperature is in the lower end of the corresponding recommended temperature range and the turbulence intensity is %.

Page of 0 0 0 RESULTS Three test cases were measured in the laboratory and simulated with CFD:. Summer case with -00 W/m cooling load (warm windows and direct solar heat load). Spring/autumn case with W/m cooling load (cold windows). Winter case with W/m (cold windows and heaters under them) In the summer case the heat load was W/m in the single-person office room case and 00 W/m in the open plan office case. The heat load levels in the test cases are presented in Table. Heat load type Persons (cylinders) Computers and displays Lights Other internal loads Heaters Solar load (panels) Windows (calculated) Total load Total load per floor area Table. Heat load levels in the test cases. Single-person office room Summer (W) 0 0 0 0 0 0 Spring/ autumn (W) 0 0 0 00 00-0 0 Open-plan office Winter (W) Summer (W) Spring/ autumn (W) 0 0 0 0 0 0 0 0-0 0 0 0 00-00 0 Winter (W) In the winter case, the heater below the window warmed up the window surface, which lowered the calculated heat loss from the room to the window surface. In the single-person office room, some extra load was added with the floor panels in the spring/autumn case in order to achieve the heat load level W/m. CFD-modeling The simplified model for the chilled beam was made of similar size as the actual device. The outlets representing room air re-circulation openings were larger than in reality covering the whole vertical surfaces of the beam. The inlet boundary conditions were determined based on the information obtained from the manufacturer. The total inlet air flow rate was l/s per unit. The dimensions of the inlet slots were.0 m x 0.0 m. This gives an inlet momentum flow rate of 0. N per unit. The air flow was blown with velocity 0. m/s to an angle of compared to the vertical direction. The inlet jet temperature was calculated from the mean room temperature and the cooling power of the chilled beam. The modeling methods and results of the open plan office case are reported in more detail by Koskela et al. (00). Single-person office room 0 0-0

Page of 0 American Society of Heating, Refrigerating and Air-Conditioning Engineers 0 0 The flow patterns and air speed distributions in the single-person room are shown in Figs. -. In the summer and winter cases, the flow direction close to the windows is upwards. This upward flow creates a circulation, which turns the inlet jet towards the corridor wall. The velocity maximum at. m level is therefore close to the corridor wall and does not cause draft to the workstation. In the spring/autumn case, the cool window surface causes a downward flow, which creates on opposite circulation compared to the summer and winter cases. Supply jet turns towards the window boosting the downward flow of cool air to the floor level.

Page of 0 0 0 Figure. Measurement results in the summer test case with W/m cooling load. Air flow pattern from smoke experiments (top), air speed distributions in the vertical measurement planes (middle) and air speed distributions at. m level (bottom left) and. m level (bottom right). Figure. Measurement results in the spring/autumn test case with W/m cooling load. Air flow pattern from smoke experiments (top), air speed distributions in the vertical measurement planes (middle) and air speed distributions at. m level (bottom left) and. m level (bottom right).

Page of 0 American Society of Heating, Refrigerating and Air-Conditioning Engineers 0 0 Figure. Measurement results in the winter test case with W/m cooling load. Air flow pattern from smoke experiments (top), air speed distributions in the vertical measurement planes (middle) and air speed distributions at. m level (bottom left) and. m level (bottom right). Effect of workstation location on the flow pattern The effect of workstation location on the flow pattern was studied with a pair of experiments in summer conditions with a heat load of W/m. The workstation was moved from the normal location in centre of the side wall to a new location close to the window. The results are shown in Figure. When the workstation is in the centre of the side wall, the downfall of the supply jet occurs closer to the corridor wall.

Page 0 of 0 0 0 This can be explained by the effect of the thermal plumes rising from the workstation heat sources to the supply jet. Also, the maximum air velocities at. m level were somewhat smaller with the workstation in the central position. However, the overall effect on the room air flow pattern was not large. Figure. Effect of the workstation location on the air speed distributions at. m level (middle) and. m level (bottom). Effect of a flow obstacle on the flow pattern The effect of the horizontal plate on the wall representing a cupboard or a shelf was also studied with a pair of experiments in the summer conditions with a heat load of W/m. In the standard configuration, the plate turned one part the supply jet from the chilled beam towards the centre of the room instead of continuing downwards along the side wall (Figure ). This caused an interaction of the two inlet jets thus widening the flow profile of the downward air flow. The maximum velocities close to the workstation were somewhat smaller with the plate installed.

Page of 0 American Society of Heating, Refrigerating and Air-Conditioning Engineers 0 0 Figure. Effect of the flow obstacle on the wall on the air speed distributions at. m level (middle) and. m level (bottom). CFD-modeling results The flow patterns obtained from the CFD-simulations in the single-person room test cases are shown in Figs. - with a comparison of air speed distributions in the measurement plane at the. m level. The mean and maximum values of air speed in the occupied zone from CFD-simulations and measurements are compared in Figure. CFD-simulations were able to correctly predict the main flow pattern in the singleperson room in all three test cases. The predicted mean air speed values were also close to the measurement results in all cases. The maximum velocities, however, were notably higher in all cases.

Page of 0 0 0 Figure. Results of the CFD-simulation in the summer test case with W/m cooling load. Air flow pattern visualized with the 0. m/s iso-surface of air speed (top) and air speed distributions at. m level from measurements (bottom left) and CFD-simulation (bottom right).

Page of 0 American Society of Heating, Refrigerating and Air-Conditioning Engineers 0 0 Figure 0. Results of the CFD-simulation in the winter test case with W/m cooling load. Air flow pattern visualized with the 0. m/s iso-surface of air speed (top) and air speed distributions at. m level from measurements (bottom left) and CFD-simulation (bottom right). Figure. Results of the CFD-simulation in the spring/autumn test case with W/m cooling load. Air flow pattern visualized with the 0. m/s iso-surface of air speed (top) and air speed distributions at. m level from measurements (bottom left) and CFD-simulation (bottom right).

Page of 0 0 0 Open-plan office Figure. Comparison of mean and maximum values of air speed in the occupied zone between CFD-simulation and measurement results The main features of the open-plan office results are shown here, but they are presented in more detail by Koskela et al. (00). The main flow patterns in the summer test case detected from the smoke experiments in the open-plan office are shown in Figure (on the left) together with the flow patterns from the CFD simulation (on the right). Figure shows the modeled distribution of air speed at two horizontal planes at 0. m and. m heights with corresponding measurement results in the central part of the office module. A strong longitudinal circulation was formed in the room due to asymmetric layout of the chilled beams compared to heat sources.

Page of 0 American Society of Heating, Refrigerating and Air-Conditioning Engineers 0 0 Figure. Flow pattern in the summer case, experimental and CFD-results Figure. Measured (top) and measured (bottom) air speed contours at 0. m (left) and. m (right) levels in the summer case. In the autumn/spring test case, the longitudinal circulation was weaker compared to the summer case (Figure on the left). This is natural, because of weaker buoyancy forces. Also the downward plumes from the cool window were weak. In the winter case, the longitudinal circulation was also weaker compared to the summer case (Figure on the right). The heaters under the windows created upward plumes as in the summer case. These plumes turned the inlet jets somewhat towards the corridor wall.

Page of 0 0 0 Figure. Flow patterns in the spring/autumn case W/m (left) and winter case W/m (right) from smoke experiments. Air speed levels in the test cases A comparison between the air speed levels in the occupied zone of the open-plan office and the singleperson room is presented in Figure. The mean and maximum air speed results are compared with the target values of ISO for summer and winter conditions. The standard does not give target values of air speed for spring or autumn conditions. Categories A, B and C correspond to draft risk levels of 0 %, 0 % and %. At. m height, the mean air speed values are at the same level (0.0 0. m/s) in both room types. Also the maximum air speed is at the same level in summer (0. 0. m/s) and spring/autumn (0.0 0. m/s) cases. In the winter case, the maximum value in the single-person room is higher (0. m/s) than in the open-plan office (0.0 m/s). At 0. m height, the velocities in the open-plan office are higher in summer (0. m/s) and winter (0. m/s) conditions than the corresponding values in the single-person room (0. m/s and 0. m/s). In spring/autumn conditions, however, the velocities are higher in the single-person room (0. m/s) than in the open-plan office (0. m/s). The main reason for the floor level air flows was the occurrence of large scale circulation in the room. The cause of this circulation was different in open-plan office and single-person room cases. In the openplan office, the main factor was the asymmetric layout of cooling units compared to heat sources and room geometry. This type of layout is typical in modern flexible offices. In the single-person room, the

Page of 0 American Society of Heating, Refrigerating and Air-Conditioning Engineers 0 0 circulation was mainly caused by the heat sources close to the window or the downdraft from cold window surface. Figure. Comparison of mean and maximum values of air speed between the open-plan office and the single-person office room compared to ISO design criteria. SUMMARY AND DISCUSSION The heat sources had a notable influence on the flow pattern causing large scale circulation and affecting the direction of inlet jets. Two main causes of draft risk were found: ) Downfall of inlet jets causing local maxima of air speed especially at the head level and ) Large scale circulation causing high air speeds especially at the floor level. In the single-person room, the direction of the circulation in the room depended on the time of the year. In the summer case it was caused by the upward plume from the warm window and in the winter case by

Page of 0 0 0 the heater below the window. In the spring/autumn case the circulation turned the supply jet towards the cool window. The combination of cool supply jet and downward flow from the widow caused relatively high velocities on the floor level. In the summer and winter cases the highest velocities were found close to the corridor wall and away from the workstation.. Findings are similar to those reported by Melikov et al. (00) and Zboril et al. (00). In the open plan office, the main circulation occurred in the longitudinal direction of the room. It was caused by the asymmetry of the heat loads and the cooling units. Also a transverse circulation was created by the heat load of the windows or the heaters below the windows. It had the effect of turning the inlet jets towards the inner wall. The transverse circulation, however, was overridden by the stronger longitudinal circulation. The collision and downfall of the inlet jets was another phenomenon causing high velocities in the open plan office. The downfall occurs locally and the position of draft risk areas can change in time and also due to changes in room heat sources. The maximum air speed in the occupied zone was in most cases relatively high compared to the recommendations, in category C or above. The mean air speed was typically in category B on the head level and in category C on the floor level. No overall difference between the velocity levels in the two room types was found. The air speed in the open plan office was highest in the summer case with high cooling load. In the single-person room, the floor level velocities were high also in the spring/autumn case. It has to be noted, however, that the ISO target values for air speed in Figure assume that the temperature is in the lower end of the corresponding recommended temperature range, which may overestimate the draft risk. The office rooms had only few flow obstacles and the space under the tables was mainly open, which made the large scale circulation flow along the floor possible. In real offices there are usually more flow obstacles, which prevent the large scale circulation and reduce the air speed. If the screens in the open plan office block the flow totally under the tables, this type of circulation is not possible. The location of the workstation had some effect on the air speed distribution. The plumes of the workstation heat sources seemed to change the location and air speed of the downfall of the inlet jet. The

Page of 0 American Society of Heating, Refrigerating and Air-Conditioning Engineers 0 0 horizontal plate on the wall representing a cupboard or a shelf also had an effect on the flow pattern and air speed distribution in the occupied zone. The CFD-simulations were able to predict the main features of the flow pattern in the room. The predicted mean air speed values were close to the measurement results in all cases. The maximum air speed values, however, were notably higher than the measured values in all cases. This is a typical result and is caused by the inability of the steady-state RANS-models to correctly predict the fluctuating room air flows. Following general conclusions can be drawn based on this study: ) Convection air flows have a notable effect on the thermal conditions in the room. The effects depend on the prevailing load conditions. Analysis of these effects requires full scale experiments or simulation with CFD-modeling techniques; analytical flow element models are not sufficient. ) Midseason conditions must be analyzed also; extreme load conditions in summer or winter do not give the whole picture. ) Current draft standards do not fully describe transient and asymmetric flow conditions such as in this study. New methods of assessing the thermal conditions in rooms should be developed. ACKNOWLEDGEMENTS The financial support of National Technology Agency of Finland is greatly appreciated. REFERENCES DIN.. DIN -; Chilled surfaces for rooms. Part. DIN, Germany. ISO. 00. ISO Standard, Ergonomics of the thermal environment - Analytical determination and interpretation of thermal comfort using calculation of the PMV and PPD indices and local thermal comfort criteria. International Organisation for Standardisation, Genève. p. Koskela, H., H. Häggblom, R., Kosonen, and M. Ruponen. 00. Air distribution in office environment with asymmetric workstation layout using chilled beams. Building and Environment (00) -. Koskela, H., J. Heikkinen, R. Niemelä, and T. Hautalampi. 00. Turbulence correction for thermal comfort calculation. Building and Environment 00; (): -. Melikov, A., B. Yordanova, L. Bozkhov, V. Zboril, and R. Kosonen. 00. Human response to thermal environment in rooms with chilled beams. Proceedings of Clima 00 Wellbeing Indoors. Finnvac ry, ISBN ----, Finland. Müller. D., I. Gores, and R. Zielinski. 00. Impact of the Thermal Load on the Room Airflow Pattern. In: Roomvent 00. Proceedings of th International Conference on Air Distribution in Rooms; 00 Sept -; Coimbra, Portugal Nielsen, P. V., F. Allard, H. B. Awbi, L. Davidson, and A. Schälin. 00. Computational fluid dynamics in ventilation design. REHVA Guidebook no 0. REHVA.

Page 0 of 0 0 0 Zboril, V., A. Melikov, B. Yordanova, L. Bozkhov, and R. Kosonen. 00. Airflow distribution in rooms with active chilled beams. Proceedings of Roomvent 00 0th international conference on air distribution in rooms. Finvac ry. ISBN ----, Finland.