info@biomedres.us   +1 (720) 414-3554
  One Westbrook Corporate Center, Suite 300, Westchester, IL 60154, USA

Biomedical Journal of Scientific & Technical Research

July, 2019, Volume 19, 4, pp 14523-14532

Research Article

Research Article

Dynamic Model-Based Monitoring of Human Thermal Comfort for Real-Time and Adaptive Control Applications

Ali Youssef, Nicolás Caballero and Jean Marie Aerts*

Author Affiliations

Division Animal and Human Health Engineering, Belgium

Received: July 05, 2019 | Published: July 17, 2019

Corresponding author: Jean Marie Aerts, Division Animal and Human Health Engineering, M3-BIORES, KU Leuven, Kasteelpark Arenberg 30 bus 2456, 3001 Heverlee, Belgium

DOI: 10.26717/BJSTR.2019.19.003350

Also View In:

Abstract

Thermal comfort and sensation are important aspects of the building design and indoor climate control as modern man spends most of the day indoors. Conventional indoor climate design and control approaches are based on static thermal comfort/sensation models that views the building occupants as passive recipients of their thermal environment. Assuming that people have relatively constant range of biological comfort requirements, and that the indoor environmental variables should be controlled to conform to that constant range. Recent advances in mobile technologies in healthcare, in particular wearable technologies (m-health) and smart clothing, have positively contributed to new possibilities in controlling and monitoring health conditions and human wellbeing in daily life applications. The wearable sensing technologies and their generated streaming data are providing a unique opportunity to understand the user’s behaviour and to predict future needs. Many advanced and accurate mechanistic thermoregulation models, such as the ‘Fiala thermal Physiology and Comfort’ model, are developed to assess the thermal strains and comfort status of humans.

However, the most reliable mechanistic models are too complex to be implemented in realtime for monitoring and control applications. Additionally, such models are using not-easily or invasively measured variables (e.g., core temperatures and metabolic rate), which are often not practical and undesirable measurements for monitoring during varied activities over prolonged periods. The main goal of this paper is to develop dynamic model-based monitoring system of the occupant’s thermal state and their thermoregulation responses under two different activity levels. In total, 25 test subjects were subjected to three different environmental temperatures, namely 5° C (cold), 20° C (moderate) and 37° C (hot) at two different activity levels (at rest and cycling). Metabolic rate, heart rate, average skin temperature, skin heat flux and aural temperature are measured continuously during the course of the experiments. The results have shown that a reduced-ordered (second order)s MISO-DTF including three input variables (wearables), namely, aural temperature, heart rate, and average skin heat flux, is best to estimate the individual’s metabolic rate (non-wearable) with mean-absolute-percentageerror of 8.7%.

A general classification model based on least-squares-support-vector-machine (LSSVM) technique is developed to predict the individual’s thermal sensation. For a 7-classes classification problem, the results have shown that the overall model accuracy of the developed classifier is 76% with a F1-score value of 84%.The developed thermal-state prediction model is promising to estimate the human’s thermal sensation/comfort status in real-time and yet reduced-order, which is suitable for wearable applications.

Keywords: Thermal Sensation; Thermal Comfort; Machine Learning; Prediction; Adaptive Controlling

Introduction

Thermal comfort (TC) is an ergonomic aspect determining the satisfaction about the surrounding environment and is defined as ‘that condition of mind which expresses satisfaction with the thermal environment and is assessed by subjective evaluation’ ASHRAE [1]. The effect of thermal environments on occupants might also be assessed in terms of thermal sensation (TS), which can be defined as ‘a conscious feeling commonly graded into the categories cold, cool, slightly cool, neutral, slightly warm, warm, and hot’ ASHRAE [1]. Thermal sensation and thermal comfort are both subjective judgements, however, thermal sensation is related to the perception of one’s thermal state, and thermal comfort to the evaluation of this perception ISO-10551 [2]. The assessment of thermal sensation has been regarded as more reliable and as such is often used to estimate thermal comfort Koelblen et al. [3]. Human thermal sensation is mainly depending on the human body temperature (core body temperature), which is a function of sets of comfort factors Enescu et al. [4,5].

