calcTNZPDF calculates the distance from the thermoneutral zone, either skin temperature or room air related. Also calculates the probability function (PDF) of the thermoneutral zone.
calcTNZPDF(ht, wt, age, gender, clo, vel, tskObs, taObs, met, rh,
fBasMet = "rosa", fSA = "duBois", percCov = 0, TcMin = 36, TcMax = 38,
plotZone = FALSE, gridTaMin = 20, gridTaMax = 30, gridTskMin = 30, gridTskMax = 42,
gridTa = 1000, gridTsk = 1000,
sa = 1.86, IbMax = 0.124, IbMin = 0.03, alphaIn = 0.08, metMin = 55.3, metMax = 57.3,
metDiff = .1, forPDF = FALSE, metAdapt = "none", trm = 15, TcPreAdapt = 37.2)a numeric value presenting body height in [cm]
a numeric value presenting body weight in [kg]
a numeric value presenting the age in [years]
a numeric value presenting sex (female = 1, male = 2)
a numeric value presenting clothing insulation level in [clo]
a numeric value presenting air velocity in [m/s]
a numeric value presenting actual mean skin temperature in [degree C]
a numeric value presenting actual air temperature in [degree C]
a numeric value presenting metabolic rate (activity related) in [met]
a numeric value presenting realtive humidity in [
a string presenting the method of calculating basal metbolic rate. Needs to be one of "rosa", "harris", "miflin", "fixed", or "direct". Fixed will result in the value of 58.2 W/m2. Direct requires definition of metMin and metMax.
a string presenting the method of calculating the surface area. Needs to be one of "duBois", "mosteller", or "direct".
a numeric value between 0 and 1 presenting the percentage of the body covered by clothes in [
a numeric value presenting the minimum allowed core temperature in [degree C].
a numeric value presenting the maximum allowed core temperature in [degree C].
a boolean variable TRUE or FALSE stating, wether TNZ should be plotted or not.
a numeric value defining the grid size in Ta dimension.
a numeric value defining the grid size in Tsk dimension.
a numeric value defining the minimum grid value for Ta, ambient temperature, in [degree C].
a numeric value defining the maximum grid value for Ta, ambient temperature, in [degree C].
a numeric value defining the minimum grid value for Tsk, skin temperature, in [degree C].
a numeric value defining the maximum grid value for Tsk, skin temperature, in [degree C].
a numeric value for surface area (only used with method fSA: direct) in [m2].
a numeric value for minimum body tissue insulation in [m2K/W].
a numeric value for maximum body tissue insulation in [m2K/W].
a numeric value for alpha (if 0, alpha will be calculated according to Fanger.
a numeric value for minimum metabolic rate (only used with method fBasMet:direct) in [W/m2].
a numeric value for maximum metabolic rate (only used with method fBasMet:direct) in [W/m2].
a numeric value for difference between minimum and maximum metabolic rate (not used with method fBasMet:direct) in [ ].
a boolean value. If TRUE, matrix for drawing of PDF will be output, if FALSE, values for dTNZ and others will be output.
a string presenting the method of calculating the surface area. Needs to be one of 'Hori', 'Q10', 'ATHB', or 'none'. NOTE: all methods applied here still in development and need further validation.
numerical value presenting the running mean outdoor temperature in [degree C]. Only used with metAdapt: Hori and ATHB.
numerical value presenting the initial core temperature before adaptation in [degree C]. Only used with metAdapt: Q10.
calcTNZPDF returns either a dataframe suitbale to draw the pdf of TNZ (by setting forPDF to TURE) or a dataframe with the columns dTNZ, dTNZTs, dTNZTa and others. Thereby
The absolute distance to the centroid of the thermoneutral zone
Relative value of distance assuming skin temperature to be dominant for sensation
Relative value of distance assuming ambient temperature to be dominant for sensation
The percentage of the body covered by clothes can be estimated e.g. based on ISO 9920 Appendix H (Figure H.1). A typical winter case leads to a value of around .86, in the summer case this goes down to values around .68.
Schweiker, M., Huebner, G. M., Kingma, B. R. M., Kramer, R., and Pallubinsky, H. Drivers of diversity in human thermal perception - A review for holistic comfort models. Temperature, 2018, 1 - 35.
Kingma, B. R., Schweiker, M., Wagner, A. and van Marken Lichtenbelt, W. D. Exploring the potential of a biophysical model to understand thermal sensation Proceedings of 9th Windsor Conference: Making Comfort Relevant Cumberland Lodge, Windsor, UK, 2016.
Kingma, B. and van Marken Lichtenbelt, W. Energy consumption in buildings and female thermal demand Nature. Clim. Change, 2015, 5, 1054 - 1056.
Kingma, B. R.; Frijns, A. J.; Schellen, L. and van Marken Lichtenbelt, W. D. Beyond the classic thermoneutral zone. Temperature, 2014, 1, 142 - 149.
see also calcdTNZ
# NOT RUN {
## Calculate and draw pdf of TNZ for a young non-obese male
longTcYoungMale <- calcTNZPDF(ht = 178, wt = 70, age = 30, gender = 2, clo = 0.5,
vel = 0.2, tskObs = 36.2, taObs = 26, met = 1,
rh = 50, fBasMet = "rosa", fSA = "duBois", percCov = 0.6,
TcMin = 36, TcMax = 38, plotZone = FALSE, gridTaMin = 20, gridTaMax = 30,
gridTskMin = 20, gridTskMax = 42, gridTa = 1000, gridTsk = 1000, sa = 2.0335, IbMax = 0.124,
IbMin = 0.03, alphaIn = 0, metMin = 55.3, metMax = 57.3, metDiff = 0.1, forPDF = TRUE)
plot(density(longTcYoungMale$X2), main="", xlim=c(14,36), ylim=c(0,.50),
xlab="Operative temperature [degree C]")
# }
Run the code above in your browser using DataLab