# Daylight Factor Calculation[LINK]

Table: Variables in Daylighting Calculations

Mathematical variable|Description|Units||FORTRAN variable ———————|———–|—–||—————- E_{h,sky}|Exterior horizontal illuminance due to light from the sky|lux|GILSK E_{h,sun}|Exterior horizontal illuminance due to light from the sun|lux|GILSU d_{sky}, d_{sun}|Interior illuminance factor due to sky, sun related light|-|DFACSK, DFACSU| w_{sky}, w_{sun}|Window luminance factor due to sky, sun related light|cd/lm|SFACSK, SFACSU b_{sky}, b_{sun}|Window background luminance factor due to sky, sun related light|cd/lm|BFACSK, BFACSU N|Number of exterior windows in a zone|-|NWD θ_{sky}, φ_{sky}|Azimuth and altitude angles of a point in the sky|radians|THSKY, PHSKY ψ_{cs}|Clear sky luminance distribution|cd/m^{2}|- ψ_{ts}|Clear turbid sky luminance distribution|cd/m^{2}|- ψ_{is}|Intermediate sky luminance distribution|cd/m^{2}|- ψ_{os}|Overcast sky luminance distribution|cd/m^{2}|- φ_{sun}|Altitude angle of the sun|radians or degrees|PHSUN γ|Angle between point in the sky and the sun; or angle between vertical and ray from reference point to window element|radians|G L_{z}|Sky zenith luminance|cd/m^{2}|ZENL m|Optical air mass of the atmosphere|m|AM h|Building altitude|m|Elevation E_{h,k}|Exterior horizontal illuminance for sky type k|lux|- N_{θ}, N_{φ}|Number of azimuth, altitude steps for sky integration||NTH, NPH |Vector from zone origin to reference point|m|RREF |Vector from zone origin to window element|m|RWIN dΩ|Solid angle subtended by window element|steradians|DOMEGA L_{w}|Luminance of a window element as seen from reference point|cd/m^{2}|WLUMSK, WLUMSU L_{w,shade}|Luminance of window element with shade in place|cd/m^{2}|WLUMSK, WLUMSU dE_{h}|Horizontal illuminance at reference point from window element|lux|- dx, dy|Size of window element|m|DWX, DWY D|Distance from reference point to window element|m|DIS B|Angle between window element’s outward normal and ray from reference point to window element|radians|- |Unit vector from reference point to window element|-|RAY |Unit vector normal to window element, pointing away from zone|-|WNORM |Unit vector along window y-axis|-|W21 |Unit vector along window x-axis|-|W23 τ_{vis}|Glass visible transmittance|-|TVISB L|Luminance of sky or obstruction |cd/m^{2}|ELUM, - Φ_{FW}|Downgoing luminous flux from a window|lm|FLFW– Φ_{CW}|Upgoing luminous flux from a window|lm|FLCW– F_{1}|First-reflected flux|lm|- ρ_{FW}|Area-weighted reflectance of floor and upper part of walls|-|SurfaceWindow%RhoFloorWall ρ_{CW}|Area-weighted reflectance of ceiling and upper part of walls|-|SurfaceWindow%RhoCeilingWall E_{r}|Average internally-reflected illuminance|lux|EINTSK, EINTSU A|Total inside surface area of a zone|m^{2}|ATOT ρ|Area-weighted average reflectance of zone interior surfaces|-|ZoneDaylight%AveVisDiffReflect θ, φ|Azimuth and altitude angle of a sky or ground element|radians|TH, PH L(θ,φ)|Luminance of sky or ground element at (θ,φ)|cd/m^{2}|HitPointLum– A_{w}|Area of glazed part of window|m^{2}|Surface%Area β|Angle of incidence, at center of window, of light from a sky or ground element|radians|- T(β)|Glazing visible transmittance at incidence angle β||TVISBR dΦ_{inc}|Luminous flux incident on window from sky or ground element|lm|- dΦ|Luminous flux from sky or ground element transmitted through window|lm|- dΦ_{FW}, dΦ_{CW}|Luminous flux from sky or ground element transmitted through window and goind downwar, upward|lm|- θ_{min}, θ_{max}|Azimuth angle integration limits|radians|THMIN, THMAX φ_{w}|Window normal altitude angle|radians|- Φ_{sh}, Φ_{unsh},|Transmitted flux through window with shade, without shade|lm|- Φ_{CW,sh}, Φ_{FW,sh}|Upgoing and downgoing portions of transmitted flux through window with shade|lm|- Φ_{CW,unsh}, Φ_{FW,unsh}|Upgoing and downgoing portions of transmitted flux through window without shade|lm|- f|Fraction of hemisphere seen by the inside of window that lies above the window midplane|-|SurfaceWindow%FractionUpgoing Φ_{inc}|Flux incident on glazing from direct sun|lm|- f_{sunlit}|Fraction of glazing that is sunlit |-|SunLitFrac Φ|Transmitted flux from direct sun|-|- L_{sh}|Luminance of window with shade|cd/m^{2}|- L_{b}|Window background luminance|cd/m^{2}|BLUM G|Discomfort glare constant|-|GTOT G_{i}|Discomfort glare constant from window i|-|- ω|Solid angle subtended by window with respect to reference point|steradians|SolidAngAtRefPt Ω|Solid angle subtended by window with respect to reference point, modified to take direction of occupant view into account|steradians|SolidAngAtRefPtWtd N_{x,} N_{y}|Number of elements in x and y direction that window is divided into for glare calculation|-|NWX, NWY p(x_{R}, y_{R})|Position factor for horizontal and vertical displacement ratios x_{R} and y_{R}|-|DayltgGlarePositionFactor p_{H}|Hopkinson position factor|-|DayltgGlarePositionFactor L_{b}|Window background luminance|cd/m^{2}|BLUM E_{b}|Illuminance on window background|lm|- E_{r}|Total internally-reflected component of daylight illuminance|lm|- E_{s}|Illuminance setpoint|lm|IllumSetPoint G_{I}|Glare index|-|GLINDX

