calcTNZPDF: Values related to TNZ approach

In marcelschweiker/comfort_R: Models and Equations for Human Comfort Research

Values related to TNZ approach

Description

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.

Usage

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)

Arguments

ht	a numeric value presenting body height in [cm].
wt	a numeric value presenting body weight in [kg].
age	a numeric value presenting the age in [years].
gender	a numeric value presenting sex (female = 1, male = 2)
clo	a numeric value presenting clothing insulation level in [clo].
vel	a numeric value presenting air velocity in [m/s].
tskObs	a numeric value presenting actual mean skin temperature in [degree C].
taObs	a numeric value presenting actual air temperature in [degree C].
met	a numeric value presenting metabolic rate (activity related) in [met].
rh	a numeric value presenting realtive humidity in [%].
fBasMet	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.
fSA	a string presenting the method of calculating the surface area. Needs to be one of "duBois", "mosteller", or "direct".
percCov	a numeric value between 0 and 1 presenting the percentage of the body covered by clothes in [%].
TcMin	a numeric value presenting the minimum allowed core temperature in [degree C].
TcMax	a numeric value presenting the maximum allowed core temperature in [degree C].
plotZone	a boolean variable TRUE or FALSE stating, wether TNZ should be plotted or not.
gridTaMin	a numeric value defining the minimum grid value for Ta, ambient temperature, in [degree C].
gridTaMax	a numeric value defining the maximum grid value for Ta, ambient temperature, in [degree C].
gridTskMin	a numeric value defining the minimum grid value for Tsk, skin temperature, in [degree C].
gridTskMax	a numeric value defining the maximum grid value for Tsk, skin temperature, in [degree C].
gridTa	a numeric value defining the grid size in Ta dimension.
gridTsk	a numeric value defining the grid size in Tsk dimension.
sa	a numeric value for surface area (only used with method fSA: direct) in [m2].
IbMax	a numeric value for maximum body tissue insulation in [m2K/W].
IbMin	a numeric value for minimum body tissue insulation in [m2K/W].
alphaIn	a numeric value for alpha (if 0, alpha will be calculated according to Fanger.
metMin	a numeric value for minimum metabolic rate (only used with method fBasMet:direct) in [W/m2].
metMax	a numeric value for maximum metabolic rate (only used with method fBasMet:direct) in [W/m2].
metDiff	a numeric value for difference between minimum and maximum metabolic rate (not used with method fBasMet:direct) in [W/m2].
forPDF	a boolean value. If TRUE, matrix for drawing of PDF will be output, if FALSE, values for dTNZ and others will be output.
metAdapt	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.
trm	numerical value presenting the running mean outdoor temperature in [degree C]. Only used with metAdapt: Hori and ATHB.
TcPreAdapt	numerical value presenting the initial core temperature before adaptation in [degree C]. Only used with metAdapt: Q10.

Details

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.

Value

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

dTNZ	The absolute distance to the centroid of the thermoneutral zone
dTNZTs	Relative value of distance assuming skin temperature to be dominant for sensation
dTNZTa	Relative value of distance assuming ambient temperature to be dominant for sensation

Note

This function was used for the review paper by Schweiker et al. (2018) (see reference above). Some of the equations implemented are still to be validated further - therefore, use this function and its parameters with great care.

This function is not (yet) implemented in calcComfInd, calcdTNZ is applied there.

Author(s)

Marcel Schweiker and Boris Kingma

References

Schweiker, Huebner, Kingma, Kramer & Pallubinsky (2018) <doi:10.1080/23328940.2018.1534490> Kingma, Schweiker, Wagner & van Marken Lichtenbelt 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 & van Marken Lichtenbelt (2015) <doi:10.1038/nclimate2741> Kingma, Frijns, Schellen & van Marken Lichtenbelt (2014) <doi:10.4161/temp.29702>

See Also

calcdTNZ

Examples

## 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]")

marcelschweiker/comfort_R documentation built on Feb. 22, 2024, 9:04 p.m.