<figure><img src="image_1.png" /></figure>

ifip

International Federation for Information Processing

This document is the original author manuscript of a paper submitted to an IFIP
conference proceedings or other IFIP publication by Springer Nature. As such, there
may be some differences in the official published version of the paper. Such
differences, if any, are usually due to reformatting during preparation for publication or
minor corrections made by the author(s) during final proofreading of the publication
manuscript.
Realtime Multi-factor Dynamic Thermal
Comfort Estimation for Indoor
Environments

Georgia Tzitziou¹, Asimina Dimara¹,²✉, Alexios Papaioannou¹,³,
Christos Tzouvaras²,⁴, Stelios Krinidis¹,³, Christos-Nikolaos Anagnostopoulos²,
Dimosthenis Ioannidis¹, and Dimitrios Tzouvaras¹

¹ Information Technologies Institute (ITI), Centre for Research Technology (CERTH), Thessaloniki, Greece
{tzitzioug,adimara,alexiopa,krinidis,djoannid,Dimitrios.Tzouvaras}@iti.gr

² Department of Cultural Technology and Communication, Intelligent Systems Lab, University of the Aegean (UoA), Mytilene, Greece
canag@aegean.gr

³ Management Science and Technology Department, International Hellenic University (IHU), Kavala, Greece

⁴ WATT AND VOLT, Anonimi Etairia Ekmetalleysis, Enallaktikon Morfon Energeias, Athens, Greece
ch.tzouvaras@watt-volt.gr

**Abstract.** Thermal comfort models are mathematical representations that simulate the thermal environment and predict human comfort based on various factors such as air temperature, air velocity, relative humidity, and radiation heat transfer. These models are used to design and evaluate heating, ventilation, and air conditioning systems, buildings, and outdoor spaces. The main issue when exploiting predicted mean vote (PMV) and predicted percentage of dissatisfied (PPD) and model for thermal comfort estimation is how to estimate clothing insulation and metabolic rate as accurately as possible. In this paper, a novel approach for calculating thermal comfort is presented that combines algorithms to enhance the precision of existing approaches. Experimental results showcase the suggested method is more accurate than other approaches.

**Keywords:** Thermal comfort · Indoor environmental conditions · Personal factors

# 1 Introduction

Looking at people's daily lives, their timeless connection with the residence they have chosen to live in may be observed. One of the main factors influencing this interaction's quality is the thermal sensation inside the residence. In recent years, the factors that influence the internal temperature of a building and how they are taken into account to give us thermal comfort have been studied to a
large extent [1]. The importance of thermal comfort is extended even in early building stages like the design phase till further stages like renovation.

Thermal comfort, as defined by ISO standard 7730 [2], is the state of being satisfied with one's thermal environment. Thermal comfort models are mathematical representations that simulate the thermal environment and predict human comfort based on various factors such as air temperature, air velocity, relative humidity, and radiation heat transfer. These models use different methods and algorithms to estimate thermal comfort and are used to design and evaluate heating, ventilation, and air conditioning (HVAC) systems, buildings, and outdoor spaces. Some common thermal comfort models include the PMV (Predicted Mean Vote) and PPD (Predicted Percentage of Dissatisfied), the Adaptive model, and the Transfer Function model [3].

The research community is continually seeking the best thermal comfort model for estimating the ideal living conditions in buildings. Fanger's *PMV model (PredictedMeanVote)* is widely considered one of the most comprehensive and accurate thermal comfort models available. There are several reasons why the *PMV* model is considered better than other models [4]:

- Accounts for many factors that affect thermal comfort, such as air temperature, air velocity, relative humidity, and clothing insulation. This makes the *PMV* model more accurate and relevant compared to models that consider only a few factors.

- It is widely used and accepted in the residential sector and is established as the standard for the prediction of thermal comfort under the International ISO 7730 standard [2].

- It is relatively easy to use, with clear instructions for estimating the *PMV* index and thermal comfort based on index values.

- Has been proven to reliably estimate human thermal comfort in a wide range of situations, making it a dependable tool for evaluating and planning HVAC systems, buildings, and outdoor spaces.

