The Doe2 Model

1.1.Weather Data

The 239 second-generation Typical Meteorological Year (TMY2) weather data files (Marion andUrban 1995), 42 first-generation Typical Meteorological Year (TMY) weather data files (NCDC1981), and four California Climate Zone (CTZ) weather data files listed in Table 2 are accessible in the Home Energy Saver. Corresponding location and climatic data (U.S. DOE 1998) are also shown for comparison. Of these parameters, only the elevation is not automatically read from weather data files by DOE-2; its value must be input explicitly, as listed in Table 1.

1.2.Materials and Constructions

Table 3 shows the physical properties of the construction materials that can be modeled in the Home Energy Saver. The constructions comprised of these materials are listed in Table 4. Of special note are interior partition walls, floors, and ceilings, which are modeled as half-constructions and doubled in area to account correctly for radiative heat transfer between interior surfaces.

Foundation constructions are handled directly in the DOE-2 BDL file. Two-dimensional response factors for all other constructions are pre-calculated using a program called WALFERF.These response factors are accessed and used by DOE-2 during simulation.Surface roughnesses are used in the calculation of exterior air film coefficients in DOE-2. Table5 shows the roughness codes (Winkelmann et al. 1993a) assumed for the exterior surfaces of the constructions listed in Table 4.Home Energy Saver users can specify the relative shades of exterior wall and roof surfaces. The solar absorptance values assigned to these shades are shown in Table 6.


The nominal door types shown in Table 7 can be modeled in the Home Energy Saver. The U-factors for these door types are drawn from ASHRAE (2001a). Alternatively, users can specify customized door types by U-factor. These U-factors incorporate the exterior air film coefficient at ASHRAE winter design conditions, while DOE-2 calculates variable exterior air film coefficients according to the hourly weather data supplied. Therefore, the input door U-factors are adjusted in DOE-2 macros according to the equation

U' = 1 / [(1 / U) – (1 / ho)](1)

U=design coefficient of heat transfer for a door (Btu/h·ft2·°F)
ho=design exterior air film coefficient (Btu/h·ft2·°F)
U'=coefficient of heat transfer for a door, excluding the effect of the exterior air film coefficient (Btu/h·ft2·°F)

In accordance with the most up-to-date calculation, the design exterior air film coefficient is assigned a value of 5.112 Btu/h·ft2·°F (Köhler 2002). All doors specified in the Home Energy Saver are assumed to be 6.67 ft high and 3 ft wide.

1.4.Windows and Skylights.

The nominal window and skylight types shown in Table 8 can be modeled in the Home Energy Saver. As with the doors, their properties are drawn from ASHRAE (2001a). Users can also specify customized window and skylight types by U-factor and solar heat gain coefficient(SHGC). The U-factors are adjusted in DOE-2 macros in the same way as those of doors.Because DOE-2 uses the shading coefficient (SC) as the measure of solar transmission through fenestration, the SHGCs must be adjusted in DOE-2 macros according to the equation

SC = SHGC / 0.87(2)


SHGC=whole-window or whole-skylight solar heat gain coefficient (dimensionless)
SC=whole-window or whole-skylight shading coefficient (dimensionless)


Slabs-on-grade, unconditioned and conditioned basements, and unvented and vented crawlspaces, with and without insulation, can be modeled in the Home Energy Saver, as shown in Table 10. Mixed or “averaged” foundation types cannot be simulated. Basements are assumed to be 8 ft high, with 1 ft above grade for standard basements and 4 ft above grade for raised basements in split-level houses. Crawlspaces are assumed to be 2.5 ft high and entirely above grade.Foundation heat transfer is simulated according to the method first devised by Huang et al.(1988) and subsequently updated (Winkelmann 1998; Huang 2003). This method utilizes effective foundation U-factors and fictitious foundation insulation layers to account for the insulating value of the surrounding soil. The effective U-factors for the foundation types modeled in the Home Energy Saver are included in Table 10.

1.6.Attics and Roofs

Unconditioned attics, conditioned attics, and cathedral or vaulted ceilings can be modeled in the Home Energy Saver. Mixed or “averaged” attic/roof types cannot be simulated.Roof apexes are always modeled parallel to the long axes of houses. Only one skylight type canbe modeled per roof.

1.7.House Geometry

Figure 1 indicates the necessary input dimensions for the eight types of house floor plans that can be modeled in the Home Energy Saver. Figure 2 provides further insight into the variety of house floor plans that can be treated.Different wall, door, and window types can be modeled on different sides of houses. However,there is no provision for modeling mixed types on the same side.Door and window areas are divided proportionately among the wall elements on a given facade. Window areas are also divided proportionately among the stories of a given façade. To avoid incorrect shading of doors and windows by one another, doors are always anchored at the left edges of walls, while windows are centered in walls. Window areas extend the full height of stories, resulting in offsetting errors in shading by other surfaces and objects.

Interior wall areas are taken to be half as great as conditioned floor areas, exclusive of basements and attics.

