AP06_6.vp 1 Introduction The provision of optimal hygrothermal conditions, i.e. above all an optimal operative temperature (calm air and air temperature reaching radiant temperature) is the principal condition for healthy human life in the interior of a building. The optimal operative temperature has in the past been calculated from the PMV (Predicted Mean Value) (see e.g. EN ISO 7730 Moderate Thermal Environment) estimated on the basis of a positive reaction from 80% of the persons presents. The feelings of human beings are very subjective values, impacted by many other factors in addition to hygro- thermal conditions, e.g. by indoor interior colors, a person’s mood, etc. In addition, due to the way in which PMV is experi- mentally estimated and proved in other experimental works (see Fishman, Pimbert 1979, Newsham, Tiller 1995), it is ap- proximately valid for the neutral zone only. The further away from the neutral zone the more the real values depart from the values calculated from PMV, see Fig. 1. What is more, the greater a person’s activity, the bigger the difference. The ap- plication of high activity values is thus impossible in practice. In Fig. 2 the mean thermal sensation vote is plotted against the operative temperature for a range of velocities. Each point represents the mean vote of thirty two subjects. The correlation between the operative temperature and the mean thermal sensation vote is high, with a correlation coeffi- cient of 0.97 (n � 5). There is no significant difference between the sexes. The solid curve is the regression line for the individ- ual vote (n � 80). For comparison, the dotted line represents © Czech Technical University Publishing House http://ctn.cvut.cz/ap/ 3 Acta Polytechnica Vol. 46 No. 6/2006 Optimal (Comfortable) Operative Temperature Estimation Based on Physiological Responses of the Human Organism M. V. Jokl, K. Kabele Problems following the application of optimal operative temperatures estimated on the basis of PMV and the necessity to apply correct values in the new Czech Government Directive No. 523/2002 Code led to experiments based on the physiological human body response instead of solely people’s feelings in a given environment. On the basis of experiments on 32 subjects (university students) it has been possible to estimate: a) the total balance of hygrothermal flows between the human body and the environment, b) the optimal operative temperature as a function of the subject’s activity, c) the thermoregulatory range for each optimal operative temperature, i.e. maximal (category Cmax) limited by the onset of sweating, minimal (category Cmin) limited by the onset of shivering (category C can be applied to naturally ventilated buildings), optimal (comfort level – category A) defined by time constant 0.368 (can be applied to air conditioned buildings), and submaximum (decreased comfort level – category B) defined by time constant 0.632 (can be applied to buildings with basic air conditioning systems). Keywords: thermal comfort, microenvironment, hygienic regulations, PMV problems, thermoregulatory ranges. 16 18 20 22 24 26 28 30 �3 �2 �1 0 1 2 3 indoor temperature (°C) co m fo rt vo te Observed vote Predicted vote a) Fig. 1: a) Comparison of mean thermal comfort votes (ASHRAE scale) with predictions by the PMV model in an English office building (Fishman and Pimbert, 1979), activity 80 W�m�2, clothing 0.64 up to 0.82 clo; b) Comparison of mean thermal comfort votes (ASHRAE scale) with predictions by the PMV model in a build- ing (Newsham and Tiller, 1995), activity 70 W�m�2, clothing 0.78�0.21 clo. �3 �2 �1 0 1 2 3 1918 20 22 24 2623 25 C o m fo rt vo te Operative temperature (°C) 21 ASHRAE PMV b) the results for 172 Japanese subjects in conditions of low air movement reported by Tanabe (1987), and the dashed line represents the calculated PMV values. Further important results are drawn from research at Reading University (Croome et al.1993). These results take into account the opening and closing of doors, i.e. the ventila- tion rate, see Fig.3. When the windows and the door were closed, the mean thermal sensation tended to be on the warm side of neutral. When the windows and doors were open, the votes were spread widely over the thermal sensation scale. However, the calculated PMV values corresponding to the tests were close to the neutral point for most of the test conditions. This suggests that, in this investigation, PMV underestimates the thermal impression for the case when the windows and doors are shut, and undervalues the change in thermal impression for the two cases. This may be due to three main reasons. The first reason is the assumed steady state laboratory conditions used in the derivation of the PMV equation. The second is the oversim- plified approach to the assessment of the metabolic rate of the occupant. The occupants rarely sat in the room for a long period, say one hour, without moving around. The third rea- son is the sensitivity of PMV to clo values (Croome et al. 1993). It can be concluded that the PMV equation overpredicts the neutral temperature by as much as 2 K and underpredicts the comfort requirement when air temperature deviates from neutrality. Humphreys and Nicol (2000) have suggested that there may be formulaic errors in such a complex index as PMV, with two contributing factors: 1) Steady state approximation. PMV, like other indices of warmth, is a steady state heat exchange equation, and therefore its application to the office environment can only be an approximation. Recent research shows that among an office population the temperature of the fin- gers varies extensively and rapidly, indicating that the thermal state of the bodies of office workers is in continual flux (Humphreys et al. 1999). This suggests that it is better to regard the people a soon as being in dynamic thermal equilibrium rather than in a steady thermal state. By extension, the same is likely to be true of other and more varied pursuits. Thus, any index built on steady state assumptions is of limited relevance to normal living. Indices that exlude thermoregulation cannot therefore simulate real life conditions. 2) Inaccurate numerical formulae for steady state. Most indices have errors in the numerical values used in the equa- tions, such as the convective and radiant heat transfer coefficients, skin temperature and sweat rates that are assumed in comfort conditions. These contribute to for- mulaic errors and additionally there are numerical errors attributable to conceptual simplifications. For example, although the calculation of PMV is based on calculated skin temperature and sweat rate but when considering external conditions to be neutral, PMV is based solely on a hypothetical heat load. This results in the same body thermal states being attributed different PMV values in different environments (Humphreys and Nicol 1996). Conceptual and numerical approximations add to the formulaic error. And, no thermoregulatory ranges can be estimated, from the PMV system. For these reasons we decided to estimate the optimal operative temperatures on the basis of the physiological re- sponse of the human organism. 4 © Czech Technical University Publishing House http://ctn.cvut.cz/ap/ Acta Polytechnica Vol. 46 No. 6/2006 Mean air velocity m/s[ ] PRESENT RESULT TANABE et al. 1987 PMV (CALCULATION) 0.52 0.29 0.22 0.16 0.1 �3 �2 �1 0 1 2 3 25 26 27 Operative temperature °C[ ] M e a n th e rm a ls e n sa tio n vo te A S H R A E Fig. 2: Mean thermal sensation vote versus operative tempera- ture for Japanese college-age subjects (Tanabe et al. 1987) 16 18 20 22 24 26 28 22.8 °C �2 �1 0 1 2 3 T h e rm a l s e n s a ti o n vo te A S H R A E a) for windows and the door closed Operative temperature [ ]°C Fo ot lev el (r= 0. 67 ) (r= 0.6 2) He ad lev el (r= 0.5 2) PM V Ov er all 22.8 °C �2 �1 0 1 2 3 16 18 20 22 24 26 28 �3T h e rm a l s e n s a ti o n vo te A S H R A E b) for windows and the door open Operative temperature [ ]°C Fo ot le ve l (r= 0. 62 ) (r= 0. 72 ) He ad le ve l ( r= 0. 70 ) PMV O ve ra ll Fig. 3: Effect of operative temperature on thermal sensation vote, activity 1.2 met, clo value results from neutral tempera- ture 22.8 °C, rh � 40�55 %, mean radiant temperature equals air temperature (Croome et al. 1993) 2 Mathematical model of the physiological body response The total heat rate production and its distribution into individual components during heat exchange between the human body and the environment are shown in Fig. 4, where qm � M � W � metabolic heat (see Jokl 1989). qres and qev, d are the components of the heat rate from the organism due to respiration and due to skin moistening (evaporation), when the human body is in the thermal neutral zone. The heat flow qdry represents the component transferred from the organism through the clothing layer with total thermal resistance Rt, wa(qdry � qc� qr). The regulatory process with- in the neutral zone is achieved mainly by vasodilation and vasoconstriction changing the body’s internal resistance into thermoregulatory and adaptational heat flux qtr and qa to the skin surface. qtr and qa is the heat flux regulating the instanta- neous value of the skin temperature during the subject’s interaction with the environment, qtr is the organism’s imme- diate response to changes in the microclimate or metabolic heat changes; qa is the reaction shift due to adaptation to heat in summer and cold in winter. qtr� qa may be negative (heat loss) or positive (heat gain). It is the transient heat flow – even in the thermal neutral zone – that is called “quasi-stationary”, to be distinct from the hyperthermia and hypothermia zone. qtr� qa represents the rates of heat storage or heat debt accumulation. When the body is in a steady-state thermal bal- ance with the environment, these therms are equal to zero. However we can consider the state of the subject in the neu- tral zone by non-steady-state conditions due to periodical changes of the metabolic heat rate, qm, or short thermal exci- tations in time followed by changes in the internal thermal resistance of the body within the neutral zone. The temporary characteristics of each non-steady process are determined, in addition to the thermal resistances Rt, i and Rt, wa, by the human body heat capacity, Ct. The values characterizing the heat exchange are: Tsk, Tcore and Tg. The internal thermal resistance, Rt, i ,also determines the changes in thermoregulation and the adaptational heat, qtr+ qa, which is necessary for maintaining the skin temperature within physiological values if the core temperature is to remain constant (Tcore � 36.7 � 0.4 °C). The heat flow balance, as presented in the model shown in Fig. 4, can be expressed by a thermal flux equation at the subject-environment boundary. Thus (if heat conduction is neglected): q R T T q q q q q q q dry t, wa g sk m res ev tr a i sw W m � � � � � � � � � � � 1 ( ) [ 2 ] (1) where qev � qev, ins � qev, sens � qev, ins� qsw [W � m �2], qm� q res� qev, ins � qi [W � m �2], qsw � 0.6(qm� 58.14) [W � m �2], � the quantity of excreted perceptible but mostly invisible sweat. This was estimated by weighing during the experi- ments as a mean value for the whole range. Heat flux within the human body can be represented as (see model in Fig. 4): q q q q G T T R T T m res tr a t,ti i sk t,ti i sk W m � � � � � � � � � ( ) ( ) [ ] 1 2 (2) where Gt, ti is total body thermal conductance, which can be expressed by Eq. (3): G q q T T q q T T G Gt,ti m res i sk tr a i sk t,m t,i� � � � � � � � [W � m�2� K�1] (3) where Gt, i is internal thermal conductance and Gt, m is meta- bolic thermal conductance. © Czech Technical University Publishing House http://ctn.cvut.cz/ap/ 5 Acta Polytechnica Vol. 46 No. 6/2006 qres Rt,wa Rt,wa Rt,i (Tsk ) CtTcore T Ti core� Tsk Tg qev qr qc qm, Rt,m METABOLIC HEAT PRODUCTION qres�qres CORE THERMOREGULATION Fig. 4: Total heat rate production and its distribution in individual components during heat exchange between the human body and the environment (qm metabolic heat, qres respiration heat, qtr thermoregulatory heat, qe evaporative heat, qc convective heat, qr radi- ant heat, Rt,wa total thermal resistance of clothing, Rt total internal thermal body resistance, Ct thermal body capacity, Ti deep body temperature, Tcore core body temperature, Tsk skin temperature, Tg globe temperature 6 © Czech Technical University Publishing House http://ctn.cvut.cz/ap/ Acta Polytechnica Vol. 46 No. 6/2006 q q � i s w [ ] W /m 2 qdry[] W/m 2 F ig . 5: G ra p h of th e re la ti on sh ip q d ry � f( q i � q s w ) fo r cl ot h in g 0. 5 cl o, p oi n ts fr om ex p er im en t, op ti m al va lu es ar e on th e li n e q d ry � q i � q s w © Czech Technical University Publishing House http://ctn.cvut.cz/ap/ 7 Acta Polytechnica Vol. 46 No. 6/2006 0 .5 c lo � � � � � q (q q ) d ry 4 .8 i sw q m / [W m ] 2 T[] g°C T g 0 .1 0 1 7 q 3 1 .7 0 8 R 0 .8 8 6 3 � � � � m 2 0 1 0 2 0 3 0 4 0 5 0 6 0 7 0 8 0 9 0 1 0 0 1 1 0 1 2 0 1 3 0 1 4 0 1 5 0 1 6 0 1 7 0 1 8 0 1 9 0 2 0 0 2 1 0 2 2 0 2 3 0 2 4 0 2 5 0 0 2 .55 7 .51 0 1 2 .51 5 1 7 .52 0 2 2 .52 5 3 0 2 7 .5 3 2 .53 5 F ig . 6: G ra p h of th e re la ti on sh ip T g � f( q m ) fo r op ti m al va lu es tr an sf er re d fr om gr ap h on F ig . 5 w it h in th e ra n ge q d ry � (q i� q s w )� � 4. 8, w h er e th e va lu e � 4. 8 W �m � 2 re p re se n ts th e m in im al th er m or eg u la to ry h ea t Thermoregulation and adaptational heat flux first affects the skin temperature, Tsk. The internal thermal resistance value, Rt, i � 1/Gt, i, characterizing the vasodilatation and vaso- constriction process, can be calculated from the equation: R T T q qt, i i sk tr a � � � [W�1� m2� K] (4) 3 Experimental estimation of mathematical model parameters An experiment over the course of several years was under- taken in a climatic chamber, leading identification of the parameters in Eq. (1) and (4). The experimental subjects were male university students. Each of them underwent six experiments lasting about three hours at four levels of activity: (1) sitting in a chair, (2) sitting on a bike-ergometer without pedaling, (3) pedaling on a bike-ergometer with a 40 W load and (4) pedaling on a bike-ergometer with a load of 1 W per kg body mass (for as long as he was able to continue). Metabolic heat production during each activity was measured by an indirect calorimetric method. Mean skin temperature, heat rate and body water loss were estimated continuously during each experiment. Two sets of clothing were used by the subjects: lightweight (pyjamas) and heavier clothing (an anti-g suit for fighter pilots). The results of the anti-g suit experiments will be pre- sented in a separate report. There were no differences between the air temperature and the surface wall temperatures – it can be assumed that the overall temperature equals the operative temperature. Six temperatures were chosen (29�3 °C and 14�3 °C, which determine the temperature ranges where some of the sub- jects started to leave the neutral zone and appeared to begin sweating or shivering). The originally chosen range of tem- peratures 8, 11, 14, 17, 20, 23, 26, 29, 32 °C was found to be excessive, and so they were reduced. Within the comfort range the relative humidity was main- tained corresponding to a partial water vapour pressure from 700 to 1850 Pa. The onset of sweating and shivering was always assessed by the same person. Experiments were car- ried out in all seasons of the year, thus reflecting the seasonal adaptation effect on maximal and minimal thermoregulatory heat, i.e. it was possible to determine adaptational heat. How- ever it became evident that the seasonal adaptation effect can be neglected (Jokl, Moos 1992), being lower than 0.2 °C (i. e. within the range of experimental error in measuring the temperatures). The same finding has been described by other authors (see Fanger 1970). The results were only ac- cepted from subjects within the thermal neutral zone with the thermoregulatory heat constant. 3.1 Graph construction of Tg � f(qm) The measured values are plotted as qdry � f ( qi� qsw) in Fig. 5, where for optimal values the linear equation: qdry � qi� qsw [W�m �2] representing equilibrium is valid. Practi- cal application of this graph is very difficult; but the relation- ship Tg � f ( qm) is useful. Therefore the linear relationship from Fig. 5 was transferred into the graph in Fig. 6 by plotting a regression line through the points limited by the equation � qdry� (qi� qsw) � � 4 .8 [W�m �2] in Fig. 6. The value of � 4 .8 W�m�2 of the regression line is the mini- mal thermoregulatory heat, i.e. it represents maximal vaso- constriction in human body, and can be estimated from the minimum value for internal thermal conductivity of the hu- man body (see Fig. 7), which equals 9.07 W�m�2�K�1 (for core body temperature Ti � 36.6 °C, skin temperature Tsk � 30.5 °C and qm � 45.7 W�m �2). q G T Ttr, t,i, i sk W mmin min ( ) . ( . . ) .� � � � � � �157 36 6 30 5 9 6 2 , i.e. �4.8 W�m�2, where G G Gtt,i, t, i, t,m,min min min . . .� � � � �9 07 7 5 1 57 W�m �2 ��, G q T Tt,m, m i sk min . . . .� � � � � 45 7 36 6 30 5 7 5 W�m�2��� 3.2 Estimation of thermoregulatory range The widest thermoregulatory range, i.e. from optimum up to the onset of visible sweating can be estimated by plotting the regression line into the points of the onset of sweating. However, to ensure comfort we need lower values, without vis- ible sweating occurring. This area is between the line of the optimum value and the tangent from the origin (which is the intersection of the line of the optimum and the regression line of the onset of sweating) to the area of beginning of shivering (see Fig. 8). These tangents are analogous to the thermo- regulatory range of category C according to CR 1752-1998. For categories A and B it must be taken into consideration that the human body is a thermoregulatory mechanism in the surrounding environment balancing the operative tempera- ture changes by thermoregulatory heat flows in the human so that equilibrium can be achieved, and this must take place at three levels (by analogy with technological mechanisms) (see Fig. 9): 8 © Czech Technical University Publishing House http://ctn.cvut.cz/ap/ Acta Polytechnica Vol. 46 No. 6/2006 0 10 20 30 40 50 60 12 16 20 24 28 32 36 0 50 100 150 200 H e a t L o s s W m [ ] � � 2 Skin Temperature °C[ ] in water in air conductance heat loss C o n d u c ta n c e W m K [ ] � � � 2 � 1 m in . Fig. 7: Human body thermal conductance Gt, ti and human body heat loss as a functionof skin temperature Tsk for a resting subject during daytime (qm � 45.6–57.4 W�m �2) (Burton, Bazett 1936, Du Bois et al. 1952, Lefevre 1898, Liebmaster 1869) (from Itoh et al. 1972) � level A, corresponding to the time constant 0.368 To, tr, max � level B, corresponding to the time constant 0.632 To, tr, max � level C, corresponding to the time constant 1.000 To, tr, max Level A is valid for building interiors with the highest requirements, and can only be attained with the use of air conditioning systems. Level C is valid for building interiors with the lowest requirements, usually only naturally venti- lated. Level B covers other buildings, where air conditioning is necessary only in some cases. The time constant according to control theory charac- terizes the system response, (the response of the human © Czech Technical University Publishing House http://ctn.cvut.cz/ap/ 9 Acta Polytechnica Vol. 46 No. 6/2006 Fig. 8: How to obtain thermoregulatory ranges (see the text for an explanation) Fig. 9: Graph of the relationship Tg � f (qm) with the regression line of onset of sweating and the thermoregulatory range for levels (cate- gories) A, B, C for warm (towards the onset of sweating) and for cold (towards the onset of shivering) organism) to the operative temperature changes, and is equal to the product of system thermal resistance R and its thermal capacity C: time constant � R�C, where R � Rti� Rtwa [W �1 �m2�K]. Thermoregulatory changes are shown in Fig. 10 as trans- ferred from Fig. 9 and also rounded up to 0.5 °C for practical application. For a complete list of values, see Table 1. 10 © Czech Technical University Publishing House http://ctn.cvut.cz/ap/ Acta Polytechnica Vol. 46 No. 6/2006 Thermoregulatory range (0.5 clo) �30 �20 �10 0 10 20 30 40 0 50 100 150 200 q tr [W /m 2 ] max C - opt max C - opt (2.5) max B - opt max B - opt (2 5). max A - opt max A - opt (2 5). min A - opt min A - opt (2 5). min B- opt min B - opt (2 5). min C - opt min C - opt (2 5). q [W/mm 2 ] T °C g ,t r [ ] �6.0 �4.0 �2.0 0 2.0 4.0 6.0 8.0 Fig. 10: Thermoregulatory changes as transferred from graph on Fig. 7 and rounded up to 0.5 °C for practical application (see also Table 1) °C 50 [W�m�2] 70 [W�m�2] 80 [W�m�2] 100 [W�m�2] 120 [W�m�2] 150 [W�m�2] 180 [W�m�2] sweat 32.2 31.4 31.1 30.3 29.6 28.5 27.4 sweat – opt 5.6 6.8 7.5 8.8 10.1 12.0 14.0 sweat – opt (0.5) 5.5 6.5 7 8.5 10 12 13.5 Max 29.3 27.9 27.2 25.8 24.4 22.3 20.2 max C – opt 2.7 3.4 3.7 4.3 4.9 5.9 6.8 max C – opt (0.5) 2.5 3.0 3.5 4.0 4.5 5.5 6.5 max B (0.632) 28.3 26.7 25.9 24.3 22.6 20.2 17.7 max B – opt 1.7 2.1 2.3 2.7 3.1 3.7 4.3 max B – opt (0.5) 1.5 2.0 2.0 2.5 3.0 3.5 4.0 max A (0.368) 27.6 25.8 24.9 23.1 21.3 18.6 15.9 max A – opt 1.0 1.2 1.4 1.6 1.8 2.2 2.5 max A – opt (0.5) 1.0 1.0 1.0 1.5 1.5 2 2.5 opt 26.6 24.6 23.6 21.5 19.5 16.5 13.4 Table 1: Optimal operative temperatures and thermoregulatory range as a function of man’s activity qm 4 A comparison between optimal values and thermoregulatory ranges with accepted values The proposed optimal temperatures and their ther- moregulatory ranges were compared with the values according to ISO 7730 (Moderate thermal environments ISO 7730-1984 (E)), CR (1752) (1998) and ISO/DIS 7730 (2003), ANSI/ASHRAE Standard 55-2004. The comparison between the above proposed opera- tive temperatures and the values according to CR 1752 and ISO/DIS 7730 is presented in Table 2 and in Fig. 11. There is agreement on operative temperatures for 50 W�m�2, 70 W�m�2 and 80 W�m�2. For higher activities the values differ: the greater the activity, the greater the operative tem- perature difference. The findings are in agreement with the experiments (sitting persons in the neutral zone) on which the PMV value is based. A comparison between proposed optimal operative tem- peratures and the values according to ISO and ANSI/ASHRAE is presented in Table 3 and Fig.12. There is evident agreement for low activities (the graph is also based on ISO 7730). 5 Discussion The optimal operative temperatures derived from PVM values (from the 1970’s) are now not fully acceptable. It is more precise to use optimal operative temperatures based on the physiological human body response and not based only on people’s feelings. This has been proved by experimental © Czech Technical University Publishing House http://ctn.cvut.cz/ap/ 11 Acta Polytechnica Vol. 46 No. 6/2006 °C 50 [W�m�2] 70 [W�m�2] 80 [W�m�2] 100 [W�m�2] 120 [W�m�2] 150 [W�m�2] 180 [W�m�2] min A (0.368) 25.9 23.7 22.6 20.4 18.2 14.9 11.6 min A – opt -0.7 �0.9 �1 �1.1 �1.3 �1.5 �1.8 min A – opt (0.5) -0.5 �0.5 �0.5 �1.0 �1.0 �1.5 �1.5 min B (0.632) 25.4 23.1 21.9 19.6 17.3 13.8 10.4 min B – opt -1.2 �1.5 �1.6 �1.9 �2.2 �2.6 �3.1 min B – opt (0.5) -1.0 �1.0 �1.5 �1.5 �2.0 �2.5 �3.0 min 24.7 22.2 21 18.5 16 12.3 8.6 min C – opt �1.9 �2.4 �2.6 �3.0 �3.5 �4.2 �4.8 min C – opt (0.5) �1.5 �2.0 �2.5 �3.0 �3.0 �4.0 �4.5 0.1 0.2 m °C/W 2 � 150 125 100 75 50 1.510.50 1.0 2.0 3.0 0 28 °C 26 °C 24 °C 22 °C 20 °C 18 °C 14 °C 16 °C 12 °C 10 °C �1 °C �2 °C �1.5 °C �2.5 °C �4 °C �3 °C �5 °C optimal �1.0 +1.5 +2.0 �1.0 +3.5 �2.5 +4.0 �3.0 +2.5 �1.5 +3.0 �2.5 26.6 °C 24.6 °C 23.6 °C 21.5 °C 19.5 °C 16.5 °C 13.4 °C +2.0 �1.5 W/m 2 met A c ti vi ty Clothing clo Fig. 11: Comparison between experimentally found optimal operative temperatures and the values in ISO 7730 (the values presented in the graph correspond to temperatures only, not to activity) works (Fishman, Pimbert 1979) and shown when ISO values have been applied in practice. The greater a person’s activ- ity, the greater the discrepancy in the optimal temperature. Because of this discrepancy, the new Czech Government Di- rective No. 523/2002 Code is based on the values presented here, and not on ISO/DIS 7730, which is based on PMV. The absence in the directive of adaptation to heat and cold, e.g. as a result of staying in a heated room in winter and in air- -conditioned cars in summer, results in the same optimal operative temperatures for winter and for summer; the tem- peratures are differentiated only by different clothing. 12 © Czech Technical University Publishing House http://ctn.cvut.cz/ap/ Acta Polytechnica Vol. 46 No. 6/2006 Category CR � ISO/DIS Jo/Ka A (air conditioning) 24.5 � 1.0 24.6 � 1.0 � 0.5 B (air conditioning and natural ventilation) 24.5 � 1.5 24.6 � 2.0 � 1.0 C (natural ventilation) 24.5 � 2.5 24.6 � 3.0 � 2.0 Table 2: A comparison between experimentally found operative temperatures and the values in CR 1752 and ISO/DIS 7730 (2003) in categories A, B, C (0.5 clo, 1.2 met) 30 20 10 % RH 40 50 60 70 80 90 Operative Temperature [°C] H u m id it y R a ti o D e w P o in t Te m p e ra tu re [°C ] Data based on ISO 7730 and ASHRAE STD 55 0.5 Clo 1.0 Clo PMW Limits 0.5 Upper Recomended Humidity Limit, 0.012 humidity ratio 10 13 16 18 21 24 27 3229 35 38 �10 �5 0 5 10 20 15 0.002 0.004 0.006 0.008 0.010 0.012 0.014 0.016 24. 0 u p t o 2 5.8 °C No Recomended Lower Humidity Limit Fig. 12: Comparison between experimentally found optimal operative temperatures and the values in ANSI/ASHRAE qm [W�m�2] 50 70 80 100 120 150 180 [met] 0.86 1.20 1.38 1.72 2.07 2.59 3.10 Jo/Ka [°C] (B) 26.6 �1.5 24.6 �2.0 23.6 �2.0 21.5 �2.5 19.5 �3.0 16.5 �3.5 13.4 �4.0 �1.0 �1.0 �1.5 �1.5 �2.0 �2.5 �3.0 ISO 7730 (1984) 26.6�1.5 24.5�1.5 23.6�2.0 22.3�2.0 20.6�2.5 18.5�2.5 16.4�2.5 CR 1752, ISO/DIS 7730 (2003) [°C] (B) – 24.5�1.5 23.5�2.0 – – – – ANSI/ASHRAE 55 (1992) To active � To sedentary� 4.5 (met � 1.2) [°C] – 24.5 23.7 22.2 20.6 18.2 16.0 Remark: (B) means category Table 3: A comparison between experimentally found operative temperatures and the values in ISO and ANSI/ASHRAE (clothing 0.5) 6 Results The mathematical model (Fig. 4) shows the role of various heat flows produced by the human body as it interacts with the environment. All the heat flows must be in mutual equilib- rium if the human body is to stay homoiotherm. Practical application of this graph is very difficult, and Tg � f (qm) is a more useful relationship (Fig. 6). This equilibrium forms the basis for optimal operative temperature estimation (Fig. 5). The experimental data on the onset of sweating and the onset of shivering enable the thermoregulatory ranges to be esti- mated (Fig. 8). The thermoregulatory area is between the line of the optimum and the tangent from pole, defined as the intersection of the line of optimum and the regression line of the onset of sweating, to the field of the onset of sweating (up- per limit, level Cmax) and to the field of the onset of shivering (lower limit, Cmin). It is interesting and in agreement with human feelings that the thermoregulatory field for cold is smaller than the thermoregulatory field of the warm area – the human body is more sensitive to temperature decreases in the cold area. The question is how to sub-divide the thermoregulatory range into categories (A, B and C). Instead of qualified assumption, it is proposed to base the categories on control theory. The human body behaves like any other system to which control theory can be applied. It is proposed that the human body time constant is used to differentiate the categories. The following values were used: time constant 0.368 T0, tr, max (A), 0.632 T0, tr, max (B) and 1.0 T0, tr, max, which correspond to categories A, B and C (Fig. 6). Category A can be applied to air conditioned buildings and category C to naturally ventilated buildings. As a result, two previously separate methods of assessment can be merged, those for air conditioned buildings, based on PMV, and these for natural ventilated buildings, based on mean monthly outdoor tem- perature. The results have been compared with the values according to ISO 7730, CR (1752) (1998), ISO/DIS 7730 (2003) (Table 2 and 3), and with the ANSI/ASHRAE Standard (Table 3). Most importantly however, it was possible to base the new Czech Government Directive No. 523/2002 Code (Table 3) on these new findings, which have been used to de- rive the compulsory microclimatic condition for workplaces in the Czech Republic. 7 Acknowledgments I would like to thank Professor D. J. Nevrala for his help with the English text. References [1] ANSI/ASHRAE Standard 55-2004. Thermal Environ- mental Conditions for Human Occupancy. [2] Croome, D. J., Gan, G., Abwi, H. B.: Evaluation of In- door Environment in Naturally Ventilated Offices. In: Research on Indoor Air Quality and Climate. CIB Pro- ceedings, Publication 163, Rotterdam 1993. [3] EN ISO 7730 Moderate Thermal Environment. [4] European technical report CR 1752-1998 “Ventilation for Buildings: Design Criteria for the Indoor Environ- ment”. [5] Fanger, P. O.: Thermal Comfort. Danish Technical Press, Copenhagen 1970. [6] Fishman, D. S., Pimbert, S. L.: Survey of the Objective Responses to the Thermal Environment in Offices. In: Indoor Climate (eds P. O. Fanger, O. Valbjorn). Copen- hagen, Danish Building Research Institute, 1979, p. 677–698. [7] Humphreyes, M. A., Nicol, J. F.: Effects of Measurement and Formulation Error on Thermal Comfort Indices in the ASHRAE Database of Field Studies. ASHRAE Trans- actions, Vol. 106 (2000), No. 2, p. 493–502. [8] Humphreyes, M. A., Nicol, J. F.: Conflicting Criteria for Thermal Sensation within the Fanger Predicted Mean Vote Equation. In: CIBSE/ASHRAE Joint National Conference Proceedings, Harrogate, UK, Vol. 2, 1996, p. 153–158. [9] Humphreyes, M. A., Nicol, J. F.: An Analysis of Some Observations of Finger-Temperature and Thermal Comfort of Office Workers. In: Indoor Air, Edinburgh, UK, 1999. [10] Itoh, S., Ogata, K., Yoshimura, H.: Advances in Climatic Physiology. IGATU SHOIN LTD. Tokyo 1972, SPRINGER VERLAG Berlin, Heidelberg, New York 1972. [11] Jokl, M. V.: Microenvironment: The Theory and Prac- tice of Indoor Climate. Thomas, Illinois, USA, 1989. [12] Jokl, M. V., Moos, P., Štverák, J.: The Human Thermo- regulatory Range within the Neutral Zone. Physiol. Res. Vol. 41 (1992), p. 227–236. [13] Jokl, M. V., Moos, P.: Die Warmeregelungsgrenze des Menschen in neutraler Zone. Bauphysik Vol. 14 (1992), No. 6, p. 175–181. [14] Nařízení vlády č. 523/2002 Sb., kterým se mění naříze- ní vlády č. 178/2001 Sb., kterým se stanoví podmínky ochrany zaměstnanců při práci (Government Directive No. 523/2002 Code., changing Government Directive No. 178/2001 Code. prescribing the conditions for em- ployees protection during the work). [15] Newsham, G. R., Tiller, D. K.: A Field Study of Office Thermal Comfort Using Questionnaire Software. Na- tional Research Council Canada, Internal report No. 708, Nov. 1995. [16] Tanabe, S. I., Kimura, K. I., Hara, T., Akimoto, T.: Effects of Air Movement on Thermal Comfort in Air- -Conditioned Spaces During Summer Season. Journal of Architecture, Planning and Environmental Engineering, Vol. 382 (1987), p. 20–30. Prof. Ing. Miloslav Jokl, DrSc. phone: +420-22435-4432 e-mail: miloslav.jokl@fsv.cvut.cz Prof. Ing. Karel Kabele phone: +420-22435-4570 e-mail: kabele@fsv.cvut.cz Dept. of Microenvironmental and Building Services Engineering Czech Technical University in Prague Faculty of Civil Engineering Thákurova 7 166 29, Prague 6, Czech Republic © Czech Technical University Publishing House http://ctn.cvut.cz/ap/ 13 Acta Polytechnica Vol. 46 No. 6/2006