Occupant Thermal Comfort[LINK]
The integration of a sophisticated building thermal analysis tool with thermal comfort models allows one to perform an energy analysis on a zone and simultaneously determine if the environmental control strategy will be sufficient for the occupants to be thermally comfortable. This section is intended to provide background on thermal comfort, present an overview of state of the art thermal comfort models, and present the mathematical models that have been incorporated into EnergyPlus.
Thermal comfort modeling in EnergyPlus is controlled primarily by the People input object. This includes input for selecting the type of thermal comfort model that is desired by the user as well as parameters that serve as inputs to the thermal comfort models. These parameters include the activity level, the work efficiency, the air velocity and the clothing insulation level for people within the space. All four of these parameters can be scheduled. More information on the People input object can be found in the EnergyPlus Input/Output Reference. More information on how each of these parameters is used and the specific modeling equations for the thermal comfort models can be found below.
Background on Thermal Comfort Models[LINK]
Throughout the last few decades, researchers have been exploring the thermal, physiological and psychological response of people in their environment in order to develop mathematical models to predict these responses. Researchers have empirically determined building occupants’ thermal responses to the combined thermal effect of the personal, environmental and physiological variables that influence the condition of thermal comfort.
There are two personal variables that influence the condition of thermal comfort: the thermal resistance of the clothing (Icl), and the metabolic rate (H/ADu).
The thermal resistance of the clothing (Icl) is measured in units of “clo.” The 1985 ASHRAE Handbook of Fundamentals (ASHRAE 1985) suggests multiplying the summation of the individual clothing items clo value by a factor of 0.82 to determine the overall clo value for clothing ensembles.
The metabolic rate (H/ADu) is a measure of the internal heat production rate of an occupant (H) in W/hr per unit of “Dubois” body surface area (ADu) in units of m2. The DuBois body surface area is given by :
ADu=0.202(weight)0.425(height)0.725
Using this equation, an area of 1.8 m2 represents the surface area of an average person of weight 70 kg. and height 1.73 m (Fanger 1967). In EnergyPlus, this average person area of 1.8 m2 is used for the body surface area all thermal comfort models. The metabolic rate is measured in mets, where 1 met = 58.2 W/m2.
The environmental variables that influence the conditions of thermal comfort include:
Air Temperature (Ta),
Mean Radiant Temperature (Tr),
Relative air velocity (v),
Water vapor pressure in ambient air (Pa)
The Air Temperature (Ta), a direct environmental index, is the drybulb temperature of the environment. It is calculated using either the mean air temperature for the zone in question or the air temperature calculated by some of the more detailed zone air models available in EnergyPlus. The Mean Radiant Temperature (Tr) is a rationally derived environmental index defined as the uniform blackbody temperature that would result in the same radiant energy exchange as in the actual environment. This value can be calculated in several ways as outlined in the Mean Radiant Temperature Calculation section later in this document. The Relative air velocity (v) a direct environmental index is a measure of the air motion obtainable via a hot wire or vane anemometers. The Water vapor pressure in ambient air (Pa) is a direct environmental index.
The physiological variables that influence the conditions of thermal comfort include:
Skin Temperature (Tsk),
Core or Internal Temperature (Tcr),
Sweat Rate,
Skin Wettedness (w),
Thermal Conductance (K) between the core and skin.
where the Skin Temperature (Tsk), the Core Temperature (Tcr) and the Sweat Rate are physiological indices. The Skin Wettedness (w) is a rationally derived physiological index defined as the ratio of the actual sweating rate to the maximum rate of sweating that would occur if the skin were completely wet.
One more consideration is important in dealing with thermal comfort  the effect of asymmetrical heating or cooling. This could occur when there is a draft or when there is a radiant flux incident on a person (which is what is of primary interest to us here). Fanger (1967) noted that the human regulatory system is quite tolerant of asymmetrical radiant flux. A reasonable upper limit on the difference in mean radiant temperature (Tr) from one direction to the opposing direction is 15°C. (ASHRAE 1984). This limit is lower if there is a high air velocity in the zone.
General Nomenclature list for Thermal Comfort Models
ADu

Dubois body surface area 
m2

 
 
H 
Internal heat production rate of an occupant per unit area = M – W 
W/m2

 
IntHeatProd 
Icl

Thermal resistance of the clothing 
clo 
 
CloUnit 
M 
Metabolic rate per unit area 
W/m2

 
ActLevel 
Pa

Water vapor pressure in ambient air 
Torr 
 
VapPress 
Ta

Air temperature 
°C 
 
AirTemp 
Tcr

Core or internal temperature 
°C 
 
CoreTemp 
Tr

Mean radiant temperature 
°C 
 
RadTemp 
Tsk

Skin temperature 
°C 
 
 
v 
Relative air velocity 
m/s 
 
AirVel 
W 
The rate of heat loss due to the performance of work 
W/m2

 
WorkEff 
w 
Skin wettedness 
 
 
 
Mathematical Models for Predicting Thermal Comfort[LINK]
Many researchers have been exploring ways to predict the thermal sensation of people in their environment based on the personal, environmental and physiological variables that influence thermal comfort. From the research done, some mathematical models that simulate occupants’ thermal response to their environment have been developed. Most thermal comfort prediction models use a seven or nine point thermal sensation scale, as in the following tables.
Seven point Thermal Sensation Scale
3 
Hot 
2 
Warm 
1 
slightly warm 
0 
neutral 
1 
slightly cool 
2 
cool 
3 
cold 
Nine point Thermal Sensation Scale
4 
very hot 
3 
hot 
2 
warm 
1 
slightly warm 
0 
neutral 
1 
slightly cool 
2 
cool 
3 
cold 
4 
very cold 
The most notable models have been developed by P.O. Fanger (the Fanger Comfort Model), the J. B. Pierce Foundation (the Pierce TwoNode Model), and researchers at Kansas State University (the KSU TwoNode Model). Berglund (1978) presents a detailed description of the theory behind these three models.
Note for all Thermal Comfort reporting: Though the published values for thermal comfort “vote” have a discrete scale (e.g. –3 to +3 or –4 to +4), the calculations in EnergyPlus are carried out on a continuous scale and, thus, reporting may be “off the scale” with specific conditions encountered in the space. This is not necessarily an error in EnergyPlus – rather a different approach that does not take the “limits” of the discrete scale values into account.
The main similarity of the three models is that all three apply an energy balance to a person and use the energy exchange mechanisms along with experimentally derived physiological parameters to predict the thermal sensation and the physiological response of a person due to their environment. The models differ somewhat in the physiological models that represent the human passive system (heat transfer through and from the body) and the human control system (the neural control of shivering, sweating and skin blood flow). The models also differ in the criteria used to predict thermal sensation. However, all three models use information from the People statement and the thermal comfort model is selected via the People statement in a user’s input file. Scheduled parameters such as the activity level, work efficiency, air velocity, and clothing insulation level all have a direct bearing on the thermal comfort models. For more information on the input of these parameters, see the People statement in the EnergyPlus Input/Output Reference. For more information on how each individual thermal comfort model uses these parameters, please consult the next several sections.
Fanger Comfort Model[LINK]
Fanger’s Comfort model was the first human comfort model that was developed. It was first published in 1967 (Fanger 1967), then updated in 1970 (Fanger 1970), and helped set the stage for the Pierce and KSU models. The mathematical model developed by P.O. Fanger is probably the most well known of the three models and is the easiest to use because it has been put in both chart and graph form.
Fanger Model Nomenclature List[LINK]
Nomenclature list for Fanger model
ADu

Dubois body surface area 
m2

 
BodySurfaceArea 
Cres

The rate of dry respiratory heat loss 
W/m2

 
DryRespHeatLoss 
Edif

The rate of heat loss from the diffusion of water vapor through the skin 
W/m2

 
EvapHeatLossDiff 
Eres

The rate of latent respiratory heat loss 
W/m2

 
LatRespHeatLoss 
Ersw,req

The rate of heat loss from the evaporation of regulatory sweating at the state of comfort 
W/m2

 
EvapHeatLossRegComf 
Esk

Total evaporative heat loss from skin 
W/m2


EvapHeatLoss 
fcl

The ratio of clothed body 
 

CloBodyRat 
feff

The fraction of surface effective for radiation ( = 0.72) 
 
 
RadSurfEff 
H 
Internal heat production rate of an occupant per unit area ( = M – W) 
W/m2

 
IntHeatProd 
hc

Convective heat transfer coefficient 
W/m2°C 
 
Hc 
L 
All the modes of energy loss from body 
W/m2

 
 
M 
Metabolic rate per unit area 
W/m2

 
ActLevel 
Pa

Water vapor pressure in ambient air 
Torr 
 
VapPress 
PMV 
Predicted Mean Vote 
 
(4)4 
PMV 
PPD 
Predicted Percentage of Dissatisfied 
 
0100 
PPD 
Qc

The rate of convective heat loss 
W/m2

 
ConvHeatLoss 
Qdry

Sensible heat flow from skin 
W/m2


DryHeatLoss 
Qr

The rate of radiative heat loss 
W/m2

 
RadHeatLoss 
Qres

The rate of respiratory heat loss 
W/m2

 
RespHeatLoss 
Ta

Air temperature 
°C 
 
AirTemp 
Tcl

Clothing surface temperature 
°C 
 
CloSurfTemp 
Tcla

Clothing surface temperature (Absolute) 
°K 
 
AbsCloSurfTemp 
Tra

Mean radiant temperature 
°K 
 
AbsRadTemp 
W 
The rate of heat loss due to the performance of work 
W/m2

 
WorkEff 
ε 
The emissivity of clothingskin surface 
 
 
SkinEmiss 
σ 
The StefanBoltzman constant ( = 5.67×10 ) 
W/m2K4

 
StefanBoltz 
Description of the model and algorithm[LINK]
Fanger developed the model based on the research he performed at Kansas State University and the Technical University of Denmark. Fanger used the sevenpoint form of a thermal sensation scale along with numerous experiments involving human subjects in various environments. He related the subjects’ response to various environmental variables which influence the condition of thermal comfort. Fanger’s model is based upon an energy analysis that takes into account all the modes of energy loss (L) from the body, including: the convection and radiant heat loss from the outer surface of the clothing, the heat loss by water vapor diffusion through the skin, the heat loss by evaporation of sweat from the skin surface, the latent and dry respiration heat loss and the heat transfer from the skin to the outer surface of the clothing. The model assumes that the person is thermally at steady state with the interior environment.
M=L
L=Qres+Qdry+Esk
Qres=Eres+Cres=0.000017M(5867−Pa)+0.0014M(34−Ta)
In the EnergyPlus source code, Qres, Eres and Cres are calculated using the following equations, respectively:
RespHeatLoss = LatRespHeatLoss + DryRespHeatLoss
LatRespHeatLoss = 0.000017*ActLevel*(5867.  VapPress)
DryRespHeatLoss = 0.0014*ActLevel*(34. AirTemp)
Qr=hr×fcl(Tcl−Tr)
hc=8.6×v0.53
Ersw=0.42(H−58.2)for(H>58.2)
In the EnergyPlus source code, Qdry, Qc and Qr are calculated using the following equations, respectively:
DryHeatLoss = ConvHeatLoss + RadHeatLoss
ConvHeatLos = CloBodyRat*Hc*(CloSurfTemp  AirTemp)
RadHeatLoss = 3.96*CloBodyRat*[(0.01*AbsCloSurfTemp)**4  (0.01*AbsRadTemp)**4]
Ersw=0for(H≤58.2)
Esk=Ersw+Ediff
PPD=100.0−95.0∗EXP(−0.03353PMV4−0.2179PMV2)
Ediff=wdiff×Emax
In the EnergyPlus source code, Ersw, Ef and Esk are calculated using the following equations, respectively:
EvapHeatLossRegComf = 0.42*(IntHeatProd  ActLevelConv) or 0.0 depending on H value
EvapHeatLossDiff = 0.00305*(5733.0  6.99*IntHeatProd  VapPress)
EvapHeatLoss = EvapHeatLossRegComf + EvapHeatLossDiff
The Fanger model was developed using the thermal sensation votes from subjects at KSU and in Denmark, and the Predicted Mean Vote (PMV) thermal sensation scale is based on how the energy loss (L) deviates from the metabolic rate (M) in the following form:
PMV=(0.303e−0.036M+0.028)(H−L)
In the EnergyPlus source code, PMV is calculated using the following equations in conjunction with the equations shown above:
PMV = ThermSensTransCoef*(IntHeatProd  EvapHeatLoss  RespHeatLoss  DryHeatLoss)
ThermSensTransCoef = 0.303*EXP(0.036*ActLevel) + 0.028
IntHeatProd = ActLevel  WorkEff
Predicted Percent of Dissatisfied (PPD) people at each PMV is calculated as follows:
SIGsk=Tsk−33.7
Pierce TwoNode Model[LINK]
The Pierce TwoNode model was developed at the John B. Pierce Foundation at Yale University. The model has been continually expanding since its first publication in 1970 (Gagge et.al. 1970). The most recent version on the model appears in the 1986 ASHRAE Transactions (Gagge et.al. 1986).
Pierce TwoNode Model Nomenclature List[LINK]
Nomenclature list for Pierce TwoNode model
Cdil

