Towards Model-Based Online Monitoring of Cyclist’s Head Thermal Comfort: Smart Helmet Concept and Prototype

Featured Application: In this work, we introduce the basis for a personalised adaptive model to predict head thermal comfort using streaming data of easily measured variables, which can be used for real-time monitoring of a cyclist’s thermal comfort and adaptive controlling of smart wearable applications. Abstract: Bicyclists can be subjected to crashes, which can cause injuries over the whole body, especially the head. Head injuries can be prevented by wearing bicycle helmets; however, bicycle helmets are frequently not worn due to a variety of reasons. One of the most common complaints about wearing bicycle helmets relates to thermal discomfort. So far, insu ﬃ cient attention has been given to the thermal performance of helmets. This paper aimed to introduce and develop an adaptive model for the online monitoring of head thermal comfort based on easily measured variables, which can be measured continuously using impeded sensors in the helmet. During the course of this work, 22 participants in total were subjected to di ﬀ erent levels of environmental conditions (air temperature, air velocity, mechanical work and helmet thermal resistance) to develop a general model to predict head thermal comfort. A reduced-order general linear regression model with three input variables, namely, temperature di ﬀ erence between ambient temperature and average under-helmet temperature, cyclist’s heart rate and the interaction between ambient temperature and helmet thermal resistance, was the most suitable to predict the cyclist’s head thermal comfort and showed maximum mean absolute percentage error (MAPE) of 8.4%. Based on the selected model variables, a smart helmet prototype (SmartHelmet) was developed using impeded sensing technology, which was used to validate the developed general model. Finally, we introduced a framework of calculation for an adaptive personalised model to predict head thermal comfort based on streaming data from the SmartHelmet prototype.


Introduction
Bicycling, for recreational, transport and sport purposes, provides health benefits for the individual as well society [1]. However, due to different reasons, bicyclists can be subjected to crashes, which can modelling body temperature response using metabolic activity alone or metabolic activity and skin temperature as inputs by means of, respectively, a single-input single-output (SISO) or a multiple-input single-output (MISO) discrete-time transfer function model.
Recent developments in compact wireless sensors allow the implementation of sensors in wearable devices such as bicycle helmets. Considering this, bicycle helmet design should be optimised for thermal comfort, so that bicycle helmets not only allow monitoring an individual's thermal comfort but also support its active control.
In this reported research, we aimed at (i) identifying a general model to estimate thermal comfort based on a few variables, the measurements of which can be integrated in helmets; (ii) developing and testing a prototype of a smart helmet based on the identified general thermal comfort model; and (iii) introducing the framework of calculation for an adaptive personalised reduced-order model to predict a cyclist's under-helmet thermal comfort using nonintrusive, easily measured variables.

Materials and Methods
The main goal of this paper is to introduce a framework for developing a personalised adaptive model for predicting a cyclist's head thermal comfort by utilising the smart helmet concept. Figure 1 presents the general framework introduced in the present paper.
Appl. Sci. 2019, 9, x FOR PEER REVIEW 3 of 19 and control applications. Data-based models, on the other hand, are less complex and thus more adequate for real-time monitoring and control purposes. Youssef et. al. [17] demonstrated that such compact data-based mechanistic models are promising for modelling body temperature response using metabolic activity alone or metabolic activity and skin temperature as inputs by means of, respectively, a single-input single-output (SISO) or a multiple-input single-output (MISO) discretetime transfer function model. Recent developments in compact wireless sensors allow the implementation of sensors in wearable devices such as bicycle helmets. Considering this, bicycle helmet design should be optimised for thermal comfort, so that bicycle helmets not only allow monitoring an individual's thermal comfort but also support its active control.
In this reported research, we aimed at (i) identifying a general model to estimate thermal comfort based on a few variables, the measurements of which can be integrated in helmets; (ii) developing and testing a prototype of a smart helmet based on the identified general thermal comfort model; and (iii) introducing the framework of calculation for an adaptive personalised reduced-order model to predict a cyclist's under-helmet thermal comfort using nonintrusive, easily measured variables.

Materials and Methods
The main goal of this paper is to introduce a framework for developing a personalised adaptive model for predicting a cyclist's head thermal comfort by utilising the smart helmet concept. Figure 1 presents the general framework introduced in the present paper.

Experimental Setup and Test Subjects
During the course of these experiments, 15 male test subjects with an average age ( ̅ ) of 22 (±1) years and an average weight ( ̅ ) of 74.3 (±9.2) kg were used in this study. The experimental protocol was approved by the Social and Societal Ethics Committee (SMEC) of KU Leuven. The experiments were conducted using a professional bicycle trainer (Tacx Ironman Smart) with a

