Article

# Experimental Investigation of Adaptive Thermal Comfort in French Healthcare Buildings

Zoubayre El Akili <fnref n="1" />*, Youcef Bouzidi <fnref n="1" />,*, Abdelatif Merabtine <fnref n="2" />, Guillaume Polidori <fnref n="3" /> and Amal Chkeir <fnref n="1" />

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

check for
updates

**Citation:** El Akili, Z.; Bouzidi, Y.; Merabtine, A.; Polidori, G.; Chkeir, A. Experimental Investigation of Adaptive Thermal Comfort in French Healthcare Buildings. *Buildings* **2021**, *11*, 551. https://doi.org/10.3390/buildings11110551

Academic Editor: Cinzia Buratti

Received: 6 October 2021

Accepted: 11 November 2021

Published: 17 November 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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

**Copyright:** © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

**Abstract:** The thermal comfort requirements of disabled people in healthcare buildings are an important research topic that concerns a specific population with medical conditions impacted by the indoor environment. This paper experimentally investigated adaptive thermal comfort in buildings belonging to the Association of Parents of Disabled Children, located in the city of Troyes, France, during the winter season. Thermal comfort was evaluated using subjective measurements and objective physical parameters. The thermal sensations of respondents were determined by questionnaires adapted to their disability. Indoor environmental parameters such as relative humidity, mean radiant temperature, air temperature, and air velocity were measured using a thermal microclimate station during winter in February and March 2020. The main results indicated a strong correlation between operative temperature, predicted mean vote, and adaptive predicted mean vote, with the adaptive temperature estimated at around 21.65 °C. These findings highlighted the need to propose an adaptive thermal comfort strategy. Thus, a new adaptive model of the predicted mean vote was proposed and discussed, with a focus on the relationship between patient sensations and the thermal environment.

**Keywords:** thermal comfort; disabled people; predicted mean vote (PMV); actual mean vote (AMV); operative temperature; adaptive predicted mean vote (aPMV)

## 1. Introduction

Indoor thermal comfort has become an important topic for sustainable building research. Thermal comfort improves the performance of the building, by reducing the energy consumption and greenhouse gas emissions [1]. Sensing, controlling, and predicting thermal comfort is required to ensure healthy indoor conditions for occupants. Providing an appropriate indoor environment, especially in medico-social institutions, is crucial since patients spend 80-90% of the day indoors and are significantly impacted by variations in hygrothermal parameters. Furthermore, such buildings may experience different hydrothermal conditions depending on patient rooms and building design, not to mention the type of disease. Ensuring a comfortable environment can positively contribute to their treatment and care [2,3].

The American Society of Heating, Refrigerating, and Air-Conditioning Engineers (ASHRAE) defined thermal comfort as the state of mind in which satisfaction is expressed with the thermal environment [4]. Over the years, an extensive amount of research has been carried out on thermal comfort but mainly focused on healthy occupants. However, few studies have explored the thermal sensation of occupants in specific buildings such as medico-social institutions due to the lack of knowledge in this area [5,6]. Thermal

<footnote n="1">¹ Interdisciplinary Research on Society-Technology-Environment Interactions (InSyTE), CREIDD, University of Technology of Troyes, 12 Rue Marie Curie, CS 42060, CEDEX, 10004 Troyes, France; youcef.bouzidi@utt.fr (Y.B.); amal.chkeir@utt.fr (A.C.)</footnote>

<footnote n="2">² EPF School of Engineering, 2 Rue Fernand Sastre, 10430 Rosières-Prés-Troyes, France; abdelatif.merabtine@epf.fr</footnote>

<footnote n="3">³ MATIM, Faculty of Exact and Natural Sciences, University of Reims Champagne Ardenne, CEDEX 2, 51687 Reims, France; guillaume.polidori@univ-reims.fr</footnote>

<footnote>* Correspondence: zoubayre.el_akili@utt.fr</footnote>
comfort requirements in medico-social institutions can differ from other types of build-
ings. Nevertheless, the current ASHRAE standards 55-2017 [7] and ISO/TS 14415 [8] lack
information on this topic. The standards outline methods to determine the thermal envi-
ronmental conditions (temperature, humidity, air speed, and radiant effects) for healthy
adults. Nevertheless, these standards are not intended to override any health or critical
process requirements (see ASHRAE Standards: 55-2017, Section 2).

For people with limited adaptive opportunities who are recognized as vulnerable people, an acceptable environment to healthy individuals may be considered unacceptable to them [9]. Assessing thermal comfort in medico-social institutions is therefore challenging. The major challenge is to best meet the thermal comfort needs of different in-situ populations: i.e., medico-social staff and patients.

Among the thermal comfort indices, Fanger's predicted mean vote (PMV) and predicted percentage of dissatisfied (PPD) are the most applicable [10]. The PMV-PPD model requires two personal factors—metabolic rate and clothing insulation—and four environmental factors—air temperature, radiant temperature, air velocity, and humidity [11]. When patients take medications that affect their thermoregulatory system [12], the predictability of comfort ratings may be less reliable. Therefore, thermal adaptation can play a key role in evaluating the thermal sensation of patients.

Recently, the adaptive approach to thermal comfort has interested researchers [13,14], because it takes into account many factors such as climate, culture, as well as psychological and behavioral adaptations, which play an important role in controlling thermal comfort. This approach could be defined as follows: "If a change occurs such as to produce discomfort, people react in ways which tend to restore their comfort" [15]. The cultural background is an important parameter in terms of AMV values (subjective approach) and in terms of PMV values (objective approach). The cultural variation and the availability of personal environmental control options, have an important influence on the PMV model accuracy [16]. Food habit and clothes style, for example, can affect internal thermal energy production and increase the clothing's factor, respectively. On the other hand, several scientific have highlighted the importance of including the culture influence on thermal comfort, this topic would be important in architect field, by improving the design possibilities [17]. Pereira et al. [18] showed that detection of occupant behavior in buildings is an important parameter, which leading to better indoor environment and thermal comfort. Pereira et al. [19] showed that occupant behavior in buildings is decisive for their hygrothermal performance. Verheyen et al. [20] investigated the thermal comfort of patients in a Belgian healthcare facility by comparing objective parameters and subjective measures of thermal comfort for different patient groups. They concluded that PMV can adequately predict mean thermal sensation for the majority of patients. On the contrary, Sattayakorn et al. [21] indicated that the PMV-PPD model was unsuitable for evaluating the thermal comfort of occupants in healthcare facilities in tropical regions. This result was later confirmed by Yau et al. [22], who indicated that the PMV-PPD model may not be suitable for tropical hospitals as occupants may be more satisfied with a hot indoor environment. In the study of Hashiguchi et al. [23], a comparison of the thermal comfort of patients and medical staff concluded that most patients were comfortable, while the medical staff were uncomfortable. However, this study did not compare subjective responses and objective measures with PMV predictions, which is of crucial importance in this field.