These elements combine to make *PMV* an effective and trustworthy tool for forecasting thermal comfort, and these are some of the reasons why it is seen to be superior to other models. Based on ASHRAE, a fast and high-accuracy function has been developed that calculates the predicted mean vote (*PMV*) and the predicted percentage of satisfied (PPD) [5,6].

The input variables of the model, such as dry bulb air temperature, mean radiant temperature, average air speed, and relative humidity, can be measured accurately through sensor installation. However, there are no standard real-time methods to calculate the metabolic rate and clothing insulation of people habiting inside the building. The temperature inside a building greatly affects the type of clothing one wears [7]. Furthermore, the variation in clothing insulation values is attributable to the fact that what we wear is influenced by factors other than temperature, such as gender, age, and cold and heat tolerance [7]. In addition, temperature plays an important role in the determination of metabolic rate [7]. The behavior of building occupants is difficult to predict, model, or calculate due to the complexity of humans. Furthermore, the absence of established standards
and protocols for data collection, as well as the precision of data, pose difficulties in the field of building occupant research [8].

As a result, the main issue when exploiting the PMV and PPD model for thermal comfort estimation is how to estimate clothing insulation and metabolic rate as accurately as possible. Within this context, the aim of this paper is to present a novel approach by fusing widely used methodologies [7], one for the determination of clothing insulation and one for the determination of metabolic rate. The main objective is to mitigate the error in estimating personal factors by taking into account as many factors as possible, such as indoor and outdoor conditions.

The remainder of this paper is organized as follows: In Sect. 2 the way the thermal comfort is estimated is suggested, highlighting a novel approach for estimating the personal factors. Section 3 presents the experiment set up along with the results. Finally, in Sect. 4, conclusions are drawn.

# 2 Methodology

In this section, a novel approach for calculating thermal comfort is presented that combines algorithms to enhance the precision of existing approaches. A tool that would provide real-time, practical, and accurate thermal comfort estimation is suggested.

## 2.1 Thermal Comfort Inference

As addressed in the Introduction, Fanger's Predicted Mean Vote (PMV) and Predicted Percentage of Dissatisfied (PPD) model is widely used to estimate thermal comfort. The models are specified through a set of equations that are outlined below [9, 10]:

$$PPD = 100 - 95 \cdot e^{-(0.03353 \cdot PMV^4 + 0.2179 \cdot PMV^2)} \quad (1)$$

$$PMV = (0.303 \cdot e^{-0.036 \cdot M} + 0.028) \cdot L \quad (2)$$

$$L = M - W \quad (3)$$

$$M - W = C + R + E_{sk} + (C_{res} + E_{res}) \quad (4)$$

$$C = f_{cl} \cdot h_{c} \cdot (T_{cl} - T_{a}) \quad (5)$$

$$R = \sigma \cdot \varepsilon_{cl} \cdot f_{cl} \cdot F_{vf} \cdot \left[ (T_{cl} + 273.15)^4 - (T_{r} + 273.15)^4 \right] \quad (6)$$

$$C_{res} + E_{res} = 0.014 \cdot M \cdot (34 - T_a) + 0.0173 \cdot M \cdot (5.87 - P_a) \quad (7)$$

$$E_{sk} = 3.05 \cdot (5.73 - 0.007 \cdot M - P_a) + 0.42 \cdot (M - 58.15)) \quad (8)$$
All the used variables are briefly explained in Table 1 [13]. Finally, to estimate the thermal comfort based on the above equations the factors needed to estimate real-time PMV are [12]:

$$PMV = f(T_a, RH, I_{cl}, M), \qquad (9)$$

where $T_a$ refers to room temperature, RH refers to room humidity, $I_{cl}$ and $M$ refer to clothing insulation and metabolic rate respectively.

Room environmental conditions (i.e., $T_a$, RH) may be retrieved from indoor temperature and humidity sensors while estimating the personal factors may be proven to be a thorny problem [12]. Even though values for the personal factors may be retrieved from the ASHRAE table, an estimation of $I_{cl}$ and $M$ values based solely on the ASHRAE table may generate a significant mistake during thermal comfort estimation [14, 15]. As a result, another approach for estimating personal factors should be followed.