These comfort factors are including indoor environmental factors, namely mean air temperature around the body, relative air velocity around the body, humidity, and mean radiant temperature to the body Parsons [5]. Additionally, some personal (individualrelated) factors, namely, metabolic rate or internal heat production in the body, which vary with the activity level and clothing thermos-physical properties (such as clothing insulation and vapour clothing resistance), are included. It should be mentioned that the individual thermal perception is deepening, as well, on psychological factors include naturalness (an environment where the people tolerate wide changes of the physical environment), expectations and short/long-term experience, which directly affect individuals’ perceptions, time of exposure, perceived control, and environmental stimulation Nikolopoulou et al. [6].

The most considered way to have an accurate assessment of TS is to ask the individuals directly about their thermal sensation perception Enescu et al. [4,5]. Thermal sensation mathematical models have been developed in order to overcome the difficulties of direct enquiry of subjects. The development of such models is mostly depending on statistical approaches that by correlating experimental conditions (i.e., environmental and person-related variables) data to thermal sensation votes obtained from human subjects Koelblen et al. [3,5]. Most of these models (e.g., PMV) are static in the sense that they predict the average vote of a large group of people based on the seven-point thermal sensation scale, instead of individual thermal comfort, they only describes the overall thermal sensation of multiple occupants in a shared thermal environment. To overcome the disadvantages of static models, adaptive thermal comfort models aims to provide insights in increasing opportunities for personal and responsive control, thermal comfort enhancement, energy consumption reduction and climatically responsive and environmentally responsible building design De Dear et al. [7,8].

The idea behind adaptive model is that occupants and individuals are no longer regarded as passive recipients of the thermal environment but rather, play an active role in creating their own thermal preferences De Dear et al. [7]. Besides regression analysis, thermal sensation prediction can also be seen as a classification problem where various classification algorithms can be implemented Lu et al. [8]. Recently, number of research work (e.g., Chaudhuri et al. [9]; Dai et al. [10]; Farhan et al. [11]; Huang, Yang, and Newman [12]; Kim et al. [13] have demonstrated the possibility of using machine learning techniques, such as support vector machine (SVM), to assess and predict human thermal sensation. It can be concluded based on the published work (see the recent literature review by Lu et al. [8] that classificationbased models have performed so well as regression models. Recent advances in mobile technologies in healthcare, in particular wearable technologies (m-health) and smart clothing, have positively contributed to new possibilities in controlling and monitoring health conditions and human wellbeing in daily life applications.

The wearable sensing technologies and their generated streaming data are providing a unique opportunity to understand the user’s behavior and to predict future needs Hussain, Kang, and Lee [14]. The generated streaming data is unique due to the personal nature of the wearable devices. However, the generated streaming data is forming a challenge related to the need of personalized adaptive models that can handle newly arrived personal data. Current HVAC control systems can be divided into two types: air temperature regulator (ATR) and thermal comfort regulator (TCR). Most TCR controllers use static models, mainly PMV, as a performance criterion. This paper is aiming to develop an adaptive model for real-time monitoring of human thermal states using personal non-intrusive sensing techniques. The developed model should be suitable for real-time adaptive controlling of indoor climate systems and smart wearable applications.

Materials and Methods

Experiments and Experimental Setup

Climate Chambers (Body & Mind Room): The “Body & Mind Room” is consisting of three climate-controlled chambers (A, B and C) designed and built to investigate the dynamic mental and physiological responses of humans to specific indoor climate conditions. The Body & Mind Rooms are experimental facilities at the M3-BIORES laboratory (Animal and Human Health Engineering Division, KU Leuven). The three rooms are dimensionally identical; however, each room is designed to provide different ranges of climate conditions as shown in (Table 1).

Table 1: The different temperature and relative humidity ranges that can be provided by the different Body & Mind (A, B and C).

Experimental Protocol: The experimental protocol used in the present study is designed in such way to investigate the subjects’ thermal and physiological responses to predefined three different temperature (low, normal and high) that under two levels of physical activities (seated = low and cycling = high). The three predefined temperatures (low = 5°C, normal = 24°C and high = 37°C) are chosen based on the thermal-comfort-chart of the ASHRAE [15] the effects on health according to the Wind Chill Chart for cold exposure (National Weather Service of the US) and for hot temperatures exposure according to Dewhirst et al. [16]. The conducted experiments are consisted of two phases (Figure 1), upper graph), namely, low activity and high activity phases. During the first experimental phase, low activity phase, the test subjects (while being seated = low activity) are exposed, during 55 minutes, to three levels of temperatures in the following order: normal, low, high and normal again (Figure 1). During the high activity phase, the test subjects is exposed to 15 minutes of light physical stress (80W of cycling on a fastened racing bicycle). During the course (75 minutes) of the active phase, each test subject is exposed to the predefined three temperature levels (Figure 1), lower graph). During each temperature level, starting from the normal level (24 °C), the test subjects are performed 15 minutes of cycling (with 80 W power) and followed 4 minutes of resting (seated). During the course of conducted experiments, the clothing insulation factor (Col) is kept constant at Col = 0.34 which accounts for a cotton short and t-shirt as a standard clothing for all test subjects. The experimental protocol is approved| by the SMEC (SociaalMaatschappelijke Ethische Comissie), on the 16th of January with number G-2018 12 1464.