Skoog et al. [24] showed little difference in the thermal sensations between staff and patients in Swedish hospitals. In another study by Hill et al. [25] on patients with physical disabilities, their most frequent request was to be warmer, whereas staff generally wanted to be cooler. In the same context, Khodakarami et al. [26] highlighted that hospital occupants have different thermal comfort requirements that are difficult to accommodate in the same space, and thus adaptive comfort is recommended for each group. Regarding the management of the indoor environment, Kim et al. [27] stated that indoor hygrothermal conditions must be carefully managed in healthcare facilities to improve staff satisfaction with their indoor environment. They concluded that the green features implemented
by leadership in energy and environmental design (LEED) in healthcare facilities had a good effect on healthcare staff's level of comfort. Pereira et al. [28] highlighted that energy efficient and smarter buildings must include the requirements of the occupants, also discussed the importance of occupant behaviour motivations for the design and implementation of smart systems in buildings. Walker et al. [29] discussed the importance of thermal comfort management in care homes, and the need to distinguish them from other seemingly similar buildings, because the thermal comfort requirements of patients can be associated with health risks.

To manage the indoor environment, Nomura et al. [30] noted that to bring about positive health outcomes for patients, there must be a minimum of six air changes per hour (ACH), although in spaces with heating, ventilation, and air conditioning (HVAC) systems, this rate may be reduced to 4 ACH. Hwang et al. [31] reported that patients are more comfortable in warmer and humid indoor environments in Taiwanese hospitals and are insensitive to indoor thermal changes. The study by Kameel et al. [32] in healthcare settings showed that hygrothermal parameters can activate or deactivate viruses and that low humidity levels can increase susceptibility to respiratory disease. Smith et al. [33] showed that patients usually preferred air temperatures between 21.5 °C and 22 °C and relative humidity between 30% and 70%. Bouzidi et al. [34] investigated adaptive thermal comfort during summer in French healthcare buildings, results shows that the adaptive temperature was 25.0 °C with upper and lower limits of 24.7 °C and 25.4 °C. Pereira et al. [35] reported that the design of the buildings can be improved by consideration of the effects that the spatial and human characteristics have on the indoor environment quality.

Despite the importance of these studies on thermal comfort, the well-being of occupants in specific situations such as medico-social institutions requires further exploration. First, in the case of healthcare buildings, the indoor temperature set-point is usually determined by healthcare staff, meaning that it is not correlated to patients' thermal sensation. Second, the adaptive approach is essential to ensure thermal comfort in healthcare buildings, although no study conducted in a French climate has determined the adaptive temperatures for patients. Finally, a new thermal comfort index is necessary to correctly assess the thermal conditions of healthcare buildings.

Hence, the main purpose of this field study is to determine the thermal sensation of occupants in order to adapt thermal approaches to vulnerable people in medico-social buildings. Thermal comfort can take on multiple meanings, being associated with vulnerability, the indoor environment, and the provision of effective care. Thus, adapting the thermal environment of occupants in medico-social institutions is crucial to maintain a good quality of service and support a population that is particularly vulnerable to illness.

To achieve this goal, we correlate thermal sensation to the indoor environment. Indeed, the HVAC system can be set based on more appropriate metering in the future to reduce energy consumption while improving the indoor climate for patients. Therefore, this study investigates the effect of indoor thermal conditions on the satisfaction of occupants in healthcare facilities. This investigation consists of two parts. The first evaluates the indoor environmental parameters of the studied buildings in terms of air temperature, operative temperature, air velocity, and humidity, as well as the activity and clothing insulation of their occupants. The second part explores the perception of patients regarding their satisfaction with the indoor thermal environment in terms of the actual mean vote (AMV). The optimal temperature is calculated based on the adaptive thermal comfort approach to develop a new PMV model that provides a relationship between the indoor parameters and the thermal sensation of occupants.

## 2. Research Methodology

The survey is based on a mixed approach consisting of both objective and subjective evaluations. These complementary methods were simultaneously implemented in the same locations. Combining the indoor environment with physiological and psychological parameters (thermal perception) provides a clearer understanding of the thermal interac-
tions between the buildings and their occupants, and a more complete understanding of how these elements can influence the overall satisfaction of occupants (Figure 1). Indeed, the objective part of the study deals with PMV and PPD indices to assess the indoor thermal environment, whereas the subjective investigations are based on a questionnaire and a ruler with pictorial representations to evaluate the actual thermal perception of occupants. We then combined these methods, first to adapt the indoor environment of patients and then to propose a better index to assess their thermal comfort.

<figure><img src="image_3.png" /><figcaption>**Figure 1.** Proposed research methodology.</figcaption></figure>

In this paper, thermal comfort parameters were experimentally investigated during the winter season in the “Gai Soleil” medico-social institution in the city of Troyes, located in eastern France (latitude 48.32°, longitude 4.08°). These medico-social buildings accommodate people with mental disabilities. They were constructed in 1963 and are managed by the “A.P.E.I. of Aube.”. The building envelope is a non-insulating brick wall’s structure (except for the roofs are insulated with mineral wool insulation), the thermal conductivity of the brick is 1.6 W·m⁻¹·K⁻¹, and the density is 2300 kg·m⁻³. Generally, standard insulation methods are not applicable. The normal working hours for this institute are weekdays, from 8 am to 6 pm. The survey was conducted in six buildings or groups of buildings as highlighted by the red color in Figure 2b from 3 February to 13 March 2020. Table 1 shows the distribution of patients in the different buildings.
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Surveyed buildings: (a) outside view of the medico-social buildings and (b) plan of the buildings.</figcaption></figure>

<table><caption>**Table 1.** Distribution of patients and staff in the buildings.</caption><thead><tr><td>Building Group</td><td>Number of Patients</td><td>Number of Staff</td></tr></thead><tbody><tr><td>1</td><td>8</td><td>2</td></tr><tr><td>2</td><td>10</td><td>2</td></tr><tr><td>3</td><td>7</td><td>2</td></tr><tr><td>4</td><td>5</td><td>1</td></tr><tr><td>5</td><td>10</td><td>2</td></tr><tr><td>6</td><td>7</td><td>2</td></tr><tr><td>Total</td><td>47</td><td>11</td></tr></tbody></table>
## 2.1. *Subjective Method*

A longitudinal thermal comfort survey (14 days over 2 months) was carried out at the “Gai Soleil” from 3 February to 13 March 2020 (Figure 3). Patients and healthcare staff completed the surveys in the shared spaces within each building, leading to a total of 423 valid questionnaires for patients (47 subjects) and 62 for staff (11 subjects).

<figure><img src="image_5.png" /><figcaption>**Figure 3.** Survey picture.</figcaption></figure>

