Window Calculation Module[LINK]
This section describes two potential modeling approaches for Windows. The first (layer by layer) is implemented. The second, simple approach, reuses the layerbylayer approach but converts an arbitrary window performance into an equivalent single layer.
The primary Window calculation is a layerbylayer approach where windows are considered to be composed of the following components, only the first of which, glazing, is required to be present:
Glazing, which consists of one or more plane/parallel glass layers. If there are two or more glass layers, the layers are separated by gaps filled with air or another gas. The glazing optical and thermal calculations are based on algorithms from the WINDOW 4 and WINDOW 5 programs [Arasteh et al., 1989], [Finlayson et al., 1993]. Glazing layers are described using te input object WindowMaterial:Glazing.
Gap, layers filled with air or another gas that separate glazing layers. Gaps are described using the input object WindowMaterial:Gas.
Frame, which surrounds the glazing on four sides. Frames are described using the input object WindowProperty:FrameAndDivider.
Divider, which consists of horizontal and/or vertical elements that divide the glazing into individual lites.
Shading device, which is a separate layer, such as drapery, roller shade or blind, on the inside or outside of the glazing, whose purpose is to reduce solar gain, reduce heat loss (movable insulation) or control daylight glare. Shading layers are described using “WindowShadingControl” input objects.
In the following, the description of the layerbylayer glazing algorithms is based on material from Finlayson et al., 1993. The frame and divider thermal model, and the shading device optical and thermal models, are new to EnergyPlus.
A second approch has been developed where windows are modeled in a simplified approach that requires minimal user input that is processed to develop and equivalent layer that then reuses much of the layerbymodel. This “Simple Window Construction: model is described below.
Optical Properties of Glazing[LINK]
The solar radiation transmitted by a system of glass layers and the solar radiation absorbed in each layer depends on the solar transmittance, reflectance and absorptance properties of the individual layers. The absorbed solar radiation enters the glazing heat balance calculation that determines the inside surface temperature and, therefore, the heat gain to the zone from the glazing (see “Window Heat Balance Calculation“). The transmitted solar radiation is absorbed by interior zone surfaces and, therefore, contributes to the zone heat balance. In addition, the visible transmittance of the glazing is an important factor in the calculation of interior daylight illuminance from the glazing.
Variables in Window Calculations
T 
Transmittance 
 
 
R 
Reflectance 
 
 
Rf, Rb

Front reflectance, back reflectance 
 
 
Ti,j

Transmittance through glass layers i to j 
 
 
Tdirgl

Direct transmittance of glazing 
 
 
Rfi,j, Rbi,j

Front reflectance, back reflectance from glass layers i to j 
 
 
Rdirgl,f, Rdirgl,b

Direct front and back reflectance of glazing 
 
 
Afi, Abi

Front absorptance, back absorptance of layer i 
 
 
N 
Number of glass layers 
 
Nlayer 
λ 
Wavelength 
microns 
Wle 
Es(λ) 
Solar spectral irradiance function 
W/m2micron 
E 
V(λ) 
Photopic response function of the eye 
 
y30 
φ 
Angle of incidence (angle between surface normal and direction of incident beam radiation) 
Rad 
Phi 
τ 
Transmittivity or transmittance 
 
tf0 
ρ 
Reflectivity or reflectance 
 
rf0, rb0 
α 
Spectral absorption coefficient 
m−1

 
d 
Glass thickness 
M 
Material.Thickness 
n 
Index of refraction 
 
ngf, ngb 
κ 
Extinction coefficient 
 
 
β 
Intermediate variable 
 
betaf, betab 
P, p 
A general property, such as transmittance 
 
 
τsh

Shade transmittance 
 
Material.Trans 
ρsh

Shade reflectance 
 
Material.ReflectShade 
αsh

Shade absorptance 
 
Material.AbsorpSolar 
τbl ρbl αbl

Blind transmittance, reflectance, absorptance 
 
 
Q, G, J 
Source, irradiance and radiosity for blind optical properties calculation 
W/m2

 
Fij

View factor between segments i and j 
 
 
fswitch

Switching factor 
 
SwitchFac 
T 
Transmittance 
 
 
R 
Reflectance 
 
 
Rf, Rb

Front reflectance, back reflectance 
 
 
Ti,j

Transmittance through glass layers i to j 
 
 
Rfi,j, Rbi,j

Front reflectance, back reflectance from glass layers i to j 
 
 
Afi, Abi

Front absorptance, back absorptance of layer i 
 
 
N 
Number of glass layers 
 
Nlayer 
λ 
Wavelength 
microns 
Wle 
Es(λ)

Solar spectral irradiance function 
W/m2micron 
E 
R(λ) 
Photopic response function of the eye 
 
y30 
φ’ 
Relative azimuth angle (angle between screen surface normal and vertical plane through sun, Ref. Figure 87) 
Rad 
SunAzimuthToScreenNormal 
α’ 
Relative altitude angle (angle between screen surface horizontal normal plane and direction of incident beam radiation, Ref. Figure 87) 
Rad 
SunAltitudeToScreenNormal 
ρsc

Beamtodiffuse solar reflectance of screen material 
 
Screens.ReflectCylinder 
γ 
Screen material aspect ratio 
 
Screens.ScreenDiameterTo SpacingRatio 
Α 
Spectral absorption coefficient 
m−1

 
D 
Glass thickness 
M 
Material.Thickness 
N 
Index of refraction 
 
ngf, ngb 
Κ 
Extinction coefficient 
 
 
Β 
Intermediate variable 
 