<table><caption>**Table 1.** PMV and PPD variables</caption><thead><tr><th>Variable</th><th>Meaning</th><th>Unit</th></tr></thead><tbody><tr><td>W</td><td>external work</td><td>W/m<sup>2</sup></td></tr><tr><td>M</td><td>metabolic rate (internal energy production)</td><td>W/m<sup>2</sup></td></tr><tr><td>C</td><td>heat loss by convection</td><td>W/m<sup>2</sup></td></tr><tr><td>R</td><td>heat loss by thermal radiation</td><td>W/m<sup>2</sup></td></tr><tr><td>E<sub>sk</sub></td><td>heat loss by evaporation from the skin</td><td>W/m<sup>2</sup></td></tr><tr><td>C<sub>res</sub></td><td>sensible heat loss due to respiration</td><td>W/m<sup>2</sup></td></tr><tr><td>E<sub>res</sub></td><td>heat loss by evaporation from the skin</td><td>W/m<sup>2</sup></td></tr><tr><td>T<sub>cl</sub></td><td>clothing surface temperature</td><td>°C</td></tr><tr><td>T<sub>a</sub></td><td>ambient air temperature (indoor)</td><td>°C</td></tr><tr><td>h<sub>c</sub></td><td>heat transfer coefficient</td><td>W/m<sup>2</sup>·K</td></tr><tr><td>V<sub>a</sub></td><td>air velocity</td><td>m/s</td></tr><tr><td>f<sub>cl</sub></td><td>clothing area factor</td><td>clo</td></tr><tr><td>ε<sub>cl</sub></td><td>emmisivity of clothing</td><td></td></tr><tr><td>F<sub>vf</sub></td><td>view factor between the body and the surrounding</td><td></td></tr><tr><td>σ</td><td>Stefan Boltzmann constant [11]</td><td>W/m<sup>2</sup>·K<sup>4</sup></td></tr><tr><td>T<sub>r</sub></td><td>radiant temperature</td><td>°C</td></tr><tr><td>P<sub>a</sub></td><td>partial vapour pressure</td><td>P<sub>a</sub></td></tr><tr><td>RH</td><td>relative humidity</td><td>%</td></tr><tr><td>I<sub>cl</sub></td><td>thermal insulation of clothing</td><td>clo</td></tr></tbody></table>
## 2.2 Clothing Insulation Estimation

Indoor temperature is an essential factor in clothing worn inside buildings, as suggested by the non-linear regression relationship between clothing insulation and indoor temperatures based on feedback observations in [7]. In this paper [7], an equation is used to calculate the insulation of clothing based on the temperature inside a room:

$$Cl_{T_{in}} = f(T_a) = 89.279(T_a)^{-1.592} \quad (10)$$

Additionally, the influence of season on clothing insulation is significant in determining thermal comfort, as individuals require different levels of insulation in their clothing to maintain comfort in different seasons [18]. During colder seasons, individuals require more insulation to keep warm, while during warmer seasons, individuals require less insulation to avoid overheating [21]. There are some standard values that correspond to each season, such as the 'Typical summer indoor clothing' equal to 0.5 and the 'Typical winter indoor clothing' equal to 1.0. A year is separated into four periods on the basis of seasonality. Each season is represented by average clothing insulation and the clothing insulation by season $Cl_s$ is depicted in Fig. 1.

<figure><img src="image_2.png" /><figcaption>Fig. 1. $Cl_s$: Average clothing insulation by season values</figcaption></figure>

The insulation value of clothing can also be expressed as representative of a particular clothing ensemble as a function of the outdoor temperature [18].

$$Cl_{T_{out}} = 2.1 \cdot 10^{-5}T_{out}^{3} + 8 \cdot 10^{-4}T_{out}^{2} - 0.0282T_{out} + 0.8167 \quad (11)$$