The subjective approach consists of a survey with questions and a ruler with pictorial representations, which was designed in collaboration with a psychologist. As shown in Figure 4, the ruler is a subjective measuring tool based on the ISO seven-point thermal sensation scale, which has the shape of a large thermometer with a pictorial representation. This adaptation of the thermal comfort scales to disabled people resulted in a good response rate [36]. To make the questionnaire easier and quickly obtain the answers, patients were questioned in two phases. First, they simply indicated their thermal sensation (Hot, Neutral, Cold) and then depending on the response, they were asked to use the pictorial representations to specify this thermal sensation as cold, cool, or slightly cool, or hot, warm, or slightly warm. This field study was conducted simultaneously with patients and healthcare staff in the same indoor environment. The aim was to understand their thermal feelings and preferences under given thermal environments. The inclusion of staff allowed us to identify possible differences between the two groups under similar conditions. The procedure took about 5-8 min per participant. During the measurements, the interior doors were left open in most of the patient areas, whereas the windows were often closed. To identify which patients would be invited to participate in our thermal sensation surveys, we selected, in collaboration with the healthcare staff, patients who were able to understand the questions and articulate their answers, either orally or by gestures. The results thus obtained indicated that only 55.3% of patients were able to respond (Table 2).

## 2.2. *Objective Method*

In this study, the physical parameters were continuously measured according to the standard ISO 7726 for seated persons [37]. The indoor environmental parameters were measured using the microclimate station HD32.3 produced by DeltaOHM (Figure 5). The Datalogger records the physical parameters included in the PMV and PPD indices such as the relative humidity (RH(%)), mean radiant temperature, air temperature, and air velocity (V(m/s)). For low air velocity and small temperature variances, it is possible to evaluate the operative temperature as the average between the indoor air temperature and the mean radiant temperature.
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Subjective measuring tool.</figcaption></figure>

<table><caption>**Table 2.** Summary of patients' ability to answer the questions.</caption><thead><tr><td>Gender</td><td>Number of Patients</td><td>Able to Respond</td><td>Percentage (%)</td></tr></thead><tbody><tr><td>Female</td><td>28</td><td>17</td><td>60.61</td></tr><tr><td>Male</td><td>57</td><td>30</td><td>52.63</td></tr><tr><td>Total</td><td>85</td><td>47</td><td>55.29</td></tr></tbody></table>

<figure><img src="image_7.png" /><figcaption>**Figure 5.** Thermal microclimate station.</figcaption></figure>

The operative temperature defined as "the uniform temperature of an imaginary black enclosure in which an occupant would exchange the same amount of heat by radiation and convection as in the actual non uniform environment" [38].
$T_{op}$, PMV, and PPD indices are thus defined as follows:

$$T_{op} = \frac{T_a + T_{rm}}{2} \qquad (1)$$

$$PMV = (0.303 \cdot e^{(-0.036 \cdot M)} + 0.028) \cdot \{(M - W) - 3.05 \cdot 10^{-3} \cdot [5733 - 6.99 \cdot (M - W) - p_a] - 0.42 \cdot [(M - W) - 58.15] - 1.7 \cdot 10^{-5} \cdot M \cdot (5867 - p_a) - 0.0014 \cdot M \cdot (34 - T_a) - 3.96 \cdot 10^{-8} \cdot f_{cl} \cdot [(T_{cl} + 273)^4 - (T_{rm} + 273)^4] - f_{cl} \cdot h_c \cdot (T_{cl} - T_a)\} \qquad (2)$$

$$PPD = 1 - 0.95 \cdot e^{(-0.03353 \cdot PMV^4 - 0.2179 \cdot PMV^2)} \qquad (3)$$

where M (W/m²), W (W/m²), and p_a (Pa) are the metabolic rate, external work, and partial vapor pressure, respectively; f-cl is the ratio of the clothed body surface to the naked body surface; h_c (W·m⁻²·K⁻¹) is the convective heat transfer coefficient; T_a (°C) is the air temperature; T_cl (°C) is the clothing surface temperature; T_rms (°C) is the mean radiant temperature; and T_op (°C) is the operative temperature.

Table 3 depicts the measured physical parameters and their accuracy. The metabolic rate and clothing insulation were estimated based on ISO 7730 [39]. In this study, the metabolic rate was set at 1.2 met, corresponding to sedentary activities, while the mean clothing insulation value was calculated as 1.1 clo. Table 4 summarizes the physiological parameters of the respondents, and Table 5 summarizes the environmental parameters.

<table><caption>**Table 3.** Physical parameters and instrument accuracy.</caption><thead><tr><th rowspan="2">Parameter</th><th colspan="2">Thermal Microclimate Station</th></tr><tr><th>Accuracy</th><th>Valid Range</th></tr></thead><tbody><tr><td>Air Velocity (m/s)</td><td>± 0.2<br/>± 0.3</td><td>0-1<br/>1-5</td></tr><tr><td>Relative Humidity (%)</td><td>± 1.5<br/>± 2.0</td><td>0-90<br/>90-100</td></tr><tr><td>Temperature (°C)</td><td>Class 1/3 DIN</td><td>-40 to +100</td></tr><tr><td>Globe Temperature (°C)</td><td>Class 1/3 DIN</td><td>-10 to +100</td></tr></tbody></table>

<table><caption>**Table 4.** Summary of the physiological parameters of patients.</caption><thead><tr><th rowspan="2">Parameter</th><th colspan="4">Male</th><th colspan="4">Female</th></tr><tr><th>Max</th><th>Min</th><th>Mean</th><th>SD</th><th>Max</th><th>Min</th><th>Mean</th><th>SD</th></tr></thead><tbody><tr><td>Age (years)</td><td>20</td><td>10</td><td>14.03</td><td>3.04</td><td>19</td><td>8</td><td>14.76</td><td>2.74</td></tr><tr><td>Weight (kg)</td><td>110</td><td>20</td><td>45.66</td><td>17.88</td><td>122</td><td>19</td><td>51</td><td>23.39</td></tr><tr><td>Height (m)</td><td>1.8</td><td>1.22</td><td>1.53</td><td>0.15</td><td>1.82</td><td>1.1</td><td>1.48</td><td>0.16</td></tr><tr><td>BMR (kcal)</td><td>2312.3</td><td>872.2</td><td>1360.6</td><td>292.5</td><td>2080</td><td>1004.4</td><td>1350.2</td><td>238.6</td></tr><tr><td>BMI (kg/m²)</td><td>38.97</td><td>13.43</td><td>19.10</td><td>4.97</td><td>35.92</td><td>14.75</td><td>22.37</td><td>6.53</td></tr></tbody></table>