Experimental Setup and Test Subjects
During the course of these experiments, 15 male test subjects with an average age (µ age ) of 22 (±1) years and an average weight (µ mass ) of 74.3 (±9.2) kg were used in this study. The experimental protocol was approved by the Social and Societal Ethics Committee (SMEC) of KU Leuven. The experiments were conducted using a professional bicycle trainer (Tacx Ironman Smart) with a fastened racing bicycle (BH L52C8 Speedrom) controlling the power delivery of the subject with a power brake.
Appl. Sci. 2019, 9, 3170 4 of 20 The power brake itself was wirelessly controlled via the Tacx Trainer software. The bicycle trainer was placed in a customised wind tunnel to simulate the wind effect on the test subjects during the course of the experiment (Figure 2). The wind tunnel was 2.1 m high, 2.3 m long and 1.5 m wide. Four rows with three fans each (Fancom type 1435/L7-588 fans) were used as the actuators for wind speed. Each fan produced a maximum ventilation rate of 3000 m 3 ·h −1 . A 50 cm long honeycomb gauze structure, placed 25 cm from the fans, was used to obtain a quasi-laminar flow within the open-loop wind tunnel (for more information about the wind tunnel, see [18]). The air speed near the test subject's head was set to 2.5 m·s −1 to simulate recreational cycling for adults and children. The wind tunnel was placed inside a climate-controlled chamber (Figure 2), the inner dimensions of which were 4 × 11 × 5 m (w × l × h). The air temperature within the climate chamber was controllable within the range of 15-35 • C. Additionally, the ventilation rate within the climate chamber was controllable within the range of 0-2700 m 3 ·h −1 (i.e., 0-11.25 volume refreshments per hour). fastened racing bicycle (BH L52C8 Speedrom) controlling the power delivery of the subject with a power brake. The power brake itself was wirelessly controlled via the Tacx Trainer software. The bicycle trainer was placed in a customised wind tunnel to simulate the wind effect on the test subjects during the course of the experiment ( Figure 2). The wind tunnel was 2.1 m high, 2.3 m long and 1.5 m wide. Four rows with three fans each (Fancom type 1435/L7-588 fans) were used as the actuators for wind speed. Each fan produced a maximum ventilation rate of 3000 m 3 ·h −1 . A 50 cm long honeycomb gauze structure, placed 25 cm from the fans, was used to obtain a quasi-laminar flow within the open-loop wind tunnel (for more information about the wind tunnel, see [18]). The air speed near the test subject's head was set to 2.5 m·s −1 to simulate recreational cycling for adults and children. The wind tunnel was placed inside a climate-controlled chamber (Figure 2), the inner dimensions of which were 4 × 11 × 5 m ( × × ℎ). The air temperature within the climate chamber was controllable within the range of 15-35 °C. Additionally, the ventilation rate within the climate chamber was controllable within the range of 0-2700 m 3 ·h −1 (i.e., 0-11.25 volume refreshments per hour).

Pretest Experiments
The pretest was a modification of the widely used physiological test protocols described by the Australian Institute of Sport [19]. The aim of the pretest was to obtain a power ( ) value that could be maintained by each of the 15 test subjects for a period of at least 20 min. The maximal lactate steady state (MLSS) and the corresponding workload steady state (WLSS) are presumed to be the maximum workload that can be maintained for endurance sports [20,21]. This lactate threshold is defined as the highest oxygen consumption rate that can be achieved during exercise without a systematic increase in blood lactate concentration [22]. A respiratory exchange ratio (RER) > 1.0 is an indication of the growing contribution of anaerobic metabolism, which causes muscle acidification and leads to muscle fatigue [23].
A bicycle incremental step test was designed in such way that the power increased 30 W every 5 min starting from 100 W. During the test, the subject's RER was measured with a spirometer (Metamax 3B) and the test was terminated when he exceeded an RER value of 1 for more than 20 s (the corresponding power, , was used further in the thermal comfort experiment). The tests were conducted at normal indoor climate conditions with 47% (±4%) relative humidity and an ambient temperature of 20 (±1) °C.

Pretest Experiments
The pretest was a modification of the widely used physiological test protocols described by the Australian Institute of Sport [19]. The aim of the pretest was to obtain a power (P) value that could be maintained by each of the 15 test subjects for a period of at least 20 min. The maximal lactate steady state (MLSS) and the corresponding workload steady state (WLSS) are presumed to be the maximum workload that can be maintained for endurance sports [20,21]. This lactate threshold is defined as the highest oxygen consumption rate that can be achieved during exercise without a systematic increase in blood lactate concentration [22]. A respiratory exchange ratio (RER) > 1.0 is an indication of the growing contribution of anaerobic metabolism, which causes muscle acidification and leads to muscle fatigue [23].
A bicycle incremental step test was designed in such way that the power increased 30 W every 5 min starting from 100 W. During the test, the subject's RER was measured with a spirometer (Metamax 3B) and the test was terminated when he exceeded an RER value of 1 for more than 20 s (the corresponding power, P RER = 1 , was used further in the thermal comfort experiment). The tests were conducted at normal indoor climate conditions with 47% (±4%) relative humidity and an ambient temperature of 20 (±1) • C.

Thermal Comfort and Variable Screening Experimental Protocol
The main objective of this stage was to screen the most suitable variables to predict the cyclist's under-helmet thermal comfort which can also be easily measured so as to be combatable for smart helmet application. During the course of these experiments, each experimental trial lasted 80 min and was divided into four consecutive timeslots of 20 min each. At each timeslot, a combination of Appl. Sci. 2019, 9, 3170 5 of 20 changes in the environmental variables, namely, relative air velocity imposed by the fan (v), ambient air temperature (T a ), thermal resistance of the scalp (R h ) and the delivered cycling power (P), was applied. The quantification of the scalp thermal resistance (R h ) was developed based on computational fluid dynamic (CFD) simulation for a bare head [24]. The thermal resistance was quantified (see Table 1) for the following cases: no-helmet wearing, where R h was 0; wearing helmet; and wearing helmet with helmet fast (aeroshell). The applied combinations of the different variables with their different levels (low, mid and high) are shown in Table 1. During the course of each trial, the heart rate (H R ), in bpm, of the test subject was measured and logged with a validated heart rate belt sensor (Zephyr TM bioharness Bt) in combination with a built-in optical heart rate sensor (PPG, Lifebeam) in the bicycle helmet (Lazer Z1 and Lazer Z1 fast = Lazer Z1 + aeroshell). The temperatures of the subject's forehead, neck, inside of the ear and the air under the bicycle helmet (at front, back, right and left) were continuously measured using calibrated thermocouples (type-T) with a sampling frequency of 1 Hz.
During the experiment, all test subjects were verbally asked about their thermal comfort every 5 min from the start (minute 0) until the end (minute 80) based on the thermal comfort scale introduced by Gagge et al. [10]. For convenience, the cold thermal sensation votes were excluded, as shown in Table 2, as the present work only focused on discomfort perception due to high temperatures. Table 2. Thermal comfort scale introduced by Gagge et al. [10], excluding the cold sensation votes.

