# Comfort modeling

Comfort modelling in semi-open spaces, with variable environmental conditions and potentially transient occupancy, requires a higher level of detail than is generally used in buildings. Conventional steady-state approaches of the PMV type are unsuitable for such environments and a transient model of human metabolism is then used. In the following article, we present a comfort indicator based on such a model, as well as an exploration of the influence of physiological variability.

## 1. Introduction

Many approaches have been developed to quantify the level of comfort of individuals in indoor spaces, often for productivity purposes: the first studies in the field consisted in comparing factories' indoor temperatures with production quality or absenteeism [H.M. Vernon et al., 1936 ].

In peculiar cases, the indoor comfort temperature can be determined by an empirical method. Thus the adaptive approach of [R. De Dear, 2001] used for naturally ventilated buildings gives a linear relationship between comfort temperature and the one-month moving average of the outside temperature, calibrated from thousands of measurements.

In indoor environments, the most commonly used indicator is [Fanger, 1970]'s Predicted Mean Vote which connects the heat flux (negative or positive) on the individual to the feeling of comfort or discomfort. This link was made possible by combining an equational approach, allowing the heat balance to be established, with a statistical regression on the comfort level of a sample of a few thousand people. However, recent work has shown the limitations of the PMV approach [van Hoof, 2008] which is based on a male morphology and causes significant shifts in the comfort range according to gender [Kingma, 2015].

The equational approach discussed here is the two-node model developed at the J.-B. Pierce Foundation (also called "Pierce's model") during the XXe century [Gagge, 1936]. The purpose was to meet military requirements, including determining the thermal stress resistance of soldiers or fighter pilots under highly variable atmospheric pressures [Nishi, 1977]. The formulation of the model is relatively simple and consists of representing the human body in the form of two concentric cylinders representing the centre of the body or "core" and the skin layer, surrounded by the clothing (see Fig. 1). In this model, metabolism is a system regulated in temperature by corrective actions (vasomotricity, sweating, perspiration) that evolve over time and according to the surrounding conditions.

Indeed, in outdoor or semi-outdoor environments, rapid variations in ambient conditions make steady-state approaches unsuitable for predicting comfort. This is all the more true when these spaces are occupied temporarily (e.g.outside passageways, very open buildings, stations): the temperature rise time of the metabolism is about one hour in summer and several hours in winter [Hoeppe, 2002].

As an illustration, Fig. 2 shows the simulation of body temperatures and skin wetness of an individual with Pierce's two-node model. From the state of thermal equilibrium, the individual is put in a warm environment with strong solar radiation for thirty minutes. The evolution of body temperatures under transient conditions (solid lines) shows an asymptotic evolution towards the steady state that would be obtained under these conditions (dotted lines). Skin wetness follows a slightly slower dynamic and the steady state is not reached in half an hour of exposure.

The physiological reactions generated by the environment on metabolism via convection, radiation and water exchanges allow to calculate the SET (Standard Effective Temperature). This is the operating temperature of a reference environment that would produce the same skin temperature and skin wetness as the actual environment being studied. This reference environment has a low air velocity, a relative humidity of 50% and the clothing is standardized with respect to the individual's activity.

The following sections provide a mathematical description of the human body and associated transfers in order to determine the SET from its physiological reactions. The effect of the environment on comfort will then be discussed, and the influence of physiological variability on the dispersion of results will be explored.

## 2. Activity and clothing characteristics

In this section, the morphological characteristics of humans and their clothing are described.

### 2.1 Metabolic activity

The metabolism is broken down into so-called "basal" metabolism, which maintains temperature, and activity metabolism, which is linked to the task performed. The term basal metabolic activity, $M$ [W], expressed as a function of height $H$ [m], of mass $m$ [kg] and of age and is calculated using the equation below from [VDI, 2008]:

$$M = a \times m^{0.75} \Big (1 +b \times (30.0 - \text{age} \big ) + c \times \big ( \frac{100 H}{m^{1/3}} - d) \Big )$$

For both male and female individuals, the constants a,b,c and d respectively take the values $a=3,45$, $b=0,004$, $c=0,01$, $d=43,4$ and $a=3,19$, $b=0,004$, $c=0,018$ and $d=42,1$.

The surface area of the human body (in m²) is determined from the empirical equation of Dubois which gives a relationship between size H and mass m :

$A = 0.203 \times m^{0.425} \times H^{0.725}$

A typical value for an average individual is $A \sim 1,8$ [m²].

The metabolic source term can be computed by two methods:

• either directly in[met] with 1[met] = 58.2[W/m²] of body surface area. A light office activity generally corresponds to 1[met],
• or as a combination of basal metabolism and activity metabolism related to the execution of the task, generally in [W] which is converted to [W/m²] using the two previous relationships.

The table below gives some orders of magnitude of the heat released as a function of activity.

Activity [W/m²] [met]
Rest, lying down 45 0,8
Rest, sitting 58 1
Standing, light office activity 70 1,2
Standing, light activity (lab work) 95 1,6
Medium standing activity (working on a machine) 115 2
Heavy work 175 3

Table: Metabolism and activity

### 2.2 Clothing properties

Clothing is taken into account via a thermal resistance $i_{clo}$ in [clo], such that 1 [clo] = 0.155 [m².K/W] (marked as $R_\text{cl}$ in this unit).

For example, the value of 1[clo] corresponds to a shirt, jacket and trousers suit. The table below gives the level of insulation in[clo] for some sets of clothing.

Clothing [clo]
Light summer clothing 0,3
Office clothing 0,7
Indoor winter clothing 1,0
Outdoor winter clothing 1,5

Correspondence between clothing and level of insulation in [clo]