<table><caption>**Table 5.** Summary of the environmental parameters.</caption><thead><tr><th>Parameter</th><th>T<sub>a</sub> (°C)</th><th>T<sub>op</sub> (°C)</th><th>T<sub>rm</sub> (°C)</th><th>RH (%)</th><th>V (m/s)</th></tr></thead><tbody><tr><td>Min</td><td>19.8</td><td>19.55</td><td>19.3</td><td>25.4</td><td>0.00</td></tr><tr><td>Max</td><td>26.9</td><td>26.9</td><td>26.9</td><td>61.5</td><td>0.11</td></tr><tr><td>Mean</td><td>23.31</td><td>23.35</td><td>23.38</td><td>38.47</td><td>0.03</td></tr><tr><td>SD</td><td>1.13</td><td>1.17</td><td>1.24</td><td>6.49</td><td>0.01</td></tr></tbody></table>
**3. Results and Discussion**

All variables were screened to ensure that there was sufficient variation to perform regression analysis to evaluate the thermal comfort conditions and obtain a patient-adapted PMV model. The results were processed, and correlations between thermal comfort indices, mean radiant temperatures, and operative temperatures were identified.

**3.1. *Longitudinal Thermal Comfort Survey***

Figure 6 depicts the results regarding the thermal sensation of patients and healthcare staff. Concerning patients, 62.16% of their thermal sensation votes ranged from -1 to -3 (slightly cool, cool, cold), while 23.16% perceived the thermal sensation as neutral. In general, the negative values obtained from the survey indicate a cooler thermal sensation from the point of view of patients. However, 85.5% of the thermal sensation votes of healthcare staff ranged from 0 to 1 (neutral, slightly warm). This indicates that most staff remained satisfied with the indoor thermal environment. This difference is essentially due to the adaptation capacities of respondents. Staff satisfaction can be explained by their ability to adapt to the indoor environment by means of their clothing, drinking, eating, activity level, and opening or closing windows. On the other hand, the adaptive opportunities of patients may be restricted by their disability and health conditions, which also differ from one patient to another. Therefore, an adaptive indoor environment needs to be created for patients.

<figure><img src="image_8.png" /><figcaption>**Figure 6.** Subjective thermal sensation of patients and healthcare staff.</figcaption></figure>

**3.2. *Relationship between PMV, PPD, and Operative Temperature***

Most indoor comfort studies consider the relationship between PMV and operative temperature. This relationship was also successfully established in the present work and is presented in Equation (4) where the high coefficient $R^2 = 0.96$ indicates the strong dependence of PMV on operative temperature (Figure 7). Our findings show that the PMV decreased at a low operative temperature (i.e., >22 °C) and increased at a high operative temperature (i.e., <22 °C).
<figure><img src="image_9.png" /><figcaption>Figure 7. Predicted mean vote versus operative temperature.</figcaption></figure>

Based on Equation (5) and as shown in Figure 8, when the indoor operative temperature is higher than 22.7 °C or lower than 20.8 °C, the PPD is outside the neutral thermal comfort zone (PPD < 6%: the range with the smallest percentage of environmental dissatisfaction) recommended for spaces occupied by fragile people with special needs [40]. To ensure a higher thermal comfort level, the optimal temperature is calculated based on the adaptive thermal comfort model.

$$PMV = 0.22 \cdot T_{op} - 4.82 \qquad (4)$$

$$PPD = 1.02 \cdot T_{op}^2 - 44.24 \cdot T_{op} + 485.34 \qquad (5)$$

<figure><img src="image_10.png" /><figcaption>Figure 8. Predicted percentage of dissatisfaction versus operative temperature.</figcaption></figure>
### 3.3. Comparison between PMV and AMV

Figure 9 illustrates the comparison between PMV and AMV values based on the measurements for men and women versus the operative temperature ranging from 19.55 °C to 26.9 °C. Female and male patients had similar results in the thermal sensation vote. The mean absolute difference in AMV during the winter season is equal to 0.14. These results also reveal that the fitted regression line for subjects' AMV is below the PMV linear curve (Figure 9). Therefore, both men and women experience the indoor environment as colder than the measurement results according to Fanger's model (note the absence of a significant physiological variance between men and women at the 95% confidence level with *p* < 0.05; see Table 6). This discrepancy may be explained by the patients' limited ability to adapt to the indoor environment (adaptive opportunities), which is not considered in Fanger's PMV model. Furthermore, the relationship significance between PMV and AMV values is significantly low (R² = 0.04; Figure 10), which accords with previous studies conducted in other types of buildings [41,42].

<table><caption>**Table 6.** Analysis of variance for physiological parameters.</caption><thead><tr><th>Source of Variation</th><th>SS</th><th>DF</th><th>MS</th><th>F-Value</th><th>p-Value</th><th>F-Crit</th></tr></thead><tbody><tr><td colspan="7">Age</td></tr><tr><td>Between groups</td><td>5.804</td><td>1.00</td><td>5.804</td><td>0.586</td><td>0.448</td><td>4.057</td></tr><tr><td>Within groups</td><td>446.025</td><td>45.00</td><td>9.912</td><td></td><td></td><td></td></tr><tr><td>Total</td><td>451.830</td><td>46.00</td><td></td><td></td><td></td><td></td></tr><tr><td colspan="7">Weight</td></tr><tr><td>Between groups</td><td>308.652</td><td>1.00</td><td>308.652</td><td>0.770</td><td>0.385</td><td>4.057</td></tr><tr><td>Within groups</td><td>18,034.667</td><td>45.00</td><td>400.770</td><td></td><td></td><td></td></tr><tr><td>Total</td><td>18,343.319</td><td>46.00</td><td></td><td></td><td></td><td></td></tr><tr><td colspan="7">Height</td></tr><tr><td>Between groups</td><td>0.019</td><td>1.00</td><td>0.019</td><td>0.749</td><td>0.391</td><td>4.057</td></tr><tr><td>Within groups</td><td>1.153</td><td>45.00</td><td>0.026</td><td></td><td></td><td></td></tr><tr><td>Total</td><td>1.172</td><td>46.00</td><td></td><td></td><td></td><td></td></tr><tr><td colspan="7">BMR</td></tr><tr><td>Between groups</td><td>1173.945</td><td>1.00</td><td>1173.945</td><td>0.016</td><td>0.901</td><td>4.057</td></tr><tr><td>Within groups</td><td>3,392,461.974</td><td>45.00</td><td>75,388.04</td><td></td><td></td><td></td></tr><tr><td>Total</td><td>3,393,635.919</td><td>46.00</td><td></td><td></td><td></td><td></td></tr><tr><td colspan="7">BMI</td></tr><tr><td>Between groups</td><td>115.644</td><td>1.00</td><td>115.644</td><td>3.72</td><td>0.06</td><td>4.057</td></tr><tr><td>Within groups</td><td>1399.057</td><td>45.00</td><td>31.090</td><td></td><td></td><td></td></tr><tr><td>Total</td><td>1514.701</td><td>46.00</td><td></td><td></td><td></td><td></td></tr></tbody></table>

