Outside Surface Heat Balance[LINK]
Figure 18. Outside Heat Balance Control Volume Diagram
The heat balance on the outside face is:
q′′αsol+q′′LWR+q′′conv−q′′ko=0
where: _{q′′αsol } = Absorbed direct and diffuse solar (short wavelength) radiation heat flux. _{q′′LWR } = Net long wavelength (thermal) radiation flux exchange with the air and surroundings. _{q′′conv } = Convective flux exchange with outside air. _{q′′ko } = Conduction heat flux (q/A) into the wall.
All terms are positive for net flux to the face except the conduction term, which is traditionally taken to be positive in the direction from outside to inside of the wall. Simplified procedures generally combine the first three terms by using the concept of a solair temperature. Each of these heat balance components is introduced briefly below.
External Shortwave Radiation[LINK]
_{q′′αsol } is calculated using procedures presented later in this manual and includes both direct and diffuse incident solar radiation absorbed by the surface face. This is influenced by location, surface facing angle and tilt, surface face material properties, weather conditions, etc.
External Longwave Radiation[LINK]
q′′LWR is a standard radiation exchange formulation between the surface, the sky, and the ground. The radiation heat flux is calculated from the surface absorptivity, surface temperature, sky and ground temperatures, and sky and ground view factors.
The longwave radiation heat exchange between surfaces is dependent on surface temperatures, spatial relationships between surfaces and surroundings, and material properties of the surfaces. The relevant material properties of the surface, emissivity e and absorptivity a, are complex functions of temperature, angle, and wavelength for each participating surface. However, it is generally agreed that reasonable assumptions for building loads calculations are (Chapman 1984; Lienhard 1981):
· each surface emits or reflects diffusely and is gray and opaque (a = e, t = 0, r = 1 e)
· each surface is at a uniform temperature
· energy flux leaving a surface is evenly distributed across the surface,
· the medium within the enclosure is nonparticipating.
These assumptions are frequently used in all but the most critical engineering applications.
Table 3. Nomenclature List of Variables.
Mathematical variable

Description

Units

Range

q“_{LWR}

Exterior surface longwave radiation flux

W/m^{2}


h_{r}

Linearized radiative heat transfer coefficient to air temperature

W/(m^{2} K)


T_{surf}

Surface Outside face temperatures

K


T_{air}

Outside air temperature

K


T_{gnd}

Environmental ground surface temperature

K


T_{sky}

Sky Effective temperature

K


F_{gnd}

view factor of wall surface to ground surface


0~1

F_{sky}

View factor of wall surface to sky


0~1

F_{air}

View factor of wall surface to air


0~1

e

Surface longwave emissivity


0~1

s

StefanBoltzmann constant

W/m^{2}K^{4}

0.0000000567

Consider an enclosure consisting of building exterior surface, surrounding ground surface, and sky. Using the assumptions above, we can determine the longwave radiative heat flux at the building exterior surface (Walton 1983; McClellan and Pedersen 1997). The total longwave radiative heat flux is the sum of components due to radiation exchange with the ground, sky, and air.
q‘LWR‘=q′′gnd+q′′sky+q′′air
Applying the StefanBoltzmann Law to each component yields:
q‘LWR‘=εσFgnd(T4gnd−T4surf)+εσFsky(T4sky−T4surf)+εσFair(T4air−T4surf)
where
e = longwave emittance of the surface
s = StefanBoltzmann constant
F_{gnd} = view factor of wall surface to ground surface temperature
F_{sky} = view factor of wall surface to sky temperature
F_{air} = view factor of wall surface to air temperature
T_{surf} = outside surface temperature
T_{gnd} = ground surface temperature
T_{sky} = sky temperature
T_{air} = air temperature
Linearized radiative heat transfer coefficients are introduced to render the above equation more compatible with the heat balance formulation,
q‘LWR‘=hr,gnd(Tgnd−Tsurf)+hr,sky(Tsky−Tsurf)+hr,air(Tair−Tsurf)
where
hr,gnd=εσFgnd(T4surf−T4gnd)Tsurf−Tgnd
hr,sky=εσFsky(T4surf−T4sky)Tsurf−Tsky
hr,air=εσFair(T4surf−T4air)Tsurf−Tair
The longwave view factors to ground and sky are calculated with the following expressions (Walton 1983):
Fground=0.5(1−cosϕ)
Fsky=0.5(1+cosϕ)
where f is the tilt angle of the surface. The view factor to the sky is further split between sky and air radiation by:
β=√0.5(1+cosϕ)
The ground surface temperature is assumed to be the same as the air temperature. The final forms of the radiative heat transfer coefficients are shown here.
hr,gnd=εσFgnd(T4surf−T4air)Tsurf−Tair
hr,sky=εσFskyβ(T4surf−T4sky)Tsurf−Tsky
hr,air=εσFsky(1−β)(T4surf−T4air)Tsurf−Tair
ASHRAE. 1993. 1993 ASHRAE Handbook – Fundamentals. Atlanta: American Society of Heating, Refrigerating, and AirConditioning Engineers, Inc.
Chapman, A. J. 1984. Heat Transfer, 4^{th} Edition, New York: Macmillan Publishing Company.
Lienhard, J. H. 1981. A Heat Transfer Textbook, Englewood Cliffs, N.J.: PrenticeHall, Inc.
McClellan, T. M., and C. O. Pedersen. 1997. Investigation of Outside Heat Balance Models for Use in a Heat Balance Cooling Load Calculation. ASHRAE Transactions, Vol. 103, Part 2, pp. 469484.
Walton, G. N. 1983. Thermal Analysis Research Program Reference Manual. NBSSIR 832655. National Bureau of Standards.
Atmospheric Variation[LINK]
All buildings are located in the troposphere, the lowest layer of the atmosphere. The troposphere extends from sea level to an altitude of 11 km. Throughout the troposphere, air temperature decreases almost linearly with altitude at a rate of approximately 1°C per 150 m. Barometric pressure decreases more slowly. Wind speed, on the other hand, increases with altitude.
Because the atmosphere changes with altitude (defined as height above ground in this case), tall buildings can experience significant differences in local atmospheric properties between the ground floor and the top floor. Buildings interact with the atmosphere through convective heat transfer between the outdoor air and the exterior surfaces of the building envelope, and through the exchange of air between the outside and inside of the building via infiltration and ventilation.
Impetus for using this modeling is illustrated in the next table. Using a 70 story (284 meters) building as an example, the atmospheric variables are significant.
Table 4. Atmospheric Variables at Two Different Altitudes above Ground Level
Variable

1.5 Meters

284 meters

Absolute Diff

Percent Diff

Air Temperature

15°C

13.15°C

1.85°C

12.3%

Barometric Pressure

101,325 Pa

97,960 Pa

3,365 Pa

3.3%

Wind Speed

2.46 m/s

7.75 m/s

5.29 m/s

215%