Scale
Thermal Comfort Perception   1  Comfortable  2  Slightly uncomfortable  3  Uncomfortable  4 Very uncomfortable The experimental design was done using JMP Pro software. A preliminary screening experiment was set up to investigate the contribution of the different variables that, potentially, have an effect on the thermal sensation and thermal comfort under the bicycle helmet. Therefore, each subject was subjected to a combination of different levels of environmental conditions during the experiment.
The experiment was designed to investigate the main effects of the defined environmental input variables, the two-variable interactions between these variables and, due to the particular interest in the effect of a bicycle helmet, the quadratic effect of R h , which can be mathematically expressed as follows: where T C is the thermal comfort and β i is the weighting factor for each variable or variable combination (i). The inclusion of the quadratic effect, which is the interaction effect of the variable with itself, was necessary to generate an experiment that has multiple levels of R h , so that analysis of a dynamic response due to the bicycle helmet was possible. With the help of the JMP Pro ® software, different combinations (referred to as runs) of the input variables were generated. In general, each participant (test subject) was subjected to four runs (combinations) of the generated ones. Table 3 shows the experimental design for test subjects (j) 1 and 8 as an example, where each time slot corresponds to one run (a combination of the four input variables).

General Linear Regression (LR) Model Identification and Offline Parameter Estimation
The main objective of this stage was to identify a general reduced-order and the most parametrically efficient (parsimonious) model structure with the minimum number of easily measured variables (based on the results of the previous stage) to predict the cyclist's under-helmet thermal comfort. For the sake of the main objective of the present work, the selected predictive model had to be suitable, concerning the computational cost, for wearable sensing technology. Due to the subjective nature of the thermal comfort data, it could not be performed in a continuous pattern, unlike the other input variables, which was a challenge for identifying the predictive model. Hence, in the present paper, we used a simple multivariate regression model with the following general form [25]: where T ci ∈ R is the response (thermal comfort) for the ith observation, α ∈ R is the regression intercept, β j ∈ R is the jth predictor's slope, u ij ∈ R is the jth predictor for the ith observation and i ∼ N 0, σ 2 is an independent and identically distributed Gaussian error term. This can be formulated in matrix form as follows: T C = Xβ + , subjected to : T C ∈ R n×1 and X ∈ R n×m where n and m are the number of samples and number of predictors (input variables), respectively. In the present paper, we used the ordinary least-squares (OLR) approach to find the regression coefficients estimates (β) that minimised the sum of the squared errors as follows:

Development of Smart Helmet Prototype
A standard cyclist helmet (312 g) was utilised for the development of the smart helmet prototype (Lazer Bullet 1.0, Lazer Sport, Antwerp, Belgium). The helmet was equipped with a Lifebeam heart rate sensor (Lazer Sport, Antwerp, Belgium; Figure 3) and a 3 × 3 mm digital humidity and temperature sensor (CJMCU-1080 HTC1080, Texas Instruments, Dallas, Texas; accuracy: ±2% for relative humidity and ±0.2 • C for temperature) to measure the surrounding air humidity and temperature. Additionally, four temperature sensors (Negative-Temperature-Coefficient "NTC" temperature sensors, 100 kΩ at 25 • C; Figure 4) were used at the front, back, right and left of the helmet inner body. The final weight of the equipped helmet was 358 g. All sensors were connected directly to a microcontroller (Adafruit Feather Appl. Sci. 2019, 9, 3170 7 of 20 32u4 Bluefruit, Adafruit Industries, New York, NY, USA) that transmitted all data from the helmet to a smartphone via Bluetooth. The Adafruit Bluefruit was chosen as it is the smallest "all-in-one" Arduino-compatible and Bluetooth Low Energy microcontroller with built-in USB and battery charging. The developed system was compatible with a 3.7 V Li-polymer rechargeable battery (LP-523450-1S-3) with the ability to power the system for up to 10 h. A circuit diagram of the used electronics and sensors is shown in Figure 4. The impeded electronics and sensor technology in the smart helmet increased the final original weight of the helmet (312 g) by only 14.7% and did not alter the geometric and aerodynamic characteristics of the original standard helmet. As such, the developed smart helmet is comparable to the original standard helmet (Lazer Bullet 1.0, Lazer Sport, Antwerp, Belgium).
Appl. Sci. 2019, 9, x FOR PEER REVIEW 7 of 19 helmet inner body. The final weight of the equipped helmet was 358 g. All sensors were connected directly to a microcontroller (Adafruit Feather 32u4 Bluefruit, Adafruit Industries, New York, NY, USA) that transmitted all data from the helmet to a smartphone via Bluetooth. The Adafruit Bluefruit was chosen as it is the smallest "all-in-one" Arduino-compatible and Bluetooth Low Energy microcontroller with built-in USB and battery charging. The developed system was compatible with a 3.7 V Li-polymer rechargeable battery (LP-523450-1S-3) with the ability to power the system for up to 10 h. A circuit diagram of the used electronics and sensors is shown in Figure 4. The impeded electronics and sensor technology in the smart helmet increased the final original weight of the helmet (312 g) by only 14.7% and did not alter the geometric and aerodynamic characteristics of the original standard helmet. As such, the developed smart helmet is comparable to the original standard helmet (Lazer Bullet 1.0, Lazer Sport, Antwerp, Belgium).   helmet inner body. The final weight of the equipped helmet was 358 g. All sensors were connected directly to a microcontroller (Adafruit Feather 32u4 Bluefruit, Adafruit Industries, New York, NY, USA) that transmitted all data from the helmet to a smartphone via Bluetooth. The Adafruit Bluefruit was chosen as it is the smallest "all-in-one" Arduino-compatible and Bluetooth Low Energy microcontroller with built-in USB and battery charging. The developed system was compatible with a 3.7 V Li-polymer rechargeable battery (LP-523450-1S-3) with the ability to power the system for up to 10 h. A circuit diagram of the used electronics and sensors is shown in Figure 4. The impeded electronics and sensor technology in the smart helmet increased the final original weight of the helmet (312 g) by only 14.7% and did not alter the geometric and aerodynamic characteristics of the original standard helmet. As such, the developed smart helmet is comparable to the original standard helmet (Lazer Bullet 1.0, Lazer Sport, Antwerp, Belgium).   An android-based application "SmartHelmet App" ( Figure 5) was developed to simultaneously communicate with both the Adafruit Feather microcontroller and the Lifebeam heart rate monitor via a Bluetooth communication protocol. The SmartHelmet App was developed using the AppyBuilder online platform (App Inventor, Massachusetts Institute of Technology, Cambridge, Massachusetts, USA). The application was designed to receive, display in real-time and store all the data from the SmartHelmet at a 0.2 Hz sampling rate. An android-based application "SmartHelmet App" ( Figure 5) was developed to simultaneously communicate with both the Adafruit Feather microcontroller and the Lifebeam heart rate monitor via a Bluetooth communication protocol. The SmartHelmet App was developed using the AppyBuilder online platform (App Inventor, Massachusetts Institute of Technology, Cambridge, Massachusetts, USA). The application was designed to receive, display in real-time and store all the data from the SmartHelmet at a 0.2 Hz sampling rate.