where: BMR: basal metabolic rate (kcal); BMI: body mass index (kg·m⁻²) and F-Crit: critical value.

### 3.4. Adaptive Thermal Comfort and Patients

The adaptive approach combines the subjective and objective approaches to create a comfortable thermal environment. In the literature, most adaptive thermal comfort studies are carried out in educational buildings as opposed to medical buildings [43,44]. A thermal model known as the adaptive predicted mean vote (aPMV) model, which considers the aforementioned factors and draws on the “black box” theory, was introduced by Yao et al. [45]. This model is more suitable to describing the indoor thermal environment in the buildings considered in this study. Figure 11 shows the flowchart of the holistic principle underlying the adaptive model. The interaction between the thermal environment and the occupants is complex and depends on many factors (thermoregulatory mechanisms, medical treatment, activity level, clothing insulation, etc.).
<figure><img src="image_11.png" /></figure>

<figure><img src="image_12.png" /><figcaption>**Figure 9.** Comparison between actual mean vote (AMV) and predicted mean vote (PMV) versus operative temperature for (a) male and (b) female patients.</figcaption></figure>
<figure><img src="image_13.png" /><figcaption>Figure 10. Predicted mean vote (PMV) versus actual mean vote (AMV).</figcaption></figure>

<figure><img src="image_14.png" /><figcaption>Figure 11. Flowchart of the holistic principle underlying the adaptive model.</figcaption></figure>

3.5. *Adaptive Opportunities*

Regarding the use of certain adaptive actions to create a comfortable thermal envi-
ronment, adaptive opportunities presented in the questionnaire focused on the patients'
daily habits. As shown in Figure 12, almost 49% of patients indicated that drinking a hot
beverage could increase their thermal sensation, an activity probably related to the winter
season of the study, when patients want to be warmer. Nevertheless, 28.37% reported doing
nothing in terms of thermal adaptation; these results are in agreement with the neutral
thermal sensations expressed by patients (Figure 6). Thus, an adaptive model in which the
indoor environment can be adapted to all patients should be developed. The proportion of
other adaptive responses used by patients to improve their thermal sensation ranged from
1.65% to 8.75%.
<figure><img src="image_15.png" /><figcaption>Figure 12. Proportion of patients' adaptive responses.</figcaption></figure>

3.6. *Adaptive Indoor Temperature*

Due to the discrepancy between PMV and AMV, the optimal temperature needs to be calculated based on the adaptive thermal comfort. Using the black box theory (a system viewed in terms of inputs and outputs without knowledge of the internal procedure, in the black box theory, the adaptive predicted mean vote aPMV uses PMV index as the input; Figure 13), aPMV can be described as follows:

$$aPMV = T \cdot E - aPMV \cdot K \cdot T \qquad (6)$$

where:

$$PMV = T \cdot E \qquad (7)$$

So:

$$aPMV = \frac{PMV}{\left(1 + \frac{K \cdot PMV}{E}\right)} \qquad (8)$$

We set:

$$\alpha = \frac{K}{E} \qquad (9)$$

Finally:

$$aPMV = \frac{PMV}{(1 + \alpha \cdot PMV)} \qquad (10)$$

where T is the transfer function (thermoregulatory system), K the feedback system (psychological and behavioral impact coefficient), and E the physical stimuli.

Here the adaptive coefficient α, representing the patient's capacity to adapt to the environment, considering parameters such as culture, climate, social, psychological, and behavioral adaptations, which have an impact on the thermal perception. The coefficient α is based on the findings of surveys of thermal comfort conducted in the field, it is therefore a coefficient linked to the type of population. It is, therefore, necessary to carry out several surveys studies to meet the different thermal comfort conditions required by different occupants in different places, e.g., in schools, hospitals..., α was calculated using
the least square method to adjust the field data sets. α can be described by the following equation [46]:

$$ \alpha = \frac{\sum_{i}^{n} y_{i} - x_{i}}{i} \qquad (11) $$

where:

$$ x = \frac{1}{\text{AMV}} \qquad y = \frac{1}{\text{PMV}} \qquad (12) $$

i: number of data.

In cooler conditions (PMV < 0), there are eight sets of data in which the value of the adaptive coefficient is calculated as follows:

$$ \alpha = \frac{\sum_{i}^{8} y_{i} - x_{i}}{8} = -7.31 \qquad (13) $$

So:

$$ aPMV = \frac{PMV}{(1 - 7.31 \cdot PMV)} \qquad (14) $$

In warmer conditions (PMV > 0), there are 14 sets of data in which the value of the adaptive coefficient is calculated as follows:

$$ \alpha = \frac{\sum_{i}^{14} y_{i} - x_{i}}{14} = 5.24 \qquad (15) $$

So:

$$ aPMV = \frac{PMV}{(1 + 5.24 \cdot PMV)} \qquad (16) $$

<figure><img src="image_16.png" /><figcaption>Figure 13. Thermal comfort adaptive model mechanism modified from [47].</figcaption></figure>

Figure 14 shows that the adaptive predicted mean vote (aPMV) varies from -0.10 to 0.16, which is within the neutral thermal comfort zone [-0.2, +0.2] recommended for spaces occupied by very sensitive and fragile people. aPMV significantly reduced the sensation of discomfort compared to the PMV model, which varies from -0.47 to 1.11. As we can see in Equations (14) and (16), the advantage of the adaptive model is that complex adaptation is represented as a single value.

Figure 15 shows the polynomial correlation equation between the calculated adaptive predicted vote and the operative temperature. The adaptive temperature was calculated so that PMV = aPMV (i.e., patient-adaptive environment). Since the objective is to ensure that patients feel warm, we choose the high temperature top = 21.7 °C based on Equation (17) (this result is limited to the temperature range in which it was carried out):

$$ -0.012 \cdot T_{\text{op}}^{2} + 0.603 \cdot T_{\text{op}} - 7.431 = 0 \qquad (17) $$
<figure><img src="image_17.png" /><figcaption>**Figure 14.** Comparison between adaptive predicted mean vote (aPMV) and predicted mean vote (PMV) versus actual mean vote (winter).</figcaption></figure>

<figure><img src="image_18.png" /><figcaption>**Figure 15.** Adaptive predicted mean vote (aPMV) versus operative temperature.</figcaption></figure>

## 4. Correcting the PMV Model for the Patients

A new design of the PMV model takes into account within-group and between-group differences. For each response, we simultaneously measured the indoor environmental parameters. Therefore, a first-order correction to the PMV model for the patients was possible. Our data are longitudinal, while our data points might not be truly independent. We have six groups (Figure 2) as well as different observations per group, while our survey study may be insufficient if we try to fit models with too many parameters. Therefore,
a linear mixed-effect model was chosen to incorporate all the data, even in the case of many covariates.

