A comparative study on cooling period thermal comfort assessment in modern open office landscape in Estonia

.


Introduction
Modern low energy office buildings require energy efficient heating, ventilation, and air conditioning (HVAC) systems which can provide comfortable and healthy indoor environment.In temperate climate countries, mechanical ventilation and active cooling systems are common practice in such buildings.However, mechanical HVAC systems do not always provide satisfactory thermal conditions [1].It is important to properly apply control strategies, design and install room cooling units and ventilation supply air elements, as well as to operate and maintain the systems to provide comfortable indoor climate without temperature fluctuations and draught risk in the cooling season [2][3][4][5][6][7][8].Office plans, in terms of occupant positions and density, can be very different from initial design and vary significantly, resulting in changing conditions and dynamic settings which makes it difficult to design the systems adequately to ensure stable thermal environment.Open office layout design is used commonly in most office buildings mainly to allow flexibility in workspaces allocation [9].This creates a difficult task for HVAC systems design, requiring careful planning to assure adequate conditions in the occupied zone in different layout cases.
As occupant satisfaction with thermal environment is dependent on many factors, such as gender, age, health, activity, mood, and other physiological and psychological factors, assessing thermal comfort (TC) based on temperature and air movement measurements is usually not sufficient for adequate estimation [6,[10][11][12][13].Thus, evaluation by questioning the occupants is usually also needed to specify the problems and get a comprehensive overview of the TC situation.Studies on office workers thermal sensation have shown that the predicted TC and actual sensation can differ significantly [12,14,15].For example, gender specific analysis indicates higher dissatisfaction rates for female occupants [11,[16][17][18][19]. Recent research has widely focused on individual perception of TC [20][21][22], developed methods to analyze the preferences for TC using machine learning algorithms [23,24] and adapt systems to provide preferable personal comfort by implementing Personalized Comfort Systems [20][21][22]25].Utilization of such systems in buildings requires paradigm shifts in occupant interaction with HVAC systems as well as system design practices, integration of advanced controls and information technologies solutions [26,27].
In addition to the individual preferences and system specific aspects influencing thermal comfort (TC), there are many building related design factors that can affect the performance of HVAC systems and in turn influence the thermal environment.Of these factors, façade design, namely window sizes, layout, and glazing parameters, can have large impact on cooling load as well as radiant temperature asymmetry and thus major influence on the overall thermal conditions in the office [28,29].Thalfeldt et al. [28] showed the importance of façade design by analyzing the effect on office buildings energy efficiency and cooling load in cold climate countries.Window-to-wall ratio (WWR) of 0.25 was found optimal for triple glazing window solutions.Larger glazing results in higher cooling loads and increase the need for larger room cooling units, higher cooled airflow rates or lower supply air temperatures to maintain the room temperature.The latter factors also increase the risk of draught in occupied spaces.In several studies, draught rate (DR) has been identified as the main cause of discomfort even if other thermal environment factors are at satisfactory levels [6,15,29,30].
Depending mainly on the cooling load, cooling plant solution, and interior design, different water based room cooling solutions are used in offices, which can be classified by supply water temperature as low temperature room cooling units e.g., fan coil units and high temperature units, such as thermally active building systems (TABS), passive cooling beams, or active cooling beams, combined with ventilation supply air terminals [31,32].In low energy buildings, high temperature cooling is usually preferred to achieve higher energy efficiency for cooled water production by cooling plants [32].The performance of these systems is extensively analyzed in various recent studies.Most of the research is based on either computer simulations, mainly computational fluid dynamics (CFD) studies or studies conducted in controlled laboratory environments [33][34][35][36][37][38][39][40][41][42][43][44][45][46][47].The research in real office settings is mainly focused on buildings located in warm and hot climate countries, dominated by cooling need [48][49][50][51].To the knowledge of the authors, only few extensive studies have been carried out in cold and temperate climates and in low energy buildings.In Germany, Pfafferott, Herkel, Kalz, and Zeuschner [14] have conducted research on summertime TC in 12 low energy office buildings which are passively cooled with local heat sink based TABS.Results showed, that 41% of occupants were dissatisfied with thermal environment in summer, but assessment, according to the standard CEN EN 15251 [52], showed measured indoor temperature-based classification relative to the indoor climate category I (highest) and II, indicating a gap between perceived and assessed TC conditions and the need for more detailed comfort assessment.Hens [15] investigated TC in two office buildings in Belgium cooled with active chilled beams and air-cooling systems.He found that the Fanger [53] predicted mean vote (PMV)/predicted percentage dissatisfied (PPD) curve underestimated the actual number of dissatisfied occupants and that standards should not be considered as absolute references.It was also concluded that one should be very careful when interpreting the results of TC studies.
The theoretical knowledge involved or access to expert engineers during office building design is feasible and implemented in practice in Estonia.However, the volatile quality of different parts of the design, building phase simplifications with budget cuts and the contradiction between the initial task and the actual situation leads to a risk of an outcome failure regarding TC for the occupant.In the Estonian construction market, great emphasis is given on diplomas, professional certificates, software for both building information modeling (BIM), and product selection programs integrating BIM solutions in the building process.In reality, during the construction process the HVAC designer and after the warranty period, the constructor retreat.Therefore, in a short period of time during the design a huge effort is invested in the project definition, while after the realization phase much of the expert advice is ignored.Otherwise, complaints regarding draught or room temperature were not topical issues.In Estonia as well, in-depth research on cooling season TC and occupant satisfaction is practically non-existent, a few studies in office buildings have been conducted with the main focus on heating season performance and mostly aimed towards energy efficiency analysis.The conducted studies indicate problems and dissatisfaction with thermal environment but lack the detail to specify the causes and details of occupants' thermal conditions and HVAC systems performance in terms of room equipment.
Regardless of the design and performance of the HVAC systems installed in actual open offices, the hypothesis of this study proposes that high proportion of occupants are dissatisfied with the TC conditions.The goal of this study is to determine, whether the thermal sensation dissatisfaction of the occupants in modern office spaces is verifiable and in accordance with the valid standard criteria.This paper aims to fill the gap of summer TC assessment by extensive field studies and thorough occupant survey in modern office buildings in Estonia, a temperate climate country.We have investigated four recently constructed and one reconstructed office buildings with open plan office layouts designed with different ventilation and cooling solutions, including mixing and displacement ventilation, TABS, radiant cooling panels, fan coil units, and active cooling beams.The on-site measurements conducted in the offices consist of high resolution and accuracy temperature and air velocity measurements with DR and TC calculations, which are described in the following chapter.

Methods
The flow chart of research methodology is shown in Figure 1.Section of methods is divided between description of reference objects, measurement set-up and equipment specifications, data analysis, and indoor climate questionnaire (ICQ).We used standard-based [52,54,55] methods in this study to measure and calculate TC parameters and to perform an online ICQ survey.The TC measuring probe and tripod mobile and flexible kit set [56] we used was designed for research and development purposes.
Atmosphere 2019, 10, x FOR PEER REVIEW 3 of 21 software for both building information modeling (BIM), and product selection programs integrating BIM solutions in the building process.In reality, during the construction process the HVAC designer and after the warranty period, the constructor retreat.Therefore, in a short period of time during the design a huge effort is invested in the project definition, while after the realization phase much of the expert advice is ignored.Otherwise, complaints regarding draught or room temperature were not topical issues.In Estonia as well, in-depth research on cooling season TC and occupant satisfaction is practically non-existent, a few studies in office buildings have been conducted with the main focus on heating season performance and mostly aimed towards energy efficiency analysis.The conducted studies indicate problems and dissatisfaction with thermal environment but lack the detail to specify the causes and details of occupants' thermal conditions and HVAC systems performance in terms of room equipment.Regardless of the design and performance of the HVAC systems installed in actual open offices, the hypothesis of this study proposes that high proportion of occupants are dissatisfied with the TC conditions.The goal of this study is to determine, whether the thermal sensation dissatisfaction of the occupants in modern office spaces is verifiable and in accordance with the valid standard criteria.This paper aims to fill the gap of summer TC assessment by extensive field studies and thorough occupant survey in modern office buildings in Estonia, a temperate climate country.We have investigated four recently constructed and one reconstructed office buildings with open plan office layouts designed with different ventilation and cooling solutions, including mixing and displacement ventilation, TABS, radiant cooling panels, fan coil units, and active cooling beams.The on-site measurements conducted in the offices consist of high resolution and accuracy temperature and air velocity measurements with DR and TC calculations, which are described in the following chapter.

Methods
The flow chart of research methodology is shown in Figure 1.Section of methods is divided between description of reference objects, measurement set-up and equipment specifications, data analysis, and indoor climate questionnaire (ICQ).We used standard-based [52,54,55] methods in this study to measure and calculate TC parameters and to perform an online ICQ survey.The TC measuring probe and tripod mobile and flexible kit set [56] we used was designed for research and development purposes.

Reference Objects
General information regarding reference objects are provided in Table 1.The scope and range of measurement points with main building envelope characteristics, such as thermal transmittance for main surfaces, such as walls, windows, floors on ground and roofs are listed with specific heat loss of external envelopes and window-to-wall ratios.

Reference Objects
General information regarding reference objects are provided in Table 1.The scope and range of measurement points with main building envelope characteristics, such as thermal transmittance for main surfaces, such as walls, windows, floors on ground and roofs are listed with specific heat loss of external envelopes and window-to-wall ratios.In Buildings A, B, C, and D, measurements were also taken on the highest floor and in Buildings B and D on the lowest floor.The temperature of slabs was considered to be close to t i and therefore the impact on operative temperature was not accounted for, as heat transmission through the building envelope in such low energy buildings is negligible compared to the heat gains through glazed surfaces and have little effect on TC.A variety of HVAC systems was involved in measurements zones (Table 2) including new and innovative solutions in the Estonian construction market.The buildings involved in this study were chosen from a range of modern office spaces in Tallinn.First criteria for reference objects was the correspondence with the Estonian energy efficiency regulations, which were first set in 2007 [57].This created the prerequisites for new buildings and HVAC systems criteria, such as envelope related parameters, such as air tightness, external wall insulation thickness, window glazing solutions, and HVAC system parameters, e.g., effectiveness of room units and energy sources, heat recovery effectiveness, specific fan power of ventilation units etc. Buildings A and B have high temperature heating systems and district heating.Building B is using low temperature heating and a ground source heat pump, Building D has high temperature heating water produced and a gas boiler and electrical heating convectors are installed in Building E. All of the studied buildings are equipped with dedicated outdoor air ventilation systems with heat recovery.Ventilation air distribution methods were classified as mixing ventilation, except for Building B, where supply air systems were built in a way to support displacement ventilation method.Buildings A and C were using active chilled beams for supply air distribution.Buildings A, C, and D are built with chillers to supply the cooling system.In all Buildings, except for E, high temperature cooling is used in room conditioning units as supply air is dehumidified in the air handling units.Multi-split fan coil units with refrigerant without the option of heating function were in operation in Building E. Room conditioning units in Buildings C and D, including the Building B with thermally active buildings system, operated both for heating and cooling purposes.

Measurement Equipment
Experimental measurements in this study were carried out with a TC measurement system Dantec Dynamics ComfortSense [56].The system is designed for high quality multi-point measurements of v a , t i , RH, and t o .The set is equipped with software, what allows easy setup for measuring sequence and positions giving researchers a comprehensive overview if the measured data.Measurement equipment probe data is described in Table 3. Ø: diameter of air diffuser connection duct.
The buildings involved in this study were chosen from a range of modern office spaces in Tallinn.First criteria for reference objects was the correspondence with the Estonian energy efficiency regulations, which were first set in 2007 [57].This created the prerequisites for new buildings and HVAC systems criteria, such as envelope related parameters, such as air tightness, external wall insulation thickness, window glazing solutions, and HVAC system parameters, e.g., effectiveness of room units and energy sources, heat recovery effectiveness, specific fan power of ventilation units etc. Buildings A and B have high temperature heating systems and district heating.Building B is using low temperature heating and a ground source heat pump, Building D has high temperature heating water produced and a gas boiler and electrical heating convectors are installed in Building E. All of the studied buildings are equipped with dedicated outdoor air ventilation systems with heat recovery.Ventilation air distribution methods were classified as mixing ventilation, except for Building B, where supply air systems were built in a way to support displacement ventilation method.Buildings A and C were using active chilled beams for supply air distribution.Buildings A, C, and D are built with chillers to supply the cooling system.In all Buildings, except for E, high temperature cooling is used in room conditioning units as supply air is dehumidified in the air handling units.Multi-split fan coil units with refrigerant without the option of heating function were in operation in Building E. Room conditioning units in Buildings C and D, including the Building B with thermally active buildings system, operated both for heating and cooling purposes.

Measurement Equipment
Experimental measurements in this study were carried out with a TC measurement system Dantec Dynamics ComfortSense [56].The system is designed for high quality multi-point measurements of va, ti, RH, and to.The set is equipped with software, what allows easy setup for measuring sequence and positions giving researchers a comprehensive overview if the measured data.Measurement equipment probe data is described in Table 3.
Table 3. Specifications of measuring equipment [56].The set is mounted on a tripod including five draft probes, one humidity and one to probe.For a sitting position, ISO standard [55] recommends measuring heights for ankle level 0.1 m, abdomen level 0.6 m, and head level 1.1 m.Conformably to Fanger and Christensen [6], mean va and standard deviation at three heights around the sitting occupant body were measured according to heights shown on Figure 2. RH probe was set at 1.0 m as a fixed height for measuring has not been fixed for front of windows.Installed room unit heating power 60 W/m².mounted in the suspended ceiling (height 2.5-2.7 m).
mounted in the suspended ceiling (height 2.7 m). Installed total cooling power 78 W/m².(Ventilation supply air is not chilled) Ø: diameter of air diffuser connection duct.
The buildings involved in this study were chosen from a range of modern office spaces in Tallinn.First criteria for reference objects was the correspondence with the Estonian energy efficiency regulations, which were first set in 2007 [57].This created the prerequisites for new buildings and HVAC systems criteria, such as envelope related parameters, such as air tightness, external wall insulation thickness, window glazing solutions, and HVAC system parameters, e.g., effectiveness of room units and energy sources, heat recovery effectiveness, specific fan power of ventilation units etc. Buildings A and B have high temperature heating systems and district heating.Building B is using low temperature heating and a ground source heat pump, Building D has high temperature heating water produced and a gas boiler and electrical heating convectors are installed in Building E. All of the studied buildings are equipped with dedicated outdoor air ventilation systems with heat recovery.Ventilation air distribution methods were classified as mixing ventilation, except for Building B, where supply air systems were built in a way to support displacement ventilation method.Buildings A and C were using active chilled beams for supply air distribution.Buildings A, C, and D are built with chillers to supply the cooling system.In all Buildings, except for E, high temperature cooling is used in room conditioning units as supply air is dehumidified in the air handling units.Multi-split fan coil units with refrigerant without the option of heating function were in operation in Building E. Room conditioning units in Buildings C and D, including the Building B with thermally active buildings system, operated both for heating and cooling purposes.

Measurement Equipment
Experimental measurements in this study were carried out with a TC measurement system Dantec Dynamics ComfortSense [56].The system is designed for high quality multi-point measurements of va, ti, RH, and to.The set is equipped with software, what allows easy setup for measuring sequence and positions giving researchers a comprehensive overview if the measured data.Measurement equipment probe data is described in Table 3.
Table 3. Specifications of measuring equipment [56].The set is mounted on a tripod including five draft probes, one humidity and one to probe.For a sitting position, ISO standard [55] recommends measuring heights for ankle level 0.1 m, abdomen level 0.6 m, and head level 1.1 m.Conformably to Fanger and Christensen [6], mean va and standard deviation at three heights around the sitting occupant body were measured according to heights shown on Figure 2. RH probe was set at 1.0 m as a fixed height for measuring has not been fixed for Ø: diameter of air diffuser connection duct.
The buildings involved in this study were chosen from a range of modern office spaces in Tallinn.First criteria for reference objects was the correspondence with the Estonian energy efficiency regulations, which were first set in 2007 [57].This created the prerequisites for new buildings and HVAC systems criteria, such as envelope related parameters, such as air tightness, external wall insulation thickness, window glazing solutions, and HVAC system parameters, e.g., effectiveness of room units and energy sources, heat recovery effectiveness, specific fan power of ventilation units etc. Buildings A and B have high temperature heating systems and district heating.Building B is using low temperature heating and a ground source heat pump, Building D has high temperature heating water produced and a gas boiler and electrical heating convectors are installed in Building E. All of the studied buildings are equipped with dedicated outdoor air ventilation systems with heat recovery.Ventilation air distribution methods were classified as mixing ventilation, except for Building B, where supply air systems were built in a way to support displacement ventilation method.Buildings A and C were using active chilled beams for supply air distribution.Buildings A, C, and D are built with chillers to supply the cooling system.In all Buildings, except for E, high temperature cooling is used in room conditioning units as supply air is dehumidified in the air handling units.Multi-split fan coil units with refrigerant without the option of heating function were in operation in Building E. Room conditioning units in Buildings C and D, including the Building B with thermally active buildings system, operated both for heating and cooling purposes.

Measurement Equipment
Experimental measurements in this study were carried out with a TC measurement system Dantec Dynamics ComfortSense [56].The system is designed for high quality multi-point measurements of va, ti, RH, and to.The set is equipped with software, what allows easy setup for measuring sequence and positions giving researchers a comprehensive overview if the measured data.Measurement equipment probe data is described in Table 3.
Table 3. Specifications of measuring equipment [56].The set is mounted on a tripod including five draft probes, one humidity and one to probe.For a sitting position, ISO standard [55] recommends measuring heights for ankle level 0.1 m, abdomen level 0.6 m, and head level 1.1 m.Conformably to Fanger and Christensen [6], mean va and standard deviation at three heights around the sitting occupant body were measured according to heights shown on Figure 2. RH probe was set at 1.0 m as a fixed height for measuring has not been fixed for The set is mounted on a tripod including five draft probes, one humidity and one t o probe.For a sitting position, ISO standard [55] recommends measuring heights for ankle level 0.1 m, abdomen level 0.6 m, and head level 1.1 m.Conformably to Fanger and Christensen [6], mean v a and standard deviation at three heights around the sitting occupant body were measured according to heights shown on Figure 2. RH probe was set at 1.0 m as a fixed height for measuring has not been fixed for measurements.The t o probe was mounted with the angle of 30 • at the height of 0.6 m as the abdomen level of a sitting person [55].Probes were connected with 54N90 ComfortSense main frame [56], using 7 channels of 16.Main frame was in turn connected with laptop computer where the measurement data was stored using ComfortSense software version 4, (Dantec Dynamics A/S, Skovlunde, Denmark) [56].Measurement period of 180 seconds as the least time recommended [59] was used.

Data Analysis
Measurement data, including ti, va, RH, and to was recorded with the sampling rate 20 Hz with ComfortSense [56] and processed in Microsoft Excel.TC parameters are calculated for each measurement positions with equations for Tu, DR, PMV, and PPD followed.To assess DR, the fluctuation rate of va is described as Tu, which is calculated by [59] where SD is standard deviation of measured local mean va (m/s) for one measurement.With ti, va, and Tu, the percentage of people predicted to be dissatisfied because of draught may be calculated as [60] DR = (34 -ti) × (va -0.05) 0.62 × (0.37 × va × Tu + 3.14 To predict the mean value of the subjective ratings of a group of people in a given environment, PMV index is used.Consisting of a set of parameters with sub-formulas, the PMV equation is given by [60] PMV The PMV index in Equation (3) was calculated using Equations ( 4)- (11).In the equations provided, M (W/m²) is metabolic rate and W (W/m²) is the effective mechanical power.Assumption of metabolic rate 1.2 met for sedentary activity for summer season provided in EN 16798-1 was used.Sedentary activity does not suppose producing effective mechanical power, therefore 0 (W/m²) was used in analysis.The next symbol Hd in Equation (3) represents dry heat loss, which is found as where mtsk is mean skin temperature [°C] and in Equation ( 4), tcl is expressed using to and calculated through iterative process, by In Equation ( 5  Probes were connected with 54N90 ComfortSense main frame [56], using 7 channels of 16.Main frame was in turn connected with laptop computer where the measurement data was stored using ComfortSense software version 4, (Dantec Dynamics A/S, Skovlunde, Denmark) [56].Measurement period of 180 seconds as the least time recommended [59] was used.

Data Analysis
Measurement data, including t i , v a , RH, and t o was recorded with the sampling rate 20 Hz with ComfortSense [56] and processed in Microsoft Excel.TC parameters are calculated for each measurement positions with equations for Tu, DR, PMV, and PPD followed.To assess DR, the fluctuation rate of v a is described as Tu, which is calculated by [59] where SD is standard deviation of measured local mean v a (m/s) for one measurement.With t i , v a , and Tu, the percentage of people predicted to be dissatisfied because of draught may be calculated as [60] DR = (34 To predict the mean value of the subjective ratings of a group of people in a given environment, PMV index is used.Consisting of a set of parameters with sub-formulas, the PMV equation is given by [60] PMV The PMV index in Equation (3) was calculated using Equations ( 4)- (11).In the equations provided, M (W/m 2 ) is metabolic rate and W (W/m 2 ) is the effective mechanical power.Assumption of metabolic rate 1.2 met for sedentary activity for summer season provided in EN 16798-1 was used.Sedentary activity does not suppose producing effective mechanical power, therefore 0 (W/m 2 ) was used in analysis.The next symbol H d in Equation (3) represents dry heat loss, which is found as where mt sk is mean skin temperature [ • C] and in Equation (4), t cl is expressed using t o and calculated through iterative process, by In Equation (5), where v ar was set equal to the v a as occupants were intended to be stationary sensing draught.Equation (5) includes f cl , which is calculated by Continuing the PMV index calculation, in Equation ( 3), evaporative heat exchange at the skin, when the person experiences a sensation of thermal neutrality E c given as where p a is the water vapor partial pressure [Pa], calculated using measured RH by In addition, Equation (3) for PMV includes respiratory convective heat exchange C res , calculated as and Equation (3) includes also respiratory evaporative heat exchange E res , given as Finally, to predict the rate of people dissatisfied in a thermal environment, the PPD index is used.Knowing PMV, PPD can be calculated as [60] Measured values are shown in Building result figures in the results chapter and used in Equations ( 1) and (2) for calculating Tu and DR, and in Equations ( 3) and ( 12) to calculate PMV and PPD.

Indoor Climate Questionnaire
To study occupant satisfaction we provided online questionnaires to the employees of the measured office spaces.As some organizations involved in this study are moving towards policy of a paperless work management, we used Google Forms [61] application.In addition to standard CEN EN 15251 [52] suggestions, we added also questions about age, gender, amount of time behind the desk during workday, and the working environment regarding cabinet or open office plan.The ICQ is presented in Appendix A.

On-Site Measurements
This section provides an overview of the TC measurement time and weather information (Table 4), followed by measurement results with calculated TC indicative parameters Tu, DR, PMV, and PPD.ICQ survey results are summarized at the end of the results sections.The experiments were carried out on regular workdays during August.Measurements were taken by two persons, by the main author of this article assisted by graduate students in different buildings.HVAC systems were in regular performance mode without disfunctions or failures recorded.Internal gains by occupants, office equipment, and lighting were in use by default as some desks were empty by unused space, duties, or vacation.No serious defects in HVAC design or construction were observed.Although, some air flow and velocity aspects were noticeable.As in Buildings A and C, active beams were in use, occupants were not always placed sitting according to rule of thumbs, according to the architectural layout, or number of persons.Possible air flow obstacles by lighting fixture (Figure 3a) were noticed with open ceiling in Building A. DR risk was also predictable in building E (Figure 3b) where some vanes were taped to closed position.DR risk was more carefully considered in Buildings B and D.
Atmosphere 2019, 10, x FOR PEER REVIEW 8 of 21 The experiments were carried out on regular workdays during August.Measurements were taken by two persons, by the main author of this article assisted by graduate students in different buildings.HVAC systems were in regular performance mode without disfunctions or failures recorded.Internal gains by occupants, office equipment, and lighting were in use by default as some desks were empty by unused space, duties, or vacation.No serious defects in HVAC design or construction were observed.Although, some air flow and velocity aspects were noticeable.As in Buildings A and C, active beams were in use, occupants were not always placed sitting according to rule of thumbs, according to the architectural layout, or number of persons.Possible air flow obstacles by lighting fixture (Figure 3a) were noticed with open ceiling in Building A. DR risk was also predictable in building E (Figure 3b) where some vanes were taped to closed position.DR risk was more carefully considered in Buildings B and D.

Results
Based on on-site measurements, the summary of va in each measurement position are shown below for each Building.According to three heights provided in Figure 2, va values during measurement period are shown with box and whiskers plot.Minimum and maximum are at the end of the whiskers, the lower and the upper line of the box are first and third quartiles, the line between is median and the cross shows mean va value of the measurement in one position.
On the box and whiskers plot, the category of the indoor climate category is colored according to the lowest criteria achieved during measurements, meaning if one of the three height is in III category, the measurement point is placed in the least, III category.In the table part on the result figures below the box and whiskers plot, measured values and calculated parameters are colored according to the category reached to be more easily distinguishable.Tu and RH measured values are not colored as being not categorized.Nonetheless, measurement point indoor climate category is defined by the inferior measured value or calculated parameter reached altogether.

Results
Based on on-site measurements, the summary of v a in each measurement position are shown below for each Building.According to three heights provided in Figure 2, v a values during measurement period are shown with box and whiskers plot.Minimum and maximum are at the end of the whiskers, the lower and the upper line of the box are first and third quartiles, the line between is median and the cross shows mean v a value of the measurement in one position.
On the box and whiskers plot, the category of the indoor climate category is colored according to the lowest criteria achieved during measurements, meaning if one of the three height is in III category, the measurement point is placed in the least, III category.In the table part on the result figures below the box and whiskers plot, measured values and calculated parameters are colored according to the category reached to be more easily distinguishable.Tu and RH measured values are not colored as being not categorized.Nonetheless, measurement point indoor climate category is defined by the inferior measured value or calculated parameter reached altogether.

Building A Results
The v a results and TC parameters in Building A equipped with open ceiling active chilled beams are provided below in Figure 4.In Building A, in 2/3 of the measured positions the v a was below the first indoor climate category threshold.Five positions met the II category requirement and in one position the v a was above the category II threshold.Measurement No 14 was taken in an office space with unusually high internal gains, where also multi-split fan coil units were additionally added to the environment due to the specifics of the lessee.The results of t i and v a including DR, PMV, and PPD are placed in the first category mainly.

Building B Results
The va results and TC parameters of Building B with slab-based TABS system are given below in Figure 5. Building B had more measured points in the second category by PMV and PPD compared to Building A. DR met the II category in four measurement positions.Positions 4-8 were in an office, where the ventilation rate had been doubled by the request of the lessee.These four measurements stand out above the others.Regarding the other four buildings observed, displacement ventilation effect can be seen, as va fluctuates more near the floor.

Building B Results
The v a results and TC parameters of Building B with slab-based TABS system are given below in Figure 5. Building B had more measured in the second category by PMV and PPD compared to Building A. DR met the II category in four measurement positions.Positions 4-8 were in an office, where the ventilation rate had been doubled by the request of the lessee.These four measurements stand out above the others.Regarding the other four buildings observed, displacement ventilation effect can be seen, as v a fluctuates more near the floor.

Building C Results
Building C was equipped with suspended ceiling active chilled beams and the results of va and parameters of TC are presented below in Figure 6.PMV, PPD, and ti were similar to Buildings A and B, at the same time va and DR were measured at two positions in the II category and three times in the III category.The va is more fluctuating on the height of the sitting person neck.m/s 0.07 0.07 0.12 0.11 0.11 0.06 0.10 0.12 0.16 0.06 0.13 0.12 0.14 0.05 0.07 0.14 0.08 0.08 0.09

Building D Results
Equipped with radiant cooling panels, results of va and parameters of TC in Building D are showed below in Figure 7. Compared to other buildings, Building D with the least number of positions had the best results on all analyzed parameters.In all cases, I category DR was achieved.At all times, mean va remained below 0.10 m/s being more fluctuating near the floor.

Building D Results
Equipped with radiant cooling panels, results of v a and parameters of TC in Building D are showed below in Figure 7. Compared to other buildings, Building D with the least number of positions had the best results on all analyzed parameters.In all cases, I category DR was achieved.At all times, mean v a remained below 0.10 m/s being more fluctuating near the floor.

Building E Results
According to the results, Building E achieved the worst TC values by categories.DR was in the II category in 4 positions of 14, ti was in III category four times.PMV and PPD second category was not reached 5 times.Fluctuations of va were random depending on the height.The va results and TC parameters in Building E, with fan coil units mounted in the suspended ceiling, are compared below in Figure 8.

Building E Results
According to the results, Building E achieved the worst TC values by categories.DR was in the II category in 4 positions of 14, t i was in III category four times.PMV and PPD second category was not reached 5 times.Fluctuations of v a were random depending on the height.The v a results and TC parameters in Building E, with fan coil units mounted in the suspended ceiling, are compared below in Figure 8.

Results of the Indoor Climate Questionnaire
Based on ICQ survey, summary of the results for thermal environment are shown below in Figure 9, the ICQ results for PMV and PPD are presented below in Figure 10.The highest number of answers were in the Building A with 36 responses divided between all age groups equally between men and women.A total of 83% were working in open office layout and 86% were spending most of the day at their workplace.For 83% of the respondents, ti was described as suitable.Meanwhile, 6 occupants found it to be warm and 7 slightly cooler.A total of 89% had not or had perceived slight odor, 72% did not find lighting fixtures or sunlight to be disturbing, and 81% found ICQ to be suitable or better.A total of 61% perceived overall acoustics and 36% perceived other noises to be disturbing.Roughly half of the respondents rarely felt eye problems, headaches, or concentration matters and

Results of the Indoor Climate Questionnaire
Based on ICQ survey, summary of the results for thermal environment are shown below in Figure 9, the ICQ results for PMV and PPD are presented below in Figure 10.The highest number of answers were in the Building A with 36 responses divided between all age groups equally between men and women.A total of 83% were working in open office layout and 86% were spending most of the day at their workplace.For 83% of the respondents, t i was described as suitable.Meanwhile, 6 occupants found it to be warm and 7 slightly cooler.A total of 89% had not or had perceived slight odor, 72% did not find lighting fixtures or sunlight to be disturbing, and 81% found ICQ to be suitable or better.A total of 61% perceived overall acoustics and 36% perceived other noises to be disturbing.Roughly half of the respondents rarely felt eye problems, headaches, or concentration matters and 64% rarely felt nasal or throat irritation.Extra comments mentioned occasional lack of ventilation and air dryness.Number of respondents of the ICQ is below the least recommended sample size [63], therefore the results of the ICQ include higher uncertainty (Figure 10).Thermal sensation voted by occupants covers significantly wider range than PMV calculated from measurements.Majority of the respondents were working in open office.The most unsatisfying ti was in the Building E and the most suitable ti was in the Building D. In general, unsuitable ti was perceived more as warmer than cooler.In Buildings A, B, and D the ti was perceived suitable for over 80% of the employees, while it was 67% in the Building E and the 60% in the Building C.

Discussion
The on-site measurement results showed, that the during cooling summertime DR risk can be stated in all observed buildings.Preconception of avoiding fan coil units for cooling does not immediately guarantee a superior thermal environment without draught.However, draught risk was the lowest in Building D with radiant cooling panels as room conditioning units.Number of respondents of the ICQ is below the least recommended sample size [63], therefore the results of the ICQ include higher uncertainty (Figure 10).Thermal sensation voted by occupants covers significantly wider range than PMV calculated from measurements.Majority of the respondents were working in open office.The most unsatisfying ti was in the Building E and the most suitable ti was in the Building D. In general, unsuitable ti was perceived more as warmer than cooler.In Buildings A, B, and D the ti was perceived suitable for over 80% of the employees, while it was 67% in the Building E and the 60% in the Building C.

Discussion
The on-site measurement results showed, that the during cooling summertime DR risk can be stated in all observed buildings.Preconception of avoiding fan coil units for cooling does not immediately guarantee a superior thermal environment without draught.However, draught risk was the lowest in Building D with radiant cooling panels as room conditioning units.Respondents in the Building B were 38% females, 2/3 aged between 26-35 or 36-45 and 1/2 spending half of the workday behind the desk.A total of 72% of them working in open office environment.Ninety percent found t i to be suitable.A total of 13 of the 29 respondents did not perceive odor.Lighting was disturbing for 21% and sunlight for 14%, meanwhile 14% were dissatisfied with ICQ.Seven percent did not find room acoustics and 17% general noise in office to be disturbing.Half of the respondents had rarely felt eye dryness or irritation, occur headaches or fatigue, and felt nasal problem or dry throat.A total of 62% had rarely felt concentration problems.
Seventy percent of the 20 ICQ respondents in Building C were women.Answers were divided between the age of 26 to 65 with the majority of them working in open office landscape, 2/3 working behind their desk most of the day.Perceived as too warm by 20%, t i was suitable by 75% of the occupants.Ninety percent had not perceived or had perceived slight odor.A total of 70% did not find lighting equipment to be disturbing and 75% was not disturbed by the sunlight.Forty percent of the respondents found air quality to be not suitable or unacceptable.A total of 85% perceived colleagues' speech and overall room acoustic to be somewhat disturbing, while 65% claimed other noises to be distracting.A total of 1/3 had rarely felt eye problems, occurred headaches, or tiredness.A total of 45% had rarely felt nasal or throat irritation and 20% had rarely had concentration issues.Extra comments mentioned lower fresh air rate in the end of the day.
Building D had only 8 responses for the online ICQ all of them working in the open office.For the majority of the answers, t i was suitable.Odor was rarely noticed, lighting or sunlight was not disturbing.ICQ was suitable or better, while room acoustics was more disturbing than other noises.Nasal issues were more often to occur compared to eye dryness or headaches and concentration issues.Extra comments noted that open office may be cheaper option for the employer being unsuitable for the employees.
A total of 2/3 of the 22 respondents in Building E were in the second age group between 26-35 years and 36% in overall were females.A total of 77% of the tenants were working in an open office environment, while 2/3 of them were spending most of their day behind the desk.One-third found t i to be suitable and 2/3 claimed the t i to be slightly warm, warm, or hot.Fifty percent perceived weak or moderate odor.Room lighting equipment did not disturb 82% and the sun did not disturb 60% of the respondents.A total of 2/3 marked ICQ suitable, good, or very good.Room acoustic level was not claimed to be disturbing for 40% and other noise for 23% of the respondents.Fifty percent had rarely felt eye dryness or irritation, 64% had rarely occur headaches or fatigue, 82% had rarely felt nasal problems or dry throat, and 50% mentioned concentration issues sometimes, often, or all the time.Extra comments noted that air quality decreases in the second phase of the day and the missing option for opening windows was also described as a disadvantage.
Number of respondents of the ICQ is below the least recommended sample size [63], therefore the results of the ICQ include higher uncertainty (Figure 10).Thermal sensation voted by occupants covers significantly wider range than PMV calculated from measurements.Majority of the respondents were working in open office.The most unsatisfying t i was in the Building E and the most suitable t i was in the Building D. In general, unsuitable t i was perceived more as warmer than cooler.In Buildings A, B, and D the t i was perceived suitable for over 80% of the employees, while it was 67% in the Building E and the 60% in the Building C.

Discussion
The on-site measurement results showed, that the during cooling summertime DR risk can be stated in all observed buildings.Preconception of avoiding fan coil units for cooling does not immediately guarantee a superior thermal environment without draught.However, draught risk was the lowest in Building D with radiant cooling panels as room conditioning units.
Possible causes, v a and DR was not significantly higher in the case of fan coil units in Building E was the taping of air distribution vanes (Figure 3b) and also positioning of the working stations was carried out avoiding direct draught from the fan coil units.This could explain the higher thermal environment temperatures.The induced airflow rate is manually adjustable for open ceiling active chilled beams in Building A and was adjusted into different positions for avoiding possible draught between two beams in various places.In Building C, few suspended ceiling active chilled beams had paper covers blocking air flow from the nozzles.These modifications were made due to the complaints, decrease in productivity or spatial plan and the layout of the workspaces.Described modifications in Buildings A, C, and E refer to possible ineffective floor space areas.Therefore, whether the design or construction may have been inaccurate or user-based thermal environment setpoints do not meet the requirements for v a and DR.
The v a limit values in EN 16798−1:2019 [54] have been calculated assuming t o +23 • C and Tu 40%.Figure 11a illustrates that the Tu is considerably higher than the default value, which increases the unsatisfaction with local TC.However, the measured t o was higher than the default value in most of the measured positions in all buildings, which decreases the number of dissatisfied.Figure 11b shows that, in general, the DR calculated based on measured t o and Tu is in the same scale with the one calculated with the default values.
paper covers blocking air flow from the nozzles.These modifications were made due to the complaints, decrease in productivity or spatial plan and the layout of the workspaces.Described modifications in Buildings A, C, and E refer to possible ineffective floor space areas.Therefore, whether the design or construction may have been inaccurate or user-based thermal environment setpoints do not meet the requirements for va and DR.
The va limit values in EN 16798−1:2019 [54] have been calculated assuming to +23 °C and Tu 40%.Figure 11a illustrates that the Tu is considerably higher than the default value, which increases the unsatisfaction with local TC.However, the measured to was higher than the default value in most of the measured positions in all buildings, which decreases the number of dissatisfied.Figure 11b shows that, in general, the DR calculated based on measured to and Tu is in the same scale with the one calculated with the default values.In further analysis of this study ti will be more deeply discussed, foreseeing to include transitional period and heating period measurements, façade inspection and ti periodical data analysis in the reference buildings.Therefore, the performance of the cooling units according to ti could be more clearly presented by period or duration curve.Periodical data analysis on ti is mandatory as ti presented in this study reflects only a fragment of the thermal environment.
Positioning TC measurement values on periodical ti data can indicate TC measurement accuracy and dispersion.IQC survey number of respondents also needs additional attention, how to achieve a higher response rate.
There are several limitations to this work.Authors had no control over the boundary conditions during measurements.This study only focuses on a few office spaces in five different building in Tallinn.More further studies of actual work environment need to be performed in order to be able to draw general conclusions about studied room conditioning solutions air distribution performance.

Conclusions
This study was based on TC measurements in open office environments in Tallinn.First or second category measured general thermal comfort in four buildings out of five were still inconvenient for significant number of occupants because of local thermal discomfort caused by In further analysis of this study t i will be more deeply discussed, foreseeing to include transitional period and heating period measurements, façade inspection and t i periodical data analysis in the reference buildings.Therefore, the performance of the cooling units according to t i could be more clearly presented by period or duration curve.Periodical data analysis on t i is mandatory as t i presented in this study reflects only a fragment of the thermal environment.Positioning TC measurement values on periodical t i data can indicate TC measurement accuracy and dispersion.IQC survey number of respondents also needs additional attention, how to achieve a higher response rate.
There are several limitations to this work.Authors had no control over the boundary conditions during measurements.This study only focuses on a few office spaces in five different building in Tallinn.More further studies of actual work environment need to be performed in order to be able to draw general conclusions about studied room conditioning solutions air distribution performance.

Conclusions
This study was based on TC measurements in open office environments in Tallinn.First or second category measured general thermal comfort in four buildings out of five were still inconvenient for significant number of occupants because of local thermal discomfort caused by draught and by some additional dissatisfaction indicated by questionnaires.Questionnaire survey showed deviation from predicted PPD in both directions.Some small occupant groups were either more satisfied or less satisfied at slightly cool or slightly warm thermal sensation, but at neutral sensation the results were more consistent.Less satisfied occupant groups exposed to higher air velocities has likely affected their satisfaction reported in thermal sensation questions because there was no specific draught question available.
Temperature measurements showed that air and operative temperature was the worst in Building E which was close to drop out from category III, while measurement results in Buildings A-D remained in between I and II category.According to the questionnaire over 80% of the employees in Buildings A, B, and C, and 75% in Building D were thermally satisfied.In the Building E, 59% of occupants found the thermal environment unsuitable or unacceptable.Generally, the average thermal satisfaction of occupants was well in line with the measurements.
Air velocity and draught rate measurements showed that modern offices do not necessarily reach to generally expected good indoor climate category II air velocity and draught rate values.A room conditioning solution with suspended ceiling active chilled beams in Building C, displacement ventilation in Building B with TABS and fan coil units in Building E showed category III performance only.Open ceiling active chilled beams in Building A corresponded to category II requirements and ceiling panels for room conditioning in Building D showed superior Category I performance.Category II and III results with active chilled beams indicate that dedicated air distribution solution together with proper design and sizing is needed to reach category II.
We found that existing standards do not provide enough detailed questionnaire for the assessment of occupant dissatisfaction.Our results suggest that questionnaire could be an easier compliance assessment method compared to measurements, which need an expensive equipment and carefully selected measurement days.For the compliance assessment with the measurement, there is more guidance needed especially how to select relevant measurement conditions and locations for draught rate measurement.Future office buildings with open-plan layouts revealed to be demanding environments where careful air distribution design is needed in order to meet comfort requirements.
Author Contributions: J.K. conceived and designed the experiments.M.K. prepared agreements with the building owners, performed the measurements and analyzed the data.M.T. and J.K. helped to perform the data analysis.M.K., R.S., M.T., and J.K. wrote this paper.All authors have read and agreed to the published version of the manuscript.

Acknowledgments:
The authors are grateful for the provided cooperation of the building owners, questionnaire respondents for their time and the valuable help from Tallinn University of Technology graduate students.

Conflicts of Interest:
The authors declare no conflict of interest.

Figure 1 .
Figure 1.Flow chart of research methodology.

Figure 1 .
Figure 1.Flow chart of research methodology.

Figure 2 .
Figure 2. Recommended air velocity probe heights behind the feet, elbow, and neck for a sitting person (a); Recommended operative temperature probe person's angle factor to their surroundings (b) [58].

Figure 2 .
Figure 2. Recommended air velocity probe heights behind the feet, elbow, and neck for a sitting person (a); Recommended operative temperature probe person's angle factor to their surroundings (b) [58].

Figure 3 .
Figure 3. Possible air flow obstacles with open active beam solution (a); modified airflow distribution with fan coil unit (b).

Figure 3 .
Figure 3. Possible air flow obstacles with open active beam solution (a); modified airflow distribution with fan coil unit (b).

Figure 4 .
Figure 4. Building A air velocity results in measurement points 1-18 and the thermal comfort parameters.

Figure 4 .
Figure 4. Building A air velocity results in measurement points 1-18 and the thermal comfort parameters.

Figure 5 .
Figure 5. Building B air velocity results in measurement points 1-16 and the thermal comfort parameters.

Figure 5 .
Figure 5. Building B air velocity results in measurement points 1-16 and the thermal comfort parameters.

Figure 6 .
Figure 6.Building C air velocity results in measurement points 1-19 and the thermal comfort parameters.

Figure 6 .
Figure 6.Building C air velocity results in measurement points 1-19 and the thermal comfort parameters.

Figure 7 .
Figure 7. Building D air velocity results in measurement points 1-11 and the thermal comfort parameters.

Figure 7 .
Figure 7. Building D air velocity results in measurement points 1-11 and the thermal comfort parameters.

Figure 8 .
Figure 8. Building E air velocity results in measurement points 1-14 and the thermal comfort parameters.

Figure 8 .
Figure 8. Building E air velocity results in measurement points 1-14 and the thermal comfort parameters.

Atmosphere 2019 , 21 Figure 9 .
Figure 9. Indoor climate questionnaire results for indoor air temperature.The descriptions of y-axis are the room air temperature sensation question (upper) and verification questions (middle and lower) from the indoor climate questionnaire (see Appendix A).

Figure 10 .
Figure 10.Indoor climate questionnaire results for predicted percentage of dissatisfied and predicted mean vote.

Figure 9 .
Figure 9. Indoor climate questionnaire results for indoor air temperature.The descriptions of y-axis are the room air temperature sensation question (upper) and verification questions (middle and lower) from the indoor climate questionnaire (see Appendix A).

Atmosphere 2019 , 21 Figure 9 .
Figure 9. Indoor climate questionnaire results for indoor air temperature.The descriptions of y-axis are the room air temperature sensation question (upper) and verification questions (middle and lower) from the indoor climate questionnaire (see Appendix A).

Figure 10 .
Figure 10.Indoor climate questionnaire results for predicted percentage of dissatisfied and predicted mean vote.

Figure 10 .
Figure 10.Indoor climate questionnaire results for predicted percentage of dissatisfied and predicted mean vote.

Figure 11 .
Figure 11.Air velocity and turbulence intensity results according to maximum air velocity categories I-III in summer (a); Draught rate correlation in measured and standard-based [54] conditions according to draught rate categories I-III (b).

Figure 11 .
Figure 11.Air velocity and turbulence intensity results according to maximum air velocity categories I-III in summer (a); Draught rate correlation in measured and standard-based [54] conditions according to draught rate categories I-III (b).

Funding:
This research was supported by the Estonian Centre of Excellence in Zero Energy and Resource Efficient Smart Buildings and Districts, ZEBE (grant 2014-2020.4.01.15-0016) funded by the European Regional Development Fund, by the programme Mobilitas Pluss (Grant No-2014-2020.4.01.16-0024,MOBTP88), by the European Commission through the H2020 project Finest Twins (grant No. 856602) and the Estonian Research Council grant (PSG409).

Table 1 .
General building information of reference objects.

Table 2 .
Heating, ventilation, and air conditioning room design solutions of reference objects.
Clothing unit 0.5 clo for summer season provided in EN 16798-1 was used in calculations.Equation (5) includes h c , which is given ash c = 2.38 × |t cl − t i | 0.25 for 2.38 × |t cl − t i | 0.25 > 12.1 × is the clothing insulation, f cl is the clothing surface area factor, v ar (m/s) is the relative air velocity, h c [W/(m 2 × K)] is the convective heat transfer coefficient, and t cl ( • C) is the clothing surface temperature.

Table 4 .
[62] of measurements and weather information from the Estonian Weather Service[62].