Wearing clothing increases the thermal resistance between the external environment and the skin surface. However, it also increases the surface area exposed to heat transfer: the factor $f_{clo}$ is introduced in order to characterize the increase in the heat exchange surface related to the clothing:

$$f_{cl} = 1 + k \times i_\text{clo}$$

In the previous equation, $k = 0,2$ si $i_{clo} < 0,5$ and $k = 0,15$ else. The surface $A \times f_{cl}$ is thus increased by 20% for the wearing of a thermal resistance of 1[clo].

The amount of clothing worn varies according to the weather and individual factors. However, the work of [Schiavon, 2013] has shown that the level of individuals' clothing can be determined from the outside temperature at 6 a.m. on the same day, according to the law shown in Fig. 3.

## 3. Equation of human metabolism

In this section, the heat balance on the two nodes of the Pierce model is presented. The calculation of heat and mass transfers to the human body and the regulation of metabolic temperatures are then detailed. Finally, the calculation of the SET will be explained.

### 3.1 Governing equations of Pierce's model

The two-node model of [Gagge, 1971] considers the human body as two concentric cylinders (nucleus and skin) regulated in temperature by several mechanisms: perspiration, sweating, shivering, vasomotricity. Metabolic heat as well as energy dissipated as a result of the activity is represented by a source term. Two balance equations are used to calculate the evolution of the temperature of the skin's core as a function of radiative, convective and evaporative transfers with the external environment (conduction is neglected for standing posture, the contact surface with the ground being about 2% of that of the body). The figure below gives a representation of the cylindrical model and the heat flux applied:

A power balance is then established on the core part, such that the algebraic sum of the metabolic heat M, of the power exchanged with the skin layer (conduction $U_\text{sk}$ and blood flow $\dot{m}_b c_b$) and respiratory exchanges $Q_r$ is equal to the storage of heat in the core $S_c$ (where the index " c " stands for core).

$$M - (U_\text{sk} + \dot{m}_b c_b )(T_\text{c}-T_\text{sk}) - Q_\text{r} = S_c$$

In the above equation, the term $M$ stands for metabolic activity in [W/m²], $U_\text{sk}$ is the thermal conductivity of the tissues between core and skin (such that $U_\text{sk} =5.28$ [W/(m².K)]). $\dot{m} c_b$ is the product of the blood mass flow rate $\dot{m}$ and the heat capacity of blood $c_b = 4180$ [J/(kg.K)]. $Q_r$ is the total heat dissipated by breathing (sensible and latent).

At the skin surface, the storage $S_{sk}$ (skin) is the sum of the heat transferred from the core by conduction and by blood flow, of convective heat $C$radiant heat $R$ and evaporative loss $E$ :

$$(U_\text{sk} + \dot{m}_b c_b) (T_\text{c}-T_\text{sk}) - C - R- E = S_{sk}$$

In the above equation, the right hand side represents the energy storage in the skin layer. To the left of equality, the term $E$ represents the evaporative exchanges on the surface of the skin related to perspiration and sweating. The terms $C$ and $R$ respectively represent convective and radiative heat exchange. The equivalent electrical diagram is given below:

Skin mass varies according to the blood flow that swells the external tissues to a greater or lesser extent. The fraction $\alpha$ of the body mass contained in the outer cylinder of the model is introduced. Variations in skin temperature $T_\text{sk}$ and core temperature $T_c$ are then described by following differential equations:

$$S_\text{sk} = \alpha \times \frac{m c_p}{A} \frac{dT_\text{sk}}{dt}$$

$$S_c = (1-\alpha)\times \frac{ m c_p}{A} \frac{dT_\text{c}}{dt}$$

In the second equation, the thermal capacity of the body is $c_p= 3492$ [J/(kg.K)] and the coefficient $\alpha$ stands for the fraction of the body mass constituted by the skin, which varies according to the blood flow. Fig. 4 shows the complete model and fluxes.

Note : The percentage of fat mass $G$ of the human body can be taken into account by weighting the thermal capacities of the tissues, such that:

$$c_p = G \times c_g + (1 - G) \times c_b$$

where respectively $c_g = 2510$ [J/kg/K] the thermal capacity of fat and $c_b = 3650$ [J/kg/K] is the one of tissues.

By setting the initial conditions for $T_\text{sk},T_\text{c}$ as the set values $T_\text{sk}^\text{set},T_\text{c}^\text{set}$, it is possible to solve the two temporal differential equations of the two-node model. The model governed by this system thus gives the evolution of body temperatures according to the exchanges with the surrounding environment. The calculation of the various terms presented in the above balance sheets is explained in the following sections.

### 3.2 Respiratory exchanges

The pulmonary ventilation rate is given as a function of metabolic activity in met such that :

$$\dot{m}_{p} = 0.3 \times \text{met} \text{ [kg/m^2/h]}$$

If we convert for a direct calculation with $M$ in [W/m²] and a flow rate in [kg/s], we get with 1 [met] = 58.2 [W/m²]:

$$\dot{m}_{p}= M \times \frac{0.3}{58.2 \times 3600} = M \times 1.432 \times 10^{-6} \text{ [kg/m^2/s]}$$

The exhaled air temperature is given by the following correlation:

$$T_\text{exp}= 32.5 + 0.066 \times T_a + 1.99 \times 10 ^{-6} \times p_a$$

where $T_a$ is the air temperature and $p_a$ is the ambient air vapour pressure in [Pa].

The sensible heat dissipated is then worth:

$$C_\text{resp} = \dot{m}_{p} c_p (T_\text{exp}- T_a)$$

With $c_a = 1006$ [J/(kg.K)] for air, we find the well known equation [Thellier, 1989], [Fanger, 1970]:

$$C_\text{resp} = 1.432 \times 10^{-3} M (T_\text{exp}- T_a) \text{ [W/m^2]}$$

The above formulation, although widely used, fails to take into account the sensible energy used by the temperature rise of water vapour $w_a$ kg/kg] contained in the exhaled air. A rigorous expression of sensible losses by respiration, taking into account the specific heat of the vapour $c_v$ is the following:

$$C_\text{resp} = \dot{m}_p \times (c_\text{a} + w_\text{a} c_\text{v})\times (T_\text{exp}- T_\text{a})$$

The latent energy dissipated is calculated from the difference in absolute lung moisture $w_\text{exp}$ and the one of ambiant air $w_a$ :

$$E_\text{resp} =\dot{m}_{p} L_v (w_\text{exp}- w_a)$$

It is more common to work with vapour pressure than with absolute humidity. We therefore convert the previous equation with psychrometric relationships. According to the law of perfect gases, we know that the moisture content of the air depends on atmospheric pressure $p_\text{atm}$ and the ratio of molar masses of $M_\text{a}$ and water vapour $M_{v}$:

$$w = \frac{M_{v}}{M_\text{a}}\times \frac{p_a}{p_\text{atm} - p_a} \simeq 0.622 \times \frac{p_a}{p_\text{atm} - p_a} \text{ [kg/kg]}$$

Thus the equation of the dissipated latent energy becomes:

$$E_{\text{resp}} = 0.622 \times \dot{m}_p L_{\text{v}} \times \Big (\frac{p_\text{a}}{p_\text{atm} - p_\text{a}} - \frac{p_\text{v exp}}{p_\text{atm} - p_\text{v exp}} \Big ) \text{ [W/m^2]}$$

The vapour pressure of the exhaled air is $p\text{v exp}$ is at saturated vapour pressure $p\text{v exp} = p_\text{vs}(T_\text{exp})$ with the following correlation for the saturated vapour pressure:

$$p_\text{vs}(T) = 10^{\frac{2.7877 + 7.625 \times T }{241 + T}}$$

In general, the above relationship is simplified by considering that the exhaled and ambient vapour pressures (in the order of $3.10^3$ [Pa]) are low in front of atmospheric pressure in the order of $10^5$ [Pa]. Hence we set down:

$$p_\text{atm} - p_\text{v} \simeq p_\text{atm} - p_\text{v exp} \simeq p_\text{atm}$$

The latent energy dissipated then becomes:

$E_\text{resp} \simeq 0.622 \times \frac{\dot{m}_p L_v }{p_\text{atm}}(p_\text{v exp}- p_a)$
$E_\text{resp} \simeq 2.14 \times 10^{-5} \times M \times (p_\text{v exp}- p_a) \text{ [W/m^2]}$