There are three types of daylight factors: interior illuminance factors, window luminance factors, and window background luminance factors. To calculate these factors the following steps are carried out for each hourly sun position on the sun paths for the design days and for representative days during the simulation run period:

Calculate exterior horizontal daylight illuminance from sky and sun for standard (CIE) clear and overcast skies.

Calculate interior illuminance, window luminance and window background luminance for each window/reference-point combination, for bare and for shaded window conditions (if a shading device has been specified), for overcast sky and for standard clear sky.

Divide by exterior horizontal illuminance to obtain daylight factors.

## Interior Illuminance Components[LINK]

To calculate daylight factors, daylight incident on a window is separated into two components: (1) light that originates from the *sky* and reaches the window directly or by reflection from exterior surfaces; and (2) light that originates from the *sun* and reaches the window directly or by reflection from exterior surfaces. Light from the window reaches the workplane directly or via reflection from the interior surfaces of the room.

For fixed sun position, sky condition (clear or overcast) and room geometry, the sky-related interior daylight will be proportional to the exterior horizontal illuminance, *E*_{h,sky}, due to light from the sky. Similarly, the sun-related interior daylight will be proportional to the exterior horizontal solar illuminance, *E*_{h,sun}.

## Daylight Factors[LINK]

The following daylight factors are calculated:

For a daylit zone with *N* windows these six daylight factors are calculated for each of the following combinations of reference point, window, sky-condition/sun-position and shading device:

## Sky Luminance Distributions[LINK]

The luminance distribution of the sky is represented as a superposition of four standard CIE skies using the approach described in (Perez et al. 1990). The standard skies are as follows.

The clear sky luminance distribution has the form (Kittler, 1965; CIE, 1973)

Here, *L*_{z} is the zenith luminance (i.e., the luminance of the sky at a point directly overhead). In the calculation of daylight factors, which are ratios of interior and exterior illumination quantities that are both proportional to *L*_{z}, the zenith luminance cancels out. For this reason we will use *L*_{z} = 1.0 for all sky luminance distributions.

The various angles, which are defined in the building coordinate system, are shown in Figure 48. The angle, *γ*, between sun and sky element is given by

The general characteristics of the clear-sky luminance distribution are a large peak near the sun; a minimum at a point on the other side of the zenith from the sun, in the vertical plane containing the sun; and an increase in luminance as the horizon is approached.

### Clear Turbid Sky[LINK]

The clear turbid sky luminance distribution has the form [Matsuura, 1987]

The intermediate sky luminance distribution has the form [Matsuura, 1987]