Constant for skin blood flow 


SkinBloodFlowConst 
Cres

The rate of dry respiratory heat loss 
W/m2

 
DryRespHeatLoss 
Csw

Proportionality constant for sweat control 
g/m2hr 

SweatContConst 
DISC 
Predicted discomfort vote 
 
(5)5 
DISC 
Edif

The rate of heat loss from the diffusion of water vapor through the skin 
W/m2

 
EvapHeatLossDiff 
Emax

Maximum evaporative heat loss 
W/m2


EvapHeatLossMax 
Esk

Total evaporative heat loss from skin 
W/m2


EvapHeatLoss 
Eres

The rate of latent respiratory heat loss 
W/m2

 
LatRespHeatLoss 
Ersw

The rate of heat loss from the evaporation of regulatory sweating 
W/m2

 
EvapHeatLossRegSweat 
Ersw,req

The rate of heat loss from the evaporation of regulatory sweating at the state of comfort 
W/m2


EvapHeatLossRegComf 
ET* 
Effective Temperature 
°C 
 
ET 
fcl

The ratio of clothed body 
 

CloBodyRat 
feff

The fraction of surface effective for radiation ( = 0.72) 
 
 
RadSurfEff 
H 
Internal heat production rate of an occupant per unit area ( = M – W) 
W/m2

 
IntHeatProd 
h 
Combined heat transfer coefficient 
W/m2°C 

H 
hc

Convective heat transfer coefficient 
W/m2°C 
 
Hc 
he’ 
Combined evaporative heat transfer coefficient 
W/(m2kPa) 

 
hr

Radiant heat transfer coefficient 
W/m2°C 
 
Hr 
Icl

Clothing insulation 
m2°C/W 

 
L 
All the modes of energy loss from body 
W/m2

 
 
LET∗

All the modes of energy loss from body at ET* 
W/m2


 
LSET∗

All the modes of energy loss from body at SET* 
W/m2


 
M 
Metabolic rate per unit area 
W/m2

 
ActLevel 
Mact

Metabolic heat production due to activity 
W/m2


 
Mshiv

Metabolic heat production due to shivering 
W/m2


ShivResponse 
Pa

Water vapor pressure in ambient air 
Torr 
 
VapPress 
PMV* 
Predicted Mean Vote modified by ET* or SET* 
 
(4)4 
PMVET PMVSET 
PSET∗

Water vapor pressure at SET* 
°C 

StdVapPressSET 
Psk

Saturated water vapor pressure at required skin temperature 
Torr 
 
SatSkinVapPress 
Qc

The rate of convective heat loss 
W/m2

 
ConvHeatLoss 
Qcrsk

Heat flow from core to skin 
W/m2


HeatFlow 
Qdry

Sensible heat flow from skin 
W/m2


DryHeatLoss 
Qr

The rate of radiative heat loss 
W/m2

 
RadHeatLoss 
Qres

The rate of respiratory heat loss 
W/m2

 
RespHeatLoss 
Scr

Heat storage in core compartment 
W/m2


CoreheatStorage 
SET* 
Standard Effective Temperature 
°C 
 
SET 
SIGb

Thermal signal of body 
°C 

BodyThermSigCold BodyThermSigWarm 
SIGcr

Thermal signal of core 
°C 

CoreThermSigCold CoreThermSigWarm 
SIGsk

Thermal signal of skin 
°C 

SkinThermSigCold SkinThermSigWarm 
SKBF 
Skin blood flow 
L/m hr 

SkinBloodFlow 
Ssk

Heat storage in skin compartment 
W/m2


SkinHeatStorage 
Str

Constriction constant of skin blood flow for average person 


Str 
SWreg

The rate of regulatory sweating 
g/m2hr 

RegSweat 
Ta

Air temperature 
°C 
 
AirTemp 
Tb

Mean body temperature 


AvgBodyTemp 
Tb−c

Mean body temperature when DISC is zero (lower limit) 
°C 

AvgBodyTempLow 
Tb−h

Mean body temperature when HSI is 100 (upper limit) 
°C 

AvgBodyTempHigh 
Tcl

Clothing surface temperature 
°C 
 
CloSurfTemp 
Tcr

Core or internal temperature 
°C 
 
CoreTemp 
Tr

Mean radiant temperature 
°C 
 
RadTemp 
TSENS 
Thermal sensation vote 
 
(5)5 
TSENS 
Tsk

Skin temperature 
°C 

SkinTemp 
W 
The rate of heat loss due to the performance of work 
W/m2

 
WorkEff 
wdif

Skin wettedness due to diffusion trough the skin 


SkinWetDiff 
wrsw

Skin wettedness due to regulatory sweating 


SkinWetSweat 
ε 
The emissivity of clothingskin surface 
 
 
SkinEmiss 
σ 
The StefanBoltzman constant ( = 5.67×10 ) 
W/m2K4

 
StefanBoltz 
Description of the model and algorithm[LINK]
The Pierce model thermally lumps the human body as two isothermal, concentric compartments, one representing the internal section or core (where all the metabolic heat is assumed to be generated) and the other representing the skin. This allows the passive heat conduction from the core compartment to the skin to be taken into account. The boundary line between the two compartments changes with the skin blood flow rate per unit skin surface area (SKBF in L/h•m2) and is described by alpha – the fraction of total body mass attributed to the skin compartment (Doherty and Arens 1988).
α=0.0417737+0.7451832/(SKBF+0.585417)
In the EnergyPlus code, α or SkinMassRat is calculated using the equivalent equation:
SkinMassRat = 0.0417737 + 0.7451832/(SkinBloodFlow + 0.585417)
Furthermore, the model takes into account the deviations of the core, skin, and mean body temperature weighted by alpha from their respective setpoints. Thermoregulatory effector mechanisms (regulatory sweating, skin blood flow, and shivering) are defined in terms of thermal signals from the core, skin and body in ∘C (Doherty and Arens 1988) by the following equations:
SIGb=Tb−36.49
SKBF=(6.3+Cdil×SIGcr)/(1+Str×(SIGsk))
SWreg=CswSIGbe(SIGsk/10.7)
In the EnergyPlus code, when SkinTemp is greater than SkinTempSet (33.7∘C), then SkinThermSigWarm = SkinTemp  SkinTempSet and SkinThermSigCold = 0.0. Otherwise, when SkinTemp is less than or equal to SkinTempSet, then SkinThermSigWarm = 0.0 and SkinThermSigCold = SkinTempSet  SkinTemp.
In the EnergyPlus code, when CoreTemp is greater than CoreTempSet (36.8∘C), then CoreThermSigWarm = CoreTemp  CoreTempSet and CoreThermSigCold = 0.0. Otherwise, when CoreTemp is less than or equal to CoreTempSet, then CoreThermSigWarm = 0.0 and CoreThermSigCold = CoreTempSet  CoreTemp.
In the EnergyPlus code, when AvgBodyTemp is greater than AvgBodyTempSet(36.49∘C), then BodyThermSigWarm = AvgBodyTemp  AvgBodyTempSet and BodyThermSigCold = 0.0. Otherwise, when AvgBodyTemp is less than or equal to AvgBodyTempSet, then BodyThermSigWarm = 0.0 and BodyThermSigCold = AvgBodyTempSetAvgBodyTemp.
The skin blood flow rate in L/hr⋅m2 is calculated using:
Mshiv=19.4(SIGcr)(SIGsk)
In the EnergyPlus code, this equation is calculated using the following equations:
SkinBloodFlow = (6.3 + VasodilationFac)/(1. + VasoconstrictFac)
VasodilationFac = SkinBloodFlowConst*CoreWarmDelTemp
VasoconstrictFac = Str*SkinColdDelTemp
The rate of regulatory sweating in g/hr•m2 is calculated using:
M=Mact+Mshiv
and the shivering response in W/m2 is calculated using:
L=Qres+Qdry+Esk+W
In the EnergyPlus code, the sweat and shivering response determined by the above two equations are calculated using the follow expressions:
RegSweat = SweatContConst*BodyThermSigWarm*EXP(SkinThermSigWarm/10.7)
ShivResponse = 19.4*SkinThermSigCold*CoreThermSigCold
The latest version of the Pierce model (Fountain and Huizenga 1997) discusses the concepts of SET* and ET*. The Pierce model converts the actual environment into a “standard environment” at a Standard Effective Temperature, SET*. SET* is the drybulb temperature of a hypothetical environment at 50% relative humidity for subjects wearing clothing that would be standard for the given activity in the real environment. Furthermore, in this standard environment, the same physiological strain, i.e. the same skin temperature, skin wettedness and heat loss to the environment, would exist as in the real environment. The Pierce model also converts the actual environment into a environment at an Effective Temperature, ET*, that is the drybulb temperature of a hypothetical environment at 50% relative humidity and uniform temperature (Ta = MRT) where the subjects would experience the same physiological strain as in the real environment.
In the latest version of the model, it is suggested that the classical Fanged PMV be modified by using ET* or SET* instead of the operative temperature. This gives a new index PMV* which is proposed for dry or humid environments. It is also suggested that PMV* is very responsive to the changes in vapor permeation efficiency of the occupants clothing.
Qres=Eres+Cres=0.0023M(44−Pa(torr))+0.0014M(34−Ta)Qres=0.017251M(5.8662−Pa(kPa))+0.0014M(34−Ta)
which in the EnergyPlus code takes on the form ActLevel = ActLevel + ActShiv.
Qdry=Qc+Qr
Qc=hc×fcl(Tcl−Ta)
In the EnergyPlus code, the calculation of Qres or RespHeatLoss is accomplished using the following equations:
RespHeatLoss = LatRespHeatLoss + DryRespHeatLoss
LatRespHeatLoss = 0.017251*ActLevel*(5.8662  VapPress)
DryRespHeatLoss = 0.0014*ActLevel*(34. AirTemp)
Qr=hr×fcl(Tcl−Tr)
hc=8.6×v0.53
hc=5.66(M/58.2−0.85)0.39
Qdry or DryHeatLoss is calculated in the EnergyPlus code using the following equation:
DryHeatLoss = CloBodyRat*(Hc*(CloSurfTemp  AirTemp) + Hr*(CloSurfTemp  RadTemp))
In Pierce model, the convective heat transfer coefficient, hc, varies with the air velocity around body and metabolic rate. The model uses the maximum value of following equations:
hr=4.×feffεσ((Tcl+Tr)/2+273.15)3
Tcl=(Tsk/Icl+fcl(hcTa+hrTr))/(1/Icl+fcl(hc+hr))
In the EnergyPlus code, the following equations are used to calculate hc or Hc:
Hc = 8.6*AirVel**0.53 where AirVel is reset to 0.137 if it drops below this value
HcAct = 5.66*(ActMet  0.85)**0.39 when ActMet is greater than 0.9
When ActMet is greater than 0.9, EnergyPlus uses the larger value of Hc and HcAct for hc.
Also, in the model, the radiant heat transfer coefficient, hr, is defined by following equation (Doherty and Arens 1988):
Esk=Ersw+Ediff
In the EnergyPlus code, hr or Hr is calculated using the following equation:
Hr = 4.*RadSurfEff*StefanBoltz*((CloSurfTemp + RadTemp)/2. + TAbsConv)**3
In the Pierce model, Tcl in ∘C is estimated by each iteration using following equation:
Ersw=0.68×SWreg
Or as shown in the EnergyPlus code:
CloSurfTemp = (CloCond*SkinTemp + CloBodyRat*(Hc*AirTemp+Hr*RadTemp))
/(CloCond + CloBodyRat*(Hc+Hr))
The total evaporative heat loss from the skin, Esk, includes evaporation of water produced by regulatory sweating, Ersw, and evaporation of water vapor that diffuses through the skin surface, Ediff.
Ediff=wdiff×Emax
Emax=he′(Psk−Pa)
Qcrsk=(5.28+1.163SKBF)(Tcr−Tsk)
Within the EnergyPlus source code, these equations are shown as:
EvapHeatLoss = EvapHeatLossRegSweat + EvapHeatLossRegDiff
RegHeatLoss = 0.68*RegSweat
DiffHeatLoss = SkinWetDiff*MaxEvapHeatLoss
where 0.68 is the passive water vapor diffusion rate in g/h·m2·Torr and
wdiff=0.06(1−wrsw)
Ssk=Qc−s−Qc−Qr−Esk
wrsw=Ersw/Emax
These equations within the EnergyPlus source code appear as:
SkinWetDiff = .06*(1.SkinWetSweat)
MaxEvapHeatLoss = (1./TotEvapHeatResist)*(SatSkinVapPress  VapPress)
SkinWetSweat = EvapHeatLossRegSweat/MaxEvapHeatLoss
The Pierce model has one additional heat flow term describing the heat transfer between the internal core compartment and the outer skin shell (Doherty and Arens 1988).
Scr=M−W−Qres−Qc−s
or as shown in the EnergyPlus source code:
HeatFlow = (CoreTempSkinTemp)*(5.28 + 1.163*SkinBloodFlow)
where 5.28 is the average body tissue conductance in W/m$^2 $ K and 1.163 is the thermal capacity of blood in W ⋅ h/L ⋅ K.
Thus, individual heat balance equations for core and skin compartments are expressed using this term, Qc−s. New temperatures of core, skin and body are calculated for each iteration from rates of heat storage in the core and skin.
Tb−c=(0.185/58.2)(M−W)+36.313
Tb−h=(0.359/58.2)(M−W)+36.664
These equations are represented by the following source code lines in EnergyPlus:
SkinHeatStorage = HeatFlow  DryHeatLoss  EvapHeatLoss
CoreHeatStorage = IntHeatProd  RespHeatLoss  HeatFlow
As a result, PMV evaluated for ET* and SET* can be calculated using:
PMVET=(0.303e−0.036M+0.028)(H−LET∗)
PMVSET=(0.303e−0.036M+0.028)(H−LSET∗)
These equations are represented by the following lines of source code in EnergyPlus:
ThermSensTransCoef = 0.303*EXP(0.036*ActLevel) + 0.028
PMVET = ThermSensTransCoef*(IntHeatProd  EvapHeatLossDiff  EvapHeatLossRegComf
 RespHeatLoss  DryHeatLossET)