Comparing the annual energy usage between 60 discretely modeled floors of a building, it turns out that the effect due to wind speed change is dominant over the first ten floors. But at floor 25, surprisingly, the effect due to air temperature has caught up and is about equal to the effect of wind speed. Above floor 25 the effect due to air temperature is now dominant. Clearly it is desirable to model air temperature variation with altitude for highrise buildings.
To accommodate atmospheric variation EnergyPlus automatically calculates the local outdoor air temperature and wind speed separately for each zone and surface that is exposed to the outdoor environment. The zone centroid or surface centroid are used to determine the height above ground. Only local outdoor air temperature and wind speed are currently calculated because they are important factors for the exterior convection calculation for surfaces (see Exterior Convection below) and can also be factors in the zone infiltration and ventilation calculations. Variation in barometric pressure, however, is considered when using the Airflow Network objects.
Local Outdoor Air Temperature Calculation[LINK]
Variation in outdoor air temperature is calculated using the U.S. Standard Atmosphere (1976). According to this model, the relationship between air temperature and altitude in a given layer of the atmosphere is:
Tz=Tb+L(Hz−Hb)
where
T_{z} = air temperature at altitude z
T_{b} = air temperature at the base of the layer, i.e., ground level for the troposphere
L = air temperature gradient, equal to –0.0065 K/m in the troposphere
H_{b} = offset equal to zero for the troposphere
H_{z} = geopotential altitude.
The variable H_{z} is defined by:
Hz=Ez(E+z)
where
E = 6,356 km, the radius of the Earth
z = altitude.
For the purpose of modeling buildings in the troposphere, altitude z refers to the height above ground level, not the height above sea level. The height above ground is calculated as the height of the centroid, or areaweighted center point, for each zone and surface.
The air temperature at ground level, T_{b}, is derived from the weather file air temperature by inverting the equation above:
Tb=Tz,met−L(EzmetE+zmet−Hb)
where
T_{z,met} = weather file air temperature (measured at the meteorological station)
z_{met} = height above ground of the air temperature sensor at the meteorological station.
The default value for z_{met} for air temperature measurement is 1.5 m above ground. This value can be overridden by using the Site:WeatherStation object.
Local Wind Speed Calculation[LINK]
Chapter 16 of the Handbook of Fundamentals (ASHRAE 2005). The wind speed measured at a meteorological station is extrapolated to other altitudes with the equation:
Vz=Vmet(δmetzmet)αmet(zδ)α
where
z = altitude, height above ground
V_{z} = wind speed at altitude z
a = wind speed profile exponent at the site
d = wind speed profile boundary layer thickness at the site
z_{met} = height above ground of the wind speed sensor at the meteorological station
V_{met} = wind speed measured at the meteorological station
a_{met} = wind speed profile exponent at the meteorological station
d_{met} = wind speed profile boundary layer thickness at the meteorological station.
The wind speed profile coefficients a, d, a_{met}, and d_{met}, are variables that depend on the roughness characteristics of the surrounding terrain. Typical values for a and d are shown in the following table:
Table 5. Wind Speed Profile Coefficients (ASHRAE Fundamentals 2005).
Terrain Description

Exponent,a

Boundary Layer Thickness,d (m)

Flat, open country

0.14

270

Rough, wooded country

0.22

370

Towns and cities

0.33

460

Ocean

0.10

210

Urban, industrial, forest

0.22

370

The terrain types above map to the options in the Terrain field of the Building object. The Terrain field can be overridden with specific values for a and d by using the Site:HeightVariation object.
The default value for z_{met} for wind speed measurement is 10 m above ground. The default values for a_{met} and d_{met} are 0.14 and 270 m, respectively, because most meteorological stations are located in an open field. These values can be overridden by using the Site:WeatherStation object.
Outdoor/Exterior Convection[LINK]
Heat transfer from surface convection is modeled using the classical formulation:
Qc=hc,extA(Tsurf−Tair)
where
Q_{c} = rate of exterior convective heat transfer
h_{c,ext} = exterior convection coefficient
A = surface area
T_{surf} = surface temperature
T_{air} = outdoor air temperature
Substantial research has gone into the formulation of models for estimating the exterior convection coefficient. Since the 1930’s there have been many different methods published for calculating this coefficient, with much disparity between them (Cole and Sturrock 1977; Yazdanian and Klems 1994). More recently Palyvos (2008) surveyed correlations cataloging some 91 different correlations into four categories based on functional form of the model equation. EnergyPlus therefore offers a wide selection of different methods for determining values for h_{c,ext}. The selection of model equations for h_{c,ext} can be made at two different levels. The first is the set of options available in the input object SurfaceConvectionAlgorithm:Outside that provides a way of broadly selecting which model equations are applied throughout the model. The input objects SurfaceProperty:ConvectionCoefficients and SurfaceProperty:ConvectionCoefficients:MultipleSurface also provide ways of selecting which model equations or values are applied for specific surfaces. These basic options are identified by the key used for input and include:
Note that when the outside environment indicates that it is raining, the exterior surfaces (exposed to wind) are assumed to be wet. The convection coefficient is set to a very high number (1000) and the outside temperature used for the surface will be the wetbulb temperature. (If you choose to report this variable, you will see 1000 as its value.)
When the AdaptiveConvectionAlgorithm is used, there is a second, deeper level of control available for selecting among a larger variety of h_{c,ext} equations and also defining custom equations using curve or table objects. These options are described in this section.
In addition to the correlation choices described below, it is also possible to override the convection coefficients on the outside of any surface by other means:
* Use the SurfaceProperty:ConvectionCoefficients object in the input file to set the convection coefficient value on either side of any surface.
* Use the SurfaceProperty:OtherSideCoefficients object in the input file to set heat transfer coefficients and temperatures on surfaces.
* Use the EnergyManagementSystem Actuators that are available for overriding h_{c} values.
These options can also use schedules to control values over time. Specific details are given in the Input Output Reference document.
Simple Combined[LINK]
The simple algorithm uses surface roughness and local surface windspeed to calculate the exterior heat transfer coefficient (key:SimpleCombined). The basic equation used is:
h=D+EVz+FVz2
where
h = heat transfer coefficient
V_{z} = local wind speed calculated at the height above ground of the surface centroid
D, E, F = material roughness coefficients
The roughness correlation is taken from Figure 1, Page 22.4, ASHRAE Handbook of Fundamentals (ASHRAE 1989). The roughness coefficients are shown in the following table:
Table 6. Roughness Coefficients D, E, and F.
Roughness Index

D

E

F

Example Material

1 (Very Rough)

11.58

5.894

0.0

Stucco

2 (Rough)

12.49

4.065

0.028

Brick

3 (Medium Rough)

10.79

4.192

0.0

Concrete

4 (Medium Smooth)

8.23

4.0

0.057

Clear pine

5 (Smooth)

10.22

3.1

0.0

Smooth Plaster

6 (Very Smooth)

8.23

3.33

0.036

Glass

Note that the simple correlation yields a combined convection and radiation heat transfer coefficient. Radiation to sky, ground, and air is included in the exterior convection coefficient for this algorithm.
All other algorithms yield a convection only heat transfer coefficient. Radiation to sky, ground, and air is calculated automatically by the program.
TARP ALGORITHM[LINK]
TARP, or Thermal Analysis Research Program, is an important predecessor of EnergyPlus (Walton 1983). Walton developed a comprehensive model for exterior convection by blending correlations from ASHRAE and flat plate experiments by Sparrow et. al. In older versions of EnergyPlus, prior to version 6, the “TARP” model was called “Detailed.” The model was reimplemented in version 6 to use Area and Perimeter values for the group of surfaces that make up a facade or roof, rather than the single surface being modeled.
Table 7. Nomenclature List of Variables.
Variable