where

The overcast sky luminance distribution has the form [Moon & Spencer, 1942]

Unlike the clear sky case, the overcast sky distribution does not depend on the solar azimuth or the sky azimuth. Note that at fixed solar altitude the zenith ( ) is three times brighter than the horizon ( ).

## Direct Normal Solar Illuminance[LINK]

For purposes of calculating daylight factors associated with beam solar illuminance, the direct normal solar illuminance is taken to be 1.0 W/m^{2}. The actual direct normal solar illuminance, determined from direct normal solar irradiance from the weather file and empirically-determined luminious efficacy, is used in the time-step calculation.

## Exterior Horizontal Illuminance[LINK]

The illuminance on an unobstructed horizontal plane due to diffuse radiation from the sky is calculated for each of the four sky types by integrating over the appropriate sky luminance distribution:

where *k* = *cs*, *ts*, *is* or *os*. The integral is evaluated as a double summation:

where

and were found to give a accuracy in the calculation of .

## Direct Component of Interior Daylight Illuminance[LINK]

The direct daylight illuminance at a reference point from a particular window is determined by dividing the window into an x-y grid and finding the flux reaching the reference point from each grid element. The geometry involved is shown in Figure 57. The horizontal illuminance at the reference point, , due to a window element is

where *L*_{w} is the luminance of the window element as seen from the reference point.

The subtended solid angle is approximated by

where

*CosB* is found from

where

Equation becomes exact as and is accurate to better than about 1% for .

The net illuminance from the window is obtained by summing the contributions from all the window elements:

In performing the summation, window elements that lie below the workplane ( ) are omitted since light from these elements cannot reach the workplane directly.

### Unshaded Window[LINK]

For the unshaded window case, the luminance of the window element is found by projecting the ray from reference point to window element and determining whether it intersects the sky or an exterior obstruction such as an overhang. If *L* is the corresponding luminance of the sky or obstruction, the window luminance is

where is the visible transmittance of the glass for incidence angle *B*.

Exterior obstructions are generally opaque (like fins, overhangs, neighboring buildings, and the building’s own wall and roof surfaces) but can be transmitting (like a tree or translucent awning). Exterior obstructions are assumed to be non-reflecting. If *L*_{sky} is the sky luminance and _{obs} is the transmittance of the obstruction (assumed independent of incidence angle), then *L* = *L*_{skyobs}. Interior obstructions are assumed to be opaque (_{obs} = 0).

### Shaded Window[LINK]

For the window-plus-shade case the shade is assumed to be a perfect diffuser, i.e., the luminance of the shade is independent of angle of emission of light, position on shade, and angle of incidence of solar radiation falling on the shade. Closely-woven drapery fabric and translucent roller shades are closer to being perfect diffusers than Venetian blinds or other slatted devices, which usually have non-uniform luminance characteristics.

The calculation of the window luminance with the shade in place, *L*_{w,sh}, is described in [Winkelmann, 1983]. The illuminance contribution at the reference point from a shaded window element is then given by Eq. (152) with .

## Internally-Reflected Component of Interior Daylight Illuminance[LINK]

Daylight reaching a reference point after reflection from interior surfaces is calculated using the *split-flux* method [Hopkinson et al., 1954], [Lynes, 1968]. In this method the daylight transmitted by the window is split into two parts—a downward-going flux, (lumens), which falls on the floor and portions of the walls below the imaginary horizontal plane passing through the center of the window (*window midplane*), and an upward-going flux, , that strikes the ceiling and portions of the walls above the window midplane. A fraction of these fluxes is absorbed by the room surfaces. The remainder, the first-reflected flux, *F*_{1}, is approximated by

where *ρ*_{FW} is the area-weighted average reflectance of the floor and those parts of the walls below the window midplane, and *ρ*_{CW} is the area-weighted average reflectance of the ceiling and those parts of the walls above the window midplane.

To find the final average internally-reflected illuminance, *E*_{r}, on the room surfaces (which in this method is uniform throughout the room) a flux balance is used. The total reflected flux absorbed by the room surfaces (or lost through the windows) is *AE*_{r}(1-ρ), where *A* is the total inside surface area of the floor, walls, ceiling and windows in the room, and *ρ* is the area-weighted average reflectance of the room surfaces, including windows. From conservation of energy

or