Linear mixed-effect models are statistical models containing both fixed and random effects. These models are useful in a wide variety of disciplines but rarely used in the thermal comfort field [38]. They are particularly useful in settings with repeated measurements. As many types of mixed-error models exist, we compared random intercept, random slope, intercept, and random slope models. Figure 16 shows that the intercept changes from one group to another, while the slope remains mostly stable. To select the best model, the Bayesian information criterion (BIC) was used, and the model with the lowest BIC was preferred. Finally, a random intercept model was selected. BIC is defined as [48]:

$$BIC = k \cdot \ln(n) - 2 \cdot \ln(\hat{L}) \qquad (18)$$

where $\hat{L}$ is the maximized value of the likelihood function of the model, n the sample size, and k the number of parameters estimated by the model.

<figure><img src="image_19.png" /><figcaption>Figure 16. Predicted mean vote (PMV) versus operative temperature in each building.</figcaption></figure>

For i = (1... n) and j = (1... m), linear mixed models were described as follows:

$$Y_{ij} = \beta_{0i} + \beta_{1i}X_{ij} + \varepsilon_{ij} \qquad (19)$$

$$\beta_{0i} = \alpha_{0} + \mu_{0i}, \beta_{1i} = \alpha_{1} + \mu_{1i} \qquad (20)$$
$$Y_{ij} = \alpha_0 + \mu_{0i} + \alpha_1 X_{ij} + \mu_{1i} X_{ij} + \varepsilon_{ij} \quad (21)$$

$$Y_{ij} = (\alpha_0 + \alpha_1 X_{ij}) + (\mu_{0i} + \mu_{1i} C_{ij} + \varepsilon_{ij}) \quad (22)$$

where $(\alpha_0 + \alpha_1 C_{ij})$ are the fixed effects, $(\mu_{0i} + \mu_{1i} C_{ij} + \varepsilon_{ij})$ the random effects, n the number of groups, m the number of repetitions per group, X the independent variable, $\mu_{0i}$ the random intercept associated with each group, $\mu_{1i}$ the random slope associated with each group, and $\varepsilon_{ij}$ the error term.

In this study, a random intercept model was selected:

$$Y_{ij} = \alpha_0 + \alpha_1 X_{ij} + \mu_{0i} + \varepsilon_{ij} \quad (23)$$

We looked for a model in which the random intercept was the same for all groups. Equation (23) thus becomes:

$$\tilde{Y}_i = \alpha_0 + \alpha_1 X_i + \mu_0 \quad (24)$$

Using R programming language, we propose an AMV-corrected PMV model in which the independent variable is the operative temperature, and the output is the predicted mean vote for patients (PMVp). Figure 17 compares the evaluations of PMVp and PMV as a function of the operative temperature. The variation of PMV and PMVp in both cases is homoscedastic, so the linearity assumption is valid. The validity of the model can be verified using the coefficient of determination ($R^2_{PMVP} = 0.99$), while the residual analysis is defined as the difference between the actual observation and the corresponding fitted value [49]. The residuals versus predicted values were randomly distributed around zero (Figure 18), with the normal probability plot of the residuals resembling a straight line (Figure 19); therefore, the validity of the model is confirmed.

$$\text{PMV}_p = -4.80952 + 0.221020 \cdot \text{T}_{\text{op}} + \mu_0 \sim N(0, \sigma_0^2 = 0.0001881) \quad (25)$$

where $\sigma$ is the standard deviation and $\mu_0$ the random intercept over the entire data set.

<figure><img src="image_20.png" /><figcaption>Figure 17. Predicted mean vote for patients (PMVp) and predicted mean vote (PMV) as a function of operative temperature.</figcaption></figure>
<figure><img src="image_21.png" /><figcaption>Figure 18. Residuals versus fitted values in the predicted mean vote for patients (PMVp).</figcaption></figure>

<figure><img src="image_22.png" /><figcaption>Figure 19. Normal probability plot of residuals.</figcaption></figure>

**5. Conclusions, Limitations, and Perspectives**

The research on thermal comfort described in this paper was carried out in a medico-social institution in the French city of Troyes in the winter season. The main outcome is that great care should be taken when interpreting the results of thermal comfort studies for vulnerable populations, including people with disabilities. The most important conclusions are as follows:

i. The comparison between patients and staff showed that thermal comfort is strongly correlated with subjective thermal perception, which is influenced by health conditions (including disease type and treatment). Regarding the thermal sensation of patients and healthcare staff. A total of 62.16% of patients thermal sensation votes
ranged from -1 to -3 (slightly cool, cool, cold), while 23.16% perceived the thermal sensation as neutral. However, 85.5% of the thermal sensation votes of healthcare staff ranged from 0 to 1 (neutral, slightly warm). This difference is essentially due to the adaptation capacities of respondents.

ii. The relationship between PMV and PPD indices and operative temperature was successfully established. However, in our case study, the PMV always overestimated the thermal sensation of patients.

iii. Concerning adaptive opportunities, 49.00% of patients indicated that drinking a hot beverage could increase their thermal sensation. Nevertheless, 28.37% reported doing nothing in terms of thermal adaptation. The proportion of other adaptive responses used by patients to improve their thermal sensation ranged from 1.65% to 8.75%. These results probably related to the effect of the winter season of the study.

iv. Due to the discrepancy between PMV and AMV, the optimal temperature needs to be calculated based on the adaptive thermal comfort. In the studied buildings, the adaptive temperature is calculated, using the black box theory and is around 21.65 °C.

v. We proposed a new patient-predicted mean vote (PMVp) that explores the relationship between patients' sensations and their thermal environment. This model takes into account within-group and between-group differences.

These results address the current concerns about assessing the thermal comfort of vulnerable individuals. In terms of study limitations, the reduced number of respondents limited the survey results, even though the percentage of participants remains acceptable for this vulnerable population (see the study of Del Ferraro et al., with only 30 subjects [50]). Given the lack of data on CO₂ concentration, CO₂ measurements analysis is also considered as a limitation of this study.

Our findings may be considered a first attempt at improving the current standards for vulnerable populations regarding their thermal comfort requirements. Further studies should be conducted to evaluate the thermal comfort conditions in healthcare buildings. In the future, we intend to link the thermo-physiological state of this vulnerable population to their adaptive thermal comfort.

**Author Contributions:** Conceptualization, Z.E.A., Y.B., A.M., G.P. and A.C.; methodology, Z.E.A., Y.B., A.M., G.P. and A.C.; software and materials, Z.E.A., Y.B. and A.M.; investigation, Z.E.A., Y.B., A.M., A.C. and G.P.; writing—original draft preparation, Z.E.A.; writing—review and editing, Z.E.A., Y.B., A.M. and G.P.; supervision and funding acquisition, Y.B., A.M. and G.P. All authors have read and agreed to the published version of the manuscript.