PMVSET = ThermSensTransCoef*(IntHeatProd  EvapHeatLossDiff
 EvapRegHeatLossRegComf  RespHeatLoss  DryHeatLossSET)
Besides PMV*, the Pierce Two Node Model uses the indices TSENS and DISC as predictors of thermal comfort. Where TSENS is the classical index used by the Pierce foundation and is a function of the mean body temperature. DISC is defined as the relative thermoregulatory strain that is needed to bring about a state of comfort and thermal equilibrium. DISC is a function of the heat stress and heat strain in hot environments and equal to TSENS in cold environments. In summary, the Pierce Model, for our purposes, uses four thermal comfort indices; PMVETa function of ET*, PMVSET a function of SET*, TSENS and DISC.
TSENSc=0.68175(Tb−Tb−c)
TSENSh=4.7(Tb−Tb−c)/(Tb−h−Tb−c)
When Tb≤Tb−c, TSENS is calculated using:
WcrCcrdTcr/dt=M−W−Qres−KS(Tcr−Tsk)
Otherwise, when Tb>Tb−c, TSENS is calculated using:
WskCskdTsk/dt=KS(Tcr−Tsk)−Qdry−Esk
Finally, DISC is calculated using the following equation:
DISC=5.(Ersw−Ersw−comf)/(Emax−Ersw−comf−Ediff)
These final equations for the Pierce Two Node Model are calculated using the following lines of source code:
AvgBodyTempLow = (0.185/ActLevelConv)*IntHeatProd + 36.313
AvgBodyTempHigh = (0.359/ActLevelConv)*IntHeatProd + 36.664
TSENS = .68175*(AvgBodyTempAvgBodyTempLow)
TSENS = 4.7*(AvgBodyTempAvgBodyTempLow)/(AvgBodyTempHighAvgBodyTempLow)
DISC = 5.*(EvapHeatLossRegSweatEvapHeatLossRegComf)
/(MaxEvapHeatLossEvapHeatLossRegComfDiffHeatLoss)
KSU TwoNode Model[LINK]
The KSU TwoNode Model, developed at Kansas State University, was published in 1977 (Azer and Hsu 1977). The KSU model is quite similar to that of the Pierce Foundation. The main difference between the two models is that the KSU model predicts thermal sensation (TSV) differently for warm and cold environment.
KSU Two Node Model Nomenclature List[LINK]
Nomenclature list for KSU TwoNode model
Ccr

Specific heat of body core 
Whr/kg°C 


Csk

Specific heat of skin 
Whr/kg°C 


Cres

The rate of dry respiratory heat loss 
W/m2

 
DryRespHeatLoss 
Edif

The rate of heat loss from the diffusion of water vapor through the skin 
W/m2

 
EvapHeatLossDiff 
Emax

Maximum evaporative heat loss 
W/m2


EvapHeatLossMax 
Esk

Total evaporative heat loss from skin 
W/m2


EvapHeatLoss 
Esw

Equivalent evaporation heat loss from the sweat secreted 
W/m2


EvapHeatLossSweat 
Esw.d

Sweat function for warm and dry skin 
W/m2


EvapHeatLossDrySweat 
Eres

The rate of latent respiratory heat loss 
W/m2

 
LatRespHeatLoss 
Fcl

The Burton thermal efficiency factor for clothing 

 
CloThermEff 
Fpcl

Permeation efficiency factor for clothing 

 
CloPermeatEff 
H 
Internal heat production rate of an occupant per unit area = M  W 
W/m2

 
IntHeatProd 
H 
Combined heat transfer coefficient 
W/m2°C 

H 
hc

Convective heat transfer coefficient 
W/m2°C 
 
Hc 
hr

Radiant heat transfer coefficient 
W/m2°C 
 
Hr 
KS 
Overall skin thermal conductance 
W/m2°C 

ThermCndct 
KSo

Skin conductance at thermal neutrality 
W/m2°C 

ThermCndctNeut 
KS−4

Skin conductance at thermal sensation very cold 
W/m2°C 

ThermCndctMin 
M 
Metabolic rate per unit area 
W/m2

 
ActLevel 
Mshiv

Metabolic heat production due to shivering 
W/m2


ShivResponse 
Pa

Water vapor pressure in ambient air 
Torr 
 
VapPress 
Psk

Saturated water vapor pressure at required skin temperature 
Torr 
 
SatSkinVapPress 
PTaccl

The pattern of acclimation 


AcclPattern 
Qc

The rate of convective heat loss 
W/m2

 
ConvHeatLoss 
Qdry

Sensible heat flow from skin 
W/m2


DryHeatLoss 
Qr

The rate of radiative heat loss 
W/m2

 
RadHeatLoss 
Qres

The rate of respiratory heat loss 
W/m2

 
RespHeatLoss 
RH 
Relative humidity 


RelHum 
Ta

Air temperature 
°C 
 
AirTemp 
Tcr

Core or internal temperature 
°C 
 
CoreTemp 
To

Operative temperature 
°C 
 
OpTemp 
Tr

Mean radiant temperature 
°C 
 
RadTemp 
Tsk

Skin temperature 
°C 

SkinTemp 
TSV 
Thermal sensation vote 

(4)4 
TSV 
V 
Relative air velocity 
m/s 
 
AirVel 
W 
The rate of heat loss due to the performance of work 
W/m2

 
WorkEff 
w 
Skin wettedness 
 
 
SkinWetTot 
Wcr

Mass of body core per unit body surface 
kg/m2


 
wrsw

Skin wettedness due to regulatory sweating 


SkinWetSweat 
wrsw−o

Skin wettedness at thermal neutrality 


SkinWetSweatNeut 
Wsk