This procedure assumes that the room behaves like an integrating sphere with perfectly diffusing interior surfaces and with no internal obstructions. It therefore works best for rooms that are close to cubical in shape, have matte surfaces (which is usually the case), and have no internal partitions. Deviations from these conditions, such as would be the case for rooms whose depth measured from the window-wall is more than three times greater than ceiling height, can lead to substantial inaccuracies in the split-flux calculation.

## Transmitted Flux from Sky and Ground[LINK]

The luminous flux incident on the center of the window from a luminous element of sky or ground at angular position , of luminance , and subtending a solid angle is

The transmitted flux is

where *T(β)* is the window transmittance for light at incidence angle *β*. This transmittance depends on whether or not the window has a shade.

For an unshaded window the total downgoing transmitted flux is obtained by integrating over the part of the exterior hemisphere seen by the window that lies above the window midplane. This gives

The upgoing flux is obtained similarly by integrating over the part of the exterior hemisphere that lies below the window midplane:

where is the angle the window outward normal makes with the horizontal plane.

For a window with a diffusing shade the total transmitted flux is

The downgoing and upgoing portions of this flux are

where *f*, the fraction of the hemisphere seen by the inside of the window that lies above the window midplane, is given by

For a vertical window ( ) the up- and down-going transmitted fluxes are equal:

.

For a horizontal skylight ( ):

The limits of integration of *θ* in Equations (153), (154) and (155) depend on . From [Figure 12 - Winkelmann, 1983] we have

which gives

Thus

## Transmitted Flux from Direct Sun[LINK]

The flux incident on the window from direct sun is

The transmitted flux is

where T is the net transmittance of the window glazing (plus shade, if present).

For an unshaded window all of the transmitted flux is downward since the sun always lies above the window midplane. Therefore

For a window with a diffusing shade

## Luminance of Shaded Window[LINK]

The luminance of a shaded window is determined at the same time that the transmitted flux is calculated. It is given by

## Daylight Discomfort Glare[LINK]

The discomfort glare at a reference point due to luminance contrast between a window and the interior surfaces surrounding the window is given by [Hopkinson, 1970] and [Hopkinson, 1972]:

where

*G* = discomfort glare constant

*L*_{w} = average luminance of the window as seen from the reference point

*Ω* = solid angle subtended by window, modified to take direction of occupant view into account

*L*_{b} = luminance of the background area surrounding the window

By dividing the window into *N*_{x} by *N*_{y} rectangular elements, as is done for calculating the direct component of interior illuminance, we have

where *L*_{w}(i,j) is the luminance of element *(i,j)* as seen from the reference point.

Similarly,

where *dω(i,j)* is the solid angle subtended by element *(i,j)* with respect to the reference point.

The modified solid angle is

where p is a “position factor” [Petherbridge & Longmore, 1954] that accounts for the decrease in visual excitation as the luminous element moves away from the line of sight. This factor depends on the horizontal and vertical displacement ratios, *x*_{R} and *y*_{R} (Figure 58), given by

where

The factor *p* can be obtained from graphs given in [Petherbridge & Longmore, 1954] or it can be calculated from tabulated values of *p*_{H}, the Hopkinson position factor [Hopkinson, 1966], since . The values resulting from the latter approach are given in Table 23. Interpolation of this table is used in EnergyPlus to evaluate *p* at intermediate values of *x*_{R} and *y*_{R}.

Position factor for glare calculation
**x**_{R}: Horizontal Displacement Factor |

**0** |

**y**_{R}: |

**0.5** |

**1.0** |

**1.5** |

**2.0** |

**>2.0** |

The background luminance is

where *ρ*_{b}~~is approximated by the average interior surface reflectance of the entire room and

where *E*_{r} is the total internally-reflected component of daylight illuminance produced by all the windows in the room and *E*_{s} is the illuminance setpoint at the reference point at which glare is being calculated. A precise calculation of *E*_{b} is not required since the glare index (see next section) is logarithmic. A factor of two variation in *E*_{b} generally produces a change of only 0.5 to 1.0 in the glare index.

The net daylight glare at a reference point due to all of the windows in a room is expressed in terms of a *glare index* given by

where *G*_{i} is the glare constant at the reference point due to the *i*^{th}^^window

## Daylight Factor Calculation[LINK]

Table: Variables in Daylighting Calculations

