Window Heat Balance Calculation[LINK]
Table 36. Fortran Variables used in Window Heat Balance Calculations
Mathematical variable

Description

Units

FORTRAN variable

N

Number of glass layers


nlayer

σ

StefanBoltzmann constant


sigma

ε_{i}

Emissivity of face i


emis

k_{i}

Conductance of glass layer i

W/m^{2}K

scon

h_{o}, h_{i}

Outside, inside air film convective conductance

W/m^{2}K

hcout, hcout

h_{j}

Conductance of gap j

W/m^{2}K

hgap

T_{o}, T_{i}

Outdoor and indoor air temperatures

K

tout, tin

E_{o}, E_{i}

Exterior, interior longwave radiation incident on window

W/m^{2}

outir, rmir

θ_{i}

Temperature of face i

K

thetas

S_{i}

Radiation (shortwave, and longwave from zone internal sources) absorbed by face i

W/m^{2}

AbsRadGlassFace

I^{ext}_{bm}

Exterior beam normal solar irradiance

W/m^{2}

BeamSolarRad

I^{ext}_{dif}

Exterior diffuse solar irradiance on glazing

W/m^{2}


I^{int}_{sw}

Interior shortwave radiation (from lights and from reflected diffuse solar) incident on glazing from inside

W/m^{2}

QS

I^{int}_{lw}

Longwave radiation from lights and equipment incident on glazing from inside

W/m^{2}

QL

φ

Angle of incidence

radians


A^{f}_{j}

Front beam solar absorptance of glass layer j



A_{j}^{f,dif}, A_{j}^{b,dif}

Front and back diffuse solar absorptance of glass layer j


AbsDiff, AbsDiffBack

A, B

Matrices used to solve glazing heat balance equations

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

Aface, Bface

h_{r,i}

Radiative conductance for face i

W/m^{2}K

hr(i)

Δθ_{i}

Difference in temperature of face i between successive iterations

K


The Glazing Heat Balance Equations[LINK]
The window glass face temperatures are determined by solving the heat balance equations on each face every time step. For a window with N glass layers there are 2N faces and therefore 2N equations to solve. Figure 95 shows the variables used for double glazing (N = 2).
Figure 95. Glazing system with two glass layers showing variables used in heat balance equations.
The following assumptions are made in deriving the heat balance equations:
1) The glass layers are thin enough (a few millimeters) that heat storage in the glass can be neglected; therefore, there are no heat capacity terms in the equations.
2) The heat flow is perpendicular to the glass faces and is one dimensional. See “Edge of Glass Corrections,” below, for adjustments to the gap conduction in multipane glazing to account for 2D conduction effects across the pane separators at the boundaries of the glazing.
3) The glass layers are opaque to IR. This is true for most glass products. For thin plastic suspended films this is not a good assumption, so the heat balance equations would have to be modified to handle this case.
4) The glass faces are isothermal. This is generally a good assumption since glass conductivity is very high.
5) The short wave radiation absorbed in a glass layer can be apportioned equally to the two faces of the layer.
The four equations for doubleglazing are as follows. The equations for single glazing (N = 1) and for N = 3 and N = 4 are analogous and are not shown.
\[{E_o}{\varepsilon_1}  {\varepsilon_1}\sigma \theta_1^4 + {k_1}\left( {{\theta_2}  {\theta_1}} \right) + {h_o}\left( {{T_o}  {\theta_1}} \right) + {S_1} = 0\]
\[{k_1}\left( {{\theta_1}  {\theta_2}} \right) + {h_1}\left( {{\theta_3}  {\theta_2}} \right) + \sigma \frac{{{\varepsilon_2}{\varepsilon_3}}}{{1  \left( {1  {\varepsilon_2}} \right)\left( {1  {\varepsilon_3}} \right)}}\left( {\theta_3^4  \theta_2^4} \right) + {S_2} = 0\]
\[{h_1}({\theta_2}  {\theta_3}) + {k_2}({\theta_4}  {\theta_3}) + \sigma \frac{{{\varepsilon_2}{\varepsilon_3}}}{{1  \left( {1  {\varepsilon_2}} \right)\left( {1  {\varepsilon_3}} \right)}}\left( {\theta_2^4  \theta_3^4} \right) + {S_3} = 0\]
\[{E_i}{\varepsilon_4}  {\varepsilon_4}\sigma \theta_4^4 + {k_2}\left( {{\theta_3}  {\theta_4}} \right) + {h_i}\left( {{T_i}  {\theta_4}} \right) + {S_4} = 0\]
Absorbed Radiation[LINK]
S_{i} in Equations to is the radiation (shortwave and longwave from zone lights and equipment) absorbed on the i^{th} face. Shortwave radiation (solar and shortwave from lights) is assumed to be absorbed uniformly along a glass layer, so for the purposes of the heat balance calculation it is split equally between the two faces of a layer. Glass layers are assumed to be opaque to IR so that the thermal radiation from lights and equipment is assigned only to the inside (roomside) face of the inside glass layer. For N glass layers S_{i} is given by
\[{S_{2j  1}} = {S_{2j}} = \frac{1}{2}\left( {I_{bm}^{ext}\cos \phi A_j^f(\phi ) + I_{dif}^{ext}A_j^{f,dif} + I_{sw}^{{\mathop{\rm int}} }A_j^{b,dif}} \right),{\rm{ }}j = 1{\rm{ to }}N\]
\[{S_{2N}} = {S_{2N}} + {\varepsilon_{2N}}I_{lw}^{{\mathop{\rm int}} }\]
Here
\(I_{bm}^{ext}\) = exterior beam normal solar irradiance
\(I_{dif}^{ext}\) = exterior diffuse solar incident on glazing from outside
\(I_{sw}^{{\mathop{\rm int}} }\) = interior shortwave radiation (from lights and from reflected diffuse solar) incident on glazing from inside
\(I_{lw}^{{\mathop{\rm int}} }\) = longwave radiation from lights and equipment incident on glazing from inside
\({\varepsilon_{2N}}\) = emissivity (thermal absorptance) of the roomside face of the inside glass layer
RoomSide Convection[LINK]
The correlation for roomside convection coefficient, \({h_i}\) , is from ISO 15099 section 8.3.2.2. (Prior to EnergyPlus version 3.1, the value for \({h_i}\) was modeled using the “Detailed” algorithm for opaque surface heat transfer, e.g. for a vertical surface \({h_i} = 1.31{\left {\Delta T} \right^{{\raise0.7ex\hbox{1} \!\mathord{\left/ {\vphantom {1 3}}\right.}\!\lower0.7ex\hbox{3}}}}\) ; see section Detailed Natural Convection Algorithm). The ISO 15099 correlation is for still room air and is determined in terms of the Nusselt number, \(Nu\) , where
\[{h_i} = Nu\left( {\frac{\lambda }{H}} \right)\]
where,
\(\lambda\) is the thermal conductivity of air, and
\(H\) is the height of the window.
The Rayleigh number based on height, \(R{a_H}\) , is calculated using,
\[R{a_H} = \frac{{{\rho ^2}{H^3}g\,{c_p}\left {{T_{surf,i}}  {T_{air}}} \right}}{{{T_{m,f}}\mu \lambda }}\]
where,
\(\rho\) is the density of air
\(g\) is the acceleration due to gravity,
\({c_p}\) is the specific heat of air,
\(\mu\) is the dynamic viscosity of air, and
\({T_{m,f}}\) is the mean film temperature in Kelvin given by,
\[{T_{m,f}} = {T_{air}} + \frac{1}{4}\left( {{T_{surf,i}}  {T_{air}}} \right)\]
There are four cases for the Nusselt correlation that vary by the tilt angle in degrees, \(\gamma\) , and are based on heating conditions. For cooling conditions (where \({T_{surf,i}} > {T_{air}}\) ) the tilt angle is complemented so that \(\gamma = 180  \gamma\)
Case A. \({0^\circ } \le \gamma < 15^\circ\)
\[Nu = 0.13Ra_H^{{\raise0.7ex\hbox{1} \!\mathord{\left/ {\vphantom {1 3}}\right.}\!\lower0.7ex\hbox{3}}}\]
Case B. \(15^\circ \le \gamma \le 90^\circ\)
\[R{a_{cv}} = 2.5 \times {10^5}{\left( {\frac{{{e^{0.72\gamma }}}}{{\sin \lambda }}} \right)^{{\raise0.7ex\hbox{1} \!\mathord{\left/ {\vphantom {1 5}}\right.}\!\lower0.7ex\hbox{5}}}}\]
\[Nu = 0.56{\left( {R{a_H}\sin \gamma } \right)^{{\raise0.7ex\hbox{1} \!\mathord{\left/ {\vphantom {1 4}}\right.}\!\lower0.7ex\hbox{4}}}};\;for\;R{a_H} \le R{a_{CV}}\]
\[Nu = 0.13\left( {Ra_H^{{\raise0.7ex\hbox{1} \!\mathord{\left/ {\vphantom {1 3}}\right.}\!\lower0.7ex\hbox{3}}}  Ra_{CV}^{{\raise0.7ex\hbox{1} \!\mathord{\left/ {\vphantom {1 3}}\right.}\!\lower0.7ex\hbox{3}}}} \right) + 0.56{\left( {R{a_{CV}}\sin \gamma } \right)^{{\raise0.7ex\hbox{1} \!\mathord{\left/ {\vphantom {1 4}}\right.}\!\lower0.7ex\hbox{4}}}};\;R{a_H} > R{a_{CV}}\]
Case C. \(90^\circ < \gamma \le 179^\circ\)
\[Nu = 0.56{\left( {R{a_H}\sin \gamma } \right)^{{\raise0.7ex\hbox{1} \!\mathord{\left/ {\vphantom {1 4}}\right.}\!\lower0.7ex\hbox{4}}}};\;{10^5} \le R{a_H}\sin \gamma < {10^{11}}\]
Case D. \(179^\circ < \gamma \le 180^\circ\)
\[Nu = 0.58Ra_H^{{\raise0.7ex\hbox{1} \!\mathord{\left/ {\vphantom {1 5}}\right.}\!\lower0.7ex\hbox{5}}};\;R{a_H} \le {10^{11}}\]
The material properties are evaluated at the mean film temperature. Standard EnergyPlus pyschrometric functions are used for \(\rho\) and \({c_p}\) . Thermal conductivity is calculated using,
\(\lambda = 2.873 \times {10^{  3}} + 7.76 \times {10^{  5}}{T_{m,f}}\) .
Kinematic viscosity is calculated using,
\(\mu = 3.723 \times {10^{  6}} + 4.94 \times {10^{  8}}{T_{m,f}}\) .
This correlation depends on the surface temperature of the roomside glazing surface and is therefore included inside the window heat balance interation loop.
Solving the Glazing Heat Balance Equations[LINK]
The equations are solved as follows:
1) Linearize the equations by defining \({h_{r,i}} = {\varepsilon_i}\sigma \theta_i^3\) . For example, Equation becomes
\[{E_o}{\varepsilon_1}  {h_{r,1}}{\theta_1} + {k_1}\left( {{\theta_2}  {\theta_1}} \right) + {h_o}\left( {{T_o}  {\theta_1}} \right) + {S_1} = 0\]
2) Write the equations in the matrix form \(A\theta = B\)
3) Use previous time step’s values of \({\theta_i}\) as initial values for the current time step. For the first time step of a design day or run period the initial values are estimated by treating the layers as a simple RC network.
4) Save the \({\theta_i}\) for use in the next iteration: \({\theta_{prev,i}} = {\theta_i}\)
5) Using \({\theta_{2N}}\) , reevaluate the roomside face surface convection coefficient \({h_i}\)
6) Using the \({\theta_i}\) to evaluate the radiative conductances \({h_{r,i}}\)
7) Find the solution \(\theta = {A^{  1}}B\) by LU decomposition
8) Perform relaxation on the the new \({\theta_i}\) : \({\theta_i} \to \left( {{\theta_i} + {\theta_{prev,i}}} \right)/2\)
9) Go to step 4
Repeat steps 4 to 9 until the difference, \(\Delta {\theta_i}\) , between values of the \({\theta_i}\) in successive iterations is less than some tolerance value. Currently, the test is
\[\frac{1}{{2N}}\sum\limits_{i = 1}^{2N} {\Delta {\theta_i}{\rm{ }} < 0.02K}\]
If this test does not pass after 100 iterations, the tolerance is increased to 0.2K. If the test still fails the program stops and an error message is issued.
The value of the inside face temperature, \({\theta_{2N}}\) , determined in this way participates in the zone heat balance solution (see Outdoor/Exterior Convection) and thermal comfort calculation (see Occupant Thermal Comfort).
EdgeOfGlass Effects[LINK]
Table 37. Fortran Variables used in Edge of Glass calculations
Mathematical variable

Description

Units

FORTRAN variable

\(\bar h\)

Areaweighted net conductance of glazing including edgeofglass effects

W/m^{2}K


A_{cg}

Area of centerofglass region

m^{2}

CenterGlArea

A_{fe}

Area of frame edge region

m^{2}

FrameEdgeArea

A_{de}

Area of divider edge region

m^{2}

DividerEdgeArea

A_{tot}

Total glazing area

m^{2}

Surface%Area

h_{cg}

Conductance of centerofglass region (without air films)

W/m^{2}K


h_{fe}

Conductance of frame edge region (without air films)

W/m^{2}K


h_{de}

Conductance of divider edge region (without air films)

W/m^{2}K


h_{ck}

Convective conductance of gap k

W/m^{2}K


h_{rk}

Radiative conductance of gap k

W/m^{2}K


η

Area ratio



α

Conductance ratio


FrEdgeToCenterGlCondRatio, DivEdgeToCenterGlCondRatio

Because of thermal bridging across the spacer separating the glass layers in multipane glazing, the conductance of the glazing near the frame and divider, where the spacers are located, is higher than it is in the center of the glass. The areaweighted net conductance (without inside and outside air films) of the glazing in this case can be written
\[\overline h = \left( {{A_{cg}}{h_{cg}} + {A_{fe}}{h_{fe}} + {A_{de}}{h_{de}}} \right)/{A_{tot}}\]
where
h_{cg} = conductance of centerofglass region (without air films)
h_{fe} = conductance of frame edge region (without air films)
h_{de} = conductance of divider edge region (without air films)
A_{cg} = area of centerofglass region
A_{fe} = area of frame edge region
A_{de} = area of divider edge region
A_{tot} = total glazing area = \({A_{cg}} + {A_{fe}} + {A_{de}}\)
The different regions are shown in Figure 96:
Figure 96: Different types of glass regions.
Equation can be rewritten as
\[\overline h = {h_{cg}}\left( {{\eta_{cg}} + {\alpha_{fe}}{\eta_{fe}} + {\alpha_{de}}{\eta_{de}}} \right)\]
where
\[{\eta_{cg}} = {A_{cg}}/{A_{tot}}\]
\[{\eta_{fe}} = {A_{fe}}/{A_{tot}}\]
\[{\eta_{de}} = {A_{de}}/{A_{tot}}\]
\[{\alpha_{fe}} = {h_{fe}}/{h_{cg}}\]
\[{\alpha_{de}} = {h_{de}}/{h_{cg}}\]
The conductance ratios \({\alpha_{fe}}\) and \({\alpha_{de}}\) are user inputs obtained from Window 5. They depend on the glazing construction as well as the spacer type, gap width, and frame and divider type.
In the EnergyPlus glazing heat balance calculation effective gap convective conductances are used to account for the edgeofglass effects. These effective conductances are determined as follows for the case with two gaps (triple glazing). The approach for other numbers of gaps is analogous.
Neglecting the very small resistance of the glass layers, the centerofglass conductance (without inside and outside air films) can be written as
\[{h_{cg}} = {\left( {{{\left( {{h_{r,1}} + {h_{c,1}}} \right)}^{  1}} + {{\left( {{h_{r,2}} + {h_{c,2}}} \right)}^{  1}}} \right)^{  1}}\]
where
${h_{c,k}} = $ convective conductance of the k^{th} gap
${h_{r,k}} = $ radiative conductance of the k^{th} gap
\[ = \frac{1}{2}\sigma \frac{{{\varepsilon_i}{\varepsilon_j}}}{{1  \left( {1  {\varepsilon_i}} \right)\left( {1  {\varepsilon_j}} \right)}}{\left( {{\theta_i} + {\theta_j}} \right)^3}\]
${_i},{_j} = $ emissivity of the faces bounding the gap
${_i},{_j} = $ temperature of faces bounding the gap (K)
Equation then becomes
\[\overline h = \left( {{\eta_{cg}} + {\alpha_{fe}}{\eta_{fe}} + {\alpha_{de}}{\eta_{de}}} \right){\left( {{{\left( {{h_{r,1}} + {h_{c,1}}} \right)}^{  1}} + {{\left( {{h_{r,2}} + {h_{c,2}}} \right)}^{  1}}} \right)^{  1}}\]
We can also write \(\overline h\) in terms of effective convective conductances of the gaps as
\[\overline h = {\left( {{{\left( {{h_{r,1}} + {{\overline h }_{c,1}}} \right)}^{  1}} + {{\left( {{h_{r,2}} + {{\overline h }_{c,2}}} \right)}^{  1}}} \right)^{  1}}\]
Comparing Eqs. and we obtain
\[{h_{r,k}} + {\overline h_{c,k}} = \left( {{\eta_{cg}} + {\alpha_{fe}}{\eta_{fe}} + {\alpha_{de}}{\eta_{de}}} \right)\left( {{h_{r,k}} + {h_{c,k}}} \right)\]
Using \({\eta_{cg}} = 1  {\eta_{fe}}  {\eta_{de}}\) gives
\[{\overline h_{c,k}} = {h_{r,k}}\left( {{\eta_{fe}}\left( {{\alpha_{fe}}  1} \right) + {\eta_{de}}\left( {{\alpha_{de}}  1} \right)} \right) + {h_{c,k}}\left( {1 + {\eta_{fe}}\left( {{\alpha_{fe}}  1} \right) + {\eta_{de}}\left( {{\alpha_{de}}  1} \right)} \right)\]
This is the expression used by EnergyPlus for the gap convective conductance when a frame or divider is present.
Apportioning of Absorbed ShortWave Radiation in Shading Device Layers[LINK]
If a shading device has a nonzero shortwave transmittance then absorption takes place throughout the shading device layer. The following algorithm is used to apportion the absorbed shortwave radiation to the two faces of the layer. Here f_{1} is the fraction assigned to the face closest to the incident radiation and f_{2} is the fraction assigned to the face furthest from the incident radiation.
\({f_1} = 1,{\rm{ }}{f_2} = 0{\rm{ if }}{\tau_{sh}} = 0\)
Otherwise
\[\begin{array}{cl}
f_1 = 0,\;f_2 = 0 &\text{if} \alpha_{sh} \le 0.01 \\
f_1 = 1,\;f_2 = 0 &\text{if} \alpha_{sh} \gt 0.999 \\
\left.
\begin{array}{l}
f_1 = \frac{1  e^{0.5 \ln (1  \alpha_{sh}}} {\alpha_{sh}} \\
f_2 = 1  {f_1}
\end{array}
\right\} &\text{if} 0.01 lt \alpha_{sh} \le 0.999
\end{array}\]
Window Frame and Divider Calculation[LINK]
For the zone heat balance calculation the inside surface temperature of the frame and that of the divider are needed. These temperatures are determined by solving the heat balance equations on the inside and outside surfaces of the frame and divider.
Table 38. Fortran Variables used in Window/Frame and Divider calculations
Mathematical variable

Description

Units

FORTRAN variable

Q_{ExtIR,abs}

IR from the exterior surround absorbed by outside frame surfaces

W


Q_{IR,emitted}

IR emitted by outside frame surfaces

W


Q_{conv}

Convection from outside air to outside frame surfaces

W


Q_{cond}

Conduction through frame from inside frame surfaces to outside frame surfaces

W


Q_{abs}

Solar radiation plus outside glass IR absorbed by outside of frame

W


Q^{dif}_{abs,sol}

Diffuse solar absorbed by outside frame surfaces, per unit frame face area

W/ m^{2}


Q^{bm}_{abs,sol}

Beam solar absorbed by outside frame surfaces, per unit frame face area

W/ m^{2}


I^{dif}_{ext}

Diffuse solar incident on window

W/ m^{2}


I^{bm}_{ext}

Direct normal solar irradiance

W/ m^{2}


α^{fr}_{sol}

Solar absorptance of frame


FrameSolAbsorp

R_{gl}^{f,dif}

Front diffuse solar reflectance of glazing



R_{gl}^{f,bm}

Front beam solar reflectance of glazing



cos(β_{face})

Cosine of angle of incidence of beam solar on frame outside face


CosIncAng

Cos(β_{h})

Cosine of angle of incidence of beam solar on frame projection parallel to window xaxis


CosIncAngHorProj

Cos(β_{v})

Cosine of angle of incidence of beam solar on frame projection parallel to window yaxis


CosIncAngVertProj

f_{sunlit}

Fraction of window that is sunlit


SunlitFrac

A_{f}

Area of frame’s outside face (same as area of frame’s inside face)

m^{2}


A_{p1}, A_{p2}

Area of frame’s outside and inside projection faces

m^{2}


F_{f}

Form factor of frame’s outside or inside face for IR



F_{p1}, F_{p2}

Form factor of frame outside projection for exterior IR; form factor of frame inside projection for interior IR



E_{o}

Exterior IR incident on window plane

W/m^{2}

outir

E_{i}

Interior IR incident on window plane

W/m^{2}

SurroundIRfromParentZone

ε_{1}, ε_{2}

Outside, inside frame surface emissivity


FrameEmis

θ_{1}, θ_{2}

Frame outside, inside surface temperature

K

FrameTempSurfOut, FrameTempSurfIn

T_{o}, T_{i}

Outdoor and indoor air temperatures

K

tout, tin

h_{o,c}, h_{i,c}

Frame outside and inside air film convective conductance

W/m2K

HOutConv, HInConv

k

Effective insidesurface to outsidesurface conductance of frame per unit area of frame projected onto window plane

W/m^{2}K

FrameConductance, FrameCon

S_{1}

Q_{abs}/A_{f}

W/m^{2}K

FrameQRadOutAbs

S_{2}

Interior shortwave radiation plus interior IR from internal sources absorbed by inside of frame divided by A_{f}

W/m^{2}K

FrameQRadInAbs

η_{1}, η_{2}

A_{p1}/A_{f}, A_{p2}/A_{f}



H

Height of glazed portion of window

m

Surface%Height

W

Width of glazed portion of window

m

Surface%Width

w_{f}, w_{d}

Frame width, divider width

m

FrameWidth, DividerWidth

p_{f1}, p_{f2}

Frame outside, inside projection

m

FrameProjectionOut, FrameProjectionIn

N_{h}, N_{v}

Number of horizontal, vertical dividers


HorDividers, VertDividers

T_{o,r}, T_{i,r}

Frame outside, inside radiative temperature

K

TOutRadFr, TInRadFr

h_{o,r}, h_{i,r}

Frame outside, inside surface radiative conductance

W/m^{2}K

HOutRad, HInRad

A

Intermediate variable in frame heat balance solution

K

Afac

C

Intermediate variable in frame heat balance solution


Efac

B, D

Intermediate variables in frame heat balance solution


Bfac, Dfac