Figure 1: Plots showing the climate chambers’ set-point temperatures programed during the 55 minutes low activity phase (upper graph) and the 75 minutes high activity phase (lower graph).

Test Subjects: In total 25 healthy participants (6 females and 19 males), between the age of 25 and 35 (average age 26 ± 4.2) years, with average weight and height of 70.90 (±12.70) kg and 1.74 (±0.10) m, respectively, are volunteered to perform the aforementioned experimental protocol.

Measurements and Gold Standards: During the course of the experiments, participants’ heart rate, metabolic rate, average skin temperature, heat flux between the skin and the ambient air, core body temperature represented by the aural temperature are measured continuously. The heart rate monitoring is performed using the Polar H7 ECG strap that is placed under the chest, with a sampling frequency of 128 Hz. The metabolic rate as metabolic equivalent tasks (METs) of each test subject is calculated based on indirect calorimetry using MetaMAX 3B spiroergometer sensor. The average skin temperature is calculated based on measurements from three body-placed, namely, scapula, chest and arm (Figure 2). The skin temperature measurements are performed using one Shimmer temperature sensor and two gSKIN ® bodyTEMP patches. Two heat flux gSKIN ® patches are placed on both the chest and the left arm (Figure 2). The skin temperatures and heat flux measurements are acquired at sampling frequency of 1 Hz. Core body temperature is estimated based on aural temperature measure measurements, which is performed using in-ear wireless (Bluetooth) temperature sensor (Cosinuss One) with a sampling rate if 1Hz. At the end of each applied temperature level during the course of both experimental phases, a thermal sensation questionnaire, based on ASHRAE 7-poins thermal scale, is performed for each test subject.

Figure 2: Sensor placement.

a) Ear channel for aural temperature measurement via the Cosinuss One,

b) Upper arm where skin temperature and heat flux are measured with the gSKIN patch,

c) Middle upper chest where skin temperature and heat flux are measured with the gSKIN patch,

d) Lower chest where heart rate is measured with the Polar H7,

e) Scapula where skin temperature is measured with the shimmer,

f) Mouth and nose where metabolic rate is measured via the MetaMAX-3B spiroergometer sensor.

Modelling and Classification