Mathematical variable|Description|Units||FORTRAN variable ———————|———–|—–||—————- E

_{h,sky}|Exterior horizontal illuminance due to light from the sky|lux|GILSK E_{h,sun}|Exterior horizontal illuminance due to light from the sun|lux|GILSU d_{sky}, d_{sun}|Interior illuminance factor due to sky, sun related light|-|DFACSK, DFACSU| w_{sky}, w_{sun}|Window luminance factor due to sky, sun related light|cd/lm|SFACSK, SFACSU b_{sky}, b_{sun}|Window background luminance factor due to sky, sun related light|cd/lm|BFACSK, BFACSU N|Number of exterior windows in a zone|-|NWD θ_{sky}, φ_{sky}|Azimuth and altitude angles of a point in the sky|radians|THSKY, PHSKY ψ_{cs}|Clear sky luminance distribution|cd/m^{2}|- ψ_{ts}|Clear turbid sky luminance distribution|cd/m^{2}|- ψ_{is}|Intermediate sky luminance distribution|cd/m^{2}|- ψ_{os}|Overcast sky luminance distribution|cd/m^{2}|- φ_{sun}|Altitude angle of the sun|radians or degrees|PHSUN γ|Angle between point in the sky and the sun; or angle between vertical and ray from reference point to window element|radians|G L_{z}|Sky zenith luminance|cd/m^{2}|ZENL m|Optical air mass of the atmosphere|m|AM h|Building altitude|m|Elevation E_{h,k}|Exterior horizontal illuminance for sky type k|lux|- N_{θ}, N_{φ}|Number of azimuth, altitude steps for sky integration||NTH, NPH |Vector from zone origin to reference point|m|RREF |Vector from zone origin to window element|m|RWIN dΩ|Solid angle subtended by window element|steradians|DOMEGA L_{w}|Luminance of a window element as seen from reference point|cd/m^{2}|WLUMSK, WLUMSU L_{w,shade}|Luminance of window element with shade in place|cd/m^{2}|WLUMSK, WLUMSU dE_{h}|Horizontal illuminance at reference point from window element|lux|- dx, dy|Size of window element|m|DWX, DWY D|Distance from reference point to window element|m|DIS B|Angle between window element’s outward normal and ray from reference point to window element|radians|- |Unit vector from reference point to window element|-|RAY |Unit vector normal to window element, pointing away from zone|-|WNORM |Unit vector along window y-axis|-|W21 |Unit vector along window x-axis|-|W23 τ_{vis}|Glass visible transmittance|-|TVISB L|Luminance of sky or obstruction |cd/m^{2}|ELUM, - Φ_{FW}|Downgoing luminous flux from a window|lm|FLFW– Φ_{CW}|Upgoing luminous flux from a window|lm|FLCW– F_{1}|First-reflected flux|lm|- ρ_{FW}|Area-weighted reflectance of floor and upper part of walls|-|SurfaceWindow%RhoFloorWall ρ_{CW}|Area-weighted reflectance of ceiling and upper part of walls|-|SurfaceWindow%RhoCeilingWall E_{r}|Average internally-reflected illuminance|lux|EINTSK, EINTSU A|Total inside surface area of a zone|m^{2}|ATOT ρ|Area-weighted average reflectance of zone interior surfaces|-|ZoneDaylight%AveVisDiffReflect θ, φ|Azimuth and altitude angle of a sky or ground element|radians|TH, PH L(θ,φ)|Luminance of sky or ground element at (θ,φ)|cd/m^{2}|HitPointLum– A_{w}|Area of glazed part of window|m^{2}|Surface%Area β|Angle of incidence, at center of window, of light from a sky or ground element|radians|- T(β)|Glazing visible transmittance at incidence angle β||TVISBR dΦ_{inc}|Luminous flux incident on window from sky or ground element|lm|- dΦ|Luminous flux from sky or ground element transmitted through window|lm|- dΦ_{FW}, dΦ_{CW}|Luminous flux from sky or ground element transmitted through window and goind downwar, upward|lm|- θ_{min}, θ_{max}|Azimuth angle integration limits|radians|THMIN, THMAX φ_{w}|Window normal altitude angle|radians|- Φ_{sh}, Φ_{unsh},|Transmitted flux through window with shade, without shade|lm|- Φ_{CW,sh}, Φ_{FW,sh}|Upgoing and downgoing portions of transmitted flux through window with shade|lm|- Φ_{CW,unsh}, Φ_{FW,unsh}|Upgoing and downgoing portions of transmitted flux through window without shade|lm|- f|Fraction of hemisphere seen by the inside of window that lies above the window midplane|-|SurfaceWindow%FractionUpgoing Φ_{inc}|Flux incident on glazing from direct sun|lm|- f_{sunlit}|Fraction of glazing that is sunlit |-|SunLitFrac Φ|Transmitted flux from direct sun|-|- L_{sh}|Luminance of window with shade|cd/m^{2}|- L_{b}|Window background luminance|cd/m^{2}|BLUM G|Discomfort glare constant|-|GTOT G_{i}|Discomfort glare constant from window i|-|- ω|Solid angle subtended by window with respect to reference point|steradians|SolidAngAtRefPt Ω|Solid angle subtended by window with respect to reference point, modified to take direction of occupant view into account|steradians|SolidAngAtRefPtWtd N_{x,}N_{y}|Number of elements in x and y direction that window is divided into for glare calculation|-|NWX, NWY p(x_{R}, y_{R})|Position factor for horizontal and vertical displacement ratios x_{R}and y_{R}|-|DayltgGlarePositionFactor p_{H}|Hopkinson position factor|-|DayltgGlarePositionFactor L_{b}|Window background luminance|cd/m^{2}|BLUM E_{b}|Illuminance on window background|lm|- E_{r}|Total internally-reflected component of daylight illuminance|lm|- E_{s}|Illuminance setpoint|lm|IllumSetPoint G_{I}|Glare index|-|GLINDX## Overview[LINK]