Frame Temperature Calculation[LINK]
Figure 97 shows a cross section through a window showing frame and divider. The outside and inside frame and divider surfaces are assumed to be isothermal. The frame and divider profiles are approximated as rectangular since this simplifies calculating heat gains and losses (see “Error Due to Assuming a Rectangular Profile,” below).
Figure 97. Cross section through a window showing frame and divider (exaggerated horizontally).
Frame Outside Surface Heat Balance[LINK]
The outside surface heat balance equation is
\[{Q_{ExtIR,abs}}  {Q_{IR,emitted}} + {Q_{conv}} + {Q_{cond}} + {Q_{abs}} = 0\]
where
\({Q_{ExtIR,abs}}\) = IR from the exterior surround (sky and ground) absorbed by outside frame surfaces
\({Q_{IR,emitted}}\) = IR emitted by outside frame surfaces
\({Q_{conv}}\) = convection from outside air to outside frame surfaces
\({Q_{cond}}\) = conduction through frame from inside frame surfaces to outside frame surfaces
\({Q_{abs}}\) = solar radiation (from sun, sky and ground) plus IR from outside window surface absorbed by outside frame surfaces (see “Calculation of Absorbed Solar Radiation,” below).
The first term can be written as the sum of the exterior IR absorbed by the outside face of the frame and the exterior IR absorbed by the frame’s outside projection surfaces.
\[{Q_{ExtIR,abs}} = {\varepsilon_1}{E_o}{A_f}{F_f} + {\varepsilon_1}{E_o}{A_{p1}}{F_{p1}}\]
where ε_{1} is the outside surface emissivity.
The exterior IR incident on the plane of the window, E_{o}, is the sum of the IR from the sky, ground and obstructions. For the purposes of the frame heat balance calculation it is assumed to be isotropic. For isotropic incident IR, F_{f} = 1.0 and F_{p1} = 0.5, which gives
\[{Q_{ExtIR,abs}} = {\varepsilon_1}{E_o}({A_f} + \frac{1}{2}{A_{p1}})\]
The IR emitted by the outside frame surfaces is
\[{Q_{ExtIR,emitted}} = {\varepsilon_1}\sigma ({A_f} + {A_{p1}})\theta_1^4\]
The convective heat flow from the outside air to the outside frame surfaces is
\[{Q_{conv}} = {h_{o,c}}({A_f} + {A_{p1}})({T_o}  {\theta_1})\]
The conduction through the frame from inside to outside is
\[{Q_{cond}} = k{A_f}({\theta_2}  {\theta_1})\]
Note that A_{f} is used here since the conductance, k, is, by definition, per unit area of frame projected onto the plane of the window.
Adding these expressions for the Q terms and dividing by* A_{f}* gives
\[{E_0}{\varepsilon_1}\left( {1 + \frac{1}{2}{\eta_1}} \right)  {\varepsilon_1}\left( {1 + {\eta_1}} \right)\theta_1^4 + {h_{o,c}}\left( {1 + {\eta_1}} \right)\left( {{T_0}  {\theta_1}} \right) + k\left( {{\theta_2}  {\theta_1}} \right) + {S_1} = 0\]
where S_{1} = Q_{abs}/A_{f} and
\[{\eta_1} = \frac{{{A_{p1}}}}{{{A_f}}} = \left( {\frac{{{p_{f,1}}}}{{{w_f}}}} \right)\frac{{H + W  ({N_h} + {N_v}){w_d}}}{{H + W + 2{w_f}}}\]
We linearize Eq. as follows.
Write the first two terms as
\[{\varepsilon_1}(1 + {\eta_1})[{E_o}(1 + \frac{1}{2}{\eta_1})/(1 + {\eta_1})  \theta_1^4]\]
and define a radiative temperature
\[{T_{o,r}} = {[{E_o}(1 + \frac{1}{2}{\eta_1})/(1 + {\eta_1})]^{1/4}}\]
This gives
\[{\varepsilon_1}(1 + {\eta_1})[T_{o,r}^4  \theta_1^4]\]
which, within a few percent, equals
\[{\varepsilon_1}(1 + {\eta_1})\frac{{{{({T_{o,r}} + {\theta_1})}^3}}}{2}({T_{o,r}}  {\theta_1})\]
Defining an outside surface radiative conductance as follows
\[{h_{o,r}} = {\varepsilon_1}(1 + {\eta_1})\frac{{{{({T_{o,r}} + {\theta_1})}^3}}}{2}\]
then gives
\[{h_{o,r}}({T_{o,r}}  {\theta_1})\]
The final outside surface heat balance equation in linearized form is then
\[{h_{o,r}}({T_{o,r}}  {\theta_1}) + {h_{o,c}}(1 + {\eta_1})({T_o}  {\theta_1}) + k({\theta_2}  {\theta_1}) + {S_1} = 0\]
Frame Inside Surface Heat Balance[LINK]
A similar approach can be used to obtain the following linearized inside surface heat balance equation:
\[{h_{i,r}}({T_{i,r}}  {\theta_2}) + {h_{i,c}}(1 + {\eta_2})({T_i}  {\theta_2}) + k({\theta_1}  {\theta_2}) + {S_2} = 0\]
where
\[{T_{i,r}} = {[{E_i}(1 + \frac{1}{2}{\eta_2})/(1 + {\eta_2})]^{1/4}}\]
\[{\eta_2} = \frac{{{A_{p2}}}}{{{A_f}}} = \left( {\frac{{{p_{f,2}}}}{{{w_f}}}} \right)\frac{{H + W  ({N_h} + {N_v}){w_d}}}{{H + W + 2{w_f}}}\]
and E_{i} is the interior IR irradiance incident on the plane of the window.
Solving Eqs. and simultaneously gives
\[{\theta_2} = \frac{{D + CA}}{{1  CB}}\]
with
\[A = \frac{{{h_{o,r}}{T_{o,r}} + {h_{o,c}}{T_o} + {S_1}}}{{{h_{o,r}} + k + {h_{o,c}}}}\]
\[B = \frac{k}{{{h_{o,r}} + k + {h_{o,c}}}}\]
\[C = \frac{k}{{{h_{i,r}} + k + {h_{i,c}}}}\]
\[D = \frac{{{h_{i,r}}{T_{i,r}} + {h_{i,c}}{T_i} + {S_2}}}{{{h_{i,r}} + k + {h_{i,c}}}}\]
Calculation of Solar Radiation Absorbed by Frame[LINK]
The frame outside face and outside projections and inside projections absorb beam solar radiation (if sunlight is striking the window) and diffuse solar radiation from the sky and ground. For the outside surfaces of the frame, the absorbed diffuse solar per unit frame face area is
\[Q_{abs,sol}^{dif} = I_{ext}^{dif}{\alpha_{fr,sol}}\left( {{A_f} + {F_{p1}}{A_{p1}}} \right)/{A_f} = I_{ext}^{dif}{\alpha_{fr,sol}}\left( {1 + 0.5\frac{{{A_{p1}}}}{{{A_f}}}} \right)\]
If there is no exterior window shade, I^{dif}_{ext} includes the effect of diffuse solar reflecting off of the glazing onto the outside frame projection, i.e.,
\[I_{ext}^{dif} \to I_{ext}^{dif}\left( {1 + R_{gl}^{f,dif}} \right)\]
The beam solar absorbed by the outside face of the frame, per unit frame face area is
\[Q_{abs,sol}^{bm,face} = I_{ext}^{bm}{\alpha_{fr,sol}}\cos {\beta_{face}}{f_{sunlit}}\]
The beam solar absorbed by the frame outside projection parallel to the window xaxis is
\[Q_{abs,sol}^{bm,h} = I_{ext}^{bm}{\alpha_{fr,sol}}\cos {\beta_h}{p_{f1}}\left( {W  {N_v}{w_d}} \right){f_{sunlit}}/{A_f}\]
Here it is assumed that the sunlit fraction, f_{sunlit}, for the window can be applied to the window frame. Note that at any given time beam solar can strike only one of the two projection surfaces that are parallel to the window xaxis. If there is no exterior window shade, I^{bm}_{ext} includes the effect of beam solar reflecting off of the glazing onto the outside frame projection, i.e.,
\[I_{ext}^{bm} \to I_{ext}^{bm}\left( {1 + R_{gl}^{f,bm}} \right)\]
The beam solar absorbed by the frame outside projection parallel to the window yaxis is
\[Q_{abs,sol}^{bm,v} = I_{ext}^{bm}{\alpha_{fr,sol}}\cos {\beta_v}{p_{f1}}\left( {H  {N_h}{w_d}} \right){f_{sunlit}}/{A_f}\]
Using a similar approach, the beam and diffuse solar absorbed by the inside frame projections is calculated, taking the transmittance of the glazing into account.
Error Due to Assuming a Rectangular Profile[LINK]
Assuming that the inside and outside frame profile is rectangular introduces an error in the surface heat transfer calculation if the profile is nonrectangular. The percent error in the calculation of convection and emitted IR is approximately 100\(\left {{L_{profile,rect}}  {L_{profile,actual}}} \right/{L_{profile,rect}}\) , where L_{profile,rect} is the profile length for a rectangular profile (w_{f} +* p_{f1}* for outside of frame or w_{f} + p_{f2}for inside of frame) and L_{profile,actual} is the actual profile length. For example, for a circular profile vs a square profile the error is about 22%. The error in the calculation of absorbed beam radiation is close to zero since the beam radiation intercepted by the profile is insensitive to the shape of the profile. The error in the absorbed diffuse radiation and absorbed IR depends on details of the shape of the profile. For example, for a circular profile vs. a square profile the error is about 15%.
Divider Temperature Calculation[LINK]
The divider inside and outside surface temperatures are determined by a heat balance calculation that is analogous to the frame heat balance calculation described above.
Beam Solar Reflection from Window Reveal Surfaces[LINK]
This section describes how beam solar radiation that is reflected from window reveal surfaces is calculated. Reflection from outside reveal surfaces—which are associated with the setback of the glazing from the outside surface of the window’s parent wall—increases the solar gain through the glazing. Reflection from inside reveal surfaces—which are associated with the setback of the glazing from the inside surface of the window’s parent wall—decreases the solar gain to the zone because some of this radiation is reflected back out of the window.
The amount of beam solar reflected from reveal surfaces depends, among other things, on the extent to which reveal surfaces are shadowed by other reveal surfaces. An example of this shadowing is shown in Figure 98. In this case the sun is positioned such that the top reveal surfaces shadow the left and bottom reveal surfaces. And the right reveal surfaces shadow the bottom reveal surfaces. The result is that the left/outside, bottom/outside, left/inside and bottom/inside reveal surfaces each have sunlit areas. Note that the top and right reveal surfaces are facing away from the sun in this example so their sunlit areas are zero.
Figure 98. Example of shadowing of reveal surfaces by other reveal surfaces.
The size of the shadowed areas, and the size of the corresponding illuminated areas, depends on the following factors:
The sun position relative to the window
The height and width of the window
The depth of the outside and inside reveal surfaces
We will assume that the reveal surfaces are perpendicular to the window plane and that the window is rectangular. Then the above factors determine a unique shadow pattern. From the geometry of the pattern the shadowed areas and corresponding illuminated areas can be determined. This calculation is done in subroutine CalcBeamSolarReflectedFromWinRevealSurface in the SolarShading module. The window reveal input data is specified in the WindowProperty:FrameAndDivider object expect for the depth of the outside reveal, which is determined from the vertex locations of the window and its parent wall.
If an exterior shading device (shade, screen or blind) is in place it is assumed that it blocks beam solar before it reaches outside or inside reveal surfaces. Correspondingly, it is assumed that an interior or betweenglass shading device blocks beam solar before it reaches inside reveal surfaces.
Representative shadow patterns are shown in Figure 99 for a window with no shading device, and without and with a frame. The case with a frame has to be considered separately because the frame can cast an additional shadow on the inside reveal surfaces.
The patterns shown apply to both vertical and horizontal reveal surfaces. It is important to keep in mind that, for a window of arbitrary tilt, if the left reveal surfaces are illuminated the right surfaces will not be, and vice versa. And if the bottom reveal surfaces are illuminated the top surfaces will not be, and vice versa. (Of course, for a vertical window, the top reveal surfaces will never be illuminated by beam solar if the reveal surfaces are perpendicular to the glazing, as is being assumed.
For each shadow pattern in Figure 99, equations are given for the shadowed areas \({A_{1,sh}}\) and \({A_{2,sh}}\) of the outside and inside reveal surfaces, respectively. The variables in these equations are the following (see also Figure 100):
\({d_1}\) = depth of outside reveal, measured from the outside plane of the glazing to the edge of the reveal, plus one half of the glazing thickness.
\({d_2}\) = depth of inside reveal (or, for illumination on bottom reveal surfaces, inside sill depth), measured from the inside plane of the glazing to the edge of the reveal or the sill, plus one half of the glazing thickness.
\(L\) = window height for vertical reveal surfaces or window width for horizontal reveal surfaces
\(\alpha\) = vertical solar profile angle for shadowing on vertical reveal surfaces or horizontal solar profile angle for shadowing on horizontal reveal surfaces.
\({p_1}({p_2})\) = distance from outside (inside) surface of frame to glazing midplane.
\({d_2}^\prime\) = depth of shadow cast by top reveal on bottom reveal, or by left reveal on right reveal, or by right reveal on left reveal.
\({d_2}^{\prime \prime }\) = depth of shadow cast by frame.
For simplicity it is assumed that, for the case without a frame, the shadowed and illuminated areas extend into the glazing region. For this reason, \({d_1}\) and \({d_2}\) are measured from the midplane of the glazing. For the case with a frame, the beam solar absorbed by the surfaces formed by the frame outside and inside projections perpendicular to the glazing is calculated as described in “Window Frame and Divider Calculation: Calculation of Solar Radiation Absorbed by Frame.”
Figure 99. Expression for area of shaded regions for different shadow patterns: (a) window without frame, (b) window with frame
Figure 100. Vertical section through a vertical window with outside and inside reveal showing calculation of the shadows cast by the top reveal onto the inside sill and by the frame onto the inside sill.
The following logic gives expressions for the shadowed areas for all possible shadow patterns. Here:
d1 = \({d_1}\)
d2 = \({d_2}\)
P1 = \({p_1}\)
P2 = \({p_2}\)
f1 = \({d_1}  {p_1}\)
f2 = \({d_2}  {p_2}\)
d2prime = \({d_2}^\prime\)
d2prime2 = \({d_2}^{\prime \prime }\)
d12 = \({d_1} + {d_2}  {d_2}^\prime\)
TanAlpha = \(\tan \alpha\)
A1sh = \({A_{1,sh}}\)
A2sh = \({A_{2,sh}}\)
L = \(L\)
L1 = average distance to frame of illuminated area of outside reveal (used to calculate view factor to frame).
L2 = average distance to frame of illuminated area of inside reveal (used to calculate view factor to frame).
IF(window does not have a frame) THEN
IF(d2prime < = d2) THEN
IF(d12*TanAlpha < = L) THEN
A1sh = 0.5*TanAlpha*d1**2
A2sh = d2prime*L + 0.5*TanAlpha*d12**2  A1sh
ELSE ! d12*TanAlpha > L
IF(d1*TanAlpha < = L) THEN
A1sh = 0.5*TanAlpha*d1**2
A2sh = d2*L  0.5*TanAlpha*(L/TanAlpha  d1)**2
ELSE ! d1*TanAlpha > L
A1sh = d1*L  (0.5/TanAlpha)*L**2
A2sh = d2*L
END IF
END IF
ELSE ! d2prime > d2
A2sh = d2*L
IF(d2prime < d1+d2) THEN
IF(d12*TanAlpha < = L) THEN
A1sh = L*(d2primed2) + 0.5*TanAlpha*d12**2
ELSE ! d12*TanAlpha > L
A1sh = d1*L  0.5*L**2/TanAlpha
END IF
ELSE ! d2prime > = d1+d2
A1sh = d1*L
END IF
END IF
ELSE ! Window has a frame
f1 = d1P1
f2 = d2P2
d2prime2 = FrameWidth/TanGamma
IF(vertical reveal) THEN ! Vertical reveal
IF(InsReveal+0.5*GlazingThickness < = P2) d2 = P2 + 0.001
ELSE ! Horizontal
IF(bottom reveal surfaces may be illuminated) THEN
! Bottom reveal surfaces may be illuminated
IF(InsSillDepth+0.5*GlazingThickness< = P2) d2 = P2 + 0.001
ELSE
! Top reveal surfaces may be illuminated
IF(InsReveal+0.5*GlazingThickness < = P2) d2 = P2 + 0.001
END IF
END IF
IF(d2prime < = f2) THEN
! Shadow from opposing reveal does not go beyond inside
! surface of frame
IF(d12*TanAlpha < = L) THEN
A1sh = 0.5*TanAlpha*f1**2
L1 = f1*(f1*TanAlpha/(6*L)+0.5)
IF(d2(d2prime+d2prime2+P2) > = 0.) THEN
A2sh = (d2prime+d2prime2)*L + &
0.5*TanAlpha*((d1+d2d2prime)**2d1+p2+d2prime2)**2)
L2 = d2prime2 + 0.5*(d2(d2prime+d2prime2+P2))
ELSE ! d2(d2prime+d2prime2+P2) < 0.
! Inside reveal is fully shadowed by frame and/or
!opposing reveal
A2sh = f2*L
L2 = f2
END IF
ELSE ! d12*TanAlpha > = L
IF((d1+P2)*TanAlpha < = L) THEN
A1sh = 0.5*TanAlpha*f1**2
L1 = f1*((f1*TanAlpha)/(6*L) + 0.5)
IF((d1+P2+d2prime2)*TanAlpha > = L) THEN
A2sh = f2*L
L2 = f2
ELSE ! (d1+P2+d2prime2)*TanAlpha < L
A2sh = f2*L  0.5*(L(d1+P2)*TanAlpha)**2/TanAlpha &
+ d2prime2*(L(d1+P2+d2prime2/2)*TanAlpha)
L2 = d2prime2 + (L/TanAlpha  (d1+P2+d2prime2))/3
END IF
ELSE ! (d1+P2)*TanAlpha > L
L2 = f2
A2sh = f2*L
IF(f1*TanAlpha < = L) THEN
A1sh = 0.5*TanAlpha*f1**2
L1 = f1*((f1*TanAlpha)/(6*L) + 0.5)
ELSE ! f1*TanAlpha > L
A1sh = f1*L  0.5*L**2/TanAlpha
L1 = f1(L/TanAlpha)/3
END IF
END IF
END IF
ELSE
! d2prime > f2  Shadow from opposing reveal goes beyond
! inside of frame
A2sh = f2*L
L2 = f2
IF(d2prime > = d1+d2) THEN
A1sh = 0.0
L1 = f1
ELSE ! d2prime < d1+d2
IF(d2prime < = d2+P1) THEN
IF(f1*TanAlpha < = L) THEN
A1sh = 0.5*TanAlpha*f1**2
L1 = f1*((f1*TanAlpha)/(6*L) + 0.5)
ELSE ! f1*TanAlpha > L
A1sh = f1*L  0.5*L**2/TanAlpha
L1 = f1  (L/TanAlpha)/3
END IF
ELSE ! d2prime > d2+P1
IF(d12*TanAlpha < = L) THEN
A1sh = L*(d2prime(d2+P1)) + 0.5*TanAlpha*d12**2
L1 = (L*(f1d12/2)d12*TanAlpha* &
(f1/2d12/3))/(Ld12*TanAlpha/2)
ELSE ! d12*TanAlpha > L
A1sh = f1*L  0.5*L**2/TanAlpha
L1 = f1  (L/TanAlpha)/3
END IF
END IF
END IF
END IF
FracToGlassOuts = 0.5*(1.0  ATAN(FrameWidth/L1)/PiOvr2)
FracToGlassIns = 0.5*(1.0  ATAN(FrameWidth/L2)/PiOvr2)
END IF ! End of check if window has frame
The beam solar reflected from a sunlit region of area \(A\) is given by
\[R = {I_B}A\cos \beta (1  a){\rm{ }}\]
where
\(R\) = reflected solar radiation [W]
\({I_B}\) = beam normal irradiance [W/m^{2}]
\(A\) = sunlit area [m^{2}]
\(\beta\) = beam solar angle of incidence on reveal surface
\(a\) = solar absorptance of reveal surface
All reflected radiation is assumed to be isotropic diffuse. For outside reveal surfaces it is assumed that \(R/2\) goes toward the window and \(R/2\) goes to the exterior environment. Of the portion that goes toward the window a fraction \({F_1}\) goes toward the frame, if present, and \(1  {F_1}\) goes toward the glazing.
The view factor \({F_1}\) to the frame calculated by assuming that the illuminated area can be considered to be a line source. Then the areaweighted average distance, \({L_1}\) , of the source to the frame is calculated from the shape of the illuminated area (see above psuedocode). Then \({F_1}\) is related as follows to the average angle subtended by the frame of width \({w_f}\) :
\[{F_1} = \frac{{{{\tan }^{  1}}({w_f}/{L_1})}}{{\pi /2}}\]
For the portion going towards the frame, \((R/2){F_1}{a_f}\) is absorbed by the frame (where \({a_f}\) is the solar absorptance of the frame) and contributes to the frame heat conduction calculation. The rest, \((R/2){F_1}(1  {a_f})\) , is assumed to be reflected to the exterior environment.
If the glazing has diffuse transmittance \({\tau_{diff}}\) , diffuse front reflectance \(\rho_{diff}^f\) , and layer front absorptance \(\alpha_{l,diff}^f\) , then, of the portion, \((R/2)(1  {F_1})\) , that goes toward the glazing, \((R/2)(1  {F_1}){\tau_{diff}}\) is transmitted to the zone, \((R/2)(1  {F_1})\alpha_{l,diff}^f\) is absorbed in glass layer \(l\) and contributes to the glazing heat balance calculation, and \((R/2)(1  {F_1})\rho_{dif{f_{}}}^f\) is reflected to the exterior environment.
The beam solar absorbed by an outside reveal surface is added to the other solar radiation absorbed by the outside of the window’s parent wall.
For inside reveal surfaces it is assumed that \(R/2\) goes towards the window and \(R/2\) goes into the zone. Of the portion that goes toward the window a fraction \((R/2){F_2}\) goes toward the frame, if present, and \((R/2)(1  {F_2})\) goes toward the glazing (\({F_2}\) is calculated using a method analogous to that used for \({F_1}\) ). For the portion going towards the frame, \((R/2){F_2}{a_f}\) is absorbed by the frame and contributes to the frame heat conduction calculation. The rest, \((R/2){F_2}(1  {a_f})\) , is assumed to be reflected back into the zone.
If the glazing has diffuse back reflectance \(\rho_{diff}^b\) , and layer back absorptance \(\alpha_{l,diff}^b\) , then, of the portion \((R/2)(1  {F_2})\) that goes toward the glazing, \((R/2)(1  {F_2}){\tau_{diff}}\) is transmitted back out the glazing, \((R/2)(1  {F_2})\alpha_{l,diff}^b\) is absorbed in glass layer \(l\) and contributes to the glazing heat balance calculation, and \((R/2)(1  {F_2})\rho_{dif{f_{}}}^b\) is reflected into the zone.
The beam solar absorbed by an inside reveal surface is added to the other solar radiation absorbed by the inside of the window’s parent wall.
Shading Device Thermal Model[LINK]
Shading devices in EnergyPlus can be on the exterior or interior sides of the window or between glass layers. The window shading device thermal model accounts for the thermal interactions between the shading layer (shade, screen or blind) and the adjacent glass, and between the shading layer and the room (for interior shading) or the shading layer and the outside surround (for exterior shading).
An important feature of the shading device thermal model is calculating the natural convection airflow between the shading device and glass. This flow affects the temperature of the shading device and glazing and, for interior shading, is a determinant of the convective heat gain from the shading layer and glazing to the zone air. The airflow model is based on one described in the ISO Standard 15099, “Thermal Performance of Windows, Doors and Shading Devices—Detailed Calculations” [ISO15099, 2001]. (Betweenglass forced airflow is also modeled; see “Airflow Windows.”)
The following effects are considered by the shading device thermal model:
For interior and exterior shading device: Longwave radiation (IR) from the surround absorbed by shading device, or transmitted by the shading device and absorbed by the adjacent glass. For interior shading the surround consists of the other zone surfaces. For exterior shading the surround is the sky and ground plus exterior shadowing surfaces and exterior building surfaces “seen” by the window.
Interreflection of IR between the shading device and adjacent glass.
Direct and diffuse solar radiation absorbed by the shading device.
Interreflection of solar radiation between shading layer and glass layers.
Convection from shading layer and glass to the air in the gap (or, for betweenglass shading, gaps) between the shading layer and adjacent glass, and convection from interior shading layer to zone air or from exterior shading layer to outside air.
Natural convection airflow in the gap (or, for betweenglass shading, gaps) between shading layer and adjacent glass induced by buoyancy effects, and the effect of this flow on the shadingtogap and glasstogap convection coefficients.
For interior shading, convective gain (or loss) to zone air from gap airflow.
In the following it is assumed that the shading device, when in place, covers the glazed part of the window (and dividers, if present) and is parallel to the glazing. For interior and exterior shading devices it is assumed that the shading layer is separated from the glazing by an air gap. A betweenglass shading layer is assumed to be centered between two glass layers and separated from the adjacent glass layers by gaps that is filled with the same gas. If the window has a frame, it is assumed that the shading device does not cover the frame.
Heat Balance Equations for Shading Device and Adjacent Glass[LINK]
If a window shading device is deployed the heat balance equations for the glass surfaces facing the shading layer are modified, and two new equations, one for each face of the shading layer, are added. Figure 101 illustrates the case of double glazing with an interior shading device.
Figure 101. Glazing system with two glass layers and an interior shading layer showing variables used in heat balance equations.
The heat balance equation for the glass surface facing the gap between glass and shading layer (called in the following, “gap”) is
\[\frac{{{E_i}{\varepsilon_4}{\tau_{sh}}}}{{1  {\rho_4}{\rho_{sh}}}} + \frac{{\sigma {\varepsilon_4}}}{{1  {\rho_4}{\rho_{sh}}}}\left[ {\theta_5^4{\varepsilon_{sh}}  \theta_4^4(1  {\rho_{sh}})} \right] + {k_2}\left( {{\theta_3}  {\theta_4}} \right) + {h_{cv}}\left( {{T_{gap}}  {\theta_4}} \right) + {S_4} = 0\]
where
τ_{sh} = IR diffuse transmittance of shading device
ε_{sh} = diffuse emissivity of shading device
ρ_{sh} = IR diffuse reflectance of shading device ( = 1  ( τ_{sh} + ε_{sh}))
θ_{5} = temperature of the surface of the shading layer that faces the gap (K).
The term 1 – ρ_{4} ρ_{sh} accounts for the interreflection of IR radiation between glass and shading layer.
The convective heat transfer from glass layer #2 to the air in the gap is
\[{q_{c,gl}} = {h_{cv}}({\theta_4}  {T_{gap}})\]
where
T_{gap} = effective mean temperature of the gap air (K).
h_{cv} = convective heat transfer coefficient from glass or shading layer to gap air (W/m^{2}K).
The corresponding heat transfer from shading layer to gap air is
\[{q_{c,sh}} = {h_{cv}}({\theta_5}  {T_{gap}})\]
The convective heat transfer coefficient is given by
\[{h_{cv}} = 2{h_c} + 4v\]
where
h_{c} = surfacetosurface heat transfer coefficient for nonvented (closed) cavities (W/m^{2}K)
v = mean air velocity in the gap (m/s).
The quantities h_{cv} and T_{gap} depend on the airflow velocity in the gap, which in turn depends on several factors, including height of shading layer, glass/shading layer separation (gap depth), zone air temperature for interior shading or outside air temperature for exterior shading, and shading layer and glass face temperatures. The calculation of h_{cv} and T_{gap} is described in the following sections.
The heat balance equation for the shading layer surface facing the gap is
\[\frac{E_i \tau_{sh} \rho_4 \varepsilon_{sh}} {1  \rho_4 \rho_{sh}} +
\frac{\sigma \varepsilon_{sh}} {1  \rho_4 \rho_{sh}} \left[
\varepsilon_4 \theta_4^4  \theta_5^4 (
1  \rho_4(
\varepsilon_{sh} + \rho_{sh}
)
)
\right] +
k_{sh} \left(
\theta_6  \theta_5
\right) +
h_{cv} \left(
T_{gap}  \theta_5
\right) +
S_{sh,1} = 0\]
where
k_{sh} = shading layer conductance (W/m^{2}K).
θ_{6} = temperature of shading layer surface facing the zone air (K).
S_{sh,1} = solar radiation plus shortwave radiation from lights plus IR radiation from lights and zone equipment absorbed by the gapside face of the shading layer (W/m^{2}K).
The heat balance equation for the shading layer surface facing the zone air is
\[{E_i}{\varepsilon_{sh}}  {\varepsilon_{sh}}\sigma \theta_6^4 + {k_{sh}}\left( {{\theta_5}  {\theta_6}} \right) + {h_i}\left( {{T_i}  {\theta_6}} \right) + {S_{sh,2}} = 0\]
where
S_{sh,2} = solar radiation plus shortwave radiation from lights plus IR radiation from lights and zone equipment absorbed by the zoneside face of the shading layer (W/m^{2}K).
Solving for Gap Airflow and Temperature[LINK]
For interior and exterior shading devices a pressurebalance equation is used to determine gap air velocity, gap air mean equivalent temperature and gap outlet air temperature given values of zone air temperature (or outside temperature for exterior shading), shading layer face temperatures and gap geometry. The pressure balance equates the buoyancy pressure acting on the gap air to the pressure losses associated with gap airflow between gap inlet and outlet [ISO15099, 2001]. The variables used in the following analysis of the interior shading case are shown in Figure 102.
Figure 102. Vertical section (a) and perspective view (b) of glass layer and interior shading layer showing variables used in the gap airflow analysis. The opening areas A_{bot}, A_{top}, A_{l}, A_{r} and A_{h} are shown schematically.
Pressure Balance Equation[LINK]
The pressure balance equation for airflow through the gap is
\[\Delta {p_T} = \Delta {p_B} + \Delta {p_{HP}} + \Delta {p_Z}\]
Here, Δp_{T}is the driving pressure difference between room air and gap air. It is given by
\[\Delta {p_T} = {\rho_0}{T_0}gH\sin \phi \frac{{\left {{T_{gap}}  {T_{gap,in}}} \right}}{{{T_{gap}}{T_{gap,in}}}}\]
where
ρ_{0} = density of air at temperature T_{0} (kg/m^{3})
T_{0} = reference temperature (283K)
g = acceleration due to gravity (m/s^{2})
H = height of shading layer (m)
φ = tilt angle of window (vertical = 90^{o})
T_{gap} = effective mean temperature of the gap air (K)
T_{gap,in} = gap inlet temperature ( = zone air temperature for interior shading) (K)
The Δp_{B} term is due to the acceleration of air to velocity v (Bernoulli’s law). It is given by
\[\Delta {p_B} = \frac{\rho }{2}{v^2}{\rm{ (Pa)}}\]
where ρis the gap air density evaluated at T_{gap} (kg/m^{3}).
The Δp_{HP} term represents the pressure drop due to friction with the shading layer and glass surfaces as the air moves through the gap. Assuming steady laminar flow, it is given by the HagenPoiseuille law for flow between parallel plates [Munson et al. 1998]:
\[\Delta {p_{HP}} = 12\mu \frac{H}{{{s^2}}}{v^{}}{\rm{ (Pa)}}\]
where μ is the viscosity of air at temperature T_{gap} (Pas).
The Δp_{Z} term is the sum of the pressure drops at the inlet and outlet openings:
\[\Delta {p_Z} = \frac{{\rho {v^2}}}{2}\left( {{Z_{in}} + {Z_{out}}} \right){\rm{ (Pa)}}\]
Here, the inlet pressure drop factor,* Z_{in}, and the outlet pressure drop factor,Z_{out}*, are given by
\[\begin{array}{rl}
Z_{in} & = \left( \frac{A_{gap}}{0.66A_{eq,in}}  1 \right)^2 \\
Z_{out} & = \left( \frac{A_{gap}}{0.60A_{eq,out}}  1 \right)^2
\end{array}\]
\[\begin{array}{l}{Z_{in}} = {\left( {\frac{{{A_{gap}}}}{{0.66{A_{eq,in}}}}  1} \right)^2}\\\{Z_{out}} = {\left( {\frac{{{A_{gap}}}}{{0.60{A_{eq,out}}}}  1} \right)^2}\end{array}\]
where
A_{eq,in} = equivalent inlet opening area (m^{2})
A_{eq,out} = equivalent outlet opening area (m^{2})
A_{gap} = crosssectional area of the gap = sW(m^{2})
If T_{gap} > T_{gap,in}
\[\begin{array}{rl}
A_{eq,in} & = A_{bot} + \frac{A_{top}}{2(A_{bot}+A_{top}} \left( A_l + A_r + A_h \right) \\
A_{eq,out} & = A_{top} + \frac{A_{bot}}{2(A_{bot}+A_{top}} \left( A_l + A_r + A_h \right)
\end{array}\]
If T_{gap} ≤ T_{gap,in}
\[\begin{array}{rl}
A_{eq,in} & = A_{top} + \frac{A_{bot}}{2(A_{bot}+A_{top}} \left( A_l + A_r + A_h \right)
A_{eq,out} & = A_{bot} + \frac{A_{top}}{2(A_{bot}+A_{top}} \left( A_l + A_r + A_h \right) \\
\end{array}\]
Here, the area of the openings through which airflow occurs (see Figure 102 and Figure 103) are defined as follows:
A_{bot} = area of the bottom opening (m^{2})
A_{top} = area of the top opening (m^{2})
A_{l} = area of the leftside opening (m^{2})
A_{r} = area of the rightside opening (m^{2})
A_{h} = air permeability of the shading device expressed as the total area of openings (“holes”) in the shade surface (these openings are assumed to be uniformly distributed over the shade) (m^{2})
Figure 103 shows examples of A_{bot}, A_{top}, A_{l} and A_{r} for different shading device configurations. These areas range from zero to a maximum value equal to the associated shade/screen/blindtoglass crosssectional area; i.e., A_{bot} and A_{top} ≤ sW, A_{l} and A_{r} ≤ sH.
Figure 103. Examples of openings for an interior shading layer covering glass of height H and width W. Not to scale. (a) Horizontal section through shading layer with openings on the left and right sides (top view). (b) Vertical section through shading layer with openings at the top and bottom (side view).
Expression for the Gap Air Velocity
Expressing Equation in terms of v yields the following quadratic equation:
\[\frac{{\rho {v^2}}}{2}\left( {1 + {Z_{in}} + {Z_{out}}} \right) + \frac{{12\mu H}}{{{s^2}}}v  {\rho_0}{T_0}gH\sin \phi \frac{{\left {{T_{gap,in}}  {T_{gap}}} \right}}{{{T_{gap,in}}{T_{gap}}}} = 0\]
Solving this gives
\[v = \frac{{{{\left[ {{{\left( {\frac{{12\mu H}}{{{s^2}}}} \right)}^2} + \frac{{2{\rho ^2}(1 + {Z_{in}} + {Z_{out}}){\rho_0}{T_0}gH\sin \phi \left {{T_{gap,in}}  {T_{gap}}} \right}}{{{T_{gap,in}}{T_{gap}}}}} \right]}^{1/2}}  \frac{{12\mu H}}{{{s^2}}}}}{{\rho (1 + {Z_{in}} + {Z_{out}})}}\]
The choice of the root of the quadratic equation is dictated by the requirement that v = 0 if T_{gap,in} = T_{gap}.
Gap Outlet Temperature and Equivalent Mean Air Temperature
The temperature of air in the gap as a function of distance, h, from the gap inlet (Figure 104) is
\[{T_{gap}}(h) = {T_{ave}}  ({T_{ave}}  {T_{gap,in}}){e^{  h/{H_0}}}\]
where
\[{T_{ave}} = \frac{{{T_{gl}} + {T_{sh}}}}{2}\]
is the average temperature of the glass and shading layer surfaces facing the gap (K).
H_{0} = characteristic height (m), given by
\[{H_0} = \frac{{\rho {C_p}s}}{{2{h_{cv}}}}v\]
where C_{p} is the heat capacity of air.
The gap outlet temperature is given by
\[{T_{gap,out}} = {T_{ave}}  ({T_{ave}}  {T_{gap,in}}){e^{  H/{H_0}}}\]
The thermal equivalent mean temperature of the gap air is
\[{T_{gap}} = \frac{1}{H}\int\limits_0^H {{T_{gap}}(h)dh = {T_{ave}}  \frac{{{H_0}}}{H}} \left( {{T_{gap,out}}  {T_{gap,in}}} \right)\]
Figure 104. Variation of gap air temperature with distance from the inlet for upward flow.
Solution Sequence for Gap Air Velocity and Outlet Temperature
The routine WinShadeGapFlow is called within the glazing heat balance iterative loop in SolveForWindowTemperatures to determine v and T_{gap,out}. The solution sequence in WinShadeGapFlow is as follows:
At start of iteration, guess T_{gap} as ((T_{gl} + T_{sh})/2 + T_{gap,in})/2. Thereafter use value from previous iteration.
Get stillair conductance, h_{c}, by calling WindowGasConductance and NusseltNumber.
Get v from Equation
Get h_{cv} from Equation
Get T_{ave}from Equation
Get T_{gap,out} from Equation
Get new value of T_{gap}from Equation
The values of h_{cv} and T_{gap} so determined are then used in the window heat balance equations to find new values of the face temperatures of the glass and shading layers. These temperatures are used in turn to get new values of h_{cv} and T_{gap} until the whole iterative process converges.
Convective Heat Gain to Zone from Gap Airflow
The heat added (or removed) from the air as it passes through the gap produces a convective gain (or loss) to the zone air given by
\[{q_v} = LW\left[ {{h_{cv}}\left( {{T_{gl}}  {T_{gap}}} \right) + {h_{cv}}\left( {{T_{sh}}  {T_{gap}}} \right)} \right] = 2{h_{cv}}LW\left( {{T_{ave}}  {T_{gap}}} \right){\rm{ (W)}}\]
This can also be expressed as
\[{q_v} = \dot m{C_p}\left( {{T_{gap,out}}  {T_{gap,in}}} \right){\rm{ (W)}}\]
where the air mass flow rate in the gap is given by
\[\dot m = \rho {A_{gap}}v{\rm{ (kg/s)}}\]
Heat Balance Equations for BetweenGlass Shading Device[LINK]
In EnergyPlus shading devices are allowed between the two glass panes of double glazing and between the two inner glass panes of triple glazing. Figure 105 shows the case of a betweenglass shading device in double glazing.
Figure 105. Glazing system with two glass layers and a betweenglass shading device showing variables used in the heat balance equations.
The heat balance equations for the two glass surfaces facing the shading device are the following.
For face #2:
\[h_{cv,1} (T_{gap,1}\theta_2) + k_1 (\theta_1\theta_2) +
\frac{\sigma \epsilon_2}{1\rho_2 R_1} \left[
\frac{\tau_{sh}}{1\rho_6\rho_3} (
\epsilon_3 \theta_3^4 + \epsilon_6\theta_6^4\rho_3
)
+ \epsilon_5\theta_5^4 + \epsilon_2\theta_2^4R_1
\right]\]
where
\[{R_1} = {\rho_5} + \frac{{\tau_{sh}^2{\rho_3}}}{{1  {\rho_6}{\rho_3}}}\]
${T_{gap,1}} = $ effective mean air temperature in gap 1 (K)
${h_{cv,1}} = $ convective heat transfer coefficient from glass or shading layer to gas in gap 1 (W/m^{2}K)
For face #3:
\[h_{cv,2} (T_{gap,2}\theta_3) + k_2 (\theta_4\theta_3) +
\frac{\sigma \epsilon_3}{1\rho_3 R_2} \left[
\frac{\tau_{sh}}{1\rho_5\rho_2} (
\epsilon_3 \theta_2^4 + \epsilon_5\theta_5^4\rho_2
)
+ \epsilon_6\theta_6^4 + \epsilon_7\theta_3^4R_2
\right]\]
where
\[{R_2} = {\rho_6} + \frac{{\tau_{sh}^2{\rho_2}}}{{1  {\rho_5}{\rho_2}}}\]
${T_{gap,2}} = $ effective mean air temperature in gap 2 (K)
${h_{cv,2}} = $ convective heat transfer coefficient from glass or shading layer to gas in gap 2 (W/m^{2}K)
The heat balance equations for the shading layer faces are:
For face #5:
\[h_{cv,1} (T_{gap,1}\theta_5) + k_sh (\theta_6\theta_5) +
\frac{\sigma \epsilon_5}{1\rho_2 R_1} \left[
\frac{\tau_{sh}\rho_2}{1\rho_5\rho_3} (
\epsilon_3 \theta_3^4 + \epsilon_6\theta_6^4\rho_3
)
+ \epsilon_2\theta_2^4 + \epsilon_5\theta_5^4\rho_2
\right]\]
For face #6:
\[h_{cv,2} (T_{gap,2}\theta_3) + k_2 (\theta_4\theta_3) +
\frac{\sigma \epsilon_3}{1\rho_3 R_2} \left[
\frac{\tau_{sh}}{1\rho_5\rho_2} (
\epsilon_2 \theta_2^4 + \epsilon_5\theta_5^4\rho_2
)
+ \epsilon_2\theta_2^4 + \epsilon_7\theta_3^4R_2
\right]\]
The convective heat transfer coefficients are given by
\[\begin{array}{l}{h_{cv,1}} = 2{h_{c,1}} + 4v\\\{h_{cv,2}} = 2{h_{c,2}} + 4v\end{array}\]
where
${h_{c,1}},{h_{c,2}} = $ surfacetosurface heat transfer coefficients for gap #1 and #2, respectively, when these gaps are nonvented (closed).
$v = $ air velocity in the gaps (m/s). It is assumed that the gap widths are equal, so that the velocity of flow in the gaps is equal and opposite, i.e., when the airflow is upward in gap #1 it is downward in gap #2 and viceversa.
In analogy to the interior or exterior shading device case, the air velocity is determined by solving the following pressure balance equation:
\[\Delta {p_{T,1,2}} = \Delta {p_{B,1}} + \Delta {p_{HP,1}} + \Delta {p_{Z,1}} + \Delta {p_{B,2}} + \Delta {p_{HP,2}} + \Delta {p_{Z,2}}\]
where the driving pressure difference between gap #1 and #2 is
\[\Delta {p_{T,1,2}} = {\rho_0}{T_0}gH\sin \phi \frac{{\left {{T_{gap,1}}  {T_{gap,2}}} \right}}{{{T_{gap,1}}{T_{gap,2}}}}{\rm{ (Pa)}}\]
The pressure drops on the righthand side of this equation are:
\[\begin{array}{l}\Delta {p_{B,i}} = \frac{{{\rho_{gap,}}_i}}{2}{v^2}\\\Delta {p_{HP,i}} = 12{\mu_{gap,}}_i\frac{H}{{{s^2}}}\\\Delta {p_{Z,i}} = \frac{{{\rho_{gap,}}_i{v^2}}}{2}({Z_{in,i}} + {Z_{out,i}})\end{array}\]
where i = gap number (1 or 2).
It can be shown that \({Z_{in,1}} + {Z_{out,1}} = {Z_{in,2}} + {Z_{out,2}}\) . Then, inserting these pressure drop expressions into , we obtain the following expression for the airflow velocity:
\[v = \frac{{{{\left[ {{{\left( {\frac{{12({\mu_{gap,}}_1 + {\mu_{gap,}}_2)H}}{{{s^2}}}} \right)}^2} + 2\Delta {p_{T,1,2}}({\rho_{gap,1}} + {\rho_{gap,2}})(1 + {Z_{in}} + {Z_{out}})} \right]}^{1/2}}  \frac{{12({\mu_{gap,}}_1 + {\mu_{gap,}}_2)H}}{{{s^2}}}}}{{({\rho_{gap,}}_1 + {\rho_{gap,}}_2)(1 + {Z_{in}} + {Z_{out}})}}\]
The choice of the sign of the square root term is dictated by the requirement that \(v = 0\) if \(\Delta {p_{T,1,2}} = 0\) , i.e., \({T_{gap,1}} = {T_{gap,2}}\) .
Given \(v\) we can now calculate \({T_{gap,1}}\) and \({T_{gap,2}}\) , which gives \(\Delta {p_{T,1,2}}\) . The procedure is as follows. We have
\[{T_{gap,1,out}} = {T_{ave,1}}  ({T_{ave,1}}  {T_{gap,1,in}}){\xi_1}\]
where \({T_{ave,1}} = ({\theta_2} + {\theta_5})/2\) and \({\xi_1} = {e^{  \frac{H}{{{H_{0,1}}}}}}\) with \({H_{0,1}} = {\rho_{gap,1}}{C_p}sv/(2{h_{cv,1}})\) . Since \({T_{gap,1,in}} = {T_{gap,2,out}}\) this gives:
\({T_{gap,1,out}} = {T_{ave,1}}  ({T_{ave,1}}  {T_{gap,2,out}}){\xi_1}\)
Similarly,
\[{T_{gap,2,out}} = {T_{ave,2}}  ({T_{ave,2}}  {T_{gap,1,out}}){\xi_2}\]
Solving these simultaneous equations gives:
\(\begin{array}{l}{T_{gap,1,out}} = \frac{{{T_{ave,1}}(1  {\xi_1}) + {\xi_1}{T_{ave,2}}(1  {\xi_2})}}{{1  {\xi_1}{\xi_2}}}\\\\\{T_{gap,2,out}} = \frac{{{T_{ave,2}}(1  {\xi_2}) + {\xi_2}{T_{ave,1}}(1  {\xi_1})}}{{1  {\xi_1}{\xi_2}}}\end{array}\)
Using these in
\[{T_{gap,1}} = {T_{ave,1}}  \frac{{{H_{0,1}}}}{H}\left( {{T_{gap,1,out}}  {T_{gap,2,out}}} \right)\]
gives
\({T_{gap,1}} = {T_{ave,1}}  \frac{{{H_{0,1}}}}{H}\xi \left( {{T_{ave,1}}  {T_{ave,2}}} \right)\)
with
\[\xi = \frac{{(1  {\xi_1})(1  {\xi_2})}}{{1  {\xi_1}{\xi_2}}}\]
Similarly, from
\[{T_{gap,2}} = {T_{ave,2}}  \frac{{{H_{0,2}}}}{H}\left( {{T_{gap,2,out}}  {T_{gap,1,out}}} \right)\]
we get
\[{T_{gap,2}} = {T_{ave,2}}  \frac{{{H_{0,2}}}}{H}\xi \left( {{T_{ave,2}}  {T_{ave,1}}} \right)\]
The overall solution sequence is as follows.
At start of iteration guess \({T_{gap,1}} = {T_{ave,1}}\) and \({T_{gap,2}} = {T_{ave,2}}\) . Then
1) Get \({\rho_{gap,1}},{\rho_{gap,2}},{\mu_{gap,1}},{\mu_{gap,2}}\) using \({T_{gap,1}},{T_{gap,2}}\) .
2) Get stillair conductances \({h_{c,1}},{h_{c,2}}\) by calling WindowGasConductance and NusseltNumber.
3) Get \(v\) from Equation
4) Get \({h_{cv,1}},{h_{cv,2}}\) from Equation
5) Get \({T_{ave,1}},{T_{ave,2}}\)
6) Get H_{o,1}, H_{o,2}, \({\xi_1}\) and \({\xi_2}\) .
7) Get \({T_{gap,1}}\) ,\({T_{gap,2}}\) from Equations and
The values \({h_{cv,1}}{\rm{, }}{h_{cv,2}}{\rm{, }}{T_{gap,1}}\) and \({T_{gap,2}}\) are then used in the face heat balance equations to find new values of the face temperatures \({\theta_2}{\rm{, }}{\theta_3},{\rm{ }}{\theta_5}\) and \({\theta_6}\) . These are used in turn to get new values of \({h_{cv,1}}{\rm{, }}{h_{cv,2}}{\rm{, }}{T_{gap,1}}\) and \({T_{gap,2}}\) until the whole iterative process converges.
Airflow Windows[LINK]
In airflow windows forced air flows in the gap between adjacent layers of glass. Such windows are also known as “heatextract windows” and “climate windows.”
Five configurations of airflow windows are modeled (Figure 106) that depend on the source and destination of forced air. The allowed combinations of Airflow Source and Airflow Destination are:
InsideAir à OutsideAir
InsideAir à InsideAir
InsideAir à ReturnAir
OutsideAir à InsideAir
OutsideAir à OutsideAir
Figure 106. Gap airflow configurations for airflow windows. From “Active facades,” Version no. 1, Belgian Building Research Institute, June 2002.
A common application of airflow windows is to reduce the zone cooling load by exhausting indoor air through the window, thereby picking up and rejecting heat from the glazing (Figure 106).
Figure 107 shows the variables used in the heat balance equations for forced airflow in a doubleglazed window.
Figure 107. Glazing system with forced airflow between two glass layers showing variables used in the heat balance equations.
The heat balance equation for the lefthand glass surface facing the gap in Figure 107 is:
\[{k_1}({\theta_1}  {\theta_2}) + {h_{cv}}({T_{gap}}  {\theta_2}) + \sigma \frac{{{\varepsilon_2}{\varepsilon_3}}}{{1  (1  {\varepsilon_2})(1  {\varepsilon_3})}}\left( {\theta_3^4  \theta_2^4} \right) + {S_2} = 0\]
The corresponding equation for the righthand glass surface facing the gap is:
\[{k_2}({\theta_4}  {\theta_3}) + {h_{cv}}({T_{gap}}  {\theta_3}) + \sigma \frac{{{\varepsilon_2}{\varepsilon_3}}}{{1  (1  {\varepsilon_2})(1  {\varepsilon_3})}}\left( {\theta_2^4  \theta_3^4} \right) + {S_3} = 0\]
Here,
T_{gap} = effective mean temperature of the gap air (K)
h_{cv} = convective heat transfer coefficient from glass to gap air (W/m2K).
The convective heat transfer coefficient is given by
\[{h_{cv}} = 2{h_c} + 4v\]
where
h_{c} = glasstoglass heat transfer coefficient for nonvented (closed) cavity (W/m^{2}K)
v = mean air velocity in the gap (m/s).
The air velocity is determined by the gap crosssectional area in the flow direction and the air flow rate, which is an input value that is constant or can vary according to a userspecified schedule:
\[v = \frac{F}{{{A_{gap}}}}{\rm{ (m/s)}}\]
where
F = airflow rate (m^{3}/s)
A_{gap} = gap crosssectional area in direction of flow (m2)
It is assumed that the airflow is uniform across the width of the window.
The mean temperature of the gap air is given by the following expression, whose derivation follows that for for the case of an interior shading device:
\[{T_{gap}} = {T_{ave}}  \frac{H}{{{H_0}}}({T_{gap,out}}  {T_{gap,in}})\]
where
\[{T_{ave}} = \frac{{{\theta_2} + {\theta_3}}}{2}\]
H = glazing height (m)
\[{H_0} = \frac{{\rho {C_p}s}}{{2{h_{cv}}}}v\]
T_{gap,in}_{ } = gap air inlet temperature (T_{i} if the airflow source is indoor air, T_{o} if the airflow source is outside air) (K)
The outlet air temperature is given by
\[{T_{gap,out}} = {T_{ave}}  ({T_{ave}}  {T_{gap,in}}){e^{  H/{H_0}}}\]
The equations for glass face #1 and #4 are the same as those for no airflow in the gap (Equations and ).
The convective heat gain to the zone air due to the gap airflow when the airflow destination is indoor air is
\[{q_v} = \dot m({C_{p,out}}{T_{gap,out}}  {C_{p,i}}{T_i}){\rm{ (W)}}\]
where
C_{p,i}_{ } = heat capacity of the indoor air (J/kgK)
C_{p,out} = heat capacity of the gap outlet air (J/kgK)
and where the air mass flow rate in the gap is
\[\dot m = \rho F{\rm{ (kg/s)}}\]
The fan energy used to move air through the gap is very small and is ignored.
Airflow Window with BetweenGlass Shading Device[LINK]
Figure 108 shows the case of a doubleglazed airflow window with a between glass shading device. The heat balance equations in this case are the same as those for the betweenglass shading device with natural convection (Figure 105 and following equations) except that now
\[v = \frac{{F/2}}{{{A_{gap}}}}{\rm{ (m/s)}}\]
where A_{gap} = sW is the crosssectional area of the gap on either side of the shading device. It is assumed that the shading device is centered between the two panes of glass so that the airflow, F, is divided equally between the two gaps.
The convective heat gain to the zone air due to the airflow through the two gaps when the airflow destination is indoor air is
\[{q_v} = \dot m({C_{p,ave,out}}{T_{gap,ave,out}}  {C_{p,i}}{T_i}{\rm{) (W)}}\]
where the average temperature of the two outlet air streams is
\[{T_{gap,ave,out}}{\rm{ = }}({T_{gap,1,out}} + {T_{gap,2,out}})/2\]
and
C_{p,ave,out}_{ } = heat capacity of the outlet air evaluated at T_{gap,ave,out} (J/kgK)
Figure 108. Airflow window with betweenglass shading device showing variables used in the heat balance equations.
Evacuated Glazing Unit (EGU)[LINK]
Evacuated glazing Units (EGU) are an emerging technology, developed as a concept some 2030 years ago, but only now approaching widespread commercialization and adoption.
Evacuated glazing unit (EGU): cavity thermal conductance[LINK]
Thermal Conductance of the space in an evacuated glazing unit (EGU) is the sum of the conductance of the low pressure gas (air) and radiation conductance.
\[{C_{gas}} = {C_{COND}} + {C_{RAD}}\]
Conductance of the lowpressure gas
Conductance of low pressure gasses is calculated using formula by Corrucini (Corruccini, R. (1959)).
\[{C_{COND}} = \alpha \cdot \left[ {\frac{{\gamma + 1}}{{\gamma  1}}} \right] \cdot {\left[ {\frac{R}{{8 \cdot \pi \cdot M \cdot T}}} \right]^{\frac{1}{2}}} \cdot P\]
\[B = \alpha \cdot \left[ {\frac{{\gamma + 1}}{{\gamma  1}}} \right] \cdot {\left[ {\frac{R}{{8 \cdot \pi \cdot M \cdot T}}} \right]^{\frac{1}{2}}}\]
\[{C_{COND}} = B \cdot P\]
Where:
\[\alpha = \frac{{{\alpha_1}{\alpha_2}}}{{{\alpha_2} + {\alpha_1}\left( {1  {\alpha_2}} \right)}}\]
a_{1}, a_{2} = Accommodation coefficients of the gas molecules with the two surfaces. These values depend on the temperature, surface conditions, etc. For the present configuration and conditions, it is expected that a is approximately 0.5. If conservative value is needed than value of 1.0 could be used. With a_{1}, a_{2} = 0.5, a = 0.333
g = Specific heat ratio, g_{air} = 1.40.
R = Universal gas constant, R = 8,314.462175 J/mol·K
M = Molecular Weight, M_{air} = 28.97 [mol/g]
T = (T_{1} + T_{2})/2 [K]
P = Pressure of the gas [N/m^{2}]
From the paper Collins and Robinson (Collins, R., & Robinson, S. (1991)), B is set at approximately 50 for Air, if pressure is given in torr. Therefore according to Collins and Robinson, for air and approximate conditions of EGU:
C_{COND} ≈ 50·P
Where P is in torr (i.e., mm Hg).
Note: Conversion from Pa to torr is accomplished by multiplying value in torr by 133.28.
Using formula 2 and assuming T_{1} to be 20 ºC and T_{2} to be 18 ºC (expected temperatures of glass surfaces in EGU, if one glass surface is lowe and unit is exposed to NFRC standard environmental conditions), and using SI system of units, the B is calculated as 54.4, which is very close to the value of 50, proposed by Collins and Robinson
C_{COND} ≈ 54.4·P
We will use exact value, calculated by the formula, so values of will be input into the calculations, which enables more flexible model that can account for special treatment of glass surfaces.
Radiation conductance[LINK]
Radiation conductance for the two parallel plates is given by:
\[{C_{RAD}} = \frac{1}{{{\varepsilon_1}^{  1} + {\varepsilon_2}^{  1}  1}}\sigma \frac{{\left( {{T_1}^4  {T_2}^4} \right)}}{{{T_1}  {T_2}}}\]
Where:
e_{1} = emissivity of the first facing glass surface, [  ]
e_{2} = emissivity of the second facing glass surface, [  ]
s = StefanBoltzmann Constant, 5.67 x 10^{8}, [W/(m^{2}·K^{4})]
T_{1} = Temperature of the first facing glass surface, [K]
T_{2} = Temperature of the second facing glass surface. [K]
Assuming glass surface temperatures of 20 ºC and 18 ºC, respectively, the following radiation conductances are obtained for the three different glass emissivities:
Clear Glass (e_{1} = e_{2} = 0.84): C_{RAD} = 3.4 W/m^{2}K
Hard Coat Lowe (e_{1} or e_{2} = 0.15): C_{RAD} = 0.68 W/m^{2}K
Soft Coat Lowe (e_{1} or e_{2} = 0.04): C_{RAD} = 0.19 W/m^{2}K
Note: C_{RAD} of 0.09 is theoretically possible using best lowe technology today (i.e., e_{1} or e_{2} = 0.02).
Note: Lowe values above are typical values, which will vary by manufacturer. Some more recent hard coat lowe values are at or below 0.1.
It should also be noted that values above are based on the fix set of temperatures, while in reality temperatures will depend on the environmental conditions and surface emissivities (e.g., it cannot be expected that clear glass will have same T_{1} and T_{2} as lowe glass).
Equation is precise formulation for two parallel plates at the constant temperature. Simplified equation under these conditions is given in the form of:
\[{C_{RAD}} = 4\sigma \frac{1}{{{\varepsilon_1}^{  1} + {\varepsilon_2}^{  1}  1}}{T_m}^3\]
Where:
T_{m} = mean temperature, [K]
\[{T_m} = \frac{{\left( {{T_1} + {T_2}} \right)}}{2}\]
Calculation of the Ufactor
\[U = \frac{1}{{{R_{tot}}}} = \frac{1}{{{R_o} + 2 \cdot {R_{glass}} + {R_{gap}} + {R_i}}}\]
\[{R_{gap}} = \frac{1}{{{C_{gap}}}} = \frac{1}{{{C_{COND}} + {C_{RAD}}}}\]
\[{R_{glass}} = \frac{{{t_{glass}}}}{{{k_{glass}}}}\]
Where:
t_{glass} = glass thickness; [m]
k_{glass} = glass conductivity; k_{glass} = 1 W/(m·K)
R_{glass} = 0.003 m^{2}K/W (for 3 mm glass pane)
R_{o} ≈ 0.033 m^{2}K/W
R_{i} ≈ 0.14 m^{2}K/W
Ufactor of EGU without any pillars (pretending that this is possible) would be calculated using C_{RAD} only. From above radiation conductance calculations:
Clear Glass: U = 2.64 W/(m^{2}·K) [0.464 Btu/(hr·ft^{2}·ºF)]
Hard Coat Lowe: U = 0.62 W/(m^{2}·K) [0.109 Btu/(hr·ft^{2}·ºF)]
Soft Coat Lowe: U = 0.19 W/(m^{2}·K) [0.034 Btu/(hr·ft^{2}·ºF)]
Adding conductance of the air at 0.001 torr (C_{COND} = 0.08 W/(m^{2}·K), these values become:
Clear Glass: U = 2.66 W/(m^{2}·K) [0.468 Btu/(hr·ft^{2}·ºF)]
Hard Coat Lowe: U = 0.68 W/(m^{2}·K) [0.120 Btu/(hr·ft^{2}·ºF)]
Soft Coat Lowe: U = 0.27 W/(m^{2}·K) [0.048 Btu/(hr·ft^{2}·ºF)]
In contrast, the Ufactor of the same configuration with the air at atmospheric pressure will be (For the space width of 50 mm, C_{COND} ≈ 450 W/m^{2}K):
Clear Glass: U = 5.52 W/(m^{2}·K) [0.468 Btu/(hr·ft^{2}·ºF)]
Hard Coat Lowe: U = 5.52 W/(m^{2}·K) [0.120 Btu/(hr·ft^{2}·ºF)]
Soft Coat Lowe: U = 5.52 W/(m^{2}·K) [0.048 Btu/(hr·ft^{2}·ºF)]
It is clear that emissivity of the glass surface makes no difference, because of the dominant conductance of the air space. Also, it is worth noting that the Ufactor of such configuration is very close to the Ufactor of single glazing.
Evacuated glazing unit (EGU): Glass support element thermal conductance[LINK]
Glass panes in the EGU are separated by an array of small support elements. Typically, these support elements have a cylindrical shape and are often referred to as “pillars.” Typical geometry of the pillar is 0.51.0 mm diameter and 0.05 mm (50 mm) height. They are typically spaced 12 in. apart in a form of square or staggered matrix.
Calculation of the Ufactor[LINK]
The conductance of these elements can be measured or numerically modeled to determine accurate thermal performance. Approximate method also exists and is based on the combination of modeling and analytical work for the conduction through small cylinders in contact with infinite parallel plates with thickness much larger than cylinder height.
The following formula can be used to determine conductance of the single pillar, C_{p} (Collins and FisherCripps 1991):
\[{C_p} = \frac{{2ka}}{{1 + \frac{{2h}}{{\pi a}}}}\]
Where:
k = conductivity of glass, W/(m·K)
a = radius of the pillar (m)
h = pillar height, m
For the square array of support pillars (Collins and FischerCripps 1991) proposes the following formula for their conductance, C_{pa}:
\[{C_{pa}} = \frac{{2ka}}{{{\lambda ^2}\left( {1 + \frac{{2h}}{{\pi a}}} \right)}}\]
Where:
l = pillar spacing, m
This formula is approximate and does not include effect of the conductivity of the pillar, but it gives good approximation for common materials used in this technology, since conductivity of the pillar does not play substantial role for noninsulating pillars (where “noninsulating” would mean that conductivity of the pillar is equal or higher than the conductivity of the glass pane.
The Ufactor of the EGU with support pillars is then:
\[U = \frac{1}{{{R_{tot}}}} = \frac{1}{{{R_o} + 2 \cdot {R_{glass}} + {R_{gap}} + {R_i}}}\]
Where:
\[{R_{gap}} = \frac{1}{{{C_{gap}}}} = \frac{1}{{{C_{COND}} + {C_{RAD}} + {C_{pa}}}}\]
Collins, R., & FischerCripps, A. 1991. “Design of Support Pillar Arrays in Flat Evacuated Windows.”. Aust. J. Phys.
Collins, R., & Robinson, S. 1991. “Evacuated Glazing”. Solar Energy. Vol. 47, No. 1, pp. 2738.
Corruccini, R. 1959. “Gaseous Heat Conduction at Low Pressures and Temperatures”. Vacuum. Vol. 78, pp.1929.
Deflection of insulated glazing unit (IGU) can result in thermal performance degradation or improvement due to the reduction or increase of gap space width. Convection of the gas fill is affected by changed gap space and due to modified convection pattern and shorter or longer thermal path at the center of the glazing unit can result in increased or decreased thermal performance. For the most part, Ufactor is mostly affected as a direct result of changed thermal performance; however note that solar heat gain through the window (SHGC) can also be affected because of the effect of inward flowing fraction of absorbed solar radiation, which is affected by thermal performance of the IGU.
Deflection in sealed IGU is caused by the difference in gas pressure in IGU gap vs. outdoor/indoor pressure. Indoor and outdoor pressure can be considered equal, since indoor building environment is in pretty good contact with outdoor environment. We will call this pressure an atmospheric pressure, P_{a}. The differences in pressure between atmospheric and gap pressure is due to several factors, listed here:
1) Difference in atmospheric pressure between IGU fabrication location and end use location
2) Difference in temperature during fabrication and actual operating conditions for the glazing. It should be noted that initial temperature can be higher than ambient temperature during fabrication process, due to elevated sealant temperatures, which can raise local temperatures within the IGU.
3) Unbalanced gas fill leakage through the sealants, resulting in lower gap pressure and inward deflection.
4) Wind or static load pressure
Effects 1 and 2 will be modeled using equations presented below, while effect 3 does not have credible mathematical model. However, cumulative deflection, resulting from all three effects can be measured in the field and its effect on thermal performance can be modeled by specifying center glazing deflection.
Wind or static load pressure effects on deflection is not included in this model at this time, but will be considered for future versions.
Recognizing that indoor and outdoor air pressure could be different, such as in hot box test environment, future plans for the extension of the model will include option to specify different values for indoor and outdoor pressure. Another future improvement to the model will also include linking certain air gaps with indoor or outdoor environment, meaning that respective pressures in linked spaces will be set to equal.
Mathematical Model[LINK]
Mathematical model described in detail here is based on the research work by Bernier and Bourret (1997) and Timoshenko and WoinowskyKrieger (1959). Bernier and Bourret (1997) of the Ecole Polytechnique Montréal adopted Timoshenko and WoinowskyKrieger (1959) model for calculating flat plate deflection subjected to the differential pressure field (static), while their original contribution was to develop correlations for changes in thermal performance, based on IGU deflection at the center of glazing location. In addition to adopting Bernier and Bourret (1997) model here, we have also developed model for calculating change in thermal performance of deflected units when this deflection is measured in the field. Therefore, the mathematical formulation, presented here is divided into two sections; 1) calculation of the deflection and resulting thermal performance caused by pressure and temperature effects and 2) calculation of the thermal performance of the IGU when the deflection is measured.
If coordinate system is set as shown in Figure 109 and Figure 110, it is possible to calculate deflection distribution at each point of pane by using following equation:
Figure 109. Deflection Coordinate System  2D
\[{L_{D\left( i \right)}}\left( {x,y} \right) = \frac{{16\cdot \Delta {P_{\left( i \right)}}}}{{{\pi ^6}\cdot {D_{\left( i \right)}}}}\mathop \sum \limits_{m = 1,3,5 \ldots }^\infty \mathop \sum \limits_{n = 1,3,5 \ldots }^\infty \frac{{\sin \frac{{m\pi x}}{W}\sin \frac{{n\pi y}}{H}}}{{mn{{\left( {{{\left( {\frac{m}{W}} \right)}^2} + {{\left( {\frac{n}{H}} \right)}^2}} \right)}^2}}}\]
Where,
\[{D_{\left( i \right)}} = \frac{{E\cdot t_{\left( i \right)}^3}}{{12\cdot \left( {1  {\nu ^2}} \right)}}\]
Where,
E = Young’s modulus (7.2 x 10^{10}) [Force per unit Area; SI: Pa, IP: psi]
t = thickness of glazing pane [Length; SI: m, IP: in.]
n = poison’s ratio (0.22 for glass) [NonDimensional]
Figure 110. Deflection Coordinate System  3D
ΔP_{i} = P_{gap(i)}  P_{gap(i1)} (for ith pane) [Force per unit Area; SI: Pa, IP: psi]
ΔP_{i} = P_{gap(1)}  P_{a} (first pane) [Force per unit Area; SI: Pa, IP: psi]
ΔP_{i} = P_{a}  P_{gap(n1)} (last pane) [Force per unit Area; SI: Pa, IP: psi]
Where,
P_{a} = atmospheric pressure. [Force per unit Area; SI: Pa, IP: psi]
\[{P_{gap\left( i \right)}} = \frac{{{P_{ini}}{V_{ini\left( i \right)}}{T_{gap\left( i \right)}}}}{{{T_{ini}}{V_{gap\left( i \right)}}}}\]
Where,
P_{ini} = Initial pressure. Applies to all gaps as a single value (input data  measured or otherwise) [Force per unit Area; SI: Pa, IP: psi]
T_{ini} = Initial temperature. Applies to all gaps as a single value (input data  measured or otherwise) [Degree Temperature; SI: K, IP: R]
V_{ini(i)} = Initial volume for ith gap. [Length*Length*Length; SI: m^{3}, IP: in^{3}]
\[{V_{ini\left( i \right)}} = {L_i}\cdot W\cdot H\]
Where,
L_{i} = nondeflected glazing gap width (for ith gap) [Length; SI: m, IP: in.]
W = IGU width [Length; SI: m, IP: in.]
H = IGU height [Length; SI: m, IP: in.]
T_{gap(i)} = temperature of the gap between two glass panes (calculated using center of glazing thermal calculation algorithm, as described in ISO 15099 (ISO 2003). This value is first calculated using nondeflected state and is recalculated after the resulting deflection is calculated. This process is repeated until temperature at next iteration does not differ by more than 0.1 ºC
V_{gap(i)} = volume of the IGU gap in deflected state [Lentgh*Length*Length; SI: m^{3}, IP: in^{3}]
\[{V_{gap\left( i \right)}} = {V_{ini\left( i \right)}} + W\cdot H\cdot \left( {\overline {{L_{D,i}}}  \overline {{L_{D,i + 1}}} } \right)\]
Where,
\(\overline {{L_{D,i}}}\) is mean deflection value for ith pane. [Length; SI: m, IP: in.]
Deflection of each pane can be positive or negative and is done solely to establish reference. Current frame of reference is that positive deflection means that pane is deflecting towards left side, while negative deflection means that pane is deflecting towards right side . Whether the deflection is in the direction of reducing the gap width or increasing it, it will be the result of pressure difference, as described in . When pressure in the glazing unit is higher than surrounding environmental pressure, the deflection will be towards increasing gap width (i.e., ballooning), while the opposite situation will result in decreasing gap width (i.e., vacuuming)
Figure 111. Deflection Direction Convention
The important part of calculating deflection of the IGU is to determine mean deflection value for each glazing pane. Mean deflection value is used to calculate gap volume in deflected state . Mean deflection of glazing pane can be calculated by integrating :
\[\overline {{L_{D\left( i \right)}}} = \mathop \smallint \limits_{x = 0}^W \mathop \smallint \limits_{y = 0}^H \frac{{16\cdot \Delta {P_{\left( i \right)}}}}{{{\pi ^6}\cdot {D_{\left( i \right)}}}}\mathop \sum \limits_{m = 1,3,5 \ldots }^\infty \mathop \sum \limits_{n = 1,3,5 \ldots }^\infty \frac{{\sin \frac{{m\pi x}}{W}\sin \frac{{n\pi y}}{H}}}{{mn{{\left( {{{\left( {\frac{m}{W}} \right)}^2} + {{\left( {\frac{n}{H}} \right)}^2}} \right)}^2}}}\]
Which is identical with the following expression:
\[\overline {{L_{D\left( i \right)}}} = \frac{{16\cdot \Delta {P_{\left( i \right)}}}}{{{\pi ^6}\cdot {D_{\left( i \right)}}}}\mathop \sum \limits_{m = 1,3,5 \ldots }^\infty \mathop \sum \limits_{n = 1,3,5 \ldots }^\infty \frac{{\mathop \smallint \nolimits_{x = 0}^W \mathop \smallint \nolimits_{y = 0}^H \sin \frac{{m\pi x}}{W}\sin \frac{{n\pi y}}{H}}}{{mn{{\left( {{{\left( {\frac{m}{W}} \right)}^2} + {{\left( {\frac{n}{H}} \right)}^2}} \right)}^2}}}\]
and because integral of sin(x) is equal with –cos(x), above equation will become:
\[\overline {{L_{D\left( i \right)}}} = \frac{{16\cdot \Delta {P_{\left( i \right)}}}}{{{\pi ^6}\cdot {D_{\left( i \right)}}}}\mathop \sum \limits_{m = 1,3,5 \ldots }^\infty \mathop \sum \limits_{n = 1,3,5 \ldots }^\infty \frac{{\left( {1  cos\left( {m\pi } \right)} \right)\left( {1  cos\left( {n\pi } \right)} \right)}}{{mn{{\left( {{{\left( {\frac{m}{W}} \right)}^2} + {{\left( {\frac{n}{H}} \right)}^2}} \right)}^2}}}\]
Finally, because \(cos\left( {m\pi } \right)\) and \(cos\left( {n\pi } \right)\) values are always equal to 1 for the given range of m and n, above equation will became:
\[\overline {{L_{D\left( i \right)}}} = \frac{{16\cdot \Delta {P_{\left( i \right)}}}}{{{\pi ^6}\cdot {D_{\left( i \right)}}}}\mathop \sum \limits_{m = 1,3,5 \ldots }^\infty \mathop \sum \limits_{n = 1,3,5 \ldots }^\infty \frac{4}{{{m^2}{n^2}{\pi ^2}{{\left( {{{\left( {\frac{m}{W}} \right)}^2} + {{\left( {\frac{n}{H}} \right)}^2}} \right)}^2}}}\]
After calculating mean pane deflection the following equation is used to calculate mean gap width:
\[{L_{r\left( i \right)}} = L + \left( {\overline {{L_{D,\left( i \right)}}}  \overline {{L_{D,\left( {i + 1} \right)}}} } \right)\]
Where,
L_{r(i)} = Mean gap “i” width after incorporating glazing deflection. This mean gap width is used to recalculate thermal performance of deflected IGU.
\(\overline {{L_{D,i}}}\) * = mean glazing deflection for each pane “i”.*
Calculation of the deflection at the center of glazing and mean glazing deflection for each pane is an iterative process, where the initial temperature distribution is calculated for nondeflected state, then deflection is calculated based on this temperature distribution, new temperature distribution is calculated for this deflected state, then temperatures from previous iteration are compared to the current iteration and the process is repeated until the difference is no larger than 0.1 ºC.
At the end of calculations, program will calculate and return maximum deflection value for each pane (i.e., center of glazing deflection). If we label maximum deflection of each pane as L_{D(i),max}, we can calculate this value by substituting x = W/2 and y = H/2 in equation to determine deflection at the center point. Therefore,
\[{L_{D\left( i \right),max}} = \frac{{16\cdot \Delta {P_{\left( i \right)}}}}{{{\pi ^6}\cdot {D_{\left( i \right)}}}}\mathop \sum \limits_{m = 1,3,5 \ldots }^\infty \mathop \sum \limits_{n = 1,3,5 \ldots }^\infty \frac{{\sin \frac{{m\pi }}{2}\sin \frac{{n\pi }}{2}}}{{mn{{\left( {{{\left( {\frac{m}{2}} \right)}^2} + {{\left( {\frac{n}{2}} \right)}^2}} \right)}^2}}}\]
For glazing systems with more than two glazing layers, meaning multiple gas filled gaps, the deflection will be calculated for each glazing pane assuming that the pressure in a gap is independent from each other and calculated separately, unless spaces are “linked” together (e.g., stretched film middle glazing that has hole for equalizing pressure).
NonLinked Gaps in 3 or more glazing layer system:[LINK]
The procedure shown above generally applies to the 3 or more layer glazing system, with the exception that neighboring pressures are no longer P_{a}, but rather could be P_{a} on one side and P_{gap} on the other, or have P_{gap} on both sides, as shown in Figure 111 for gap “i”. Center of glazing thermal calculation will determine new temperature distribution, after deflection is calculated for each glazing and will be used to determine new P_{gap}, as per the procedure above.
Linked Gaps in 3 or more glazing layer system:[LINK]
When one or more gaps are linked together, their pressure is assumed to be identical (e.g., in triple glazing IGU P_{gap,1} = P_{gap,2}.) This pressure is calculated from temperatures of bounding glazing for linked gaps (e.g., for triple glazing IGU, glazing 1 and 3) and using neighboring pressures outside of those bounding glazing (e.g., for triple glazed IGU, P_{a} on both sides).
Note: This feature is not implemented in WINDOW 7.1. It is considered for future enhancements to the program.
Gap(s) Linked to Indoor or Outdoor Environment:[LINK]
If one or more glazing gaps are linked to either indoor or outdoor environment its pressure is fixed to P_{a}. In combination situations, such as two or more gaps linked together with one of them being linked to indoor or outdoor environment, they will all have fixed pressure of P_{a}.
When deflection is measured, it is normally measured at the point of maximum deflection. Maximum deflection occurs at center of the IGU (at W/2 and H/2).
Measured value is typically gap width at the point of maximum deflection, which we can label L_{G(i)}. For ith measured gap the width is equal to:
\[{L_{G\left( i \right)}} = {L_{\left( i \right)}} + \left( {{L_{D\left( i \right),max}}  {L_{D\left( {i + 1} \right),max}}} \right)\]
Figure 112. Sketch of the nonsymetrically Deflected Glazing Panes
If we label ratio of mean deflection and maximum deflection as R_{(i)}, then:
\[{R_{(i)}} = \frac{{\overline {{L_{D(i)}}} }}{{{L_{D(i)}}_{,\max }}} = \sum\limits_{m = 1,3,5,...}^\infty {\sum\limits_{n = 1,3,5,...}^\infty {\frac{4}{{{m^2}{n^2}{\pi ^2}{{\left( {{{(\frac{m}{W})}^2} + {{(\frac{n}{H})}^2}} \right)}^2}}}} }\]
\($ Important thing to note is that ratios (\){R_{( i )}}$ ) for all gaps in glazing system are equal.
\[{R_{\left( 1 \right)}} = {R_{\left( 2 \right)}} = \ldots = {R_{\left( i \right)}} = \ldots = {R_{\left( {n  1} \right)}} = R\]
Replacing and into the following equation is obtained:
\[{L_{r\left( i \right)}} = {L_{\left( i \right)}} + R\left( {{L_{D\left( i \right),max}}  {L_{D\left( {i + 1} \right),max}}} \right)\]
Combining with we get the following expression for the mean gap width:
\[{L_{r\left( i \right)}} = {L_{\left( i \right)}} + R\left( {{L_{G\left( i \right)}}  {L_{\left( i \right)}}} \right)\]
Number of equations given in expression is equal to n1, where n is number of panes. Therefore, we need one more equation to complete the system of equations that would allow us to solve for all independent variables. To get the last equation we can rewrite in slightly different manner:
\[{L_{D\left( i \right),max}} = \frac{{\Delta {P_{\left( i \right)}}}}{{{D_{\left( i \right)}}}}\cdot K\]
Where coefficient K combines all constant terms, while \({D_{\left( i \right)}}\) is given by and \(\Delta {P_{\left( i \right)}}\) is calculated by, and . Summing over all deflections, \({L_{D\left( i \right),max}}\) the following equation is obtained:
\[\mathop \sum \limits_{i = 1}^n \frac{{{D_{\left( i \right)}}}}{K}\cdot {L_{D\left( i \right),max}} = \mathop \sum \limits_{i = 1}^n \Delta {P_{\left( i \right)}} = 0\]
Note that sum of all \(\Delta {P_{\left( i \right)}}\) is equal to zero since outside pressure is equal to inside. Therefore, the remaining equation that completes the set of equations is:
\[\mathop \sum \limits_{i = 1}^n {D_{\left( i \right)}}\cdot {L_{D\left( i \right),max}} = 0\]
Solving system of equations[LINK]
In order to solve system of equations we will present in slightly different manner:
\[0 = {L_{\left( i \right)}}  {L_{G\left( i \right)}} + {L_{D\left( i \right),max}}  {L_{D\left( {i + 1} \right),max}}\]
Which in developed form will look like this:
\[0 = {L_{\left( 1 \right)}}  {L_{G\left( 1 \right)}} + {L_{D\left( 1 \right),max}}  {L_{D\left( 2 \right),max}}\]
\[0 = {L_{\left( 2 \right)}}  {L_{G\left( 2 \right)}} + {L_{D\left( 2 \right),max}}  {L_{D\left( 3 \right),max}}\]
\[0 = {L_{\left( i \right)}}  {L_{G\left( i \right)}} + {L_{D\left( i \right),max}}  {L_{D\left( {i + 1} \right),max}}\]
\[0 = {L_{\left( {n  1} \right)}}  {L_{G\left( {n  1} \right)}} + {L_{D\left( {n  1} \right),max}}  {L_{D\left( n \right),max}}\]
In order to express each \({L_{D\left( i \right),max}}\) as dependence from \({L_{D\left( n \right),max}}\) (deflection of inside/last pane) we will need to make sum from first to last, then from second to last, third to last and so on. This procedure will create following set of equations:
\[{L_{D\left( 1 \right),max}} = \mathop \sum \limits_{k = 1}^{k = n  1} \left( {{L_{G\left( k \right)}}  {L_{\left( k \right)}}} \right) + {L_{D\left( n \right),max}}\]
\[{L_{D\left( 2 \right),max}} = \mathop \sum \limits_{k = 2}^{k = n  1} \left( {{L_{G\left( k \right)}}  {L_{\left( k \right)}}} \right) + {L_{D\left( n \right),max}}\]
\[{L_{D\left( i \right),max}} = \mathop \sum \limits_{k = i}^{k = n  1} \left( {{L_{G\left( k \right)}}  {L_{\left( k \right)}}} \right) + {L_{D\left( n \right),max}}\]
\[{L_{D\left( {n  1} \right),max}} = \mathop \sum \limits_{k = n  1}^{k = n  1} \left( {{L_{G\left( k \right)}}  {L_{\left( k \right)}}} \right) + {L_{D\left( n \right),max}}\]
Now replacing this set of equations back to :
\[\mathop \sum \limits_{i = 1}^{n  1} {D_{\left( i \right)}}\cdot \left( {\mathop \sum \limits_{k = i}^{k = n  1} \left( {{L_{G\left( k \right)}}  {L_{\left( k \right)}}} \right) + {L_{D\left( n \right),max}}} \right) + {D_{\left( n \right)}}\cdot {L_{D\left( n \right),max}} = 0\]
Which solving by \({L_{D\left( n \right),max}}\) leads to the following equation:
\[{L_{D\left( n \right),max}} = \frac{{\mathop \sum \nolimits_{i = 1}^{n  1} ({D_{\left( i \right)}}\cdot \mathop \sum \nolimits_{k = i}^{k = n  1} \left( {{L_{\left( k \right)}}  {L_{G\left( k \right)}}} \right))}}{{\mathop \sum \nolimits_{i = 1}^n {D_{\left( i \right)}}}}\]
Calculating \({L_{D\left( n \right),max}}\) value from this equation and substituting it in will enable calculation of the deflection of remaining panes.
Arasteh, D.K., J.C. Kohler and B.T. Griffith. Draft 2009. Modeling Windows in EnergyPlus with only U, SHGC, and optionally VT. LBNL report. Full reference to be determined.
Arasteh, D. J. 2009. Modeling Windows in EnergyPlus with only U, SHGC, and optionally VT. Lawrence Berkeley National Laboratory.
Arasteh, D.K., M.S. Reilly and M.D. Rubin. 1989. A versatile procedure for calculating heat transfer through windows. American Society of Heating, Refrigeration and AirConditioning Engineers, ASHRAE Transactions, Vol. 95, Pt. 2.
Bernier, M., & Bourret, B. January 1997. “Effects of Glass Plate Curvature on the UFactor of Sealed Insulated Glazing Units”. Atlanta, GA: ASHRAE Transactions. Vol. 103, Pt. 1. American Society for Heating, Refrigerating and AirConditioning Engineers.
Collins, R., & FischerCripps, A. 1991. “Design of Support Pillar Arrays in Flat Evacuated Windows.”. Aust. J. Phys.
Collins, R., & Robinson, S. 1991. “Evacuated Glazing”. Solar Energy. Vol. 47, No. 1, pp. 2738.
Corruccini, R. (1959). “Gaseous Heat Conduction at Low Pressures and Temperatures”. Vacuum. Vol. 78, pp.1929.
Finlayson, E.U., D.K. Arasteh, C. Huizenga, M.D. Rubin and M.S. Reilly. 1993. WINDOW 4.0: documentation of calculation procedures. Lawrence Berkeley National Laboratory report no. LBL33943.
ISO. 2003. ISO 15099:2003. Thermal performance of windows, doors, and shading devices – Detailed calculations. International Organization for Standardization.
Klems, J. H. 1994A. “A New Method for Predicting the Solar Heat Gain of Complex Fenestration Systems: I. Overview and Derivation of the Matrix Layer Calculation.”. ASHRAE Transactions. 100 (pt.1): 10731086.
Klems, J. H. 1994B. “A New Method for Predicting the Solar Heat Gain of Complex Fenestration Systems: II. Detailed Description of the Matrix Layer Calculation.”. ASHRAE Transactions. 100 (pt.1): 10731086.
Klems, J. H. 1995. “Measurements of Bidirectional Optical Properties of Complex Shading Devices.”. ASHRAE Transactions. 101 (pt 1; Symposium Paper CH9581 (RP548)): 791801.
Klems, J. H. 1996. “A Comparison between Calculated and Measured SHGC for Complex Glazing Systems.”. ASHRAE Transactions. 102 (Pt. 1; Symposium Paper AT96161): 931939.
Klems, J. H. 1996. “Calorimetric Measurements of InwardFlowing Fraction for Complex Glazing and Shading Systems.”. ASHRAE Trans. 102(Pt. 1; Symposium Paper AT96163): 947954.
Munson, B.R, D.F. Young and T.H. Okiishi. 1998. “Fundamentals of Fluid Mechanics,” Third Edition Update, John Wiley & Sons, Inc.
Papamichael, K. J. 1998. “Determination and Application of Bidirectional SolarOptical Properties of Fenestration Systems.”. Cambridge, MA: 13th National Passive Solar Conference.
Simmler, H., U. Fischer and F. Winkelmann. 1996. SolarThermal Window Blind Model for DOE2. Lawrence Berkeley National Laboratory, Simulation Research Group internal report, (unpublished).
Timoshenko, S., & KreigerWoinowsky, S. 1959. “Theory of Plates and Shells” 2nd Edition. McGrawHill Company.
Equivalent Layer Fenestration Model[LINK]
The section describes the equivalent layer fenestration optical and thermal model. The Equivalent Layer fenestration model can have four types of attachments: drapes, venetian blinds, roller blinds and insect screens. In this model shading layers are assumed to be uniform and can be represented by an equivalent homogenous layer that has spatiallyaveraged “effective” optical and thermal properties (ASHRAE 1311RP). Likewise, venetian blinds can be characterized using effective optical and thermal properties. When solar radiation strikes a window surface some fraction of the incident solar radiation passes unobstructed through openings in a shading layer and the remaining fraction is intercepted by the structure of the layer. The intercepted radiation is partly absorbed, partly reflected and partly transmitted. These reflected and transmitted components of the scattered solar radiation are assumed to be uniformly diffuse. Shading layers, because of their openness, generally transmit longwave radiation, and the effective infrared properties of shades account for that. Using effective optical properties and a beam/diffuse split of solar radiation at each layer, the equivalent layer approach can represent multilayer systems. This representation provides virtually unlimited flexibility to combine different types of shading layers in a fenestration. The equivalent layer window model requires a few set of optical data to characterize a particular layer and this set of data is used to calculate effective layer properties. For instance, the effective solar optical properties of a venetian blind can be calculated as a function of slats optical properties and geometry. Also, it is possible to adjust slat angle at each time step in response to the changing angular position of the sun. Moreover, the model provides control strategies as a function of slat angle that can be changed at each time step as needed. Likewise, effective properties of a pleated drape are calculated as a function of fabric properties and a specified value of fullness. The only input data needed to fully characterize drapery fabrics, roller blinds and insect screens are material openness as area fraction, and the transmittance and reflectance at normal incidence. Shade openness area fraction is the same as the beambeam transmittance at normal incidence. In multilayer fenestration, each layer is separated by a gap. A gap in equivalent layer model is defined by specifying the fill gas and the gap spacing. Currently five gas types are allowed: Air, Argon, Xenon, Krypton and Custom. The convective heat transfer coefficient in a gap is calculated depending on the spacing, the temperatures of the layers and the fill gas properties. Equivalentlayer concept – offers wide range of multiple glazing and shading layers combination and can simulate multilayer complex fenestration systems. The effective layer properties of venetian blinds, pleated drapes, roller blinds, and insect screens are calculated from geometric layer models and material properties. A set of empirical correlations for estimating offnormal material properties were developed under ASHRAE research project (ASHRAE 1311RP).
The Equivalent Layer Analysis[LINK]
The equivalent layer windows system is treated as a series of parallel layers separated by gaps as shown in Figure 113. This multilayer structure has been used in several computer programs and the underlying theory has been documented (ASHRAE 1311RP).
Figure 113. Multilayer fenestration analysis structure (ASHRAE 1311RP)
Equivalent layer fenestration model uses twostep analysis. First, the flux of absorbed solar radiation at each layer, S_{i}, caused by the incident radiation flux, I_{sol}, is determined using net radiation analysis. Second, an energy balance is applied at each layer, accounting for heat transfer and the known set of absorbed solar radiations S_{i} values, in order to solve for the set of layer temperatures, T_{i}, and the corresponding heat fluxes. The fenestration model also accounts for the diathermanous shade layers in the longwave radiant exchange. The latter can be significant for shading layers. Glass is considered opaque with respect to longwave radiation.
Equivalent Layer Optical Model[LINK]
The multilayer optical model is based on an algorithm originally developed by Edwards (1977) and extended by Wright and Kotey (2006). The algorithm models the interaction of incident solar radiation with a glazing system composed of any number of parallel, planar, specular glazing layers. The shading layers scatter portion of the incident solar radiation diffusely, and the model tracks the beam and diffuse components of solar radiation as they interact with a multilayer system of glazing and shading layers. The conceptual arrangement for tracking beam and diffuse solar flux components is illustrated in Figure 114.
Analysis yields beambeam, beamdiffuse and diffusediffuse fluxes, providing full detail concerning the quantities of reflected, transmitted and absorbed radiation.
Figure 114. Solar analysis of the multilayer glazing/shading system showing beam and diffuse fluxes (ASHRAE 1311RP)
The beam and diffuse characterization of solar radiation demands an expanded set of solar optical properties (Wright and Kotey 2006). The quantities of interest for single layer are shown inFigure 115.
Where,
t = transmittance of a glazing or a shading layer, ()
r = reflectance of a glazing or a shading layer, ()
e = thermal emissivity of a glazing or a shading layer, ()
q = incident angle, (degree)
f = slat angle, (degree)
W = profile angle, (degree)
J = radiosity of a layer surface, (W/m2)
T = temperature of a layer, (C)
Subscripts
f = front side of a layer
b = back side of a layer
bb = beambeam optical property. Represents a fraction of the beam radiation incident at a given layer and leaves the layer without being scattered.
bd = beamdiffuse optical property. Represents a fraction of the beam radiation incident at a given layer and leaves the layer diffusely
dd = diffusediffuse optical property. Represents a diffuse radiation incident at a given layer and leaves the layer as diffuse
i = layer index, ()
Superscripts
m = represents material property (e.g., fabric material)
str = represents an apparent material structure property of roller blinds
w = represents an apparent wire property
Figure 115. Twelve solar properties assigned at each layer (ASHRAE 1311RP)
Each glazing or shading layer, in general, require eleven set of solar properties (Wright and Kotey 2006). Only beambeam properties are needed for common glazing layers. Each shading layer is represented by a few set of component properties and geometries that are converted to the full set of layer optical and thermal properties by equivalent layer models (ASHRAE 1311RP). The models also make account for offnormal incidence of solar radiation and can calculate the diffuse component from normal incidence values using numerical integration.
Equivalent Layer Thermal Model[LINK]
A surface energy balance is applied at each layer and the resulting set of equation is solved for layer temperatures and heat transfer fluxes. A schematic drawing of the multilayer illustration and variables designation is shown in Figure 116. The net radiation formulation based on the radiosities, Jf,i and Jb,i the radiant flux leaving the front and back surfaces of the i^{th} layer, respectively, is the used as a solution technique. The net radiant heat flux across a gap can be expressed as the difference between the radiosities of the bounding surfaces. The net radiation model analysis yields the layer temperatures and corresponding heat transfer coefficients that are used to determine Ufactor and SHGC. The thermal analysis is done in “ASHWAT_Thermal” routine.
Figure 116. Radiosity model used in thermal analysis of the multilayer glazing/shading system (ASHRAE 1311RP)
Glazing Layer Analysis[LINK]
The EquivalentLayer window model allows a wide range of options regarding the selection of glazing layers and gaps combinations. It also allows mixing glazing and shading layers in any sequence. For instance, it allows placing a venetian blind or a roller blind in between glazing layers. Only beambeam solar properties at normal incidence, infrared transmittance and infrared emissivity are required for common glazing types. Offnormal solar properties at a given sun position are estimated by adjusting the normal incidence values in “ASHWAT_OffNormalProperties” routine. The model assumes that the ratio between normal and offnormal transmittance is the same for the glazing layer in question and a reference piece of uncoated 6 mm glass with a moderate tint (ASHRAE 1311RP). A similar procedure is used to convert solar reflectance from normal to offnormal.
Gap Layer Analysis[LINK]
Gaps can be specified as sealed, vented indoors or vented outdoors. The sealed gap is considered as enclosed cavity. Vented gaps is specified for inner and outer most gaps only and when the gas type is AIR. The vented gap model assumes that air flow at the perimeter of the window is not restricted. Any fill gas can be specified by molecular mass and thermophysical properties of viscosity, specific heat and thermal conductivity. This is done by providing “a”, “b” and “c” coefficients of the quadratic equation of the form, p = a + bT + cT^{2}, where the T (K) is temperature and p is the property being evaluated. EnergyPlus has builtindata available for common fill gas components including Air, Argon, Krypton and Xenon. Users specify CUSTOM gas by defining the coefficients as an input.
Shade Layer Analysis[LINK]
Equivalent layer fenestration model has a complete set of solar and longwave models for the four shading layer types: drapes, venetian blinds, roller blinds and insect screens developed under ASHRAE 1311RP. And semiempirical models are used to evaluate the offnormal properties of drape, roller blind and insect screen materials. The effective layer properties of venetian blinds and the effect of pleating in drapes are determined using a more fundamental net radiation scheme. The openness fraction, A_{o} and beambeam solar transmittance at normal incidence, t_{bb} (q = 0) really represent a geometric quantity and it has been confirmed that they can be used interchangeably (ASHRAE 1311RP). Openness is simply the fraction of a material, by area, that is open. In equivalent layer fenestration model, beambeam transmittance at normal incidence (i.e., openness fraction), beamdiffuse transmittance and reflectance at normal incidence are required to characterize drapery fabric, roller blind and insect screen material. A conventional venetian blind in equivalent layer model can be characterized by specifying the geometry, solar reflectance and transmittance, and emissivity of the slats. The offnormal solar properties of drape, roller blind and insect screen materials were formulated based on measurements (ASHRAE 1311RP). The longwave properties of the of the drape fabrics, roller blinds and insect screens can be calculated from the material emissivity and openness of the fabric (ASHRE 1311RP). The optical model development for venetian blinds is presented by Yahoda and Wright (2004 and 2005), and Kotey et al. (2008).
Drapes and Curtains[LINK]
The offnormal optical properties of drapery fabric is determined from user specified optical properties at normal incidence (q = 0) using empirical correlations (Kotey et al., 2009a). The input data required to characterize drapery fabric are the beambeam transmittance of the material at normal incidence t_{bb} (q = 0) = Ao, the beamdiffuse transmittance of the material at normal incidencet_{bd} (q = 0), and the beamdiffuse (total) reflectance of the materialr_{bt} (q = 0). The offnormal properties are calculated as follows.
Offnormal Transmittance:
\[{\rm{\tau }}_{{\rm{bb}}}^{\rm{m}}\left( {\rm{\theta }} \right) = {\rm{\tau }}_{{\rm{bb}}}^{\rm{m}}{\rm{(\theta }} = {\rm{0)}} \cdot {\rm{co}}{{\rm{s}}^{\rm{b}}}\left( {\rm{\theta }} \right)\]
\[{\rm{\tau }}_{{\rm{bt}}}^{\rm{m}}\left( {\rm{\theta }} \right) = {\rm{\tau }}_{{\rm{bt}}}^{\rm{m}}{\rm{(\theta }} = {\rm{0)}} \cdot {\rm{co}}{{\rm{s}}^{\rm{b}}}\left( {\rm{\theta }} \right)\]
\[{\rm{b}} = {\left[ {\;  0.5 \cdot \ln \;\left( {\;{{\left[ {{\rm{\tau }}_{{\rm{bb}}}^{\rm{m}}{\rm{(\theta }} = {\rm{0)}}\;{\rm{,}}\;\;{\rm{0}}{\rm{.01}}} \right]}^{\;\max }}\;} \right)\;\;,\;\;0.35\;} \right]^{\;\max }}\]
\[{\rm{b}} = {\left[ {\;  0.5 \cdot \ln \;\left( {\;{{\left[ {{\rm{\tau }}_{{\rm{bt}}}^{\rm{m}}{\rm{(\theta }} = {\rm{0)}}\;{\rm{,}}\;\;{\rm{0}}{\rm{.01}}} \right]}^{\;\max }}\;} \right)\;\;,\;\;0.35\;} \right]^{\;\max }}\]
\[{\rm{\tau }}_{{\rm{bt}}}^{\rm{m}}\left( {\rm{\theta }} \right) = {\rm{\tau }}_{{\rm{bd}}}^{\rm{m}}\left( {\rm{\theta }} \right) + {\rm{\tau }}_{{\rm{bb}}}^{\rm{m}}\left( {\rm{\theta }} \right)\]
The offnormal reflectance:
\[{\rm{\rho }}_{{\rm{bt}}}^{\rm{m}}\left( {\rm{\theta }} \right) = {\rm{\rho }}_{{\rm{bd}}}^{\rm{m}}\left( {\rm{\theta }} \right) = {\rm{\rho }}_{{\rm{bt}}}^{\rm{m}}\left( {{\rm{\theta }} = {\rm{0}}} \right) + \left( {{\rm{\rho }}_{{\rm{bt}}}^{\rm{m}}\left( {{\rm{\theta }} = 9{\rm{0}}^\circ } \right)  {\rm{\rho }}_{{\rm{bt}}}^{\rm{m}}\left( {{\rm{\theta }} = {\rm{0}}} \right)} \right) \cdot \left( {1  {\rm{co}}{{\rm{s}}^{{\rm{0}}{\rm{.6}}}}\left( {\rm{\theta }} \right)} \right)\]
\[{\rm{\rho }}_{{\rm{bt}}}^{\rm{m}}\left( {{\rm{\theta }} = {\rm{90}}^\circ } \right) = {\rm{\rho }}_{{\rm{bt}}}^{\rm{m}}\left( {{\rm{\theta }} = {\rm{0}}} \right) + \left( {{\rm{1}}  {\rm{\rho }}_{{\rm{bt}}}^{\rm{m}}\left( {{\rm{\theta }} = {\rm{0}}} \right)} \right)\;\left( {\,0.7{{\left( {\;{\rm{\rho }}_{}^{\rm{y}}} \right)}^{\;0.7}}} \right)\]
The apparent yarn reflectance is given by:
\[{\rho ^y} = \frac{{\rho_{{\rm{bt}}}^{\rm{m}}\left( {{\rm{\theta }} = {\rm{0}}} \right)}}{{{\rm{1}}  \tau_{{\rm{bb}}}^{\rm{m}}\left( {{\rm{\theta }} = {\rm{0}}} \right)}} = \frac{{\rho_{{\rm{bt}}}^{\rm{m}}\left( {{\rm{\theta }} = {\rm{0}}} \right)}}{{{\rm{1}}  {{\rm{A}}_{\rm{o}}}}}\]
The above set of equations for drapery fabrics are used subject to the condition that the solar absorptance of the fabric, at normal incidence, is not less than 1% (ASHRE 1311RP). The diffusediffuse material properties, for Equivalent layer window model, are determined using Rhomberg integration with 32 panels covering the range from q = 0 to q = 90° (ASHRAE 1311RP). The subscript “X” stands for either front or back side of the layer.
\[\tau_{{\rm{x,dd}}}^{\rm{m}} = \;\;{\rm{2}}\;\int\limits_{{\rm{\theta }} = {\rm{0}}}^{{\raise0.5ex\hbox{\scriptstyle {\rm{\pi }}}\kern0.1em/\kern0.15em\lower0.25ex\hbox{\scriptstyle {\rm{2}}}}} {\tau_{{\rm{x,bt}}}^{\rm{m}}\left( {\rm{\theta }} \right) \cdot {\rm{cos}}\left( {\rm{\theta }} \right) \cdot {\rm{sin}}\left( {\rm{\theta }} \right)} \;{\rm{d\theta }}\]
\[\rho_{{\rm{x,dd}}}^{\rm{m}} = \;\;{\rm{2}}\;\int\limits_{{\rm{\theta }} = {\rm{0}}}^{{\raise0.5ex\hbox{\scriptstyle {\rm{\pi }}}\kern0.1em/\kern0.15em\lower0.25ex\hbox{\scriptstyle {\rm{2}}}}} {\rho_{{\rm{x,bt}}}^{\rm{m}}\left( {\rm{\theta }} \right) \cdot {\rm{cos}}\left( {\rm{\theta }} \right) \cdot {\rm{sin}}\left( {\rm{\theta }} \right)} \;{\rm{d\theta }}\]
The above set of equations for drapery fabrics apply to the full range of A_{o}, fabric transmittance and fabric reflectance including that falls within the bounds of Keyes’ (1967) fabric chart plus sheer fabrics (ASHRAE 1311RP). The longwave thermal emissivity and thermal transmittances of drapery fabric are calculated using the following correlations and fabric openness fraction (Kotey et al. 2008).
\[{\varepsilon ^{\rm{m}}} = {\rm{0}}{\rm{.87}} \cdot \left( {{\rm{1  }}{{\rm{A}}_{\rm{o}}}} \right)\]
\[{\tau ^{\rm{m}}} = {\rm{1}}  {\rm{0}}{\rm{.95}} \cdot \left( {{\rm{1}}  {{\rm{A}}_{\rm{o}}}} \right)\]
The optical and thermal properties determined using the above same sets of equations are equally valid for pleated drape shades (Kotey, et. al., 2009a). For pleated drape, the effective beambeam and beamdiffuse solar properties are determined by tracking both radiation components, for a given incident angle, and interaction with a fabric pleated rectangular geometry shown in Figure 117. The solar optical properties of the two different pleat planes are evaluated on the basis of the local solar incidence angle. Therefore, the effective layer properties are influenced not just by horizontal solar profile angle, W_{H}, but also by incidence angle (ASHRAE 1311RP).
Figure 117. Geometry used for Pleated Drape Analysis
The solar diffusediffuse and longwave effective properties of the pleated drape are evaluated using a much simpler netradiation analysis using conventional shape factors (Kotey, et. al., 2009a). Users can chose and apply the pleated drape model to any fabric and any degree of fullness (ASHRAE 1311RP).
Venetian Blinds[LINK]
The effective shortwave optical and longwave optical properties of venetian blind layer is estimated using analytical models (Yahoda and Wright 2004, 2005; Kotey et al. 2008). The model requires properties of venetian blind slats and geometry of the slats shown in Figure 1118.
Figure 1118. Likewise, the effective longwave properties are obtained for the layer knowing longwave properties of the slats.
Figure 118. Geometry and properties used for venetian blind analysis
The model assumes that venetian blind slats reflect and transmit solar radiation diffusely (ASHRAE 1311RP). The same assumption is made regarding thermal radiation. The input data required to characterize a venetian blind are: front and back side reflectance and transmittance of the slat, geometry and infrared emissivity and transmittance of the slate. The effective optical properties of the venetian blind are determined by tracking beam and diffuse solar radiation through various interactions with slats positioned at a given slat angle. The model uses simple foursurface model if the slats are fully sunlit and a sixsurface model if the slats are partially sunlit (ASHRAE 1311RP). Slats are assumed to be thin and flat but a correction is applied to account for slat curvature effect (Kotey et al. 2008).
Roller Blinds[LINK]
The offnormal properties of rollerblind are determined from solar properties of roller blind fabric measured at normal incidence (q = 0) using correlations (Kotey, et. al., 2009b). The offnormal properties for roller blind shades are calculated using the set equations given below:
\[{{\rm{\tau }}_{bb}}\left( {\rm{\theta }} \right) = \left\{ \begin{array}{l}{{\rm{\tau }}_{bb}}{\rm{(\theta }} = {\rm{0)}} \cdot {\rm{co}}{{\rm{s}}^{\rm{b}}}{\kern 1pt} \left( {\frac{{\rm{\theta }}}{{{{\rm{\theta }}_{{\rm{cutoff}}}}}} \cdot \frac{\pi }{2}} \right)\quad \theta < {\theta_{CutOff}}\\0\quad \quad \quad \theta \ge {\theta_{CutOff}}\end{array} \right.\]
\[{\rm{b}} = {\rm{0}}{\rm{.6}} \cdot {\rm{co}}{{\rm{s}}^{{\rm{0}}{\rm{.3}}}}\left( {{{\rm{A}}_{\rm{o}}} \cdot \frac{{\rm{\pi }}}{{\rm{2}}}} \right)\]
\[{{\rm{\theta }}_{{\rm{cutoff}}}} = {\rm{65}}^\circ + \;\left( {{\rm{90}}^\circ  {\rm{65}}^\circ } \right) \cdot \left( {{\rm{1}}  {\rm{cos}}\left( {{{\rm{A}}_{\rm{o}}} \cdot \frac{{\rm{\pi }}}{{\rm{2}}}} \right)} \right)\]
\[{\rm{\tau }}_{{\rm{bt}}}^{}\left( {\rm{\theta }} \right) = {\rm{\tau }}_{{\rm{bt}}}^{}{\rm{(\theta }} = {\rm{0)}} \cdot {\rm{co}}{{\rm{s}}^{\rm{b}}}{\kern 1pt} \left( {\rm{\theta }} \right)\]
\[b = \left\{
\begin{array}{lc}
0.133 \cdot \left( \rm{\tau}^{str} + 0.003 \right)^{0.467} &\text{if} \left( 0 \le \rm{\tau}^{str} \le 0.33 \right) \\
0.33 \cdot \left( 1  \rm{\tau}^{str} \right) &\text{if} \left( 0.33 < \rm{\tau}^{str} \le 1 \right)
\end{array}
\right.\]
\[{\rm{\tau }}_{}^{{\rm{str}}} = \frac{{{\rm{\tau }}_{{\rm{bt}}}^{}\left( {{\rm{\theta }} = {\rm{0}}} \right)  {{\rm{\tau }}_{{\rm{bb}}}}\left( {{\rm{\theta }} = {\rm{0}}} \right)}}{{{\rm{1}}  {{\rm{\tau }}_{{\rm{bb}}}}\left( {{\rm{\theta }} = {\rm{0}}} \right)}}\]
\[{\rm{\tau bt}}\left( {\rm{\theta }} \right) = {\rm{\tau bd}}\left( {\rm{\theta }} \right) + {\rm{\tau bb}}\left( {\rm{\theta }} \right)\]
The offnormal solar property calculation of roller blind is based on a set of correlations developed from measurement data using samples of commonly used commercially produced roller blind material openness range of 0.0 – 0.14. Thus, these correlations are not applicable for shades with materials openness fraction exceeding 0.20. The mean solar reflectance of a roller blind material was found to be purely diffuse and unaffected by incidence angle and is given by:
\[{\rm{\rho }}_{{\rm{bt}}}^{}\left( {{\rm{\theta }} = {\rm{0}}} \right) = {\rm{\rho }}_{{\rm{bd}}}^{}\left( {{\rm{\theta }} = {\rm{0}}} \right) = {\rm{\rho }}_{{\rm{bt}}}^{}\left( {\rm{\theta }} \right) = {\rm{\rho }}_{{\rm{bd}}}^{}\left( {\rm{\theta }} \right)\]
The diffusediffuse transmittance and reflectance are obtained by Rhomberg numerical integration. The longwave properties of roller blind material determined using the material property and the openness fraction (Kotey et al. 2008) as shown below:
\[{\varepsilon ^{{\rm{rb}}}} = {\rm{0}}{\rm{.91}} \cdot \left( {{\rm{1  }}{{\rm{A}}_{\rm{o}}}} \right)\]
\[{\tau ^{{\rm{rb}}}} = {\rm{1}}  {\rm{0}}{\rm{.95}} \cdot \left( {{\rm{1}}  {{\rm{A}}_{\rm{o}}}} \right)\]
Insect Screens[LINK]
The empirical correlations formulated to obtain the effective offnormal solar and longwave properties of insect screens were based on measurements (Kotey et al. (2009a). Insect screen geometry is shown in Figure 119. The calculation of effective solar properties requires a set of properties measured at normal incidence: t_{bb} (q = 0), t_{bd} (q = 0) and r_{bt} (q = 0).
Figure 119. Geometry used for insect screen analysis
Openness can be determined by optical measurement at normal incidence, A_{o} = t_{bb} (q = 0), but in the case of insect screens A_{o} can reliably be calculated knowing wire diameter (d), and wire spacing (s) as follows:
\[{A_o} = {\left( {1  {d \mathord{\left/ {\vphantom {d s}} \right. } s}} \right)^2}\]
The incidence angle beyond which direct beam transmission is cut off, q_{CutOff}, can also be estimated purely from geometry and is given by:
\[{\theta_{CutOff}} = co{s^{  1}}\left( {{d \mathord{\left/ {\vphantom {d S}} \right. } S}} \right)\]
The offnormal properties are calculated as follows.
\[{\tau_{bb}}\left( {\rm{\theta }} \right) = \left\{ \begin{array}{l}{\tau_{bb}}{\rm{(\theta }} = {\rm{0)}} \cdot {\rm{co}}{{\rm{s}}^{\rm{b}}}\left( {\frac{{\rm{\theta }}}{{{{\rm{\theta }}_{{\rm{cutoff}}}}}} \cdot \frac{\pi }{{\rm{2}}}} \right)\quad {\rm{\theta }} < {{\rm{\theta }}_{{\rm{cutoff}}}}\\0\quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \;\,{\rm{\theta }} \ge {{\rm{\theta }}_{{\rm{cutoff}}}}\end{array} \right.\]
\[{\rm{b}} =  {\rm{0}}{\rm{.45}}\;{\rm{ln}}\;\left( {{\rm{MAX}}\left( {{{\rm{\tau }}_{{\rm{bb}}}}\left( {{\rm{\theta }} = {\rm{0}}} \right){\rm{,}}\;{\rm{0}}{\rm{.01}}} \right)} \right)\; + \;{\rm{0}}{\rm{.1}}\]
\[{\rm{\tau }}_{{\rm{bt}}}^{}\left( {\rm{\theta }} \right) = {\rm{\tau }}_{{\rm{bt}}}^{}{\rm{(\theta }} = {\rm{0)}} \cdot {\rm{co}}{{\rm{s}}^{\rm{b}}}\left( {\rm{\theta }} \right)\]
\[{\rm{b}} =  {\rm{0}}{\rm{.65}}\;{\rm{ln}}\;\left( {{\rm{MAX}}\left( {{{\rm{\tau }}_{{\rm{bt}}}}\left( {{\rm{\theta }} = {\rm{0}}} \right){\rm{,}}\;{\rm{0}}{\rm{.01}}} \right)} \right)\; + \;{\rm{0}}{\rm{.1}}\]
\[{\rm{\tau }}_{{\rm{bd}}}^{}\left( {\rm{\theta }} \right) = {\rm{\tau }}_{{\rm{bt}}}^{}\left( {\rm{\theta }} \right)  {\rm{\tau }}_{{\rm{bb}}}^{}\left( {\rm{\theta }} \right)\]
\[{\rm{\rho }}_{{\rm{bt}}}^{}\left( {\rm{\theta }} \right) = {\rm{\rho }}_{{\rm{bt}}}^{}\left( {{\rm{\theta }} = {\rm{0}}} \right) + \left( {{\rm{\rho }}_{{\rm{bt}}}^{}\left( {{\rm{\theta }} = 9{\rm{0}}^\circ } \right)  {\rm{\rho }}_{{\rm{bt}}}^{}\left( {{\rm{\theta }} = {\rm{0}}} \right)} \right) \cdot \left( {1  {\rm{co}}{{\rm{s}}^{\rm{b}}}\left( {\rm{\theta }} \right)} \right)\]
\[{\rm{b}} =  {\rm{0}}{\rm{.45}}\;{\rm{ln}}\;{\left( {{{\rm{\rho }}^{\rm{w}}}} \right)^{}}\]
\[{\rm{\rho }}_{}^{\rm{w}} = \frac{{{\rm{\rho }}_{{\rm{bt}}}^{}\left( {{\rm{\theta }} = {\rm{0}}} \right)}}{{{\rm{1}}  {\rm{\tau }}_{{\rm{bb}}}^{}\left( {{\rm{\theta }} = {\rm{0}}} \right)}}\]
\[{\rho_{{\rm{bt}}}}\left( {{\rm{\theta }} = {\rm{90}}^\circ } \right) = {\rho_{{\rm{bt}}}}\left( {{\rm{\theta }} = {\rm{0}}} \right) + \left( {1  {\rho_{{\rm{bt}}}}\left( {{\rm{\theta }} = {\rm{0}}} \right)} \right)\;\left( {{\rm{0}}{\rm{.35}}{\rho ^{\rm{w}}}} \right)\]
The diffusediffuse material properties are obtained by Rhomberg numerical integration. The longwave properties of insect screen are given by expressions similar those formulated for drapery fabrics, and roller blinds (Kotey et al. 2008) is given by:
\[{\varepsilon ^{screen}} = {\varepsilon ^{\rm{w}}} \cdot \left( {{\rm{1  }}{{\rm{A}}_{\rm{o}}}} \right)\]
\[{\tau ^{{\rm{screen}}}} = {\rm{1  }}{\tau ^{\rm{w}}} \cdot \left( {{\rm{1  }}{{\rm{A}}_{\rm{o}}}} \right)\]
The apparent wire material emissivity can be taken as e^{w} = 0.93 for common insect screens  screens with dark, rough, nonmetallic wire surfaces. The corresponding infrared wire material transmittance is t^{w} = 0.98. A lower infrared emissivity can be used for screens constructed with shiny metallic wire. For example, to model stainlesssteel wire mesh use e^{w} = 0.32 and t^{w} = 0.81.
Integration with Heat Balance Method[LINK]
The solar model calculations are performed surface by surface following the existing structure of EnergyPlus’s heat balance algorithm. The optical and thermal calculations are performed for each surface at each time step. The thermal model needs to be updated at each time step during the surface heat balance iteration. Thus, window thermal model is invoked during the surface heat balance but only once. This is consistent with the current EnergyPlus’s window model inside surface heat balance iteration scheme. The equivalent layer window thermal model also has internal iteration scheme. The thermal model routine also requires the solar and shortwave radiation flux absorbed by each layer of a fenestration at each time step. The calculation of the absorbed radiation flux is performed using the existing scheme except that the optical properties are calculated by the equivalent layer optical model. The equivalent layer thermal model returns the temperatures and the fluxes at each layer.
Equivalent Layer Window Solar Model[LINK]
The equivalent layer window solar model calculates the transmittance of the window and absorptance of each layer. Separate optical properties calculations are performed for exterior beam and diffuse radiations, and interior diffuse radiation (e.g., lighting and interreflected solar gain). The fractions for beam radiations depend on solar incident angle and hence updated for each time step. Diffuse radiation fractions are the same for all time steps unless the shade characteristics are altered; for instance, when the venetian blind slat angle is controlled. Otherwise, the diffuse fractions are constant for a given configuration and time step. The hub of the optical properties calculation routine for Equivalent Layer window model is the routine “ASHWAT_Solar”. For beam radiation, at each time step, first the normal incidence optical properties are modified for current incident angle. Using the offnormal properties, then the ASHWAT_Solar routine sets up the coefficient matrix based on the net radiation concept to determine the effective absorptance of each layer and the window transmittance. The total intensity of shortwave radiation absorbed at each layer is determined by multiplying the incident solar radiation, reflected internal solar radiations, and internal shortwave source components by the appropriate absorbed fractions and summing them.
Equivalent Layer Window Thermal Model[LINK]
The equivalent layer thermal model is calculated only once for each time step. But the thermal model has internal iterative solutions scheme. During each time step, the procedure initially assumes that room air and means radiant temperatures are known. The surface heat balance loops through all zone surfaces while invoking the Equivalent Layer thermal model only once for each window surface at the first iteration. The window surfaces temperatures from the first iteration are used to complete the heat balances for the indoor (and implicitly the outdoor) face of each surface iteratively. Once indoor surface temperatures are calculated using the surface heat balance, the zone air temperature can be updated and the loads are predicted. Shaded fenestration in general do not have single inside temperature by virture their longwave radiation transmittance. The equivalent layer window model accounts for this using effective emissivity of the composite layers derived for each fenestration (ASHRAE 1311RP) as shown below:
\[{\varepsilon ^{eff}} = \sum\limits_{j = 0}^{nl} {\left[ {{\varepsilon_j} \cdot \prod\limits_{k = j{ + 1}}^{nl + {1}} {{\tau_k}} } \right]}\]
where
e^{eff} = composite indoor (roomside) longwave emissivity, ()
e_{j} = effective emissivity of layer j, ()
t_{k} = effective infrared transmittance of layer k (t_{nl+1} = 1)
nl = number of layers in fenestration system (glazing and shade). Layers are numbered outside to inside (layer 1 is outermost, layer nl is innermost).
Each equivalent layer window surafce yields net longwave radiant flux exchanged with the zone surfaces. Net longwave radiation exchange from the window to the zone is recast for a composite surface temperature calculation as follows:
\[{T^{eff}} = \sqrt[4]{{\frac{{{Q_{lw}}}}{{\sigma {\varepsilon ^{eff}}}}}} + {T_0}\]
where
T^{eff} = inside surface temperature of the composite layer, °C (°F)
Q_{lw} = infrared radiant gain from zone, W/m2 (Btu/hft2)
s = StefanBoltzmann constant, W/m2K4 (Btu/hft2R4)
T_{0} = Temperature of absolute zero, 273.15 °C (459.67 °F)
Recalculating effective inside face window temperature may result in extra convection flux. The “extra” (Other) convective flux is computed; this is the gain in excess of that resulting from the standard surface heat balance model. The net other convection term QX calculated below is added to the zone air heat balance (ASHRAE 1311RP).
\[QX = {Q_{conv}}  hc \cdot ({T^{eff}}  {T_a})\]
Where,
h_{c} = inside convective coefficient of the fenestration, W/m^{2}K (Btu/hft^{2}F)
Q_{conv} = total convective heat flux to zone from equivalent layer window thermal model, W/m^{2} (Btu/hft^{2}); includes openchannel gains and impact of inside surface convective coefficient
T_{a} = zone air temperature, °C (°F)
Edwards, D.K. 1977. Solar absorption by each element in an absorbercoverglass array, Technical Note, Solar Energy, Vol. 19, pp. 401402.
Parmelee, G. V., and W. W. Aubele. 1952. The shading of sunlit glass: an analysis of the effect of uniformly spaced flat opaque slats, ASHVE Transactions, Vol. 58, pp. 377398.
Farber, Erich A.; William A. Smith, C.W. Pennington, John C. Reed. 1963. Theoretical analysis of solar heat gain through insulating glass with inside shading. ASHRAE Transactions, Vol. 69, pp.393405.
Rheault, S., and E. Bilgen. 1989. Heat transfer analysis in an automated venetian blind system, Journal of Solar Energy, Vol. 111 (Feb.), pp. 8995.
Pfrommer, P., K. J. Lomas, and C. Kupke. 1996. “Solar Radiation Transport through SlatType Blinds: a New Model and its Application for Thermal Simulation of Buildings,” Solar Energy, Vol. 57, No. 2, pp. 7791.
Rosenfeld, J.L.J., W. J. Platzer, H. Van Dijk, and A. Maccari. 2000. “Modelling the Optical and Thermal Properties of Complex Glazing: Overview of Recent Developments”, Solar Energy, Vol. 69 Supplement, No. 16, pp.113.
Yahoda, D. S. and J. L. Wright. 2004. “Methods for Calculating the Effective Longwave Radiative Properties of a Venetian Blind Layer,” ASHRAE Transactions, Vol. 110, Pt. 1., pp. 463473.
Yahoda, D. S. and J. L. Wright. 2005. “Methods for Calculating the Effective SolarOptical Properties of a Venetian Blind Layer,” ASHRAE Transactions, Vol. 111, Pt. 1, pp. 572586.
Yahoda, D. S. and J. L. Wright. 2004. “Heat Transfer Analysis of a BetweenPanes Venetian Blind Using Effective Longwave Radiative Properties,” ASHRAE Transactions, Vol. 110, Pt. 1., pp. 455462.
Huang, N.Y.T., J. L. Wright, M. R. Collins. 2006. “Thermal Resistance of a Window with an Enclosed Venetian Blind: Guarded Heater Plate Measurements,” ASHRAE Transactions, Vol. 112, Pt. 2. pp. 1321.
Wright, J. L. 2008. “Calculating CentreGlass Performance Indices of Glazing Systems with Shading Devices,” ASHRAE Transactions, Vol. 114, Pt. 2.
Wright, J. L., N. Y. T. Huang, and M. R. Collins. 2008. “Thermal Resistance of a Window with an Enclosed Venetian Blind: A Simplified Model,” ASHRAE Transactions, Vol. 114, Pt. 1.
Kotey, N. A., J. L. Wright, and M. R. Collins. 2008. “Determining Longwave Radiative Properties of Flat Shading Materials,” 33rd Annual SESCI / 3rd CSBC Conference Proceedings, Fredericton, NB.
Kotey, N.A., Wright, J.L., M. R. Collins. 2009a. “Determination of AngleDependent Solar Optical Properties of Roller Blind Materials,” drafted for submission to ASHRAE Transactions, Vol. 115, Pt. 1.
Kotey, N.A., Wright, J.L., M. R. Collins. 2009b. “Determination of AngleDependent Solar Optical Properties of Drapery Fabrics,” in review, ASHRAE Transactions, Vol. 115, Pt. 2.
Wright, John L., Charles S. Barnaby, Michael R. Collins, and Nathan A. Kotey. Improving Cooling Load Calculations for Fenestration with Shading Devices . ASHRAE 1311RP, Final Report, February 11, 2009.
Window Heat Balance Calculation[LINK]
Table 36. Fortran Variables used in Window Heat Balance CalculationsThe Glazing Heat Balance Equations[LINK]
The window glass face temperatures are determined by solving the heat balance equations on each face every time step. For a window with N glass layers there are 2N faces and therefore 2N equations to solve. Figure 95 shows the variables used for double glazing (N = 2).
Figure 95. Glazing system with two glass layers showing variables used in heat balance equations.
The following assumptions are made in deriving the heat balance equations:
1) The glass layers are thin enough (a few millimeters) that heat storage in the glass can be neglected; therefore, there are no heat capacity terms in the equations.
2) The heat flow is perpendicular to the glass faces and is one dimensional. See “Edge of Glass Corrections,” below, for adjustments to the gap conduction in multipane glazing to account for 2D conduction effects across the pane separators at the boundaries of the glazing.
3) The glass layers are opaque to IR. This is true for most glass products. For thin plastic suspended films this is not a good assumption, so the heat balance equations would have to be modified to handle this case.
4) The glass faces are isothermal. This is generally a good assumption since glass conductivity is very high.
5) The short wave radiation absorbed in a glass layer can be apportioned equally to the two faces of the layer.
The four equations for doubleglazing are as follows. The equations for single glazing (N = 1) and for N = 3 and N = 4 are analogous and are not shown.
\[{E_o}{\varepsilon_1}  {\varepsilon_1}\sigma \theta_1^4 + {k_1}\left( {{\theta_2}  {\theta_1}} \right) + {h_o}\left( {{T_o}  {\theta_1}} \right) + {S_1} = 0\]
\[{k_1}\left( {{\theta_1}  {\theta_2}} \right) + {h_1}\left( {{\theta_3}  {\theta_2}} \right) + \sigma \frac{{{\varepsilon_2}{\varepsilon_3}}}{{1  \left( {1  {\varepsilon_2}} \right)\left( {1  {\varepsilon_3}} \right)}}\left( {\theta_3^4  \theta_2^4} \right) + {S_2} = 0\]
\[{h_1}({\theta_2}  {\theta_3}) + {k_2}({\theta_4}  {\theta_3}) + \sigma \frac{{{\varepsilon_2}{\varepsilon_3}}}{{1  \left( {1  {\varepsilon_2}} \right)\left( {1  {\varepsilon_3}} \right)}}\left( {\theta_2^4  \theta_3^4} \right) + {S_3} = 0\]
\[{E_i}{\varepsilon_4}  {\varepsilon_4}\sigma \theta_4^4 + {k_2}\left( {{\theta_3}  {\theta_4}} \right) + {h_i}\left( {{T_i}  {\theta_4}} \right) + {S_4} = 0\]
Absorbed Radiation[LINK]
S_{i} in Equations to is the radiation (shortwave and longwave from zone lights and equipment) absorbed on the i^{th} face. Shortwave radiation (solar and shortwave from lights) is assumed to be absorbed uniformly along a glass layer, so for the purposes of the heat balance calculation it is split equally between the two faces of a layer. Glass layers are assumed to be opaque to IR so that the thermal radiation from lights and equipment is assigned only to the inside (roomside) face of the inside glass layer. For N glass layers S_{i} is given by
\[{S_{2j  1}} = {S_{2j}} = \frac{1}{2}\left( {I_{bm}^{ext}\cos \phi A_j^f(\phi ) + I_{dif}^{ext}A_j^{f,dif} + I_{sw}^{{\mathop{\rm int}} }A_j^{b,dif}} \right),{\rm{ }}j = 1{\rm{ to }}N\]
\[{S_{2N}} = {S_{2N}} + {\varepsilon_{2N}}I_{lw}^{{\mathop{\rm int}} }\]
Here
\(I_{bm}^{ext}\) = exterior beam normal solar irradiance
\(I_{dif}^{ext}\) = exterior diffuse solar incident on glazing from outside
\(I_{sw}^{{\mathop{\rm int}} }\) = interior shortwave radiation (from lights and from reflected diffuse solar) incident on glazing from inside
\(I_{lw}^{{\mathop{\rm int}} }\) = longwave radiation from lights and equipment incident on glazing from inside
\({\varepsilon_{2N}}\) = emissivity (thermal absorptance) of the roomside face of the inside glass layer
RoomSide Convection[LINK]
The correlation for roomside convection coefficient, \({h_i}\) , is from ISO 15099 section 8.3.2.2. (Prior to EnergyPlus version 3.1, the value for \({h_i}\) was modeled using the “Detailed” algorithm for opaque surface heat transfer, e.g. for a vertical surface \({h_i} = 1.31{\left {\Delta T} \right^{{\raise0.7ex\hbox{1} \!\mathord{\left/ {\vphantom {1 3}}\right.}\!\lower0.7ex\hbox{3}}}}\) ; see section Detailed Natural Convection Algorithm). The ISO 15099 correlation is for still room air and is determined in terms of the Nusselt number, \(Nu\) , where
\[{h_i} = Nu\left( {\frac{\lambda }{H}} \right)\]
where,
\(\lambda\) is the thermal conductivity of air, and
\(H\) is the height of the window.
The Rayleigh number based on height, \(R{a_H}\) , is calculated using,
\[R{a_H} = \frac{{{\rho ^2}{H^3}g\,{c_p}\left {{T_{surf,i}}  {T_{air}}} \right}}{{{T_{m,f}}\mu \lambda }}\]
where,
\(\rho\) is the density of air
\(g\) is the acceleration due to gravity,
\({c_p}\) is the specific heat of air,
\(\mu\) is the dynamic viscosity of air, and
\({T_{m,f}}\) is the mean film temperature in Kelvin given by,
\[{T_{m,f}} = {T_{air}} + \frac{1}{4}\left( {{T_{surf,i}}  {T_{air}}} \right)\]
There are four cases for the Nusselt correlation that vary by the tilt angle in degrees, \(\gamma\) , and are based on heating conditions. For cooling conditions (where \({T_{surf,i}} > {T_{air}}\) ) the tilt angle is complemented so that \(\gamma = 180  \gamma\)
Case A. \({0^\circ } \le \gamma < 15^\circ\)
\[Nu = 0.13Ra_H^{{\raise0.7ex\hbox{1} \!\mathord{\left/ {\vphantom {1 3}}\right.}\!\lower0.7ex\hbox{3}}}\]
Case B. \(15^\circ \le \gamma \le 90^\circ\)
\[R{a_{cv}} = 2.5 \times {10^5}{\left( {\frac{{{e^{0.72\gamma }}}}{{\sin \lambda }}} \right)^{{\raise0.7ex\hbox{1} \!\mathord{\left/ {\vphantom {1 5}}\right.}\!\lower0.7ex\hbox{5}}}}\]
\[Nu = 0.56{\left( {R{a_H}\sin \gamma } \right)^{{\raise0.7ex\hbox{1} \!\mathord{\left/ {\vphantom {1 4}}\right.}\!\lower0.7ex\hbox{4}}}};\;for\;R{a_H} \le R{a_{CV}}\]
\[Nu = 0.13\left( {Ra_H^{{\raise0.7ex\hbox{1} \!\mathord{\left/ {\vphantom {1 3}}\right.}\!\lower0.7ex\hbox{3}}}  Ra_{CV}^{{\raise0.7ex\hbox{1} \!\mathord{\left/ {\vphantom {1 3}}\right.}\!\lower0.7ex\hbox{3}}}} \right) + 0.56{\left( {R{a_{CV}}\sin \gamma } \right)^{{\raise0.7ex\hbox{1} \!\mathord{\left/ {\vphantom {1 4}}\right.}\!\lower0.7ex\hbox{4}}}};\;R{a_H} > R{a_{CV}}\]
Case C. \(90^\circ < \gamma \le 179^\circ\)
\[Nu = 0.56{\left( {R{a_H}\sin \gamma } \right)^{{\raise0.7ex\hbox{1} \!\mathord{\left/ {\vphantom {1 4}}\right.}\!\lower0.7ex\hbox{4}}}};\;{10^5} \le R{a_H}\sin \gamma < {10^{11}}\]
Case D. \(179^\circ < \gamma \le 180^\circ\)
\[Nu = 0.58Ra_H^{{\raise0.7ex\hbox{1} \!\mathord{\left/ {\vphantom {1 5}}\right.}\!\lower0.7ex\hbox{5}}};\;R{a_H} \le {10^{11}}\]
The material properties are evaluated at the mean film temperature. Standard EnergyPlus pyschrometric functions are used for \(\rho\) and \({c_p}\) . Thermal conductivity is calculated using,
\(\lambda = 2.873 \times {10^{  3}} + 7.76 \times {10^{  5}}{T_{m,f}}\) .
Kinematic viscosity is calculated using,
\(\mu = 3.723 \times {10^{  6}} + 4.94 \times {10^{  8}}{T_{m,f}}\) .
This correlation depends on the surface temperature of the roomside glazing surface and is therefore included inside the window heat balance interation loop.
Solving the Glazing Heat Balance Equations[LINK]
The equations are solved as follows:
1) Linearize the equations by defining \({h_{r,i}} = {\varepsilon_i}\sigma \theta_i^3\) . For example, Equation becomes
\[{E_o}{\varepsilon_1}  {h_{r,1}}{\theta_1} + {k_1}\left( {{\theta_2}  {\theta_1}} \right) + {h_o}\left( {{T_o}  {\theta_1}} \right) + {S_1} = 0\]
2) Write the equations in the matrix form \(A\theta = B\)
3) Use previous time step’s values of \({\theta_i}\) as initial values for the current time step. For the first time step of a design day or run period the initial values are estimated by treating the layers as a simple RC network.
4) Save the \({\theta_i}\) for use in the next iteration: \({\theta_{prev,i}} = {\theta_i}\)
5) Using \({\theta_{2N}}\) , reevaluate the roomside face surface convection coefficient \({h_i}\)
6) Using the \({\theta_i}\) to evaluate the radiative conductances \({h_{r,i}}\)
7) Find the solution \(\theta = {A^{  1}}B\) by LU decomposition
8) Perform relaxation on the the new \({\theta_i}\) : \({\theta_i} \to \left( {{\theta_i} + {\theta_{prev,i}}} \right)/2\)
9) Go to step 4
Repeat steps 4 to 9 until the difference, \(\Delta {\theta_i}\) , between values of the \({\theta_i}\) in successive iterations is less than some tolerance value. Currently, the test is
\[\frac{1}{{2N}}\sum\limits_{i = 1}^{2N} {\Delta {\theta_i}{\rm{ }} < 0.02K}\]
If this test does not pass after 100 iterations, the tolerance is increased to 0.2K. If the test still fails the program stops and an error message is issued.
The value of the inside face temperature, \({\theta_{2N}}\) , determined in this way participates in the zone heat balance solution (see Outdoor/Exterior Convection) and thermal comfort calculation (see Occupant Thermal Comfort).
EdgeOfGlass Effects[LINK]
Table 37. Fortran Variables used in Edge of Glass calculationsBecause of thermal bridging across the spacer separating the glass layers in multipane glazing, the conductance of the glazing near the frame and divider, where the spacers are located, is higher than it is in the center of the glass. The areaweighted net conductance (without inside and outside air films) of the glazing in this case can be written
\[\overline h = \left( {{A_{cg}}{h_{cg}} + {A_{fe}}{h_{fe}} + {A_{de}}{h_{de}}} \right)/{A_{tot}}\]
where
h_{cg} = conductance of centerofglass region (without air films)
h_{fe} = conductance of frame edge region (without air films)
h_{de} = conductance of divider edge region (without air films)
A_{cg} = area of centerofglass region
A_{fe} = area of frame edge region
A_{de} = area of divider edge region
A_{tot} = total glazing area = \({A_{cg}} + {A_{fe}} + {A_{de}}\)
The different regions are shown in Figure 96:
Figure 96: Different types of glass regions.
Equation can be rewritten as
\[\overline h = {h_{cg}}\left( {{\eta_{cg}} + {\alpha_{fe}}{\eta_{fe}} + {\alpha_{de}}{\eta_{de}}} \right)\]
where
\[{\eta_{cg}} = {A_{cg}}/{A_{tot}}\]
\[{\eta_{fe}} = {A_{fe}}/{A_{tot}}\]
\[{\eta_{de}} = {A_{de}}/{A_{tot}}\]
\[{\alpha_{fe}} = {h_{fe}}/{h_{cg}}\]
\[{\alpha_{de}} = {h_{de}}/{h_{cg}}\]
The conductance ratios \({\alpha_{fe}}\) and \({\alpha_{de}}\) are user inputs obtained from Window 5. They depend on the glazing construction as well as the spacer type, gap width, and frame and divider type.
In the EnergyPlus glazing heat balance calculation effective gap convective conductances are used to account for the edgeofglass effects. These effective conductances are determined as follows for the case with two gaps (triple glazing). The approach for other numbers of gaps is analogous.
Neglecting the very small resistance of the glass layers, the centerofglass conductance (without inside and outside air films) can be written as
\[{h_{cg}} = {\left( {{{\left( {{h_{r,1}} + {h_{c,1}}} \right)}^{  1}} + {{\left( {{h_{r,2}} + {h_{c,2}}} \right)}^{  1}}} \right)^{  1}}\]
where
${h_{c,k}} = $ convective conductance of the k^{th} gap
${h_{r,k}} = $ radiative conductance of the k^{th} gap
\[ = \frac{1}{2}\sigma \frac{{{\varepsilon_i}{\varepsilon_j}}}{{1  \left( {1  {\varepsilon_i}} \right)\left( {1  {\varepsilon_j}} \right)}}{\left( {{\theta_i} + {\theta_j}} \right)^3}\]
${_i},{_j} = $ emissivity of the faces bounding the gap
${_i},{_j} = $ temperature of faces bounding the gap (K)
Equation then becomes
\[\overline h = \left( {{\eta_{cg}} + {\alpha_{fe}}{\eta_{fe}} + {\alpha_{de}}{\eta_{de}}} \right){\left( {{{\left( {{h_{r,1}} + {h_{c,1}}} \right)}^{  1}} + {{\left( {{h_{r,2}} + {h_{c,2}}} \right)}^{  1}}} \right)^{  1}}\]
We can also write \(\overline h\) in terms of effective convective conductances of the gaps as
\[\overline h = {\left( {{{\left( {{h_{r,1}} + {{\overline h }_{c,1}}} \right)}^{  1}} + {{\left( {{h_{r,2}} + {{\overline h }_{c,2}}} \right)}^{  1}}} \right)^{  1}}\]
Comparing Eqs. and we obtain
\[{h_{r,k}} + {\overline h_{c,k}} = \left( {{\eta_{cg}} + {\alpha_{fe}}{\eta_{fe}} + {\alpha_{de}}{\eta_{de}}} \right)\left( {{h_{r,k}} + {h_{c,k}}} \right)\]
Using \({\eta_{cg}} = 1  {\eta_{fe}}  {\eta_{de}}\) gives
\[{\overline h_{c,k}} = {h_{r,k}}\left( {{\eta_{fe}}\left( {{\alpha_{fe}}  1} \right) + {\eta_{de}}\left( {{\alpha_{de}}  1} \right)} \right) + {h_{c,k}}\left( {1 + {\eta_{fe}}\left( {{\alpha_{fe}}  1} \right) + {\eta_{de}}\left( {{\alpha_{de}}  1} \right)} \right)\]
This is the expression used by EnergyPlus for the gap convective conductance when a frame or divider is present.
Apportioning of Absorbed ShortWave Radiation in Shading Device Layers[LINK]
If a shading device has a nonzero shortwave transmittance then absorption takes place throughout the shading device layer. The following algorithm is used to apportion the absorbed shortwave radiation to the two faces of the layer. Here f_{1} is the fraction assigned to the face closest to the incident radiation and f_{2} is the fraction assigned to the face furthest from the incident radiation.
\({f_1} = 1,{\rm{ }}{f_2} = 0{\rm{ if }}{\tau_{sh}} = 0\)
Otherwise
\[\begin{array}{cl} f_1 = 0,\;f_2 = 0 &\text{if} \alpha_{sh} \le 0.01 \\ f_1 = 1,\;f_2 = 0 &\text{if} \alpha_{sh} \gt 0.999 \\ \left. \begin{array}{l} f_1 = \frac{1  e^{0.5 \ln (1  \alpha_{sh}}} {\alpha_{sh}} \\ f_2 = 1  {f_1} \end{array} \right\} &\text{if} 0.01 lt \alpha_{sh} \le 0.999 \end{array}\]
Window Frame and Divider Calculation[LINK]
For the zone heat balance calculation the inside surface temperature of the frame and that of the divider are needed. These temperatures are determined by solving the heat balance equations on the inside and outside surfaces of the frame and divider.
Table 38. Fortran Variables used in Window/Frame and Divider calculationsFrame Temperature Calculation[LINK]
Figure 97 shows a cross section through a window showing frame and divider. The outside and inside frame and divider surfaces are assumed to be isothermal. The frame and divider profiles are approximated as rectangular since this simplifies calculating heat gains and losses (see “Error Due to Assuming a Rectangular Profile,” below).
Figure 97. Cross section through a window showing frame and divider (exaggerated horizontally).
Frame Outside Surface Heat Balance[LINK]
The outside surface heat balance equation is
\[{Q_{ExtIR,abs}}  {Q_{IR,emitted}} + {Q_{conv}} + {Q_{cond}} + {Q_{abs}} = 0\]
where
\({Q_{ExtIR,abs}}\) = IR from the exterior surround (sky and ground) absorbed by outside frame surfaces
\({Q_{IR,emitted}}\) = IR emitted by outside frame surfaces
\({Q_{conv}}\) = convection from outside air to outside frame surfaces
\({Q_{cond}}\) = conduction through frame from inside frame surfaces to outside frame surfaces
\({Q_{abs}}\) = solar radiation (from sun, sky and ground) plus IR from outside window surface absorbed by outside frame surfaces (see “Calculation of Absorbed Solar Radiation,” below).
The first term can be written as the sum of the exterior IR absorbed by the outside face of the frame and the exterior IR absorbed by the frame’s outside projection surfaces.
\[{Q_{ExtIR,abs}} = {\varepsilon_1}{E_o}{A_f}{F_f} + {\varepsilon_1}{E_o}{A_{p1}}{F_{p1}}\]
where ε_{1} is the outside surface emissivity.
The exterior IR incident on the plane of the window, E_{o}, is the sum of the IR from the sky, ground and obstructions. For the purposes of the frame heat balance calculation it is assumed to be isotropic. For isotropic incident IR, F_{f} = 1.0 and F_{p1} = 0.5, which gives
\[{Q_{ExtIR,abs}} = {\varepsilon_1}{E_o}({A_f} + \frac{1}{2}{A_{p1}})\]
The IR emitted by the outside frame surfaces is
\[{Q_{ExtIR,emitted}} = {\varepsilon_1}\sigma ({A_f} + {A_{p1}})\theta_1^4\]
The convective heat flow from the outside air to the outside frame surfaces is
\[{Q_{conv}} = {h_{o,c}}({A_f} + {A_{p1}})({T_o}  {\theta_1})\]
The conduction through the frame from inside to outside is
\[{Q_{cond}} = k{A_f}({\theta_2}  {\theta_1})\]
Note that A_{f} is used here since the conductance, k, is, by definition, per unit area of frame projected onto the plane of the window.
Adding these expressions for the Q terms and dividing by* A_{f}* gives
\[{E_0}{\varepsilon_1}\left( {1 + \frac{1}{2}{\eta_1}} \right)  {\varepsilon_1}\left( {1 + {\eta_1}} \right)\theta_1^4 + {h_{o,c}}\left( {1 + {\eta_1}} \right)\left( {{T_0}  {\theta_1}} \right) + k\left( {{\theta_2}  {\theta_1}} \right) + {S_1} = 0\]
where S_{1} = Q_{abs}/A_{f} and
\[{\eta_1} = \frac{{{A_{p1}}}}{{{A_f}}} = \left( {\frac{{{p_{f,1}}}}{{{w_f}}}} \right)\frac{{H + W  ({N_h} + {N_v}){w_d}}}{{H + W + 2{w_f}}}\]
We linearize Eq. as follows.
Write the first two terms as
\[{\varepsilon_1}(1 + {\eta_1})[{E_o}(1 + \frac{1}{2}{\eta_1})/(1 + {\eta_1})  \theta_1^4]\]
and define a radiative temperature
\[{T_{o,r}} = {[{E_o}(1 + \frac{1}{2}{\eta_1})/(1 + {\eta_1})]^{1/4}}\]
This gives
\[{\varepsilon_1}(1 + {\eta_1})[T_{o,r}^4  \theta_1^4]\]
which, within a few percent, equals
\[{\varepsilon_1}(1 + {\eta_1})\frac{{{{({T_{o,r}} + {\theta_1})}^3}}}{2}({T_{o,r}}  {\theta_1})\]
Defining an outside surface radiative conductance as follows
\[{h_{o,r}} = {\varepsilon_1}(1 + {\eta_1})\frac{{{{({T_{o,r}} + {\theta_1})}^3}}}{2}\]
then gives
\[{h_{o,r}}({T_{o,r}}  {\theta_1})\]
The final outside surface heat balance equation in linearized form is then
\[{h_{o,r}}({T_{o,r}}  {\theta_1}) + {h_{o,c}}(1 + {\eta_1})({T_o}  {\theta_1}) + k({\theta_2}  {\theta_1}) + {S_1} = 0\]
Frame Inside Surface Heat Balance[LINK]
A similar approach can be used to obtain the following linearized inside surface heat balance equation:
\[{h_{i,r}}({T_{i,r}}  {\theta_2}) + {h_{i,c}}(1 + {\eta_2})({T_i}  {\theta_2}) + k({\theta_1}  {\theta_2}) + {S_2} = 0\]
where
\[{T_{i,r}} = {[{E_i}(1 + \frac{1}{2}{\eta_2})/(1 + {\eta_2})]^{1/4}}\]
\[{\eta_2} = \frac{{{A_{p2}}}}{{{A_f}}} = \left( {\frac{{{p_{f,2}}}}{{{w_f}}}} \right)\frac{{H + W  ({N_h} + {N_v}){w_d}}}{{H + W + 2{w_f}}}\]
and E_{i} is the interior IR irradiance incident on the plane of the window.
Solving Eqs. and simultaneously gives
\[{\theta_2} = \frac{{D + CA}}{{1  CB}}\]
with
\[A = \frac{{{h_{o,r}}{T_{o,r}} + {h_{o,c}}{T_o} + {S_1}}}{{{h_{o,r}} + k + {h_{o,c}}}}\]
\[B = \frac{k}{{{h_{o,r}} + k + {h_{o,c}}}}\]
\[C = \frac{k}{{{h_{i,r}} + k + {h_{i,c}}}}\]
\[D = \frac{{{h_{i,r}}{T_{i,r}} + {h_{i,c}}{T_i} + {S_2}}}{{{h_{i,r}} + k + {h_{i,c}}}}\]
Calculation of Solar Radiation Absorbed by Frame[LINK]
The frame outside face and outside projections and inside projections absorb beam solar radiation (if sunlight is striking the window) and diffuse solar radiation from the sky and ground. For the outside surfaces of the frame, the absorbed diffuse solar per unit frame face area is
\[Q_{abs,sol}^{dif} = I_{ext}^{dif}{\alpha_{fr,sol}}\left( {{A_f} + {F_{p1}}{A_{p1}}} \right)/{A_f} = I_{ext}^{dif}{\alpha_{fr,sol}}\left( {1 + 0.5\frac{{{A_{p1}}}}{{{A_f}}}} \right)\]
If there is no exterior window shade, I^{dif}_{ext} includes the effect of diffuse solar reflecting off of the glazing onto the outside frame projection, i.e.,
\[I_{ext}^{dif} \to I_{ext}^{dif}\left( {1 + R_{gl}^{f,dif}} \right)\]
The beam solar absorbed by the outside face of the frame, per unit frame face area is
\[Q_{abs,sol}^{bm,face} = I_{ext}^{bm}{\alpha_{fr,sol}}\cos {\beta_{face}}{f_{sunlit}}\]
The beam solar absorbed by the frame outside projection parallel to the window xaxis is
\[Q_{abs,sol}^{bm,h} = I_{ext}^{bm}{\alpha_{fr,sol}}\cos {\beta_h}{p_{f1}}\left( {W  {N_v}{w_d}} \right){f_{sunlit}}/{A_f}\]
Here it is assumed that the sunlit fraction, f_{sunlit}, for the window can be applied to the window frame. Note that at any given time beam solar can strike only one of the two projection surfaces that are parallel to the window xaxis. If there is no exterior window shade, I^{bm}_{ext} includes the effect of beam solar reflecting off of the glazing onto the outside frame projection, i.e.,
\[I_{ext}^{bm} \to I_{ext}^{bm}\left( {1 + R_{gl}^{f,bm}} \right)\]
The beam solar absorbed by the frame outside projection parallel to the window yaxis is
\[Q_{abs,sol}^{bm,v} = I_{ext}^{bm}{\alpha_{fr,sol}}\cos {\beta_v}{p_{f1}}\left( {H  {N_h}{w_d}} \right){f_{sunlit}}/{A_f}\]
Using a similar approach, the beam and diffuse solar absorbed by the inside frame projections is calculated, taking the transmittance of the glazing into account.
Error Due to Assuming a Rectangular Profile[LINK]
Assuming that the inside and outside frame profile is rectangular introduces an error in the surface heat transfer calculation if the profile is nonrectangular. The percent error in the calculation of convection and emitted IR is approximately 100\(\left {{L_{profile,rect}}  {L_{profile,actual}}} \right/{L_{profile,rect}}\) , where L_{profile,rect} is the profile length for a rectangular profile (w_{f} +* p_{f1}* for outside of frame or w_{f} + p_{f2}for inside of frame) and L_{profile,actual} is the actual profile length. For example, for a circular profile vs a square profile the error is about 22%. The error in the calculation of absorbed beam radiation is close to zero since the beam radiation intercepted by the profile is insensitive to the shape of the profile. The error in the absorbed diffuse radiation and absorbed IR depends on details of the shape of the profile. For example, for a circular profile vs. a square profile the error is about 15%.
Divider Temperature Calculation[LINK]
The divider inside and outside surface temperatures are determined by a heat balance calculation that is analogous to the frame heat balance calculation described above.
Beam Solar Reflection from Window Reveal Surfaces[LINK]
This section describes how beam solar radiation that is reflected from window reveal surfaces is calculated. Reflection from outside reveal surfaces—which are associated with the setback of the glazing from the outside surface of the window’s parent wall—increases the solar gain through the glazing. Reflection from inside reveal surfaces—which are associated with the setback of the glazing from the inside surface of the window’s parent wall—decreases the solar gain to the zone because some of this radiation is reflected back out of the window.
The amount of beam solar reflected from reveal surfaces depends, among other things, on the extent to which reveal surfaces are shadowed by other reveal surfaces. An example of this shadowing is shown in Figure 98. In this case the sun is positioned such that the top reveal surfaces shadow the left and bottom reveal surfaces. And the right reveal surfaces shadow the bottom reveal surfaces. The result is that the left/outside, bottom/outside, left/inside and bottom/inside reveal surfaces each have sunlit areas. Note that the top and right reveal surfaces are facing away from the sun in this example so their sunlit areas are zero.
Figure 98. Example of shadowing of reveal surfaces by other reveal surfaces.
The size of the shadowed areas, and the size of the corresponding illuminated areas, depends on the following factors:
The sun position relative to the window
The height and width of the window
The depth of the outside and inside reveal surfaces
We will assume that the reveal surfaces are perpendicular to the window plane and that the window is rectangular. Then the above factors determine a unique shadow pattern. From the geometry of the pattern the shadowed areas and corresponding illuminated areas can be determined. This calculation is done in subroutine CalcBeamSolarReflectedFromWinRevealSurface in the SolarShading module. The window reveal input data is specified in the WindowProperty:FrameAndDivider object expect for the depth of the outside reveal, which is determined from the vertex locations of the window and its parent wall.
If an exterior shading device (shade, screen or blind) is in place it is assumed that it blocks beam solar before it reaches outside or inside reveal surfaces. Correspondingly, it is assumed that an interior or betweenglass shading device blocks beam solar before it reaches inside reveal surfaces.
Representative shadow patterns are shown in Figure 99 for a window with no shading device, and without and with a frame. The case with a frame has to be considered separately because the frame can cast an additional shadow on the inside reveal surfaces.
The patterns shown apply to both vertical and horizontal reveal surfaces. It is important to keep in mind that, for a window of arbitrary tilt, if the left reveal surfaces are illuminated the right surfaces will not be, and vice versa. And if the bottom reveal surfaces are illuminated the top surfaces will not be, and vice versa. (Of course, for a vertical window, the top reveal surfaces will never be illuminated by beam solar if the reveal surfaces are perpendicular to the glazing, as is being assumed.
For each shadow pattern in Figure 99, equations are given for the shadowed areas \({A_{1,sh}}\) and \({A_{2,sh}}\) of the outside and inside reveal surfaces, respectively. The variables in these equations are the following (see also Figure 100):
\({d_1}\) = depth of outside reveal, measured from the outside plane of the glazing to the edge of the reveal, plus one half of the glazing thickness.
\({d_2}\) = depth of inside reveal (or, for illumination on bottom reveal surfaces, inside sill depth), measured from the inside plane of the glazing to the edge of the reveal or the sill, plus one half of the glazing thickness.
\(L\) = window height for vertical reveal surfaces or window width for horizontal reveal surfaces
\(\alpha\) = vertical solar profile angle for shadowing on vertical reveal surfaces or horizontal solar profile angle for shadowing on horizontal reveal surfaces.
\({p_1}({p_2})\) = distance from outside (inside) surface of frame to glazing midplane.
\({d_2}^\prime\) = depth of shadow cast by top reveal on bottom reveal, or by left reveal on right reveal, or by right reveal on left reveal.
\({d_2}^{\prime \prime }\) = depth of shadow cast by frame.
For simplicity it is assumed that, for the case without a frame, the shadowed and illuminated areas extend into the glazing region. For this reason, \({d_1}\) and \({d_2}\) are measured from the midplane of the glazing. For the case with a frame, the beam solar absorbed by the surfaces formed by the frame outside and inside projections perpendicular to the glazing is calculated as described in “Window Frame and Divider Calculation: Calculation of Solar Radiation Absorbed by Frame.”
Figure 99. Expression for area of shaded regions for different shadow patterns: (a) window without frame, (b) window with frame
Figure 100. Vertical section through a vertical window with outside and inside reveal showing calculation of the shadows cast by the top reveal onto the inside sill and by the frame onto the inside sill.
The following logic gives expressions for the shadowed areas for all possible shadow patterns. Here:
d1 = \({d_1}\)
d2 = \({d_2}\)
P1 = \({p_1}\)
P2 = \({p_2}\)
f1 = \({d_1}  {p_1}\)
f2 = \({d_2}  {p_2}\)
d2prime = \({d_2}^\prime\)
d2prime2 = \({d_2}^{\prime \prime }\)
d12 = \({d_1} + {d_2}  {d_2}^\prime\)
TanAlpha = \(\tan \alpha\)
A1sh = \({A_{1,sh}}\)
A2sh = \({A_{2,sh}}\)
L = \(L\)
L1 = average distance to frame of illuminated area of outside reveal (used to calculate view factor to frame).
L2 = average distance to frame of illuminated area of inside reveal (used to calculate view factor to frame).
The beam solar reflected from a sunlit region of area \(A\) is given by
\[R = {I_B}A\cos \beta (1  a){\rm{ }}\]
where
\(R\) = reflected solar radiation [W]
\({I_B}\) = beam normal irradiance [W/m^{2}]
\(A\) = sunlit area [m^{2}]
\(\beta\) = beam solar angle of incidence on reveal surface
\(a\) = solar absorptance of reveal surface
All reflected radiation is assumed to be isotropic diffuse. For outside reveal surfaces it is assumed that \(R/2\) goes toward the window and \(R/2\) goes to the exterior environment. Of the portion that goes toward the window a fraction \({F_1}\) goes toward the frame, if present, and \(1  {F_1}\) goes toward the glazing.
The view factor \({F_1}\) to the frame calculated by assuming that the illuminated area can be considered to be a line source. Then the areaweighted average distance, \({L_1}\) , of the source to the frame is calculated from the shape of the illuminated area (see above psuedocode). Then \({F_1}\) is related as follows to the average angle subtended by the frame of width \({w_f}\) :
\[{F_1} = \frac{{{{\tan }^{  1}}({w_f}/{L_1})}}{{\pi /2}}\]
For the portion going towards the frame, \((R/2){F_1}{a_f}\) is absorbed by the frame (where \({a_f}\) is the solar absorptance of the frame) and contributes to the frame heat conduction calculation. The rest, \((R/2){F_1}(1  {a_f})\) , is assumed to be reflected to the exterior environment.
If the glazing has diffuse transmittance \({\tau_{diff}}\) , diffuse front reflectance \(\rho_{diff}^f\) , and layer front absorptance \(\alpha_{l,diff}^f\) , then, of the portion, \((R/2)(1  {F_1})\) , that goes toward the glazing, \((R/2)(1  {F_1}){\tau_{diff}}\) is transmitted to the zone, \((R/2)(1  {F_1})\alpha_{l,diff}^f\) is absorbed in glass layer \(l\) and contributes to the glazing heat balance calculation, and \((R/2)(1  {F_1})\rho_{dif{f_{}}}^f\) is reflected to the exterior environment.
The beam solar absorbed by an outside reveal surface is added to the other solar radiation absorbed by the outside of the window’s parent wall.
For inside reveal surfaces it is assumed that \(R/2\) goes towards the window and \(R/2\) goes into the zone. Of the portion that goes toward the window a fraction \((R/2){F_2}\) goes toward the frame, if present, and \((R/2)(1  {F_2})\) goes toward the glazing (\({F_2}\) is calculated using a method analogous to that used for \({F_1}\) ). For the portion going towards the frame, \((R/2){F_2}{a_f}\) is absorbed by the frame and contributes to the frame heat conduction calculation. The rest, \((R/2){F_2}(1  {a_f})\) , is assumed to be reflected back into the zone.
If the glazing has diffuse back reflectance \(\rho_{diff}^b\) , and layer back absorptance \(\alpha_{l,diff}^b\) , then, of the portion \((R/2)(1  {F_2})\) that goes toward the glazing, \((R/2)(1  {F_2}){\tau_{diff}}\) is transmitted back out the glazing, \((R/2)(1  {F_2})\alpha_{l,diff}^b\) is absorbed in glass layer \(l\) and contributes to the glazing heat balance calculation, and \((R/2)(1  {F_2})\rho_{dif{f_{}}}^b\) is reflected into the zone.
The beam solar absorbed by an inside reveal surface is added to the other solar radiation absorbed by the inside of the window’s parent wall.
Shading Device Thermal Model[LINK]
Shading devices in EnergyPlus can be on the exterior or interior sides of the window or between glass layers. The window shading device thermal model accounts for the thermal interactions between the shading layer (shade, screen or blind) and the adjacent glass, and between the shading layer and the room (for interior shading) or the shading layer and the outside surround (for exterior shading).
An important feature of the shading device thermal model is calculating the natural convection airflow between the shading device and glass. This flow affects the temperature of the shading device and glazing and, for interior shading, is a determinant of the convective heat gain from the shading layer and glazing to the zone air. The airflow model is based on one described in the ISO Standard 15099, “Thermal Performance of Windows, Doors and Shading Devices—Detailed Calculations” [ISO15099, 2001]. (Betweenglass forced airflow is also modeled; see “Airflow Windows.”)
The following effects are considered by the shading device thermal model:
For interior and exterior shading device: Longwave radiation (IR) from the surround absorbed by shading device, or transmitted by the shading device and absorbed by the adjacent glass. For interior shading the surround consists of the other zone surfaces. For exterior shading the surround is the sky and ground plus exterior shadowing surfaces and exterior building surfaces “seen” by the window.
Interreflection of IR between the shading device and adjacent glass.
Direct and diffuse solar radiation absorbed by the shading device.
Interreflection of solar radiation between shading layer and glass layers.
Convection from shading layer and glass to the air in the gap (or, for betweenglass shading, gaps) between the shading layer and adjacent glass, and convection from interior shading layer to zone air or from exterior shading layer to outside air.
Natural convection airflow in the gap (or, for betweenglass shading, gaps) between shading layer and adjacent glass induced by buoyancy effects, and the effect of this flow on the shadingtogap and glasstogap convection coefficients.
For interior shading, convective gain (or loss) to zone air from gap airflow.
In the following it is assumed that the shading device, when in place, covers the glazed part of the window (and dividers, if present) and is parallel to the glazing. For interior and exterior shading devices it is assumed that the shading layer is separated from the glazing by an air gap. A betweenglass shading layer is assumed to be centered between two glass layers and separated from the adjacent glass layers by gaps that is filled with the same gas. If the window has a frame, it is assumed that the shading device does not cover the frame.
Heat Balance Equations for Shading Device and Adjacent Glass[LINK]
If a window shading device is deployed the heat balance equations for the glass surfaces facing the shading layer are modified, and two new equations, one for each face of the shading layer, are added. Figure 101 illustrates the case of double glazing with an interior shading device.
Figure 101. Glazing system with two glass layers and an interior shading layer showing variables used in heat balance equations.
The heat balance equation for the glass surface facing the gap between glass and shading layer (called in the following, “gap”) is
\[\frac{{{E_i}{\varepsilon_4}{\tau_{sh}}}}{{1  {\rho_4}{\rho_{sh}}}} + \frac{{\sigma {\varepsilon_4}}}{{1  {\rho_4}{\rho_{sh}}}}\left[ {\theta_5^4{\varepsilon_{sh}}  \theta_4^4(1  {\rho_{sh}})} \right] + {k_2}\left( {{\theta_3}  {\theta_4}} \right) + {h_{cv}}\left( {{T_{gap}}  {\theta_4}} \right) + {S_4} = 0\]
where
τ_{sh} = IR diffuse transmittance of shading device
ε_{sh} = diffuse emissivity of shading device
ρ_{sh} = IR diffuse reflectance of shading device ( = 1  ( τ_{sh} + ε_{sh}))
θ_{5} = temperature of the surface of the shading layer that faces the gap (K).
The term 1 – ρ_{4} ρ_{sh} accounts for the interreflection of IR radiation between glass and shading layer.
The convective heat transfer from glass layer #2 to the air in the gap is
\[{q_{c,gl}} = {h_{cv}}({\theta_4}  {T_{gap}})\]
where
T_{gap} = effective mean temperature of the gap air (K).
h_{cv} = convective heat transfer coefficient from glass or shading layer to gap air (W/m^{2}K).
The corresponding heat transfer from shading layer to gap air is
\[{q_{c,sh}} = {h_{cv}}({\theta_5}  {T_{gap}})\]
The convective heat transfer coefficient is given by
\[{h_{cv}} = 2{h_c} + 4v\]
where
h_{c} = surfacetosurface heat transfer coefficient for nonvented (closed) cavities (W/m^{2}K)
v = mean air velocity in the gap (m/s).
The quantities h_{cv} and T_{gap} depend on the airflow velocity in the gap, which in turn depends on several factors, including height of shading layer, glass/shading layer separation (gap depth), zone air temperature for interior shading or outside air temperature for exterior shading, and shading layer and glass face temperatures. The calculation of h_{cv} and T_{gap} is described in the following sections.
The heat balance equation for the shading layer surface facing the gap is
\[\frac{E_i \tau_{sh} \rho_4 \varepsilon_{sh}} {1  \rho_4 \rho_{sh}} + \frac{\sigma \varepsilon_{sh}} {1  \rho_4 \rho_{sh}} \left[ \varepsilon_4 \theta_4^4  \theta_5^4 ( 1  \rho_4( \varepsilon_{sh} + \rho_{sh} ) ) \right] + k_{sh} \left( \theta_6  \theta_5 \right) + h_{cv} \left( T_{gap}  \theta_5 \right) + S_{sh,1} = 0\]
where
k_{sh} = shading layer conductance (W/m^{2}K).
θ_{6} = temperature of shading layer surface facing the zone air (K).
S_{sh,1} = solar radiation plus shortwave radiation from lights plus IR radiation from lights and zone equipment absorbed by the gapside face of the shading layer (W/m^{2}K).
The heat balance equation for the shading layer surface facing the zone air is
\[{E_i}{\varepsilon_{sh}}  {\varepsilon_{sh}}\sigma \theta_6^4 + {k_{sh}}\left( {{\theta_5}  {\theta_6}} \right) + {h_i}\left( {{T_i}  {\theta_6}} \right) + {S_{sh,2}} = 0\]
where
S_{sh,2} = solar radiation plus shortwave radiation from lights plus IR radiation from lights and zone equipment absorbed by the zoneside face of the shading layer (W/m^{2}K).
Solving for Gap Airflow and Temperature[LINK]
For interior and exterior shading devices a pressurebalance equation is used to determine gap air velocity, gap air mean equivalent temperature and gap outlet air temperature given values of zone air temperature (or outside temperature for exterior shading), shading layer face temperatures and gap geometry. The pressure balance equates the buoyancy pressure acting on the gap air to the pressure losses associated with gap airflow between gap inlet and outlet [ISO15099, 2001]. The variables used in the following analysis of the interior shading case are shown in Figure 102.
Figure 102. Vertical section (a) and perspective view (b) of glass layer and interior shading layer showing variables used in the gap airflow analysis. The opening areas A_{bot}, A_{top}, A_{l}, A_{r} and A_{h} are shown schematically.
Pressure Balance Equation[LINK]
The pressure balance equation for airflow through the gap is
\[\Delta {p_T} = \Delta {p_B} + \Delta {p_{HP}} + \Delta {p_Z}\]
Here, Δp_{T}is the driving pressure difference between room air and gap air. It is given by
\[\Delta {p_T} = {\rho_0}{T_0}gH\sin \phi \frac{{\left {{T_{gap}}  {T_{gap,in}}} \right}}{{{T_{gap}}{T_{gap,in}}}}\]
where
ρ_{0} = density of air at temperature T_{0} (kg/m^{3})
T_{0} = reference temperature (283K)
g = acceleration due to gravity (m/s^{2})
H = height of shading layer (m)
φ = tilt angle of window (vertical = 90^{o})
T_{gap} = effective mean temperature of the gap air (K)
T_{gap,in} = gap inlet temperature ( = zone air temperature for interior shading) (K)
The Δp_{B} term is due to the acceleration of air to velocity v (Bernoulli’s law). It is given by
\[\Delta {p_B} = \frac{\rho }{2}{v^2}{\rm{ (Pa)}}\]
where ρis the gap air density evaluated at T_{gap} (kg/m^{3}).
The Δp_{HP} term represents the pressure drop due to friction with the shading layer and glass surfaces as the air moves through the gap. Assuming steady laminar flow, it is given by the HagenPoiseuille law for flow between parallel plates [Munson et al. 1998]:
\[\Delta {p_{HP}} = 12\mu \frac{H}{{{s^2}}}{v^{}}{\rm{ (Pa)}}\]
where μ is the viscosity of air at temperature T_{gap} (Pas).
The Δp_{Z} term is the sum of the pressure drops at the inlet and outlet openings:
\[\Delta {p_Z} = \frac{{\rho {v^2}}}{2}\left( {{Z_{in}} + {Z_{out}}} \right){\rm{ (Pa)}}\]
Here, the inlet pressure drop factor,* Z_{in}, and the outlet pressure drop factor,Z_{out}*, are given by
\[\begin{array}{rl} Z_{in} & = \left( \frac{A_{gap}}{0.66A_{eq,in}}  1 \right)^2 \\ Z_{out} & = \left( \frac{A_{gap}}{0.60A_{eq,out}}  1 \right)^2 \end{array}\]
\[\begin{array}{l}{Z_{in}} = {\left( {\frac{{{A_{gap}}}}{{0.66{A_{eq,in}}}}  1} \right)^2}\\\{Z_{out}} = {\left( {\frac{{{A_{gap}}}}{{0.60{A_{eq,out}}}}  1} \right)^2}\end{array}\]
where
A_{eq,in} = equivalent inlet opening area (m^{2})
A_{eq,out} = equivalent outlet opening area (m^{2})
A_{gap} = crosssectional area of the gap = sW(m^{2})
If T_{gap} > T_{gap,in}
\[\begin{array}{rl} A_{eq,in} & = A_{bot} + \frac{A_{top}}{2(A_{bot}+A_{top}} \left( A_l + A_r + A_h \right) \\ A_{eq,out} & = A_{top} + \frac{A_{bot}}{2(A_{bot}+A_{top}} \left( A_l + A_r + A_h \right) \end{array}\]
If T_{gap} ≤ T_{gap,in}
\[\begin{array}{rl} A_{eq,in} & = A_{top} + \frac{A_{bot}}{2(A_{bot}+A_{top}} \left( A_l + A_r + A_h \right) A_{eq,out} & = A_{bot} + \frac{A_{top}}{2(A_{bot}+A_{top}} \left( A_l + A_r + A_h \right) \\ \end{array}\]
Here, the area of the openings through which airflow occurs (see Figure 102 and Figure 103) are defined as follows:
A_{bot} = area of the bottom opening (m^{2})
A_{top} = area of the top opening (m^{2})
A_{l} = area of the leftside opening (m^{2})
A_{r} = area of the rightside opening (m^{2})
A_{h} = air permeability of the shading device expressed as the total area of openings (“holes”) in the shade surface (these openings are assumed to be uniformly distributed over the shade) (m^{2})
Figure 103 shows examples of A_{bot}, A_{top}, A_{l} and A_{r} for different shading device configurations. These areas range from zero to a maximum value equal to the associated shade/screen/blindtoglass crosssectional area; i.e., A_{bot} and A_{top} ≤ sW, A_{l} and A_{r} ≤ sH.
Figure 103. Examples of openings for an interior shading layer covering glass of height H and width W. Not to scale. (a) Horizontal section through shading layer with openings on the left and right sides (top view). (b) Vertical section through shading layer with openings at the top and bottom (side view).
Expression for the Gap Air Velocity
Expressing Equation in terms of v yields the following quadratic equation:
\[\frac{{\rho {v^2}}}{2}\left( {1 + {Z_{in}} + {Z_{out}}} \right) + \frac{{12\mu H}}{{{s^2}}}v  {\rho_0}{T_0}gH\sin \phi \frac{{\left {{T_{gap,in}}  {T_{gap}}} \right}}{{{T_{gap,in}}{T_{gap}}}} = 0\]
Solving this gives
\[v = \frac{{{{\left[ {{{\left( {\frac{{12\mu H}}{{{s^2}}}} \right)}^2} + \frac{{2{\rho ^2}(1 + {Z_{in}} + {Z_{out}}){\rho_0}{T_0}gH\sin \phi \left {{T_{gap,in}}  {T_{gap}}} \right}}{{{T_{gap,in}}{T_{gap}}}}} \right]}^{1/2}}  \frac{{12\mu H}}{{{s^2}}}}}{{\rho (1 + {Z_{in}} + {Z_{out}})}}\]
The choice of the root of the quadratic equation is dictated by the requirement that v = 0 if T_{gap,in} = T_{gap}.
Gap Outlet Temperature and Equivalent Mean Air Temperature
The temperature of air in the gap as a function of distance, h, from the gap inlet (Figure 104) is
\[{T_{gap}}(h) = {T_{ave}}  ({T_{ave}}  {T_{gap,in}}){e^{  h/{H_0}}}\]
where
\[{T_{ave}} = \frac{{{T_{gl}} + {T_{sh}}}}{2}\]
is the average temperature of the glass and shading layer surfaces facing the gap (K).
H_{0} = characteristic height (m), given by
\[{H_0} = \frac{{\rho {C_p}s}}{{2{h_{cv}}}}v\]
where C_{p} is the heat capacity of air.
The gap outlet temperature is given by
\[{T_{gap,out}} = {T_{ave}}  ({T_{ave}}  {T_{gap,in}}){e^{  H/{H_0}}}\]
The thermal equivalent mean temperature of the gap air is
\[{T_{gap}} = \frac{1}{H}\int\limits_0^H {{T_{gap}}(h)dh = {T_{ave}}  \frac{{{H_0}}}{H}} \left( {{T_{gap,out}}  {T_{gap,in}}} \right)\]
Figure 104. Variation of gap air temperature with distance from the inlet for upward flow.
Solution Sequence for Gap Air Velocity and Outlet Temperature
The routine WinShadeGapFlow is called within the glazing heat balance iterative loop in SolveForWindowTemperatures to determine v and T_{gap,out}. The solution sequence in WinShadeGapFlow is as follows:
At start of iteration, guess T_{gap} as ((T_{gl} + T_{sh})/2 + T_{gap,in})/2. Thereafter use value from previous iteration.
Get stillair conductance, h_{c}, by calling WindowGasConductance and NusseltNumber.
Get v from Equation
Get h_{cv} from Equation
Get T_{ave}from Equation
Get T_{gap,out} from Equation
Get new value of T_{gap}from Equation
The values of h_{cv} and T_{gap} so determined are then used in the window heat balance equations to find new values of the face temperatures of the glass and shading layers. These temperatures are used in turn to get new values of h_{cv} and T_{gap} until the whole iterative process converges.
Convective Heat Gain to Zone from Gap Airflow
The heat added (or removed) from the air as it passes through the gap produces a convective gain (or loss) to the zone air given by
\[{q_v} = LW\left[ {{h_{cv}}\left( {{T_{gl}}  {T_{gap}}} \right) + {h_{cv}}\left( {{T_{sh}}  {T_{gap}}} \right)} \right] = 2{h_{cv}}LW\left( {{T_{ave}}  {T_{gap}}} \right){\rm{ (W)}}\]
This can also be expressed as
\[{q_v} = \dot m{C_p}\left( {{T_{gap,out}}  {T_{gap,in}}} \right){\rm{ (W)}}\]
where the air mass flow rate in the gap is given by
\[\dot m = \rho {A_{gap}}v{\rm{ (kg/s)}}\]
Heat Balance Equations for BetweenGlass Shading Device[LINK]
In EnergyPlus shading devices are allowed between the two glass panes of double glazing and between the two inner glass panes of triple glazing. Figure 105 shows the case of a betweenglass shading device in double glazing.
Figure 105. Glazing system with two glass layers and a betweenglass shading device showing variables used in the heat balance equations.
The heat balance equations for the two glass surfaces facing the shading device are the following.
For face #2:
\[h_{cv,1} (T_{gap,1}\theta_2) + k_1 (\theta_1\theta_2) + \frac{\sigma \epsilon_2}{1\rho_2 R_1} \left[ \frac{\tau_{sh}}{1\rho_6\rho_3} ( \epsilon_3 \theta_3^4 + \epsilon_6\theta_6^4\rho_3 ) + \epsilon_5\theta_5^4 + \epsilon_2\theta_2^4R_1 \right]\]
where
\[{R_1} = {\rho_5} + \frac{{\tau_{sh}^2{\rho_3}}}{{1  {\rho_6}{\rho_3}}}\]
${T_{gap,1}} = $ effective mean air temperature in gap 1 (K)
${h_{cv,1}} = $ convective heat transfer coefficient from glass or shading layer to gas in gap 1 (W/m^{2}K)
For face #3:
\[h_{cv,2} (T_{gap,2}\theta_3) + k_2 (\theta_4\theta_3) + \frac{\sigma \epsilon_3}{1\rho_3 R_2} \left[ \frac{\tau_{sh}}{1\rho_5\rho_2} ( \epsilon_3 \theta_2^4 + \epsilon_5\theta_5^4\rho_2 ) + \epsilon_6\theta_6^4 + \epsilon_7\theta_3^4R_2 \right]\]
where
\[{R_2} = {\rho_6} + \frac{{\tau_{sh}^2{\rho_2}}}{{1  {\rho_5}{\rho_2}}}\]
${T_{gap,2}} = $ effective mean air temperature in gap 2 (K)
${h_{cv,2}} = $ convective heat transfer coefficient from glass or shading layer to gas in gap 2 (W/m^{2}K)
The heat balance equations for the shading layer faces are:
For face #5:
\[h_{cv,1} (T_{gap,1}\theta_5) + k_sh (\theta_6\theta_5) + \frac{\sigma \epsilon_5}{1\rho_2 R_1} \left[ \frac{\tau_{sh}\rho_2}{1\rho_5\rho_3} ( \epsilon_3 \theta_3^4 + \epsilon_6\theta_6^4\rho_3 ) + \epsilon_2\theta_2^4 + \epsilon_5\theta_5^4\rho_2 \right]\]
For face #6:
\[h_{cv,2} (T_{gap,2}\theta_3) + k_2 (\theta_4\theta_3) + \frac{\sigma \epsilon_3}{1\rho_3 R_2} \left[ \frac{\tau_{sh}}{1\rho_5\rho_2} ( \epsilon_2 \theta_2^4 + \epsilon_5\theta_5^4\rho_2 ) + \epsilon_2\theta_2^4 + \epsilon_7\theta_3^4R_2 \right]\]
The convective heat transfer coefficients are given by
\[\begin{array}{l}{h_{cv,1}} = 2{h_{c,1}} + 4v\\\{h_{cv,2}} = 2{h_{c,2}} + 4v\end{array}\]
where
${h_{c,1}},{h_{c,2}} = $ surfacetosurface heat transfer coefficients for gap #1 and #2, respectively, when these gaps are nonvented (closed).
$v = $ air velocity in the gaps (m/s). It is assumed that the gap widths are equal, so that the velocity of flow in the gaps is equal and opposite, i.e., when the airflow is upward in gap #1 it is downward in gap #2 and viceversa.
In analogy to the interior or exterior shading device case, the air velocity is determined by solving the following pressure balance equation:
\[\Delta {p_{T,1,2}} = \Delta {p_{B,1}} + \Delta {p_{HP,1}} + \Delta {p_{Z,1}} + \Delta {p_{B,2}} + \Delta {p_{HP,2}} + \Delta {p_{Z,2}}\]
where the driving pressure difference between gap #1 and #2 is
\[\Delta {p_{T,1,2}} = {\rho_0}{T_0}gH\sin \phi \frac{{\left {{T_{gap,1}}  {T_{gap,2}}} \right}}{{{T_{gap,1}}{T_{gap,2}}}}{\rm{ (Pa)}}\]
The pressure drops on the righthand side of this equation are:
\[\begin{array}{l}\Delta {p_{B,i}} = \frac{{{\rho_{gap,}}_i}}{2}{v^2}\\\Delta {p_{HP,i}} = 12{\mu_{gap,}}_i\frac{H}{{{s^2}}}\\\Delta {p_{Z,i}} = \frac{{{\rho_{gap,}}_i{v^2}}}{2}({Z_{in,i}} + {Z_{out,i}})\end{array}\]
where i = gap number (1 or 2).
It can be shown that \({Z_{in,1}} + {Z_{out,1}} = {Z_{in,2}} + {Z_{out,2}}\) . Then, inserting these pressure drop expressions into , we obtain the following expression for the airflow velocity:
\[v = \frac{{{{\left[ {{{\left( {\frac{{12({\mu_{gap,}}_1 + {\mu_{gap,}}_2)H}}{{{s^2}}}} \right)}^2} + 2\Delta {p_{T,1,2}}({\rho_{gap,1}} + {\rho_{gap,2}})(1 + {Z_{in}} + {Z_{out}})} \right]}^{1/2}}  \frac{{12({\mu_{gap,}}_1 + {\mu_{gap,}}_2)H}}{{{s^2}}}}}{{({\rho_{gap,}}_1 + {\rho_{gap,}}_2)(1 + {Z_{in}} + {Z_{out}})}}\]
The choice of the sign of the square root term is dictated by the requirement that \(v = 0\) if \(\Delta {p_{T,1,2}} = 0\) , i.e., \({T_{gap,1}} = {T_{gap,2}}\) .
Given \(v\) we can now calculate \({T_{gap,1}}\) and \({T_{gap,2}}\) , which gives \(\Delta {p_{T,1,2}}\) . The procedure is as follows. We have
\[{T_{gap,1,out}} = {T_{ave,1}}  ({T_{ave,1}}  {T_{gap,1,in}}){\xi_1}\]
where \({T_{ave,1}} = ({\theta_2} + {\theta_5})/2\) and \({\xi_1} = {e^{  \frac{H}{{{H_{0,1}}}}}}\) with \({H_{0,1}} = {\rho_{gap,1}}{C_p}sv/(2{h_{cv,1}})\) . Since \({T_{gap,1,in}} = {T_{gap,2,out}}\) this gives:
\({T_{gap,1,out}} = {T_{ave,1}}  ({T_{ave,1}}  {T_{gap,2,out}}){\xi_1}\)
Similarly,
\[{T_{gap,2,out}} = {T_{ave,2}}  ({T_{ave,2}}  {T_{gap,1,out}}){\xi_2}\]
Solving these simultaneous equations gives:
\(\begin{array}{l}{T_{gap,1,out}} = \frac{{{T_{ave,1}}(1  {\xi_1}) + {\xi_1}{T_{ave,2}}(1  {\xi_2})}}{{1  {\xi_1}{\xi_2}}}\\\\\{T_{gap,2,out}} = \frac{{{T_{ave,2}}(1  {\xi_2}) + {\xi_2}{T_{ave,1}}(1  {\xi_1})}}{{1  {\xi_1}{\xi_2}}}\end{array}\)
Using these in
\[{T_{gap,1}} = {T_{ave,1}}  \frac{{{H_{0,1}}}}{H}\left( {{T_{gap,1,out}}  {T_{gap,2,out}}} \right)\]
gives
\({T_{gap,1}} = {T_{ave,1}}  \frac{{{H_{0,1}}}}{H}\xi \left( {{T_{ave,1}}  {T_{ave,2}}} \right)\)
with
\[\xi = \frac{{(1  {\xi_1})(1  {\xi_2})}}{{1  {\xi_1}{\xi_2}}}\]
Similarly, from
\[{T_{gap,2}} = {T_{ave,2}}  \frac{{{H_{0,2}}}}{H}\left( {{T_{gap,2,out}}  {T_{gap,1,out}}} \right)\]
we get
\[{T_{gap,2}} = {T_{ave,2}}  \frac{{{H_{0,2}}}}{H}\xi \left( {{T_{ave,2}}  {T_{ave,1}}} \right)\]
The overall solution sequence is as follows.
At start of iteration guess \({T_{gap,1}} = {T_{ave,1}}\) and \({T_{gap,2}} = {T_{ave,2}}\) . Then
1) Get \({\rho_{gap,1}},{\rho_{gap,2}},{\mu_{gap,1}},{\mu_{gap,2}}\) using \({T_{gap,1}},{T_{gap,2}}\) .
2) Get stillair conductances \({h_{c,1}},{h_{c,2}}\) by calling WindowGasConductance and NusseltNumber.
3) Get \(v\) from Equation
4) Get \({h_{cv,1}},{h_{cv,2}}\) from Equation
5) Get \({T_{ave,1}},{T_{ave,2}}\)
6) Get H_{o,1}, H_{o,2}, \({\xi_1}\) and \({\xi_2}\) .
7) Get \({T_{gap,1}}\) ,\({T_{gap,2}}\) from Equations and
The values \({h_{cv,1}}{\rm{, }}{h_{cv,2}}{\rm{, }}{T_{gap,1}}\) and \({T_{gap,2}}\) are then used in the face heat balance equations to find new values of the face temperatures \({\theta_2}{\rm{, }}{\theta_3},{\rm{ }}{\theta_5}\) and \({\theta_6}\) . These are used in turn to get new values of \({h_{cv,1}}{\rm{, }}{h_{cv,2}}{\rm{, }}{T_{gap,1}}\) and \({T_{gap,2}}\) until the whole iterative process converges.
Airflow Windows[LINK]
In airflow windows forced air flows in the gap between adjacent layers of glass. Such windows are also known as “heatextract windows” and “climate windows.”
Five configurations of airflow windows are modeled (Figure 106) that depend on the source and destination of forced air. The allowed combinations of Airflow Source and Airflow Destination are:
InsideAir à OutsideAir
InsideAir à InsideAir
InsideAir à ReturnAir
OutsideAir à InsideAir
OutsideAir à OutsideAir
Figure 106. Gap airflow configurations for airflow windows. From “Active facades,” Version no. 1, Belgian Building Research Institute, June 2002.
A common application of airflow windows is to reduce the zone cooling load by exhausting indoor air through the window, thereby picking up and rejecting heat from the glazing (Figure 106).
Figure 107 shows the variables used in the heat balance equations for forced airflow in a doubleglazed window.
Figure 107. Glazing system with forced airflow between two glass layers showing variables used in the heat balance equations.
The heat balance equation for the lefthand glass surface facing the gap in Figure 107 is:
\[{k_1}({\theta_1}  {\theta_2}) + {h_{cv}}({T_{gap}}  {\theta_2}) + \sigma \frac{{{\varepsilon_2}{\varepsilon_3}}}{{1  (1  {\varepsilon_2})(1  {\varepsilon_3})}}\left( {\theta_3^4  \theta_2^4} \right) + {S_2} = 0\]
The corresponding equation for the righthand glass surface facing the gap is:
\[{k_2}({\theta_4}  {\theta_3}) + {h_{cv}}({T_{gap}}  {\theta_3}) + \sigma \frac{{{\varepsilon_2}{\varepsilon_3}}}{{1  (1  {\varepsilon_2})(1  {\varepsilon_3})}}\left( {\theta_2^4  \theta_3^4} \right) + {S_3} = 0\]
Here,
T_{gap} = effective mean temperature of the gap air (K)
h_{cv} = convective heat transfer coefficient from glass to gap air (W/m2K).
The convective heat transfer coefficient is given by
\[{h_{cv}} = 2{h_c} + 4v\]
where
h_{c} = glasstoglass heat transfer coefficient for nonvented (closed) cavity (W/m^{2}K)
v = mean air velocity in the gap (m/s).
The air velocity is determined by the gap crosssectional area in the flow direction and the air flow rate, which is an input value that is constant or can vary according to a userspecified schedule:
\[v = \frac{F}{{{A_{gap}}}}{\rm{ (m/s)}}\]
where
F = airflow rate (m^{3}/s)
A_{gap} = gap crosssectional area in direction of flow (m2)
It is assumed that the airflow is uniform across the width of the window.
The mean temperature of the gap air is given by the following expression, whose derivation follows that for for the case of an interior shading device:
\[{T_{gap}} = {T_{ave}}  \frac{H}{{{H_0}}}({T_{gap,out}}  {T_{gap,in}})\]
where
\[{T_{ave}} = \frac{{{\theta_2} + {\theta_3}}}{2}\]
H = glazing height (m)
\[{H_0} = \frac{{\rho {C_p}s}}{{2{h_{cv}}}}v\]
T_{gap,in}_{ } = gap air inlet temperature (T_{i} if the airflow source is indoor air, T_{o} if the airflow source is outside air) (K)
The outlet air temperature is given by
\[{T_{gap,out}} = {T_{ave}}  ({T_{ave}}  {T_{gap,in}}){e^{  H/{H_0}}}\]
The equations for glass face #1 and #4 are the same as those for no airflow in the gap (Equations and ).
The convective heat gain to the zone air due to the gap airflow when the airflow destination is indoor air is
\[{q_v} = \dot m({C_{p,out}}{T_{gap,out}}  {C_{p,i}}{T_i}){\rm{ (W)}}\]
where
C_{p,i}_{ } = heat capacity of the indoor air (J/kgK)
C_{p,out} = heat capacity of the gap outlet air (J/kgK)
and where the air mass flow rate in the gap is
\[\dot m = \rho F{\rm{ (kg/s)}}\]
Fan Energy[LINK]
The fan energy used to move air through the gap is very small and is ignored.
Airflow Window with BetweenGlass Shading Device[LINK]
Figure 108 shows the case of a doubleglazed airflow window with a between glass shading device. The heat balance equations in this case are the same as those for the betweenglass shading device with natural convection (Figure 105 and following equations) except that now
\[v = \frac{{F/2}}{{{A_{gap}}}}{\rm{ (m/s)}}\]
where A_{gap} = sW is the crosssectional area of the gap on either side of the shading device. It is assumed that the shading device is centered between the two panes of glass so that the airflow, F, is divided equally between the two gaps.
The convective heat gain to the zone air due to the airflow through the two gaps when the airflow destination is indoor air is
\[{q_v} = \dot m({C_{p,ave,out}}{T_{gap,ave,out}}  {C_{p,i}}{T_i}{\rm{) (W)}}\]
where the average temperature of the two outlet air streams is
\[{T_{gap,ave,out}}{\rm{ = }}({T_{gap,1,out}} + {T_{gap,2,out}})/2\]
and
C_{p,ave,out}_{ } = heat capacity of the outlet air evaluated at T_{gap,ave,out} (J/kgK)
Figure 108. Airflow window with betweenglass shading device showing variables used in the heat balance equations.
Evacuated Glazing Unit (EGU)[LINK]
Evacuated glazing Units (EGU) are an emerging technology, developed as a concept some 2030 years ago, but only now approaching widespread commercialization and adoption.
Evacuated glazing unit (EGU): cavity thermal conductance[LINK]
Thermal Conductance of the space in an evacuated glazing unit (EGU) is the sum of the conductance of the low pressure gas (air) and radiation conductance.
\[{C_{gas}} = {C_{COND}} + {C_{RAD}}\]
Conductance of the lowpressure gas
Conductance of low pressure gasses is calculated using formula by Corrucini (Corruccini, R. (1959)).
\[{C_{COND}} = \alpha \cdot \left[ {\frac{{\gamma + 1}}{{\gamma  1}}} \right] \cdot {\left[ {\frac{R}{{8 \cdot \pi \cdot M \cdot T}}} \right]^{\frac{1}{2}}} \cdot P\]
\[B = \alpha \cdot \left[ {\frac{{\gamma + 1}}{{\gamma  1}}} \right] \cdot {\left[ {\frac{R}{{8 \cdot \pi \cdot M \cdot T}}} \right]^{\frac{1}{2}}}\]
\[{C_{COND}} = B \cdot P\]
Where:
\[\alpha = \frac{{{\alpha_1}{\alpha_2}}}{{{\alpha_2} + {\alpha_1}\left( {1  {\alpha_2}} \right)}}\]
a_{1}, a_{2} = Accommodation coefficients of the gas molecules with the two surfaces. These values depend on the temperature, surface conditions, etc. For the present configuration and conditions, it is expected that a is approximately 0.5. If conservative value is needed than value of 1.0 could be used. With a_{1}, a_{2} = 0.5, a = 0.333
g = Specific heat ratio, g_{air} = 1.40.
R = Universal gas constant, R = 8,314.462175 J/mol·K
M = Molecular Weight, M_{air} = 28.97 [mol/g]
T = (T_{1} + T_{2})/2 [K]
P = Pressure of the gas [N/m^{2}]
From the paper Collins and Robinson (Collins, R., & Robinson, S. (1991)), B is set at approximately 50 for Air, if pressure is given in torr. Therefore according to Collins and Robinson, for air and approximate conditions of EGU:
C_{COND} ≈ 50·P
Where P is in torr (i.e., mm Hg).
Note: Conversion from Pa to torr is accomplished by multiplying value in torr by 133.28.
Using formula 2 and assuming T_{1} to be 20 ºC and T_{2} to be 18 ºC (expected temperatures of glass surfaces in EGU, if one glass surface is lowe and unit is exposed to NFRC standard environmental conditions), and using SI system of units, the B is calculated as 54.4, which is very close to the value of 50, proposed by Collins and Robinson
C_{COND} ≈ 54.4·P
We will use exact value, calculated by the formula, so values of will be input into the calculations, which enables more flexible model that can account for special treatment of glass surfaces.
Radiation conductance[LINK]
Radiation conductance for the two parallel plates is given by:
\[{C_{RAD}} = \frac{1}{{{\varepsilon_1}^{  1} + {\varepsilon_2}^{  1}  1}}\sigma \frac{{\left( {{T_1}^4  {T_2}^4} \right)}}{{{T_1}  {T_2}}}\]
Where:
e_{1} = emissivity of the first facing glass surface, [  ]
e_{2} = emissivity of the second facing glass surface, [  ]
s = StefanBoltzmann Constant, 5.67 x 10^{8}, [W/(m^{2}·K^{4})]
T_{1} = Temperature of the first facing glass surface, [K]
T_{2} = Temperature of the second facing glass surface. [K]
Assuming glass surface temperatures of 20 ºC and 18 ºC, respectively, the following radiation conductances are obtained for the three different glass emissivities:
Clear Glass (e_{1} = e_{2} = 0.84): C_{RAD} = 3.4 W/m^{2}K
Hard Coat Lowe (e_{1} or e_{2} = 0.15): C_{RAD} = 0.68 W/m^{2}K
Soft Coat Lowe (e_{1} or e_{2} = 0.04): C_{RAD} = 0.19 W/m^{2}K
Note: C_{RAD} of 0.09 is theoretically possible using best lowe technology today (i.e., e_{1} or e_{2} = 0.02).
Note: Lowe values above are typical values, which will vary by manufacturer. Some more recent hard coat lowe values are at or below 0.1.
It should also be noted that values above are based on the fix set of temperatures, while in reality temperatures will depend on the environmental conditions and surface emissivities (e.g., it cannot be expected that clear glass will have same T_{1} and T_{2} as lowe glass).
Equation is precise formulation for two parallel plates at the constant temperature. Simplified equation under these conditions is given in the form of:
\[{C_{RAD}} = 4\sigma \frac{1}{{{\varepsilon_1}^{  1} + {\varepsilon_2}^{  1}  1}}{T_m}^3\]
Where:
T_{m} = mean temperature, [K]
\[{T_m} = \frac{{\left( {{T_1} + {T_2}} \right)}}{2}\]
Calculation of the Ufactor
\[U = \frac{1}{{{R_{tot}}}} = \frac{1}{{{R_o} + 2 \cdot {R_{glass}} + {R_{gap}} + {R_i}}}\]
\[{R_{gap}} = \frac{1}{{{C_{gap}}}} = \frac{1}{{{C_{COND}} + {C_{RAD}}}}\]
\[{R_{glass}} = \frac{{{t_{glass}}}}{{{k_{glass}}}}\]
Where:
t_{glass} = glass thickness; [m]
k_{glass} = glass conductivity; k_{glass} = 1 W/(m·K)
R_{glass} = 0.003 m^{2}K/W (for 3 mm glass pane)
R_{o} ≈ 0.033 m^{2}K/W
R_{i} ≈ 0.14 m^{2}K/W
Ufactor of EGU without any pillars (pretending that this is possible) would be calculated using C_{RAD} only. From above radiation conductance calculations:
Clear Glass: U = 2.64 W/(m^{2}·K) [0.464 Btu/(hr·ft^{2}·ºF)]
Hard Coat Lowe: U = 0.62 W/(m^{2}·K) [0.109 Btu/(hr·ft^{2}·ºF)]
Soft Coat Lowe: U = 0.19 W/(m^{2}·K) [0.034 Btu/(hr·ft^{2}·ºF)]
Adding conductance of the air at 0.001 torr (C_{COND} = 0.08 W/(m^{2}·K), these values become:
Clear Glass: U = 2.66 W/(m^{2}·K) [0.468 Btu/(hr·ft^{2}·ºF)]
Hard Coat Lowe: U = 0.68 W/(m^{2}·K) [0.120 Btu/(hr·ft^{2}·ºF)]
Soft Coat Lowe: U = 0.27 W/(m^{2}·K) [0.048 Btu/(hr·ft^{2}·ºF)]
In contrast, the Ufactor of the same configuration with the air at atmospheric pressure will be (For the space width of 50 mm, C_{COND} ≈ 450 W/m^{2}K):
Clear Glass: U = 5.52 W/(m^{2}·K) [0.468 Btu/(hr·ft^{2}·ºF)]
Hard Coat Lowe: U = 5.52 W/(m^{2}·K) [0.120 Btu/(hr·ft^{2}·ºF)]
Soft Coat Lowe: U = 5.52 W/(m^{2}·K) [0.048 Btu/(hr·ft^{2}·ºF)]
It is clear that emissivity of the glass surface makes no difference, because of the dominant conductance of the air space. Also, it is worth noting that the Ufactor of such configuration is very close to the Ufactor of single glazing.
Evacuated glazing unit (EGU): Glass support element thermal conductance[LINK]
Glass panes in the EGU are separated by an array of small support elements. Typically, these support elements have a cylindrical shape and are often referred to as “pillars.” Typical geometry of the pillar is 0.51.0 mm diameter and 0.05 mm (50 mm) height. They are typically spaced 12 in. apart in a form of square or staggered matrix.
Calculation of the Ufactor[LINK]
The conductance of these elements can be measured or numerically modeled to determine accurate thermal performance. Approximate method also exists and is based on the combination of modeling and analytical work for the conduction through small cylinders in contact with infinite parallel plates with thickness much larger than cylinder height.
The following formula can be used to determine conductance of the single pillar, C_{p} (Collins and FisherCripps 1991):
\[{C_p} = \frac{{2ka}}{{1 + \frac{{2h}}{{\pi a}}}}\]
Where:
k = conductivity of glass, W/(m·K)
a = radius of the pillar (m)
h = pillar height, m
For the square array of support pillars (Collins and FischerCripps 1991) proposes the following formula for their conductance, C_{pa}:
\[{C_{pa}} = \frac{{2ka}}{{{\lambda ^2}\left( {1 + \frac{{2h}}{{\pi a}}} \right)}}\]
Where:
l = pillar spacing, m
This formula is approximate and does not include effect of the conductivity of the pillar, but it gives good approximation for common materials used in this technology, since conductivity of the pillar does not play substantial role for noninsulating pillars (where “noninsulating” would mean that conductivity of the pillar is equal or higher than the conductivity of the glass pane.
The Ufactor of the EGU with support pillars is then:
\[U = \frac{1}{{{R_{tot}}}} = \frac{1}{{{R_o} + 2 \cdot {R_{glass}} + {R_{gap}} + {R_i}}}\]
Where:
\[{R_{gap}} = \frac{1}{{{C_{gap}}}} = \frac{1}{{{C_{COND}} + {C_{RAD}} + {C_{pa}}}}\]
References[LINK]
Collins, R., & FischerCripps, A. 1991. “Design of Support Pillar Arrays in Flat Evacuated Windows.”. Aust. J. Phys.
Collins, R., & Robinson, S. 1991. “Evacuated Glazing”. Solar Energy. Vol. 47, No. 1, pp. 2738.
Corruccini, R. 1959. “Gaseous Heat Conduction at Low Pressures and Temperatures”. Vacuum. Vol. 78, pp.1929.
Thermal Performance of Deflected Insulated Glazing Unit (IGU)[LINK]
Introduction[LINK]
Deflection of insulated glazing unit (IGU) can result in thermal performance degradation or improvement due to the reduction or increase of gap space width. Convection of the gas fill is affected by changed gap space and due to modified convection pattern and shorter or longer thermal path at the center of the glazing unit can result in increased or decreased thermal performance. For the most part, Ufactor is mostly affected as a direct result of changed thermal performance; however note that solar heat gain through the window (SHGC) can also be affected because of the effect of inward flowing fraction of absorbed solar radiation, which is affected by thermal performance of the IGU.
Deflection in sealed IGU is caused by the difference in gas pressure in IGU gap vs. outdoor/indoor pressure. Indoor and outdoor pressure can be considered equal, since indoor building environment is in pretty good contact with outdoor environment. We will call this pressure an atmospheric pressure, P_{a}. The differences in pressure between atmospheric and gap pressure is due to several factors, listed here:
1) Difference in atmospheric pressure between IGU fabrication location and end use location
2) Difference in temperature during fabrication and actual operating conditions for the glazing. It should be noted that initial temperature can be higher than ambient temperature during fabrication process, due to elevated sealant temperatures, which can raise local temperatures within the IGU.
3) Unbalanced gas fill leakage through the sealants, resulting in lower gap pressure and inward deflection.
4) Wind or static load pressure
Effects 1 and 2 will be modeled using equations presented below, while effect 3 does not have credible mathematical model. However, cumulative deflection, resulting from all three effects can be measured in the field and its effect on thermal performance can be modeled by specifying center glazing deflection.
Wind or static load pressure effects on deflection is not included in this model at this time, but will be considered for future versions.
Recognizing that indoor and outdoor air pressure could be different, such as in hot box test environment, future plans for the extension of the model will include option to specify different values for indoor and outdoor pressure. Another future improvement to the model will also include linking certain air gaps with indoor or outdoor environment, meaning that respective pressures in linked spaces will be set to equal.
Mathematical Model[LINK]
Mathematical model described in detail here is based on the research work by Bernier and Bourret (1997) and Timoshenko and WoinowskyKrieger (1959). Bernier and Bourret (1997) of the Ecole Polytechnique Montréal adopted Timoshenko and WoinowskyKrieger (1959) model for calculating flat plate deflection subjected to the differential pressure field (static), while their original contribution was to develop correlations for changes in thermal performance, based on IGU deflection at the center of glazing location. In addition to adopting Bernier and Bourret (1997) model here, we have also developed model for calculating change in thermal performance of deflected units when this deflection is measured in the field. Therefore, the mathematical formulation, presented here is divided into two sections; 1) calculation of the deflection and resulting thermal performance caused by pressure and temperature effects and 2) calculation of the thermal performance of the IGU when the deflection is measured.
Calculation of the deflection and thermal performance caused by pressure and temperature effects[LINK]
If coordinate system is set as shown in Figure 109 and Figure 110, it is possible to calculate deflection distribution at each point of pane by using following equation:
Figure 109. Deflection Coordinate System  2D
\[{L_{D\left( i \right)}}\left( {x,y} \right) = \frac{{16\cdot \Delta {P_{\left( i \right)}}}}{{{\pi ^6}\cdot {D_{\left( i \right)}}}}\mathop \sum \limits_{m = 1,3,5 \ldots }^\infty \mathop \sum \limits_{n = 1,3,5 \ldots }^\infty \frac{{\sin \frac{{m\pi x}}{W}\sin \frac{{n\pi y}}{H}}}{{mn{{\left( {{{\left( {\frac{m}{W}} \right)}^2} + {{\left( {\frac{n}{H}} \right)}^2}} \right)}^2}}}\]
Where,
\[{D_{\left( i \right)}} = \frac{{E\cdot t_{\left( i \right)}^3}}{{12\cdot \left( {1  {\nu ^2}} \right)}}\]
Where,
E = Young’s modulus (7.2 x 10^{10}) [Force per unit Area; SI: Pa, IP: psi]
t = thickness of glazing pane [Length; SI: m, IP: in.]
n = poison’s ratio (0.22 for glass) [NonDimensional]
Figure 110. Deflection Coordinate System  3D
ΔP_{i} = P_{gap(i)}  P_{gap(i1)} (for ith pane) [Force per unit Area; SI: Pa, IP: psi]
ΔP_{i} = P_{gap(1)}  P_{a} (first pane) [Force per unit Area; SI: Pa, IP: psi]
ΔP_{i} = P_{a}  P_{gap(n1)} (last pane) [Force per unit Area; SI: Pa, IP: psi]
Where,
P_{a} = atmospheric pressure. [Force per unit Area; SI: Pa, IP: psi]
\[{P_{gap\left( i \right)}} = \frac{{{P_{ini}}{V_{ini\left( i \right)}}{T_{gap\left( i \right)}}}}{{{T_{ini}}{V_{gap\left( i \right)}}}}\]
Where,
P_{ini} = Initial pressure. Applies to all gaps as a single value (input data  measured or otherwise) [Force per unit Area; SI: Pa, IP: psi]
T_{ini} = Initial temperature. Applies to all gaps as a single value (input data  measured or otherwise) [Degree Temperature; SI: K, IP: R]
V_{ini(i)} = Initial volume for ith gap. [Length*Length*Length; SI: m^{3}, IP: in^{3}]
\[{V_{ini\left( i \right)}} = {L_i}\cdot W\cdot H\]
Where,
L_{i} = nondeflected glazing gap width (for ith gap) [Length; SI: m, IP: in.]
W = IGU width [Length; SI: m, IP: in.]
H = IGU height [Length; SI: m, IP: in.]
T_{gap(i)} = temperature of the gap between two glass panes (calculated using center of glazing thermal calculation algorithm, as described in ISO 15099 (ISO 2003). This value is first calculated using nondeflected state and is recalculated after the resulting deflection is calculated. This process is repeated until temperature at next iteration does not differ by more than 0.1 ºC
V_{gap(i)} = volume of the IGU gap in deflected state [Lentgh*Length*Length; SI: m^{3}, IP: in^{3}]
\[{V_{gap\left( i \right)}} = {V_{ini\left( i \right)}} + W\cdot H\cdot \left( {\overline {{L_{D,i}}}  \overline {{L_{D,i + 1}}} } \right)\]
Where,
\(\overline {{L_{D,i}}}\) is mean deflection value for ith pane. [Length; SI: m, IP: in.]
Deflection of each pane can be positive or negative and is done solely to establish reference. Current frame of reference is that positive deflection means that pane is deflecting towards left side, while negative deflection means that pane is deflecting towards right side . Whether the deflection is in the direction of reducing the gap width or increasing it, it will be the result of pressure difference, as described in . When pressure in the glazing unit is higher than surrounding environmental pressure, the deflection will be towards increasing gap width (i.e., ballooning), while the opposite situation will result in decreasing gap width (i.e., vacuuming)
Figure 111. Deflection Direction Convention
The important part of calculating deflection of the IGU is to determine mean deflection value for each glazing pane. Mean deflection value is used to calculate gap volume in deflected state . Mean deflection of glazing pane can be calculated by integrating :
\[\overline {{L_{D\left( i \right)}}} = \mathop \smallint \limits_{x = 0}^W \mathop \smallint \limits_{y = 0}^H \frac{{16\cdot \Delta {P_{\left( i \right)}}}}{{{\pi ^6}\cdot {D_{\left( i \right)}}}}\mathop \sum \limits_{m = 1,3,5 \ldots }^\infty \mathop \sum \limits_{n = 1,3,5 \ldots }^\infty \frac{{\sin \frac{{m\pi x}}{W}\sin \frac{{n\pi y}}{H}}}{{mn{{\left( {{{\left( {\frac{m}{W}} \right)}^2} + {{\left( {\frac{n}{H}} \right)}^2}} \right)}^2}}}\]
Which is identical with the following expression:
\[\overline {{L_{D\left( i \right)}}} = \frac{{16\cdot \Delta {P_{\left( i \right)}}}}{{{\pi ^6}\cdot {D_{\left( i \right)}}}}\mathop \sum \limits_{m = 1,3,5 \ldots }^\infty \mathop \sum \limits_{n = 1,3,5 \ldots }^\infty \frac{{\mathop \smallint \nolimits_{x = 0}^W \mathop \smallint \nolimits_{y = 0}^H \sin \frac{{m\pi x}}{W}\sin \frac{{n\pi y}}{H}}}{{mn{{\left( {{{\left( {\frac{m}{W}} \right)}^2} + {{\left( {\frac{n}{H}} \right)}^2}} \right)}^2}}}\]
and because integral of sin(x) is equal with –cos(x), above equation will become:
\[\overline {{L_{D\left( i \right)}}} = \frac{{16\cdot \Delta {P_{\left( i \right)}}}}{{{\pi ^6}\cdot {D_{\left( i \right)}}}}\mathop \sum \limits_{m = 1,3,5 \ldots }^\infty \mathop \sum \limits_{n = 1,3,5 \ldots }^\infty \frac{{\left( {1  cos\left( {m\pi } \right)} \right)\left( {1  cos\left( {n\pi } \right)} \right)}}{{mn{{\left( {{{\left( {\frac{m}{W}} \right)}^2} + {{\left( {\frac{n}{H}} \right)}^2}} \right)}^2}}}\]
Finally, because \(cos\left( {m\pi } \right)\) and \(cos\left( {n\pi } \right)\) values are always equal to 1 for the given range of m and n, above equation will became:
\[\overline {{L_{D\left( i \right)}}} = \frac{{16\cdot \Delta {P_{\left( i \right)}}}}{{{\pi ^6}\cdot {D_{\left( i \right)}}}}\mathop \sum \limits_{m = 1,3,5 \ldots }^\infty \mathop \sum \limits_{n = 1,3,5 \ldots }^\infty \frac{4}{{{m^2}{n^2}{\pi ^2}{{\left( {{{\left( {\frac{m}{W}} \right)}^2} + {{\left( {\frac{n}{H}} \right)}^2}} \right)}^2}}}\]
After calculating mean pane deflection the following equation is used to calculate mean gap width:
\[{L_{r\left( i \right)}} = L + \left( {\overline {{L_{D,\left( i \right)}}}  \overline {{L_{D,\left( {i + 1} \right)}}} } \right)\]
Where,
L_{r(i)} = Mean gap “i” width after incorporating glazing deflection. This mean gap width is used to recalculate thermal performance of deflected IGU.
\(\overline {{L_{D,i}}}\) * = mean glazing deflection for each pane “i”.*
Calculation of the deflection at the center of glazing and mean glazing deflection for each pane is an iterative process, where the initial temperature distribution is calculated for nondeflected state, then deflection is calculated based on this temperature distribution, new temperature distribution is calculated for this deflected state, then temperatures from previous iteration are compared to the current iteration and the process is repeated until the difference is no larger than 0.1 ºC.
At the end of calculations, program will calculate and return maximum deflection value for each pane (i.e., center of glazing deflection). If we label maximum deflection of each pane as L_{D(i),max}, we can calculate this value by substituting x = W/2 and y = H/2 in equation to determine deflection at the center point. Therefore,
\[{L_{D\left( i \right),max}} = \frac{{16\cdot \Delta {P_{\left( i \right)}}}}{{{\pi ^6}\cdot {D_{\left( i \right)}}}}\mathop \sum \limits_{m = 1,3,5 \ldots }^\infty \mathop \sum \limits_{n = 1,3,5 \ldots }^\infty \frac{{\sin \frac{{m\pi }}{2}\sin \frac{{n\pi }}{2}}}{{mn{{\left( {{{\left( {\frac{m}{2}} \right)}^2} + {{\left( {\frac{n}{2}} \right)}^2}} \right)}^2}}}\]
For glazing systems with more than two glazing layers, meaning multiple gas filled gaps, the deflection will be calculated for each glazing pane assuming that the pressure in a gap is independent from each other and calculated separately, unless spaces are “linked” together (e.g., stretched film middle glazing that has hole for equalizing pressure).
NonLinked Gaps in 3 or more glazing layer system:[LINK]
The procedure shown above generally applies to the 3 or more layer glazing system, with the exception that neighboring pressures are no longer P_{a}, but rather could be P_{a} on one side and P_{gap} on the other, or have P_{gap} on both sides, as shown in Figure 111 for gap “i”. Center of glazing thermal calculation will determine new temperature distribution, after deflection is calculated for each glazing and will be used to determine new P_{gap}, as per the procedure above.
Linked Gaps in 3 or more glazing layer system:[LINK]
When one or more gaps are linked together, their pressure is assumed to be identical (e.g., in triple glazing IGU P_{gap,1} = P_{gap,2}.) This pressure is calculated from temperatures of bounding glazing for linked gaps (e.g., for triple glazing IGU, glazing 1 and 3) and using neighboring pressures outside of those bounding glazing (e.g., for triple glazed IGU, P_{a} on both sides).
Note: This feature is not implemented in WINDOW 7.1. It is considered for future enhancements to the program.
Gap(s) Linked to Indoor or Outdoor Environment:[LINK]
If one or more glazing gaps are linked to either indoor or outdoor environment its pressure is fixed to P_{a}. In combination situations, such as two or more gaps linked together with one of them being linked to indoor or outdoor environment, they will all have fixed pressure of P_{a}.
Calculation of the thermal performance caused by measured deflection[LINK]
When deflection is measured, it is normally measured at the point of maximum deflection. Maximum deflection occurs at center of the IGU (at W/2 and H/2).
Measured value is typically gap width at the point of maximum deflection, which we can label L_{G(i)}. For ith measured gap the width is equal to:
\[{L_{G\left( i \right)}} = {L_{\left( i \right)}} + \left( {{L_{D\left( i \right),max}}  {L_{D\left( {i + 1} \right),max}}} \right)\]
Figure 112. Sketch of the nonsymetrically Deflected Glazing Panes
If we label ratio of mean deflection and maximum deflection as R_{(i)}, then:
\[{R_{(i)}} = \frac{{\overline {{L_{D(i)}}} }}{{{L_{D(i)}}_{,\max }}} = \sum\limits_{m = 1,3,5,...}^\infty {\sum\limits_{n = 1,3,5,...}^\infty {\frac{4}{{{m^2}{n^2}{\pi ^2}{{\left( {{{(\frac{m}{W})}^2} + {{(\frac{n}{H})}^2}} \right)}^2}}}} }\]
\($ Important thing to note is that ratios (\){R_{( i )}}$ ) for all gaps in glazing system are equal.
\[{R_{\left( 1 \right)}} = {R_{\left( 2 \right)}} = \ldots = {R_{\left( i \right)}} = \ldots = {R_{\left( {n  1} \right)}} = R\]
Replacing and into the following equation is obtained:
\[{L_{r\left( i \right)}} = {L_{\left( i \right)}} + R\left( {{L_{D\left( i \right),max}}  {L_{D\left( {i + 1} \right),max}}} \right)\]
Combining with we get the following expression for the mean gap width:
\[{L_{r\left( i \right)}} = {L_{\left( i \right)}} + R\left( {{L_{G\left( i \right)}}  {L_{\left( i \right)}}} \right)\]
Number of equations given in expression is equal to n1, where n is number of panes. Therefore, we need one more equation to complete the system of equations that would allow us to solve for all independent variables. To get the last equation we can rewrite in slightly different manner:
\[{L_{D\left( i \right),max}} = \frac{{\Delta {P_{\left( i \right)}}}}{{{D_{\left( i \right)}}}}\cdot K\]
Where coefficient K combines all constant terms, while \({D_{\left( i \right)}}\) is given by and \(\Delta {P_{\left( i \right)}}\) is calculated by, and . Summing over all deflections, \({L_{D\left( i \right),max}}\) the following equation is obtained:
\[\mathop \sum \limits_{i = 1}^n \frac{{{D_{\left( i \right)}}}}{K}\cdot {L_{D\left( i \right),max}} = \mathop \sum \limits_{i = 1}^n \Delta {P_{\left( i \right)}} = 0\]
Note that sum of all \(\Delta {P_{\left( i \right)}}\) is equal to zero since outside pressure is equal to inside. Therefore, the remaining equation that completes the set of equations is:
\[\mathop \sum \limits_{i = 1}^n {D_{\left( i \right)}}\cdot {L_{D\left( i \right),max}} = 0\]
Solving system of equations[LINK]
In order to solve system of equations we will present in slightly different manner:
\[0 = {L_{\left( i \right)}}  {L_{G\left( i \right)}} + {L_{D\left( i \right),max}}  {L_{D\left( {i + 1} \right),max}}\]
Which in developed form will look like this:
\[0 = {L_{\left( 1 \right)}}  {L_{G\left( 1 \right)}} + {L_{D\left( 1 \right),max}}  {L_{D\left( 2 \right),max}}\]
\[0 = {L_{\left( 2 \right)}}  {L_{G\left( 2 \right)}} + {L_{D\left( 2 \right),max}}  {L_{D\left( 3 \right),max}}\]
\[0 = {L_{\left( i \right)}}  {L_{G\left( i \right)}} + {L_{D\left( i \right),max}}  {L_{D\left( {i + 1} \right),max}}\]
\[0 = {L_{\left( {n  1} \right)}}  {L_{G\left( {n  1} \right)}} + {L_{D\left( {n  1} \right),max}}  {L_{D\left( n \right),max}}\]
In order to express each \({L_{D\left( i \right),max}}\) as dependence from \({L_{D\left( n \right),max}}\) (deflection of inside/last pane) we will need to make sum from first to last, then from second to last, third to last and so on. This procedure will create following set of equations:
\[{L_{D\left( 1 \right),max}} = \mathop \sum \limits_{k = 1}^{k = n  1} \left( {{L_{G\left( k \right)}}  {L_{\left( k \right)}}} \right) + {L_{D\left( n \right),max}}\]
\[{L_{D\left( 2 \right),max}} = \mathop \sum \limits_{k = 2}^{k = n  1} \left( {{L_{G\left( k \right)}}  {L_{\left( k \right)}}} \right) + {L_{D\left( n \right),max}}\]
\[{L_{D\left( i \right),max}} = \mathop \sum \limits_{k = i}^{k = n  1} \left( {{L_{G\left( k \right)}}  {L_{\left( k \right)}}} \right) + {L_{D\left( n \right),max}}\]
\[{L_{D\left( {n  1} \right),max}} = \mathop \sum \limits_{k = n  1}^{k = n  1} \left( {{L_{G\left( k \right)}}  {L_{\left( k \right)}}} \right) + {L_{D\left( n \right),max}}\]
Now replacing this set of equations back to :
\[\mathop \sum \limits_{i = 1}^{n  1} {D_{\left( i \right)}}\cdot \left( {\mathop \sum \limits_{k = i}^{k = n  1} \left( {{L_{G\left( k \right)}}  {L_{\left( k \right)}}} \right) + {L_{D\left( n \right),max}}} \right) + {D_{\left( n \right)}}\cdot {L_{D\left( n \right),max}} = 0\]
Which solving by \({L_{D\left( n \right),max}}\) leads to the following equation:
\[{L_{D\left( n \right),max}} = \frac{{\mathop \sum \nolimits_{i = 1}^{n  1} ({D_{\left( i \right)}}\cdot \mathop \sum \nolimits_{k = i}^{k = n  1} \left( {{L_{\left( k \right)}}  {L_{G\left( k \right)}}} \right))}}{{\mathop \sum \nolimits_{i = 1}^n {D_{\left( i \right)}}}}\]
Calculating \({L_{D\left( n \right),max}}\) value from this equation and substituting it in will enable calculation of the deflection of remaining panes.
References[LINK]
Arasteh, D.K., J.C. Kohler and B.T. Griffith. Draft 2009. Modeling Windows in EnergyPlus with only U, SHGC, and optionally VT. LBNL report. Full reference to be determined.
Arasteh, D. J. 2009. Modeling Windows in EnergyPlus with only U, SHGC, and optionally VT. Lawrence Berkeley National Laboratory.
Arasteh, D.K., M.S. Reilly and M.D. Rubin. 1989. A versatile procedure for calculating heat transfer through windows. American Society of Heating, Refrigeration and AirConditioning Engineers, ASHRAE Transactions, Vol. 95, Pt. 2.
Bernier, M., & Bourret, B. January 1997. “Effects of Glass Plate Curvature on the UFactor of Sealed Insulated Glazing Units”. Atlanta, GA: ASHRAE Transactions. Vol. 103, Pt. 1. American Society for Heating, Refrigerating and AirConditioning Engineers.
Collins, R., & FischerCripps, A. 1991. “Design of Support Pillar Arrays in Flat Evacuated Windows.”. Aust. J. Phys.
Collins, R., & Robinson, S. 1991. “Evacuated Glazing”. Solar Energy. Vol. 47, No. 1, pp. 2738.
Corruccini, R. (1959). “Gaseous Heat Conduction at Low Pressures and Temperatures”. Vacuum. Vol. 78, pp.1929.
Finlayson, E.U., D.K. Arasteh, C. Huizenga, M.D. Rubin and M.S. Reilly. 1993. WINDOW 4.0: documentation of calculation procedures. Lawrence Berkeley National Laboratory report no. LBL33943.
ISO. 2003. ISO 15099:2003. Thermal performance of windows, doors, and shading devices – Detailed calculations. International Organization for Standardization.
Klems, J. H. 1994A. “A New Method for Predicting the Solar Heat Gain of Complex Fenestration Systems: I. Overview and Derivation of the Matrix Layer Calculation.”. ASHRAE Transactions. 100 (pt.1): 10731086.
Klems, J. H. 1994B. “A New Method for Predicting the Solar Heat Gain of Complex Fenestration Systems: II. Detailed Description of the Matrix Layer Calculation.”. ASHRAE Transactions. 100 (pt.1): 10731086.
Klems, J. H. 1995. “Measurements of Bidirectional Optical Properties of Complex Shading Devices.”. ASHRAE Transactions. 101 (pt 1; Symposium Paper CH9581 (RP548)): 791801.
Klems, J. H. 1996. “A Comparison between Calculated and Measured SHGC for Complex Glazing Systems.”. ASHRAE Transactions. 102 (Pt. 1; Symposium Paper AT96161): 931939.
Klems, J. H. 1996. “Calorimetric Measurements of InwardFlowing Fraction for Complex Glazing and Shading Systems.”. ASHRAE Trans. 102(Pt. 1; Symposium Paper AT96163): 947954.
Munson, B.R, D.F. Young and T.H. Okiishi. 1998. “Fundamentals of Fluid Mechanics,” Third Edition Update, John Wiley & Sons, Inc.
Papamichael, K. J. 1998. “Determination and Application of Bidirectional SolarOptical Properties of Fenestration Systems.”. Cambridge, MA: 13th National Passive Solar Conference.
Simmler, H., U. Fischer and F. Winkelmann. 1996. SolarThermal Window Blind Model for DOE2. Lawrence Berkeley National Laboratory, Simulation Research Group internal report, (unpublished).
Timoshenko, S., & KreigerWoinowsky, S. 1959. “Theory of Plates and Shells” 2nd Edition. McGrawHill Company.
Equivalent Layer Fenestration Model[LINK]
The section describes the equivalent layer fenestration optical and thermal model. The Equivalent Layer fenestration model can have four types of attachments: drapes, venetian blinds, roller blinds and insect screens. In this model shading layers are assumed to be uniform and can be represented by an equivalent homogenous layer that has spatiallyaveraged “effective” optical and thermal properties (ASHRAE 1311RP). Likewise, venetian blinds can be characterized using effective optical and thermal properties. When solar radiation strikes a window surface some fraction of the incident solar radiation passes unobstructed through openings in a shading layer and the remaining fraction is intercepted by the structure of the layer. The intercepted radiation is partly absorbed, partly reflected and partly transmitted. These reflected and transmitted components of the scattered solar radiation are assumed to be uniformly diffuse. Shading layers, because of their openness, generally transmit longwave radiation, and the effective infrared properties of shades account for that. Using effective optical properties and a beam/diffuse split of solar radiation at each layer, the equivalent layer approach can represent multilayer systems. This representation provides virtually unlimited flexibility to combine different types of shading layers in a fenestration. The equivalent layer window model requires a few set of optical data to characterize a particular layer and this set of data is used to calculate effective layer properties. For instance, the effective solar optical properties of a venetian blind can be calculated as a function of slats optical properties and geometry. Also, it is possible to adjust slat angle at each time step in response to the changing angular position of the sun. Moreover, the model provides control strategies as a function of slat angle that can be changed at each time step as needed. Likewise, effective properties of a pleated drape are calculated as a function of fabric properties and a specified value of fullness. The only input data needed to fully characterize drapery fabrics, roller blinds and insect screens are material openness as area fraction, and the transmittance and reflectance at normal incidence. Shade openness area fraction is the same as the beambeam transmittance at normal incidence. In multilayer fenestration, each layer is separated by a gap. A gap in equivalent layer model is defined by specifying the fill gas and the gap spacing. Currently five gas types are allowed: Air, Argon, Xenon, Krypton and Custom. The convective heat transfer coefficient in a gap is calculated depending on the spacing, the temperatures of the layers and the fill gas properties. Equivalentlayer concept – offers wide range of multiple glazing and shading layers combination and can simulate multilayer complex fenestration systems. The effective layer properties of venetian blinds, pleated drapes, roller blinds, and insect screens are calculated from geometric layer models and material properties. A set of empirical correlations for estimating offnormal material properties were developed under ASHRAE research project (ASHRAE 1311RP).
The Equivalent Layer Analysis[LINK]
The equivalent layer windows system is treated as a series of parallel layers separated by gaps as shown in Figure 113. This multilayer structure has been used in several computer programs and the underlying theory has been documented (ASHRAE 1311RP).
Figure 113. Multilayer fenestration analysis structure (ASHRAE 1311RP)
Equivalent layer fenestration model uses twostep analysis. First, the flux of absorbed solar radiation at each layer, S_{i}, caused by the incident radiation flux, I_{sol}, is determined using net radiation analysis. Second, an energy balance is applied at each layer, accounting for heat transfer and the known set of absorbed solar radiations S_{i} values, in order to solve for the set of layer temperatures, T_{i}, and the corresponding heat fluxes. The fenestration model also accounts for the diathermanous shade layers in the longwave radiant exchange. The latter can be significant for shading layers. Glass is considered opaque with respect to longwave radiation.
Equivalent Layer Optical Model[LINK]
The multilayer optical model is based on an algorithm originally developed by Edwards (1977) and extended by Wright and Kotey (2006). The algorithm models the interaction of incident solar radiation with a glazing system composed of any number of parallel, planar, specular glazing layers. The shading layers scatter portion of the incident solar radiation diffusely, and the model tracks the beam and diffuse components of solar radiation as they interact with a multilayer system of glazing and shading layers. The conceptual arrangement for tracking beam and diffuse solar flux components is illustrated in Figure 114.
Analysis yields beambeam, beamdiffuse and diffusediffuse fluxes, providing full detail concerning the quantities of reflected, transmitted and absorbed radiation.
Figure 114. Solar analysis of the multilayer glazing/shading system showing beam and diffuse fluxes (ASHRAE 1311RP)
The beam and diffuse characterization of solar radiation demands an expanded set of solar optical properties (Wright and Kotey 2006). The quantities of interest for single layer are shown inFigure 115.
Where,
t = transmittance of a glazing or a shading layer, ()
r = reflectance of a glazing or a shading layer, ()
e = thermal emissivity of a glazing or a shading layer, ()
q = incident angle, (degree)
f = slat angle, (degree)
W = profile angle, (degree)
J = radiosity of a layer surface, (W/m2)
T = temperature of a layer, (C)
Subscripts
f = front side of a layer
b = back side of a layer
bb = beambeam optical property. Represents a fraction of the beam radiation incident at a given layer and leaves the layer without being scattered.
bd = beamdiffuse optical property. Represents a fraction of the beam radiation incident at a given layer and leaves the layer diffusely
dd = diffusediffuse optical property. Represents a diffuse radiation incident at a given layer and leaves the layer as diffuse
i = layer index, ()
Superscripts
m = represents material property (e.g., fabric material)
str = represents an apparent material structure property of roller blinds
w = represents an apparent wire property
Figure 115. Twelve solar properties assigned at each layer (ASHRAE 1311RP)
Each glazing or shading layer, in general, require eleven set of solar properties (Wright and Kotey 2006). Only beambeam properties are needed for common glazing layers. Each shading layer is represented by a few set of component properties and geometries that are converted to the full set of layer optical and thermal properties by equivalent layer models (ASHRAE 1311RP). The models also make account for offnormal incidence of solar radiation and can calculate the diffuse component from normal incidence values using numerical integration.
Equivalent Layer Thermal Model[LINK]
A surface energy balance is applied at each layer and the resulting set of equation is solved for layer temperatures and heat transfer fluxes. A schematic drawing of the multilayer illustration and variables designation is shown in Figure 116. The net radiation formulation based on the radiosities, Jf,i and Jb,i the radiant flux leaving the front and back surfaces of the i^{th} layer, respectively, is the used as a solution technique. The net radiant heat flux across a gap can be expressed as the difference between the radiosities of the bounding surfaces. The net radiation model analysis yields the layer temperatures and corresponding heat transfer coefficients that are used to determine Ufactor and SHGC. The thermal analysis is done in “ASHWAT_Thermal” routine.
Figure 116. Radiosity model used in thermal analysis of the multilayer glazing/shading system (ASHRAE 1311RP)
Glazing Layer Analysis[LINK]
The EquivalentLayer window model allows a wide range of options regarding the selection of glazing layers and gaps combinations. It also allows mixing glazing and shading layers in any sequence. For instance, it allows placing a venetian blind or a roller blind in between glazing layers. Only beambeam solar properties at normal incidence, infrared transmittance and infrared emissivity are required for common glazing types. Offnormal solar properties at a given sun position are estimated by adjusting the normal incidence values in “ASHWAT_OffNormalProperties” routine. The model assumes that the ratio between normal and offnormal transmittance is the same for the glazing layer in question and a reference piece of uncoated 6 mm glass with a moderate tint (ASHRAE 1311RP). A similar procedure is used to convert solar reflectance from normal to offnormal.
Gap Layer Analysis[LINK]
Gaps can be specified as sealed, vented indoors or vented outdoors. The sealed gap is considered as enclosed cavity. Vented gaps is specified for inner and outer most gaps only and when the gas type is AIR. The vented gap model assumes that air flow at the perimeter of the window is not restricted. Any fill gas can be specified by molecular mass and thermophysical properties of viscosity, specific heat and thermal conductivity. This is done by providing “a”, “b” and “c” coefficients of the quadratic equation of the form, p = a + bT + cT^{2}, where the T (K) is temperature and p is the property being evaluated. EnergyPlus has builtindata available for common fill gas components including Air, Argon, Krypton and Xenon. Users specify CUSTOM gas by defining the coefficients as an input.
Shade Layer Analysis[LINK]
Equivalent layer fenestration model has a complete set of solar and longwave models for the four shading layer types: drapes, venetian blinds, roller blinds and insect screens developed under ASHRAE 1311RP. And semiempirical models are used to evaluate the offnormal properties of drape, roller blind and insect screen materials. The effective layer properties of venetian blinds and the effect of pleating in drapes are determined using a more fundamental net radiation scheme. The openness fraction, A_{o} and beambeam solar transmittance at normal incidence, t_{bb} (q = 0) really represent a geometric quantity and it has been confirmed that they can be used interchangeably (ASHRAE 1311RP). Openness is simply the fraction of a material, by area, that is open. In equivalent layer fenestration model, beambeam transmittance at normal incidence (i.e., openness fraction), beamdiffuse transmittance and reflectance at normal incidence are required to characterize drapery fabric, roller blind and insect screen material. A conventional venetian blind in equivalent layer model can be characterized by specifying the geometry, solar reflectance and transmittance, and emissivity of the slats. The offnormal solar properties of drape, roller blind and insect screen materials were formulated based on measurements (ASHRAE 1311RP). The longwave properties of the of the drape fabrics, roller blinds and insect screens can be calculated from the material emissivity and openness of the fabric (ASHRE 1311RP). The optical model development for venetian blinds is presented by Yahoda and Wright (2004 and 2005), and Kotey et al. (2008).
Drapes and Curtains[LINK]
The offnormal optical properties of drapery fabric is determined from user specified optical properties at normal incidence (q = 0) using empirical correlations (Kotey et al., 2009a). The input data required to characterize drapery fabric are the beambeam transmittance of the material at normal incidence t_{bb} (q = 0) = Ao, the beamdiffuse transmittance of the material at normal incidencet_{bd} (q = 0), and the beamdiffuse (total) reflectance of the materialr_{bt} (q = 0). The offnormal properties are calculated as follows.
Offnormal Transmittance:
\[{\rm{\tau }}_{{\rm{bb}}}^{\rm{m}}\left( {\rm{\theta }} \right) = {\rm{\tau }}_{{\rm{bb}}}^{\rm{m}}{\rm{(\theta }} = {\rm{0)}} \cdot {\rm{co}}{{\rm{s}}^{\rm{b}}}\left( {\rm{\theta }} \right)\]
\[{\rm{\tau }}_{{\rm{bt}}}^{\rm{m}}\left( {\rm{\theta }} \right) = {\rm{\tau }}_{{\rm{bt}}}^{\rm{m}}{\rm{(\theta }} = {\rm{0)}} \cdot {\rm{co}}{{\rm{s}}^{\rm{b}}}\left( {\rm{\theta }} \right)\]
\[{\rm{b}} = {\left[ {\;  0.5 \cdot \ln \;\left( {\;{{\left[ {{\rm{\tau }}_{{\rm{bb}}}^{\rm{m}}{\rm{(\theta }} = {\rm{0)}}\;{\rm{,}}\;\;{\rm{0}}{\rm{.01}}} \right]}^{\;\max }}\;} \right)\;\;,\;\;0.35\;} \right]^{\;\max }}\]
\[{\rm{b}} = {\left[ {\;  0.5 \cdot \ln \;\left( {\;{{\left[ {{\rm{\tau }}_{{\rm{bt}}}^{\rm{m}}{\rm{(\theta }} = {\rm{0)}}\;{\rm{,}}\;\;{\rm{0}}{\rm{.01}}} \right]}^{\;\max }}\;} \right)\;\;,\;\;0.35\;} \right]^{\;\max }}\]
\[{\rm{\tau }}_{{\rm{bt}}}^{\rm{m}}\left( {\rm{\theta }} \right) = {\rm{\tau }}_{{\rm{bd}}}^{\rm{m}}\left( {\rm{\theta }} \right) + {\rm{\tau }}_{{\rm{bb}}}^{\rm{m}}\left( {\rm{\theta }} \right)\]
The offnormal reflectance:
\[{\rm{\rho }}_{{\rm{bt}}}^{\rm{m}}\left( {\rm{\theta }} \right) = {\rm{\rho }}_{{\rm{bd}}}^{\rm{m}}\left( {\rm{\theta }} \right) = {\rm{\rho }}_{{\rm{bt}}}^{\rm{m}}\left( {{\rm{\theta }} = {\rm{0}}} \right) + \left( {{\rm{\rho }}_{{\rm{bt}}}^{\rm{m}}\left( {{\rm{\theta }} = 9{\rm{0}}^\circ } \right)  {\rm{\rho }}_{{\rm{bt}}}^{\rm{m}}\left( {{\rm{\theta }} = {\rm{0}}} \right)} \right) \cdot \left( {1  {\rm{co}}{{\rm{s}}^{{\rm{0}}{\rm{.6}}}}\left( {\rm{\theta }} \right)} \right)\]
\[{\rm{\rho }}_{{\rm{bt}}}^{\rm{m}}\left( {{\rm{\theta }} = {\rm{90}}^\circ } \right) = {\rm{\rho }}_{{\rm{bt}}}^{\rm{m}}\left( {{\rm{\theta }} = {\rm{0}}} \right) + \left( {{\rm{1}}  {\rm{\rho }}_{{\rm{bt}}}^{\rm{m}}\left( {{\rm{\theta }} = {\rm{0}}} \right)} \right)\;\left( {\,0.7{{\left( {\;{\rm{\rho }}_{}^{\rm{y}}} \right)}^{\;0.7}}} \right)\]
The apparent yarn reflectance is given by:
\[{\rho ^y} = \frac{{\rho_{{\rm{bt}}}^{\rm{m}}\left( {{\rm{\theta }} = {\rm{0}}} \right)}}{{{\rm{1}}  \tau_{{\rm{bb}}}^{\rm{m}}\left( {{\rm{\theta }} = {\rm{0}}} \right)}} = \frac{{\rho_{{\rm{bt}}}^{\rm{m}}\left( {{\rm{\theta }} = {\rm{0}}} \right)}}{{{\rm{1}}  {{\rm{A}}_{\rm{o}}}}}\]
The above set of equations for drapery fabrics are used subject to the condition that the solar absorptance of the fabric, at normal incidence, is not less than 1% (ASHRE 1311RP). The diffusediffuse material properties, for Equivalent layer window model, are determined using Rhomberg integration with 32 panels covering the range from q = 0 to q = 90° (ASHRAE 1311RP). The subscript “X” stands for either front or back side of the layer.
\[\tau_{{\rm{x,dd}}}^{\rm{m}} = \;\;{\rm{2}}\;\int\limits_{{\rm{\theta }} = {\rm{0}}}^{{\raise0.5ex\hbox{\scriptstyle {\rm{\pi }}}\kern0.1em/\kern0.15em\lower0.25ex\hbox{\scriptstyle {\rm{2}}}}} {\tau_{{\rm{x,bt}}}^{\rm{m}}\left( {\rm{\theta }} \right) \cdot {\rm{cos}}\left( {\rm{\theta }} \right) \cdot {\rm{sin}}\left( {\rm{\theta }} \right)} \;{\rm{d\theta }}\]
\[\rho_{{\rm{x,dd}}}^{\rm{m}} = \;\;{\rm{2}}\;\int\limits_{{\rm{\theta }} = {\rm{0}}}^{{\raise0.5ex\hbox{\scriptstyle {\rm{\pi }}}\kern0.1em/\kern0.15em\lower0.25ex\hbox{\scriptstyle {\rm{2}}}}} {\rho_{{\rm{x,bt}}}^{\rm{m}}\left( {\rm{\theta }} \right) \cdot {\rm{cos}}\left( {\rm{\theta }} \right) \cdot {\rm{sin}}\left( {\rm{\theta }} \right)} \;{\rm{d\theta }}\]
The above set of equations for drapery fabrics apply to the full range of A_{o}, fabric transmittance and fabric reflectance including that falls within the bounds of Keyes’ (1967) fabric chart plus sheer fabrics (ASHRAE 1311RP). The longwave thermal emissivity and thermal transmittances of drapery fabric are calculated using the following correlations and fabric openness fraction (Kotey et al. 2008).
\[{\varepsilon ^{\rm{m}}} = {\rm{0}}{\rm{.87}} \cdot \left( {{\rm{1  }}{{\rm{A}}_{\rm{o}}}} \right)\]
\[{\tau ^{\rm{m}}} = {\rm{1}}  {\rm{0}}{\rm{.95}} \cdot \left( {{\rm{1}}  {{\rm{A}}_{\rm{o}}}} \right)\]
The optical and thermal properties determined using the above same sets of equations are equally valid for pleated drape shades (Kotey, et. al., 2009a). For pleated drape, the effective beambeam and beamdiffuse solar properties are determined by tracking both radiation components, for a given incident angle, and interaction with a fabric pleated rectangular geometry shown in Figure 117. The solar optical properties of the two different pleat planes are evaluated on the basis of the local solar incidence angle. Therefore, the effective layer properties are influenced not just by horizontal solar profile angle, W_{H}, but also by incidence angle (ASHRAE 1311RP).
Figure 117. Geometry used for Pleated Drape Analysis
The solar diffusediffuse and longwave effective properties of the pleated drape are evaluated using a much simpler netradiation analysis using conventional shape factors (Kotey, et. al., 2009a). Users can chose and apply the pleated drape model to any fabric and any degree of fullness (ASHRAE 1311RP).
Venetian Blinds[LINK]
The effective shortwave optical and longwave optical properties of venetian blind layer is estimated using analytical models (Yahoda and Wright 2004, 2005; Kotey et al. 2008). The model requires properties of venetian blind slats and geometry of the slats shown in Figure 1118.
Figure 1118. Likewise, the effective longwave properties are obtained for the layer knowing longwave properties of the slats.
Figure 118. Geometry and properties used for venetian blind analysis
The model assumes that venetian blind slats reflect and transmit solar radiation diffusely (ASHRAE 1311RP). The same assumption is made regarding thermal radiation. The input data required to characterize a venetian blind are: front and back side reflectance and transmittance of the slat, geometry and infrared emissivity and transmittance of the slate. The effective optical properties of the venetian blind are determined by tracking beam and diffuse solar radiation through various interactions with slats positioned at a given slat angle. The model uses simple foursurface model if the slats are fully sunlit and a sixsurface model if the slats are partially sunlit (ASHRAE 1311RP). Slats are assumed to be thin and flat but a correction is applied to account for slat curvature effect (Kotey et al. 2008).
Roller Blinds[LINK]
The offnormal properties of rollerblind are determined from solar properties of roller blind fabric measured at normal incidence (q = 0) using correlations (Kotey, et. al., 2009b). The offnormal properties for roller blind shades are calculated using the set equations given below:
\[{{\rm{\tau }}_{bb}}\left( {\rm{\theta }} \right) = \left\{ \begin{array}{l}{{\rm{\tau }}_{bb}}{\rm{(\theta }} = {\rm{0)}} \cdot {\rm{co}}{{\rm{s}}^{\rm{b}}}{\kern 1pt} \left( {\frac{{\rm{\theta }}}{{{{\rm{\theta }}_{{\rm{cutoff}}}}}} \cdot \frac{\pi }{2}} \right)\quad \theta < {\theta_{CutOff}}\\0\quad \quad \quad \theta \ge {\theta_{CutOff}}\end{array} \right.\]
\[{\rm{b}} = {\rm{0}}{\rm{.6}} \cdot {\rm{co}}{{\rm{s}}^{{\rm{0}}{\rm{.3}}}}\left( {{{\rm{A}}_{\rm{o}}} \cdot \frac{{\rm{\pi }}}{{\rm{2}}}} \right)\]
\[{{\rm{\theta }}_{{\rm{cutoff}}}} = {\rm{65}}^\circ + \;\left( {{\rm{90}}^\circ  {\rm{65}}^\circ } \right) \cdot \left( {{\rm{1}}  {\rm{cos}}\left( {{{\rm{A}}_{\rm{o}}} \cdot \frac{{\rm{\pi }}}{{\rm{2}}}} \right)} \right)\]
\[{\rm{\tau }}_{{\rm{bt}}}^{}\left( {\rm{\theta }} \right) = {\rm{\tau }}_{{\rm{bt}}}^{}{\rm{(\theta }} = {\rm{0)}} \cdot {\rm{co}}{{\rm{s}}^{\rm{b}}}{\kern 1pt} \left( {\rm{\theta }} \right)\]
\[b = \left\{ \begin{array}{lc} 0.133 \cdot \left( \rm{\tau}^{str} + 0.003 \right)^{0.467} &\text{if} \left( 0 \le \rm{\tau}^{str} \le 0.33 \right) \\ 0.33 \cdot \left( 1  \rm{\tau}^{str} \right) &\text{if} \left( 0.33 < \rm{\tau}^{str} \le 1 \right) \end{array} \right.\]
\[{\rm{\tau }}_{}^{{\rm{str}}} = \frac{{{\rm{\tau }}_{{\rm{bt}}}^{}\left( {{\rm{\theta }} = {\rm{0}}} \right)  {{\rm{\tau }}_{{\rm{bb}}}}\left( {{\rm{\theta }} = {\rm{0}}} \right)}}{{{\rm{1}}  {{\rm{\tau }}_{{\rm{bb}}}}\left( {{\rm{\theta }} = {\rm{0}}} \right)}}\]
\[{\rm{\tau bt}}\left( {\rm{\theta }} \right) = {\rm{\tau bd}}\left( {\rm{\theta }} \right) + {\rm{\tau bb}}\left( {\rm{\theta }} \right)\]
The offnormal solar property calculation of roller blind is based on a set of correlations developed from measurement data using samples of commonly used commercially produced roller blind material openness range of 0.0 – 0.14. Thus, these correlations are not applicable for shades with materials openness fraction exceeding 0.20. The mean solar reflectance of a roller blind material was found to be purely diffuse and unaffected by incidence angle and is given by:
\[{\rm{\rho }}_{{\rm{bt}}}^{}\left( {{\rm{\theta }} = {\rm{0}}} \right) = {\rm{\rho }}_{{\rm{bd}}}^{}\left( {{\rm{\theta }} = {\rm{0}}} \right) = {\rm{\rho }}_{{\rm{bt}}}^{}\left( {\rm{\theta }} \right) = {\rm{\rho }}_{{\rm{bd}}}^{}\left( {\rm{\theta }} \right)\]
The diffusediffuse transmittance and reflectance are obtained by Rhomberg numerical integration. The longwave properties of roller blind material determined using the material property and the openness fraction (Kotey et al. 2008) as shown below:
\[{\varepsilon ^{{\rm{rb}}}} = {\rm{0}}{\rm{.91}} \cdot \left( {{\rm{1  }}{{\rm{A}}_{\rm{o}}}} \right)\]
\[{\tau ^{{\rm{rb}}}} = {\rm{1}}  {\rm{0}}{\rm{.95}} \cdot \left( {{\rm{1}}  {{\rm{A}}_{\rm{o}}}} \right)\]
Insect Screens[LINK]
The empirical correlations formulated to obtain the effective offnormal solar and longwave properties of insect screens were based on measurements (Kotey et al. (2009a). Insect screen geometry is shown in Figure 119. The calculation of effective solar properties requires a set of properties measured at normal incidence: t_{bb} (q = 0), t_{bd} (q = 0) and r_{bt} (q = 0).
Figure 119. Geometry used for insect screen analysis
Openness can be determined by optical measurement at normal incidence, A_{o} = t_{bb} (q = 0), but in the case of insect screens A_{o} can reliably be calculated knowing wire diameter (d), and wire spacing (s) as follows:
\[{A_o} = {\left( {1  {d \mathord{\left/ {\vphantom {d s}} \right. } s}} \right)^2}\]
The incidence angle beyond which direct beam transmission is cut off, q_{CutOff}, can also be estimated purely from geometry and is given by:
\[{\theta_{CutOff}} = co{s^{  1}}\left( {{d \mathord{\left/ {\vphantom {d S}} \right. } S}} \right)\]
The offnormal properties are calculated as follows.
\[{\tau_{bb}}\left( {\rm{\theta }} \right) = \left\{ \begin{array}{l}{\tau_{bb}}{\rm{(\theta }} = {\rm{0)}} \cdot {\rm{co}}{{\rm{s}}^{\rm{b}}}\left( {\frac{{\rm{\theta }}}{{{{\rm{\theta }}_{{\rm{cutoff}}}}}} \cdot \frac{\pi }{{\rm{2}}}} \right)\quad {\rm{\theta }} < {{\rm{\theta }}_{{\rm{cutoff}}}}\\0\quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \;\,{\rm{\theta }} \ge {{\rm{\theta }}_{{\rm{cutoff}}}}\end{array} \right.\]
\[{\rm{b}} =  {\rm{0}}{\rm{.45}}\;{\rm{ln}}\;\left( {{\rm{MAX}}\left( {{{\rm{\tau }}_{{\rm{bb}}}}\left( {{\rm{\theta }} = {\rm{0}}} \right){\rm{,}}\;{\rm{0}}{\rm{.01}}} \right)} \right)\; + \;{\rm{0}}{\rm{.1}}\]
\[{\rm{\tau }}_{{\rm{bt}}}^{}\left( {\rm{\theta }} \right) = {\rm{\tau }}_{{\rm{bt}}}^{}{\rm{(\theta }} = {\rm{0)}} \cdot {\rm{co}}{{\rm{s}}^{\rm{b}}}\left( {\rm{\theta }} \right)\]
\[{\rm{b}} =  {\rm{0}}{\rm{.65}}\;{\rm{ln}}\;\left( {{\rm{MAX}}\left( {{{\rm{\tau }}_{{\rm{bt}}}}\left( {{\rm{\theta }} = {\rm{0}}} \right){\rm{,}}\;{\rm{0}}{\rm{.01}}} \right)} \right)\; + \;{\rm{0}}{\rm{.1}}\]
\[{\rm{\tau }}_{{\rm{bd}}}^{}\left( {\rm{\theta }} \right) = {\rm{\tau }}_{{\rm{bt}}}^{}\left( {\rm{\theta }} \right)  {\rm{\tau }}_{{\rm{bb}}}^{}\left( {\rm{\theta }} \right)\]
\[{\rm{\rho }}_{{\rm{bt}}}^{}\left( {\rm{\theta }} \right) = {\rm{\rho }}_{{\rm{bt}}}^{}\left( {{\rm{\theta }} = {\rm{0}}} \right) + \left( {{\rm{\rho }}_{{\rm{bt}}}^{}\left( {{\rm{\theta }} = 9{\rm{0}}^\circ } \right)  {\rm{\rho }}_{{\rm{bt}}}^{}\left( {{\rm{\theta }} = {\rm{0}}} \right)} \right) \cdot \left( {1  {\rm{co}}{{\rm{s}}^{\rm{b}}}\left( {\rm{\theta }} \right)} \right)\]
\[{\rm{b}} =  {\rm{0}}{\rm{.45}}\;{\rm{ln}}\;{\left( {{{\rm{\rho }}^{\rm{w}}}} \right)^{}}\]
\[{\rm{\rho }}_{}^{\rm{w}} = \frac{{{\rm{\rho }}_{{\rm{bt}}}^{}\left( {{\rm{\theta }} = {\rm{0}}} \right)}}{{{\rm{1}}  {\rm{\tau }}_{{\rm{bb}}}^{}\left( {{\rm{\theta }} = {\rm{0}}} \right)}}\]
\[{\rho_{{\rm{bt}}}}\left( {{\rm{\theta }} = {\rm{90}}^\circ } \right) = {\rho_{{\rm{bt}}}}\left( {{\rm{\theta }} = {\rm{0}}} \right) + \left( {1  {\rho_{{\rm{bt}}}}\left( {{\rm{\theta }} = {\rm{0}}} \right)} \right)\;\left( {{\rm{0}}{\rm{.35}}{\rho ^{\rm{w}}}} \right)\]
The diffusediffuse material properties are obtained by Rhomberg numerical integration. The longwave properties of insect screen are given by expressions similar those formulated for drapery fabrics, and roller blinds (Kotey et al. 2008) is given by:
\[{\varepsilon ^{screen}} = {\varepsilon ^{\rm{w}}} \cdot \left( {{\rm{1  }}{{\rm{A}}_{\rm{o}}}} \right)\]
\[{\tau ^{{\rm{screen}}}} = {\rm{1  }}{\tau ^{\rm{w}}} \cdot \left( {{\rm{1  }}{{\rm{A}}_{\rm{o}}}} \right)\]
The apparent wire material emissivity can be taken as e^{w} = 0.93 for common insect screens  screens with dark, rough, nonmetallic wire surfaces. The corresponding infrared wire material transmittance is t^{w} = 0.98. A lower infrared emissivity can be used for screens constructed with shiny metallic wire. For example, to model stainlesssteel wire mesh use e^{w} = 0.32 and t^{w} = 0.81.
Integration with Heat Balance Method[LINK]
The solar model calculations are performed surface by surface following the existing structure of EnergyPlus’s heat balance algorithm. The optical and thermal calculations are performed for each surface at each time step. The thermal model needs to be updated at each time step during the surface heat balance iteration. Thus, window thermal model is invoked during the surface heat balance but only once. This is consistent with the current EnergyPlus’s window model inside surface heat balance iteration scheme. The equivalent layer window thermal model also has internal iteration scheme. The thermal model routine also requires the solar and shortwave radiation flux absorbed by each layer of a fenestration at each time step. The calculation of the absorbed radiation flux is performed using the existing scheme except that the optical properties are calculated by the equivalent layer optical model. The equivalent layer thermal model returns the temperatures and the fluxes at each layer.
Equivalent Layer Window Solar Model[LINK]
The equivalent layer window solar model calculates the transmittance of the window and absorptance of each layer. Separate optical properties calculations are performed for exterior beam and diffuse radiations, and interior diffuse radiation (e.g., lighting and interreflected solar gain). The fractions for beam radiations depend on solar incident angle and hence updated for each time step. Diffuse radiation fractions are the same for all time steps unless the shade characteristics are altered; for instance, when the venetian blind slat angle is controlled. Otherwise, the diffuse fractions are constant for a given configuration and time step. The hub of the optical properties calculation routine for Equivalent Layer window model is the routine “ASHWAT_Solar”. For beam radiation, at each time step, first the normal incidence optical properties are modified for current incident angle. Using the offnormal properties, then the ASHWAT_Solar routine sets up the coefficient matrix based on the net radiation concept to determine the effective absorptance of each layer and the window transmittance. The total intensity of shortwave radiation absorbed at each layer is determined by multiplying the incident solar radiation, reflected internal solar radiations, and internal shortwave source components by the appropriate absorbed fractions and summing them.
Equivalent Layer Window Thermal Model[LINK]
The equivalent layer thermal model is calculated only once for each time step. But the thermal model has internal iterative solutions scheme. During each time step, the procedure initially assumes that room air and means radiant temperatures are known. The surface heat balance loops through all zone surfaces while invoking the Equivalent Layer thermal model only once for each window surface at the first iteration. The window surfaces temperatures from the first iteration are used to complete the heat balances for the indoor (and implicitly the outdoor) face of each surface iteratively. Once indoor surface temperatures are calculated using the surface heat balance, the zone air temperature can be updated and the loads are predicted. Shaded fenestration in general do not have single inside temperature by virture their longwave radiation transmittance. The equivalent layer window model accounts for this using effective emissivity of the composite layers derived for each fenestration (ASHRAE 1311RP) as shown below:
\[{\varepsilon ^{eff}} = \sum\limits_{j = 0}^{nl} {\left[ {{\varepsilon_j} \cdot \prod\limits_{k = j{ + 1}}^{nl + {1}} {{\tau_k}} } \right]}\]
where
e^{eff} = composite indoor (roomside) longwave emissivity, ()
e_{j} = effective emissivity of layer j, ()
t_{k} = effective infrared transmittance of layer k (t_{nl+1} = 1)
nl = number of layers in fenestration system (glazing and shade). Layers are numbered outside to inside (layer 1 is outermost, layer nl is innermost).
Each equivalent layer window surafce yields net longwave radiant flux exchanged with the zone surfaces. Net longwave radiation exchange from the window to the zone is recast for a composite surface temperature calculation as follows:
\[{T^{eff}} = \sqrt[4]{{\frac{{{Q_{lw}}}}{{\sigma {\varepsilon ^{eff}}}}}} + {T_0}\]
where
T^{eff} = inside surface temperature of the composite layer, °C (°F)
Q_{lw} = infrared radiant gain from zone, W/m2 (Btu/hft2)
s = StefanBoltzmann constant, W/m2K4 (Btu/hft2R4)
T_{0} = Temperature of absolute zero, 273.15 °C (459.67 °F)
Recalculating effective inside face window temperature may result in extra convection flux. The “extra” (Other) convective flux is computed; this is the gain in excess of that resulting from the standard surface heat balance model. The net other convection term QX calculated below is added to the zone air heat balance (ASHRAE 1311RP).
\[QX = {Q_{conv}}  hc \cdot ({T^{eff}}  {T_a})\]
Where,
h_{c} = inside convective coefficient of the fenestration, W/m^{2}K (Btu/hft^{2}F)
Q_{conv} = total convective heat flux to zone from equivalent layer window thermal model, W/m^{2} (Btu/hft^{2}); includes openchannel gains and impact of inside surface convective coefficient
T_{a} = zone air temperature, °C (°F)
References:[LINK]
Edwards, D.K. 1977. Solar absorption by each element in an absorbercoverglass array, Technical Note, Solar Energy, Vol. 19, pp. 401402.
Parmelee, G. V., and W. W. Aubele. 1952. The shading of sunlit glass: an analysis of the effect of uniformly spaced flat opaque slats, ASHVE Transactions, Vol. 58, pp. 377398.
Farber, Erich A.; William A. Smith, C.W. Pennington, John C. Reed. 1963. Theoretical analysis of solar heat gain through insulating glass with inside shading. ASHRAE Transactions, Vol. 69, pp.393405.
Rheault, S., and E. Bilgen. 1989. Heat transfer analysis in an automated venetian blind system, Journal of Solar Energy, Vol. 111 (Feb.), pp. 8995.
Pfrommer, P., K. J. Lomas, and C. Kupke. 1996. “Solar Radiation Transport through SlatType Blinds: a New Model and its Application for Thermal Simulation of Buildings,” Solar Energy, Vol. 57, No. 2, pp. 7791.
Rosenfeld, J.L.J., W. J. Platzer, H. Van Dijk, and A. Maccari. 2000. “Modelling the Optical and Thermal Properties of Complex Glazing: Overview of Recent Developments”, Solar Energy, Vol. 69 Supplement, No. 16, pp.113.
Yahoda, D. S. and J. L. Wright. 2004. “Methods for Calculating the Effective Longwave Radiative Properties of a Venetian Blind Layer,” ASHRAE Transactions, Vol. 110, Pt. 1., pp. 463473.
Yahoda, D. S. and J. L. Wright. 2005. “Methods for Calculating the Effective SolarOptical Properties of a Venetian Blind Layer,” ASHRAE Transactions, Vol. 111, Pt. 1, pp. 572586.
Yahoda, D. S. and J. L. Wright. 2004. “Heat Transfer Analysis of a BetweenPanes Venetian Blind Using Effective Longwave Radiative Properties,” ASHRAE Transactions, Vol. 110, Pt. 1., pp. 455462.
Huang, N.Y.T., J. L. Wright, M. R. Collins. 2006. “Thermal Resistance of a Window with an Enclosed Venetian Blind: Guarded Heater Plate Measurements,” ASHRAE Transactions, Vol. 112, Pt. 2. pp. 1321.
Wright, J. L. 2008. “Calculating CentreGlass Performance Indices of Glazing Systems with Shading Devices,” ASHRAE Transactions, Vol. 114, Pt. 2.
Wright, J. L., N. Y. T. Huang, and M. R. Collins. 2008. “Thermal Resistance of a Window with an Enclosed Venetian Blind: A Simplified Model,” ASHRAE Transactions, Vol. 114, Pt. 1.
Kotey, N. A., J. L. Wright, and M. R. Collins. 2008. “Determining Longwave Radiative Properties of Flat Shading Materials,” 33rd Annual SESCI / 3rd CSBC Conference Proceedings, Fredericton, NB.
Kotey, N.A., Wright, J.L., M. R. Collins. 2009a. “Determination of AngleDependent Solar Optical Properties of Roller Blind Materials,” drafted for submission to ASHRAE Transactions, Vol. 115, Pt. 1.
Kotey, N.A., Wright, J.L., M. R. Collins. 2009b. “Determination of AngleDependent Solar Optical Properties of Drapery Fabrics,” in review, ASHRAE Transactions, Vol. 115, Pt. 2.
Wright, John L., Charles S. Barnaby, Michael R. Collins, and Nathan A. Kotey. Improving Cooling Load Calculations for Fenestration with Shading Devices . ASHRAE 1311RP, Final Report, February 11, 2009.
Documentation content copyright © 19962015 The Board of Trustees of the University of Illinois and the Regents of the University of California through the Ernest Orlando Lawrence Berkeley National Laboratory. All rights reserved. EnergyPlus is a trademark of the US Department of Energy.
This documentation is made available under the EnergyPlus Open Source License v1.0.