Mass of skin per unit body surface 
kg/m2


 
Description of the model and algorithm[LINK]
The KSU TwoNode Model is based on the changes that occur in the thermal conductance between the core and the skin temperature in cold environments. In warm environments, the thermal conductance of the model is based on changes in the skin wettedness.
In this model, metabolic heat production is generated in the core which exchanges energy with the environment by respiration, and the skin exchanges energy with the interior environment by convection and radiation. In addition, body heat is dissipated through the evaporation of sweat and/or water vapor diffusion through the skin. These principles are outlined in the following system of equations.
Qres=Eres+Cres=0.0023M(44−Pa(torr))+0.0014M(34−Ta)
Qdry=Qc+Qr=hfclFcl(Tsk−To)
where Qres and Qdry are calculated using:
h=hc+hr
hc=8.3√v
The calculation of Qres (RespHeatLoss) and Qdry (DryHeatLoss) in the source code are accomplished through the following equations:
RespHeatLoss = LatRespHeatLoss + DryRespHeatLoss
LatRespHeatLoss = 0.0023*ActLevelTot*(44.  VapPress)
DryRespHeatLoss = 0.0014*ActLevelTot*(34.  AirTemp)
DryHeatLoss = H*CloBodyRat*CloThermEff*(SkinTemp  OpTemp)
In the equation for Qdry, h and To are calculated using:
hr=ϵavg(3.87+0.031Tr)
To=(hcTa+hrTr)/(hc+hr)
For(Esw≤Emax),Esk=Esw+(1−wrsw)Ediff
For(Esw>Emax),Esk=Emax
These equations are represented by the following source code equations:
H = Hc + Hr
Hc = 8.3*SQRT(AirVel) where AirVel is reset to 0.137 if it drops below this value
Hr = EmissAvg*(3.87 + 0.031*RadTemp)
OpTemp = (Hc*AirTemp + Hr*RadTemp)/H
In the KSU TwoNode Model, evaporative heat losses are defined by these equations:
Ediff=0.408(Psk−Pa)
Emax=2.2hcFpcl(Psk−Pa)
KS=5.3+[6.75+42.45(Tcr−36.98)+8.15(Tcr−35.15)0.8(Tsk−33.8)]1+0.4(32.1−Tsk)
Mshiv=20(36.9−Tcr)(32.5−Tsk)+5(32.5−Tsk)
These equations are shown in the EnergyPlus source code as:
EvapHeatLoss = SkinWetSweat*EvapHeatLossMax+(1.  SkinWetSweat)*EvapHeatLossDiff
SkinWetSweat = EvapHeatLossDrySweat/EvapHeatLossMax
EvapHeatLossDiff = 0.408*(SkinVapPress  VapPress)
EvapHeatLossMax = 2.2*Hc*(SkinVapPress  VapPress)*CloPermeatEff
Here, control signals, based on setpoint temperatures in the skin and core, are introduced into the system of equations, and these equations are integrated numerically for small time increments or small increments in core and skin temperature. The control signals modulate the thermoregulatory mechanism and regulate the peripheral blood flow, the sweat rate, and the increase of metabolic heat by active muscle shivering. The development of the controlling functions of skin conductance (KS), sweat rate (Esw), and shivering (Mshiv) is based on their correlation with the deviations in skin and core temperatures from their setpoints.
ϕ=0.5+0.5exp[−5.6(w−0.4)]
Esw=ϕ[260(Tcr−36.9)+26(Tsk−33.8)]exp[(Tsk−33.8)/8.5]1.0+0.05(33.37−Tsk)2.4
wrsw−o=0.2+0.4{1.0−exp[−0.6(H/58.2−1.0)]}
where ϕ=1.0 for w≤0.4 and when w>0.4:
ϕ=0.5+0.5exp[−5.6(w−0.4)]
Within EnergyPlus, these variables are calculated using the following source code equations:
SkinCndctDilation = 42.45*CoreSignalWarmMax+8.15*CoreSignalSkinSens**0.8*SkinSignalWarmMax
SkinCndctConstriction = 1.0 + 0.4*SkinSignalColdMax
ThermCndct = 5.3+(6.75+SkinCndctDilation)/SkinCndctConstriction
WeighFac = 260.+70.*AcclPattern where the pattern of acclimation is assumed to be zero in the equation for Esw above
SweatCtrlFac = 1. + 0.05*SkinSignalSweatColdMax**2.4
EvapHeatLossDrySweat = ((WeighFac*CoreSignalSweatMax+0.1*WeighFac*SkinSignalSweatMax)
*EXP(SkinSignalSweatMax/8.5))/SweatCtrlFac
SweatSuppFac = 1. for w≤0.4
SweatSuppFac = 0.5 + 0.5*EXP(5.6*SkinWetSignal) for w>0.4
ShivResponse = 20.*CoreSignalShivMax*SkinSignalShivMax + 5.*SkinSignalShivMax
In the KSU model, two new parameters are introduced and used in correlating thermal sensations with their associated physiological responses. Instead of correlating warm thermal sensations with skin wettedness, it is correlated with a wettedness factor defined by
εwsw=(wrsw−wrsw−o)/(1.0−wrsw−o)
where:
εvc=(KSo−KS)/(KSo−KS(−4))
wrsw−o=0.2+0.4{1.0−exp[−0.6(H/58.2−1.0)]}
The wettedness factor is calculated using the following EnergyPlus equations:
SkinWetFac = (SkinWetSweat  SkinWetNeut)/(1.  SkinWetNeut)
SkinWetSweat = DrySweatRate/EvapHeatLossMax
SkinWetNeut = 0.02 + 0.4*(1.EXP(0.6*(IntHeatProdMetMax  1.)))
In addition, instead of correlating cold thermal sensation with the skin temperature in the KSU model, it is correlated with a factor identified as vasoconstriction factor defined by:
εvc=(KSo−KS)/(KSo−KS(−4))
or in the source code:
VasoconstrictFac=(ThermCndctNeutThermCndct)/(ThermCndctNeutThermCndctMin)
Thus, for a cold environment, TSV is a function of a vasoconstriction factor (εvc):
TSV=−1.46×εvc+3.75×εvc2−6.17×εvc3
and for the warm environments, TSV is defined as:
TSV=[5.0−6.56(RH−0.5)]×εwsw
In EnergyPlus, these two equations for TSV in cold and warm environments, respectively, show up as:
TSV = 1.46153*VasoconstrictFac + 3.74721*VasoconstrictFac**2
 6.168856*VasoconstrictFac**3
TSV = (5.  6.56*(RelHum  0.50))*SkinWetFac
The KSU model’s TSV was developed from experimental conditions in all temperature ranges and for clo levels between 0.05 clo and 0.7 clo and for activities levels ranging from 1 to 6 mets (Berglund 1978).
Adaptive Comfort Model
Adaptive comfort model, intended for use in naturally ventilated buildings, determines the acceptability of indoor conditions given the monthly mean outdoor air temperature and the indoor operative temperature. This is used as an index for occupant adaptation to outdoor conditions, and determines the acceptability of indoor conditions. The model also accounts for people’s clothing adaptation in naturally conditioned spaces by relating the acceptable range of indoor temperatures to the outdoor climate, so it is not necessary to estimate the clothing values for the space. No humidity or airspeed limits are required when this option is used. This section summarizes the adaptive comfort models based on the ASHRAE Standard 552010 and CEN 15251. Details are available in the two standards.
Adaptive Comfort Model Based on ASHRAE Standard 552010
In ASHRAE Standard 55, the monthly mean outdoor air temperature, used in the adaptive comfort model, is defined as the simple running average of the previous thirty daily average outdoor air temperatures.
The model defines two comfort regions: 80% Acceptability, and 90% Acceptability. If the monthly mean outdoor air temperature is not within the specified domain of 10.0 to 33.5C, the model is not applicable.
The central line of the model (shown in red), or comfort temperature, is defined as
Tot=0.31⋅To+17.8
Where
Tot – operative temperature (°C), calculated as the average of the indoor air drybulb temperature and the mean radiant temperature of zone inside surfaces
To – monthly mean outdoor air drybulb temperature (°C).
If the .stat file is provided for the simulation, To is drawn directly from the daily average temperatures in the .stat file, which provides a value for each month. If no .stat file is provided, the monthly mean outdoor temperature is a simple running average of the previous thirty daily average temperatures, calculated directly from the weather file (.epw):
α=2930
To=(1−α)⋅Tod−1+α⋅Tod−1
Tr=(Tr−avg+Tsurf)/2
Tod−i is defined as the daily average temperature of the ith previous day.
Note that the weather file must be a standard .epw containing a full year of data.
The comfort regions for 80% and 90% acceptability are symmetric about the central line.
90% Acceptability Limits: Tot = 0.31* To + 17.8± 2.5
80% Acceptability Limits: Tot = 0.31* To + 17.8± 3.5
If, using either method, To is less than 10°(C) or greater than 33.5°(C), the model is not applicable.
For a detailed description of this model, please see ASHRAE Standard 552010, Thermal Environmental Conditions for Human Occupancy.
Adaptive Comfort Model Based on European Standard EN152512007[LINK]
The EN152512007 is similar to ASHRAE 552010, but with slightly different curves of the indoor operative temperature and acceptability limits (Fig. 2). The model, intended for use in naturally ventilated buildings, determines the acceptability of indoor conditions given the 7day weighted mean outdoor air temperature and the indoor operative temperature. The 7day weighted mean outdoor air temperature (Trm) is defined as the weighted running average of the previous 7 daily average outdoor air temperatures.
This weighted running average is calculated from a full annual weather file that must be specified for the simulation. This is used as an index for occupant adaptation to outdoor conditions, and determines the acceptability of indoor conditions. The model also accounts for people’s clothing adaptation in naturally conditioned spaces by relating the acceptable range of indoor temperatures to the outdoor climate, so it is not necessary to estimate the clothing values for the space. No humidity or airspeed limits are required when this option is used. The model defines three comfort regions: Category I (90%) Acceptability, Category II (80%) Acceptability, and Category III (65%) Acceptability. If Trm is not within the specified domain of 10.0 to 30.0C, the model is not applicable.
Central line (shown as redFigure): Tot = 0.33*To + 18.8
Category I, 90% Acceptability Limits: Tot = 0.33*To + 18.8± 2.0
Category II, 80% Acceptability Limits: Tot = 0.33*To + 18.8± 3.0
Category III, 65% Acceptability Limits: Tot = 0.33*To + 18.8± 4.0
For 10°(C) < Trm < 15°(C), the comfort temperature of the lower boundaries of the comfort regions is Tcomf = 23.75°(C). That is, the lower boundaries are constant according to the same ranges above:
Category I, 90% Acceptability Limits: Tot = 23.75 2.0
Category II, 80% Acceptability Limits: Tot = 23.75  3.0
Category III, 65% Acceptability Limits: Tot = 23.75 4.0
Where
Tot – operative temperature (°C) , calculated as the average of the indoor air drybulb temperature and the mean radiant temperature of zone inside surfaces
To – mean outdoor air drybulb temperature (°C), calculated as the weighted mean of the previous 7day daily mean outdoor air drybulb temperature (Tod):
To = (1  α)*{Tod−1 + α *Tod−2 + α2 *Tod−3 + α3 *Tod−4 + α4 *Tod−5 + α5 *Tod−6 + α6 *Tod−7}
To = (1  α)*Tod−1 + α *To−1
α = 0.8
Dynamic Clothing Model[LINK]
In most building energy simulations, thermal comfort condition is calculated based on the assumption that the clothing insulation is equal to a constant value of 0.5 clo during the cooling season and 1.0 clo during heating season. Usually those two values are used and the change from 0.5 to 1 or viceversa is made suddenly from one day to another. In addition, there is no standardized guideline on how to set clothing insulation schedules in the international standards. This simplified assumption may lead to systems that are incorrectly sized and operated and to the incorrect assessment of comfort conditions. In reality, occupants frequently adjust their clothing depending on the thermal conditions around them, as opposed to the assumption of constant clothing values. Therefore, the clothing insulation variation should be captured during the building simulation to realistically model HVAC systems. In order to overcome the limitations of the constant clothing insulation assumption, three new predictive clothing insulation models were developed by Schiavon and Lee (2012) based on 6,333 selected observations taken from ASHRAE RP884 and RP921 databases. The first and third models vary the clothing insulation as a function of outdoor air temperature measured at 6 o’clock and the second model takes into account both 6 o’clock outdoor air temperature and indoor operative temperature when adjusting the clothing insulation. The dynamic clothing models should be implemented in dynamic building energy simulation.
The model proposed to ASHRAE 55 is described below.
For ta(out,6) < 5°C, Icl=1 .00
For 5°C ≤ ta(out,6) < 5°C, Icl=0.818−0.0364∗ta(out,6)
For 5°C ≤ ta(out,6) < 26°C Icl=10(−0.1635−0.0066∗ta(out,6))
or ta(out,6) ≥ 26°C Icl=0.46
Where, Icl is the clothing insulation value, ta(out,6) is the outdoor air temperature measured at 6 o’clock in the morning. The following figure illustrates the proposed clothing insulation model.
In the figure above, clothing insulation schedule for a fixed model (blue) typically used in energy simulation software and for the clothing model based on outdoor air temperature measured at 6 o’clock. Climate data for Chicago O’hare International Airport has been used.
The following figure illustrates the new clothing insulation model.
The dynamic predictive clothing insulation model is implemented into EnergyPlus for realistic energy simulation. Addenda A to ASHRAE 55 with the clothing model has been approved by the ASHRAE committee and the chance to be included in ASHRAE 552013 is high.
Schiavon S, Lee KH. 2013. Dynamic predictive clothing insulation models based on outdoor air and indoor operative temperatures. Building and Environment. Volume 59, 250260. http://dx.doi.org/10.1016/j.buildenv.2012.08.024 (link to the journal) http://escholarship.org/uc/item/3338m9qf (link to the freely available preprint version)
Lee KH, Schiavon S. 2013. Influence of three dynamic predictive clothing insulation models on building energy use, HVAC sizing and thermal comfort. Submitted to Journal. http://escholarship.org/uc/item/3sx6n876 (link to the freely available preprint version)
Mean Radiant Temperature Calculation[LINK]
There are three options to calculate mean radiant temperature in the thermal comfort models. One is the zone averaged MRT, another is the surface weighted MRT, and the other is angle factor MRT. The zone averaged MRT is calculated on the assumption that a person is in the center of a space, whereas the surface weighted MRT is calculated in consideration of the surface that a person is closest to, and the angle factor MRT is calculated based on angle factors between a person and the different surfaces in a space. Here, the surface weighted MRT is the average temperature of the selected surface and zone averaged MRT and is intended to represent conditions in the limit as a person gets closer and closer to a particular surface. In that limit, half of the person’s radiant field will be dominated by that surface and the other half will be exposed to the rest of the zone. Note that the surface weighted MRT is only an approximation. The angle factor MRT is the mean temperature of the surrounding surface temperatures weighted according to the magnitude of the respective angle factors and allows the user to more accurately predict thermal comfort at a particular location within a space.
Nomenclature and variable list for MRT calculation
Tr

