Simulating the daylight performance of complex fenestration systems using bidirectional scattering distribution functions within radiance.
As energy efficiency and indoor environmental quality objectives become increasingly aggressive in this charged environment to mitigate climate change, it becomes more important that window performance be quantified accurately. A routine, efficient method is urgently needed to model complex fenestration systems such as sunlight-redirecting mirrored louvers, holographic optical elements, and even conventional roller shades.
To address this problem, a method was defined to quantify solar gains through windows with nonspecular shading devices, termed complex fenestration systems (CFS), that bridged the gap between expensive solar calorimeter measurements and first-principles analytical models [Klems 1994a]. The method requires bidirectional optical measurements of nonspecular layers such as shades or blinds, use of optical data for specular glass layers, and matrix layer calculations to combine these individual layer properties to produce overall window system transmittances and layer absorptances from which solar gains can be computed as a function of incident direction. This method is also applicable to the assessment of daylighting.
Bidirectional transmittance and reflectance or scattering function (BSDF) data for nonspecular devices is critical for routine implementation of this method and obtaining such data has been and continues to be a nontrivial task. A detailed review of the methods used to characterize nonspecular materials is given by Andersen and deBoer : starting with the first scanning goniophotometer built in 1988 and validation of the calculation method [Klems and others, 1997] to today's heterogeneous methods encompassing measurement and simulation techniques [Jonsson and others, 2009; Andersen and others, 2009]. Modifications to the Window software were made to enable combination of layer data into any arbitrarily configured system of layers [Mitchell and others, 2008] ; its database now contains BSDF data for some common shading systems. Other databases are also publicly available [DeBoer 2006], but differ in format.
In parallel, BSDF-enabled simulation tools have been and are currently under development. The DELight subroutine was implemented within EnergyPlus to enable daylight modeling of CFS using radiosity methods [Carroll and Hitchcock 2005]. Further modifications to EnergyPlus are underway to enable calculation of solar gains. Kaemph  adapted the prism2 function within Radiance to model CFS, deriving two principle emerging directions from measured BTDF data. DeBoer  implemented a method within Adeline which treats the CFS as a light-emitting surface.
This paper describes modifications made to Radiance [Ward 1998] which enable both visualizations and annual daylighting evaluations of CFS. A brief overview of the fundamental underlying method is given, then the methods of implementation within Radiance are described. Validation and limitations of the methods are discussed.
Klems showed that with a few simplifying assumptions, it is possible to derive the overall system transmittance matrix of a set of layers, including all the effects of multiple interreflections between layers, from the optical properties of individual layers [Klems 1994a, 1994b]. Two approximations were made:
1. Optical properties are spatially averaged over a suitablysized area so that spatially inhomogeneous systems, such as venetian blinds with repetitions of unit slats, are treated as homogeneous layers. Direct radiation through the blind would be a uniformly lit patch rather than a series of stripes, which was deemed sufficient for the purpose of solar gain calculation. Spatial averaging for two parallel layers with spatial inhomogeneities of comparable size will not be correctly treated: for example, two parallel, partially-closed blinds positioned such that the opening between slats on one blind is aligned with the closed slat of the second blind will not result in zero direct transmittance.
[FIGURE 1 OMITTED]
2. The original method described in 1994 assumed properties that were measured using a source that very roughly approximated the solar spectrum after transmittance through a layer of glass; that is, spectrally averaged and unpolarized. However, Window 6 implemented the Klems method and applies the matrix layer calculation method at each wavelength, if detailed spectral information are available for the individual layers, to produce spectralangular system properties which can then be used directly or appropriately averaged. Interference and polarization effects are at present not accounted for.
With these approximations, the solar optical properties of a system can be defined in terms of its bidirectional transmittance and reflectance distribution functions [Nicodemus 1965], which are in effect a set of hemispherical luminous coefficients defined by paired incident and outgoing angles. The incoming and outgoing hemispheres are subdivided into a grid of elements (referred to as a "basis") and the radiance is averaged over each solid angle element or "patch". To obtain total transmitted radiation for a given incident angle, [[theta].sub.1], [[phi].sub.1], the luminous coefficients are multiplied by the incident irradiance and summed for all patches of the hemispherical basis using a matrix calculation, as in (1),
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (1)
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (2) where, [theta] and [phi] (Fig. 1) define the boundaries of each patch of the basis and (2) shows integration of flux over the Klems full 145 x 145 basis, which is the default basis used in Window 6. (1)
For multilayer calculations, output flux can be used as input for the next layer, enabling one to create a coplanar window system from an arbitrary configuration of individual layers (for example, glass layer, between pane shade, glass). The method assumes that the layer planes are infinitely long compared to the depth of the system (1/day > 5), ignoring the edge effects of the window frame. Similarly, total reflected radiation for a multilayer system can be computed. The outgoing total radiation is then distributed to interior room surfaces to calculate solar gains.
When this method is applied to daylighting assessments, spatial averaging over a patch results is perhaps the most limiting aspect of this method. The spatial resolution of outgoing flux is as coarse as the coordinate system or basis and hence, one cannot resolve features smaller than that. Radiance provides work-arounds for some of these issues, which is discussed in the following sections.
3 CFS MODELING IN RADIANCE
Until recently, there was no way to effectively utilize general BSDF data in Radiance. Such data were difficult to obtain, and representation standards were lacking. With the development of Window 6, these barriers are lifting.
The backwards ray-tracing approach in Radiance needs to know a priori where to look for light sources and where to find peaks in the BSDF, so we know where to send our samples for optimal results. The standard approach for CFS in Radiance is to precompute the light passing through each window or skylight, treating it as a "white box" in the final rendering, otherwise known as a "secondary light source." The original mkillum program was written to convert a CFS, described as geometry and materials, into a secondary light source with a luminous intensity distribution corresponding to a particular exterior condition (that is, sun, sky, obstructions, reflections, etc.). In this mode, mkillum is computing all the light interactions of whatever CFS might be placed in the opening. During rendering, thus-converted windows and skylights will usually have a dual nature, applying the mkillum-computed light distribution for illuminating interior surfaces, while applying the original CFS geometry and materials for direct viewing and beam penetration. Computing light at a point on an interior surface then involves tracing "shadow rays" to all light sources, including secondary sources that turn into their original geometry when intersected. In the case of venetian blinds, this process provides a representative view out the window, since the actual blinds geometry is represented, plus striped shadows on the floor when sunlight penetrates.
While the standard mkillum technique just described works well most of the time, there are important cases where it fails. Using a backwards ray-tracing approach, it is nearly impossible to account for sunlight transmitting through a system with curved, specular elements. However, if we are given the BSDF as input, we can apply it directly to convert the exterior illumination arriving at the window or skylight into an interior luminous intensity distribution, as we describe in the following section. A photon-mapping extension (that is, forward ray-tracing) may also be used for such systems, although it fundamentally changes the way Radiance works [Shregle & Weinhold 2004]. Neither approach is practical for annual simulation, since a separate precomputation step is required for every solar condition. A daylight coefficient (DC) method is needed.
The Daysim version of Radiance was developed specifically for annual simulation using daylight coefficients [Reinhart and Herkel 2000], but it cannot take advantage of mkillum in its calculations. All light passing into the space must be accounted by the single-step backwards ray-tracing approach in Radiance, making venetian blinds and other CFS a serious challenge. Particularly if the CFS is operable, the basic DC approach can become quite expensive. In Section 3.2, we describe a new DC method that separates the effect of the BSDF, making annual simulations truly practical, even with complex, operable fenestration.
In the following sections, we describe two techniques for incorporating BSDF data in Radiance. The first technique enhances mkillum so it can precompute light output from CFS portals (windows or skylights) using BSDF data, permitting a broader array of systems to be modeled. The second technique extends the DC approach using the rtcontrib program to compute coefficients for incoming and outgoing directions on CFS with BSDF data connecting the two. This offers advantages for annual simulation and operable shading devices, as changing the control setting is as simple as swapping matrices in an inner time-step loop.
3.1 SINGLE TIME-STEP EVALUATIONS AND RENDERINGS OF CFS IN RADIANCE USING MKILLUM
The main challenge with CFS is that their output is a) difficult to compute, and b) varies over their surface. While it is possible to rely upon the default Radiance interreflection calculation to integrate light arriving through the windows and skylights, the more complex their behavior, the more time it takes to get a reasonable result. Since we know that windows are important contributors to illumination, we can turn them into secondary light sources, precomputing and smoothing their output distributions for efficient rendering. This is the purpose of mkillum.
Since Radiance 3.9 was released in 2008, mkillum has supported BSDF descriptions of window and skylights, using XML files exported from Window 6. A key advantage is that mkillum does not need to model how light gets through the window, since the BSDF informs us directly. All mkillum needs to determine is how light is impinging on the exterior opening.
The obvious way to compute the external incident radiation is to send sample rays out from the window into the exterior environment. Obstructions, interre-flections, overhangs, buildings and all will be accounted for by this external sampling. However, we need to make special allowances for light sources such as the sun, which constitute important, localized contributors to illumination. Random sampling is no way to integrate the sun, even if we can account for any complex redirection to the interior with a BSDF. We therefore send a user-selected number of samples from random points on each opening's exterior towards the sun and any other localized sources described in the scene model. The radiation from these samples will be accounted for separately and passed to the appropriate incoming direction of the BSDF. (2) What happens in the final stage depends on the setting of the mkillum '1' flag.
The '1' flag indicates whether the window being computed is a secondary source ('1-') or a standard light source ('1+'). In the latter case, all contributions are fed through the BSDF and resampled to create a luminous intensity distribution in a mkillum-specific polar coordinate system for treating the window as if it were a luminaire. This is not normally done for windows, since a direct view will show only a uniform output corresponding to the average behavior in the viewer's direction. Also, for any window system with an unobstructed beam component, the resolution of the BSDF and the recorded luminous intensity distribution will tend to smear the light path excessively, resulting in blurry and featureless shadows. This is why the default method computes a secondary light source (12), leaving the final rendering phase to create a realistic view out the window. Similarly, direct contributions passing through the window are also rendered, such as the detailed shadows cast by venetian blinds, which would be nearly impossible to produce using a BSDF [de Boer 2006].
In the secondary source ('1-') method just described, a direct beam through a window is accounted in the final rendering phase, so we must discount this contribution in the luminous intensity table. To accomplish this, we make a separate tally of all source contributions when accumulating our incoming radiation. In the final mkillum stage, as we pass this light through our BSDF and redistribute the output into a luminous intensity table, we subtract source light that passes straight through the BSDF and continues in the same direction. Light that is redirected by the BSDF is unchanged. The user then creates a geometric representation of the CFS for proper view and beam shading. The rendered result is an accurate depiction of the window geometry, the indirect contribution from the window, and the direct contribution passing through the window (Fig. 2).
[FIGURE 2 OMITTED]
While the calculation of luminance distributions for CFS characterized with BSDF data is fairly efficient using mkillum, it still takes some time to sample the exterior environment, account for interreflections and so on, then considerably longer to render a final image using the computed secondary light sources. The time is well-spent if the desired output is a high-quality rendering at a particular time and day, but may be impractical in an annual simulation context. A DC approach is required.
3.2 ANNUAL TIME-STEP CALCULATIONS OF CFS IN RADIANCE USING RTCONTRIB
The DC approach is a proven method for annual simulation that subdivides the sky into patches whose partial contributions are computed independently [Tregenza 1983]. Any general lighting simulation method may be used to compute the daylight coefficients in advance of the annual loop. The inner time-step then assigns a set of values to the sky patches for a particular day/time/weather condition, multiplying this against the precomputed daylight coefficients and summing the partial contributions. This vector dot product is quick to evaluate, making the DC precomputation well worth the effort, but the standard approach must be altered substantially to handle BSDF data and operable systems.
We have the same issues in a DC calculation as we do in a single time-step rendering, which is the complex interaction at the windows. Unfortunately, mkillum cannot help us here, because it only computes the output distribution for a particular daylight condition. We want to know what happens under any condition, which would require a separate run of mkillum for every sky patch. While certainly possible, such a calculation is very inefficient and would take too long to complete. This is where we introduce the rtcontrib tool.
The rtcontrib tool can compute light flux transfer between almost arbitrary sources and destinations using a general geometric ray-tracing method. In the context of daylight coefficients and CFS, we use rtcontrib to calculate a matrix linking sky patches to incident directions on a particular window or skylight group. We call this matrix D for "daylight," and we compute separate D matrices for each fenestration group, where groups are divided by elevation (orientation) and proximity to external obstructions, an issue we explore more fully in Section 4. From matrix D and a set of sky patch luminances, we can compute the luminance arriving from specified directions on each window group, as all daylight illumination ultimately comes from the sky. (3) This is illustrated in Fig. 3b-3d, which show a Tregenza sky subdivision outside the window shaded pink in Fig. 3a through a Klems basis. BSDF patches that see only sky will get a contribution from a few sky patches, whereas BSDF patches that see other surfaces in the scene will receive an integration of many sky patches. The D matrix thereby incorporates all daylight contributions, direct and indirect, accounting for diffuse and specular reflections, obstructions, and so on In addition, we provide a single ground patch for the distant horizon (shown in Fig. 3d), since we cannot specify an infinite ground plane as Radiance geometry.
For a given sky condition, we could quickly compute the outgoing luminance on the window interior from a BSDF matrix. We could use these luminous intensity distributions to render an image directly, but as this is also part of our time-step calculation, it would be much better to precompute this transfer function as well. A second run of rtcontrib is used to relate the exiting directions on our window interiors to the desired image or sensor values, and we call this the V or "view" matrix. This matrix is illustrated in the case of an interior view, shown as projected colors in Fig. 4.
[FIGURE 4 OMITTED]
Identifying the BSDF matrix as T, we can write our combined matrix formula for a single window group thus:
i =VTDs (3)
i=result vector (luminances, illuminances, image, etc.)
V=view matrix, relating outgoing directions on window to desired results at interior
T=transmission matrix, relating incident window directions to exiting directions (BSDF)
D=daylight matrix, relating sky patches to incident directions on window
s = sky vector, assigning luminance values to patches representing sky directions.
The dimensions of the vectors and matrices will depend on the desired resolution in the calculation. If we are computing 100 interior illuminances from a Tregenza sky subdivision with 146 patches (including ground), then the s vector will be 146x1, the D matrix will be 145x146, the T matrix using a full Klems BSDF basis will be 145x145, the V matrix will be 100x145, and the resultant i vector will be 100x1.
When multiple window groups are involved, we evaluate vector i for each group and sum the results together. For operable fenestration systems, the above formula permits us to quickly swap the T matrix for an alternate position (for example, slat angle) and iterate until we arrive at the desired result vector i. In cases where the CFS is not operable, we can collapse the V, T, and D matrices into a single rectangular daylight coefficient matrix for the sake of efficiency. The important point is that we have eliminated all ray-tracing calculations from the inner time-step loop, converting it into a straightforward matrix-vector multiplication. This technique has been used successfully in a pilot study of 60 real buildings [Saxena and others, 2010].
Figure 5 shows our window office again, this time rendered using rtcontrib. Precalculation of the D and V matrices took about an hour on a laptop computer, followed by 12 seconds to evaluate the i vector (1.5 megapixels long) for a particular date, time, and slat angle. Once again, we have lost information for the slat shadows, but in this case, we have no convenient remedy as we are relying completely on the BSDF to transfer external radiation to the interior. (The issue of BSDF resolution is addressed in the following section.)
[FIGURE 5 OMITTED]
4 VALIDATION AND APPLICABILITY OF SIMULATION METHODS
4.1 ASSIGNMENT OF FLUX IN EXITING DIRECTION WITHIN MKILLUM AND EFFECT ON RENDERINGS
To ensure that the Window 6 BSDF data coordinate system was properly interpreted within Radiance, we ran a simple check to validate input and output directions by reproducing BSDF output using mkillum. A simple test of an incident normal light ray that redirects light into specific Klems patches verified that zone directions are being correctly addressed by mkillum (and therefore by rtcontrib).
In evaluating the single patch output, which considered output in patches numbered 1, 5, 25, 30, 63, 72, 96, 107, and 138 (Fig. 6), it is clear that light expected to be redirected into a single patch is distributed over multiple patches (Fig. 7). This is due to the differences between the Klems basis and polar representation of the distribution that is created and applied by mkillum: incoming flux from a Klems patch must be split up and assigned to outgoing patches of the polar basis, which differ slightly in angular position. Neighboring patches which had BTDF values of zero are receiving some of the flux from adjacent patches that contained individual spikes in the distribution.
[FIGURE 6 OMITTED]
The '1+' and '1-' options also provide different results. As described above, with the mkillum '1+' option, the specular component is determined from the BSDF data. Using the '1-' option, the specular component is removed and is addressed using a description of the material's geometry which is placed within a zone of designated thickness adjacent to the planar BSDF material. Figure 7 shows the differences in performance for a small rectangular aperture facing the inside of a hemisphere. In this example, the eight patches with nonzero values cause some light to be distributed beyond their boundaries. The central patch of the polar representation is much larger than the Klems patch 1.
[FIGURE 7 OMITTED]
[FIGURE 8 OMITTED]
To address performance with real materials, validation tests were conducted using mkillum with the '1+' and '1-' settings for a horizontal blind (calculated in Window 6) and a fabric shade (measured with a goniophotometer). Both shading devices permit direct sunlight to penetrate at the specular angle. The results of these tests show that mkillum is correctly able to track the amount of flux admitted to the space quite accurately under both the '1+' and '1-' approaches, but is unable to properly track the flux direction when there is a sharp spike in the distribution and the '1+' option is applied. The flux that passes directly through the material and exits within one of the 145 patches is often distributed across multiple patches.
In the case of the horizontal blinds (1+), the flux in patch 4 is approximately 21 percent of its theoretical value, while in patch 1-3 and 5-9 it is roughly seven to twelve times the theoretical value determined using the BTDF, and in patches 13, 14 and 15, it is 37 to 87 times higher (Fig. 8a). In all cases, the total flux emitted by the material over all directions is still very close to the theoretical value. When the direct flux is removed from the BSDF data using the '1-' option, the total emitted flux and the directional information are quite accurate.
Examples of the theoretical and mkillum derived center zone luminous intensity and total flux for a horizontal blind and fabric shade material with open weave are shown in Fig. 8 for the two different mkillum options. Separate tests at angles where a concentrated direct source beam is completely blocked by the shade, as well as a test where a uniformly luminous hemisphere illuminated the BTDF material, showed excellent agreement across all patches. Since the basis definition for incoming and exiting flux is the same in rtcontrib, the spread of flux to adjacent patches seen in mkillum is not a concern in rtcontrib.
[FIGURE 9 OMITTED]
[FIGURE 10 OMITTED]
The results of these tests show that BSDFs can be applied with mkillum, but beam contributions that pass directly through shading devices should be separated out from the BSDFs using the '1-' option to permit this flux to be properly addressed. Very narrow spikes in the output distribution resulting from specular reflections within a CFS are unlikely to be treated accurately by either approach, since flux will be distributed across multiple patches as shown in Fig. 7.
4.2 RESOLUTION OF THE BSDF AND SKY SUBDIVISION AND EFFECTS ON ACCURACY
Up to now, we have discussed hemispherical subdivision schemes with 145 patches (full Klems BSDF basis and Tregenza for sky subdivisions), but because values are averaged over patches, the resolution of the hemispherical division schemes affects accuracy. A resolution of 145 patches equates to an average solid angle of 0.0433 sr per division (cone with 13.5[degrees] apex angle). For comparison, the sun has a 0.5[degrees] apex angle.
For the sky subdivision, the Radiance program used to create the s vector spreads the radiance of the sun over three Tregenza sky divisions (a 0.5[degrees] source spread over three 13.5[degrees] patches) to maintain the appropriate centroid. In addition, the sky divisions do not align with Klems divisions causing the solar radiance to be spread over six or more Klems patches. If the simulated CFS has an output peak corresponding to sunlight transmission or re-direction than what may in reality be a focused patch of light may be spread up to 20[degrees] off axis. To reduce this effect, a higher resolution subdivision of the sky such as the Reinhart scheme, can be used which further subdivides the Tregenza patches while preserving a zenith patch (Fig. 9).
[FIGURE 11 OMITTED]
Incorrect flux placement can also occur when a low-resolution BSDF is used to represent specularly-reflecting CFS with peaky output. Averaging over the Klems 145-patch basis reduces the intensity and spreads the narrow peak over a larger angle. The redirected sunlight would cover a larger portion of the ceiling and be less bright at any single location.
To better understand impacts, high-resolution angle bases were created to better understand impacts on simulation accuracy (Fig. 10). A variable resolution angle basis would most efficiently characterize a peaky optical redirection system. The division schemes depicted are fixed resolution schemes, but might represent the finest level of resolution for a variable division scheme. Images in Fig. 11 illustrate the size of BSDF patches for three bases resolutions.
[FIGURE 12 OMITTED]
[FIGURE 13 OMITTED]
Images in Fig. 12 demonstrate the spread of direct flux (sunlight) for a space with clear glass using various resolutions of BSDF and sky divisions for the mkillum and rtcortt rib tools. The rpict rendering represents "ground truth" or the reference condition against which the percentage difference in luminance levels are indicated in Fig. 13. (4) Note how the differences are greatest for both the view portion of the window wall and interior space for the low-resolution cases and decreases as resolution is increased. Resolution of small-area glare sources and the luminance patterns of the windows are also affected, skewing assessments of visual discomfort. The apparent window luminance is averaged over outgoing Klems patches, effectively replacing small intense sources like the sun orb with a larger, less intense glare source (Fig. 14).
[FIGURE 14 OMITTED]
In some cases, an averaged result may be more desirable than a discontinuous pattern. For example, when evaluating work plane illuminance, the real world result depends on the position of the simulation point within the shadow pattern. Whether a high sunlight illuminance or a low shadow illuminance is reported is left to chance. The result from an averaged continuous BSDF output more closely resembles the average illuminance on the work plane. An averaged result may also be more desirable for a viewpoint-specific glare assessment. In the nonBSDF case, the uncertainty of reporting either a higher degree of discomfort or no discomfort is based on the chance of the viewpoint being in sun or shadow. The average BSDF result will report discomfort less than the sun condition but more than the shadow condition. Clearly, simulations using higher-resolution datasets will yield greater accuracy, particularly for peaky systems, but impacts on annual energy and comfort assessments will need to be further investigated.
4.3 SUBDIVISION OF WINDOW GROUPS AND WINDOW
Use of rtcontrib requires intelligent partitioning of windows into groups to account for nearby obstructions such as an overhang, operable shades whose height varies on a time-step basis, and interactions both inside and outside a building. Inaccuracies can result if this subdivision is not well considered.
The D matrix contains DC contributions either for a single point on the facade or integrated over an area of the facade. To obtain a more accurate value of incident flux on the window, a simulation may need to be split into many parts with more than one D matrix to reproduce nonuniform shading effects caused by external obstructions. The number of D matrices required depends on the proximity and size of external obstructions (Fig. 15-17). Similarly, an attached building obstruction may require subdivision of individual windows. Subdividing a window with an overhang into horizontal bands, for example, reproduces the variation in shading over the height of the window provided by the overhang (Fig. 18-19).
[FIGURE 15 OMITTED]
[FIGURE 16 OMITTED]
[FIGURE 17 OMITTED]
[FIGURE 18 OMITTED]
[FIGURE 19 OMITTED]
Window groupings and subdivisions used for the V matrix should mirror those used for the D matrix. For some cases (that is, variable height shading systems), it is desirable to divide the windows further for the V matrix computation. Further division for V matrix computation is permissible as long as D matrix groups are reflected in V matrix subgroups. For example, variable-height shading systems, such as a roller shade, require at least two BSDFs to represent optical properties above and below the bottom of the shade - Fig. 20 shows a window divided into 4 equal bands enabling simulation of discrete shade height positions in 25 percent increments.
Currently, Radiance does not perform automatic grouping or subdivision of windows. The user is left to determine when window grouping and subdivision is necessary and to determine the appropriate level of subdivision for their simulation needs.
[FIGURE 20 OMITTED]
A new method of modeling the daylighting performance of complex fenestration systems (CFS) was implemented in Radiance. The method relies on bidirectional scattering distribution functions data output from Window 6, which enables end users to describe any arbitrary assembly of glazing, shading, and other optically-complex coplanar layers. Modifications to two of Radiance's existing tools, mkillum and rtcontrib, enable end users to accurately render an image of an interior space, including shadows cast by an inhomogenous system such as a Venetian blind, and perform efficient annual simulations. The modified tools will require judicious use in order to obtain accurate results. Further work is required to determine the sensitivity of annual and time-step evaluations of peaky systems when resolution of the BSDF basis and sky subdivisions is varied.
We gratefully acknowledge the contributions of Marilyne Andersen, Ecole Poly-technique Federale de Lausanne, toward the validation of this new capability. This work was supported by the Assistant Secretary for Energy Efficiency and Renewable Energy, Building Technologies Program, of the U.S. Department of Energy under Contract No. DE-ACO2-05CH11231. This research used the Lawrencium computational cluster resource provided by the IT Division at the Lawrence Berkeley National Laboratory (Supported by the Director, Office of Science, Office of Basic Energy Sciences, of the U.S. Department of Energy under Contract No. DE-ACO2-05CH11231).
[c] 2011 The Illuminating Engineering Society of North America doi:
Andersen M, Stokes E, Gayeski N, Browne C. 2009. Using digital imaging to assess spectral solar-optical properties of complex fenestration materials: A new approach in videogoniophotometry. Sol Energy. 84(4): 549-562. doi: 10.1016/j.solener.2009.02.005.
Andersen M, de Boer J. 2006. Goniophotometry and assessment of bi-directional photometric properties of complex fenestration systems, Energy Build. 38(7): 836-848.
Carroll WL, Hitchcock RJ. 2005. DElight2 daylighting analysis in EnergyPlus: Integration and preliminary user results. In: Proceedings of building simulation 2005. Montreal, Canada.
de Boer J. 2006. Modelling indoor illumination by complex fenestration systems based on bidirectional photometric data. Energy Build. 38(7): 849-868.
Jonsson JC, Rubin M, Nilsson AM, Jonsson A, Roos A. 2009. Optical characterization of fritted glass for architectural applications. Opt Mater. 31(6): 949-958.
Kaempf JH. 2004. Integration of BT(R)DF data into Radiance lighting simulation programme. Technical report CTI Project 4881.1 Bidrectional Goniophotometer. Ecole Polytechnique Federale de Lausanne.
Klems JH. 1994a. A new method for predicting the solar heat gain of complex fenestration systems: I. Overview and derivation of the matrix layer calculation. ASHRAE Trans. 100(1): 1065-1072.
Klems JH. 1994b. A new method for predicting the solar heat gain of complex fenestration systems: II. Detailed description of the matrix layer calculation. ASHRAE Trans. 100(1): 1073-1086.
Klems JH, Warner JL, Kelley GO. 1997. A comparison between calculated and measured SHGC for complex glazing systems. ASHRAE Trans. 102(1): 931-939.
Mitchell R, Kohler C, Klems J, Rubin M, Arasteh D, Huizenga C, Yu T, Curcijaf D, editors. 2008. Window 6.2/ Therm 6.2 Research Version User Manual. Lawrence Berkeley National Laboratory. LBNL-941. Berkeley, CA 94720.
Nicodemus FE. 1965. Directional Reflectance and Emissivity of an Opaque Surface. Appl Opt. 4: 767-773.
Reinhart CF. 2001. Daylight availability and manual lighting control in office buildings - Simulation studies and analysis of measurements. Thesis. University of Karlsruhe.
Saxena M, Ward G, Perry T, Heschong L, Higa R. 2010. Dynamic Radiance - Predicting annual daylighting with variable fenestration optics using BSDFs. In: Proceedings of SimBuild 2010. 4th National Conference of IBPSA-USA. 2010 August 11-13. New York.
Schregle R, Weinold J. 2004. Physical validation of global illumination methods: Measurement and error analysis. Comput Graph Forum. 23(4): 761-81.
Tregenza PR. 1983. Daylight coefficients. Light Res Tech. 15(2): 65-71.
Larson, Ward G, Shakespeare R. 1998. Rendering with radiance: The art and science of lighting visualization. San Francisco. Morgan Kaufmann. 656 p.
(1.) Klems modified the Tregenza hemispherical subdivision to give higher resolution in incident angle and a weighting of the patches proportional to their solid angle and projected area. The angular resolution of this basis is approximately -[+ or -] 5[degrees]in incident angle and much coarser in azimuth.
(2.) If a source is large or local, its energy may enter multiple patches in the BSDF basis. If one of our previous random samples happens to hit a light source, it is given zero energy to avoid double-counting.
(3.) For the time being, the sun is assumed to have its energy distributed to a few sky patches in our analysis.
(4.) These examples use clear glass which can be defined natively by Radiance and many other simulation software packages. Using a BSDF to simulate clear glass is unnecessary, however is done here to illustrate errors associated with a peaky BSDF that can be understood intuitively.
G. Ward, Anyhere Software, 950 Creston Road, Berkeley, CA, 94708 USA; R. Mistrick, Ph D., Penn State University, 104 Engineering A, University Park, PA 16802 USA; E. S. Lee, A. McNeil, J. Jonsson, Ph.D., Building Technologies Program, Environmental Energy Technologies Division, Lawrence Berkeley National Laboratory, Mailstop 90-3111, 1 Cyclotron Road, Berkeley, CA 94 720 USA
|Printer friendly Cite/link Email Feedback|
|Author:||Ward, G; Mistrick, R.; Lee, E.S.; McNeil, A.; Jonsson, J.|
|Date:||Apr 1, 2011|
|Previous Article:||X-Ray computer tomography in end-of-life investigations of HID lamps.|
|Next Article:||Vehicle kinematics in turns and the role of cornering lamps in driver vision.|