There are three types of daylight factors: interior illuminance factors, window luminance factors, and window background luminance factors. To calculate these factors the following steps are carried out for each hourly sun position on the sun paths for the design days and for representative days during the simulation run period:

Calculate exterior horizontal daylight illuminance from sky and sun for standard (CIE) clear and overcast skies.

Calculate interior illuminance, window luminance and window background luminance for each window/reference-point combination, for bare and for shaded window conditions (if a shading device has been specified), for overcast sky and for standard clear sky.

Divide by exterior horizontal illuminance to obtain daylight factors.

## Interior Illuminance Components[LINK]

To calculate daylight factors, daylight incident on a window is separated into two components: (1) light that originates from the

skyand reaches the window directly or by reflection from exterior surfaces; and (2) light that originates from thesunand reaches the window directly or by reflection from exterior surfaces. Light from the window reaches the workplane directly or via reflection from the interior surfaces of the room.For fixed sun position, sky condition (clear or overcast) and room geometry, the sky-related interior daylight will be proportional to the exterior horizontal illuminance,

E, due to light from the sky. Similarly, the sun-related interior daylight will be proportional to the exterior horizontal solar illuminance,_{h,sky}E._{h,sun}## Daylight Factors[LINK]

The following daylight factors are calculated:

For a daylit zone with

Nwindows these six daylight factors are calculated for each of the following combinations of reference point, window, sky-condition/sun-position and shading device:## Sky Luminance Distributions[LINK]

The luminance distribution of the sky is represented as a superposition of four standard CIE skies using the approach described in (Perez et al. 1990). The standard skies are as follows.

## Clear Sky[LINK]

The clear sky luminance distribution has the form (Kittler, 1965; CIE, 1973)

Here,

Lis the zenith luminance (i.e., the luminance of the sky at a point directly overhead). In the calculation of daylight factors, which are ratios of interior and exterior illumination quantities that are both proportional to_{z}L, the zenith luminance cancels out. For this reason we will use_{z}L= 1.0 for all sky luminance distributions._{z}The various angles, which are defined in the building coordinate system, are shown in Figure 48. The angle,

γ, between sun and sky element is given byThe general characteristics of the clear-sky luminance distribution are a large peak near the sun; a minimum at a point on the other side of the zenith from the sun, in the vertical plane containing the sun; and an increase in luminance as the horizon is approached.

## Clear Turbid Sky[LINK]

The clear turbid sky luminance distribution has the form [Matsuura, 1987]

## Intermediate Sky[LINK]

The intermediate sky luminance distribution has the form [Matsuura, 1987]

where