where $T_{out}$ is the outdoor temperature expressed in °C. This equation is determining the Icl insulation of clothing based on real-time values.

There is a method that combines the $Cl_{T_{in}}$ and the $Cl_s$ to estimate total clothing insulation as [12]:

$$I_{cl} = wCl_s + (1 - w)Cl_{T_{in}}, \quad (12)$$

this method was tested in real-life environments and has proven to be accurate enough. Nonetheless, this method may produce a significant error during periods of time when the outdoor temperature is extremely high or low than the average
seasonal outdoor temperatures (e.g., during spring outdoor temperatures can be more than 25 °C or less than 10 °C). As a result, in this paper, we suggest a new method that will “correct” this specific error by combining all temperature factors. The insulation value of clothing is:

$$I_{cl} = aCl_s + bCl_{T_{in}} + (1 - a - b)Cl_{T_{out}}. \qquad (13)$$

During an experimental phase, the values of a and b were set to a = 0.5 and b = 0.25.

Finally, when estimating the clothing insulation of a person in motion, it is important to account for the dynamic nature of the insulation, which is affected by both the individual’s activity level and the air speed around them. As per the ISO 7730 standard, it is necessary to correct for these factors [2]. Similarly, the ASHRAE 55 Standard provides a correction equation for body movement for activities with a metabolic rate of 1.2 met or higher, expressed as:

$$I_{cl} = I_{cl}(0.6 + \frac{0.4}{met}), \qquad (14)$$

## 2.3 Metabolic Rate

<figure><img src="image_3.png" /><figcaption>Fig. 2. Metabolic rates for typical tasks ASHRAE Standard 55 [17]</figcaption></figure>

The metabolic rate that corresponds to the ideal level of comfort ranges from 84.8 W/m² to 89.9 W/m². This range is determined based on occupants’ feedback who reported feeling thermally comfortable during a survey [16]. According to the ASHRAE metabolic rate table in Fig. 2, this range of activity is classified as low-level activity, similar to activities such as standing or relaxing. It appears that low-level activities are more feasible in a household setting, while more strenuous activities such as cleaning the house are not so frequent [16]. As a result, the metabolic rate could be set to low-level activities (that is, 70–80 W/m² or 1.2 met), such as standing and relaxed activities. During the night, the activity can be set to sleep (that is, 40 W/m² or 0.7 met).
## 2.4 Overall System

The overall conceptual architecture of the system suggested in Sect. 2 is depicted in Fig. 3. The sensors push real-time indoor environmental conditions to be utilized for both thermal comfort and clothing insulation estimation. The date and time are utilized for the estimation of metabolic rate and clothing estimation. In addition, outdoor conditions are utilized for clothing insulation. It may be observed how clothing insulation is estimated through a multifactor process.

<figure><img src="image_4.png" /><figcaption>Fig. 3. Thermal comfort estimation conceptual architecture</figcaption></figure>

# 3 Results

In this section, the results of the suggested methodology of thermal comfort estimation will be presented. Initially, how the experiment was set up will be discussed. Finally, experimental results will be showcased in the last subsection.

## 3.1 Experiment Set up

To test the suggested methodology as described in Sect. 2, data from the CERT.ITI smart-home were exploited [19]. Data from a temperature humidity sensor (i.e., [20]) were retrieved. Values from a specific room for a year (i.e., temperature, humidity) were utilized for the tests (Fig. 4).

The average values during the tested period are presented in Table 2.
<figure><img src="image_5.png" /><figcaption>Fig. 4. Thermal comfort set up experiment</figcaption></figure>

