modelFanger

Fangers thermal comfort model

Extends from IDEAS.Buildings.Components.Comfort.BaseClasses.PartialComfort (Partial for comfort models).

Information

This model is based on the buildings library.

Thermal comfort model according to Fanger, as described in the ASHRAE Fundamentals (1997).

The thermal sensation of a human being is mainly related to the thermal balance of its body as a whole. This balance is influenced by two groups of factors, personal and physical. The activity level and clothing thermal insulation of the subject form the group of personal factors, while the environmental parameters: air temperature, mean radiant temperature, air velocity, and air humidity compose the group of physical factors. When the personal factors have been estimated and the physical factors have been measured the thermal sensation for the body as a whole (general thermal comfort) can be predicted by calculating the PMV index. The PPD index, obtained from the PMV index, provides information on thermal discomfort (thermal dissatisfaction) by predicting the percentage of people likely to feel too hot or too cold in the given thermal environment.

The Predicted Mean Vote (PMV) model combines four physical variables (air temperature, air velocity, mean radiant temperature, and relative humidity), and two personal variables (clothing insulation and activity level) into an index that can be used to predict the average thermal sensation of a large group of people.

To determine appropriate thermal conditions, practitioners refer to standards such as ASHRAE Standard 55 (ASHRAE, 1992) and ISO Standard 7730 (ISO, 1994). These standards define temperature ranges that should result in thermal satisfaction for at least 80% of occupants in a space.

PMV thermal sensation scale

The PMV index predicts the mean value of the votes of a large group of people on the following 7-point thermal sensation scale:

Cold Cool Slightly cool Neutral Slightly warm Warm Hot
-3 -2 -1 0 +1 +2 +3

Operative temperature

For a given space there exists an optimum operative temperature corresponding to PMV=0 (neutral). The operative temperature is defined as: The uniform temperature of an imaginary black enclosure in which an occupant would exchange the same amount of heat by radiation plus convection as in the actual nonuniform environment. The operative temperature is computed as the average of the air temperature and the mean radiant temperature, weighted by their respective heat transfer coefficients (see ASHRAE Fundamentals, 1997, page 8.3, eq (8)).

Optimum operative temperatures

Winter: activity 1.2 met,
clothing = 0.9 clo (sweater, long sleeve shirt, heavy pants),
air flow = 30 fpm (0.15 m/sec),
mean radiant temperature equal to air temperature,
Optimum Operative Temperature (top) = 22.7° C (71° F)

Summer: clothing = 0.5 clo,
air flow = 50 fpm (0.25 m/sec),
Optimum Operative Temperature (top) = 24.4° C (76° F).

All equation numbers in the model refer to the ASHRAE Handbook Fundamentals, Chapter 8, Thermal Comfort, 1997.

Usual ranges of variables (ISO)

M = 46 to 232 W/m^2 (0.8 to 4 met)
ICl = 0 to 2 clo (0 to 0.310 m^2*K/W)
TAir_degC = 10 to 30° C
TRad_degC = 10 to 40° C
vAir = 0 to 1 m/s
pSte = 0 to 2700 Pa

Insulation for clothing ensembles

Clothing is defined in terms of clo units. Clo is a unit used to express the thermal insulation provided by garments and clothing ensembles, where 1 clo = 0.155 (m^2*K/W) (ASHRAE 55-92).

The following table is obtained from ASHRAE page 8.8

Clothing ensembleclo
ASHRAE Standard 55 Winter0.90
ASHRAE Standard 55 Summer0.50
Walking shorts, short-sleeve shirt 0.36
Trousers, long-sleeve shirt 0.61
Trousers, long-sleeve shirt, suit jacket 0.96
Trousers, long-sleeve shirt, suit jacket, T-shirt 1.14
Trousers, long-sleeve shirt, long-sleeve sweater, T-shirt 1.01
Same as above + suit jacket, long underwear bottoms 1.30
Sweat pants, sweat shirt 0.74
Knee-length skirt, short-sleeve shirt, panty hose, sandals 0.54
Knee-length skirt, long-sleeve shirt, full slip, panty hose 0.67
Knee-length skirt, long-sleeve shirt, half slip, panty hose, long sleeve sweater 1.10
Long-sleeve coveralls, T-shirt 0.72
Insulated coveralls, long-sleeve, thermal underwear, long underwear bottoms 1.37

Metabolic rates

One met is defined as 58.2 Watts per square meter which is equal to the energy produced per unit surface area of a seated person at rest.

The following table is obtained from ASHRAE page 8.6.

ActivityW/m2 body surface area
ASHRAE Standard 5558.2
reclining 45
seated and quiet 60
sedentary activity (reading, writing) 60
standing, relaxed 70
office (filling while standing)80
office (walking)100
Sleeping 40
Seated quiet 60
Standing Relaxed 70
Walking 3.2 - 6.4km/h 115-220
Reading 55
Writing 60
Typing 65
Lifting/packing 120
Driving Car 60-115
Driving Heavy vehicle 185
Cooking 95-115
Housecleaning 115-200
Machine work 105-235
Pick and shovel work 235-280
Dancing-Social 140-225
Calisthenics 175-235
Basketball 290-440
Wrestling 410-505