Description

Units

Range

A

Surface area of the surface

m^{2}

/ = 0

h_{c}

Surface exterior convective heat transfer coefficient

W/(m^{2}K)


h_{f}

Forced convective heat transfer coefficient

W/(m^{2}K)


h_{n}

Natural convective heat transfer coefficient

W/(m^{2}K)


P

Perimeter of surface

m


R_{f}

Surface roughness multiplier



T_{air}

Local outdoor air temperature calculated at the height above ground of the surface centroid

°C


T_{so}

Outside surface temperature

°C


DT

Temperature difference between the surface and air,

°C


V_{z}

Local wind speed calculated at the height above ground of the surface centroid

m/s


W_{f}

Wind direction modifier



f

Angle between the ground outward normal and the surface outward normal

degree


Roughness Index

Surface roughness index (6 = very smooth, 5 = smooth, 4 = medium smooth, 3 = medium rough, 2 = rough, 1 = very rough)


1~6

The Detailed, BLAST, and TARP convection models are very similar. In all three models, convection is split into forced and natural components (Walton 1981). The total convection coefficient is the sum of these components.
hc=hf+hn
The forced convection component is based on a correlation by Sparrow, Ramsey, and Mass (1979):
hf=2.537WfRf(PVzA)1/2
where
W_{f} = 1.0 for windward surfaces
or
W_{f} = 0.5 for leeward surfaces
Leeward is defined as greater than 100 degrees from normal incidence (Walton 1981).
The surface roughness multiplier Rf is based on the ASHRAE graph of surface conductance (ASHRAE 1981) and may be obtained from the following table:
Table 8. Surface Roughness Multipliers (Walton 1981).
Roughness Index

Rf

Example Material

1 (Very Rough)

2.17

Stucco

2 (Rough)

1.67

Brick

3 (Medium Rough)

1.52

Concrete

4 (Medium Smooth)

1.13

Clear pine

5 (Smooth)

1.11

Smooth Plaster

6 (Very Smooth)

1.00

Glass

The natural convection component h_{n} is calculated in the same way as the interior “Detailed” model. The detailed natural convection model correlates the convective heat transfer coefficient to the surface orientation and the difference between the surface and zone air temperatures (where DT = Air Temperature  Surface Temperature). The algorithm is taken directly from Walton (1983). Walton derived his algorithm from the ASHRAE Handbook (2001), Table 5 on p. 3.12, which gives equations for natural convection heat transfer coefficients in the turbulent range for large, vertical plates and for large, horizontal plates facing upward when heated (or downward when cooled). A note in the text also gives an approximation for large, horizontal plates facing downward when heated (or upward when cooled) recommending that it should be half of the facing upward value. Walton adds a curve fit as a function of the cosine of the tilt angle to provide intermediate values between vertical and horizontal. The curve fit values at the extremes match the ASHRAE values very well.
For no temperature difference OR a vertical surface the following correlation is used:
h=1.31ΔT13
For (DT < 0.0 AND an upward facing surface) OR (DT > 0.0 AND an downward facing surface) an enhanced convection correlation is used:
h=9.482ΔT137.283−cosΣ
where S is the surface tilt angle.
For (DT > 0.0 AND an upward facing surface) OR (DT < 0.0 AND an downward facing surface) a reduced convection correlation is used:
h=1.810ΔT131.382+cosΣ
where S is the surface tilt angle.
MoWiTT Algorithm[LINK]
Table 9. Nomenclature List of Variables.
Variable

Description

Units

Range

A

Constant