1.8.House Shading and Buffering

House shading by roof eaves, patios, carports, trellises, and similar overhangs can be modeled in the Home Energy Saver. Different extensions can be specified on different sides of houses. All overhangs are positioned at the roof line and assumed to run the lengths of the façades from which they extend. They are assigned a uniform solar transmittance of 0.DOE-2 test simulations were performed to ascertain the buffering effects of a garage on different sides of a small, single-story house in different climates. Reductions in heating and cooling energy consumption of 1.5% were consistently obtained. These effects should be even less forlarger houses. Because estimates of energy consumption from the Home Energy Saver are intended to be conservative, garages are not considered in the model.Shade trees and neighboring houses can be specified independently on different sides of housesin the Home Energy Saver. Both are centered with respect to house façades. Trees are considered to have 6-ft-tall trunks, assigned widths of 15 ft, and positioned 10 ft away from the houses they shade. Neighboring houses are assigned the same widths as the houses they shade and positioned20 ft away. Trees are assigned a solar transmittance of 0.70.

1.9.Internal Loads

House occupants are assumed to release sensible and latent loads of 230 Btu/h and 190 Btu/h,respectively, in the Home Energy Saver. Internal gains from appliances and lights are passed toDOE-2 from elsewhere in the Homer Energy Saver (Pinckard et al. 2003). Gains from appliances are assumed to be 80% sensible and 20% latent. The typical load profiles plotted in Figure 3apply to all 365 days of the simulated year.1.10.Infiltration and Natural Ventilation The Sherman-Grimsrud Method (Sherman 1980; Winkelmann et al. 1993b) is the specified infiltration calculation option in the Home Energy Saver. Terrain parameters for suburban environments are used. Fractional leakage areas (FLAs) for conditioned spaces above grade areselected outside of DOE-2 according to user-specified house airtightness, vintage, and number of stories (Pinckard et al. 2003) and passed to DOE-2 through a macro input variable. Alternatively,users can specify measured or estimated air leakage rates. These leakage rates are converted to SI units in a DOE-2 macro according to the simple relationship

Q50 = 0.0004719474 Q50ip(3)


Q50ip=conditioned house envelope air leakage rate at 50 Pa (ft3/min)
Q50=conditioned house envelope air leakage rate at 50 Pa (m3/s)

Conditioned house envelope air leakage rates are then found using the equation

Q4 = Q50 (4 / 50)0.65(4)

whereQ50=conditioned house envelope air leakage rate at 50 Pa (m3/s)
Q4=conditioned house envelope air leakage rate at 4 Pa (m3/s)

Effective leakage areas (ELAs) of conditioned house envelopes are determined using the equation

ELA = Q4 / [4 (2 / ρ)]0.5(5)

whereQ4=conditioned house envelope air leakage rate at 4 Pa (m3/s)
ρ=air density (1.2 kg/m3)
ELA=effective leakage area of a conditioned house envelope (dimensionless)

Lastly, ELAs are converted to FLAs, including a unit conversion, using the relationship

FLA = (ELA / 0.09290304) / Af(6)


ELA=effective leakage area of a conditioned house envelope (dimensionless)
Af=conditioned house floor area (ft2)
FLA=fractional leakage area of a conditioned house envelope (dimensionless)

For simplicity, basements are assumed to have no air exchange with the exterior. Unvented and vented crawlspaces are assigned FLAs of 0.0015 and 0.0030, respectively. Unconditioned atticsare assigned FLAs of 0.00242 (Treidler 1993).The opening of windows for natural ventilation is simulated whenever the exterior temperature and humidity would result in a cooling effect. One quarter of the total user-specified window area is assumed open, and a discharge coefficient of 0.6 is applied.

1.11.HVAC Systems

The HVAC system types shown in Table 11 can be modeled in the Home Energy Saver.Performance curves for equipment efficiency as a function of part-load ratio are drawn from Henderson et al. (1999). Curves for equipment capacity and efficiency as functions of outdoor temperature are unpublished LBNL derivations based on commercially available equipment performance data.Hot water heating is simulated in DOE-2 for tankless and indirect boiler types. The results for system types simulated as other system types are converted after capture in the Web programming as necessary.Window-mounted air conditioners that are controlled manually (i.e., by on-off switches rather than thermostats), whole-house fans, portable fans, and portable electric resistance heaters aretreated separately from DOE-2 as appliances (Pinckard et al. 2003).

1.12.Ducts and Pipes

A simplified estimate of duct delivery efficiency, based on draft ASHRAE Standard 152P(ASHRAE 2001b), is currently used in the Home Energy Saver. This single numerical value is a weighted average corresponding to the pre-estimated annual heating and cooling loads for the house to be modeled. The development of a DOE-2 function for the hourly calculation of ductv efficiency, based on the same Standard, is proposed.Boiler pipes are simply treated. Uninsulated and insulated pipes are assigned delivery efficiencies of 0.90 and 0.95, respectively.