betaf, betab 
P, p 
A general property, such as transmittance 
 
 
Glass Layer Properties[LINK]
In EnergyPlus, the optical properties of individual glass layers are given by the following quantities at normal incidence as a function of wavelength:
Transmittance, T
Front reflectance, Rf
Back reflectance, Rb
Here “front” refers to radiation incident on the side of the glass closest to the outside environment, and “back” refers to radiant incident on the side of the glass closest to the inside environment. For glazing in exterior walls, “front” is therefore the side closest to the outside air and “back” is the side closest to the zone air. For glazing in interior (i.e., interzone) walls, “back” is the side closest to the zone in which the wall is defined in and “front” is the side closest to the adjacent zone.
Glass Optical Properties Conversion[LINK]
Conversion from Glass Optical Properties Specified as Index of Refraction and Transmittance at Normal Incidence[LINK]
The optical properties of uncoated glass are sometimes specified by index of refraction, n,* * and transmittance at normal incidence, T.
The following equations show how to convert from this set of values to the transmittance and reflectance values required by Material:WindowGlass. These equations apply only to uncoated glass, and can be used to convert either spectralaverage solar properties or spectralaverage visible properties (in general, n and T are different for the solar and visible). Note that since the glass is uncoated, the front and back reflectances are the same and equal to the R that is solved for in the following equations.
Given n and T, find R:
r=(n−1n+1)2τ=[(1−r)4+4r2T2]1/2−(1−r)22r2TR=r+(1−r)2rτ21−r2τ2
Example:
T=0.86156n=1.526r=(1.526−11.526+1)2τ=0.93974R=0.07846
Simple Window Model[LINK]
EnergyPlus includes an alternate model that allows users to enter in simplified window performance indices. This model is accessed through the WindowMaterial:SimpleGlazingSystem input object and converts the simple indices into an equivalent single layer window. (In addition a special model is used to determine the angular properties of the system – described below). Once the model generates the properties for the layer, the program reuses the bulk of the layerbylayer model for subsequent calculations. The properties of the equivalent layer are determined using the step by step method outlined by Arasteh, Kohler, and Griffith (2009). The core equations are documented here. The reference contains additional information.
The simplified window model accepts U and SHGC indices and is useful for several reasons:
1) Sometimes, the only thing that is known about the window are its U and SHGC;
2) Codes, standards, and voluntary programs are developed in these terms;
3) A singlelayer calculation is faster than multilayer calculations.
Note: This use of U and SHGC to describe the thermal properties of windows is only appropriate for specular glazings.
While it is important to include the ability to model windows with only Uvalue and SHGC, we note that any method to use U and SHGC alone in building simulation software will inherently be approximate. This is due primarily to the following factors:
SHGC combines directly transmitted solar radiation and radiation absorbed by the glass which flows inward. These have different implications for space heating/cooling. Different windows with the same SHGC often have different ratios of transmitted to absorbed solar radiation.
SHGC is determined at normal incidence; angular properties of glazings vary with number of layers, tints, coatings. So products which have the same SHGC, can have different angular properties.
Window Ufactors vary with temperatures.
Thus, for modeling specific windows, we recommend using more detailed data than just the U and SHGC, if at all possible.
The simplified window model determines the properties of an equivalent layer in the following steps.
Step 1. Determine glasstoglass Resistance.[LINK]
Window Uvalues include interior and exterior surface heat transfer coefficients. The resistance of the bare window product, or glasstoglass resistance is augmented by these film coefficients so that,
1U=Ri,w+Ro,w+Rl,w
Where,
Ri,w is the resistance of the interior film coefficient under standard winter conditions in units of m2·K/W,
Ro,w is the resistance of the exterior film coefficient under standard winter conditions in units of m2·K/W, and
Rl,w is the resisance of the bare window under winter conditions (without the film coefficients) in units of m2·K/W.
The values for Ri,w and Ro,w depend on U and are calculated using the following correlations.
Ri,w=1(0.359073Ln(U)+6.949915);forU<5.85
Ri,w=1(1.788041U−2.886625);forU≥5.85
Ro,w=1(0.025342U+29.163853)
So that the glasstoglass resistance is calculated using:
Rl,w=1U−Ri,w−Ro,w
Because the window model in EnergyPlus is for flat geometries, the models are not necessarily applicable to lowperformance projecting products, such as skylights with unisulated curbs. The model cannot support glazing systems with a U higher than 7.0 because the thermal resistance of the film coefficients alone can provide this level of performance and none of the various resistances can be negative.
Step 2. Determine Layer Thickness.[LINK]
The thickness of the equivalent layer in units of meters is calculated using,
Thickness=⎧⎪⎨⎪⎩0.002for 1Rl,w>7.00.05914−0.00714Rl,wfor 1Rl,w≤7.0
Step 3. Determine Layer Thermal Conductivity[LINK]
The effective thermal conductivity, λeff, of the equivalent layer is calculated using,
λeff=ThicknessRl,w
Step 4. Determine Layer Solar Transmittance[LINK]
The layer’s solar transmittance at normal incidence, Tsol, is calculated using correlations that are a function of SHGC and UFactor.
Tsol=0.939998SHGC2+0.20332SHGC;U>4.5;SHGC<0.7206
Tsol=1.30415SHGC−0.30515;U>4.5;SHGC≥0.7206
Tsol=0.41040SHGC;U<3.4;SHGC≤0.15
Tsol=0.085775SHGC2+0.963954SHGC−0.084958;U<3.4;SHGC>0.15
And for Uvalues between 3.4 and 4.5, the value for Tsol is interpolated using results of the equations for both ranges.
Step 5. Determine Layer Solar Reflectance[LINK]
The layer’s solar reflectance is calculated by first determining the inward flowing fraction which requires values for the resistance of the inside and outside film coefficients under summer conditions, Ri,s and Ro,s, respectively. The correlations are:
Ri,s=1(29.436546(SHGC−TSol)3−21.943415(SHGC−TSol)2+9.945872(SHGC−TSol)+7.426151);U>4.5Ri,s=1(199.8208128(SHGC−TSol)3−90.639733(SHGC−TSol)2+19.737055(SHGC−TSol)+6.766575);U<3.4Ro,s=1(2.225824(SHGC−TSol)+20.57708);U>4.5Ro,s=1(5.763355(SHGC−TSol)+20.541528);U<3.4
And for Uvalues between 3.4 and 4.5, the values are interpolated using results from both sets of equations.
The inward flowing fraction, Fracinward, is then calculated using:
Fracinward=(Ro,s+0.5Rl,w)(Ro,s+Rl,w+Ri,s)
Then, the solar reflectances of the front face, Rs,f, and back face, Rs,b, are calculated using:
Rs,f=Rs,b=1−TSol−(SHGC−TSol)Fracinward
The thermal absorptance, or emittance, is taken as 0.84 for both the front and back and the longwave transmittance is 0.0.
Step 6. Determine Layer Visible Properties[LINK]
The user has the option of entering a value for visible transmittance as one of the simple performance indices. If the user does not enter a value, then the visible properties are the same as the solar properties. If the user does enter a value then layer’s visible transmittance at normal incidence, TVis, is set to that value. The visible light reflectance for the back surface is calculated using:
RVis,b=−0.7409T3Vis+1.6531T2Vis−1.2299TVis+0.4547
The visible light reflectance for the front surface is calculated using:
RVis,f=−0.0622T3Vis+0.4277T2Vis−0.4169TVis+0.2399
The angular properties of windows are important because during energy modeling, the solar incidence angles are usually fairly high. Angles of incidence are defined as angles from the normal direction extending out from the window. The simple glazing system model includes a range of correlations that are selected based on the values for U and SHGC. These were chosen to match the types of windows likely to have such performance levels. The matrix of possible combinations of U and SHGC values have been mapped to set of 28 bins shown in the following figure.
There are ten different correlations, A thru J, for both transmission and reflectance. The correlations are used in various weighting and interpolation schemes according the figure above. The correlations are normalized against the performance at normal incidence. EnergyPlus uses these correlations to store the glazing system’s angular performance at 10 degree increments and interpolates between them during simulations. The model equations use the cosine of the incidence angle, cos(φ), as the independent variable. The correlations have the form:
TorR=acos(φ)4+bcos(φ)3+ccos(φ)2+dcos(φ)+e
The coefficient values for a, b, c, d, and e are listed in the following tables for each of the curves.
Application Issues[LINK]
EnergyPlus’s normal process of running the detailed layerbylayer model, with the equivalent layer produced by this model, creates reports (sent to the EIO file) of the overall performance indices and the properties of the equivalent layer. Both of these raise issues that may be confusing.
The simplified window model does not reuse all aspects of the detailed layerbylayer model, in that the angular solar transmission properties use a different model when the simple window model is in effect. If the user takes the material properties of an equivalent glazing layer from the simple window model and then reenters them into just the detailed model, then the performance will not be the same because of the angular transmission model will have changed. It is not proper use of the model to reenter the equivalent layer’s properties and expect the exact level of performance.
There may not be exact agreement between the performance indices echoed out and those input in the model. This is expected with the model and the result of a number of factors. For example, although input is allowed to go up to U7 W/m2⋅K, the actual outcome is limited to no higher than about 5.8W/m2⋅K. This is because the thermal resistance to heat transfer at the surfaces is already enough resistance to provide an upper limit to the conductance of a planar surface. Sometimes there is conflict between the SHGC and the U that are not physical and compromises need to be made. In general, the simple window model is intended to generate a physicallyreasonable glazing that approximates the input entered as well as possible. But the model is not always be able to do exactly what is specified when the specifications are not physical.
Arasteh, D., J.C. Kohler, B. Griffith, Modeling Windows in EnergyPlus with Simple Performance Indices. Lawrence Berkeley National Laboratory. In Draft. Available at
Glazing System Properties[LINK]
The optical properties of a glazing system consisting of N glass layers separated by nonabsorbing gas layers (see Figure 4) are determined by solving the following recursion relations for Ti,j, the transmittance through layers i to j; Rfi,j and Rbi,j, the front and back reflectance, respectively, from layers i to j; and Aj, the absorption in layer j. Here layer 1 is the outermost layer and layer N is the innermost layer. These relations account for multiple internal reflections within the glazing system. Each of the variables is a function of wavelength.
Ti,j=Ti,j−1Tj,j1−Rfj,jRbj−1,i
Rfi,j=Rfi,j−1+T2i,j−1Rfj,j1−Rfj,jRbj−1,i
Rbj,i=Rbj,j+T2j,jRbj−1,i1−Rbj−1,iRfj,j
Afj=T1,j−1(1−Tj,j−Rfj,j)1−Rfj,NRbj−1,1+T1,jRfj+1,N(1−Tj,j−Rbj,j)1−Rfj,NRbj−1,1
In Equation [eq:Ajtothefequation], Ti,j = 1 and Ri,j = 0 if i<0 or j>N.
As an example, for double glazing (N=2), these equations reduce to:
T1,2=T1,1T2,21−Rf2,2Rb1,1
Rf1,2=Rf1,1+T21,1Rf2,21−Rf2,2Rb1,1
Rb2,1=Rb2,2+T22,2Rb1,11−Rb1,1Rf2,2
Af1=(1−T1,1−Rf1,1)+T1,1Rf2,2(1−T1,1−Rb1,1)1−Rf2,2Rb1,1
Af2=T1,1(1−T2,2−Rf2,2)1−Rf2,2Rb1,1
If the above transmittance and reflectance properties are input as a function of wavelength, EnergyPlus calculates “spectral average” values of the above glazing system properties by integrating over wavelength.
The spectralaverage solar property is:
Ps=∫P(λ)Es(λ)dλ∫Es(λ)dλ
The spectralaverage visible property is:
Pv=∫P(λ)Es(λ)V(λ)dλ∫Es(λ)V(λ)dλ
where Es(λ) is the solar spectral irradiance function and V(λ) is the photopic response function of the eye. The default functions are shown in Table 2 and Table 3. They can be overwritten by user defined solar and/or visible spectrum using the objects Site:SolarAndVisibleSpectrum and Site:SpectrumData. They are expressed as a set of values followed by the corresponding wavelengths for values.
When a choice of Spectral is entered as the optical data type, the correlations to store the glazing system’s angular performance are generated based on angular performance at 10 degree increments. When a choice of SpectralAndAngle is entered as the optical data type, the correlations for the glazing system will be generated using 10 degree increments or more if the SpectralAndAngle properties include data for more angles. For each incident angle, the properties of the SpectralAndAngle layer(s) is calculated by linear interpolation, and then the performance of the entire glazing system is calculated for that angle. The glazing system properties at each angle are used to generate polynomial curve fits with 6 coefficients as a function of cosine of incident angle. The polynomial curves are then used in the simulation to calculate optical properties at each timestep.
If a glazing layer has optical properties that are roughly constant with wavelength, the wavelengthdependent values of Ti,i, Rfi,i and Rbi,i in Equations [eq:Tijequation] to [eq:Ajtothefequation] can be replaced with constant values for that layer.
0.0, 9.5, 42.3, 107.8, 181.0, 246.0, 395.3, 390.1, 435.3, 438.9,
483.7, 520.3, 666.2, 712.5, 720.7, 1013.1, 1158.2, 1184.0, 1071.9, 1302.0,
1526.0, 1599.6, 1581.0, 1628.3, 1539.2, 1548.7, 1586.5, 1484.9, 1572.4, 1550.7,
1561.5, 1501.5, 1395.5, 1485.3, 1434.1, 1419.9, 1392.3, 1130.0, 1316.7, 1010.3,
1043.2, 1211.2, 1193.9, 1175.5, 643.1, 1030.7, 1131.1, 1081.6, 849.2, 785.0,
916.4, 959.9, 978.9, 933.2, 748.5, 667.5, 690.3, 403.6, 258.3, 313.6,
526.8, 646.4, 746.8, 690.5, 637.5, 412.6, 108.9, 189.1, 132.2, 339.0,
460.0, 423.6, 480.5, 413.1, 250.2, 32.5, 1.6, 55.7, 105.1, 105.5,
182.1, 262.2, 274.2, 275.0, 244.6, 247.4, 228.7, 244.5, 234.8, 220.5,
171.5, 30.7, 2.0, 1.2, 21.2, 91.1, 26.8, 99.5, 60.4, 89.1,
82.2, 71.5, 70.2, 62.0, 21.2, 18.5, 3.2
0.3000, 0.3050, 0.3100, 0.3150, 0.3200, 0.3250, 0.3300, 0.3350, 0.3400, 0.3450,
0.3500, 0.3600, 0.3700, 0.3800, 0.3900, 0.4000, 0.4100, 0.4200, 0.4300, 0.4400,
0.4500, 0.4600, 0.4700, 0.4800, 0.4900, 0.5000, 0.5100, 0.5200, 0.5300, 0.5400,
0.5500, 0.5700, 0.5900, 0.6100, 0.6300, 0.6500, 0.6700, 0.6900, 0.7100, 0.7180,
0.7244, 0.7400, 0.7525, 0.7575, 0.7625, 0.7675, 0.7800, 0.8000, 0.8160, 0.8237,
0.8315, 0.8400, 0.8600, 0.8800, 0.9050, 0.9150, 0.9250, 0.9300, 0.9370, 0.9480,
0.9650, 0.9800, 0.9935, 1.0400, 1.0700, 1.1000, 1.1200, 1.1300, 1.1370, 1.1610,
1.1800, 1.2000, 1.2350, 1.2900, 1.3200, 1.3500, 1.3950, 1.4425, 1.4625, 1.4770,
1.4970, 1.5200, 1.5390, 1.5580, 1.5780, 1.5920, 1.6100, 1.6300, 1.6460, 1.6780,
1.7400, 1.8000, 1.8600, 1.9200, 1.9600, 1.9850, 2.0050, 2.0350, 2.0650, 2.1000,
2.1480, 2.1980, 2.2700, 2.3600, 2.4500, 2.4940, 2.5370
: Solar spectral irradiance function.
[table:solarspectralirradiancefunction.]
0.0000, 0.0001, 0.0001, 0.0002, 0.0004, 0.0006, 0.0012, 0.0022, 0.0040, 0.0073,
0.0116, 0.0168, 0.0230, 0.0298, 0.0380, 0.0480, 0.0600, 0.0739, 0.0910, 0.1126,
0.1390, 0.1693, 0.2080, 0.2586, 0.3230, 0.4073, 0.5030, 0.6082, 0.7100, 0.7932,
0.8620, 0.9149, 0.9540, 0.9803, 0.9950, 1.0000, 0.9950, 0.9786, 0.9520, 0.9154,
0.8700, 0.8163, 0.7570, 0.6949, 0.6310, 0.5668, 0.5030, 0.4412, 0.3810, 0.3210,
0.2650, 0.2170, 0.1750, 0.1382, 0.1070, 0.0816, 0.0610, 0.0446, 0.0320, 0.0232,
0.0170, 0.0119, 0.0082, 0.0158, 0.0041, 0.0029, 0.0021, 0.0015, 0.0010, 0.0007,
0.0005, 0.0004, 0.0002, 0.0002, 0.0001, 0.0001, 0.0001, 0.0000, 0.0000, 0.0000,
0.0000
.380, .385, .390, .395, .400, .405, .410, .415, .420, .425,
.430, .435, .440, .445, .450, .455, .460, .465, .470, .475,
.480, .485, .490, .495, .500, .505, .510, .515, .520, .525,
.530, .535, .540, .545, .550, .555, .560, .565, .570, .575,
.580, .585, .590, .595, .600, .605, .610, .615, .620, .625,
.630, .635, .640, .645, .650, .655, .660, .665, .670, .675,
.680, .685, .690, .695, .700, .705, .710, .715, .720, .725,
.730, .735, .740, .745, .750, .755, .760, .765, .770, .775,
.780
: Photopic response function.
[table:photopicresponsefunction.]
Calculation of Angular Properties[LINK]
Calculation of optical properties is divided into two categories: uncoated glass and coated glass.
Angular Properties for Uncoated Glass[LINK]
The following discussion assumes that optical quantities such as transmissivity, reflectvity, absorptivity, and index of refraction are a function of wavelength, λ. If there are no spectral data the angular dependence is calculated based on the single values for transmittance and reflectance in the visible and solar range. In the visible range an average wavelength of 0.575 microns is used in the calculations. In the solar range an average wavelength of 0.898 microns is used.
The spectral data include the transmittance, T, and the reflectance, R. For uncoated glass the reflectance is the same for the front and back surfaces. For angle of incidence, ϕ , the transmittance and reflectance are related to the transmissivity, τ, and reflectivity, ρ, by the following relationships:
T(ϕ)=τ(ϕ)2e−αd/cosϕ′1−ρ(ϕ)2e−2αd/cosϕ′
R(ϕ)=ρ(ϕ)(1+T(ϕ)e−αd/cosϕ′)
The spectral reflectivity is calculated from Fresnel’s equation assuming unpolarized incident radiation:
ρ(ϕ)=12((ncosϕ−cosϕ′ncosϕ+cosϕ′)2+(ncosϕ′−cosϕncosϕ′+cosϕ)2)
The spectral transmittivity is given by:
τ(ϕ)=1−ρ(ϕ)
The spectral absorption coefficient is defined as:
α=4πκλ
where κ is the dimensionless spectrallydependent extinction coefficient and λ is the wavelength expressed in the same units as the sample thickness.
Solving Equation [eq:RhoofPhi] at normal incidence gives:
n=1+√ρ(0)1−√ρ(0)
Evaluating Equation [eq:RofPhiEquation] at normal incidence gives the following expression for κ:
κ=−λ4πdlnR(0)−ρ(0)ρ(0)T(0)
Eliminating the exponential in Equations [eq:TofPhiEquation] and [eq:RofPhiEquation] gives the reflectivity at normal incidence:
ρ(0)=β−√β2−4(2−R(0))R(0)2(2−R(0))
where
β=T(0)2−R(0)2+2R(0)+1
The value for the reflectivity, ρ(0), from Equation [eq:Rho0Equation] is substituted into Equations [eq:NAsFunctionOfRhoEquation] and [eq:KappaAsFunctionOfLambdaRTEquation]. The result from Equation [eq:KappaAsFunctionOfLambdaRTEquation] is used to calculate the absorption coefficient in Equation [eq:AlphaAsFunctionOfKappaLambda]. The index of refraction is used to calculate the reflectivity in Equation [eq:RhoofPhi] which is then used to calculate the transmittivity in Equation [eq:TauofPhiEquation]. The reflectivity, transmissivity and absorption coefficient are then substituted into Equations [eq:TofPhiEquation] and [eq:RofPhiEquation] to obtain the angular values of the reflectance and transmittance.
Angular Properties for Coated Glass[LINK]
A regression fit is used to calculate the angular properties of coated glass from properties at normal incidence. If the transmittance of the coated glass is > 0.645, the angular dependence of uncoated clear glass is used. If the transmittance of the coated glass is ≤ 0.645, the angular dependence of uncoated bronze glass is used. The values for the angular functions for the transmittance and reflectance of both clear glass (¯τclr,¯ρclr) and bronze glass (¯τbnz,¯ρbnz) are determined from a fourthorder polynomial regression:
¯τ(ϕ)=¯τ0+¯τ1cos(ϕ)+¯τ2cos2(ϕ)+¯τ3cos3(ϕ)+¯τ4cos4(ϕ)
and
¯ρ(ϕ)=¯ρ0+¯ρ1cos(ϕ)+¯ρ2cos2(ϕ)+¯ρ3cos3(ϕ)+¯ρ4cos4(ϕ)−¯τ(ϕ)
The polynomial coefficients are given in Table 4.
Polynomial coefficients used to determine angular properties of coated glass.
¯τclr 
0.0015 
3.355 
3.840 
1.460 
0.0288 
¯ρclr 
0.999 
0.563 
2.043 
2.532 
1.054 
¯τbnz 
0.002 
2.813 
2.341 
0.05725 
0.599 
¯ρbnz 
0.997 
1.868 
6.513 
7.862 
3.225 
These factors are used as follows to calculate the angular transmittance and reflectance:
For T(0) > 0.645:
T(ϕ)=T(0)¯τclr(ϕ)
R(ϕ)=R(0)(1−¯ρclr(ϕ))+¯ρclr(ϕ)
For T(0) ≤ 0.645:
T(ϕ)=T(0)¯τbnz(ϕ)
R(ϕ)=R(0)(1−¯ρbnz(ϕ))+¯ρbnz(ϕ)
Angular Properties for Simple Glazing Systems[LINK]
When the glazing system is modeled using the simplified method, an alternate method is used to determine the angular properties. The equation for solar transmittance as a function of incidence angle, T(ϕ), is:
T(ϕ)=T(ϕ=0)cos(ϕ)(1+(0.768+0.817SHGC4)sin3(ϕ))
where,
T(ϕ=0) is the normal incidence solar transmittance, TSol.
The equation for solar reflectance as a function of incidence angle, R(ϕ), is:
R(ϕ)=R(ϕ=0)(f1(ϕ)+f2(ϕ)√SHGC)Rfit,o
where,
f1(ϕ)=(((2.403cos(ϕ)−6.192)cos(ϕ)+5.625)cos(ϕ)−2.095)cos(ϕ)+1
f2(ϕ)=(((−1.188cos(ϕ)+2.022)cos(ϕ)+0.137)cos(ϕ)−1.720)cos(ϕ)
Rfit,o=0.7413−(0.7396√SHGC)
Calculation of Hemispherical Values[LINK]
The hemispherical value of a property is determined from the following integral:
Phemispherical=2∫π20P(ϕ)cos(ϕ)sin(ϕ)dϕ
The integral is evaluated by Simpson’s rule for property values at angles of incidence from 0 to 90 degrees in 10degree increments.
Optical Properties of Window Shading Devices[LINK]
Shading devices affect the system transmittance and glass layer absorptance for shortwave radiation and for longwave (thermal) radiation. The effect depends on the shade position (interior, exterior or betweenglass), its transmittance, and the amount of interreflection between the shading device and the glazing. Also of interest is the amount of radiation absorbed by the shading device.
In EnergyPlus, shading devices are divided into four categories, “shades,” “blinds,” “screens,” and “switchable glazing.” “Shades” are assumed to be perfect diffusers. This means that direct radiation incident on the shade is reflected and transmitted as hemispherically uniform diffuse radiation: there is no direct component of transmitted radiation. It is also assumed that the transmittance, τsh, reflectance, ρsh, and absorptance, αsh, are the same for the front and back of the shade and are independent of angle of incidence. Many types of drapery and pulldown roller devices are close to being perfect diffusers and can be categorized as “shades.”
“Blinds” in EnergyPlus are slattype devices such as venetian blinds. Unlike shades, the optical properties of blinds are strongly dependent on angle of incidence. Also, depending on slat angle and the profile angle of incident direct radiation, some of the direct radiation may pass between the slats, giving a direct component of transmitted radiation.
“Screens” are debris or insect protection devices made up of metallic or nonmetallic materials. Screens may also be used as shading devices for large glazing areas where excessive solar gain is an issue. The EnergyPlus window screen model assumes the screen is composed of intersecting orthogonallycrossed cylinders, with the surface of the cylinders assumed to be diffusely reflecting. Screens may only be used on the exterior surface of a window construction. As with blinds, the optical properties affecting the direct component of transmitted radiation are dependent on the angle of incident direct radiation.
With “Switchable glazing,” shading is achieved making the glazing more absorbing or more reflecting, usually by an electrical or chemical mechanism. An example is electrochromic glazing where the application of an electrical voltage or current causes the glazing to switch from light to dark.
Shades and blinds can be either fixed or moveable. If moveable, they can be deployed according to a schedule or according to a trigger variable, such as solar radiation incident on the window. Screens can be either fixed or moveable according to a schedule.
Shade/Glazing System Properties for ShortWave Radiation[LINK]
Shortwave radiation includes:
Beam solar radiation from the sun and diffuse solar radiation from the sky and ground incident on the outside of the window,
Beam and/or diffuse radiation reflected from exterior obstructions or the building itself,
Solar radiation reflected from the inside zone surfaces and incident as diffuse radiation on the inside of the window,
Beam solar radiation from one exterior window incident on the inside of another window in the same zone, and
Shortwave radiation from electric lights incident as diffuse radiation on the inside of the window.
Exterior Shade
For an exterior shade we have the following expressions for the system transmittance, the effective system glass layer absorptance, and the system shade absorptance, taking interreflection between shade and glazing into account. Here, “system” refers to the combination of glazing and shade. The system properties are given in terms of the isolated shade properties (i.e., shade properties in the absence of the glazing) and the isolated glazing properties (i.e., glazing properties in the absence of the shade).
Tsys(ϕ)=Tdif1,Nτsh1−Rdiffρsh
Tdifsys=Tdif1,Nτsh1−Rdiffρsh
Asysj,f(ϕ)=Adifj,fτsh1−Rfρsh,j=1 to N
Adif,sysj,f=Adifj,fτsh1−Rfρsh,j=1 to N
Adif,sysj,b=Adifj,bTdif1,Nρsh1−Rfρsh,j=1 to N
αsyssh=αsh(1+τshRf1−Rfρsh)
Interior Shade
The system properties when an interior shade is in place are the following:
Tsys(ϕ)=T1,N(ϕ)τsh1−Rdifbρsh
Tdifsys=Tdif1,Nτsh1−Rdifbρsh
Asysj,f(ϕ)=Aj,f(ϕ)+T1,N(ϕ)ρsh1−RdifbρshAdifj,b,j=1 to N
Adif,sysj,f=Adifj,f+Tdif1,Nρsh1−RdifbρshAdifj,b,j=1 to N
Adif,sysj,b=τsh1−RdifbρshAdifj,b, j=1 to N
αsyssh(ϕ)=T1,N(ϕ)αsh1−Rdifbρsh
αdif,syssh=Tdif1,Nαsh1−Rdifbρsh
LongWave Radiation Properties of Window Shades[LINK]
Longwave radiation includes:
Thermal radiation from the sky, ground and exterior obstructions incident on the outside of the window,
Thermal radiation from other room surfaces incident on the inside of the window, and
Thermal radiation from internal sources, such as equipment and electric lights, incident on the inside of the window.
The program calculates how much longwave radiation is absorbed by the shade and by the adjacent glass surface. The system emissivity (thermal absorptance) for an interior or exterior shade, taking into account reflection of longwave radiation between the glass and shade, is given by:
εlw,syssh=εlwsh⎛⎝1+τlwshρlwgl1−ρlwshρlwgl⎞⎠
where ρlwgl is the longwave reflectance of the outermost glass surface for an exterior shade or the innermost glass surface for an interior shade, and it is assumed that the longwave transmittance of the glass is zero.
The innermost (for interior shade) or outermost (for exterior shade) glass surface emissivity when the shade is present is:
εlw,sysgl=εlwglτlwsh1−ρlwshρlwgl
Switchable Glazing[LINK]
For switchable glazing, such as electrochromics, the solar and visible optical properties of the glazing can switch from a light state to a dark state. The switching factor, fswitch, determines what state the glazing is in. An optical property, p, such as transmittance or glass layer absorptance, for this state is given by:
p=(1−fswitch)plight+fswitchpdark
where
plight is the property value for the unswitched, or light state, and pdark is the property value for the fully switched, or dark state.
The value of the switching factor in a particular time step depends on what type of switching control has been specified: “schedule,” “trigger,” or “daylighting.” If “schedule,” fswitch = schedule value, which can be 0 or 1.
Thermochromic Windows[LINK]
Thermochromic (TC) materials have active, reversible optical properties that vary with temperature. Thermochromic windows are adaptive window systems for incorporation into building envelopes. Thermochromic windows respond by absorbing sunlight and turning the sunlight energy into heat. As the thermochromic film warms it changes its light transmission level from less absorbing to more absorbing. The more sunlight it absorbs the lower the light level going through it. Figure 5 shows the variations of window properties with the temperature of the thermochromic glazing layer. By using the suns own energy the window adapts based solely on the directness and amount of sunlight. Thermochromic materials will normally reduce optical transparency by absorption and/or reflection, and are specular (maintaining vision).
On cloudy days the window is at full transmission and letting in diffuse daylighting. On sunny days the window maximizes diffuse daylighting and tints based on the angle of the sun relative to the window. For a south facing window (northern hemisphere) the daylight early and late in the day is maximized and the direct sun at mid day is minimized.
The active thermochromic material can be embodied within a laminate layer or a surface film. The overall optical state of the window at a given time is a function primarily of:
thermochromic material properties
solar energy incident on the window
construction of the window system that incorporates the thermochromic layer
environmental conditions (interior, exterior, air temperature, wind, etc).
The tinted film, in combination with a heat reflecting, lowe layer allows the window to reject most of the absorbed radiation thus reducing undesirable heat load in a building. In the absence of direct sunlight the window cools and clears and again allows lower intensity diffuse radiation into a building. TC windows can be designed in several ways (Figure 6), with the most common being a triple pane windows with the TC glass layer in the middle a double pane windows with the TC layer on the inner surface of the outer pane or for sloped glazing a double pane with the laminate layer on the inner pane with a lowe layer toward the interior. The TC glass layer has variable optical properties depending on its temperature, with a lower temperature at which the optical change is initiated, and an upper temperature at which a minimum transmittance is reached. TC windows act as passive solar shading devices without the need for sensors, controls and power supplies but their optical performance is dependent on varying solar and other environmental conditions at the location of the window.
EnergyPlus describes a thermochromic window with a Construction object which references a special layer defined with a WindowMaterial:GlazingGroup:Thermochromic object. The WindowMaterial:GlazingGroup:Thermochromic object further references a series of WindowMaterial:Glazing objects corresponding to each specification temperature of the TC layer. During EnergyPlus run time, a series of TC windows corresponding to each specification temperature is created once. At the beginning of a particular time step calculations, the temperature of the TC glass layer from the previous time step is used to look up the most closed specification temperature whose corresponding TC window construction will be used for the current time step calculations. The current time step calculated temperature of the TC glass layer can be different from the previous time step, but no iterations are done in the current time step for the new TC glass layer temperature. This is an approximation that considers the reaction time of the TC glass layer can be close to EnergyPlus simulation time step say 10 to 15 minutes.
Window blinds in EnergyPlus are defined as a series of equidistant slats that are oriented horizontally or vertically. All of the slats are assumed to have the same optical properties. The overall optical properties of the blind are determined by the slat geometry (width, separation and angle) and the slat optical properties (frontside and backside transmittance and reflectance). Blind properties for direct radiation are also sensitive to the “profile angle,” which is the angle of incidence in a plane that is perpendicular to the window plane and to the direction of the slats. The blind optical model in EnergyPlus is based on Simmler, Fischer and Winkelmann, 1996; however, that document has numerous typographical errors and should be used with caution.
The following assumptions are made in calculating the blind optical properties:
The slats are flat.
The spectral dependence of interreflections between slats and glazing is ignored; spectralaverage slat optical properties are used.
The slats are perfect diffusers. They have a perfectly matte finish so that reflection from a slat is isotropic (hemispherically uniform) and independent of angle of incidence, i.e., the reflection has no specular component. This also means that absorption by the slats is hemispherically uniform with no incidence angle dependence. If the transmittance of a slat is nonzero, the transmitted radiation is isotropic and the transmittance is independent of angle of incidence.
Interreflection between the blind and wall elements near the periphery of the blind is ignored.
If the slats have holes through which support strings pass, the holes and strings are ignored. Any other structures that support or move the slats are ignored.
Slat Optical Properties[LINK]
The slat optical properties used by EnergyPlus are shown in the following table.
Slat Optical Properties [table:slatopticalproperties]
τdir,dif 
Directtodiffuse transmittance (same for front and back of slat) 
τdif,dif 
Diffusetodiffuse transmittance (same for front and back of slat) 
ρfdir,dif, ρbdir,dif