<table><caption>Table 2. Data-set average values</caption><thead><tr><td>Month</td><td>Average Indoor Temperature °C</td><td>Average Indoor Humidity %</td><td>Average Outdoor Temperature °C</td></tr></thead><tbody><tr><td>January</td><td>17.96</td><td>28</td><td>8.5</td></tr><tr><td>February</td><td>18.23</td><td>29.11</td><td>10.7</td></tr><tr><td>March</td><td>20.23</td><td>30.3</td><td>14.3</td></tr><tr><td>April</td><td>22.5</td><td>27.5</td><td>18.6</td></tr><tr><td>May</td><td>22.5</td><td>32.3</td><td>17.3</td></tr><tr><td>June</td><td>25.1</td><td>34.5</td><td>17.9</td></tr><tr><td>July</td><td>26.2</td><td>33.6</td><td>20.1</td></tr><tr><td>August</td><td>28.1</td><td>38.8</td><td>20.8</td></tr><tr><td>September</td><td>27</td><td>23.8</td><td>18.5</td></tr><tr><td>October</td><td>24.3</td><td>20.5</td><td>17.3</td></tr><tr><td>November</td><td>22.96</td><td>27.1</td><td>15.7</td></tr><tr><td>December</td><td>23.4</td><td>28.6</td><td>10.1</td></tr></tbody></table>

## 3.2 Experimental Results

In this section, some indicative results will be presented, one for each season. In the following graphs, real-time values of indoor and outdoor temperatures from the CERTH.ITI smart-home [19], along with the values of the PMV calculated with the two methods will be presented. PMVold value based on Equation (12) and the PMVnew value based on Eq. (13). The graph for each of the four seasons (i.e., Spring, Summer, Fall, and Winter), depicts a representative day from each season. When analyzing the following graphs two facts must be considered. First, the indoor temperature tends to stay more constant than the outdoor temperature and therefore closer to the average seasonal temperature,
for reasons such as insulation of the building, heating and cooling systems, and solar radiation. Second, in terms of the thermal sensation scale in Fig. 5, as the PMV value is closer to 0 there is a better thermal sensation, and it is expected that in the summer season, it tends to go higher, where in the winter season it is lower.

<figure><img src="image_6.png" /><figcaption>Fig. 5. Thermal Sensation Scale</figcaption></figure>

Figure 6 depicts values from April 13th, 2022 during the spring season. Outdoor temperatures may be observed to fall below average spring temperatures (e.g., average below 12 °C). As a result, the line of the PMV value calculated with the new equation is lower than the other, meaning that the resident is feeling colder. While the outdoor and indoor temperature increases, the deviation between the two PMV values is smaller.

<figure><img src="image_7.png" /><figcaption>Fig. 6. Thermal comfort comparison for Spring</figcaption></figure>

Figure 7 shows values from July 13th, 2022 during the summer season. On this day, outdoor temperatures are very close to indoor ones. So, here the fact that the PMV values calculated with the new equation are closer to 0, means that there is a better thermal sensation compared to the old one, as expected.

Following up, in the winter season, Fig. 8 depicts the data values for the day of December 28, 2022. The first thing that becomes apparent is that the line based on the old equation takes values higher than 0, which is in contrast to previous acknowledgments.
<figure><img src="image_8.png" /><figcaption>Fig. 7. Thermal comfort comparison for Summer</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>Fig. 8. Thermal comfort comparison for Winter</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>Fig. 9. Thermal comfort comparison for Fall</figcaption></figure>
Finally, during the fall season, Fig. 9 depicts the values from October 24th, 2022. In this case, outdoor temperature values are far below the indoor ones and this leads to the conclusion that the line base on the new equation is more accurate.

# 4 Conclusions

In general, thermal comfort is an important factor in building occupants' behav-
ior, which is difficult to predict, model, or calculate due to the complexity of
humans. In this paper, a novel approach for calculating thermal comfort was pre-
sented that combines algorithms to enhance the precision of existing approaches.
The proposed methodology was tested in real-time environments and has proven
to be accurate enough.

The analyses of the graphs lead to the conclusion that when the outdoor and
indoor temperatures are close, the new PMV values are more comparable to the
optimum thermal sensation. Furthermore, the old method may produce a sig-
nificant error in the thermal comfort calculation when the outdoor temperature
is extremely higher or lower than the average seasonal outdoor temperature and
consequently the indoor temperatures.

In summary, the proposed approach has shown promising results and could
have significant implications for improving building occupant comfort and energy
efficiency. However, to further improve the accuracy and applicability of the
method, additional steps such as conducting validation tests in different settings,
integrating the methodology with building automation systems, incorporating
other factors that impact thermal comfort, and considering user feedback and
perception should be taken into account.