For the sake of present study, the measured variables are divided into wearables, which are easily measured variables using wearable sensors and gold standards (reference) variables, which are not suitable for wearable technologies. The wearables are including heart rate biomedres-openaccess-journal-bjstr , aural temperature biomedres-openaccess-journal-bjstr , average skin temperature biomedres-openaccess-journal-bjstr , skin heat flux biomedres-openaccess-journal-bjstr and ambient air temperature biomedres-openaccess-journal-bjstr. On the other hand, the gold standards are consisting of the core temperature biomedres-openaccess-journal-bjstr, which is driven from the aural temperature biomedres-openaccess-journal-bjstr, metabolic rate biomedres-openaccess-journal-bjstr and personal thermal sensation votes TS. The ultimate goal of this work is to develop an adaptive classification model to predict the individual thermal sensation depending, solely, on the wearables or estimated variables. Hence, both of the metabolic rate and core body temperature are estimated using online dynamic modelling approach (Figure 3). Then, the individual thermal sensation is predicted using a classification model (classifier) whose inputs are the wearables and estimated the metabolic rate and core body temperature (Figure 3).

Figure 3: Overview of the main steps to predict the individual thermal sensation.

Dynamic Modelling: Although the system under study (occupant’s thermoregulation) is inherently a non-linear system, the essential perturbation behaviour can often be approximated well by simple linearized Transfer Function (TF) models Young [17,18]. For the purposes of the present paper, therefore, the following liner, Multi-input, single-output (MISO) discrete-timesystems are considered to estimate metabolic rate and core body temperature Young [19] ,

biomedres-openaccess-journal-bjstr

where denotes the value of the asscociated variable at the th sampling instant; is the output variable; biomedres-openaccess-journal-bjstr are input variables, while biomedres-openaccess-journal-bjstr and biomedres-openaccess-journal-bjstr are appropriately defined polynomials in the backshift operator biomedres-openaccess-journal-bjstr , i.e., biomedres-openaccess-journal-bjstr and biomedres-openaccess-journal-bjstr is additive noise, a serially uncorrelated sequence of random variables with variancethat accounts for measurement noise. The Simplified Refined Instrumental Variable (SRIV) algorithm was utilised in the identification and estimation of the models (model parameters and model structure) Young and Jakeman [20]. Two main statistical measures were employed to determine the most appropriate model structure. Namely, the coefficient of determination , based on the response error; and YIC (Young’s Information Criterion), which provides a combined measure of model fit and parametric efficiency, with large negative values indicating a model which explains the output data well and yet avoids over-parameterisation Young, Chotai and Tych [21] . Additionally, the estimation performance of the selected models is evaluated used the mean-absolute-error (MAE) value.

Classification Model: To predict the individual thermal sensation, a classification model (classifier) is developed and trained based on the wearables and estimated variables (metabolic rate and core body temperature) that together with the thermals sensation votes (gold standard). A modified support vector machine (SVM) technique, namely, least-squares support vector machine (LS-SVM) is used to develop and train the thermal sensation classifier Suykens et al. [22]; Suykens and Vandewalle [23]. SVMs are originally presented as binary classifiers Suykens et al. [22] that assign each data instance biomedres-openaccess-journal-bjstr to one of two classes described by a class label biomedres-openaccess-journal-bjstr based on the decision boundary that maximises the margin biomedres-openaccess-journal-bjstr between the two classes. Generally, a feature map biomedres-openaccess-journal-bjstr is used to transform the geometric boundary between the two classes to a linear boundary biomedres-openaccess-journal-bjstr in feature space, for some weight vector biomedres-openaccess-journal-bjstr and biomedres-openaccess-journal-bjstr. The class of each instance can then be found by biomedres-openaccess-journal-bjstr where sign refers to the sign function.

Due to some computational complexities of standard SVM because of the quadratic programming problem, least squares support vector machines (LS-SVM) is presented to overcome such problem. LS-SVM, in contrast with standard SVM, relies on least squares cost function as follows:

biomedres-openaccess-journal-bjstr

such that; biomedres-openaccess-journal-bjstr where biomedres-openaccess-journal-bjstr errors such that biomedres-openaccess-journal-bjstr is proportional to the signed distance of x_i from the decision boundary, and biomedres-openaccess-journal-bjstr represents the regularisation constant. In LSSVM, instead of solving the quadratic programming problem a set of linear equations to be solved is sufficient to find the optimal solution of the classifier. The LS-SVMlab (Least Squares Support Vector Machine lab) Matlab-based toolbox is used to implement the LS-SVM classification algorithm Suykens et al. [22].

The performance of the classification model is determined based on accuracy, sensitivity, precision and F1-score as follows:

biomedres-openaccess-journal-bjstr

where TP, TN, FP and FN are the true positive, true negative, false positive and false negative, respectively.

Results and Discussion

Dynamic Modelling and Estimation of Individual’s Metabolic Rate

The average metabolic rate obtained from the 25 participants

at the temperature levels (24, 5 and 37 °C) during low and high activity phases are presented in (Table 2). Different combinations of inputs variables (wearables) are tested for best estimation of individual’s metabolic rate. The SRIV algorithm, combined with and selection criteria, suggested that a second-order MISO discrete-time TF with heart rate ( biomedres-openaccess-journal-bjstr ), average skin heat flux ( biomedres-openaccess-journal-bjstr ) and aural temperature (biomedres-openaccess-journal-bjstr ) as input variables is the best (with average = 0.89 ±0.04 and = − 13.62 ±2.33) to describe and estimate the dynamic behaviour of the individual’s metabolic rate. More specifically, the SRIV algorithm identified the following general MISO discrete-time TF model structure

Table 2: Average (±standard deviation) of the measured metabolic rate obtained from the 25 test subjects during low and high activity phases.

Note: 1MET = 1.163 W.kg-1

biomedres-openaccess-journal-bjstr

where biomedres-openaccess-journal-bjstr is the estimated metabolic rate and the numerator polynomials , and are of the following orders (number of zeros) 2, 3 and 2, respectively. While the system delays , and are varied from person to another (inter-personal) with average values of 1.4, 0.20 and 0.21 minutes, respectively. A simulation example of the developed estimation model (2) for one test subject during low activity experimental phase at normal temperature (24°C) depicted in (Figure 4).

Figure 4: A simulation example of the developed MISO-DTF model (2) to estimate the metabolic rate during the low activity experimental phase at normal temperature (24oC).

The estimation performance of the selected general MISO-DTF (2) is evaluated based on the mean-absolute-percentage-error

biomedres-openaccess-journal-bjstr

The results have shown that the developed general model is shown, for all test subjects, a higher average MAPE value (10 ±2.2 %) during the low activity phases than the average MAPE value (7.6 ±2.6 %) resulted during the high activity phases. The METs (metabolic equivalent tasks) are a measure, which accounts for a normalized form of energy expenditure per kilogram of mass. There are consensus that the measurement of metabolic rate might vary amongindividuals (interpersonal) up to 75% Byrne et al. [24] and even within the same day from morning to afternoon for the same subject (intrapersonal) up to 6%, though measurements on different days might be comparable on fasted subjects Haugen et al. [25] . Hence, a general estimation model of individual metabolic rate will not be efficient in this case. However, the general estimation performance of the suggested general MISO model can be enhanced by using the online adaptive form of the SRIV algorithm Garnier, Young and Gilson [26]. The online adaptive (closed-loop) SRIV algorithm is providing the possibility to personalise the developed general model by retuning the model parameters and model delays based on the streaming data acquired from the wearable sensors (Table 3).

Table 3: Average RT2, YIC model delays and MAPE for the selected MISO-DTF model to estimate the individual’s metabolic rate obtained from the 25 test subjects during low and high activity phases.