References

  • ASHRAE Handbook, Fundamentals (SI Edition). American Society of Heating, Refrigerating and Air-Conditioning Engineers, Chapter 8, Thermal Comfort; pages 8.1-8.26; Atlanta, USA, 1997.
  • International Standards Organization (ISO). Moderate Thermal Environments: Determination of the PMV and PPD Indices and Specification of the Conditions for Thermal Comfort (ISO 7730). Geneva, Switzerland: ISO. 1994.
  • Charles, K.E. Fanger Thermal Comfort and Draught Models. Institute for Research in Construction National Research Council of Canada, Ottawa, K1A 0R6, Canada. IRC Research Report RR-162. October 2003. http://irc.nrc-cnrc.gc.ca/ircpubs.
  • Data, References and Links at: Thermal Comfort; Dr. Sam C M Hui Department of Mechanical Engineering The University of Hong Kong MEBS6006 Environmental Services I; http://me.hku.hk/msc-courses/MEBS6006/index.html

Parameters

TypeNameDefaultDescription
Booleanuse_phi_in (from PartialComfort)trueGet the relative humidity from the input connector
Realphi (from PartialComfort)0.4Default value of relative humidity
IDEAS.Buildings.Components.OccupancyType.OfficeWorkoccupancyType (from PartialComfort)
Modelica.Units.SI.HeatFluxW0Rate of mechanical work accomplished (must be non-positive, typically equal to 0)
Modelica.Units.SI.CoefficientOfHeatTransferhRad0.8*4.7Radiative heat transfer coefficient
Conditional inputs
Booleanuse_vAir_infalseGet the air velocity from the input connector
Booleanuse_M_infalseGet the metabolic rate from the input connector
Booleanuse_ICl_infalseGet the clothing insulation from the input connector
Booleanuse_pAir_infalseGet the air pressure from the input connector
Modelica.Units.SI.VelocityvAir0.05Fixed value for air velocity
Modelica.Units.SI.HeatFluxM(occupancyType.QsenPp + occupancyType.QlatPp)/1.9Fixed value for metabolic rate, assumes skin surface area of 1.9 m2
RealICloccupancyType.IClFixed value for clothing insulation in units of clo (summer=0.5; winter=0.9)
Modelica.Units.SI.PressurepAir101325Fixed value for air pressure

Connectors

TypeNameDefaultDescription
Modelica.Blocks.Interfaces.RealInputphi_in (from PartialComfort)Relative humidity
Modelica.Blocks.Interfaces.RealInputTRad (from PartialComfort)Radiation temperature
Modelica.Blocks.Interfaces.RealInputTAir (from PartialComfort)Air temperature
Modelica.Blocks.Interfaces.RealOutputPMVPMV
Modelica.Blocks.Interfaces.RealOutputPPDPPD [0.05...1]
Modelica.Blocks.Interfaces.RealInputpAir_inAir pressure
Modelica.Blocks.Interfaces.RealInputICl_inClothing thermal resistance in clo
Modelica.Blocks.Interfaces.RealInputvAir_inAir velocity
Modelica.Blocks.Interfaces.RealInputM_inMetabolic heat generation in W/m2 (not in met)

Components

TypeNameDefaultDescription
Modelica.Units.SI.TemperatureTOpeOperative temperature
Modelica.Units.SI.TemperatureTCloSurface temperature of clothing
Modelica.Units.SI.TemperatureTSkiSkin temperature
Modelica.Units.SI.CoefficientOfHeatTransferhComCombined heat transfer coefficient
Modelica.Units.SI.CoefficientOfHeatTransferhConConvective heat transfer coefficient
Modelica.Units.SI.PressurepStePartial pressure of water vapor in ambient air
Modelica.Units.SI.HeatFluxLThermal load of the body
RealfClClothing area factor (61)
Modelica.Units.SI.ThermalInsulanceRClThermal resistance of clothing (10)

Revisions

  • April 11, 2019 by Filip Jorissen:
    Revised implementation such that default value of relative humidity is used when using a dry air medium. See #1011.
  • Jul 18, 2016, by Filip Jorissen:
    Copied and changed implementation from Buildings library.
  • May 30, 2014, by Michael Wetter:
    Removed undesirable annotation Evaluate=true. Added default value for ICl.
  • October 9, 2013, by Michael Wetter:
    Corrected use of conditional connector.
  • July 14, 2010, by Michael Wetter:
    Merged models into Buildings library. Fixed bug in computation of lower value for hCon. The original implementation lead to too high a lower value.
  • August 2, 2005 Revised model, fixed bug in computing clothing surface temperature, changed various other computations, changed parameter and input to model, set clothing insulation value as input rather than computing it in model, added model to UTC library.
  • June, 2005 Michael Wetter and Sorin Costiner:
    Improved version, added PPDDraft, TOpe, performed studies
  • March 03, 2005 Michael Wetter and Sorin Costiner:
    First implementation.