Front and back directtodiffuse reflectance 
ρfdif,dif, ρbdif,dif

Front and back diffusetodiffuse reflectance 
It is assumed that there is no directtodirect transmission or reflection, so that τdir,dir=0, ρfdir,dir=0, and ρbdir,dir=0. It is further assumed that the slats are perfect diffusers, so that τdir,dif, ρfdir,dif and ρbdir,dif are independent of angle of incidence. Until the EnergyPlus model is improved to take into account the angleofincidence dependence of slat transmission and reflection, it is assumed that τdir,dif = τdif,dif, ρfdir,dif = ρfdif,dif, and ρbdir,dif = ρbdif,dif.
Direct Transmittance of Blind[LINK]
The directtodirect and directtodiffuse transmittance of a blind is calculated using the slat geometry shown in Figure 7(a), which shows the side view of one of the cells of the blind. For the case shown, each slat is divided into two segments, so that the cell is bounded by a total of six segments, denoted by s1 through s6 (note in the following that si refers to both segment i and the length of segment i).The lengths of s1 and s2 are equal to the slat separation, h, which is the distance between adjacent slat faces. s3 and s4 are the segments illuminated by direct radiation. In the case shown in Figure 7(a) the cell receives radiation by reflection of the direct radiation incident on s4 and, if the slats have nonzero transmittance, by transmission through s3, which is illuminated from above.
The goal of the blind direct transmission calculation is to determine the direct and diffuse radiation leaving the cell through s2 for unit direct radiation entering the cell through s1.
DirecttoDirect Blind Transmittance[LINK]
Figure 7(b) shows the case where some of the direct radiation passes through the cell without hitting the slats. From the geometry in this figure we see that
τdir,dirbl,f=1−wh,w≤h
where
w=scos(φb−φs)cosφs
Note that we are assuming that the slat thickness is zero. A correction for nonzero slat thickness is described later.
DirecttoDiffuse Blind Transmittance, Reflectance and Absorptance[LINK]
The directtodiffuse and transmittance and reflectance of the blind are calculated using a radiosity method that involves the following three vector quantities:
Ji = the radiosity of segment si, i.e., the total radiant flux into the cell from si
Gi = the irradiance on the cell side of si
Qi = the source flux from the cell side of si
Based on these definitions we have the following equations that relate J, G and Q for the different segments:
J1=Q1J2=Q2J3=Q3+ρbdif,difG3+τdif,difG4J4=Q4+τdif,difG3+ρfdif,difG4J5=Q5+ρbdif,difG5+τdif,difG6J6=Q6+τdif,difG5+ρfdif,difG6
In addition we have the following equation relating G and J:
Gi=6∑j=1JjFji, i=1,6
where Fji is the view factor between sj and si, i.e., Fji is the fraction of radiation leaving sj that is intercepted by si.
Using J1=Q1=0 and J2=Q2=0 and combining the above equations gives the following equation set relating J and Q:
J3−ρbdif,dif6∑j=3JjFj3−τdif,dif6∑j=3JjFj4=Q3
J4−τdif,dif6∑j=3JjFj3−ρfdif,dif6∑j=3JjFj4=Q4
J5−ρbdif,dif6∑j=3JjFj5−τdif,dif6∑j=3JjFj6=Q5
J6−τdif,dif6∑j=3JjFj3−ρfdif,dif6∑j=3JjFj6=Q6
This can be written in the form:
Q′=XJ′
where X is a 4x4 matrix and
J′=⎡⎢
⎢
⎢⎣J3J4J5J6⎤⎥
⎥
⎥⎦
Q′=⎡⎢
⎢
⎢
⎢⎣Q3Q4Q5Q6⎤⎥
⎥
⎥
⎥⎦
We then obtain J′ from:
J′=X−1Q′
The view factors, Fij, are obtained as follows. The cell we are dealing with is a convex polygon with n sides. In such a polygon the view factors must satisfy the following constraints:
n∑j=1Fij=1, i=1,n
siFij=sjFji,i=1,n; j=1,n
Fii=0, i=1,n
These constraints lead to simple equations for the view factors for n = 3 and 4. For n = 3, we have the following geometry and view factor expression:
For n = 4 we have:
Applying these to the slat cell shown in Figure 8 we have the following:
F12=d1+d2−2s2h
F13=h+s3−d32h,etc.
The sources for the directtodiffuse transmittance calculation are:
Q1=Q2=Q5=Q6=0(and therefore J1=J2=0)
Q3=τdir,difQ4=ρfdir,dif}φb≤φs+π2(beam hits front of slats)
Q3=ρbdir,difQ4=τdir,dif}φb>φs+π2(beam hits back of slats)
For unit incident direct flux, the front directtodiffuse transmittance and reflectance of the blind are:
τdir,difbl,f=G2ρdir,difbl,f=G1
where
G2=∑6j=3JjFj2G1=∑6j=3JjFj1
and J3 to J6 are given by Equation [eq:QequalsXJprime].
The front direct absorptance of the blind is then:
αdirbl,f=1−τdir,difbl,f−τdir,dirbl,f−ρdir,difbl,f
The directtodiffuse calculations are performed separately for solar and visible slat properties to get the corresponding solar and visible blind properties.
Dependence on Profile Angle[LINK]
The directtodirect and directtodiffuse blind properties are calculated for direct radiation profile angles (see Figure 7) ranging from –90o to +90o in 5o increments. (The “profile angle” is the angle of incidence in a plane that is perpendicular to the window and perpendicular to the slat direction.) In the time step loop the blind properties for a particular profile angle are obtained by interpolation.
Dependence on Slat Angle[LINK]
All blind properties are calculated for slat angles ranging from –90o to +90o in 10o increments. In the timestep loop the slat angle is determined by the slatangle control mechanism and then the blind properties at that slat angle are determined by interpolation. Three slatangle controls are available: (1) slat angle is adjusted to just block beam solar incident on the window; (2) slat angle is determined by a schedule; and (3) slat angle is fixed.
DiffusetoDiffuse Transmittance and Reflectance of Blind[LINK]
To calculate the diffusetodiffuse properties, assuming uniformly distributed incident diffuse radiation, each slat bounding the cell is divided into two segments of equal length (Figure 9), i.e., s3=s4 and s5=s6. For frontside properties we have a unit source, Q1=1. All the other Qi are zero. Using this source value, we apply the methodology described above to obtain G2 and G1. We then have:
τdif,difbl,f=G2ρdif,difbl,f=G1αdifbl,f=1−τdif,difbl,f−ρdif,difbl,f
The backside properties are calculated in a similar way by setting Q2 = 1 with the other Qi equal to zero.
The diffusetodiffuse calculations are performed separately for solar, visible and IR slat properties to get the corresponding solar, visible and IR blind properties.
Blind properties for sky and ground diffuse radiation[LINK]
For horizontal slats on a vertical window (the most common configuration) the blind diffusetodiffuse properties will be sensitve to whether the radiation is incident upward from the ground or downward from the sky (Figure 10). For this reason we also calculate the following solar properties for a blind consisting of horizontal slats in a vertical plane:
τgnd−dif,difbl,f= front transmittance for ground diffuse solar
τsky−dif,difbl,f= front transmittance for sky diffuse solar
ρgnd−dif,difbl,f= front reflectance for ground diffuse solar
ρsky−dif,difbl,f= front reflectance for sky diffuse solar
αgnd−dif,difbl,f= front absorptance for ground diffuse solar
αsky−dif,difbl,f= front absorptance for sky diffuse solar
These are obtained by integrating over sky and ground elements, as shown in Figure 10, treating each element as a source of direct radiation of irradiance I(ϕs) incident on the blind at profile angle ϕs. This gives:
τsky−dif,difbl,f=π/2∫0[τdir,dirbl,f(ϕs)+τdir,difbl,f(ϕs)]Isky(ϕs)cosϕsdϕsπ/2∫0Isky(ϕs)cosϕsdϕs
ρsky−dif,difbl,f=π/2∫0ρdir,difbl,fIsky(ϕs)cosϕsdϕsπ/2∫0Isky(ϕs)cosϕsdϕs
αsky−difbl,f=π/2∫0αdirbl,fIsky(ϕs)cosϕsdϕsπ/2∫0Isky(ϕs)cosϕsdϕs
We assume that the sky radiance is uniform. This means that Isky is independent of ϕs, giving:
τsky−dif,difbl,f=π/2∫0[τdir,dirbl,f+τdir,difbl,f]cosϕsdϕs
ρsky−dif,difbl,f=π/2∫0ρdir,difbl,fcosϕsdϕs
αsky−difbl,f=π/2∫0αdirbl,fcosϕsdϕs
The corresponding ground diffuse quantities are obtained by integrating ϕs from −π/2 to 0.
An improvement to this calculation would be to allow the sky radiance distribution to be nonuniform, i.e., to depend on sun position and sky conditions, as is done in the detailed daylighting calculation (see “Sky Luminance Distributions” under “Daylight Factor Calculation”).
Correction Factor for Slat Thickness[LINK]
A correction has to be made to the blind transmittance, reflectance and absorptance properties to account for the amount of radiation incident on a blind that is reflected and absorbed by the slat edges (the slats are assumed to be opaque to radiation striking the slat edges). This is illustrated in Figure 11 for the case of direct radiation incident on the blind. The slat crosssection is assumed to be rectangular. The quantity of interest is the fraction, fedge, of direct radiation incident on the blind that strikes the slat edges. Based on the geometry shown in Figure 11 we see that
fedge=tcosγ(h+tcosξ)cosφs=tcos(φs−ξ)(h+tcosξ)cosφs=tsin(φb−φs)(h+tsinφb)cosφs
The edge correction factor for diffuse incident radiation is calculated by averaging this value of fedge over profile angles, φs, from 90o to +90o.
As an example of how the edge correction factor is applied, the following two equations show how blind front diffuse transmittance and reflectance calculated assuming zero slat thickness are modified by the edge correction factor. It is assumed that the edge transmittance is zero and that the edge reflectance is the same as the slat front reflectance, ρf.
τdif,difbl,f→τdif,difbl,f(1−fedge)ρdifbl,f→ρdifbl,f(1−fedge)+fedgeρf
Window Calculation Module[LINK]
This section describes two potential modeling approaches for Windows. The first (layer by layer) is implemented. The second, simple approach, reuses the layerbylayer approach but converts an arbitrary window performance into an equivalent single layer.
The primary Window calculation is a layerbylayer approach where windows are considered to be composed of the following components, only the first of which, glazing, is required to be present:
Glazing, which consists of one or more plane/parallel glass layers. If there are two or more glass layers, the layers are separated by gaps filled with air or another gas. The glazing optical and thermal calculations are based on algorithms from the WINDOW 4 and WINDOW 5 programs [Arasteh et al., 1989], [Finlayson et al., 1993]. Glazing layers are described using te input object WindowMaterial:Glazing.
Gap, layers filled with air or another gas that separate glazing layers. Gaps are described using the input object WindowMaterial:Gas.
Frame, which surrounds the glazing on four sides. Frames are described using the input object WindowProperty:FrameAndDivider.
Divider, which consists of horizontal and/or vertical elements that divide the glazing into individual lites.
Shading device, which is a separate layer, such as drapery, roller shade or blind, on the inside or outside of the glazing, whose purpose is to reduce solar gain, reduce heat loss (movable insulation) or control daylight glare. Shading layers are described using “WindowShadingControl” input objects.
In the following, the description of the layerbylayer glazing algorithms is based on material from Finlayson et al., 1993. The frame and divider thermal model, and the shading device optical and thermal models, are new to EnergyPlus.
A second approch has been developed where windows are modeled in a simplified approach that requires minimal user input that is processed to develop and equivalent layer that then reuses much of the layerbymodel. This “Simple Window Construction: model is described below.
Optical Properties of Glazing[LINK]
The solar radiation transmitted by a system of glass layers and the solar radiation absorbed in each layer depends on the solar transmittance, reflectance and absorptance properties of the individual layers. The absorbed solar radiation enters the glazing heat balance calculation that determines the inside surface temperature and, therefore, the heat gain to the zone from the glazing (see “Window Heat Balance Calculation“). The transmitted solar radiation is absorbed by interior zone surfaces and, therefore, contributes to the zone heat balance. In addition, the visible transmittance of the glazing is an important factor in the calculation of interior daylight illuminance from the glazing.
Glass Layer Properties[LINK]
In EnergyPlus, the optical properties of individual glass layers are given by the following quantities at normal incidence as a function of wavelength:
Transmittance, T
Front reflectance, Rf
Back reflectance, Rb
Here “front” refers to radiation incident on the side of the glass closest to the outside environment, and “back” refers to radiant incident on the side of the glass closest to the inside environment. For glazing in exterior walls, “front” is therefore the side closest to the outside air and “back” is the side closest to the zone air. For glazing in interior (i.e., interzone) walls, “back” is the side closest to the zone in which the wall is defined in and “front” is the side closest to the adjacent zone.
Glass Optical Properties Conversion[LINK]
Conversion from Glass Optical Properties Specified as Index of Refraction and Transmittance at Normal Incidence[LINK]
The optical properties of uncoated glass are sometimes specified by index of refraction, n,* * and transmittance at normal incidence, T.
The following equations show how to convert from this set of values to the transmittance and reflectance values required by Material:WindowGlass. These equations apply only to uncoated glass, and can be used to convert either spectralaverage solar properties or spectralaverage visible properties (in general, n and T are different for the solar and visible). Note that since the glass is uncoated, the front and back reflectances are the same and equal to the R that is solved for in the following equations.
Given n and T, find R:
r=(n−1n+1)2τ=[(1−r)4+4r2T2]1/2−(1−r)22r2TR=r+(1−r)2rτ21−r2τ2
Example:
T=0.86156n=1.526r=(1.526−11.526+1)2τ=0.93974R=0.07846
Simple Window Model[LINK]
EnergyPlus includes an alternate model that allows users to enter in simplified window performance indices. This model is accessed through the WindowMaterial:SimpleGlazingSystem input object and converts the simple indices into an equivalent single layer window. (In addition a special model is used to determine the angular properties of the system – described below). Once the model generates the properties for the layer, the program reuses the bulk of the layerbylayer model for subsequent calculations. The properties of the equivalent layer are determined using the step by step method outlined by Arasteh, Kohler, and Griffith (2009). The core equations are documented here. The reference contains additional information.
The simplified window model accepts U and SHGC indices and is useful for several reasons:
1) Sometimes, the only thing that is known about the window are its U and SHGC;
2) Codes, standards, and voluntary programs are developed in these terms;
3) A singlelayer calculation is faster than multilayer calculations.
Note: This use of U and SHGC to describe the thermal properties of windows is only appropriate for specular glazings.
While it is important to include the ability to model windows with only Uvalue and SHGC, we note that any method to use U and SHGC alone in building simulation software will inherently be approximate. This is due primarily to the following factors:
SHGC combines directly transmitted solar radiation and radiation absorbed by the glass which flows inward. These have different implications for space heating/cooling. Different windows with the same SHGC often have different ratios of transmitted to absorbed solar radiation.
SHGC is determined at normal incidence; angular properties of glazings vary with number of layers, tints, coatings. So products which have the same SHGC, can have different angular properties.
Window Ufactors vary with temperatures.
Thus, for modeling specific windows, we recommend using more detailed data than just the U and SHGC, if at all possible.
The simplified window model determines the properties of an equivalent layer in the following steps.
Step 1. Determine glasstoglass Resistance.[LINK]
Window Uvalues include interior and exterior surface heat transfer coefficients. The resistance of the bare window product, or glasstoglass resistance is augmented by these film coefficients so that,
1U=Ri,w+Ro,w+Rl,w
Where,
Ri,w is the resistance of the interior film coefficient under standard winter conditions in units of m2·K/W,
Ro,w is the resistance of the exterior film coefficient under standard winter conditions in units of m2·K/W, and
Rl,w is the resisance of the bare window under winter conditions (without the film coefficients) in units of m2·K/W.
The values for Ri,w and Ro,w depend on U and are calculated using the following correlations.
Ri,w=1(0.359073Ln(U)+6.949915);forU<5.85
Ri,w=1(1.788041U−2.886625);forU≥5.85
Ro,w=1(0.025342U+29.163853)
So that the glasstoglass resistance is calculated using:
Rl,w=1U−Ri,w−Ro,w
Because the window model in EnergyPlus is for flat geometries, the models are not necessarily applicable to lowperformance projecting products, such as skylights with unisulated curbs. The model cannot support glazing systems with a U higher than 7.0 because the thermal resistance of the film coefficients alone can provide this level of performance and none of the various resistances can be negative.
Step 2. Determine Layer Thickness.[LINK]
The thickness of the equivalent layer in units of meters is calculated using,
Thickness=⎧⎪⎨⎪⎩0.002for 1Rl,w>7.00.05914−0.00714Rl,wfor 1Rl,w≤7.0
Step 3. Determine Layer Thermal Conductivity[LINK]
The effective thermal conductivity, λeff, of the equivalent layer is calculated using,
λeff=ThicknessRl,w
Step 4. Determine Layer Solar Transmittance[LINK]
The layer’s solar transmittance at normal incidence, Tsol, is calculated using correlations that are a function of SHGC and UFactor.
Tsol=0.939998SHGC2+0.20332SHGC;U>4.5;SHGC<0.7206
Tsol=1.30415SHGC−0.30515;U>4.5;SHGC≥0.7206
Tsol=0.41040SHGC;U<3.4;SHGC≤0.15
Tsol=0.085775SHGC2+0.963954SHGC−0.084958;U<3.4;SHGC>0.15
And for Uvalues between 3.4 and 4.5, the value for Tsol is interpolated using results of the equations for both ranges.
Step 5. Determine Layer Solar Reflectance[LINK]
The layer’s solar reflectance is calculated by first determining the inward flowing fraction which requires values for the resistance of the inside and outside film coefficients under summer conditions, Ri,s and Ro,s, respectively. The correlations are:
Ri,s=1(29.436546(SHGC−TSol)3−21.943415(SHGC−TSol)2+9.945872(SHGC−TSol)+7.426151);U>4.5Ri,s=1(199.8208128(SHGC−TSol)3−90.639733(SHGC−TSol)2+19.737055(SHGC−TSol)+6.766575);U<3.4Ro,s=1(2.225824(SHGC−TSol)+20.57708);U>4.5Ro,s=1(5.763355(SHGC−TSol)+20.541528);U<3.4
And for Uvalues between 3.4 and 4.5, the values are interpolated using results from both sets of equations.
The inward flowing fraction, Fracinward, is then calculated using:
Fracinward=(Ro,s+0.5Rl,w)(Ro,s+Rl,w+Ri,s)
Then, the solar reflectances of the front face, Rs,f, and back face, Rs,b, are calculated using:
Rs,f=Rs,b=1−TSol−(SHGC−TSol)Fracinward
The thermal absorptance, or emittance, is taken as 0.84 for both the front and back and the longwave transmittance is 0.0.
Step 6. Determine Layer Visible Properties[LINK]
The user has the option of entering a value for visible transmittance as one of the simple performance indices. If the user does not enter a value, then the visible properties are the same as the solar properties. If the user does enter a value then layer’s visible transmittance at normal incidence, TVis, is set to that value. The visible light reflectance for the back surface is calculated using:
RVis,b=−0.7409T3Vis+1.6531T2Vis−1.2299TVis+0.4547
The visible light reflectance for the front surface is calculated using:
RVis,f=−0.0622T3Vis+0.4277T2Vis−0.4169TVis+0.2399
Step 7. Determine Angular Performance[LINK]
The angular properties of windows are important because during energy modeling, the solar incidence angles are usually fairly high. Angles of incidence are defined as angles from the normal direction extending out from the window. The simple glazing system model includes a range of correlations that are selected based on the values for U and SHGC. These were chosen to match the types of windows likely to have such performance levels. The matrix of possible combinations of U and SHGC values have been mapped to set of 28 bins shown in the following figure.
Diagram of Transmittance and Reflectance Correlations Used based on U and SHGC [fig:diagramoftransmittanceand]
There are ten different correlations, A thru J, for both transmission and reflectance. The correlations are used in various weighting and interpolation schemes according the figure above. The correlations are normalized against the performance at normal incidence. EnergyPlus uses these correlations to store the glazing system’s angular performance at 10 degree increments and interpolates between them during simulations. The model equations use the cosine of the incidence angle, cos(φ), as the independent variable. The correlations have the form:
TorR=acos(φ)4+bcos(φ)3+ccos(φ)2+dcos(φ)+e
The coefficient values for a, b, c, d, and e are listed in the following tables for each of the curves.
Normalized Transmittance Correlations for Angular Performance [fig:normalizedtransmittancecorrelationsfor]
Normalized Reflectanct Correlations for Angular Performance [fig:normalizedreflectanctcorrelationsfor]
Application Issues[LINK]
EnergyPlus’s normal process of running the detailed layerbylayer model, with the equivalent layer produced by this model, creates reports (sent to the EIO file) of the overall performance indices and the properties of the equivalent layer. Both of these raise issues that may be confusing.
The simplified window model does not reuse all aspects of the detailed layerbylayer model, in that the angular solar transmission properties use a different model when the simple window model is in effect. If the user takes the material properties of an equivalent glazing layer from the simple window model and then reenters them into just the detailed model, then the performance will not be the same because of the angular transmission model will have changed. It is not proper use of the model to reenter the equivalent layer’s properties and expect the exact level of performance.
There may not be exact agreement between the performance indices echoed out and those input in the model. This is expected with the model and the result of a number of factors. For example, although input is allowed to go up to U7 W/m2⋅K, the actual outcome is limited to no higher than about 5.8W/m2⋅K. This is because the thermal resistance to heat transfer at the surfaces is already enough resistance to provide an upper limit to the conductance of a planar surface. Sometimes there is conflict between the SHGC and the U that are not physical and compromises need to be made. In general, the simple window model is intended to generate a physicallyreasonable glazing that approximates the input entered as well as possible. But the model is not always be able to do exactly what is specified when the specifications are not physical.
References[LINK]
Arasteh, D., J.C. Kohler, B. Griffith, Modeling Windows in EnergyPlus with Simple Performance Indices. Lawrence Berkeley National Laboratory. In Draft. Available at
Glazing System Properties[LINK]
The optical properties of a glazing system consisting of N glass layers separated by nonabsorbing gas layers (see Figure 4) are determined by solving the following recursion relations for Ti,j, the transmittance through layers i to j; Rfi,j and Rbi,j, the front and back reflectance, respectively, from layers i to j; and Aj, the absorption in layer j. Here layer 1 is the outermost layer and layer N is the innermost layer. These relations account for multiple internal reflections within the glazing system. Each of the variables is a function of wavelength.
Ti,j=Ti,j−1Tj,j1−Rfj,jRbj−1,i
Rfi,j=Rfi,j−1+T2i,j−1Rfj,j1−Rfj,jRbj−1,i
Rbj,i=Rbj,j+T2j,jRbj−1,i1−Rbj−1,iRfj,j
Afj=T1,j−1(1−Tj,j−Rfj,j)1−Rfj,NRbj−1,1+T1,jRfj+1,N(1−Tj,j−Rbj,j)1−Rfj,NRbj−1,1
In Equation [eq:Ajtothefequation], Ti,j = 1 and Ri,j = 0 if i<0 or j>N.
Schematic of transmission, reflection and absorption of solar radiation within a multilayer glazing system. [fig:schematicoftransmissionreflection]
As an example, for double glazing (N=2), these equations reduce to:
T1,2=T1,1T2,21−Rf2,2Rb1,1
Rf1,2=Rf1,1+T21,1Rf2,21−Rf2,2Rb1,1
Rb2,1=Rb2,2+T22,2Rb1,11−Rb1,1Rf2,2
Af1=(1−T1,1−Rf1,1)+T1,1Rf2,2(1−T1,1−Rb1,1)1−Rf2,2Rb1,1
Af2=T1,1(1−T2,2−Rf2,2)1−Rf2,2Rb1,1
If the above transmittance and reflectance properties are input as a function of wavelength, EnergyPlus calculates “spectral average” values of the above glazing system properties by integrating over wavelength.
The spectralaverage solar property is:
Ps=∫P(λ)Es(λ)dλ∫Es(λ)dλ
The spectralaverage visible property is:
Pv=∫P(λ)Es(λ)V(λ)dλ∫Es(λ)V(λ)dλ
where Es(λ) is the solar spectral irradiance function and V(λ) is the photopic response function of the eye. The default functions are shown in Table 2 and Table 3. They can be overwritten by user defined solar and/or visible spectrum using the objects Site:SolarAndVisibleSpectrum and Site:SpectrumData. They are expressed as a set of values followed by the corresponding wavelengths for values.
When a choice of Spectral is entered as the optical data type, the correlations to store the glazing system’s angular performance are generated based on angular performance at 10 degree increments. When a choice of SpectralAndAngle is entered as the optical data type, the correlations for the glazing system will be generated using 10 degree increments or more if the SpectralAndAngle properties include data for more angles. For each incident angle, the properties of the SpectralAndAngle layer(s) is calculated by linear interpolation, and then the performance of the entire glazing system is calculated for that angle. The glazing system properties at each angle are used to generate polynomial curve fits with 6 coefficients as a function of cosine of incident angle. The polynomial curves are then used in the simulation to calculate optical properties at each timestep.
If a glazing layer has optical properties that are roughly constant with wavelength, the wavelengthdependent values of Ti,i, Rfi,i and Rbi,i in Equations [eq:Tijequation] to [eq:Ajtothefequation] can be replaced with constant values for that layer.
: Solar spectral irradiance function.
[table:solarspectralirradiancefunction.]
: Photopic response function.
[table:photopicresponsefunction.]
Calculation of Angular Properties[LINK]
Calculation of optical properties is divided into two categories: uncoated glass and coated glass.
Angular Properties for Uncoated Glass[LINK]
The following discussion assumes that optical quantities such as transmissivity, reflectvity, absorptivity, and index of refraction are a function of wavelength, λ. If there are no spectral data the angular dependence is calculated based on the single values for transmittance and reflectance in the visible and solar range. In the visible range an average wavelength of 0.575 microns is used in the calculations. In the solar range an average wavelength of 0.898 microns is used.
The spectral data include the transmittance, T, and the reflectance, R. For uncoated glass the reflectance is the same for the front and back surfaces. For angle of incidence, ϕ , the transmittance and reflectance are related to the transmissivity, τ, and reflectivity, ρ, by the following relationships:
T(ϕ)=τ(ϕ)2e−αd/cosϕ′1−ρ(ϕ)2e−2αd/cosϕ′
R(ϕ)=ρ(ϕ)(1+T(ϕ)e−αd/cosϕ′)
The spectral reflectivity is calculated from Fresnel’s equation assuming unpolarized incident radiation:
ρ(ϕ)=12((ncosϕ−cosϕ′ncosϕ+cosϕ′)2+(ncosϕ′−cosϕncosϕ′+cosϕ)2)
The spectral transmittivity is given by:
τ(ϕ)=1−ρ(ϕ)
The spectral absorption coefficient is defined as:
α=4πκλ
where κ is the dimensionless spectrallydependent extinction coefficient and λ is the wavelength expressed in the same units as the sample thickness.
Solving Equation [eq:RhoofPhi] at normal incidence gives:
n=1+√ρ(0)1−√ρ(0)
Evaluating Equation [eq:RofPhiEquation] at normal incidence gives the following expression for κ:
κ=−λ4πdlnR(0)−ρ(0)ρ(0)T(0)
Eliminating the exponential in Equations [eq:TofPhiEquation] and [eq:RofPhiEquation] gives the reflectivity at normal incidence:
ρ(0)=β−√β2−4(2−R(0))R(0)2(2−R(0))
where
β=T(0)2−R(0)2+2R(0)+1
The value for the reflectivity, ρ(0), from Equation [eq:Rho0Equation] is substituted into Equations [eq:NAsFunctionOfRhoEquation] and [eq:KappaAsFunctionOfLambdaRTEquation]. The result from Equation [eq:KappaAsFunctionOfLambdaRTEquation] is used to calculate the absorption coefficient in Equation [eq:AlphaAsFunctionOfKappaLambda]. The index of refraction is used to calculate the reflectivity in Equation [eq:RhoofPhi] which is then used to calculate the transmittivity in Equation [eq:TauofPhiEquation]. The reflectivity, transmissivity and absorption coefficient are then substituted into Equations [eq:TofPhiEquation] and [eq:RofPhiEquation] to obtain the angular values of the reflectance and transmittance.
Angular Properties for Coated Glass[LINK]
A regression fit is used to calculate the angular properties of coated glass from properties at normal incidence. If the transmittance of the coated glass is > 0.645, the angular dependence of uncoated clear glass is used. If the transmittance of the coated glass is ≤ 0.645, the angular dependence of uncoated bronze glass is used. The values for the angular functions for the transmittance and reflectance of both clear glass (¯τclr,¯ρclr) and bronze glass (¯τbnz,¯ρbnz) are determined from a fourthorder polynomial regression:
¯τ(ϕ)=¯τ0+¯τ1cos(ϕ)+¯τ2cos2(ϕ)+¯τ3cos3(ϕ)+¯τ4cos4(ϕ)
and
¯ρ(ϕ)=¯ρ0+¯ρ1cos(ϕ)+¯ρ2cos2(ϕ)+¯ρ3cos3(ϕ)+¯ρ4cos4(ϕ)−¯τ(ϕ)
The polynomial coefficients are given in Table 4.
These factors are used as follows to calculate the angular transmittance and reflectance:
For T(0) > 0.645:
T(ϕ)=T(0)¯τclr(ϕ)
R(ϕ)=R(0)(1−¯ρclr(ϕ))+¯ρclr(ϕ)
For T(0) ≤ 0.645:
T(ϕ)=T(0)¯τbnz(ϕ)
R(ϕ)=R(0)(1−¯ρbnz(ϕ))+¯ρbnz(ϕ)
Angular Properties for Simple Glazing Systems[LINK]
When the glazing system is modeled using the simplified method, an alternate method is used to determine the angular properties. The equation for solar transmittance as a function of incidence angle, T(ϕ), is:
T(ϕ)=T(ϕ=0)cos(ϕ)(1+(0.768+0.817SHGC4)sin3(ϕ))
where,
T(ϕ=0) is the normal incidence solar transmittance, TSol.
The equation for solar reflectance as a function of incidence angle, R(ϕ), is:
R(ϕ)=R(ϕ=0)(f1(ϕ)+f2(ϕ)√SHGC)Rfit,o
where,
f1(ϕ)=(((2.403cos(ϕ)−6.192)cos(ϕ)+5.625)cos(ϕ)−2.095)cos(ϕ)+1
f2(ϕ)=(((−1.188cos(ϕ)+2.022)cos(ϕ)+0.137)cos(ϕ)−1.720)cos(ϕ)
Rfit,o=0.7413−(0.7396√SHGC)
Calculation of Hemispherical Values[LINK]
The hemispherical value of a property is determined from the following integral:
Phemispherical=2∫π20P(ϕ)cos(ϕ)sin(ϕ)dϕ
The integral is evaluated by Simpson’s rule for property values at angles of incidence from 0 to 90 degrees in 10degree increments.
Optical Properties of Window Shading Devices[LINK]
Shading devices affect the system transmittance and glass layer absorptance for shortwave radiation and for longwave (thermal) radiation. The effect depends on the shade position (interior, exterior or betweenglass), its transmittance, and the amount of interreflection between the shading device and the glazing. Also of interest is the amount of radiation absorbed by the shading device.
In EnergyPlus, shading devices are divided into four categories, “shades,” “blinds,” “screens,” and “switchable glazing.” “Shades” are assumed to be perfect diffusers. This means that direct radiation incident on the shade is reflected and transmitted as hemispherically uniform diffuse radiation: there is no direct component of transmitted radiation. It is also assumed that the transmittance, τsh, reflectance, ρsh, and absorptance, αsh, are the same for the front and back of the shade and are independent of angle of incidence. Many types of drapery and pulldown roller devices are close to being perfect diffusers and can be categorized as “shades.”
“Blinds” in EnergyPlus are slattype devices such as venetian blinds. Unlike shades, the optical properties of blinds are strongly dependent on angle of incidence. Also, depending on slat angle and the profile angle of incident direct radiation, some of the direct radiation may pass between the slats, giving a direct component of transmitted radiation.
“Screens” are debris or insect protection devices made up of metallic or nonmetallic materials. Screens may also be used as shading devices for large glazing areas where excessive solar gain is an issue. The EnergyPlus window screen model assumes the screen is composed of intersecting orthogonallycrossed cylinders, with the surface of the cylinders assumed to be diffusely reflecting. Screens may only be used on the exterior surface of a window construction. As with blinds, the optical properties affecting the direct component of transmitted radiation are dependent on the angle of incident direct radiation.
With “Switchable glazing,” shading is achieved making the glazing more absorbing or more reflecting, usually by an electrical or chemical mechanism. An example is electrochromic glazing where the application of an electrical voltage or current causes the glazing to switch from light to dark.
Shades and blinds can be either fixed or moveable. If moveable, they can be deployed according to a schedule or according to a trigger variable, such as solar radiation incident on the window. Screens can be either fixed or moveable according to a schedule.
Shades[LINK]
Shade/Glazing System Properties for ShortWave Radiation[LINK]
Shortwave radiation includes:
Beam solar radiation from the sun and diffuse solar radiation from the sky and ground incident on the outside of the window,
Beam and/or diffuse radiation reflected from exterior obstructions or the building itself,
Solar radiation reflected from the inside zone surfaces and incident as diffuse radiation on the inside of the window,
Beam solar radiation from one exterior window incident on the inside of another window in the same zone, and
Shortwave radiation from electric lights incident as diffuse radiation on the inside of the window.
Exterior Shade
For an exterior shade we have the following expressions for the system transmittance, the effective system glass layer absorptance, and the system shade absorptance, taking interreflection between shade and glazing into account. Here, “system” refers to the combination of glazing and shade. The system properties are given in terms of the isolated shade properties (i.e., shade properties in the absence of the glazing) and the isolated glazing properties (i.e., glazing properties in the absence of the shade).
Tsys(ϕ)=Tdif1,Nτsh1−Rdiffρsh
Tdifsys=Tdif1,Nτsh1−Rdiffρsh
Asysj,f(ϕ)=Adifj,fτsh1−Rfρsh,j=1 to N
Adif,sysj,f=Adifj,fτsh1−Rfρsh,j=1 to N
Adif,sysj,b=Adifj,bTdif1,Nρsh1−Rfρsh,j=1 to N
αsyssh=αsh(1+τshRf1−Rfρsh)
Interior Shade
The system properties when an interior shade is in place are the following:
Tsys(ϕ)=T1,N(ϕ)τsh1−Rdifbρsh
Tdifsys=Tdif1,Nτsh1−Rdifbρsh
Asysj,f(ϕ)=Aj,f(ϕ)+T1,N(ϕ)ρsh1−RdifbρshAdifj,b,j=1 to N
Adif,sysj,f=Adifj,f+Tdif1,Nρsh1−RdifbρshAdifj,b,j=1 to N
Adif,sysj,b=τsh1−RdifbρshAdifj,b, j=1 to N
αsyssh(ϕ)=T1,N(ϕ)αsh1−Rdifbρsh
αdif,syssh=Tdif1,Nαsh1−Rdifbρsh
LongWave Radiation Properties of Window Shades[LINK]
Longwave radiation includes:
Thermal radiation from the sky, ground and exterior obstructions incident on the outside of the window,
Thermal radiation from other room surfaces incident on the inside of the window, and
Thermal radiation from internal sources, such as equipment and electric lights, incident on the inside of the window.
The program calculates how much longwave radiation is absorbed by the shade and by the adjacent glass surface. The system emissivity (thermal absorptance) for an interior or exterior shade, taking into account reflection of longwave radiation between the glass and shade, is given by:
εlw,syssh=εlwsh⎛⎝1+τlwshρlwgl1−ρlwshρlwgl⎞⎠
where ρlwgl is the longwave reflectance of the outermost glass surface for an exterior shade or the innermost glass surface for an interior shade, and it is assumed that the longwave transmittance of the glass is zero.
The innermost (for interior shade) or outermost (for exterior shade) glass surface emissivity when the shade is present is:
εlw,sysgl=εlwglτlwsh1−ρlwshρlwgl
Switchable Glazing[LINK]
For switchable glazing, such as electrochromics, the solar and visible optical properties of the glazing can switch from a light state to a dark state. The switching factor, fswitch, determines what state the glazing is in. An optical property, p, such as transmittance or glass layer absorptance, for this state is given by:
p=(1−fswitch)plight+fswitchpdark
where
plight is the property value for the unswitched, or light state, and pdark is the property value for the fully switched, or dark state.
The value of the switching factor in a particular time step depends on what type of switching control has been specified: “schedule,” “trigger,” or “daylighting.” If “schedule,” fswitch = schedule value, which can be 0 or 1.
Thermochromic Windows[LINK]
Thermochromic (TC) materials have active, reversible optical properties that vary with temperature. Thermochromic windows are adaptive window systems for incorporation into building envelopes. Thermochromic windows respond by absorbing sunlight and turning the sunlight energy into heat. As the thermochromic film warms it changes its light transmission level from less absorbing to more absorbing. The more sunlight it absorbs the lower the light level going through it. Figure 5 shows the variations of window properties with the temperature of the thermochromic glazing layer. By using the suns own energy the window adapts based solely on the directness and amount of sunlight. Thermochromic materials will normally reduce optical transparency by absorption and/or reflection, and are specular (maintaining vision).
Variations of Window Properties with the Temperature of the Thermochromic Glazing Layer [fig:variationsofwindowpropertieswith]
On cloudy days the window is at full transmission and letting in diffuse daylighting. On sunny days the window maximizes diffuse daylighting and tints based on the angle of the sun relative to the window. For a south facing window (northern hemisphere) the daylight early and late in the day is maximized and the direct sun at mid day is minimized.
The active thermochromic material can be embodied within a laminate layer or a surface film. The overall optical state of the window at a given time is a function primarily of:
thermochromic material properties
solar energy incident on the window
construction of the window system that incorporates the thermochromic layer
environmental conditions (interior, exterior, air temperature, wind, etc).
The tinted film, in combination with a heat reflecting, lowe layer allows the window to reject most of the absorbed radiation thus reducing undesirable heat load in a building. In the absence of direct sunlight the window cools and clears and again allows lower intensity diffuse radiation into a building. TC windows can be designed in several ways (Figure 6), with the most common being a triple pane windows with the TC glass layer in the middle a double pane windows with the TC layer on the inner surface of the outer pane or for sloped glazing a double pane with the laminate layer on the inner pane with a lowe layer toward the interior. The TC glass layer has variable optical properties depending on its temperature, with a lower temperature at which the optical change is initiated, and an upper temperature at which a minimum transmittance is reached. TC windows act as passive solar shading devices without the need for sensors, controls and power supplies but their optical performance is dependent on varying solar and other environmental conditions at the location of the window.
Configurations of Thermochromic Windows [fig:configurationsofthermochromicwindows]
EnergyPlus describes a thermochromic window with a Construction object which references a special layer defined with a WindowMaterial:GlazingGroup:Thermochromic object. The WindowMaterial:GlazingGroup:Thermochromic object further references a series of WindowMaterial:Glazing objects corresponding to each specification temperature of the TC layer. During EnergyPlus run time, a series of TC windows corresponding to each specification temperature is created once. At the beginning of a particular time step calculations, the temperature of the TC glass layer from the previous time step is used to look up the most closed specification temperature whose corresponding TC window construction will be used for the current time step calculations. The current time step calculated temperature of the TC glass layer can be different from the previous time step, but no iterations are done in the current time step for the new TC glass layer temperature. This is an approximation that considers the reaction time of the TC glass layer can be close to EnergyPlus simulation time step say 10 to 15 minutes.
Blinds[LINK]
Window blinds in EnergyPlus are defined as a series of equidistant slats that are oriented horizontally or vertically. All of the slats are assumed to have the same optical properties. The overall optical properties of the blind are determined by the slat geometry (width, separation and angle) and the slat optical properties (frontside and backside transmittance and reflectance). Blind properties for direct radiation are also sensitive to the “profile angle,” which is the angle of incidence in a plane that is perpendicular to the window plane and to the direction of the slats. The blind optical model in EnergyPlus is based on Simmler, Fischer and Winkelmann, 1996; however, that document has numerous typographical errors and should be used with caution.
The following assumptions are made in calculating the blind optical properties:
The slats are flat.
The spectral dependence of interreflections between slats and glazing is ignored; spectralaverage slat optical properties are used.
The slats are perfect diffusers. They have a perfectly matte finish so that reflection from a slat is isotropic (hemispherically uniform) and independent of angle of incidence, i.e., the reflection has no specular component. This also means that absorption by the slats is hemispherically uniform with no incidence angle dependence. If the transmittance of a slat is nonzero, the transmitted radiation is isotropic and the transmittance is independent of angle of incidence.
Interreflection between the blind and wall elements near the periphery of the blind is ignored.
If the slats have holes through which support strings pass, the holes and strings are ignored. Any other structures that support or move the slats are ignored.
Slat Optical Properties[LINK]
The slat optical properties used by EnergyPlus are shown in the following table.
It is assumed that there is no directtodirect transmission or reflection, so that τdir,dir=0, ρfdir,dir=0, and ρbdir,dir=0. It is further assumed that the slats are perfect diffusers, so that τdir,dif, ρfdir,dif and ρbdir,dif are independent of angle of incidence. Until the EnergyPlus model is improved to take into account the angleofincidence dependence of slat transmission and reflection, it is assumed that τdir,dif = τdif,dif, ρfdir,dif = ρfdif,dif, and ρbdir,dif = ρbdif,dif.
Direct Transmittance of Blind[LINK]
The directtodirect and directtodiffuse transmittance of a blind is calculated using the slat geometry shown in Figure 7(a), which shows the side view of one of the cells of the blind. For the case shown, each slat is divided into two segments, so that the cell is bounded by a total of six segments, denoted by s1 through s6 (note in the following that si refers to both segment i and the length of segment i).The lengths of s1 and s2 are equal to the slat separation, h, which is the distance between adjacent slat faces. s3 and s4 are the segments illuminated by direct radiation. In the case shown in Figure 7(a) the cell receives radiation by reflection of the direct radiation incident on s4 and, if the slats have nonzero transmittance, by transmission through s3, which is illuminated from above.
The goal of the blind direct transmission calculation is to determine the direct and diffuse radiation leaving the cell through s2 for unit direct radiation entering the cell through s1.
(a) Side view of a cell formed by adjacent slats showing how the cell is divided into segments, si, for the calculation of direct solar transmittance; (b) side view of a cell showing case where some of the direct solar passes between adjacent slats without touching either of them. In this figure ϕs is the profile angle and ϕb is the slat angle. [fig:asideviewofacellformedbyadjacent]
DirecttoDirect Blind Transmittance[LINK]
Figure 7(b) shows the case where some of the direct radiation passes through the cell without hitting the slats. From the geometry in this figure we see that
τdir,dirbl,f=1−wh,w≤h
where
w=scos(φb−φs)cosφs
Note that we are assuming that the slat thickness is zero. A correction for nonzero slat thickness is described later.
DirecttoDiffuse Blind Transmittance, Reflectance and Absorptance[LINK]
The directtodiffuse and transmittance and reflectance of the blind are calculated using a radiosity method that involves the following three vector quantities:
Ji = the radiosity of segment si, i.e., the total radiant flux into the cell from si
Gi = the irradiance on the cell side of si
Qi = the source flux from the cell side of si
Based on these definitions we have the following equations that relate J, G and Q for the different segments:
J1=Q1J2=Q2J3=Q3+ρbdif,difG3+τdif,difG4J4=Q4+τdif,difG3+ρfdif,difG4J5=Q5+ρbdif,difG5+τdif,difG6J6=Q6+τdif,difG5+ρfdif,difG6
In addition we have the following equation relating G and J:
Gi=6∑j=1JjFji, i=1,6
where Fji is the view factor between sj and si, i.e., Fji is the fraction of radiation leaving sj that is intercepted by si.
Using J1=Q1=0 and J2=Q2=0 and combining the above equations gives the following equation set relating J and Q:
J3−ρbdif,dif6∑j=3JjFj3−τdif,dif6∑j=3JjFj4=Q3
J4−τdif,dif6∑j=3JjFj3−ρfdif,dif6∑j=3JjFj4=Q4
J5−ρbdif,dif6∑j=3JjFj5−τdif,dif6∑j=3JjFj6=Q5
J6−τdif,dif6∑j=3JjFj3−ρfdif,dif6∑j=3JjFj6=Q6
This can be written in the form:
Q′=XJ′
where X is a 4x4 matrix and
J′=⎡⎢ ⎢ ⎢⎣J3J4J5J6⎤⎥ ⎥ ⎥⎦
Q′=⎡⎢ ⎢ ⎢ ⎢⎣Q3Q4Q5Q6⎤⎥ ⎥ ⎥ ⎥⎦
We then obtain J′ from:
J′=X−1Q′
The view factors, Fij, are obtained as follows. The cell we are dealing with is a convex polygon with n sides. In such a polygon the view factors must satisfy the following constraints:
n∑j=1Fij=1, i=1,n
siFij=sjFji,i=1,n; j=1,n
Fii=0, i=1,n
These constraints lead to simple equations for the view factors for n = 3 and 4. For n = 3, we have the following geometry and view factor expression:
View Factor for Three Surfaces
For n = 4 we have:
View Factor for Four Surfaces
Applying these to the slat cell shown in Figure 8 we have the following:
F12=d1+d2−2s2h
F13=h+s3−d32h,etc.
Slat cell showing geometry for calculation of view factors between the segments of the cell. [fig:slatcellshowinggeometryforcalculationof]
The sources for the directtodiffuse transmittance calculation are:
Q1=Q2=Q5=Q6=0(and therefore J1=J2=0)
Q3=τdir,difQ4=ρfdir,dif}φb≤φs+π2(beam hits front of slats)
Q3=ρbdir,difQ4=τdir,dif}φb>φs+π2(beam hits back of slats)
For unit incident direct flux, the front directtodiffuse transmittance and reflectance of the blind are:
τdir,difbl,f=G2ρdir,difbl,f=G1
where
G2=∑6j=3JjFj2G1=∑6j=3JjFj1
and J3 to J6 are given by Equation [eq:QequalsXJprime].
The front direct absorptance of the blind is then:
αdirbl,f=1−τdir,difbl,f−τdir,dirbl,f−ρdir,difbl,f
The directtodiffuse calculations are performed separately for solar and visible slat properties to get the corresponding solar and visible blind properties.
Dependence on Profile Angle[LINK]
The directtodirect and directtodiffuse blind properties are calculated for direct radiation profile angles (see Figure 7) ranging from –90o to +90o in 5o increments. (The “profile angle” is the angle of incidence in a plane that is perpendicular to the window and perpendicular to the slat direction.) In the time step loop the blind properties for a particular profile angle are obtained by interpolation.
Dependence on Slat Angle[LINK]
All blind properties are calculated for slat angles ranging from –90o to +90o in 10o increments. In the timestep loop the slat angle is determined by the slatangle control mechanism and then the blind properties at that slat angle are determined by interpolation. Three slatangle controls are available: (1) slat angle is adjusted to just block beam solar incident on the window; (2) slat angle is determined by a schedule; and (3) slat angle is fixed.
DiffusetoDiffuse Transmittance and Reflectance of Blind[LINK]
To calculate the diffusetodiffuse properties, assuming uniformly distributed incident diffuse radiation, each slat bounding the cell is divided into two segments of equal length (Figure 9), i.e., s3=s4 and s5=s6. For frontside properties we have a unit source, Q1=1. All the other Qi are zero. Using this source value, we apply the methodology described above to obtain G2 and G1. We then have:
τdif,difbl,f=G2ρdif,difbl,f=G1αdifbl,f=1−τdif,difbl,f−ρdif,difbl,f
The backside properties are calculated in a similar way by setting Q2 = 1 with the other Qi equal to zero.
The diffusetodiffuse calculations are performed separately for solar, visible and IR slat properties to get the corresponding solar, visible and IR blind properties.
Slat cell showing arrangement of segments and location of source for calculation of diffusetodiffuse optical properties. [fig:slatcellshowingarrangementofsegments]
Blind properties for sky and ground diffuse radiation[LINK]
For horizontal slats on a vertical window (the most common configuration) the blind diffusetodiffuse properties will be sensitve to whether the radiation is incident upward from the ground or downward from the sky (Figure 10). For this reason we also calculate the following solar properties for a blind consisting of horizontal slats in a vertical plane:
τgnd−dif,difbl,f= front transmittance for ground diffuse solar
τsky−dif,difbl,f= front transmittance for sky diffuse solar
ρgnd−dif,difbl,f= front reflectance for ground diffuse solar
ρsky−dif,difbl,f= front reflectance for sky diffuse solar
αgnd−dif,difbl,f= front absorptance for ground diffuse solar
αsky−dif,difbl,f= front absorptance for sky diffuse solar
These are obtained by integrating over sky and ground elements, as shown in Figure 10, treating each element as a source of direct radiation of irradiance I(ϕs) incident on the blind at profile angle ϕs. This gives:
τsky−dif,difbl,f=π/2∫0[τdir,dirbl,f(ϕs)+τdir,difbl,f(ϕs)]Isky(ϕs)cosϕsdϕsπ/2∫0Isky(ϕs)cosϕsdϕs
ρsky−dif,difbl,f=π/2∫0ρdir,difbl,fIsky(ϕs)cosϕsdϕsπ/2∫0Isky(ϕs)cosϕsdϕs
αsky−difbl,f=π/2∫0αdirbl,fIsky(ϕs)cosϕsdϕsπ/2∫0Isky(ϕs)cosϕsdϕs
Side view of horizontal slats in a vertical blind showing geometry for calculating blind transmission, reflection and absorption properties for sky and ground diffuse radiation. [fig:sideviewofhorizontalslatsinavertical]
We assume that the sky radiance is uniform. This means that Isky is independent of ϕs, giving:
τsky−dif,difbl,f=π/2∫0[τdir,dirbl,f+τdir,difbl,f]cosϕsdϕs
ρsky−dif,difbl,f=π/2∫0ρdir,difbl,fcosϕsdϕs
αsky−difbl,f=π/2∫0αdirbl,fcosϕsdϕs
The corresponding ground diffuse quantities are obtained by integrating ϕs from −π/2 to 0.
An improvement to this calculation would be to allow the sky radiance distribution to be nonuniform, i.e., to depend on sun position and sky conditions, as is done in the detailed daylighting calculation (see “Sky Luminance Distributions” under “Daylight Factor Calculation”).
Correction Factor for Slat Thickness[LINK]
A correction has to be made to the blind transmittance, reflectance and absorptance properties to account for the amount of radiation incident on a blind that is reflected and absorbed by the slat edges (the slats are assumed to be opaque to radiation striking the slat edges). This is illustrated in Figure 11 for the case of direct radiation incident on the blind. The slat crosssection is assumed to be rectangular. The quantity of interest is the fraction, fedge, of direct radiation incident on the blind that strikes the slat edges. Based on the geometry shown in Figure 11 we see that
fedge=tcosγ(h+tcosξ)cosφs=tcos(φs−ξ)(h+tcosξ)cosφs=tsin(φb−φs)(h+tsinφb)cosφs
The edge correction factor for diffuse incident radiation is calculated by averaging this value of fedge over profile angles, φs, from 90o to +90o.
As an example of how the edge correction factor is applied, the following two equations show how blind front diffuse transmittance and reflectance calculated assuming zero slat thickness are modified by the edge correction factor. It is assumed that the edge transmittance is zero and that the edge reflectance is the same as the slat front reflectance, ρf.
τdif,difbl,f→τdif,difbl,f(1−fedge)ρdifbl,f→ρdifbl,f(1−fedge)+fedgeρf