Classification Model and Prediction of Individual’s Thermal Sensation

In order to give an idea about the interaction relationship between considered variables, the correlation between all measured variables are calculated and represented in a colourmap correlation coefficient (r) matrix as shown in (Figure 5). High positive or negative correlation coefficient values, such as that between heat flux and skin temperature, are reflected strong interaction between these variables, which can affect the feature selection of the classification model. The classification model for predicting the individual’s thermal sensation is developed, based on LS-SVM approach, by training the classifier on 80% of the data points, while the rest of the data (20%) is used for testing. The model accuracy, sensitivity F1-score and overall confusion matrix are computed to evaluate the performance of the developed classifier. A feature space including all the measured and estimated input variables, namely, biomedres-openaccess-journal-bjstr and biomedres-openaccess-journal-bjstr . Additionally, other features are extracted by computing the variance, min, max, root mean squares (RMS) and first derivative (biomedres-openaccess-journal-bjstr where is the variable) of the aforementioned measured and estimated variables. The age and gender of the test subjects are also included in the feature spaces.

Figure 5: Colour-map of the correlation matrix representing the correlation levels (r∈[-1,1]) between the mean values of all measured variables, namely, heart rate (HR), aural (core) temperature, arm temperature, chest temperature, scapula temperature, heat flux from the arm skin (Arm HF), heat flux from the chest skin (Chest HF) and metabolic rate (MET’s).

(Figure 6) is showing the distribution of the participants’ thermal sensation votes at the three environment temperatures (24, 5 and 37 °C). The aforementioned figure is showing the ‘confusion space’, or the area in which the reported thermal sensation votes at the three environmental temperatures are overlapping. Such observation shows that the thermal perception may overlap even with big difference in the surrounding environmental temperatures. Such confusion space is a great challenge for any predictive model of thermal sensation especially for static models such as PMV. For the sake of the main objective of the present work, the computational cost of the developed algorithm should be low enough to be compatible with wearable technology and online modelling. Hence, a feature selection procedure is employed to obtain the most reduced-dimension model yet with best error performance. Feature selection here is based on evaluating all possible feature combinations and selecting the combination with best error performance. The feature selection step is resulted in a feature space including 25 features as shown in (Table 4).

Table 4: An overview of the selected feature space including the measured and estimated variables (six variables) and some operations on these variables (× = selected).

Figure 6: Distribution of the participants (25) thermal sensation votes at the three environment temperatures (24, 5 and 37 °C) showing the confusion space.

A classification model is developed based on the selected 25 input features and trained using the LS-SVM approach. The resulted confusion matrix from the developed classification model based on the selected feature space is shown in (Figure 6) and Figure 7. The overall error performance results of the developed classification model are presented in (Table 4). For a 7-classes classificatison problem, the developed classifier has shown an overall accuracy of 76% to predict the individual’s thermal sensation. The developed classifier has shown a high (84%) F1-score, which reflects low false positive and negative.

Figure 7: The resulted confusion matrix from testing the developed LS-SVM classifier. The diagonal represents the correctly classified data points.

Table 5 The error performance results of the developed general classification for each class separately are shown in (Table 6). The results showed that the error performances of classes 1, 2, 6 and 7 are very low see (Table 6), which can be attributed to the low number (0, 2, 4, 2, respectively) of obtained votes for these classes or in other word due to the uneven class distribution. Therefore, the overall F1-score is more reliable and efficient measure of performance than accuracy in this case. SVM is used in recent studies to assess occupant’s thermal demands Dai et al. (2017), and to predict thermal comfort/sensation Farhan et al. (2015).

Table 5: Overall error performance (accuracy, sensitivity, precision and F1-score) of the developed general LS-SVM classification model.

Table 6: The error performance (precision, sensitivity, and F1-score) of the developed general LS-SVM model for the 7-classes classification problem.