W/(m^{2}K(m/s)^{b}


B

Constant



C_{t}

Turbulent natural convection constant

W/(m^{2}K^{4/3})


h_{c}

Surface exterior convective heat transfer coefficient

W/(m^{2}K)


T_{so}

Outside surface temperature

°C/K


DT

Temperature difference between the surface and air

°C/K


The MoWiTT model is based on measurements taken at the Mobile Window Thermal Test (MoWiTT) facility (Yazdanian and Klems 1994). The correlation applies to very smooth, vertical surfaces (e.g. window glass) in lowrise buildings and has the form:
hc=√[Ct(ΔT)13]2+[aVbz]2
Constants a, b and turbulent natural convection constant C_{t} are given in Table 10. The original MoWiTT model has been modified for use in EnergyPlus so that it is sensitive to the local suface’s wind speed which varies with the height above ground. The original MoWiTT model was formulated for use with the air velocity at the location of the weather station. As of Version 7.2, EnergyPlus uses the “a” model coefficients derived by Booten et al. (2012) rather than the original values from Yazdanian and Klems (1994).
NOTE: The MoWiTT algorithm may not be appropriate for rough surfaces, highrise surfaces, or surfaces that employ movable insulation.
Table 10. MoWiTT Coefficients (Yazdanian and Klems 1994, Booten et al. 2012)
Wind Direction

Ct

a

b

(Units)

W/m^{2}K^{4/3}

W/m^{2}K(m/s)^{b}


Windward

0.84

3.26

0.89

Leeward

0.84

3.55

0.617

Table 11. Nomenclature List of Variables.
Variable

Description

Units

Range

a

Constant

W/(m^{2}K(m/s)^{b}


b

Constant



h_{c}

Surface exterior convective heat transfer coefficient

W/(m^{2}K)


h_{c,glass}

Convective heat transfer coefficient for very smooth surfaces (glass)

W/(m^{2}K)


h_{n}

Natural convective heat transfer coefficient

W/(m^{2}K)


R_{f}

Surface roughness multiplier



T_{so}

Outside surface temperature

°C/K


DT

Temperature difference between the surface and air,

°C/K


f

Angle between the ground outward normal and the surface outward normal

radian


The DOE2 convection model is a combination of the MoWiTT and BLAST Detailed convection models (LBL 1994). The convection coefficient for very smooth surfaces (e.g. glass) is calculated as:
hc,glass=√h2n+[aVbz]2
h_{n} is calculated using Equation or Equation . Constants a and b are given in Table 10.
For less smooth surfaces, the convection coefficient is modified according to the equation
hc=hn+Rf(hc,glass−hn)
where R_{f} is the roughness multiplier given by Table 8.
Adaptive Convection Algorithm[LINK]
This algorithm has a structure that allows for finer control over the models used for particular surfaces. The algorithm for the outside face was developed for EnergyPlus but it borrows concepts and its name from the research done by BeausoleilMorrison (2000, 2002) for convection at the inside face (see the description below for interior convection).
The adaptive convection algorithm implemented in EnergyPlus for the outside face is much simpler than that for the inside face. The surface classification system has a total of 4 different categories for surfaces that depend on current wind direction and heat flow directions. However it is more complex in that the h_{c} equation is split into two parts and there are separate model equation selections for forced convection, h_{f}, and natural convection, h_{n}. The following table summarizes the categories and the default assignments for h_{c} equations. The individual h_{c} equations are documented below.
Table 12. Adaptive Convection Algorithm Details
#

Surface Classification

Heat Flow Direction

Wind Direction

h_{f} Models

h_{n} Models

1

Roof Stable

Down

Any

TARPWindward MoWiTTWindward DOE2Windward NusseltJurges BlockenWindward EmmelRoof ClearRoof

WaltonStableHorizontalOrTilt AlamdariStableHorizontal

2

Roof Unstable

Up

Any

TARPWindward MoWiTTWindward DOE2Windward NusseltJurges BlockenWindward EmmelRoof ClearRoof

WaltonUnstableHorizontalOrTilt AlamdariUnstableHorizontal

3

Vertical Wall Windward

Any

Windward

TARPWindward DOE2Windward MoWiTTWindward NusseltJurges McAdams Mitchell BlockenWindward EmmelVertical

ASHRAEVerticalWall AlamdariHammondVerticalWall FohannoPolidoriVerticalWall ISO15099Windows

4

Vertical Wall Leeward

Any

Leeward

TARPLeeward MoWiTTLeeward DOE2Leeward EmmelVertical NusseltJurges McAdams Mitchell

ASHRAEVerticalWall AlamdariHammondVerticalWall FohannoPolidoriVerticalWall ISO15099Windows

Outside Face Surface Classification[LINK]
During an initial setup phase, all the heat transfer surfaces in the input file are analyzed in groups to determine appropriate values for geometry scales used in many of the convection correlations. Eight separate groups are assembled for nominally vertical exterior surfaces for eight bins of azimuth: north, northeast, east, southeast, south, southwest, west, northwest. Surfaces with the same range of azimuth are grouped together and analyzed for overall geometry parameters. A ninth group is assembled for nominally horizontal exterior surfaces for a roof bin that is also analyzed for geometry. These geometry routines find bounds and limits of all the surfaces in the group and then model geometric parameters from these limits.
Sparrow Windward[LINK]
As discussed above for the TARP algorithm, a Sparrow et al. (1979) conducted flat plate measurements and develop the following correlation for finitesize flat plates oriented to windward.
hf=2.53Rf(PVzA)1/122
Sparrow Leeward[LINK]
Sparrow et al. (1979) conducted flat plate measurements and develop the following correlation for finitesize flat plates oriented to leeward.
hf=2.532Rf(PVzA)1/122
MoWITT Windward[LINK]
As discussed above, Yazdanian and Klems (1994) used outdoor laboratory measurements to develop the following correlation for smooth surfaces oriented to windward. Booten et al. (2012) developed revised coefficients for use with local surface wind speeds.
hc=√[0.84ΔT1/133]2+[2.38V0.89z]2
This model equation is for the total film coefficient and includes the natural convection portion. Therefore it should not be used in conjunction with a second natural convection model equation.
MoWITT Leeward[LINK]
Yazdanian and Klems (1994) used outdoor laboratory measurements to develop the following correlation for smooth surfaces oriented to leeward. Booten et al. (2012) developed revised coefficients for use with local surface wind speeds.
hc=√[0.84ΔT1/133]2+[2.86V0.617z]2
This model equation is for the total film coefficient and includes the natural convection portion. Therefore it should not be used in conjunction with a second natural convection model equation.
Blocken et al. (2009) developed a set of correlations for windward facing outdoor surfaces using numerical methods (key: BlockenWindward).
hf=4.6V0.8910m:θ≤11.25hf=5.0V0.8010m:11.25<θ≤33.75hf=4.6V0.8410m:33.75<θ≤56.25hf=4.5V0.8110m:56.25<θ≤100.0
Where V_{10m} is the air velocity at the location of the weather station and θ is the angle of incidence between the wind and the surface in degrees. This model is only applicable to windward surfaces and lacks a natural convection component and therefore cannot be used on its own but only within the adaptive convection algorithm for the outside face.
Clear et al. (2003) developed correlations from measurements for horizontal roofs on two commercial buildings. In EnergyPlus the implementation uses the model for natural convection plus turbulent forced convection (eq. 8A in the reference) and applies it to the center point of each surface section that makes up the roof.
hc=ηkLn0.15Ra1/133Ln+kxRf0.0296Re4/455x1/133Pr
Where
x is the distance to the surface centroid from where the wind begins to intersect the roof. In EnergyPlus this is currently simplified to half the square root of the roof surface.
Ln=AreaPerimeter of overall roof
k is the thermal conductivity of air
$= is the weighting factor for natural convection (suppressed at high forced convection rates)
R{a_{{L_n}}} = G{r_{{L_n}}}$ is the Rayleigh number
GrLn=gρ2Ln3ΔTTfμ2 is the Grashof number
Rex=Vzρxμ is the Reynolds number at x
Pr is the Prandtl number
This model only claims to be applicable to horizontal roof surfaces so it may not be applicable to tilted roofs. It combines natural and forced convection and therefore should not be used in conjunction with yet another natural convection model.
Emmel et al. (2007) developed a set of correlations for outdoor surfaces using numerical methods. The following equations are for vertical surfaces (key: EmmelVertical):
hf=5.15V0.8110m:θ≤22.5hf=3.34V0.8410m:22.5<θ≤67.5hf=4.78V0.7110m:67.5<θ≤112.5hf=4.05V0.7710m:112.5<θ≤157.5hf=3.54V0.7610m:157.5<θ≤180.0
Where V_{10m} is the air velocity at the location of the weather station and θ is the angle of incidence between the wind and the surface in degrees. The following equations are used for horizontal (roof) surfaces (key: EmmelRoof):
hf=5.11V0.7810m:θ≤22.5hf=4.60V0.7910m:22.5<θ≤67.5hf=3.67V0.8510m:67.5<θ≤90
Where θ is the angle of incidence between the wind and the longest edge of the roof surface in degrees.
This model is for all wind directions but lacks a natural convection component. The model was developed for simple, rectangular lowrise buildings. It is available only within the adaptive convection algorithm for the outside face
Nusselt Jurges[LINK]
Perhaps the oldest equation for winddriven convection was developed by Nusselt and Jurges (1922). Palyvos (2008) casts their model in simplified form in SI units as:
hc=5.8+3.94Vz
Where V_{z} is the wind velocity in m/s, in EnergyPlus that velocity is adjusted for height above ground using the z axis coordinate of the surface’s centroid and the site wind model. This model can be applied to all surfaces and the relatively large constant is assumed to represent the natural convection portion of a total convection coefficient. The model is not sensitive to wind direction nor surface roughness.
A venerable equation for winddriven convection was developed by McAdams (1954) which Palyvos (2008) casts in SI units as:
hc=5.7+3.8Vz
Where V_{z} is the wind velocity in m/s that has been adjusted for height above ground using the z axis coordinate of the surface’s centroid. This model can be applied to all surfaces and the relatively large constant is assumed to represent the natural convection portion of a total convection coefficient. The model is not sensitive to wind direction nor surface roughness.
A useful geometric scale based on building volume is used in an equation developed by Mitchell (1976). The winddriven convection equation is cast by Palyvos as:
hf=8.6V0.6zL0.4
Where V_{z} is the wind velocity in m/s that has been adjusted for height above ground using the z axis coordinate of the surface’s centroid and L is the cube root of the building’s total volume. EnergyPlus interprets this as the sum of the volume of all the zones in the input file.
Exterior/External Conduction[LINK]
The conduction term, q′′ko _{,} can in theory be calculated using a wide variety of heat conduction formulations. Typically in EnergyPlus, the Conduction Transfer Function (CTF) method is used. The available models are described in this section: Conduction Through The Walls.
ASHRAE. 1981. 1981 ASHRAE Handbook – Fundamentals, Atlanta: American Society of Heating, Refrigerating, and AirConditioning Engineers, Inc.
ASHRAE. 1989. 1989 ASHRAE Handbook – Fundamentals, Atlanta: American Society of Heating, Refrigerating, and AirConditioning Engineers, Inc.
ASHRAE. 1993. 1993 ASHRAE Handbook – Fundamentals, Chapter 3, Heat Transfer, IP & SI Editions, Atlanta: American Society of Heating, Refrigerating, and AirConditioning Engineers, Inc.
ASHRAE. 2001. 2001 ASHRAE Handbook – Fundamentals, Atlanta: American Society of Heating, Refrigerating, and AirConditioning Engineers, Inc.
ASHRAE. 2005. 2005 ASHRAE Handbook – Fundamentals, Chapter 16, Air Flow Around Buildings, Atlanta: American Society of Heating, Refrigerating, and AirConditioning Engineers, Inc.
Booten, C., N. Kruis, and C. Christensen. 2012. Identifying and Resolving Issues in EnergyPlus and DOE2 Window Heat Transfer Calculations. National Renewable Energy Laboratory. NREL/TP550055787. Golden, CO.
Cole, R. J., and N. S. Sturrock. 1977. The Convective Heat Exchange at the External Surface of Buildings. Building and Environment, Vol. 12, p. 207.
Ellis, P.G., and P.A. Torcellini. 2005. “Simulating Tall Buildings Using EnergyPlus”, Proceedings of the Ninth International IBPSA Conference, Building Simulation 2005, Montreal, Canada, August 1518, 2005.
Lawrence Berkeley Laboratory (LBL). 1994. DOE2.1E053 source code.
Sparrow, E. M., J. W. Ramsey, and E. A. Mass. 1979. Effect of Finite Width on Heat Transfer and Fluid Flow about an Inclined Rectangular Plate. Journal of Heat Transfer, Vol. 101, p. 204.
U.S. Standard Atmosphere. 1976. U.S. Government Printing Office, Washington, D.C.
Walton, G. N. 1981. Passive Solar Extension of the Building Loads Analysis and System Thermodynamics (BLAST) Program, Technical Report, United States Army Construction Engineering Research Laboratory, Champaign, IL.
Walton, G. N. 1983. Thermal Analysis Research Program Reference Manual. NBSSIR 832655. National Bureau of Standards.
Yazdanian, M. and J. H. Klems. 1994. Measurement of the Exterior Convective Film Coefficient for Windows in LowRise Buildings. ASHRAE Transactions, Vol. 100, Part 1, p. 1087.
Outside Surface Heat Balance[LINK]
Figure 18. Outside Heat Balance Control Volume Diagram
The heat balance on the outside face is:
q′′αsol+q′′LWR+q′′conv−q′′ko=0
where: _{q′′αsol } = Absorbed direct and diffuse solar (short wavelength) radiation heat flux. _{q′′LWR } = Net long wavelength (thermal) radiation flux exchange with the air and surroundings. _{q′′conv } = Convective flux exchange with outside air. _{q′′ko } = Conduction heat flux (q/A) into the wall.
All terms are positive for net flux to the face except the conduction term, which is traditionally taken to be positive in the direction from outside to inside of the wall. Simplified procedures generally combine the first three terms by using the concept of a solair temperature. Each of these heat balance components is introduced briefly below.
External Shortwave Radiation[LINK]
_{q′′αsol } is calculated using procedures presented later in this manual and includes both direct and diffuse incident solar radiation absorbed by the surface face. This is influenced by location, surface facing angle and tilt, surface face material properties, weather conditions, etc.
External Longwave Radiation[LINK]
q′′LWR is a standard radiation exchange formulation between the surface, the sky, and the ground. The radiation heat flux is calculated from the surface absorptivity, surface temperature, sky and ground temperatures, and sky and ground view factors.
The longwave radiation heat exchange between surfaces is dependent on surface temperatures, spatial relationships between surfaces and surroundings, and material properties of the surfaces. The relevant material properties of the surface, emissivity e and absorptivity a, are complex functions of temperature, angle, and wavelength for each participating surface. However, it is generally agreed that reasonable assumptions for building loads calculations are (Chapman 1984; Lienhard 1981):
· each surface emits or reflects diffusely and is gray and opaque (a = e, t = 0, r = 1 e)
· each surface is at a uniform temperature
· energy flux leaving a surface is evenly distributed across the surface,
· the medium within the enclosure is nonparticipating.
These assumptions are frequently used in all but the most critical engineering applications.
Table 3. Nomenclature List of Variables.Consider an enclosure consisting of building exterior surface, surrounding ground surface, and sky. Using the assumptions above, we can determine the longwave radiative heat flux at the building exterior surface (Walton 1983; McClellan and Pedersen 1997). The total longwave radiative heat flux is the sum of components due to radiation exchange with the ground, sky, and air.
q‘LWR‘=q′′gnd+q′′sky+q′′air
Applying the StefanBoltzmann Law to each component yields:
q‘LWR‘=εσFgnd(T4gnd−T4surf)+εσFsky(T4sky−T4surf)+εσFair(T4air−T4surf)
where
e = longwave emittance of the surface
s = StefanBoltzmann constant
F_{gnd} = view factor of wall surface to ground surface temperature
F_{sky} = view factor of wall surface to sky temperature
F_{air} = view factor of wall surface to air temperature
T_{surf} = outside surface temperature
T_{gnd} = ground surface temperature
T_{sky} = sky temperature
T_{air} = air temperature
Linearized radiative heat transfer coefficients are introduced to render the above equation more compatible with the heat balance formulation,
q‘LWR‘=hr,gnd(Tgnd−Tsurf)+hr,sky(Tsky−Tsurf)+hr,air(Tair−Tsurf)
where
hr,gnd=εσFgnd(T4surf−T4gnd)Tsurf−Tgnd
hr,sky=εσFsky(T4surf−T4sky)Tsurf−Tsky
hr,air=εσFair(T4surf−T4air)Tsurf−Tair
The longwave view factors to ground and sky are calculated with the following expressions (Walton 1983):
Fground=0.5(1−cosϕ)
Fsky=0.5(1+cosϕ)
where f is the tilt angle of the surface. The view factor to the sky is further split between sky and air radiation by:
β=√0.5(1+cosϕ)
The ground surface temperature is assumed to be the same as the air temperature. The final forms of the radiative heat transfer coefficients are shown here.
hr,gnd=εσFgnd(T4surf−T4air)Tsurf−Tair
hr,sky=εσFskyβ(T4surf−T4sky)Tsurf−Tsky
hr,air=εσFsky(1−β)(T4surf−T4air)Tsurf−Tair
References[LINK]
ASHRAE. 1993. 1993 ASHRAE Handbook – Fundamentals. Atlanta: American Society of Heating, Refrigerating, and AirConditioning Engineers, Inc.
Chapman, A. J. 1984. Heat Transfer, 4^{th} Edition, New York: Macmillan Publishing Company.
Lienhard, J. H. 1981. A Heat Transfer Textbook, Englewood Cliffs, N.J.: PrenticeHall, Inc.
McClellan, T. M., and C. O. Pedersen. 1997. Investigation of Outside Heat Balance Models for Use in a Heat Balance Cooling Load Calculation. ASHRAE Transactions, Vol. 103, Part 2, pp. 469484.
Walton, G. N. 1983. Thermal Analysis Research Program Reference Manual. NBSSIR 832655. National Bureau of Standards.
Atmospheric Variation[LINK]
All buildings are located in the troposphere, the lowest layer of the atmosphere. The troposphere extends from sea level to an altitude of 11 km. Throughout the troposphere, air temperature decreases almost linearly with altitude at a rate of approximately 1°C per 150 m. Barometric pressure decreases more slowly. Wind speed, on the other hand, increases with altitude.
Because the atmosphere changes with altitude (defined as height above ground in this case), tall buildings can experience significant differences in local atmospheric properties between the ground floor and the top floor. Buildings interact with the atmosphere through convective heat transfer between the outdoor air and the exterior surfaces of the building envelope, and through the exchange of air between the outside and inside of the building via infiltration and ventilation.
Impetus for using this modeling is illustrated in the next table. Using a 70 story (284 meters) building as an example, the atmospheric variables are significant.
Table 4. Atmospheric Variables at Two Different Altitudes above Ground LevelComparing the annual energy usage between 60 discretely modeled floors of a building, it turns out that the effect due to wind speed change is dominant over the first ten floors. But at floor 25, surprisingly, the effect due to air temperature has caught up and is about equal to the effect of wind speed. Above floor 25 the effect due to air temperature is now dominant. Clearly it is desirable to model air temperature variation with altitude for highrise buildings.
To accommodate atmospheric variation EnergyPlus automatically calculates the local outdoor air temperature and wind speed separately for each zone and surface that is exposed to the outdoor environment. The zone centroid or surface centroid are used to determine the height above ground. Only local outdoor air temperature and wind speed are currently calculated because they are important factors for the exterior convection calculation for surfaces (see Exterior Convection below) and can also be factors in the zone infiltration and ventilation calculations. Variation in barometric pressure, however, is considered when using the Airflow Network objects.
Local Outdoor Air Temperature Calculation[LINK]
Variation in outdoor air temperature is calculated using the U.S. Standard Atmosphere (1976). According to this model, the relationship between air temperature and altitude in a given layer of the atmosphere is:
Tz=Tb+L(Hz−Hb)
where
T_{z} = air temperature at altitude z
T_{b} = air temperature at the base of the layer, i.e., ground level for the troposphere
L = air temperature gradient, equal to –0.0065 K/m in the troposphere
H_{b} = offset equal to zero for the troposphere
H_{z} = geopotential altitude.
The variable H_{z} is defined by:
Hz=Ez(E+z)
where
E = 6,356 km, the radius of the Earth
z = altitude.
For the purpose of modeling buildings in the troposphere, altitude z refers to the height above ground level, not the height above sea level. The height above ground is calculated as the height of the centroid, or areaweighted center point, for each zone and surface.
The air temperature at ground level, T_{b}, is derived from the weather file air temperature by inverting the equation above:
Tb=Tz,met−L(EzmetE+zmet−Hb)
where
T_{z,met} = weather file air temperature (measured at the meteorological station)
z_{met} = height above ground of the air temperature sensor at the meteorological station.
The default value for z_{met} for air temperature measurement is 1.5 m above ground. This value can be overridden by using the Site:WeatherStation object.
Local Wind Speed Calculation[LINK]
Chapter 16 of the Handbook of Fundamentals (ASHRAE 2005). The wind speed measured at a meteorological station is extrapolated to other altitudes with the equation:
Vz=Vmet(δmetzmet)αmet(zδ)α
where
z = altitude, height above ground
V_{z} = wind speed at altitude z
a = wind speed profile exponent at the site
d = wind speed profile boundary layer thickness at the site
z_{met} = height above ground of the wind speed sensor at the meteorological station
V_{met} = wind speed measured at the meteorological station
a_{met} = wind speed profile exponent at the meteorological station
d_{met} = wind speed profile boundary layer thickness at the meteorological station.
The wind speed profile coefficients a, d, a_{met}, and d_{met}, are variables that depend on the roughness characteristics of the surrounding terrain. Typical values for a and d are shown in the following table:
Table 5. Wind Speed Profile Coefficients (ASHRAE Fundamentals 2005).The terrain types above map to the options in the Terrain field of the Building object. The Terrain field can be overridden with specific values for a and d by using the Site:HeightVariation object.
The default value for z_{met} for wind speed measurement is 10 m above ground. The default values for a_{met} and d_{met} are 0.14 and 270 m, respectively, because most meteorological stations are located in an open field. These values can be overridden by using the Site:WeatherStation object.
Outdoor/Exterior Convection[LINK]
Heat transfer from surface convection is modeled using the classical formulation:
Qc=hc,extA(Tsurf−Tair)
where
Q_{c} = rate of exterior convective heat transfer
h_{c,ext} = exterior convection coefficient
A = surface area
T_{surf} = surface temperature
T_{air} = outdoor air temperature
Substantial research has gone into the formulation of models for estimating the exterior convection coefficient. Since the 1930’s there have been many different methods published for calculating this coefficient, with much disparity between them (Cole and Sturrock 1977; Yazdanian and Klems 1994). More recently Palyvos (2008) surveyed correlations cataloging some 91 different correlations into four categories based on functional form of the model equation. EnergyPlus therefore offers a wide selection of different methods for determining values for h_{c,ext}. The selection of model equations for h_{c,ext} can be made at two different levels. The first is the set of options available in the input object SurfaceConvectionAlgorithm:Outside that provides a way of broadly selecting which model equations are applied throughout the model. The input objects SurfaceProperty:ConvectionCoefficients and SurfaceProperty:ConvectionCoefficients:MultipleSurface also provide ways of selecting which model equations or values are applied for specific surfaces. These basic options are identified by the key used for input and include:
SimpleCombined
TARP
MoWiTT
DOE2
AdaptiveConvectionAlgorithm
Note that when the outside environment indicates that it is raining, the exterior surfaces (exposed to wind) are assumed to be wet. The convection coefficient is set to a very high number (1000) and the outside temperature used for the surface will be the wetbulb temperature. (If you choose to report this variable, you will see 1000 as its value.)
When the AdaptiveConvectionAlgorithm is used, there is a second, deeper level of control available for selecting among a larger variety of h_{c,ext} equations and also defining custom equations using curve or table objects. These options are described in this section.
In addition to the correlation choices described below, it is also possible to override the convection coefficients on the outside of any surface by other means:
* Use the SurfaceProperty:ConvectionCoefficients object in the input file to set the convection coefficient value on either side of any surface.
* Use the SurfaceProperty:OtherSideCoefficients object in the input file to set heat transfer coefficients and temperatures on surfaces.
* Use the EnergyManagementSystem Actuators that are available for overriding h_{c} values.
These options can also use schedules to control values over time. Specific details are given in the Input Output Reference document.
Simple Combined[LINK]
The simple algorithm uses surface roughness and local surface windspeed to calculate the exterior heat transfer coefficient (key:SimpleCombined). The basic equation used is:
h=D+EVz+FVz2
where
h = heat transfer coefficient
V_{z} = local wind speed calculated at the height above ground of the surface centroid
D, E, F = material roughness coefficients
The roughness correlation is taken from Figure 1, Page 22.4, ASHRAE Handbook of Fundamentals (ASHRAE 1989). The roughness coefficients are shown in the following table:
Table 6. Roughness Coefficients D, E, and F.Note that the simple correlation yields a combined convection and radiation heat transfer coefficient. Radiation to sky, ground, and air is included in the exterior convection coefficient for this algorithm.
All other algorithms yield a convection only heat transfer coefficient. Radiation to sky, ground, and air is calculated automatically by the program.
TARP ALGORITHM[LINK]
TARP, or Thermal Analysis Research Program, is an important predecessor of EnergyPlus (Walton 1983). Walton developed a comprehensive model for exterior convection by blending correlations from ASHRAE and flat plate experiments by Sparrow et. al. In older versions of EnergyPlus, prior to version 6, the “TARP” model was called “Detailed.” The model was reimplemented in version 6 to use Area and Perimeter values for the group of surfaces that make up a facade or roof, rather than the single surface being modeled.
Table 7. Nomenclature List of Variables.The Detailed, BLAST, and TARP convection models are very similar. In all three models, convection is split into forced and natural components (Walton 1981). The total convection coefficient is the sum of these components.
hc=hf+hn
The forced convection component is based on a correlation by Sparrow, Ramsey, and Mass (1979):
hf=2.537WfRf(PVzA)1/2
where
W_{f} = 1.0 for windward surfaces
or
W_{f} = 0.5 for leeward surfaces
Leeward is defined as greater than 100 degrees from normal incidence (Walton 1981).
The surface roughness multiplier Rf is based on the ASHRAE graph of surface conductance (ASHRAE 1981) and may be obtained from the following table:
Table 8. Surface Roughness Multipliers (Walton 1981).The natural convection component h_{n} is calculated in the same way as the interior “Detailed” model. The detailed natural convection model correlates the convective heat transfer coefficient to the surface orientation and the difference between the surface and zone air temperatures (where DT = Air Temperature  Surface Temperature). The algorithm is taken directly from Walton (1983). Walton derived his algorithm from the ASHRAE Handbook (2001), Table 5 on p. 3.12, which gives equations for natural convection heat transfer coefficients in the turbulent range for large, vertical plates and for large, horizontal plates facing upward when heated (or downward when cooled). A note in the text also gives an approximation for large, horizontal plates facing downward when heated (or upward when cooled) recommending that it should be half of the facing upward value. Walton adds a curve fit as a function of the cosine of the tilt angle to provide intermediate values between vertical and horizontal. The curve fit values at the extremes match the ASHRAE values very well.
For no temperature difference OR a vertical surface the following correlation is used:
h=1.31ΔT13
For (DT < 0.0 AND an upward facing surface) OR (DT > 0.0 AND an downward facing surface) an enhanced convection correlation is used:
h=9.482ΔT137.283−cosΣ
where S is the surface tilt angle.
For (DT > 0.0 AND an upward facing surface) OR (DT < 0.0 AND an downward facing surface) a reduced convection correlation is used:
h=1.810ΔT131.382+cosΣ
where S is the surface tilt angle.
MoWiTT Algorithm[LINK]
Table 9. Nomenclature List of Variables.The MoWiTT model is based on measurements taken at the Mobile Window Thermal Test (MoWiTT) facility (Yazdanian and Klems 1994). The correlation applies to very smooth, vertical surfaces (e.g. window glass) in lowrise buildings and has the form:
hc=√[Ct(ΔT)13]2+[aVbz]2
Constants a, b and turbulent natural convection constant C_{t} are given in Table 10. The original MoWiTT model has been modified for use in EnergyPlus so that it is sensitive to the local suface’s wind speed which varies with the height above ground. The original MoWiTT model was formulated for use with the air velocity at the location of the weather station. As of Version 7.2, EnergyPlus uses the “a” model coefficients derived by Booten et al. (2012) rather than the original values from Yazdanian and Klems (1994).
NOTE: The MoWiTT algorithm may not be appropriate for rough surfaces, highrise surfaces, or surfaces that employ movable insulation.
Table 10. MoWiTT Coefficients (Yazdanian and Klems 1994, Booten et al. 2012)DOE2 Model[LINK]
Table 11. Nomenclature List of Variables.[LINK]
The DOE2 convection model is a combination of the MoWiTT and BLAST Detailed convection models (LBL 1994). The convection coefficient for very smooth surfaces (e.g. glass) is calculated as:
hc,glass=√h2n+[aVbz]2
h_{n} is calculated using Equation or Equation . Constants a and b are given in Table 10.
For less smooth surfaces, the convection coefficient is modified according to the equation
hc=hn+Rf(hc,glass−hn)
where R_{f} is the roughness multiplier given by Table 8.
Adaptive Convection Algorithm[LINK]
This algorithm has a structure that allows for finer control over the models used for particular surfaces. The algorithm for the outside face was developed for EnergyPlus but it borrows concepts and its name from the research done by BeausoleilMorrison (2000, 2002) for convection at the inside face (see the description below for interior convection).
The adaptive convection algorithm implemented in EnergyPlus for the outside face is much simpler than that for the inside face. The surface classification system has a total of 4 different categories for surfaces that depend on current wind direction and heat flow directions. However it is more complex in that the h_{c} equation is split into two parts and there are separate model equation selections for forced convection, h_{f}, and natural convection, h_{n}. The following table summarizes the categories and the default assignments for h_{c} equations. The individual h_{c} equations are documented below.
Table 12. Adaptive Convection Algorithm DetailsOutside Face Surface Classification[LINK]
During an initial setup phase, all the heat transfer surfaces in the input file are analyzed in groups to determine appropriate values for geometry scales used in many of the convection correlations. Eight separate groups are assembled for nominally vertical exterior surfaces for eight bins of azimuth: north, northeast, east, southeast, south, southwest, west, northwest. Surfaces with the same range of azimuth are grouped together and analyzed for overall geometry parameters. A ninth group is assembled for nominally horizontal exterior surfaces for a roof bin that is also analyzed for geometry. These geometry routines find bounds and limits of all the surfaces in the group and then model geometric parameters from these limits.
Sparrow Windward[LINK]
As discussed above for the TARP algorithm, a Sparrow et al. (1979) conducted flat plate measurements and develop the following correlation for finitesize flat plates oriented to windward.
hf=2.53Rf(PVzA)1/122
Sparrow Leeward[LINK]
Sparrow et al. (1979) conducted flat plate measurements and develop the following correlation for finitesize flat plates oriented to leeward.
hf=2.532Rf(PVzA)1/122
MoWITT Windward[LINK]
As discussed above, Yazdanian and Klems (1994) used outdoor laboratory measurements to develop the following correlation for smooth surfaces oriented to windward. Booten et al. (2012) developed revised coefficients for use with local surface wind speeds.
hc=√[0.84ΔT1/133]2+[2.38V0.89z]2
This model equation is for the total film coefficient and includes the natural convection portion. Therefore it should not be used in conjunction with a second natural convection model equation.
MoWITT Leeward[LINK]
Yazdanian and Klems (1994) used outdoor laboratory measurements to develop the following correlation for smooth surfaces oriented to leeward. Booten et al. (2012) developed revised coefficients for use with local surface wind speeds.
hc=√[0.84ΔT1/133]2+[2.86V0.617z]2
This model equation is for the total film coefficient and includes the natural convection portion. Therefore it should not be used in conjunction with a second natural convection model equation.
Blocken[LINK]
Blocken et al. (2009) developed a set of correlations for windward facing outdoor surfaces using numerical methods (key: BlockenWindward).
hf=4.6V0.8910m:θ≤11.25hf=5.0V0.8010m:11.25<θ≤33.75hf=4.6V0.8410m:33.75<θ≤56.25hf=4.5V0.8110m:56.25<θ≤100.0
Where V_{10m} is the air velocity at the location of the weather station and θ is the angle of incidence between the wind and the surface in degrees. This model is only applicable to windward surfaces and lacks a natural convection component and therefore cannot be used on its own but only within the adaptive convection algorithm for the outside face.
Clear[LINK]
Clear et al. (2003) developed correlations from measurements for horizontal roofs on two commercial buildings. In EnergyPlus the implementation uses the model for natural convection plus turbulent forced convection (eq. 8A in the reference) and applies it to the center point of each surface section that makes up the roof.
hc=ηkLn0.15Ra1/133Ln+kxRf0.0296Re4/455x1/133Pr
Where
x is the distance to the surface centroid from where the wind begins to intersect the roof. In EnergyPlus this is currently simplified to half the square root of the roof surface.
Ln=AreaPerimeter of overall roof
k is the thermal conductivity of air
$= is the weighting factor for natural convection (suppressed at high forced convection rates)
R{a_{{L_n}}} = G{r_{{L_n}}}$ is the Rayleigh number
GrLn=gρ2Ln3ΔTTfμ2 is the Grashof number
Rex=Vzρxμ is the Reynolds number at x
Pr is the Prandtl number
This model only claims to be applicable to horizontal roof surfaces so it may not be applicable to tilted roofs. It combines natural and forced convection and therefore should not be used in conjunction with yet another natural convection model.
Emmel[LINK]
Emmel et al. (2007) developed a set of correlations for outdoor surfaces using numerical methods. The following equations are for vertical surfaces (key: EmmelVertical):
hf=5.15V0.8110m:θ≤22.5hf=3.34V0.8410m:22.5<θ≤67.5hf=4.78V0.7110m:67.5<θ≤112.5hf=4.05V0.7710m:112.5<θ≤157.5hf=3.54V0.7610m:157.5<θ≤180.0
Where V_{10m} is the air velocity at the location of the weather station and θ is the angle of incidence between the wind and the surface in degrees. The following equations are used for horizontal (roof) surfaces (key: EmmelRoof):
hf=5.11V0.7810m:θ≤22.5hf=4.60V0.7910m:22.5<θ≤67.5hf=3.67V0.8510m:67.5<θ≤90
Where θ is the angle of incidence between the wind and the longest edge of the roof surface in degrees.
This model is for all wind directions but lacks a natural convection component. The model was developed for simple, rectangular lowrise buildings. It is available only within the adaptive convection algorithm for the outside face
Nusselt Jurges[LINK]
Perhaps the oldest equation for winddriven convection was developed by Nusselt and Jurges (1922). Palyvos (2008) casts their model in simplified form in SI units as:
hc=5.8+3.94Vz
Where V_{z} is the wind velocity in m/s, in EnergyPlus that velocity is adjusted for height above ground using the z axis coordinate of the surface’s centroid and the site wind model. This model can be applied to all surfaces and the relatively large constant is assumed to represent the natural convection portion of a total convection coefficient. The model is not sensitive to wind direction nor surface roughness.
McAdams[LINK]
A venerable equation for winddriven convection was developed by McAdams (1954) which Palyvos (2008) casts in SI units as:
hc=5.7+3.8Vz
Where V_{z} is the wind velocity in m/s that has been adjusted for height above ground using the z axis coordinate of the surface’s centroid. This model can be applied to all surfaces and the relatively large constant is assumed to represent the natural convection portion of a total convection coefficient. The model is not sensitive to wind direction nor surface roughness.
Mitchell[LINK]
A useful geometric scale based on building volume is used in an equation developed by Mitchell (1976). The winddriven convection equation is cast by Palyvos as:
hf=8.6V0.6zL0.4
Where V_{z} is the wind velocity in m/s that has been adjusted for height above ground using the z axis coordinate of the surface’s centroid and L is the cube root of the building’s total volume. EnergyPlus interprets this as the sum of the volume of all the zones in the input file.
Exterior/External Conduction[LINK]
The conduction term, q′′ko _{,} can in theory be calculated using a wide variety of heat conduction formulations. Typically in EnergyPlus, the Conduction Transfer Function (CTF) method is used. The available models are described in this section: Conduction Through The Walls.
References[LINK]
ASHRAE. 1981. 1981 ASHRAE Handbook – Fundamentals, Atlanta: American Society of Heating, Refrigerating, and AirConditioning Engineers, Inc.
ASHRAE. 1989. 1989 ASHRAE Handbook – Fundamentals, Atlanta: American Society of Heating, Refrigerating, and AirConditioning Engineers, Inc.
ASHRAE. 1993. 1993 ASHRAE Handbook – Fundamentals, Chapter 3, Heat Transfer, IP & SI Editions, Atlanta: American Society of Heating, Refrigerating, and AirConditioning Engineers, Inc.
ASHRAE. 2001. 2001 ASHRAE Handbook – Fundamentals, Atlanta: American Society of Heating, Refrigerating, and AirConditioning Engineers, Inc.
ASHRAE. 2005. 2005 ASHRAE Handbook – Fundamentals, Chapter 16, Air Flow Around Buildings, Atlanta: American Society of Heating, Refrigerating, and AirConditioning Engineers, Inc.
Booten, C., N. Kruis, and C. Christensen. 2012. Identifying and Resolving Issues in EnergyPlus and DOE2 Window Heat Transfer Calculations. National Renewable Energy Laboratory. NREL/TP550055787. Golden, CO.
Cole, R. J., and N. S. Sturrock. 1977. The Convective Heat Exchange at the External Surface of Buildings. Building and Environment, Vol. 12, p. 207.
Ellis, P.G., and P.A. Torcellini. 2005. “Simulating Tall Buildings Using EnergyPlus”, Proceedings of the Ninth International IBPSA Conference, Building Simulation 2005, Montreal, Canada, August 1518, 2005.
Lawrence Berkeley Laboratory (LBL). 1994. DOE2.1E053 source code.
Sparrow, E. M., J. W. Ramsey, and E. A. Mass. 1979. Effect of Finite Width on Heat Transfer and Fluid Flow about an Inclined Rectangular Plate. Journal of Heat Transfer, Vol. 101, p. 204.
U.S. Standard Atmosphere. 1976. U.S. Government Printing Office, Washington, D.C.
Walton, G. N. 1981. Passive Solar Extension of the Building Loads Analysis and System Thermodynamics (BLAST) Program, Technical Report, United States Army Construction Engineering Research Laboratory, Champaign, IL.
Walton, G. N. 1983. Thermal Analysis Research Program Reference Manual. NBSSIR 832655. National Bureau of Standards.
Yazdanian, M. and J. H. Klems. 1994. Measurement of the Exterior Convective Film Coefficient for Windows in LowRise Buildings. ASHRAE Transactions, Vol. 100, Part 1, p. 1087.
Documentation content copyright © 19962015 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.