Test Subjects
In total, seven well-trained male cyclists were recruited for the course of this experiment. Their average physical characteristics were as follows: age-34.5 (±5) years; body mass-74.5 (±7.3) kg; body height-177.6 (±5.4) cm; body mass index (BMI)-23.6 (±1.8) kg·m −2 ; and body surface area-1.9 (±0.1) m 2 . Prior to the trial, a signed written consent form was obtained from all participants after a detailed description of the protocol, discomforts and benefits. The experimental protocol was approved by the ethical review board at the University of Thessaly, School of Exercise Science in accordance with the recommendations of the Declaration of Helsinki.

Experimental Design and Protocol
Participants were exposed to a hot (34 °C and 56% relative humidity) environment and completed a 30 km cycling time-trial (TT) inside an environmental chamber. In addition, exposure to 800 W of solar radiation was simulated using compact source iodide (CSI) lamps, while a constant wind speed of 5.1 m·s −1 was provided with a large 80 cm diameter industrial fan positioned in front of the participant at a distance of 140 cm from the bicycle saddle. All participants were instructed to abstain from vigorous physical activity 24 h prior the experimental trial and consume at least 500 mL of water and a light meal 2 h before arrival at the laboratory.

Test Subjects
In total, seven well-trained male cyclists were recruited for the course of this experiment. Their average physical characteristics were as follows: age-34.5 (±5) years; body mass-74.5 (±7.3) kg; body height-177.6 (±5.4) cm; body mass index (BMI)-23.6 (±1.8) kg·m −2 ; and body surface area-1.9 (±0.1) m 2 . Prior to the trial, a signed written consent form was obtained from all participants after a detailed description of the protocol, discomforts and benefits. The experimental protocol was approved by the ethical review board at the University of Thessaly, School of Exercise Science in accordance with the recommendations of the Declaration of Helsinki.

Experimental Design and Protocol
Participants were exposed to a hot (34 • C and 56% relative humidity) environment and completed a 30 km cycling time-trial (TT) inside an environmental chamber. In addition, exposure to 800 W of solar radiation was simulated using compact source iodide (CSI) lamps, while a constant wind speed of 5.1 m·s −1 was provided with a large 80 cm diameter industrial fan positioned in front of the participant at a distance of 140 cm from the bicycle saddle. All participants were instructed to abstain from vigorous physical activity 24 h prior the experimental trial and consume at least 500 mL of water and a light meal 2 h before arrival at the laboratory.
Upon arrival at the laboratory, participants changed into their standardised cycling apparel and underwent basic anthropometric measurements. Body height was measured using a stadiometer (Seca 213; Seca GmbH & Co. KG; Hamburg, Germany), while body mass was determined with a digital weighing scale (Version 5.3 KERN & Sohn GmbH). BMI and body surface area were calculated from the measurements of body height and mass. After instrumentation, participants wore the SmartHelmet, entered the controlled environmental chamber and sat on the cycle for 10 min for a baseline period. Thereafter, they performed a 15 min warm-up followed by the 30 km TT. Participants were allowed to drink water ad libitum throughout the TT. No verbal encouragement was provided during the TT.
Cyclists performed the TT on an adjustable friction-braked cycle ergometer (CycleOps 400 Pro Serie Indoor Cycle, Fitchburg, MA, USA), which was combined with the commercially available software Rouvy (VirtualTraining, Vimperk, Czech Republic), allowing simulation of a route on a computer screen. During the 30 km TT, all cyclists were instructed to complete the race as fast as possible with free access to controlling power (W) and cadence (rpm). To simulate real cycling, participants could see their power, cadence and covered distance throughout the TT.
Ratings of perceived exertion (RPE) were reported with the 6-20 point Borg scale [26] before the baseline period, at the beginning of the warm-up period as well as at the start and end times of the TT. Thermal comfort (T C ) and thermal sensation (T S ) were measured at the same time points using 7-and 9-point scales, respectively [10].
The average power output, pedalling cadence and 30 km TT duration of all the participants (test subjects) are shown in Table 4.  Figure 6 shows the resulting P RER = 1 value for each test subject and the corresponding low and high levels of power. The corresponding low (P = 50% of P RER = 1 ) and high (P = 90% of P RER = 1 ) levels for each test subject were used in the screening experiments.