Neglecting the sensible heat of water vapour and assuming low exhaled pressures in front of atmospheric pressure, has a limited influence: depending on the temperature and humidity conditions considered, the bias introduced is in the order of 2% to 7% between exact and approximate formulations of the sensible and latent losses (the distribution of losses is given in Fig. 6. This deviation may not be negligible insofar as the losses by breathing reach between $\sim 25$ and 30% of the metabolism at rest (see the experiment of [Cain, 1989]). Fig. 6 showing the error between formulations can be used to determine whether it is necessary to implement the exact equations.

### 3.3 Heat Transfer

This section presents the body's convective and radiative exchanges with the environment.

#### 3.3.1 Convective heat transfer

The natural convection exchange is calculated from the following correlation, weighting the convection coefficient as a function of ambient pressure:

$$h_c=3\times p_\text{atm}^{0.53}\text{ [W/(m^2.K)]}$$

When the air velocity $v$ is known, the following correlation is used, depending on the velocity and atmospheric pressure shown in Fig. 7:

$$h_c = 8.6 \times (v \times p_\text{atm}\times 10^{-5})^{0.53} \text{ [W/(m^2.K)]}$$

In practice, the maximum value between the two previous equations is taken as the exchange coefficient. The literature is prolific in correlations for the calculation of the convective exchange coefficient on the human being: on this topic we recommand the synthesis made by [Thellier, 1989].

Long wave (infra-red) radiant heat transfer between the surface at temperature $T$ and the surrounding walls at temperature $T_r$ are modelled via a linearized radiative exchange coefficient $h_r$ determined as follows:

$$R = \sigma \varepsilon f (T^4 - T_r^4)$$
$$R = \sigma \varepsilon f (T^2 - T_r^2) (T^2 + T_r^2)$$
$$R =\sigma \varepsilon f (T - T_r)(T + T_r) (T^2 + T_r^2)$$

Where $\varepsilon$ is the emissivity of the clothing and $f$ is a radiation efficiency factor taking into account the fact that some body surfaces do not face the surrounding walls (under the arms, inside the legs).

One then writes out $T_m = (T+T_r)/2$, yielding:

$$R \simeq \sigma \varepsilon f (T - T_r) \times 2 T_m \times 2 T_m^2$$
$$R\simeq 4 \sigma \varepsilon f T_m^3 \times (T - T_r) = h_r \times (T - T_r)$$

Thus we obtain the linearized radiative exchange coefficient $h_r$:

$$h_r = 4 \sigma \varepsilon f T_m^3 \text{ [W/(m^2.K)] }$$

At room temperature, with $\varepsilon=0.97$ and $f = 0.72$, we obtain $h_r \sim 4.7$ [W/(m².K)].

In outdoor spaces, the individual is regularly exposed to solar radiation. A simple way to integrate solar fluxes into the heat balance is to convert them into an increase in the average radiant temperature. The flux absorbed by the human body is therefore assumed to be equal to an average radiant temperature $T_r^*$ greater than $T_r$.

$$\alpha f_{p} \varphi_0 = \sigma \varepsilon (T_r^{*4}-T_r^4)$$
$$T_r^{*} = \sqrt[4]{ T_r^4 + \frac{\alpha f_{p} \varphi_0}{\sigma \varepsilon} }$$

Where $\varphi_0$ is the flux on a surface orthogonal to radiation, $\alpha \sim 0.7$ is the absorption coefficient of solar radiation, $f_p$ is a solar radiation projection factor that depends on the solar height $h$ (a good approximation is given by $f_p \simeq 0.233 \times cos(h) +0.067$ after [Helbig, 2013]).

#### 3.3.3 Total sensible heat transfer

Convective and radiative transfer are grouped into a global heat exchange term $h_g$ :

$$h_g = h_r + h_c$$

Thus the overall thermal resistance taking into account the exchange surface increment related to the wearing of clothing is:

$$R_\text{g} = \frac{1}{f_\text{cl} h_g}$$

The operating temperature is then calculated $T_\text{op}$, i.e. the weighted average of air and wall temperatures such that :

$$T_\text{op} = \frac{h_c T_c + h_r T_r}{h_c + h_g}$$

The operative temperature can be interpreted as the temperature of an enclosure where the air and walls are at the same temperature ( $T_a = T_r$) and where the convective and radiative exchanges with the individual would be equal to those of the real environment (in which $T_a \ne T_r$).

The clothing temperature $T_\text{cl}$ can then be determined by a heat budget equation at the skin surface and at the clothing surface:

$$T_\text{cl} = \frac{R_\text{g} T_\text{sk} + R_\text{cl} T_\text{op}} {R_\text{g}+ R_\text{cl}}$$

It should be noted that the form of the previous equation is implicit: the exchange coefficient $h_r$ depends on the clothing surface temperature $T_\text{cl}$ (see the equation of the radiative transfer coefficient). Similarly, the operative temperature is also calculated from the exchange coefficient $h_r$ (see two equations above), it is therefore necessary to iterate to obtain $T_\text{cl}$.

### 3.4 Mass transfer

The transfer of steam between the skin and the ambient air, via the layer of clothing is governed by the following equality, with $p_\text{sk} = p_\text{vs}(T_\text{sk})$ being the saturation vapour pressure at skin temperature:

$$E_\text{max} = \frac{p_\text{sk} - p_\text{a}}{R_e}$$

where $R_e$ is the total water vapour permeation resistance between surface and environment, calculated from the following formula:

$$R_e = \frac{R_g + R_\text{cl}}{LR \times i_m}$$

where $i_m$ is the moisture permeability index which by definition is the ratio between thermal resistance and vapour transfer resistance multiplied by the Lewis relationship $LR$. A commonly used value is $i_m = 0.38$However, for specific garments this value may be lower. For example, for soft PVC or neoprene clothing $i_m \sim 0.15$.

The Lewis expression characterizes the relationship between mass transfer and heat transfer. It is used to calculate the vapour transfer coefficient from the convective coefficient. In practice, $LR \sim 17 \times 10^{-3}$ [K/Pa], however the value may be affected by particular ambient conditions. Since Lewis' relationship has a significant impact on the calculation of SET*, its calculation is presented in the following section.

#### 3.4.1 Computing the Lewis' relation LR

The analogy between heat and mass transfer allows the following correlations to be written for the numbers of Nusselt $Nu$ (ratio of convective and conductive heat transfer at an interface) and Sherwood $Sh$ (ratio of convective and diffuse mass exchanges at an interface):

$$Nu = C \times Re^m \times Pr ^n= \frac{h_c L_c}{\lambda }$$
$$Sh = C \times Re^m \times Sc ^n = \frac{\beta L_c}{D}$$

Where $L_c$ is the characteristic length of the phenomenon, $\beta$ [m/s] is the mass transfer velocity, $D$ [m²/s] is the mass diffusion coefficient and $\lambda$ [W/(m.K)] is the fluid's conductivity. Prandtl's number $Pr$ is the ratio of momentum quantity diffusion to thermal diffusivity and Schmidt number $Sc$ is its equivalent for mass transfer. The ratio of the two numbers is then performed:

$$\frac{Sh}{Nu}= \Big( \frac{Sc}{Pr} \Big)^n$$

Hence replacing the numbers of Sherwood and Nusselt by their expressions we obtain the number of Lewis $Le$, ratio of thermal diffusivity $\alpha =\frac{\lambda }{\rho_a c_a}$ to the mass diffusivity of vapour in air $D$ :

$$\frac{\beta L_c}{D} \frac{\lambda }{h_c L_c}= \Big( \frac{Sc}{Pr} \Big)^n = Le^n$$

Reformulating the previous expression, we can isolate the transfer velocity $\beta$ in [m/s], with generally $n=1/3$ :

$$\beta = \frac{h_c \times Le^{n-1}}{\rho_a c_a}$$

The vapour flow rate $\dot{m}_v$ is proportional to the mass transfer rate $\beta$ and to the difference in vapour density $\rho_\text{v}$ between skin surface and air:

$\dot{m}_v = \beta \times (\rho_\text{v sk}-\rho_\text{v a}) = \beta \frac{M_\text{v}}{RT} (p_\text{sk}-p_\text{a}) \text{ [kg/(m^2.s)]}$

Where $M_\text{v} = 18$ [g/mol] is the molecular weight of water vapour, $R = 8.32$ [J/(mol.K)] being the constant of perfect gases.

The heat exchanged by vapour transfer on the bare skin is then worth, with $\rho_a$ the density of air [kg/m3]:

$$E = L_v \frac{h_c}{\rho_a c_a Le^{2/3}} \frac{M_\text{v}}{R T}(p_\text{sk}-p_\text{a}) \text{ [W/m^2]}$$

In this configuration, it is convenient to have a simple relationship between the convective heat exchange coefficient and the mass transfer related to the vapour pressure gradient. Then we put Lewis' relationship $LR$ such that:

$$LR =\frac{L_v}{\rho_a c_a Le^{\frac{2}{3}} } \frac{M_\text{v}}{RT} = \frac{L_v}{\rho_a c_a} \frac{D^{\frac{2}{3}}}{\alpha ^{\frac{2}{3}}} \times \frac{M_\text{v}}{RT}$$

The expression of the power exchanged by water vapour transfer $E$ becomes as a result:

$$E = LR \times h_c \times (p_\text{sk}-p_a) \text{ [W/m^2]}$$

À température ambiante une valeur couramment acceptée du nombre de Lewis at the surface of skin is $Le = 0.83$Mass diffusivity $D$ is calculated with following relation:

$$D = D_0 \times \bigg(\frac{T_\text{sk}+273.15}{T_0} \bigg)^{1.81} \times \frac{p_\text{atm}}{p_0} \text{ [m^2/s]},$$
$$D_0 = 2.26 \times 10^{-5} \text{ [m^2/s]},$$
$$T_0 = 273.15 \text{ [K]},$$
$$p_0 = 101325 \text{ [Pa]}.$$

A typical value for $D$ is, at 20°C:

$$D = 2.26 \times 10^{-5} \times \frac{293.15}{273.15} = 2.425 \times 10^{-5} \text{[m^2/s]}$$

A correction factor $K$ applies to compensate for the dry air vapour pressure gradient going in the opposite direction to the water vapour pressure gradient, i.e. the dry air also diffuses towards the skin surface [Fobelets,1988]. The factor $K$ is calculated as follows (for more details see [ASHRAE, Chapter6]):

$$K = \rho_\text{a} \times \frac{p_0}{p_\text{atm}} \times \frac{\ln(\rho_\text{a} / \rho_\text{v})}{\rho_\text{a} - \rho_\text{v}} \text{ [kg/m}^3\text{]}$$

For a skin temperature of 35°C / 100% humidity and air at 25°C / 50% humidity, this factor is $K=1.04$In such conditions, one obtains $LR = 17.4 \times 10 ^{-3}$ [K/Pa].

#### 3.4.2. Evaporative heat transfer

When $E_\text{max}$ is known, the heat transfer related to sweating writes:

$$E_\text{sw} = L_v \times \dot{m}_\text{sw}$$

Where $L_v$ is the latent heat of vaporization and $\dot{m}_\text{sw}$ is the sweat rate (detailed in section 3.5).

The skin wetness of sweating $\omega_\text{sw}$ is then introduced. It can be seen as the equivalent fraction of wet skin to achieve the required level of evaporation. It is defined by the ratio of the power evaporated by sweating to the surface of the skin $E_\text{sw}$ and the maximum evaporation heat transfer:

$$\omega_\text{sw} = \frac{E_\text{sw}} {E_\text{max}}$$

The phenomenon of diffusion through the skin tissues is also taken into account. In the model of Pierce, we consider that 6% of the total skin wetness $\omega$ relate to diffusion:

$$\omega = 0.06 + 0.94 \times \omega_\text{sw}$$

The diffusion heat loss $E_d$ is defined as:

$$E_d = \omega E_\text{max} - E_\text{sw}$$

The total latent transfer is then calculated $E$ as the sum of the losses by diffusion and sweating on the surface of the skin :

$$E = E_\text{sw} + E_\text{d}$$

In addition, the critical wetness, determining the skin surface for which the effectiveness of regulation by sweating and perspiration declines, is expressed as a function of air velocity:

$$\omega_c = 0.59 \times v^{-0.08}$$

When the wetness is greater than the critical wetness, the sweat can no longer be completely evaporated and drops are forming. The skin wetness is then capped at $\omega = \omega_c$ and we calculated the latent heat transfer as:

$$\omega_\text{sw} = \omega_c / 0.94$$
$$E_\text{sw} = \omega_\text{sw} E_\text{max}$$

Diffusion exchanges are still considered to occur on 6% of the skin surface:

$$E_\text{d} = 0.06 \times (1-\omega_\text{sw}) E_\text{max}$$

Latent heat transfer $E$ is then computed again after the equation of E \eqref{eq:E_total}

### 3.5 Temperature control of the metabolism

In Pierce's model, metabolism is a system regulated around the set core temperatures $T_\text{c}^\text{set}$ and the set skin temperature $T_\text{sk}^\text{set}$. Various regulatory actions are used to maintain the body at a constant temperature or to prevent it from heating up.

#### 3.5.1 Blood vessel expansion/striction

Blood flow rate variation $q_b$ L/(m².h)] is described as a function of the difference between core and skin temperatures and their respective setpoints:

$$q_\text{b}= \frac{q_\text{b}^\text{set}+ C_\text{d} (T_\text{c} - T_\text{c}^{\text{set}})}{1 +C_s(T_\text{sk}^{\text{set}} - T_\text{sk}) }) \mathrm{ [L/(m^2.h)]}$$

Where the dilation and striction constants are valid $C_\text{d} = 200$ [L/(m².h.K)] et $C_\text{s} = 0.5$ [-]. The set core and skin temperatures are respectively $T_\text{c}^{\text{set}}=36.8$ [°C] et $T_{\text{sk}}^{\text{set}}=34$ [°C]. The blood flow rate is between a minimum of 0.5[L/(m².h.K)] and a maximum of 90[L/(m².h.K)].

If the temperature differences in the above equation become negative, their value is assigned to zero in order to avoid inconsistent operation of the temperature control. In it were not the case, a vasoconstriction action could take place when the skin temperature is higher than the set point, whereas on the contrary, the blood flow should be increased in order to reduce the deviation from the set point.

There is a great disparity between the values of expansion constants in the literature: thus [Hoeppe, 1993] proposes $C_\text{d}=75$ while [ASHRAE_2014] gives $C_\text{d}=120$ L/(m².h.K)] and [Doherty Arens, 1988] propose $C_\text{d}=200$. This illustrates the variability in nature, the influence of which is discussed in Section 5.