**Funding:** This work was supported by European Regional Development Fund (ERDF) for the Doctoral research programme "Retcli" agreement number: CA0023755.

**Institutional Review Board Statement:** This study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Ethics committee of APEI Troyes (February 2021).

**Informed Consent Statement:** Informed consent was obtained from all subjects involved in the study and in their absence by the authorized legal representative.

**Data Availability Statement:** The data presented in this study are available on request from the corresponding author.

**Acknowledgments:** This work was conducted in collaboration with the Association of Parents of Inadapted Children (APEI) and was supported by the Regional Council of Champagne-Ardenne (CRCA) and the European Fund for Economic and Regional Development (FEDER).

**Conflicts of Interest:** The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.
Abbreviations

<table><tr><td>A.P.E.I</td><td>Association of Parents of Inadapted Children</td></tr><tr><td>AMV</td><td>Actual Mean Vote</td></tr><tr><td>PMV</td><td>Predicted Mean Vote</td></tr><tr><td>aPMV</td><td>Adaptive Predicted Mean Vote</td></tr><tr><td>BIC</td><td>Bayesian Information Criterion</td></tr><tr><td>PPD</td><td>Percentage People Dissatisfied</td></tr><tr><td>HVAC</td><td>Heating, Ventilation, and Air Conditioning systems</td></tr><tr><td>DF</td><td>Degrees of Freedom</td></tr><tr><td>SS</td><td>Sum of Squares</td></tr><tr><td>MS</td><td>Mean Squares</td></tr></table>

References

1. Rocha, A.; Pinto, D.; Ramos, N.M.; Almeida, R.M.; Barreira, E.; Simões, M.L.; Martins, J.P.; Pereira, P.F.; Sanhudo, L. A case study to improve the winter thermal comfort of an existing bus station. *J. Build. Eng.* **2020**, *29*, 101123. [CrossRef]

2. Azizpour, F.; Moghimi, S.; Lim, C.; Mat, S.; Zaharim, A.; Sopian, K. Thermal comfort assessment in large scale hospital: Case study in Malaysia. In Proceedings of the 4th WSEAS International Conference on Energy and Development-Environment-Biomedicine, Corfu Island, Greece, July 14–16 2011; Mastorakis, N., Mladenov, V., Bojkovic, Z., Topalis, F., Psarris, K., Eds.; World Scientific and Engineering Academy and Society (WSEAS): Stevens Point, WI, USA, 2011; pp. 171–174.

3. Dijkstra, K.; Pieterse, M.E.; Pruyn, A.T. Physical environmental stimuli that turn healthcare facilities into healing environments through psychologically mediated effects: Systematic review. *J. Adv. Nurs.* **2006**, *56*, 166–181. [CrossRef] [PubMed]

4. ANSI/ASHRAE. *Thermal Environmental Conditions for Human Occupancy Standard 55-2013*; ASHRAE: Peachtree Corners, GA, USA, 2013; ISSN 1041-2336.

5. Khodakarami, J.; Nasrollahi, N. Thermal comfort in hospitals–A literature review. *Renew. Sustain. Energy Rev.* **2012**, *16*, 4071–4077. [CrossRef]

6. Pereira, P.F.D.C.; Broday, E.E.; Xavier, A.A.D.P. Thermal Comfort Applied in Hospital Environments: A Literature Review. *Appl. Sci.* **2020**, *10*, 7030. [CrossRef]

7. ANSI/ASHRAE. *Standard 55: 2017, Thermal Environmental Conditions for Human Occupancy*; ASHRAE: Atlanta, GA, USA, 2017.

8. ISO/TS 14415. *Ergonomics of the Thermal Environment-Application of International Standards to People with Physical Disabilities*; International Organization for Standardization: Geneva, Switzerland, 2005.

9. Parsons, K. *Human Thermal Environments: The Effects of Hot, Moderate, and Cold Environments on Human Health, Comfort and Performance*, 2nd ed.; CRC: Boca Raton, FL, USA, 2001.

10. Fanger, P.O. *Thermal Comfort: Analysis and Applications in Environmental Engineering*; Danish Technical Press: Copenhagen, Denmark, 1970.

11. Mora, R.; Meteyer, M. Using thermal comfort models in health care settings: A review. *ASHRAE Trans.* **2018**, *124*, 11–23.

12. Westaway, K.; Frank, O.; Husband, A.; McClure, A.; Shute, R.; Edwards, S.; Curtis, J.; Rowett, D. Medicines can affect thermoregulation and accentuate the risk of dehydration and heat-related illness during hot weather. *J. Clin. Pharm. Ther.* **2015**, *40*, 363–367. [CrossRef]

13. de Dear, R.; Akimoto, T.; Arens, E.A.; Brager, G.; Candido, C.; Cheong, D.; Li, B.; Nishihara, N.; Sekhar, S.C.; Tanabe, S.-I.; et al. Progress in thermal comfort research over the last twenty years. *Indoor Air* **2013**, *23*, 442–461. [CrossRef]

14. Carlucci, S.; Bai, L.; de Dear, R.; Yang, L. Review of adaptive thermal comfort models in built environmental regulatory documents. *Build. Environ.* **2018**, *137*, 73–89. [CrossRef]

15. Nicol, J.F. Adaptive comfort. *Build. Res. Inf.** *Inf. 2011*, *39*, 105–107. [CrossRef]

16. Cheung, T.; Schiavon, S.; Parkinson, T.; Li, P.; Brager, G. Analysis of the accuracy on PMV-PPD model using the ASHRAE Global Thermal Comfort Database II. *Build. Environ.* **2019**, *153*, 205–217. [CrossRef]

17. Kenawy, I.; Elkadi, H. The impact of cultural and climatic background on thermal sensation votes. In Proceedings of the 29th international PLEA conference, Munich, Germany, 10–12 September 2013; pp. 1–6.

18. Pereira, P.F.; Ramos, N.M. Detection of occupant actions in buildings through change point analysis of in-situ measurements. *Energy Build.* **2018**, *173*, 365–377. [CrossRef]

19. Pereira, P.F.; Ramos, N.M.; Almeida, R.M.; Simões, M.L. Methodology for detection of occupant actions in residential buildings using indoor environment monitoring systems. *Build. Environ.* **2018**, *146*, 107–118. [CrossRef]

20. Verheyen, J.; Theys, N.; Allonsius, L.; Descamps, F. Thermal comfort of patients: Objective and subjective measurements in patient rooms of a Belgian healthcare facility. *Build. Environ.* **2011**, *46*, 1195–1204. [CrossRef]