Mean radiant temperature 
°C 
 
RadTemp 
Tr−avg

Zone averaged radiant temperature 
°C 
 
ZoneRadTemp 
Tsurf

Surface temperature 
°C 
 
SurfaceTemp 
Fsurf

Angle factor between person and surface 
 
01 
AngleFactor 
Description of the model and algorithm[LINK]
The zone averaged MRT is calculated without weighting any surface temperature of the space.
Tr=Tsurf−1Fsurf−1+Tsurf−2Fsurf−2+.........+Tsurf−nFsurf−n
RadTemp = MRT(ZoneNum)
The surface weighted MRT is the average temperature of the zone averaged MRT and the temperature of the surface that a person is closest to.
Tr=(Tr−avg+Tsurf)/2
ZoneRadTemp = MRT(ZoneNum)
SurfaceTemp = GetSurfaceTemp(People(PeopleNum)%SurfacePtr)
RadTemp = (ZoneRadTemp + SurfaceTemp)/2.0
The angle factor MRT is the mean value of surrounding surface temperatures weighted by the size of the respective angle factors between a person and each surface.
Tr=Tsurf−1Fsurf−1+Tsurf−2Fsurf−2+.........+Tsurf−nFsurf−n
SurfTempAngleFacSummed = SurfTempAngleFacSummed &
 SurfaceTemp * AngleFactorList(AngleFacNum)%AngleFactor(SurfNum)