Angles appearing in the expression for the clear-sky luminance distribution.

## Overcast Sky[LINK]

The overcast sky luminance distribution has the form [Moon & Spencer, 1942]

Unlike the clear sky case, the overcast sky distribution does not depend on the solar azimuth or the sky azimuth. Note that at fixed solar altitude the zenith ( ) is three times brighter than the horizon ( ).

## Direct Normal Solar Illuminance[LINK]

For purposes of calculating daylight factors associated with beam solar illuminance, the direct normal solar illuminance is taken to be 1.0 W/m

^{2}. The actual direct normal solar illuminance, determined from direct normal solar irradiance from the weather file and empirically-determined luminious efficacy, is used in the time-step calculation.## Exterior Horizontal Illuminance[LINK]

The illuminance on an unobstructed horizontal plane due to diffuse radiation from the sky is calculated for each of the four sky types by integrating over the appropriate sky luminance distribution:

where

k=cs,ts,isoros. The integral is evaluated as a double summation:where

and were found to give a accuracy in the calculation of .

## Direct Component of Interior Daylight Illuminance[LINK]

The direct daylight illuminance at a reference point from a particular window is determined by dividing the window into an x-y grid and finding the flux reaching the reference point from each grid element. The geometry involved is shown in Figure 57. The horizontal illuminance at the reference point, , due to a window element is

where

Lis the luminance of the window element as seen from the reference point._{w}The subtended solid angle is approximated by

where

CosBis found fromwhere

Equation becomes exact as and is accurate to better than about 1% for .

The net illuminance from the window is obtained by summing the contributions from all the window elements:

In performing the summation, window elements that lie below the workplane ( ) are omitted since light from these elements cannot reach the workplane directly.

Geometry for calculation of direct component of daylight illuminance at a reference point. Vectors R

_{ref}, W_{1}, W_{2}, W_{3}and R_{win}are in the building coordinate system.## Unshaded Window[LINK]

For the unshaded window case, the luminance of the window element is found by projecting the ray from reference point to window element and determining whether it intersects the sky or an exterior obstruction such as an overhang. If

Lis the corresponding luminance of the sky or obstruction, the window luminance iswhere is the visible transmittance of the glass for incidence angle

B.Exterior obstructions are generally opaque (like fins, overhangs, neighboring buildings, and the building’s own wall and roof surfaces) but can be transmitting (like a tree or translucent awning). Exterior obstructions are assumed to be non-reflecting. If

Lis the sky luminance and_{sky}is the transmittance of the obstruction (assumed independent of incidence angle), then_{obs}L=L. Interior obstructions are assumed to be opaque (_{skyobs}= 0)._{obs}## Shaded Window[LINK]

For the window-plus-shade case the shade is assumed to be a perfect diffuser, i.e., the luminance of the shade is independent of angle of emission of light, position on shade, and angle of incidence of solar radiation falling on the shade. Closely-woven drapery fabric and translucent roller shades are closer to being perfect diffusers than Venetian blinds or other slatted devices, which usually have non-uniform luminance characteristics.

The calculation of the window luminance with the shade in place,

L, is described in [Winkelmann, 1983]. The illuminance contribution at the reference point from a shaded window element is then given by Eq. (152) with ._{w,sh}## Internally-Reflected Component of Interior Daylight Illuminance[LINK]

Daylight reaching a reference point after reflection from interior surfaces is calculated using the

split-fluxmethod [Hopkinson et al., 1954], [Lynes, 1968]. In this method the daylight transmitted by the window is split into two parts—a downward-going flux, (lumens), which falls on the floor and portions of the walls below the imaginary horizontal plane passing through the center of the window (window midplane), and an upward-going flux, , that strikes the ceiling and portions of the walls above the window midplane. A fraction of these fluxes is absorbed by the room surfaces. The remainder, the first-reflected flux,F, is approximated by_{1}where

ρis the area-weighted average reflectance of the floor and those parts of the walls below the window midplane, and_{FW}ρis the area-weighted average reflectance of the ceiling and those parts of the walls above the window midplane._{CW}To find the final average internally-reflected illuminance,

E, on the room surfaces (which in this method is uniform throughout the room) a flux balance is used. The total reflected flux absorbed by the room surfaces (or lost through the windows) is_{r}AE, where_{r}(1-ρ)Ais the total inside surface area of the floor, walls, ceiling and windows in the room, andρis the area-weighted average reflectance of the room surfaces, including windows. From conservation of energyor