Pretest Experiments
Appl. Sci. 2019, 9, x FOR PEER REVIEW 9 of 19 Upon arrival at the laboratory, participants changed into their standardised cycling apparel and underwent basic anthropometric measurements. Body height was measured using a stadiometer (Seca 213; Seca GmbH & Co. KG; Hamburg, Germany), while body mass was determined with a digital weighing scale (Version 5.3 KERN & Sohn GmbH). BMI and body surface area were calculated from the measurements of body height and mass. After instrumentation, participants wore the SmartHelmet, entered the controlled environmental chamber and sat on the cycle for 10 min for a baseline period. Thereafter, they performed a 15 min warm-up followed by the 30 km TT. Participants were allowed to drink water ad libitum throughout the TT. No verbal encouragement was provided during the TT.
Cyclists performed the TT on an adjustable friction-braked cycle ergometer (CycleOps 400 Pro Serie Indoor Cycle, Fitchburg, MA, USA), which was combined with the commercially available software Rouvy (VirtualTraining, Vimperk, Czech Republic), allowing simulation of a route on a computer screen. During the 30 km TT, all cyclists were instructed to complete the race as fast as possible with free access to controlling power (W) and cadence (rpm). To simulate real cycling, participants could see their power, cadence and covered distance throughout the TT.
Ratings of perceived exertion (RPE) were reported with the 6-20 point Borg scale [26] before the baseline period, at the beginning of the warm-up period as well as at the start and end times of the TT. Thermal comfort ( ) and thermal sensation ( ) were measured at the same time points using 7and 9-point scales, respectively [10].
The average power output, pedalling cadence and 30 km TT duration of all the participants (test subjects) are shown in Table 4.  Figure 6 shows the resulting value for each test subject and the corresponding low and high levels of power. The corresponding low ( = 50% of ) and high ( = 90% of ) levels for each test subject were used in the screening experiments.  Figure 7 shows the acquired measurements from test subject 1, including environment-related variables, namely, ambient temperature (T a ), air velocity (ν), helmet thermal resistance and applied power level (R h ) and the applied mechanical work rate (P). Bioresponse-related variables, including heart rate (H R ), the temperature difference (∆T) between the average temperature beneath the helmet and the ambient air temperature, the temperature difference (∆T ear ) between the ear temperature and the ambient air temperature as well the thermal comfort (T C ), were considered.

Development of Offline (General) Thermal Comfort Model
Appl. Sci. 2019, 9, x FOR PEER REVIEW 10 of 19 Figure 6. Obtained power values of the pretest. These power values correspond to the power value when they exceeded a respiratory exchange ratio (RER) of one. Figure 7 shows the acquired measurements from test subject 1, including environment-related variables, namely, ambient temperature ( ), air velocity (ν), helmet thermal resistance and applied power level ( ) and the applied mechanical work rate ( ). Bioresponse-related variables, including heart rate ( ), the temperature difference (∆ ) between the average temperature beneath the helmet and the ambient air temperature, the temperature difference (∆ ) between the ear temperature and the ambient air temperature as well the thermal comfort ( ), were considered. The graphs show the environmental variables (left graphs), including the ambient air temperature ( , °C), fan set-points ( , 1 = 4 ms −1 ), the helmet wearing level (0 = no helmet, 0.5 = helmet and 1 = helmet + aeroshell) and the applied mechanical work rate (power) level ( , W). The measured variables related to the bioresponses of the test subject (right graphs) were heart rate ( , bpm), the temperature difference ( ∆ , °C) between the average temperature beneath the helmet and the ambient air temperature, the temperature difference (∆ , °C) between the ear temperature and the ambient air temperature and the thermal comfort (red line) and sensation (blue line) scores.