21. Sattayakorn, S.; Ichinose, M.; Sasaki, R. Clarifying thermal comfort of healthcare occupants in tropical region: A case of indoor environment in Thai hospitals. *Energy Build.* **2017**, *149*, 45–57. [CrossRef]

22. Yau, Y.H.; Chew, B.T. Thermal comfort study of hospital workers in Malaysia. *Indoor Air* **2009**, *19*, 500–510. [CrossRef]

23. Hashiguchi, N.; Hirakawa, M.; Tochihara, Y.; Kaji, Y.; Karaki, C. Thermal Environment and Subjective Responses of Patients and Staff in a Hospital during Winter. *J. Physiol. Anthr. Appl. Hum. Sci.* **2005**, *24*, 111–115. [CrossRef]
24. Skoog, J.; Fransson, N.; Jagemar, L. Thermal environment in Swedish hospitals: Summer and winter measurements. *Energy Build.* **2005**, *37*, 872-877. [CrossRef]

25. Hill, L.D.; Webb, L.H.; Parsons, K.C. Carers' Views of the Thermal Comfort Requirements of People with Physical Disabilities. *Proc. Hum. Factors Ergon. Soc. Annu. Meet.* **2000**, *44*, 716-719. [CrossRef]

26. Khodakarami, J.; Knight, I. Measured thermal comfort conditions in Iranian hospitals for patients and staff. In *Proceedings of the Clima 2007 WellBeing Indoors, Helsinki, Finland, 10-14 June 2007; Seppaenen, O., Saeteri, J., Eds.; FINVAC: Helsinki, Finland, 2007*.

27. Kim, S.-K.; Hwang, Y.; Lee, Y.S.; Corser, W. Occupant Comfort and Satisfaction in Green Healthcare Environments: A Survey Study Focusing on Healthcare Staff. *J. Sustain. Dev.* **2015**, *8*, 156. [CrossRef]

28. Pereira, P.F.; Ramos, N.M.M. Occupant behaviour motivations in the residential context—An investigation of variation patterns and seasonality effect. *Build. Environ.* **2019**, *148*, 535-546. [CrossRef]

29. Walker, G.; Brown, S.; Neven, L. Thermal comfort in care homes: Vulnerability, responsibility and 'thermal care'. *Build. Res. Inf.* **2016**, *44*, 135-146. [CrossRef]

30. Ninomura, P.; Bartley, J. New ventilation guidelines for health-care facilities. *ASHRAE J.* **2001**, *43*, 29+30+32-33.

31. Hwang, R.-L.; Lin, T.-P.; Cheng, M.-J.; Chien, J.-H. Patient thermal comfort requirement for hospital environments in Taiwan. *Build. Environ.* **2007**, *42*, 2980-2987. [CrossRef]

32. Kameel, R.; Khalil, E. Thermal comfort vs air quality in air-conditioned healthcare applications, In *Proceedings of the 36th AIAA Thermophysics Conference, Orlando, FL, USA, 23-26 June 2003*.

33. Smith, R.M.; Rae, A. Thermal comfort of patients in hospital ward areas. *J. Hyg.* **1977**, *78*, 17-26. [CrossRef] [PubMed]

34. Bouzidi, Y.; El Akili, Z.; Gademer, A.; Tazi, N.; Chahboun, A. How Can We Adapt Thermal Comfort for Disabled Patients? A Case Study of French Healthcare Buildings in Summer. *Energies* **2021**, *14*, 4530. [CrossRef]

35. Pereira, P.F.; Ramos, N.M.; Ferreira, A. Room-scale analysis of spatial and human factors affecting indoor environmental quality in Porto residential flats. *Build. Environ.* **2020**, *186*, 107376. [CrossRef]

36. Mishra, A.K.; Ramgopal, M. Field studies on human thermal comfort—An overview. *Build. Environ.* **2013**, *64*, 94-106. [CrossRef]

37. EN 15251. *Indoor Environmental Input Parameters for Design and Assessment of Energy Performance of Buildings Addressing Indoor Air Quality, Thermal Environment, Lighting and Acoustics*; BSI: London, UK, 2006.

38. de Dear, R.J. A global database of thermal comfort field experiments. *ASHRAE Trans.* **1998**, *104*, 1141.

39. ISO 7726. *Ergonomics of the Thermal Environment-Instruments for Measuring Physical Quantities*; International Standardization Organization: Geneva, Switzerland, 1998.

40. ISO 7730. *Ergonomics of Thermal Environment-Analytical Determination and Interpretation of Thermal Comfort Using Calculation of PMV and PPD Indices and Local Thermal Comfort Criteria*; International Organization for Standardization: Geneva, Switzerland, 2006.

41. Atmaca, A.B.; Zorer Gedik, G. Determination of thermal comfort of religious buildings by measurement and survey methods: Examples of mosques in a temperate-humid climate. *J. Build. Eng.* **2020**, *30*, 101246. [CrossRef]

42. Ioannou, A.; Itard, L. In-situ and real time measurements of thermal comfort and its determinants in thirty residential dwellings in the Netherlands. *Energy Build.* **2017**, *139*, 487-505. [CrossRef]

43. Fabbri, K.; Gaspari, J.; Vandi, L. Indoor Thermal Comfort of Pregnant Women in Hospital: A Case Study Evidence. *Sustainability* **2019**, *11*, 6664. [CrossRef]

44. Zomorodian, Z.S.; Tahsildoost, M.; Hafezi, M. Thermal comfort in educational buildings: A review article. *Renew. Sustain. Energy Rev.* **2016**, *59*, 895-906. [CrossRef]

45. Yao, R.; Li, B.; Liu, J. A theoretical adaptive model of thermal comfort-Adaptive Predicted Mean Vote (aPMV). *Build. Environ.* **2009**, *44*, 2089-2096. [CrossRef]

46. Kim, J.T.; Lim, J.H.; Cho, S.H.; Yun, G.Y. Development of the adaptive PMV model for improving prediction performances. *Energy Build.* **2015**, *98*, 100-105. [CrossRef]

47. Hughes, C.; Natarajan, S.; Liu, C.; Chung, W.J.; Herrera, M. Winter thermal comfort and health in the elderly. *Energy Policy* **2019**, *134*, 110954. [CrossRef]

48. Neath, A.A.; Cavanaugh, J.E. The Bayesian information criterion: Background, derivation, and applications. *Wiley Interdiscip. Rev. Comput. Stat.* **2012**, *4*, 199-203. [CrossRef]

49. Montgomery, D.C. *Introduction to Statistical Quality Control*, 3rd ed.; John Wiley & Sons: New York, NY, USA, 1996.

50. Liu, C.; Zhou, G.; Li, H. Analysis of Thermal Environment in a Hospital Operating Room. *Procedia Eng.* **2015**, *121*, 735-742. [CrossRef]