**Acknowledgements.** This work is partially supported by the PRECEPT project funded by the EU H2020 under Grant Agreement No. 958284.

# References

1. Sansaniwal, S.K., Mathur, J., Mathur, S.: Review of practices for human ther-
mal comfort in buildings: present and future perspectives. Int. J. Ambient Energy
**43**(1), 2097–2123 (2022)

2. ISO, ISO7730: 7730: 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. Management 3(605), e615 (2005)

3. OpenAI: Thermal Comfort Models. OpenAI (2023). openai.com/thermal-comfort-models/

4. Lourenço Niza, I., Broday, E.E.: Development of thermal comfort models over the past years: a systematic literature review. Int. J. Ambient Energy **43**(1), 8830–8846 (2022)

5. ANSI and ASHRAE, Thermal Environmental Conditions for Human Occupancy. Atlanta (2020)
6. Tartarini, F., Schiavon, S.: Pythermalcomfort: a Python package for thermal comfort research. SoftwareX **12**, 100578. Elsevier BV (2020). Crossref, https://doi.org/10.1016/j.softx.2020.100578

7. Dimara, A., Timplalexis, C., Krinidis, S., Tzovaras, D.: A dynamic convergence algorithm for thermal comfort modelling. In: Tzovaras, D., Giakoumis, D., Vincze, M., Argyros, A. (eds.) ICVS 2019. LNCS, vol. 11754, pp. 680–689. Springer, Cham (2019). https://doi.org/10.1007/978-3-030-34995-0_62

8. Harputlugil, T., de Wilde, P.: The interaction between humans and buildings for energy efficiency: a critical review. Energy Res. Soc. Sci. **71**, 101828 (2021)

9. Dimara, A., Krinidis, S., Tzovaras, D.: Comfit: a novel indoor comfort inference tool, 165–170 (2021)

10. Krinidis, S., et al.: Multi-criteria HVAC control optimization. In: 2018 IEEE International Energy Conference (ENERGYCON). IEEE (2018)

11. Peeters, L., et al.: Thermal comfort in residential buildings: comfort values and scales for building energy simulation. Appl. Energy **86**(5), 772–780 (2009)

12. Dimara, A., et al.: NRG4-U: a novel home energy management system for a unique loadprofile. Energy Sources Part A: Recov. Utilization Environ. Effects **44**(1), 353–378 (2022)

13. Oğulata, R.T.: The effect of thermal insulation of clothing on human thermal comfort. Fibres Textiles Eastern Europe **15**(2), 67–72 (2007)

14. Havenith, G., Holmér, I., Parsons, K.: Personal factors in thermal comfort assessment: clothing properties and metabolic heat production. Energy Build. **34**(6), 581–591 (2002)

15. Luo, M., et al.: Human metabolic rate and thermal comfort in buildings: the problem and challenge. Build. Environ. **131**, 44–52 (2018)

16. Dimara, A., et al.: Optimal comfort conditions in residential houses. In: 2020 5th International Conference on Smart and Sustainable Technologies (SpliTech). IEEE (2020)

17. Arendt, K.: Influence of external walls' thermal capacitance on indoor thermal comfort (2013)

18. Schiavon, S., Lee, K.H.: Dynamic predictive clothing insulation models based on outdoor air and indoor operative temperatures. Build. Environ. **59**, 250–260 (2013). https://doi.org/10.1016/j.buildenv.2012.08.024

19. ITI Smart Home CERTH.ITI. https://smarthome.iti.gr/

20. Tuya ZigBee Temperature and Humidity Sensor. https://www.expert4house.com/en/smart-home/sensors-and-detectors/tuya-zigbee-temperature-and-humidity-sensor

21. Kwon, J.Y., Choi, J.: Clothing insulation and temperature, layer and mass of clothing under comfortable environmental conditions. J. Physiol. Anthropol. **32**(1) (2013). https://doi.org/10.1186/1880-6805-32-11