RadTemp = SurfTempAngleFacSummed
ASHRAE. 1984. “High Intensity Infrared Radiant Heating”, 1984 Handbook of Systems and Equipment, American Society of Heating, Refrigerating and Air Conditioning Engineers, Atlanta, GA, Chapter 18.
ASHRAE. 1985. “Physiological Principles for Comfort and Health,” 1985 Handbook of Fundamentals, American Society of Heating, Refrigerating and Air Conditioning Engineers, Atlanta, GA, Chapter 8.
ASHRAE. 1993. “Physiological Principles and Thermal Comfort”, 1993 ASHRAE Handbook of Fundamentals, American Society of Heating, Refrigerating and Air Conditioning Engineers, Atlanta, GA, Chapter 8.
ASHRAE. 2010. “Standard 552010 – Thermal Environmental Conditions for Human Occupancy (ANSI approved)”, American Society of Heating, Refrigerating and Air Conditioning Engineers, Atlanta, GA.
Azer, N.Z., Hsu, S. 1977. “The prediction of Thermal Sensation from Simple model of Human Physiological Regulatory Response”, ASHRAE Trans., Vol.83, Pt 1.
Berglund, Larry. 1978. “Mathematical Models for Predicting the Thermal Comfort Response of Building Occupants”, ASHRAE Trans., Vol.84.
Doherty, T.J., Arens, E. 1988. “Evaluation of the Physiological Bases of Thermal Comfort Models”, ASHRAE Trans., Vol.94, Pt 1.
Du Bois, D. and E.F. 1916. “A Formula to Estimate Approximate Surface Area, if Height and Weight are Known”, Archives of internal Medicine, Vol.17.
CEN. 2007. “Standard EN15251 Indoor environmental input parameters for design and assessment of energy performance of buildings addressing indoor air quality, thermal environment, lighting and acoustics”. Bruxelles: European committee for Standardisation.
Fanger, P.O. 1970. Thermal ComfortAnalysis and Applications in Environmental Engineering, Danish Technical Press, Copenhagen.
Fanger, P.O. 1986. “Radiation and Discomfort”, ASHRAE Journal. February 1986.
Fanger P.O. 1967. “Calculation of Thermal Comfort: Introduction of a Basic Comfort Equation”, ASHRE Trans., Vol.73, Pt 2.
Fountain, Marc.E., Huizenga, Charlie. 1997 “A Thermal Sensation Prediction Tool for Use by the Profession”, ASHRAE Trans., Vol.103, Pt 2.
Gagge, A.P., Stolwijk, J. A. J., Nishi, Y. 1970. “An Effective Temperature Scale Based on a Simple Model of Human Physiological Regulatory Response”, ASHRAE Trans., Vol.70, Pt 1.
Gagge, A.P., Fobelets, A.P., Berglund, L. G. 1986. “A Standard Predictive Index of Human Response to the Thermal Environment”, ASHRAE Trans., Vol.92, Pt 2.
Hsu, S. 1977. “A Thermoregulatory Model for Heat Acclimation and Some of its Application”, Ph. D. Dissertation, Kansas State University.
IntHout, D. 1990. “Thermal Comfort Calculation / A Computer Model”, ASHRAE Trans., Vol.96, Pt 1.
ISO. 1983. “Determination of the PMV and PPD Indices and Specification of the Conditions for Thermal Comfort”, DIS 7730, Moderate Thermal Environment, 1983.
Occupant Thermal Comfort[LINK]
The integration of a sophisticated building thermal analysis tool with thermal comfort models allows one to perform an energy analysis on a zone and simultaneously determine if the environmental control strategy will be sufficient for the occupants to be thermally comfortable. This section is intended to provide background on thermal comfort, present an overview of state of the art thermal comfort models, and present the mathematical models that have been incorporated into EnergyPlus.
Thermal comfort modeling in EnergyPlus is controlled primarily by the People input object. This includes input for selecting the type of thermal comfort model that is desired by the user as well as parameters that serve as inputs to the thermal comfort models. These parameters include the activity level, the work efficiency, the air velocity and the clothing insulation level for people within the space. All four of these parameters can be scheduled. More information on the People input object can be found in the EnergyPlus Input/Output Reference. More information on how each of these parameters is used and the specific modeling equations for the thermal comfort models can be found below.
Background on Thermal Comfort Models[LINK]
Throughout the last few decades, researchers have been exploring the thermal, physiological and psychological response of people in their environment in order to develop mathematical models to predict these responses. Researchers have empirically determined building occupants’ thermal responses to the combined thermal effect of the personal, environmental and physiological variables that influence the condition of thermal comfort.
There are two personal variables that influence the condition of thermal comfort: the thermal resistance of the clothing (Icl), and the metabolic rate (H/ADu).
The thermal resistance of the clothing (Icl) is measured in units of “clo.” The 1985 ASHRAE Handbook of Fundamentals (ASHRAE 1985) suggests multiplying the summation of the individual clothing items clo value by a factor of 0.82 to determine the overall clo value for clothing ensembles.
The metabolic rate (H/ADu) is a measure of the internal heat production rate of an occupant (H) in W/hr per unit of “Dubois” body surface area (ADu) in units of m2. The DuBois body surface area is given by :
ADu=0.202(weight)0.425(height)0.725
Using this equation, an area of 1.8 m2 represents the surface area of an average person of weight 70 kg. and height 1.73 m (Fanger 1967). In EnergyPlus, this average person area of 1.8 m2 is used for the body surface area all thermal comfort models. The metabolic rate is measured in mets, where 1 met = 58.2 W/m2.
The environmental variables that influence the conditions of thermal comfort include:
Air Temperature (Ta),
Mean Radiant Temperature (Tr),
Relative air velocity (v),
Water vapor pressure in ambient air (Pa)
The Air Temperature (Ta), a direct environmental index, is the drybulb temperature of the environment. It is calculated using either the mean air temperature for the zone in question or the air temperature calculated by some of the more detailed zone air models available in EnergyPlus. The Mean Radiant Temperature (Tr) is a rationally derived environmental index defined as the uniform blackbody temperature that would result in the same radiant energy exchange as in the actual environment. This value can be calculated in several ways as outlined in the Mean Radiant Temperature Calculation section later in this document. The Relative air velocity (v) a direct environmental index is a measure of the air motion obtainable via a hot wire or vane anemometers. The Water vapor pressure in ambient air (Pa) is a direct environmental index.
The physiological variables that influence the conditions of thermal comfort include:
Skin Temperature (Tsk),
Core or Internal Temperature (Tcr),
Sweat Rate,
Skin Wettedness (w),
Thermal Conductance (K) between the core and skin.
where the Skin Temperature (Tsk), the Core Temperature (Tcr) and the Sweat Rate are physiological indices. The Skin Wettedness (w) is a rationally derived physiological index defined as the ratio of the actual sweating rate to the maximum rate of sweating that would occur if the skin were completely wet.
One more consideration is important in dealing with thermal comfort  the effect of asymmetrical heating or cooling. This could occur when there is a draft or when there is a radiant flux incident on a person (which is what is of primary interest to us here). Fanger (1967) noted that the human regulatory system is quite tolerant of asymmetrical radiant flux. A reasonable upper limit on the difference in mean radiant temperature (Tr) from one direction to the opposing direction is 15°C. (ASHRAE 1984). This limit is lower if there is a high air velocity in the zone.
Mathematical Models for Predicting Thermal Comfort[LINK]
Many researchers have been exploring ways to predict the thermal sensation of people in their environment based on the personal, environmental and physiological variables that influence thermal comfort. From the research done, some mathematical models that simulate occupants’ thermal response to their environment have been developed. Most thermal comfort prediction models use a seven or nine point thermal sensation scale, as in the following tables.
The most notable models have been developed by P.O. Fanger (the Fanger Comfort Model), the J. B. Pierce Foundation (the Pierce TwoNode Model), and researchers at Kansas State University (the KSU TwoNode Model). Berglund (1978) presents a detailed description of the theory behind these three models.
Note for all Thermal Comfort reporting: Though the published values for thermal comfort “vote” have a discrete scale (e.g. –3 to +3 or –4 to +4), the calculations in EnergyPlus are carried out on a continuous scale and, thus, reporting may be “off the scale” with specific conditions encountered in the space. This is not necessarily an error in EnergyPlus – rather a different approach that does not take the “limits” of the discrete scale values into account.
The main similarity of the three models is that all three apply an energy balance to a person and use the energy exchange mechanisms along with experimentally derived physiological parameters to predict the thermal sensation and the physiological response of a person due to their environment. The models differ somewhat in the physiological models that represent the human passive system (heat transfer through and from the body) and the human control system (the neural control of shivering, sweating and skin blood flow). The models also differ in the criteria used to predict thermal sensation. However, all three models use information from the People statement and the thermal comfort model is selected via the People statement in a user’s input file. Scheduled parameters such as the activity level, work efficiency, air velocity, and clothing insulation level all have a direct bearing on the thermal comfort models. For more information on the input of these parameters, see the People statement in the EnergyPlus Input/Output Reference. For more information on how each individual thermal comfort model uses these parameters, please consult the next several sections.
Fanger Comfort Model[LINK]
Fanger’s Comfort model was the first human comfort model that was developed. It was first published in 1967 (Fanger 1967), then updated in 1970 (Fanger 1970), and helped set the stage for the Pierce and KSU models. The mathematical model developed by P.O. Fanger is probably the most well known of the three models and is the easiest to use because it has been put in both chart and graph form.
Fanger Model Nomenclature List[LINK]
Description of the model and algorithm[LINK]
Fanger developed the model based on the research he performed at Kansas State University and the Technical University of Denmark. Fanger used the sevenpoint form of a thermal sensation scale along with numerous experiments involving human subjects in various environments. He related the subjects’ response to various environmental variables which influence the condition of thermal comfort. Fanger’s model is based upon an energy analysis that takes into account all the modes of energy loss (L) from the body, including: the convection and radiant heat loss from the outer surface of the clothing, the heat loss by water vapor diffusion through the skin, the heat loss by evaporation of sweat from the skin surface, the latent and dry respiration heat loss and the heat transfer from the skin to the outer surface of the clothing. The model assumes that the person is thermally at steady state with the interior environment.
M=L
L=Qres+Qdry+Esk
Qres=Eres+Cres=0.000017M(5867−Pa)+0.0014M(34−Ta)
In the EnergyPlus source code, Qres, Eres and Cres are calculated using the following equations, respectively:
RespHeatLoss = LatRespHeatLoss + DryRespHeatLoss
LatRespHeatLoss = 0.000017*ActLevel*(5867.  VapPress)
DryRespHeatLoss = 0.0014*ActLevel*(34. AirTemp)
Qr=hr×fcl(Tcl−Tr)
hc=8.6×v0.53
Ersw=0.42(H−58.2)for(H>58.2)
In the EnergyPlus source code, Qdry, Qc and Qr are calculated using the following equations, respectively:
DryHeatLoss = ConvHeatLoss + RadHeatLoss
ConvHeatLos = CloBodyRat*Hc*(CloSurfTemp  AirTemp)
RadHeatLoss = 3.96*CloBodyRat*[(0.01*AbsCloSurfTemp)**4  (0.01*AbsRadTemp)**4]
Ersw=0for(H≤58.2)
Esk=Ersw+Ediff
PPD=100.0−95.0∗EXP(−0.03353PMV4−0.2179PMV2)
Ediff=wdiff×Emax
In the EnergyPlus source code, Ersw, Ef and Esk are calculated using the following equations, respectively:
EvapHeatLossRegComf = 0.42*(IntHeatProd  ActLevelConv) or 0.0 depending on H value
EvapHeatLossDiff = 0.00305*(5733.0  6.99*IntHeatProd  VapPress)
EvapHeatLoss = EvapHeatLossRegComf + EvapHeatLossDiff
The Fanger model was developed using the thermal sensation votes from subjects at KSU and in Denmark, and the Predicted Mean Vote (PMV) thermal sensation scale is based on how the energy loss (L) deviates from the metabolic rate (M) in the following form:
PMV=(0.303e−0.036M+0.028)(H−L)
In the EnergyPlus source code, PMV is calculated using the following equations in conjunction with the equations shown above:
PMV = ThermSensTransCoef*(IntHeatProd  EvapHeatLoss  RespHeatLoss  DryHeatLoss)
ThermSensTransCoef = 0.303*EXP(0.036*ActLevel) + 0.028
IntHeatProd = ActLevel  WorkEff
Predicted Percent of Dissatisfied (PPD) people at each PMV is calculated as follows:
SIGsk=Tsk−33.7
Pierce TwoNode Model[LINK]
The Pierce TwoNode model was developed at the John B. Pierce Foundation at Yale University. The model has been continually expanding since its first publication in 1970 (Gagge et.al. 1970). The most recent version on the model appears in the 1986 ASHRAE Transactions (Gagge et.al. 1986).
Pierce TwoNode Model Nomenclature List[LINK]
Description of the model and algorithm[LINK]
The Pierce model thermally lumps the human body as two isothermal, concentric compartments, one representing the internal section or core (where all the metabolic heat is assumed to be generated) and the other representing the skin. This allows the passive heat conduction from the core compartment to the skin to be taken into account. The boundary line between the two compartments changes with the skin blood flow rate per unit skin surface area (SKBF in L/h•m2) and is described by alpha – the fraction of total body mass attributed to the skin compartment (Doherty and Arens 1988).
α=0.0417737+0.7451832/(SKBF+0.585417)
In the EnergyPlus code, α or SkinMassRat is calculated using the equivalent equation:
SkinMassRat = 0.0417737 + 0.7451832/(SkinBloodFlow + 0.585417)
Furthermore, the model takes into account the deviations of the core, skin, and mean body temperature weighted by alpha from their respective setpoints. Thermoregulatory effector mechanisms (regulatory sweating, skin blood flow, and shivering) are defined in terms of thermal signals from the core, skin and body in ∘C (Doherty and Arens 1988) by the following equations:
SIGb=Tb−36.49
SKBF=(6.3+Cdil×SIGcr)/(1+Str×(SIGsk))
SWreg=CswSIGbe(SIGsk/10.7)
In the EnergyPlus code, when SkinTemp is greater than SkinTempSet (33.7∘C), then SkinThermSigWarm = SkinTemp  SkinTempSet and SkinThermSigCold = 0.0. Otherwise, when SkinTemp is less than or equal to SkinTempSet, then SkinThermSigWarm = 0.0 and SkinThermSigCold = SkinTempSet  SkinTemp.
In the EnergyPlus code, when CoreTemp is greater than CoreTempSet (36.8∘C), then CoreThermSigWarm = CoreTemp  CoreTempSet and CoreThermSigCold = 0.0. Otherwise, when CoreTemp is less than or equal to CoreTempSet, then CoreThermSigWarm = 0.0 and CoreThermSigCold = CoreTempSet  CoreTemp.
In the EnergyPlus code, when AvgBodyTemp is greater than AvgBodyTempSet(36.49∘C), then BodyThermSigWarm = AvgBodyTemp  AvgBodyTempSet and BodyThermSigCold = 0.0. Otherwise, when AvgBodyTemp is less than or equal to AvgBodyTempSet, then BodyThermSigWarm = 0.0 and BodyThermSigCold = AvgBodyTempSetAvgBodyTemp.
The skin blood flow rate in L/hr⋅m2 is calculated using:
Mshiv=19.4(SIGcr)(SIGsk)
In the EnergyPlus code, this equation is calculated using the following equations:
SkinBloodFlow = (6.3 + VasodilationFac)/(1. + VasoconstrictFac)
VasodilationFac = SkinBloodFlowConst*CoreWarmDelTemp
VasoconstrictFac = Str*SkinColdDelTemp
The rate of regulatory sweating in g/hr•m2 is calculated using:
M=Mact+Mshiv
and the shivering response in W/m2 is calculated using:
L=Qres+Qdry+Esk+W
In the EnergyPlus code, the sweat and shivering response determined by the above two equations are calculated using the follow expressions:
RegSweat = SweatContConst*BodyThermSigWarm*EXP(SkinThermSigWarm/10.7)
ShivResponse = 19.4*SkinThermSigCold*CoreThermSigCold
The latest version of the Pierce model (Fountain and Huizenga 1997) discusses the concepts of SET* and ET*. The Pierce model converts the actual environment into a “standard environment” at a Standard Effective Temperature, SET*. SET* is the drybulb temperature of a hypothetical environment at 50% relative humidity for subjects wearing clothing that would be standard for the given activity in the real environment. Furthermore, in this standard environment, the same physiological strain, i.e. the same skin temperature, skin wettedness and heat loss to the environment, would exist as in the real environment. The Pierce model also converts the actual environment into a environment at an Effective Temperature, ET*, that is the drybulb temperature of a hypothetical environment at 50% relative humidity and uniform temperature (Ta = MRT) where the subjects would experience the same physiological strain as in the real environment.
In the latest version of the model, it is suggested that the classical Fanged PMV be modified by using ET* or SET* instead of the operative temperature. This gives a new index PMV* which is proposed for dry or humid environments. It is also suggested that PMV* is very responsive to the changes in vapor permeation efficiency of the occupants clothing.
Qres=Eres+Cres=0.0023M(44−Pa(torr))+0.0014M(34−Ta)Qres=0.017251M(5.8662−Pa(kPa))+0.0014M(34−Ta)
which in the EnergyPlus code takes on the form ActLevel = ActLevel + ActShiv.
Qdry=Qc+Qr
Qc=hc×fcl(Tcl−Ta)
In the EnergyPlus code, the calculation of Qres or RespHeatLoss is accomplished using the following equations:
RespHeatLoss = LatRespHeatLoss + DryRespHeatLoss
LatRespHeatLoss = 0.017251*ActLevel*(5.8662  VapPress)
DryRespHeatLoss = 0.0014*ActLevel*(34. AirTemp)
Qr=hr×fcl(Tcl−Tr)
hc=8.6×v0.53
hc=5.66(M/58.2−0.85)0.39
Qdry or DryHeatLoss is calculated in the EnergyPlus code using the following equation:
DryHeatLoss = CloBodyRat*(Hc*(CloSurfTemp  AirTemp) + Hr*(CloSurfTemp  RadTemp))
In Pierce model, the convective heat transfer coefficient, hc, varies with the air velocity around body and metabolic rate. The model uses the maximum value of following equations:
hr=4.×feffεσ((Tcl+Tr)/2+273.15)3
Tcl=(Tsk/Icl+fcl(hcTa+hrTr))/(1/Icl+fcl(hc+hr))
In the EnergyPlus code, the following equations are used to calculate hc or Hc:
Hc = 8.6*AirVel**0.53 where AirVel is reset to 0.137 if it drops below this value
HcAct = 5.66*(ActMet  0.85)**0.39 when ActMet is greater than 0.9
When ActMet is greater than 0.9, EnergyPlus uses the larger value of Hc and HcAct for hc.
Also, in the model, the radiant heat transfer coefficient, hr, is defined by following equation (Doherty and Arens 1988):
Esk=Ersw+Ediff
In the EnergyPlus code, hr or Hr is calculated using the following equation:
Hr = 4.*RadSurfEff*StefanBoltz*((CloSurfTemp + RadTemp)/2. + TAbsConv)**3
In the Pierce model, Tcl in ∘C is estimated by each iteration using following equation:
Ersw=0.68×SWreg
Or as shown in the EnergyPlus code:
CloSurfTemp = (CloCond*SkinTemp + CloBodyRat*(Hc*AirTemp+Hr*RadTemp))
/(CloCond + CloBodyRat*(Hc+Hr))
The total evaporative heat loss from the skin, Esk, includes evaporation of water produced by regulatory sweating, Ersw, and evaporation of water vapor that diffuses through the skin surface, Ediff.
Ediff=wdiff×Emax
Emax=he′(Psk−Pa)
Qcrsk=(5.28+1.163SKBF)(Tcr−Tsk)
Within the EnergyPlus source code, these equations are shown as:
EvapHeatLoss = EvapHeatLossRegSweat + EvapHeatLossRegDiff
RegHeatLoss = 0.68*RegSweat
DiffHeatLoss = SkinWetDiff*MaxEvapHeatLoss
where 0.68 is the passive water vapor diffusion rate in g/h·m2·Torr and
wdiff=0.06(1−wrsw)
Ssk=Qc−s−Qc−Qr−Esk
wrsw=Ersw/Emax
These equations within the EnergyPlus source code appear as:
SkinWetDiff = .06*(1.SkinWetSweat)
MaxEvapHeatLoss = (1./TotEvapHeatResist)*(SatSkinVapPress  VapPress)
SkinWetSweat = EvapHeatLossRegSweat/MaxEvapHeatLoss
The Pierce model has one additional heat flow term describing the heat transfer between the internal core compartment and the outer skin shell (Doherty and Arens 1988).
Scr=M−W−Qres−Qc−s
or as shown in the EnergyPlus source code:
HeatFlow = (CoreTempSkinTemp)*(5.28 + 1.163*SkinBloodFlow)
where 5.28 is the average body tissue conductance in W/m$^2 $ K and 1.163 is the thermal capacity of blood in W ⋅ h/L ⋅ K.
Thus, individual heat balance equations for core and skin compartments are expressed using this term, Qc−s. New temperatures of core, skin and body are calculated for each iteration from rates of heat storage in the core and skin.
Tb−c=(0.185/58.2)(M−W)+36.313
Tb−h=(0.359/58.2)(M−W)+36.664
These equations are represented by the following source code lines in EnergyPlus:
SkinHeatStorage = HeatFlow  DryHeatLoss  EvapHeatLoss
CoreHeatStorage = IntHeatProd  RespHeatLoss  HeatFlow
As a result, PMV evaluated for ET* and SET* can be calculated using:
PMVET=(0.303e−0.036M+0.028)(H−LET∗)
PMVSET=(0.303e−0.036M+0.028)(H−LSET∗)
These equations are represented by the following lines of source code in EnergyPlus:
ThermSensTransCoef = 0.303*EXP(0.036*ActLevel) + 0.028
PMVET = ThermSensTransCoef*(IntHeatProd  EvapHeatLossDiff  EvapHeatLossRegComf
 RespHeatLoss  DryHeatLossET)