In this present paper, 25 participants are subjected to three different environmental temperatures, namely 5°C (cold), 20 °C (moderate) and 37 °C (hot) at two different activity levels, namely, at low level (rest) and high level (cycling at 80 W power). Metabolic rate, heart rate, average skin temperature (from three different body locations), heat flux and aural temperature are measured continuously during the course of the experiments. The thermal sensation votes are collected from each test subject based on ASHRAE 7-points questioner. The results have shown that a reduced-ordered (second order) MISO-DTF including three input variables (wearables), namely, aural temperature, heart rate, and average heat flux, is best to estimate the individual’s metabolic rate (non-wearable) with average MAPE of 8.7%. A general classification model based on LS-SVM technique is developed to predict the individual’s thermal sensation. For a 7-classes classification problem, the results have shown that the overall model accuracy and F1-score of the developed classifier are 76% and 84%, respectively. It is suggested in this paper that the model overall performance of the model can be enhanced by using a personalised adaptive classification algorithm based on streaming data from wearable sensors.

References

Research Article

Dynamic Model-Based Monitoring of Human Thermal Comfort for Real-Time and Adaptive Control Applications

Ali Youssef, Nicolás Caballero and Jean Marie Aerts*

Author Affiliations

Division Animal and Human Health Engineering, Belgium

Received: July 05, 2019 | Published: July 17, 2019

Corresponding author: Jean Marie Aerts, Division Animal and Human Health Engineering, M3-BIORES, KU Leuven, Kasteelpark Arenberg 30 bus 2456, 3001 Heverlee, Belgium

DOI: 10.26717/BJSTR.2019.19.003350

Also View In:

Abstract

Thermal comfort and sensation are important aspects of the building design and indoor climate control as modern man spends most of the day indoors. Conventional indoor climate design and control approaches are based on static thermal comfort/sensation models that views the building occupants as passive recipients of their thermal environment. Assuming that people have relatively constant range of biological comfort requirements, and that the indoor environmental variables should be controlled to conform to that constant range. Recent advances in mobile technologies in healthcare, in particular wearable technologies (m-health) and smart clothing, have positively contributed to new possibilities in controlling and monitoring health conditions and human wellbeing in daily life applications. The wearable sensing technologies and their generated streaming data are providing a unique opportunity to understand the user’s behaviour and to predict future needs. Many advanced and accurate mechanistic thermoregulation models, such as the ‘Fiala thermal Physiology and Comfort’ model, are developed to assess the thermal strains and comfort status of humans.

However, the most reliable mechanistic models are too complex to be implemented in realtime for monitoring and control applications. Additionally, such models are using not-easily or invasively measured variables (e.g., core temperatures and metabolic rate), which are often not practical and undesirable measurements for monitoring during varied activities over prolonged periods. The main goal of this paper is to develop dynamic model-based monitoring system of the occupant’s thermal state and their thermoregulation responses under two different activity levels. In total, 25 test subjects were subjected to three different environmental temperatures, namely 5° C (cold), 20° C (moderate) and 37° C (hot) at two different activity levels (at rest and cycling). Metabolic rate, heart rate, average skin temperature, skin heat flux and aural temperature are measured continuously during the course of the experiments. The results have shown that a reduced-ordered (second order)s MISO-DTF including three input variables (wearables), namely, aural temperature, heart rate, and average skin heat flux, is best to estimate the individual’s metabolic rate (non-wearable) with mean-absolute-percentageerror of 8.7%.

A general classification model based on least-squares-support-vector-machine (LSSVM) technique is developed to predict the individual’s thermal sensation. For a 7-classes classification problem, the results have shown that the overall model accuracy of the developed classifier is 76% with a F1-score value of 84%.The developed thermal-state prediction model is promising to estimate the human’s thermal sensation/comfort status in real-time and yet reduced-order, which is suitable for wearable applications.

Keywords: Thermal Sensation; Thermal Comfort; Machine Learning; Prediction; Adaptive Controlling