#### 3.5.2 Sweating

Sweating is triggered according to global body temperature $T_b$the latter being an average of the core and skin temperatures weighted by the mass fraction of skin $\alpha$ :

$$T_b = \alpha T_\text{sk} + (1-\alpha) T_{c}$$

where $\alpha$ depends on the blood flow that modifies the mass of the outer cylinder [Gagge, 1986]:

$$\alpha = \frac{0.0417737 + 0.7451832 }{q_\text{b} + 0.585417} \text{ [-]}$$

The sweat rate is then described as a function of the deviation of body and skin temperatures from their set points:

$$\dot{m}_\text{sw} = C_\text{sw} \times (T_\text{b} - T_\text{b}^\text{set}) \times e^{ \frac{(T_\text{sk} - T_\text{sk}^{\text{set}})}{10.7}} \text{ [kg/(m^2.h)]}$$

With $C_\text{sw}=0.170$ kg/(m².h.K)] from [Gagge, 1986]. The body setpoint temperature is defined as follows $T_b^\text{set} = 0.1 \times T_\text{sk}^\text{set} + 0.9 \times T_\text{b}^\text{set} = 36.49$°C.

If the temperature differences in the above equation become negative, they are set to zero to prevent sweating when body temperatures are below their set point, i.e. when the body or skin is colder than their set point temperature.

#### 3.5.3 Shivering

When the environment is cold, the shivering reaction allows heat to be generated by muscle contraction. In Pierce's model, this corresponds to the addition of a source term $M_\text{sh}$ adding to metabolic activity $M$defined as:

$$M_\text{sh} = C_\text{sh} \times (T_\text{sk}^{\text{set}} - T_\text{sk}) (T_\text{c}^{\text{set}} - T_\text{c})$$

The shivering coefficient is $C_\text{sh} = 19.4$ [W/(m².K²)]. Temperature differences are considered as zero when they become negative. A different formula taking into account the effects of fat percentage can also be applied [ASHRAE, Chapter9].

### 3.6 Computing the SETStandard Effective Temperature)

The Standard Effective Temperature (SET) is an indicator derived from Pierce's model that is used to quantify the comfort level of a real environment based on an individual's physiological reactions. This is an operative temperature.

First, Pierce's model is used to calculate the skin temperature $T_\text{sk}$ and the skin wetness $\omega$ that result from an individual's exposure to the real environment. The calculation is performed for an exposure time of 60 minutes with an integration time step of one minute. The skin temperature and skin wetness obtained are then stored.

The individual is then placed in a reference environment under controlled environmental conditions and the SET is calculated: it is the operative temperature of this environment that would produce the same physiological reactions ($T_\text{sk},\omega$) as the real environment.

The reference environment corresponds to an indoor office environment and is defined with :

• a low convective exchange coefficient, related only to the subject's activity,
• a standardized clothing, adapted to the level of metabolic activity,
• $\varphi = 50\%$ of relative humidity, which affects the vapour pressure $p_a$ of the reference environment,
• the radiant temperature is equal to the air temperature: this is the SET that is to be computed according to the skin temperature and skin wetness initially obtained in the real environment.

The convective heat exchange coefficient $h_c$ is defined as a function of metabolism, such that:

$$h_c =3\text{ pour } met < 0.85$$
$$h_c=5.66\times(met -0.85)^{0.39}\text{ pour } met \ge 0.85$$

The clothing is adapted to the level of metabolic activity in met such that:

$$i_\text{clo}=\frac{1.52}{\frac{M}{58.2} + 0.6944} - 0.1835 \text{ [clo]}$$

The heat flux exchanged with the environment are then recalculated according to the definitions given in sections 3.3 and 3.4.2, taking into account the changes made in the clothing resistances $R_\text{cl},R_e$ and in the convective heat transfer coefficient $h_c$. SET is obtained by solving the following non-linear equation, where $T_\text{sk}$ and $\omega$ are the physiological reactions calculated in the real environment:

$$C+R + E = \frac{1}{h_g + R_\text{cl}} \times (T_\text{sk} - SET^* ) + \frac{\omega }{R_e}\times \Big (p_\text{sk} - 50\% p_\text{vs}(SET^*) \Big)$$

The representation of SET as a function of humidity and temperature is given in following figure. One can note a dependence of comfort on humidity: at the same operating temperature, SET is greater for a higher humidity than for a low humidity. The comfortable SET range was determined by comparing the model with an experimental campaign and ranges from 22.2 to 25.6°C [Nishi, 1977].

## 4. Effect of environment on comfort

Semi-open spaces are subject to wide variations in ambient conditions, including air velocity and average radiant temperature. This section therefore focuses on their effects on simulated physiological reactions.

### 4.1 Air velocity and clothing properties

Thermal resistance and vapour permeation resistance are actually affected by the "pumping" effects associated with the wearer's movement, which promote air intake into clothing, and by wind-related air movements. This reduces their convective properties and their resistance to water vapour transfer. The following paragraphs describe the calculation of the decrease in these resistances according to ISO 9920, the genesis of which is described in [Holmer, 1999; Havenith, 1999].