PMVSET = ThermSensTransCoef*(IntHeatProd  EvapHeatLossDiff
 EvapRegHeatLossRegComf  RespHeatLoss  DryHeatLossSET)
Besides PMV*, the Pierce Two Node Model uses the indices TSENS and DISC as predictors of thermal comfort. Where TSENS is the classical index used by the Pierce foundation and is a function of the mean body temperature. DISC is defined as the relative thermoregulatory strain that is needed to bring about a state of comfort and thermal equilibrium. DISC is a function of the heat stress and heat strain in hot environments and equal to TSENS in cold environments. In summary, the Pierce Model, for our purposes, uses four thermal comfort indices; PMVETa function of ET*, PMVSET a function of SET*, TSENS and DISC.
TSENSc=0.68175(Tb−Tb−c)
TSENSh=4.7(Tb−Tb−c)/(Tb−h−Tb−c)
When Tb≤Tb−c, TSENS is calculated using:
WcrCcrdTcr/dt=M−W−Qres−KS(Tcr−Tsk)
Otherwise, when Tb>Tb−c, TSENS is calculated using:
WskCskdTsk/dt=KS(Tcr−Tsk)−Qdry−Esk
Finally, DISC is calculated using the following equation:
DISC=5.(Ersw−Ersw−comf)/(Emax−Ersw−comf−Ediff)
These final equations for the Pierce Two Node Model are calculated using the following lines of source code:
AvgBodyTempLow = (0.185/ActLevelConv)*IntHeatProd + 36.313
AvgBodyTempHigh = (0.359/ActLevelConv)*IntHeatProd + 36.664
TSENS = .68175*(AvgBodyTempAvgBodyTempLow)
TSENS = 4.7*(AvgBodyTempAvgBodyTempLow)/(AvgBodyTempHighAvgBodyTempLow)
DISC = 5.*(EvapHeatLossRegSweatEvapHeatLossRegComf)
/(MaxEvapHeatLossEvapHeatLossRegComfDiffHeatLoss)
KSU TwoNode Model[LINK]
The KSU TwoNode Model, developed at Kansas State University, was published in 1977 (Azer and Hsu 1977). The KSU model is quite similar to that of the Pierce Foundation. The main difference between the two models is that the KSU model predicts thermal sensation (TSV) differently for warm and cold environment.
KSU Two Node Model Nomenclature List[LINK]
Description of the model and algorithm[LINK]
The KSU TwoNode Model is based on the changes that occur in the thermal conductance between the core and the skin temperature in cold environments. In warm environments, the thermal conductance of the model is based on changes in the skin wettedness.
In this model, metabolic heat production is generated in the core which exchanges energy with the environment by respiration, and the skin exchanges energy with the interior environment by convection and radiation. In addition, body heat is dissipated through the evaporation of sweat and/or water vapor diffusion through the skin. These principles are outlined in the following system of equations.
Qres=Eres+Cres=0.0023M(44−Pa(torr))+0.0014M(34−Ta)
Qdry=Qc+Qr=hfclFcl(Tsk−To)
where Qres and Qdry are calculated using:
h=hc+hr
hc=8.3√v
The calculation of Qres (RespHeatLoss) and Qdry (DryHeatLoss) in the source code are accomplished through the following equations:
RespHeatLoss = LatRespHeatLoss + DryRespHeatLoss
LatRespHeatLoss = 0.0023*ActLevelTot*(44.  VapPress)
DryRespHeatLoss = 0.0014*ActLevelTot*(34.  AirTemp)
DryHeatLoss = H*CloBodyRat*CloThermEff*(SkinTemp  OpTemp)
In the equation for Qdry, h and To are calculated using:
hr=ϵavg(3.87+0.031Tr)
To=(hcTa+hrTr)/(hc+hr)
For(Esw≤Emax),Esk=Esw+(1−wrsw)Ediff
For(Esw>Emax),Esk=Emax
These equations are represented by the following source code equations:
H = Hc + Hr
Hc = 8.3*SQRT(AirVel) where AirVel is reset to 0.137 if it drops below this value
Hr = EmissAvg*(3.87 + 0.031*RadTemp)
OpTemp = (Hc*AirTemp + Hr*RadTemp)/H
In the KSU TwoNode Model, evaporative heat losses are defined by these equations:
Ediff=0.408(Psk−Pa)
Emax=2.2hcFpcl(Psk−Pa)
KS=5.3+[6.75+42.45(Tcr−36.98)+8.15(Tcr−35.15)0.8(Tsk−33.8)]1+0.4(32.1−Tsk)
Mshiv=20(36.9−Tcr)(32.5−Tsk)+5(32.5−Tsk)
These equations are shown in the EnergyPlus source code as:
EvapHeatLoss = SkinWetSweat*EvapHeatLossMax+(1.  SkinWetSweat)*EvapHeatLossDiff
SkinWetSweat = EvapHeatLossDrySweat/EvapHeatLossMax
EvapHeatLossDiff = 0.408*(SkinVapPress  VapPress)
EvapHeatLossMax = 2.2*Hc*(SkinVapPress  VapPress)*CloPermeatEff
Here, control signals, based on setpoint temperatures in the skin and core, are introduced into the system of equations, and these equations are integrated numerically for small time increments or small increments in core and skin temperature. The control signals modulate the thermoregulatory mechanism and regulate the peripheral blood flow, the sweat rate, and the increase of metabolic heat by active muscle shivering. The development of the controlling functions of skin conductance (KS), sweat rate (Esw), and shivering (Mshiv) is based on their correlation with the deviations in skin and core temperatures from their setpoints.
ϕ=0.5+0.5exp[−5.6(w−0.4)]
Esw=ϕ[260(Tcr−36.9)+26(Tsk−33.8)]exp[(Tsk−33.8)/8.5]1.0+0.05(33.37−Tsk)2.4
wrsw−o=0.2+0.4{1.0−exp[−0.6(H/58.2−1.0)]}
where ϕ=1.0 for w≤0.4 and when w>0.4:
ϕ=0.5+0.5exp[−5.6(w−0.4)]
Within EnergyPlus, these variables are calculated using the following source code equations:
SkinCndctDilation = 42.45*CoreSignalWarmMax+8.15*CoreSignalSkinSens**0.8*SkinSignalWarmMax
SkinCndctConstriction = 1.0 + 0.4*SkinSignalColdMax
ThermCndct = 5.3+(6.75+SkinCndctDilation)/SkinCndctConstriction
WeighFac = 260.+70.*AcclPattern where the pattern of acclimation is assumed to be zero in the equation for Esw above
SweatCtrlFac = 1. + 0.05*SkinSignalSweatColdMax**2.4
EvapHeatLossDrySweat = ((WeighFac*CoreSignalSweatMax+0.1*WeighFac*SkinSignalSweatMax)
*EXP(SkinSignalSweatMax/8.5))/SweatCtrlFac
SweatSuppFac = 1. for w≤0.4
SweatSuppFac = 0.5 + 0.5*EXP(5.6*SkinWetSignal) for w>0.4
ShivResponse = 20.*CoreSignalShivMax*SkinSignalShivMax + 5.*SkinSignalShivMax
In the KSU model, two new parameters are introduced and used in correlating thermal sensations with their associated physiological responses. Instead of correlating warm thermal sensations with skin wettedness, it is correlated with a wettedness factor defined by
εwsw=(wrsw−wrsw−o)/(1.0−wrsw−o)
where:
εvc=(KSo−KS)/(KSo−KS(−4))
wrsw−o=0.2+0.4{1.0−exp[−0.6(H/58.2−1.0)]}
The wettedness factor is calculated using the following EnergyPlus equations:
SkinWetFac = (SkinWetSweat  SkinWetNeut)/(1.  SkinWetNeut)
SkinWetSweat = DrySweatRate/EvapHeatLossMax
SkinWetNeut = 0.02 + 0.4*(1.EXP(0.6*(IntHeatProdMetMax  1.)))
In addition, instead of correlating cold thermal sensation with the skin temperature in the KSU model, it is correlated with a factor identified as vasoconstriction factor defined by:
εvc=(KSo−KS)/(KSo−KS(−4))
or in the source code:
VasoconstrictFac=(ThermCndctNeutThermCndct)/(ThermCndctNeutThermCndctMin)
Thus, for a cold environment, TSV is a function of a vasoconstriction factor (εvc):
TSV=−1.46×εvc+3.75×εvc2−6.17×εvc3
and for the warm environments, TSV is defined as:
TSV=[5.0−6.56(RH−0.5)]×εwsw
In EnergyPlus, these two equations for TSV in cold and warm environments, respectively, show up as:
TSV = 1.46153*VasoconstrictFac + 3.74721*VasoconstrictFac**2
 6.168856*VasoconstrictFac**3