Development of Offline (General) Thermal Comfort Model
To investigate the effect of the different inputs on thermal comfort, different linear regression models (general models) were identified to estimate and predict the perceived under-helmet thermal comfort (output) using continuously measured variables (inputs), including the aforementioned environmental and bioresponse-related variables. The most suitable combination of input variables was selected by retaining only the input variables with a significant (p < 0.05) effect on thermal comfort. Additionally, the best model structure was selected based on two main selection criteria, namely, the goodness of fit ( ) and Akaike information criterion (AIC). The results showed that the most suitable LR model structure, with the highest goodness of fit (average = 0.87 ± 0.05) and The graphs show the environmental variables (left graphs), including the ambient air temperature (T a , • C), fan set-points (ν, 1 = 4 ms −1 ), the helmet wearing level (0 = no helmet, 0.5 = helmet and 1 = helmet + aeroshell) and the applied mechanical work rate (power) level (P, W). The measured variables related to the bioresponses of the test subject (right graphs) were heart rate (H R , bpm), the temperature difference (∆T, • C) between the average temperature beneath the helmet and the ambient air temperature, the temperature difference (∆T ear , • C) between the ear temperature and the ambient air temperature and the thermal comfort (red line) and sensation (blue line) scores.
To investigate the effect of the different inputs on thermal comfort, different linear regression models (general models) were identified to estimate and predict the perceived under-helmet thermal comfort T C (output) using continuously measured variables (inputs), including the aforementioned environmental and bioresponse-related variables. The most suitable combination of input variables was selected by retaining only the input variables with a significant (p < 0.05) effect on thermal comfort. Additionally, the best model structure was selected based on two main selection criteria, namely, the goodness of fit (R 2 ) and Akaike information criterion (AIC). The results showed that the most suitable LR model structure, with the highest goodness of fit (average R 2 = 0.87 ± 0.05) and lowest Akaike information criterion (average AIC = 138 ± 12), to predict the thermal comfort for all test subjects was as follows: The average parameter estimates, t-ratio and p-value of P > |t| for each selected input variable are given in Table 5. The results showed that the main effect of the thermal resistance R h was not significant (p > 0.05); however, the variable interaction of R h with T a showed a significant (p = 0.015) effect on the prediction of the under-helmet thermal comfort. To understand the interaction effect of R h and T a on the prediction of thermal comfort, a prediction trace analysis of the model [27] was employed using prediction the JMP ® profiler tool [28], as visualised in Figure 8. For convenience of this analysis, the values of each input variable were scaled (normalised) in such a way to lie in the closed interval [−1, +1], where -1 indicates the variable's low level and +1 indicates its high level (Figure 8). The scaling of each variable value i(k) was done according to the following formula:  The prediction trace analysis [28] of the developed model (2) was based on computing the predicted response as one variable was changing while the others were held constant at certain values. The results showed that the effect of R h was dependent on the level of T a . At a low level (−1) of ambient air temperature (T a = 20 • C), for a change in thermal resistance R h from a low level (−1) (i.e., no-bicycle helmet) to a high level (1) (i.e., using the Lazer-Z1 Fast), the predicted thermal comfort scale (Table 2) decreased by 0.5 thermal comfort units but was perceived as comfortable. However, at a high level (1) of ambient air temperature (T a = 30 • C), the comfort level increased by 0.5 thermal comfort units. This information is important for actively controlling under-helmet thermal comfort, which can be done by manipulating the helmet thermal resistance via, for instance, opening/closing some of the helmet's holes.
(MAPEs) of 8.4% and 11%, respectively. The MAPE is given by where is the number of data points and is the predicted thermal comfort. The values of the input variables were normalised in the range between -1 and 1, which correspond to low and high levels, respectively. Table 6 shows the average parameter estimates of the developed compact regression model (3) for the 15 test subjects. It can be noticed that both the mechanical work rate ( ) and air velocity ( ) disappeared from the compact model (3). The heart rate (HR) variable included in the compact model (3) directly linked to the applied mechanical work rate ( ), hence the effect of , included in model (2), translated by the bioresponse represented by HR (e.g., [29]) included in model (3). According to Newton's law of Figure 8. Visualisation of the model prediction traces showing the interaction effect of the thermal resistance (R h ) and ambient temperature (T a ) on the predicted thermal comfort. (a) When the temperature was low (20 • C), additional thermal resistance was perceived as comfortable. However, (b) when the temperature was high (30 • C), additional thermal resistance was perceived as uncomfortable. The values of the input variables were normalised in the range between -1 and 1, which correspond to low and high levels, respectively. Table 6 shows the average parameter estimates of the developed compact regression model (3) for the 15 test subjects.
As expected, the heart rate (H R ) of the test subjects was found to be highly correlated (Pearson's correlation coefficient, r = 0.85) with the power (P). Additionally, the heart rate was significantly correlated (r = 0.68) with the recorded thermal comfort for all 15 test subjects.
As expected, the temperature difference (∆T) between the average air temperature beneath the helmet (T h ) and the ambient air temperature (T a ) was correlated with both relative air velocity (ν) and helmet thermal resistance (R h ), with r = 0.82 and 0.78, respectively.
By employing both heart rate (H R ) and the temperature difference (∆T) as input variables to the linear regression model, the best model structure that gave the highest average goodness of fit (with average R 2 = 0.89 ± 0.04) and lowest Akaike information criterion (average AIC = 123 ± 7) was as follows: where ∆T = T h − T a and T h is the average air temperature under the helmet, which is calculated from the four temperature sensors located under the helmet. It can be noticed that the structure of model (3) is more compact, consisting of three input variables, compared with the structure of model (2), which consisted of five input variables. Model (3) showed better prediction performance for the thermal comfort level than model (2), which had maximum mean absolute percentage errors (MAPEs) of 8.4% and 11%, respectively. The MAPE is given by where N is the number of data points andT C is the predicted thermal comfort. It can be noticed that both the mechanical work rate (P) and air velocity (ν) disappeared from the compact model (3). The heart rate (H R ) variable included in the compact model (3) directly linked to the applied mechanical work rate (P), hence the effect of P, included in model (2), translated by the bioresponse represented by H R (e.g., [29]) included in model (3). According to Newton's law of cooling, temperature difference (∆T) is the driving force for the convective heat transfer (Q h ) between the cyclist's head and the ambient air. The heat flux (q) is proportional to ∆T and the convective heat transfer coefficient (h c ) links both variables as follows: The heat transfer coefficient (h c ,W·m 2 · • C) is a combination of the heat transfer coefficient of the air (h air ) and that of the helmet (h H = 1 R h ); hence, The heat transfer coefficients of the air (h air ) and the bicycle helmet ( 1 R h ) are dependent on air velocity (ν). Hence, it is clear that the effect of ∆T is inherently connected to the effect of both ν and helmet thermal resistance (R h ).
It can be concluded from the presented results that the input variables included in model (3), namely, temperature difference (∆T), heart rate (H R ) of the cyclist and the interaction variable [T a R h ] between ambient temperature (T a ) and helmet thermal resistance (R h ), were suitable enough to estimate the cyclist's thermal comfort (T C ) under the bicycle helmet. These selected variables were the basis for developing a reduced-order personalised model for real-time monitoring of a cyclist's thermal comfort under the helmet. Additionally, from a practical point of view, these three variables were suitable to be measured using integrated sensors in the cyclist's helmet, as is shown in the following subsection.

Testing the SmartHelmet Prototype and Validation of the Developed General Model
In Figure 9, the average ratings of perceived exertion (RPE), thermal comfort (T C ) and thermal sensation (T S ) values at the start and end times of the TT are presented for all seven test subjects. The average values (±standard deviation) of all used subjective ratings showed a significant (p < 0.05) increase at the end of the TT (RPE = 17.6 ± 0.5, T C = 2.6 ± 0.5 and T S = 4.4 ± 0.6) compared with their values at the start of the trial.
Appl. Sci. 2019, 9, x FOR PEER REVIEW 13 of 19 cooling, temperature difference (∆ ) is the driving force for the convective heat transfer ( ) between the cyclist's head and the ambient air. The heat flux ( ) is proportional to ∆ and the convective heat transfer coefficient (ℎ ) links both variables as follows: The heat transfer coefficient (ℎ ,W·m 2 ·°C) is a combination of the heat transfer coefficient of the air (ℎ ) and that of the helmet (ℎ = ); hence, The heat transfer coefficients of the air (ℎ ) and the bicycle helmet ( 1 ) are dependent on air velocity ( ). Hence, it is clear that the effect of Δ is inherently connected to the effect of both and helmet thermal resistance ( ).
It can be concluded from the presented results that the input variables included in model (3), namely, temperature difference (Δ ), heart rate (HR) of the cyclist and the interaction variable [ ] between ambient temperature ( ) and helmet thermal resistance ( ), were suitable enough to estimate the cyclist's thermal comfort ( ) under the bicycle helmet. These selected variables were the basis for developing a reduced-order personalised model for real-time monitoring of a cyclist's thermal comfort under the helmet. Additionally, from a practical point of view, these three variables were suitable to be measured using integrated sensors in the cyclist's helmet, as is shown in the following subsection.

Testing the SmartHelmet Prototype and Validation of the Developed General Model
In Figure 9, the average ratings of perceived exertion ( ), thermal comfort ( ) and thermal sensation ( ) values at the start and end times of the TT are presented for all seven test subjects. The average values (±standard deviation) of all used subjective ratings showed a significant (p < 0.05) increase at the end of the TT ( = 17.6 ± 0.5, = 2.6 ± 0.5 and = 4.4 ± 0.6) compared with their values at the start of the trial.  Figure 10 shows the real-time measured average temperatures ( ) under the helmet, average temperature difference ( ) between the average temperature under the helmet and the ambient air temperature and the average heart rate (HR) obtained during the TT from all seven test subjects using the developed prototype smart helmet.  The developed offline liner regression model (3) was used to estimate the thermal comfort ( ) of all seven test subjects based on the measurements acquired from the SmartHelmet prototype and for comparison with the thermal comfort subjective rating. The model was able to estimate the thermal comfort from all test subjects and revealed an average of 0.84 (±0.03). Model (3) was able to predict the cyclist's thermal comfort under the helmet and had a maximum MAPE of 10%. However, by retuning the model parameters using the data obtained from the TT experiment, the maximum MAPE was reduced to 7.8%.
The main advantage of the proposed model is that it is a conceptually simple yet very effective tool to explore linear relationships between a response variable (output) and a set of explanatory variables (input variables), which can be easily used for wearable technology such as the SmartHelmet. On the other hand, the disadvantage of such a model is the absence of the time component; in other words, the model is not able to explain the transient response of the output. Additionally, in practice, many factors can affect and change the relationship represented by the proposed model. These factors include helmet-related factors (e.g., helmet weight), other environmental conditions (e.g., wind direction) and personal-related factors, which were not included in the model (e.g., the surface area and contour of the cyclist's head). Hence, it is clear that such general models need to be adapted to new data (personal data) and different conditions for better performance. With the help of wearable sensing technologies (SmartHelmet) and streaming modelling algorithms, an adaptive personalised model can be developed for real-time monitoring of a cyclist's head thermal comfort. The developed offline liner regression model (3) was used to estimate the thermal comfort (T C ) of all seven test subjects based on the measurements acquired from the SmartHelmet prototype and for comparison with the thermal comfort subjective rating. The model was able to estimate the thermal comfort from all test subjects and revealed an average R 2 of 0.84 (±0.03). Model (3) was able to predict the cyclist's thermal comfort under the helmet and had a maximum MAPE of 10%. However, by retuning the model parameters using the data obtained from the TT experiment, the maximum MAPE was reduced to 7.8%.
The main advantage of the proposed model is that it is a conceptually simple yet very effective tool to explore linear relationships between a response variable (output) and a set of explanatory variables (input variables), which can be easily used for wearable technology such as the SmartHelmet. On the other hand, the disadvantage of such a model is the absence of the time component; in other words, the model is not able to explain the transient response of the output. Additionally, in practice, many factors can affect and change the relationship represented by the proposed model. These factors include helmet-related factors (e.g., helmet weight), other environmental conditions (e.g., wind direction) and personal-related factors, which were not included in the model (e.g., the surface area and contour of the cyclist's head). Hence, it is clear that such general models need to be adapted to new data (personal data) and different conditions for better performance. With the help of wearable sensing technologies (SmartHelmet) and streaming modelling algorithms, an adaptive personalised model can be developed for real-time monitoring of a cyclist's head thermal comfort.
In the following subsection, we introduce the framework of online model adaptation and personalisation (streaming algorithm) based on the easily measured variables obtained from the wearable sensors impeded in the SmartHelmet.

Introduction of Online Personalisation and Adaptive Modelling Algorithm
Most of the available thermal sensation and comfort predictive models (e.g., [30][31][32][33][34][35][36][37]) are static models. That is, they predict the average vote of a large group of people based on, for example, the 7-point thermal sensation scale instead of individual thermal comfort, and they only describe the overall thermal sensation/comfort of multiple occupants in a shared thermal environment. To overcome the disadvantages of static models, adaptive thermal comfort models aim to provide insights and opportunities to personalise the thermal comfort prediction of individuals [38]. The idea behind adaptive models is that occupants and individuals are no longer regarded as passive recipients of the thermal environment, but rather, they play an active role in creating their own thermal preferences [39]. The suggested linear regression model, represented by (3), in the present paper is considered as a global model, also called an offline model [40,41], for an adaptive personalised model to assess and predict individual thermal comfort under a cyclist's helmet. Figure 11 depicts the proposed steps for retuning and personalising the offline regression model (3). The suggested personalised adaptive tuning algorithm consists of the main components shown in Figure 11.
Appl. Sci. 2019, 9, x FOR PEER REVIEW 15 of 19 In the following subsection, we introduce the framework of online model adaptation and personalisation (streaming algorithm) based on the easily measured variables obtained from the wearable sensors impeded in the SmartHelmet.

Introduction of Online Personalisation and Adaptive Modelling Algorithm
Most of the available thermal sensation and comfort predictive models (e.g., [30][31][32][33][34][35][36][37]) are static models. That is, they predict the average vote of a large group of people based on, for example, the 7-point thermal sensation scale instead of individual thermal comfort, and they only describe the overall thermal sensation/comfort of multiple occupants in a shared thermal environment. To overcome the disadvantages of static models, adaptive thermal comfort models aim to provide insights and opportunities to personalise the thermal comfort prediction of individuals [38]. The idea behind adaptive models is that occupants and individuals are no longer regarded as passive recipients of the thermal environment, but rather, they play an active role in creating their own thermal preferences [39]. The suggested linear regression model, represented by (3), in the present paper is considered as a global model, also called an offline model [40,41], for an adaptive personalised model to assess and predict individual thermal comfort under a cyclist's helmet. Figure 11 depicts the proposed steps for retuning and personalising the offline regression model (3). The suggested personalised adaptive tuning algorithm consists of the main components shown in Figure 11.

Offline Linear Regression Model
As mentioned earlier, the linear regression model (3), which was developed based on the data obtained from the 15 test subjects, is the offline base model for online prediction of personal under-helmet thermal comfort. The general form of the offline linear regression model (3) is as follows: and is subjected to T C ∈ R n×1 and X ∈ R n×3 ,where T C is the output vector (n samples of thermal comfort votes); is the model residual vector, which consists of independent and Gaussian-distributed entries; and β and X are the regression vector (of the size 3) and predictor matrix (of the size n × 3), respectively, given by

Streaming Data
The availability of real-time sensor data from the developed SmartHelmet prototype allows for streaming data, which is processed via an online algorithm (stream processing) to adapt the offline model [42]. The streaming data includes newň samples of measured sensor data (new input matrixX) acquired from the SmartHelmet sensors, and new personal thermal comfort votes (new output vectoř T C ) acquired through an interactive query provided by the developed SmartHelmet App.

Online Parameter Estimation Algorithm
As explained earlier (Section 2.1.3), the general setting of regression analysis is to identify a relationship between a response variable (Y) and one or several explanatory variables (predictors) (X) by using a learning sample [43]. In a prediction framework, the main assumption for predicting Y on a new sample of X observations is that the regression model (with the general form Y = f (X) + , where represents the model residuals) is still valid. Unfortunately, this assumption is not valid in the present case, where the thermal comfort of the individual cyclist is strongly dependent on many personal-and time-dependent factors [11]. Therefore, in this study, we adapted the original regression model (3) to a new sample (observations) by estimating a transformation [41,43] between the original regression function ( f (X)) and the new one (f (X)) while still using the same model variables and structure. Ordinary least squares (LS) is one of the most popular regression techniques, which was used here for parameter estimation of the developed offline model. However, for the online parameter estimation and in the presence of unknown parameter changes, its adaptive versions-the sliding or moving LS, recursive least squares (RLS) and recursive partial least squares (PLS)-are widely used [41,44,45].
In the present paper, the RLS algorithm is suggested for online personalisation and adaptive modelling of under-helmet thermal comfort. The suggested RLS algorithm has the advantage of being simple and computationally efficient for wearable and adaptive sensing, which was the case in the present work. In the RLS algorithm, the new regression vectorβ, as in Equation (4), can be estimated recursively as follows [25,46,47]: P n+1 = P n − P n (X n+1 ) T X n+1 P n 1 + X n+1 P n (X n+1 ) T , β n+1 =β n + P n+1 (X n+1 ) T Ť C n+1 − X n+1βn where P n = X n T X n −1 . This recursive algorithm is efficient for cases where the regression vectorβ is a function of time (time varying). However, in the case of adaptive modelling with streaming data, due to the arrival of new samples, the influence of new observations decreases gradually and the ability to track the changes inβ will be lost. Hence, to mitigate this, the widely used and popular forgetting factor approach [48] is proposed in this paper. The approach of forgetting here is based on gradually discarding older data in favour of more recent information. In the least-squares method, forgetting can be viewed as giving less weight to older data and more weight to recent data [48,49]. Hence, the forgetting factor, λ, was introduced to (5) as follows [41]: where λ ∈ (0, 1]. The forgettiing factor λ operates as a weight, which diminishes for more remote data and expands for more recent data [48,49]. The main difference here between (5) and (7) is that in conventional RLS (5), the covariance vanishes to zero with time, losing its capability to keep track of changes in the regression vectorβ. In (7), however, the covariance matrix is divided by 0 ≤ λ < 1 at each update. This slows down the fading out of the covariance matrix [49].

Conclusions
In the present work, we aimed to develop a general model approach to predict a cyclist's head thermal comfort using nonintrusive and easily measured variables, which can be measured using impeded sensors in a bicycle helmet. During the first experimental stage, 15 participants were exposed to different levels of mechanical activity, ambient temperatures, helmet thermal resistance and wind velocities in order to develop a general model to predict a cyclist's head thermal comfort. The results showed that ambient temperature, average air temperature under the helmet, cyclist heart rate, cyclist mechanical work and helmet thermal resistance significantly influenced the cyclist's head thermal comfort. A reduced-order general linear regression model with three input variables, namely, temperature difference between ambient temperature and average under-helmet temperature, cyclist's heart rate and the interaction between ambient temperature and helmet thermal resistance, was the most suitable to predict the cyclist's head thermal comfort, showing a maximum MAPE of 8.4%. The developed general model structure was based on easily measured variables that can be measured continuously using impeded sensors in the bicycle helmet but is still of reduced order and low computational cost, which is suitable for streaming and adaptive modelling. Based on the selected model variables, a smart helmet prototype (SmartHelmet) was developed using impeded sensing technology as a proof of concept. The developed general model was validated using the developed SmartHelmet prototype. During the validation experimental phase, seven well-trained male cyclists were exposed to a hot (34 • C and 56% relative humidity) environment and completed a 30 km cycling TT inside an environmental chamber. The validation results showed that the developed general model was able to predict the thermal comfort of the seven participants and had a maximum MAPE of 10%. By retuning the model parameters, the maximum MAPE decreased to 7.8%. Finally, we introduced a calculation framework of an adaptive personalised model based on the developed general model to predict a cyclist's head thermal comfort based on streaming data from the SmartHelmet prototype.