This procedure assumes that the room behaves like an integrating sphere with perfectly diffusing interior surfaces and with no internal obstructions. It therefore works best for rooms that are close to cubical in shape, have matte surfaces (which is usually the case), and have no internal partitions. Deviations from these conditions, such as would be the case for rooms whose depth measured from the window-wall is more than three times greater than ceiling height, can lead to substantial inaccuracies in the split-flux calculation.

## Transmitted Flux from Sky and Ground[LINK]

The luminous flux incident on the center of the window from a luminous element of sky or ground at angular position , of luminance , and subtending a solid angle is

The transmitted flux is

where

T(β)is the window transmittance for light at incidence angleβ. This transmittance depends on whether or not the window has a shade.For an unshaded window the total downgoing transmitted flux is obtained by integrating over the part of the exterior hemisphere seen by the window that lies above the window midplane. This gives

The upgoing flux is obtained similarly by integrating over the part of the exterior hemisphere that lies below the window midplane:

where is the angle the window outward normal makes with the horizontal plane.

For a window with a diffusing shade the total transmitted flux is

The downgoing and upgoing portions of this flux are

where

f, the fraction of the hemisphere seen by the inside of the window that lies above the window midplane, is given byFor a vertical window ( ) the up- and down-going transmitted fluxes are equal:

.

For a horizontal skylight ( ):

The limits of integration of

θin Equations (153), (154) and (155) depend on . From [Figure 12 - Winkelmann, 1983] we havewhich gives

Thus

## Transmitted Flux from Direct Sun[LINK]

The flux incident on the window from direct sun is

The transmitted flux is

where T is the net transmittance of the window glazing (plus shade, if present).

For an unshaded window all of the transmitted flux is downward since the sun always lies above the window midplane. Therefore

For a window with a diffusing shade

## Luminance of Shaded Window[LINK]

The luminance of a shaded window is determined at the same time that the transmitted flux is calculated. It is given by

## Daylight Discomfort Glare[LINK]

The discomfort glare at a reference point due to luminance contrast between a window and the interior surfaces surrounding the window is given by [Hopkinson, 1970] and [Hopkinson, 1972]:

where

G= discomfort glare constantL= average luminance of the window as seen from the reference point_{w}Ω= solid angle subtended by window, modified to take direction of occupant view into accountL= luminance of the background area surrounding the window_{b}By dividing the window into

Nby_{x}Nrectangular elements, as is done for calculating the direct component of interior illuminance, we have_{y}where

Lis the luminance of element_{w}(i,j)(i,j)as seen from the reference point.Similarly,

where

dω(i,j)is the solid angle subtended by element(i,j)with respect to the reference point.The modified solid angle is

where p is a “position factor” [Petherbridge & Longmore, 1954] that accounts for the decrease in visual excitation as the luminous element moves away from the line of sight. This factor depends on the horizontal and vertical displacement ratios,

xand_{R}y(Figure 58), given by_{R}where

Geometry for calculation of displacement ratios used in the glare formula.

The factor

pcan be obtained from graphs given in [Petherbridge & Longmore, 1954] or it can be calculated from tabulated values ofp, the Hopkinson position factor [Hopkinson, 1966], since . The values resulting from the latter approach are given in Table 23. Interpolation of this table is used in EnergyPlus to evaluate_{H}pat intermediate values ofxand_{R}y._{R}x_{R}: Horizontal Displacement Factor0y_{R}:0.51.01.52.0>2.0The background luminance is

where

ρ~~is approximated by the average interior surface reflectance of the entire room and_{b}where

Eis the total internally-reflected component of daylight illuminance produced by all the windows in the room and_{r}Eis the illuminance setpoint at the reference point at which glare is being calculated. A precise calculation of_{s}Eis not required since the glare index (see next section) is logarithmic. A factor of two variation in_{b}Egenerally produces a change of only 0.5 to 1.0 in the glare index._{b}## Glare Index[LINK]

The net daylight glare at a reference point due to all of the windows in a room is expressed in terms of a

glare indexgiven bywhere

Gis the glare constant at the reference point due to the_{i}i^^window^{th}Documentation content copyright © 1996-2014 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.