#### 4.1.1 Variation in thermal resistance

The decrease in resistance to heat transfer is defined in [Holmer, 1999]. For a walking velocity smaller than $v_w < 0.7$ m/s] (i. e. 2.5[km/h]) an equivalent walking speed is defined $v_w$ depending on the activity level $M$ :

$$v_w = 0.0052 \times (M - 58)$$

For the convective heat transfer coefficient $h_c$ the corrected velocity $v_c$ should be used, taking into account the air velocity $v$ and the walking velocity $v_w$ :

$$v_c = \sqrt{ v^2 + v_w ^2}$$

The corrected total convective thermal insulation and insulation of clothing $R_T= R_g + R_\text{cl}$, which is valid between 0.2 and 3 [m/s] and $clo \ge 0.6$ is provided from following correlation:

$$\frac{R_T^{c}}{R_T} = e^{0.043 - 0.398 \times v_c + 0.066\times v_c^2 - 0.378 \times v_w + 0.094 \times v_w^2}$$

Where $R_T$ is the thermal insulation of clothing and the ambient air layer, with the thermal resistance value of clothing at zero air velocity $R_\text{cl}$. The convective thermal resistance is calculated with $v \sim 0$ [m/s].

An illustration of the ratio $R_T^{c}/R_T$ is shown in Fig. 9 as a function of air and walking velocities. It can be seen that an air velocity of $\sim 1$ m/s] at zero walking speed results in a significant reduction of $\sim 25 \%$ resistance to thermal transfer. It is therefore important to recalculate the heat transfer taking this reduction into account in order to obtain an accurate heat balance.

#### 4.1.2 Variation of the vapour permeation resistance

Corrected vapour transfer permeability $i_m^c$ is given by [Havenith, 1999]:

$$\frac{i_{m}^{c}}{i_{m}} = 4.9 - 6.5 \Big( \frac{R_T^{c}}{R_T} \Big) + 2.6 \Big ( \frac{R_T^{c}}{R_T} \Big)^2$$

The variation in the permeability index $i_m$ depending on $\frac{R_T^{c}}{R_T}$ is shown in Fig. 10. This illustration shows that for a 70% reduction in heat transfer, vapour permeability $i_m^c$ is multiplied by $\sim 1.6$ compared to $i_m$ in a quiet atmosphere. The equation giving the resistance to steam transfer $R_e$ is modified and has a significant impact on the latent energy balance.

Fig. 11 shows the difference between SET calculated with constant (solid line) and variable (dotted line) clothing properties calculated from the equations \eqref{eq:corr_T} and \eqref{eq:im_corr} for an air velocity of $v= 1$ [m/s]. A shift from the comfort zone to higher temperatures occurs: clothing properties are reduced, sensitive and latent transfers are increased and the operating temperature must be higher in order to obtain the same flow lost on the skin surface.

### 4.2 Influence of the average radiant temperature

The equation \eqref{eq:T_mrt_star} gives the average radiant temperature increment $\Delta T_r$ depending on the solar flux incident on the person. Thus, for 600[W/m²] of incident flux, we observe an increment of $\sim 26$ [K] with respect to the radiant temperature without sunlight.

An illustration of the influence of the average radiant temperature and air velocity is given in Fig. 13, where the SET comfort indicator is given spatially. The radiant temperature is calculated from the ground temperatures and the incident solar flux. For this example, air velocities were calculated by the mechanics of the computational fluids. It has been assumed that air temperature and humidity are not affected by the urban environment. The disparity of SET shows the important effect of air velocities and average radiant temperatures.

## 5. Influence of physiological variability on SET

As we have seen in the previous sections, SET is calculated from an individual's physiological reactions. In the standard [ASHRAE, 2013], these reactions are those of a male individual. As the morphological parameters are adjustable, it is possible to use the model for other phenotypes and to evaluate the comfort resulting from their physiological reactions.

The model presented in section 3.5 has a large number of parameters. In order to show the sensitivity of the model to them, the effects of the variability of physiological parameters on the SET comfort indicator evaluated from random draws are presented here. The dispersion of the results will be assessed in terms of standard deviation: it should be kept in mind that $\pm 2$ standard deviations around an average value correspond to 95% of the population assessed.

### 5.1 Presentations of variable physiological parameters

Many physiological characteristics can be accurately represented by Gaussian statistical distributions (e.g. size, mass). Assuming that this is also the case for other model parameters, the Monte Carlo method was used to assess the dispersion of SET resulting from physiological variability.

The variable parameters chosen for the variability assessment are as follows:

• size, mass and percentage of fat mass,
• the coefficients $C_\text{d}$, $C_\text{s}$, $C_\text{sw}$, $C_\text{sh}$ respectively for dilation, tightening, sweating and tremor,
• the conductance of fabrics $U_\text{sk}$ and basal blood flow $q_\text{b}^\text{set}$,
• the setpoint core and skin temperatures $T_\text{c}^\text{set}$, $T_\text{sk}^\text{set}$.

The dispersion around mass and size was chosen from data from INSEE. The average value of the other parameters is that of the original code. For the target skin and core temperatures, the standard deviation is considered to be 0.5% of the mean value in order to respect the minimum and maximum internal temperature values recorded by [Mackowiak, 1992].

For the other parameters, the standard deviation is taken at 5% of the mean value, which is a conservative assumption considering the disparities observed for example for the expansion constant (see box 3.5.1). The conductance of fabrics can also vary more significantly [Schweiker, 2017]. The table below gives a summary of the numerical values used.