TSV = (5.  6.56*(RelHum  0.50))*SkinWetFac
The KSU model’s TSV was developed from experimental conditions in all temperature ranges and for clo levels between 0.05 clo and 0.7 clo and for activities levels ranging from 1 to 6 mets (Berglund 1978).
Adaptive Comfort Model
Adaptive comfort model, intended for use in naturally ventilated buildings, determines the acceptability of indoor conditions given the monthly mean outdoor air temperature and the indoor operative temperature. This is used as an index for occupant adaptation to outdoor conditions, and determines the acceptability of indoor conditions. The model also accounts for people’s clothing adaptation in naturally conditioned spaces by relating the acceptable range of indoor temperatures to the outdoor climate, so it is not necessary to estimate the clothing values for the space. No humidity or airspeed limits are required when this option is used. This section summarizes the adaptive comfort models based on the ASHRAE Standard 552010 and CEN 15251. Details are available in the two standards.
Adaptive Comfort Model Based on ASHRAE Standard 552010
In ASHRAE Standard 55, the monthly mean outdoor air temperature, used in the adaptive comfort model, is defined as the simple running average of the previous thirty daily average outdoor air temperatures.
The model defines two comfort regions: 80% Acceptability, and 90% Acceptability. If the monthly mean outdoor air temperature is not within the specified domain of 10.0 to 33.5C, the model is not applicable.
Acceptable operative temperature ranges for naturally conditioned spaces (ASHRAE Standard 552010)
The central line of the model (shown in red), or comfort temperature, is defined as
Tot=0.31⋅To+17.8
Where
Tot – operative temperature (°C), calculated as the average of the indoor air drybulb temperature and the mean radiant temperature of zone inside surfaces
To – monthly mean outdoor air drybulb temperature (°C).
If the .stat file is provided for the simulation, To is drawn directly from the daily average temperatures in the .stat file, which provides a value for each month. If no .stat file is provided, the monthly mean outdoor temperature is a simple running average of the previous thirty daily average temperatures, calculated directly from the weather file (.epw):
α=2930
To=(1−α)⋅Tod−1+α⋅Tod−1
Tr=(Tr−avg+Tsurf)/2
Tod−i is defined as the daily average temperature of the ith previous day.
Note that the weather file must be a standard .epw containing a full year of data.
The comfort regions for 80% and 90% acceptability are symmetric about the central line.
90% Acceptability Limits: Tot = 0.31* To + 17.8± 2.5
80% Acceptability Limits: Tot = 0.31* To + 17.8± 3.5
If, using either method, To is less than 10°(C) or greater than 33.5°(C), the model is not applicable.
For a detailed description of this model, please see ASHRAE Standard 552010, Thermal Environmental Conditions for Human Occupancy.
Adaptive Comfort Model Based on European Standard EN152512007[LINK]
The EN152512007 is similar to ASHRAE 552010, but with slightly different curves of the indoor operative temperature and acceptability limits (Fig. 2). The model, intended for use in naturally ventilated buildings, determines the acceptability of indoor conditions given the 7day weighted mean outdoor air temperature and the indoor operative temperature. The 7day weighted mean outdoor air temperature (Trm) is defined as the weighted running average of the previous 7 daily average outdoor air temperatures.
This weighted running average is calculated from a full annual weather file that must be specified for the simulation. This is used as an index for occupant adaptation to outdoor conditions, and determines the acceptability of indoor conditions. The model also accounts for people’s clothing adaptation in naturally conditioned spaces by relating the acceptable range of indoor temperatures to the outdoor climate, so it is not necessary to estimate the clothing values for the space. No humidity or airspeed limits are required when this option is used. The model defines three comfort regions: Category I (90%) Acceptability, Category II (80%) Acceptability, and Category III (65%) Acceptability. If Trm is not within the specified domain of 10.0 to 30.0C, the model is not applicable.
Categories for European Standard EN152512007
Acceptable operative temperature ranges for naturally conditioned spaces (CEN EN152512007)
Central line (shown as redFigure): Tot = 0.33*To + 18.8
Category I, 90% Acceptability Limits: Tot = 0.33*To + 18.8± 2.0
Category II, 80% Acceptability Limits: Tot = 0.33*To + 18.8± 3.0
Category III, 65% Acceptability Limits: Tot = 0.33*To + 18.8± 4.0
For 10°(C) < Trm < 15°(C), the comfort temperature of the lower boundaries of the comfort regions is Tcomf = 23.75°(C). That is, the lower boundaries are constant according to the same ranges above:
Category I, 90% Acceptability Limits: Tot = 23.75 2.0
Category II, 80% Acceptability Limits: Tot = 23.75  3.0
Category III, 65% Acceptability Limits: Tot = 23.75 4.0
Where
Tot – operative temperature (°C) , calculated as the average of the indoor air drybulb temperature and the mean radiant temperature of zone inside surfaces
To – mean outdoor air drybulb temperature (°C), calculated as the weighted mean of the previous 7day daily mean outdoor air drybulb temperature (Tod):
To = (1  α)*{Tod−1 + α *Tod−2 + α2 *Tod−3 + α3 *Tod−4 + α4 *Tod−5 + α5 *Tod−6 + α6 *Tod−7}
To = (1  α)*Tod−1 + α *To−1
α = 0.8
Dynamic Clothing Model[LINK]
In most building energy simulations, thermal comfort condition is calculated based on the assumption that the clothing insulation is equal to a constant value of 0.5 clo during the cooling season and 1.0 clo during heating season. Usually those two values are used and the change from 0.5 to 1 or viceversa is made suddenly from one day to another. In addition, there is no standardized guideline on how to set clothing insulation schedules in the international standards. This simplified assumption may lead to systems that are incorrectly sized and operated and to the incorrect assessment of comfort conditions. In reality, occupants frequently adjust their clothing depending on the thermal conditions around them, as opposed to the assumption of constant clothing values. Therefore, the clothing insulation variation should be captured during the building simulation to realistically model HVAC systems. In order to overcome the limitations of the constant clothing insulation assumption, three new predictive clothing insulation models were developed by Schiavon and Lee (2012) based on 6,333 selected observations taken from ASHRAE RP884 and RP921 databases. The first and third models vary the clothing insulation as a function of outdoor air temperature measured at 6 o’clock and the second model takes into account both 6 o’clock outdoor air temperature and indoor operative temperature when adjusting the clothing insulation. The dynamic clothing models should be implemented in dynamic building energy simulation.
The model proposed to ASHRAE 55 is described below.
For ta(out,6) < 5°C, Icl=1 .00
For 5°C ≤ ta(out,6) < 5°C, Icl=0.818−0.0364∗ta(out,6)
For 5°C ≤ ta(out,6) < 26°C Icl=10(−0.1635−0.0066∗ta(out,6))
or ta(out,6) ≥ 26°C Icl=0.46
Where, Icl is the clothing insulation value, ta(out,6) is the outdoor air temperature measured at 6 o’clock in the morning. The following figure illustrates the proposed clothing insulation model.
Clothing Schedule Illustration
In the figure above, clothing insulation schedule for a fixed model (blue) typically used in energy simulation software and for the clothing model based on outdoor air temperature measured at 6 o’clock. Climate data for Chicago O’hare International Airport has been used.
The following figure illustrates the new clothing insulation model.
Graphical representation of the proposed clothing insulation model
The dynamic predictive clothing insulation model is implemented into EnergyPlus for realistic energy simulation. Addenda A to ASHRAE 55 with the clothing model has been approved by the ASHRAE committee and the chance to be included in ASHRAE 552013 is high.
References[LINK]
Schiavon S, Lee KH. 2013. Dynamic predictive clothing insulation models based on outdoor air and indoor operative temperatures. Building and Environment. Volume 59, 250260. http://dx.doi.org/10.1016/j.buildenv.2012.08.024 (link to the journal) http://escholarship.org/uc/item/3338m9qf (link to the freely available preprint version)
Lee KH, Schiavon S. 2013. Influence of three dynamic predictive clothing insulation models on building energy use, HVAC sizing and thermal comfort. Submitted to Journal. http://escholarship.org/uc/item/3sx6n876 (link to the freely available preprint version)
Mean Radiant Temperature Calculation[LINK]
There are three options to calculate mean radiant temperature in the thermal comfort models. One is the zone averaged MRT, another is the surface weighted MRT, and the other is angle factor MRT. The zone averaged MRT is calculated on the assumption that a person is in the center of a space, whereas the surface weighted MRT is calculated in consideration of the surface that a person is closest to, and the angle factor MRT is calculated based on angle factors between a person and the different surfaces in a space. Here, the surface weighted MRT is the average temperature of the selected surface and zone averaged MRT and is intended to represent conditions in the limit as a person gets closer and closer to a particular surface. In that limit, half of the person’s radiant field will be dominated by that surface and the other half will be exposed to the rest of the zone. Note that the surface weighted MRT is only an approximation. The angle factor MRT is the mean temperature of the surrounding surface temperatures weighted according to the magnitude of the respective angle factors and allows the user to more accurately predict thermal comfort at a particular location within a space.
Description of the model and algorithm[LINK]
The zone averaged MRT is calculated without weighting any surface temperature of the space.
Tr=Tsurf−1Fsurf−1+Tsurf−2Fsurf−2+.........+Tsurf−nFsurf−n
RadTemp = MRT(ZoneNum)
The surface weighted MRT is the average temperature of the zone averaged MRT and the temperature of the surface that a person is closest to.
Tr=(Tr−avg+Tsurf)/2
ZoneRadTemp = MRT(ZoneNum)
SurfaceTemp = GetSurfaceTemp(People(PeopleNum)%SurfacePtr)
RadTemp = (ZoneRadTemp + SurfaceTemp)/2.0
The angle factor MRT is the mean value of surrounding surface temperatures weighted by the size of the respective angle factors between a person and each surface.
Tr=Tsurf−1Fsurf−1+Tsurf−2Fsurf−2+.........+Tsurf−nFsurf−n
SurfTempAngleFacSummed = SurfTempAngleFacSummed &
RadTemp = SurfTempAngleFacSummed
References[LINK]
ASHRAE. 1984. “High Intensity Infrared Radiant Heating”, 1984 Handbook of Systems and Equipment, American Society of Heating, Refrigerating and Air Conditioning Engineers, Atlanta, GA, Chapter 18.
ASHRAE. 1985. “Physiological Principles for Comfort and Health,” 1985 Handbook of Fundamentals, American Society of Heating, Refrigerating and Air Conditioning Engineers, Atlanta, GA, Chapter 8.
ASHRAE. 1993. “Physiological Principles and Thermal Comfort”, 1993 ASHRAE Handbook of Fundamentals, American Society of Heating, Refrigerating and Air Conditioning Engineers, Atlanta, GA, Chapter 8.
ASHRAE. 2010. “Standard 552010 – Thermal Environmental Conditions for Human Occupancy (ANSI approved)”, American Society of Heating, Refrigerating and Air Conditioning Engineers, Atlanta, GA.
Azer, N.Z., Hsu, S. 1977. “The prediction of Thermal Sensation from Simple model of Human Physiological Regulatory Response”, ASHRAE Trans., Vol.83, Pt 1.
Berglund, Larry. 1978. “Mathematical Models for Predicting the Thermal Comfort Response of Building Occupants”, ASHRAE Trans., Vol.84.
Doherty, T.J., Arens, E. 1988. “Evaluation of the Physiological Bases of Thermal Comfort Models”, ASHRAE Trans., Vol.94, Pt 1.
Du Bois, D. and E.F. 1916. “A Formula to Estimate Approximate Surface Area, if Height and Weight are Known”, Archives of internal Medicine, Vol.17.
CEN. 2007. “Standard EN15251 Indoor environmental input parameters for design and assessment of energy performance of buildings addressing indoor air quality, thermal environment, lighting and acoustics”. Bruxelles: European committee for Standardisation.
Fanger, P.O. 1970. Thermal ComfortAnalysis and Applications in Environmental Engineering, Danish Technical Press, Copenhagen.
Fanger, P.O. 1986. “Radiation and Discomfort”, ASHRAE Journal. February 1986.
Fanger P.O. 1967. “Calculation of Thermal Comfort: Introduction of a Basic Comfort Equation”, ASHRE Trans., Vol.73, Pt 2.
Fountain, Marc.E., Huizenga, Charlie. 1997 “A Thermal Sensation Prediction Tool for Use by the Profession”, ASHRAE Trans., Vol.103, Pt 2.
Gagge, A.P., Stolwijk, J. A. J., Nishi, Y. 1970. “An Effective Temperature Scale Based on a Simple Model of Human Physiological Regulatory Response”, ASHRAE Trans., Vol.70, Pt 1.
Gagge, A.P., Fobelets, A.P., Berglund, L. G. 1986. “A Standard Predictive Index of Human Response to the Thermal Environment”, ASHRAE Trans., Vol.92, Pt 2.
Hsu, S. 1977. “A Thermoregulatory Model for Heat Acclimation and Some of its Application”, Ph. D. Dissertation, Kansas State University.
IntHout, D. 1990. “Thermal Comfort Calculation / A Computer Model”, ASHRAE Trans., Vol.96, Pt 1.
ISO. 1983. “Determination of the PMV and PPD Indices and Specification of the Conditions for Thermal Comfort”, DIS 7730, Moderate Thermal Environment, 1983.
Documentation content copyright © 19962017 The Board of Trustees of the University of Illinois and the Regents of the University of California through the Ernest Orlando Lawrence Berkeley National Laboratory. All rights reserved. EnergyPlus is a trademark of the US Department of Energy.
This documentation is made available under the EnergyPlus Open Source License v1.0.