Paramètre Valeur moyenne µ Écart-type sigma
Height [m] 1.75 0.065
Weight [kg] 71.4 8.9
Fat [%] 23 7.5
$C_d$ [L/m².h.k] 200 10
$C_s$ [1/K] 0.5 0.025
$C_{sw}$ [kg/(m².h.K)] 170 8.5
$C_{sh}$ [kg/(m².K)] 19.4 0.97
$U_{sk}$ [W/(m².K)] 5.28 0.264
$q_b^{\text{set}}$[L/(m².h)] 6.3 0.315
$T_{sk}^{\text{set}}$ [°C] 33.7 0.169
$T_{c}^{\text{set}}$ [°C] 36.8 0.184

Parameters considered as variables in Pierce's model

The resulting SET differences are presented in the following sections.

### 5.2 Gender effect

The difference in comfort temperature between women and men has been studied experimentally and numerically[Kingma, 2015; Hoof, 2008]. Here, the average SET for male and female individuals is calculated from random draws on parameters according to Gaussian laws described in section 5.1. Basal metabolism is calculated according to the very first equation based on the size and mass of each individual and the level of clothing chosen is 1[clo].

The average SET difference between the two genders is shown in Fig. 14 for operating conditions, i. e. $T_a = T_r$at a low air velocity of $v=0.1$ [m/s]. The difference between average SET ranges from -4 to -1.25 [K] for the conditions studied.

In the range between 19 and 25°C of operating temperature which contains the comfort zone defined according to[Nishi, 1977] - see Fig. 8 - the average gender difference varies from -3.5 to -2 [K]. The average SET calculated for men and women therefore differ significantly and in some cases this difference is larger than the comfort zone, which is 3.4 [K].

It should also be noted that according to the results of the model, the SET obtained for women is higher than that of men: this implies that a higher operating temperature of the reference environment is required to produce the same physiological reaction as in men.

### 5.3 Environment effect

In this section, the effect of the variability of a population with both sexes equiprobably distributed is presented. Physiological parameters are distributed according to Gaussian distributions as described in 5.1. The clothing is 1[clo] and the activity is 1[met].

The standard deviation of the SET distributions obtained at any point in the psychrometric diagram is plotted for high wind conditions ($v=1$ [m/s]) on Fig. 15. It can be observed that the standard deviation can vary from $\sim 0.4$ à $2$ [K]. This implies that under these conditions, the temperature felt by 95% of the population varies between 1.6 and 8 [K], or four standard deviations.

The standard deviation obtained for high solar flux conditions, i. e. $T_r=T_a + 30$and a low air velocity $v=0.1$ [m/s] are shown in Fig. 16. We observe that it varies between $\sim 0.2$ and 1.4 [K], which implies that 95% of the population is in a range between 0.8 and 5.6 [K] under conditions of high radiant temperature.

While these ranges may seem extensive, the literature review of [Hoof, 2008] mentions experimental campaigns where deviations of $\sim 10$ [K] were recorded between temperatures declared comfortable. In indoor environments, it is therefore difficult to define a setpoint temperature adapted to each individual. In outdoor space, it is appropriate to integrate variability into the estimation of comfort, in the spirit of [Bruse, 2007].

## 6. Conclusion

This article exposes the fact that the determination of comfort in environments with variable conditions requires an appropriate indicator, allowing the calculation of metabolic reactions in transient conditions.

Clothing also plays a significant role in the heat and moisture balance of the human body. Its transfer properties being highly variable with wind speed, it is necessary to obtain an appropriate estimate of air velocity around the body in (semi-) outdoor spaces.

Morphological characteristics such as gender, mass, height also have a significant influence on the numerical evaluation of comfort by the SET model. A significant variation in perceived comfort can be introduced by physiological variability, especially in outdoor conditions where high air velocities and strong solar radiation are often encountered.

In this article, a dynamic model of the human body for the calculation of comfort has been presented. Other models exist and can be used, including the Physiological Equivalent Temperature (PET). The latter uses many of the equations of the two-nodes model presented here, as well as a similar reference environment. However, PET is defined as the temperature of an environment that would produce the same skin and core temperatures as the environment being studied.

Other more complex models have been developed, ranging up to 340 nodes [Fiala, 1999] leading to the comfort index named Universal Thermal Comfort Index (UTCI). However, the complexity of the model does not make it of the practical nor provides an easy calibration of physiological information, particularly if variability is to be incorporated. A correlation based on multiple linear regressions has nevertheless been derived and allows us to evaluate the comfort of a pedestrian walking at a speed of 4[km/h].

The individual's thermal history also has an influence on his or her feelings: for example, one appreciates an atmosphere at 16°C for a few minutes after making an intense effort, while the same atmosphere would be considered too cold when the body is in thermal balance. The approach to alliesthesia developed by [Parkinson, 2012] proposes a model for this phenomenon. The notion of "pleasure" felt during alliesthesic phenomena also goes beyond the usual approach, which focuses on thermal neutrality.

The psychological aspects of comfort should also not be neglected: it has been proven that under identical conditions, the feeling of comfort is better in a furnished room (furniture, decoration) than in a bare room [Rohles, 2007]. Similarly, the placebo effect has an impact: when the user thinks he or she has control over the set temperature, his or her feeling is positively affected, whether or not this control is effective [Rohles, 2007].

In this sector, local comfort features such as heaters, heated seats, office fans or seat-mounted fans provide an interesting response to the variability in comfort temperatures. Considerable energy savings have been achieved at a pilot site [Brager, 2015]: by providing each individual with a power of around 20 Watts, the decrease in the overall setpoint temperature has resulted in a power reduction of 500 W/person.

The source code of the SET and PET comfort indexes can be dowloaded on Sharing area in this site