LICASO and DAYSIM

Ian Ashdown1, Chris Jackson2, Joel Spahn3, Todd Saemisch3

1. SunTracker Technologies Ltd., Victoria, Canada. 2. Lighting Analysts Ltd., London, UK. 3. Lighting Analysts Inc., Littleton, CO.

LICASO™ and DAYSIM are daylighting analysis software programs that perform climate-based annual daylight simulations, including the calculation of annual daylight metrics. This white paper compares their performance in terms of accuracy and calculation times. Average illuminance and spatial daylight autonomy values agree to within 8 percent, while LICASO is approximately 250 times faster than DAYSIM for a benchmark model.

1.  Climate-Based Daylight Metrics

Quoting from the IES Lighting Handbook, Tenth Edition,5 the definition of daylight autonomy seems simple enough: “… the measure of the percentage of the operating period (or number of hours) that a particular daylight level is exceeded throughout the year.” It includes spatial daylight autonomy (sDA)6, which “… reports the percentage of sensors (or building area) that achieves a minimum daylight illuminance level (typically 300 lux) for a minimum percent of the analysis year (time).” Other dynamic daylight metrics include continuous daylight autonomy (cDA)6, maximum daylight autonomy (mDA)6, useful daylight illuminance (UDI)9 and more, all with seemingly simple definitions.

Simple to say, yes, but calculating these metrics is another matter entirely. Behind the scene are devilishly complex algorithms that require massive amounts of computation. Until recently, the only options for professional lighting designers and architects have been based on the justly acclaimed Radiance suite of lighting simulation tools. A typical example is DAYSIM, which is described as “… validated, Radiance-based daylighting analysis software that models the annual amount of daylight in and around buildings.”

LICASO is the first climate-based annual daylight simulation software program that is not based on Radiance. More than this, it is not even based on the Radiance computational model of ray tracing. Rather, it relies on proven radiosity methods,1 and in particular the algorithms that have been driving Lighting Analysts’ AGi32 and ElumTools lighting design and analysis software products for nearly two decades. (See Lighting Analysts’ blog article Climate-Based Daylight Modeling for further details.)

… but enough gratuitous advertising. The intent of this white paper is to compare the performance of LICASO and DAYSIM in terms of accuracy and calculation times. This is not a case of which program is better suited for any given application, but simply to see whether the two programs are indeed comparable.

2.  Benchmark Model

The spatial daylight autonomy metric has been adopted for use with green building certification by the US Green Building Council10 and the International WELL Building Institute.3 However, the United Kingdom Education Funding Agency mandates the use of both spatial daylight autonomy and useful daylight illuminance daylight metrics for its Priority School Building Programme.4 The benchmark model used in this white paper is therefore based on a typical 55-m2 classroom in accordance with the PSBP baseline design.

The benchmark model consists of four identical rooms facing north, south, east, and west, with each room having two glazed windows (FIG. 1). Each room measures 7.5 meters long by 7.0 meters wide by 3.2 meters high.

Each window measures 1.68 meters wide by 2.0 meters high, is positioned 0.85 meters above the floor and 0.5 meters from the closest wall, and has a transmittance of 70 percent.

The floor reflectance is 20 percent, the wall reflectance is 60 percent, and the ceiling reflectance is 80 percent.

A virtual ground plane with 18 percent reflectance is assumed for LICASO. The equivalent for Radiance (and hence DAYSIM) is a 180-degree glow source (essentially an upside-down sky) with uniform luminance that is the horizontal illuminance due to the diffuse skylight and direct sunlight multiplied by the ground plane reflectance.8 DAYSIM allows the user to specify the ground plane reflectance, with a default value of 20 percent.

Figure 1 - Benchmark Model

Figure 1 – Benchmark Model

A grid of 15 by 15 virtual photometers, spaced at 0.5-meter intervals, is centered in each room, with a mounting height of 0.75 meters (FIG. 2).

Figure 2 - Virtual photometer layout

Figure 2 – Virtual photometer layout

3.  Simulation Parameters

The TMY3 weather file is LONDON/GATWICK-GBR (37760), with a site location of 51.15 degrees north and 0.18 degrees west. All simulations were run in one-hour increments for the entire year, with occupied hours of 8:00 AM to 4:00 PM for a total of 2,920 hours, and Daylight Saving Time from March 29th to October 25th.

3.1    DAYSIM

DAYSIM uses a modified version of the Radiance utility program rtrace called rtrace_dc., where the “dc” suffix represents “daylight coefficients.” Of the 46 user-specified parameters available for rtrace, rtrace_dc provides user access to 13 of them with default values (Table 1).

Table 1 – DAYSIM User-Specified Parameters

Parameter Name Default Value
-aa acc ambient accuracy 0.10
-ab N ambient bounces 5
-ad N ambient divisions 1000
-ar res ambient resolution 300
-as N ambient super-samples 20
-dj frac source jitter 0.0
-dp D direct pretest density 512
-dr N direct relays 2
-ds frac source sub-structuring 0.20
-lr N limit reflection 6
-lw frac limit weight 0.004
-sj frac specular jitter 1.0
-st frac specular threshold 0.15

It is not obvious to DAYSIM users what effect these parameters will have on the calculations. However, the same parameters are available for the Radiance utility program rpict, and so reference can be made to http://radsite.lbl.gov/radiance/refer/Notes/rpict_options.html, with rendering artifacts related to the relevant parameters enumerated in Table 2.

Table 2 – Artifacts Associated with DAYSIM Parameters

Parameter Artifact Solution
-aa uneven shading boundaries in shadows decrease value by 25%
-ab lighting in shadows too flat increment value
-ad “splotches” of light double value
-ar shading wrong in some areas double or quadruple value
-as “splotches” of light increase to half of -ad setting
-dj shadows are unnaturally sharp increase value to 0.7
-dp incorrect mirror reflections double value
-dr missing multiple mirror reflections increment value
-ds large sources cast unnatural shadows decrease value by 50%
-lr some multiple specular reflections gone increment value
-lw some specular reflections gone decrease value by 50%

Some of these parameters are problematic, of course, in that their effects can only be seen in the renderings generated by rpict. Without access to these renderings, DAYSIM users have little choice but to accept its default values.

More obvious are the effects of the parameter values on the calculation times. These are (again from the rpict documentation) enumerated in Table 3.

Table 3 – Calculation Times Associated with DAYSIM Parameters

Parameter Execution Time Effect
-aa direct, doubling this value approximately quadruples rendering time
-ab direct, doubling this value can double rendering time
-ad direct, doubling value may double rendering time
-ar direct, effect depends on scene, can quadruple time for double value
-as direct, effectively adds to -ad parameter and its cost
-dj indirect, increasing value requires -ps parameter to be reduced
-dp minor, affects start-up time only, higher values take longer
-dr direct, depending on the scene each new reflection can double time
-ds inverse, halving value causes rendering time to approximately double
-lr minor, increase causes very slightly longer rendering time
-lw minor, decrease causes very slightly longer rendering time

 (The -ps parameter of rtrace is not accessible to the user, so presumably rtrace_dc modifies this parameter accordingly when the -dj parameter is changed from its default value.)

What is clear is that the -a* parameters can be very expensive in terms of calculation time, and should therefore be changed with considerable caution. In the absence of rpict renderings, however, the only indication of the effect of these parameters is on the uniformity of the virtual photometer readings.

To illustrate this point, FIG. 3 shows the isolux distribution of photometer readings for two sets of -a* parameter values, with all other parameters being set to their DAYSIM default values.

Figure 3 – Visualization of DAYSIM –a* parameter values effect on photometer uniformity

Figure 3 – Visualization of DAYSIM –a* parameter values effect on photometer uniformity

Simply by looking at the isolux distributions, it is evident that setting -ad to 1000 and -as to 20 results in “splotches” of light. Raising these parameter values to 5,000 and 500, respectively, appears to resolve this issue, but at the expense of increasing the calculation time. (In this particular example, the DAYSIM calculation time increased by a factor of 6.5 times.)

It should also be noted that the optimal -a* parameter values are scene-dependent. Mastery of these parameters requires an in-depth understanding of how Radiance interpolates its cached irradiance values.

DAYSIM further offers three options for daylight coefficients:

  • Original with 65 representative direct solar positions (e.g., FIG. 4)
  • DDS (Dynamic Daylight Simulations) with 2,305 representative direct solar positions
  • Shadow testing with hourly direct solar positions taken directly from TMY3 weather file records

In the third option, the actual solar position is bi-linearly interpolated from the representative direct solar positions for the first two options. (The third option is reportedly rarely used because it is very expensive in terms of calculation time.)

FIG. 4 –Annual solar path (65 positions) for Freiberg, Germany. (Ref. 2)

FIG. 4 –Annual solar path (65 positions) for Freiberg, Germany. (Ref. 2)

3.2    LICASO

In accordance with radiosity methods, LICASO subdivides each surface into a two-level hierarchy of patches and elements (Fig. 5). For each interreflection, light is received by the elements of a patch and then reflected (and transmitted for translucent surfaces) from the center of the patch.1 The combined direct sunlight and diffuse daylight are interreflected (or “bounced”) between elements and patches in this way until mostly absorbed.

Compared to DAYSIM, LICASO has only four user-specified parameters for climate-based annual daylight modeling:

  • Maximum surface patch area
  • Maximum window patch area
  • Number of elements per surface patch
  • Stopping criterion for absorbed light

For the benchmark model, the maximum patch area is 1.0 m2, the number of elements per patch is four, and the stopping criterion is 99 percent. (That is, the bounces of light stop when 99 percent of the interreflected light is absorbed.)

Direct sunlight and diffuse skylight incident upon the windows are received by each window patch and transmitted into the room interior from the center of each patch. For the benchmark model, each window patch has an area of approximately 0.3 m2.

Figure 5 – LICASO surface discretization into patches (blue lines) and elements (red lines).

Figure 5 – LICASO surface discretization into patches (blue lines) and elements (red lines).

LICASO defines 120 representative direct solar positions, as shown in FIG. 6, where the representative positions are calculated for each hour on the specified dates. The actual solar position for any given hour and date is then linearly interpolated from the representative direct solar positions for the same hour. (See Lighting Analysts’ blog article Climate-Based Daylight Modeling for further details.)

Figure 6 – LICASO representative direct solar positions.

Figure 6 – LICASO representative direct solar positions.

LICASO generates a variety of daylight metrics:

  • Illuminance
  • Basic Daylight Autonomy (DA)
  • Continuous Daylight Autonomy (cDA)
  • Maximum Daylight Autonomy (mDA / maxDA)
  • Minimum Daylight Autonomy (minDA)
  • Spatial Daylight Autonomy (sDA)
  • Useful Daylight Illuminance (UDI)
  • Annual Daylight Exposure (ADE)
  • Annual Sunlight Exposure (ASE)
  • Spatial Annual Sunlight Exposure (sASE)

It also provides three-dimensional rendered views (e.g., FIG. 7) and animations of single days and the entire year. This enables the user to both analyze and visualize the distribution of daylight throughout the year.

Figure 7 – LICASO Daylight Autonomy.

Figure 7 – LICASO Daylight Autonomy.

4.  Results

All tests were performed on a Windows 10 desktop computer with an Intel Core i7-4770K quad core CPU (3.5 GHz overclocked at 4.1 GHz) and 32 GB of random access memory

As previously noted, both the accuracy and execution time of DAYSIM is strongly dependent on the user-specified parameter values, as is evident from FIG. 3. Consequently, eight separate simulations were performed with different parameter settings for -ab, -ad, -ar, and -as, with each simulation being compared with the LICASO results. (Default values were used for all other DAYSIM parameters.)

Two metrics were chosen for comparison purposes: average illuminance as measured by the virtual photometers, and spatial daylight autonomy for 300 lux and 50 percent minimum time (designated as sDA300/50% by IES LM-83-12).6 The benchmark results are presented in Appendix A.

The DAYSIM versus LICASO average illuminance differences are plotted in Figure 8. Assuming that the DAYSIM simulations represent increasing accuracy with each simulation, it is evident that LICASO underestimates the average illuminance of the south room by 5 percent and the west room by 3 percent, and overestimates the average illuminance of the north room by 2 percent and the east room by 1 percent. Considering that DAYSIM and LICASO use completely different computational models, these differences are remarkably small.

Figure 8 - DAYSIM versus LICASO average illuminance differences

Figure 8 – DAYSIM versus LICASO average illuminance differences

The DAYSIM versus LICASO sDA300/50% differences are plotted in Figure 9. Again assuming that the DAYSIM simulations represent increasing accuracy with each simulation, it is evident that LICASO underestimates the sDA of the west room by 8 percent, the north room by 4 percent, and the east room by 8 percent.

Figure 9 - DAYSIM versus LICASO sDA300/50% differences

Figure 9 – DAYSIM versus LICASO sDA300/50% differences

The DAYSIM sDA value for the east room appears to be a calculation anomaly, possibly due to a remaining “splotch” of light. Ray tracing in Radiance is a stochastic (i.e., random) process, and so this anomaly may not occur if the benchmark is executed on a different machine.

The execution times for the different simulations are summarized in Table 13.

Table 13 – DAYSIM / LICASO execution times

Simulation DAYSIM Execution Time

(minutes)

LICASO Comparison

(45 seconds)

1   28   37 times
2   69   92 times
3 103 137 times
4 143 190 times
5 176 234 times
6 182 242 times
7 221 294 times
8 261 348 times

The differences in execution time between LICASO and DAYSIM are perhaps surprising, but they are typical due to differences between the ray tracing and radiosity calculation models. Simply put, radiosity methods are better able to take advantage of scene redundancy between hourly calculations.

5.  Conclusions

For both the average illuminance and spatial daylight autonomy metrics used in this benchmark comparison, it must be implicitly assumed that DAYSIM generates correct values. DAYSIM is described as “validated, Radiance-based daylighting analysis software,” but the simulations show that the two metrics converge to constant values only for Simulations 6 through 8. It is true that Radiance has been validated by a number of studies, but the accuracy of its photometric predictions is highly dependent on the parameters chosen for rtrace and, by extension, rtrace_dc.

Simulation 1 generates average illuminance results that differ by up to 6 percent from the converged values of Simulations 6 through 8. Similarly, the sDA results differ by up to 13 percent. Given this, the differences in results between DAYSIM and LICASO for Simulations 6 through 8 are arguably acceptable. (As an aside, differences of ±10 percent between predicted and measured illuminances are considered quite acceptable in electric lighting calculations.)

Regarding the difference in calculation times – LICASO is hundreds of times faster than DAYSIM – this must be put into perspective. For the past three decades, Radiance has been the gold standard for electric lighting and daylighting research, and DAYSIM has built upon this foundation by offering lighting researchers open-source software for dynamic daylight metrics, annual visual glare analysis, and electric lighting control. The innumerable user-specified parameters of rtrace and other Radiance tools (including DAYSIM) may make them difficult to master, but they are essential for lighting research.

Lighting Analysts’ LICASO, by comparison, is a commercial product that is powered by proprietary software licensed from SunTracker Technologies, and which relies on patented and patent-pending algorithms. It is intended for use as a climate-based daylighting simulation and analysis tool for professional lighting designers and architects. It further does not support glare analysis, electric lighting control, or bidirectional scattering-distribution functions (BSDFs), although these features are currently under development.

In summary, then, this benchmark analysis has shown that DAYSIM and LICASO generate comparable results in terms of dynamic daylight metrics such as spatial daylight autonomy. LICASO is clearly faster, but this comes at a cost for daylighting research, as there are fewer parameters to experiment with. Which software to choose depends, as always, on the user’s requirements.

References

  1. Ashdown, I. 1994. Radiosity: A Programmer’s Perspective. New York, NY: John Wiley & Sons.
  2. Bourgeois, D., C. F. Reinhart, and G. Ward. 2008. “Standard Daylight Coefficient Model for Dynamic Daylighting Simulations,” Building Research & Information 36(1):68-82.
  3. Delos Living. 2017. The WELL Building Standard v1 with January 2017 addenda. New York, NY: Delos Living LLC.
  4. 2014. EFA Daylight Design Guide, Lighting Strategy, Version 2. London, UK: The National Archives.
  5. 2010. The Lighting Handbook, Tenth Edition. New York, NY: Illuminating Engineering Society.
  6. 2013. IES RP-5-13, Recommended Practice for Daylighting Buildings. New York, NY: Illuminating Engineering Society.
  7. 2012. LM-83-12, IES Spatial Daylight Autonomy (sDA) and Annual Sunlight Exposure (ASE). New York, NY: Illuminating Engineering Society.
  8. Mardaljevic, J. 1998. “Daylight Simulation,” Chapter 6, in Rendering with Radiance, G. W. Larson and R. Shakespeare, Eds. San Francisco, CA: Morgan Kaufmann Publishers.
  9. Nabil, A., and J. Mardaljevic. 2006. “Useful Daylight Illuminances: A Replacement for Daylight Factors,” Energy and Buildings 38(7):905-913.
  10. 2013. LEED v4 BD+C: Schools – Daylight. Washington, DC: U.S. Green Building Council (www.usgbc.org).

Appendix A – Benchmark Results

Table A1 – DAYSIM  / LICASO Simulation 1

Time
DAYSIM 28 minutes
LICASO 45 seconds

 

Parameter Value
-ab 5
-ad 1000
-ar 300
-as 20

 

Room Average Illuminance Difference sDA300/50% Difference
DAYSIM LICASO DAYSIM LICASO
West 987.6 983 -0.5% 70.9% 71.0% +0.1%
South 2060.3 1991 -3.3% 100.0% 100.0% 0.0%
North 612.2 648 +5.8% 62.1% 65.0% +2.9%
East 1219.1 1253 +2.7% 82.4% 87.0% +4.5%

Table A2 – DAYSIM  / LICASO Simulation 2

Time
DAYSIM 69 minutes
LICASO 45 seconds

 

Parameter Value
-ab         6
-ad 2000
-ar    300
-as    200

 

Room Average Illuminance Difference sDA300/50% Difference
DAYSIM LICASO DAYSIM LICASO
West 1002.0    983 -1.8%   73.1%   71.0% -2.8%
South 2081.6 1991 -4.3% 100.0% 100.0%   0.0%
North    625.2    648 +3.6%   64.8%   65.0% +0.3%
East 1236.0 1253 +1.3%   93.4%   87.0% -6.8%

Table A3 – DAYSIM  / LICASO Simulation 3

Time
DAYSIM 103 minutes
LICASO 45 seconds

 

Parameter Value
-ab         6
-ad 3000
-ar    300
-as    300

 

Room Average Illuminance Difference sDA300/50% Difference
DAYSIM LICASO DAYSIM LICASO
West 1009.6    983 -2.6%   75.8%   71.0% -6.3%
South 2085.7 1991 -4.5% 100.0% 100.0%   0.0%
North    632.6    648 +2.4%   67.0%   65.0% -2.9%
East 1240.3 1253 +1.0%   92.3%   87.0% -5.7%

Table A4 – DAYSIM  / LICASO Simulation 4

Time
DAYSIM 143 minutes
LICASO 45 seconds

 

Parameter Value
-ab         6
-ad 4000
-ar    300
-as    400

 

Room Average Illuminance Difference sDA300/50% Difference
DAYSIM LICASO DAYSIM LICASO
West 1011.4    983 -2.8%   75.3%   71.0% -5.7%
South 2092.0 1991 -4.8% 100.0% 100.0%   0.0%
North    631.3    648 +2.6%   65.9%   65.0% -1.3%
East 1243.2 1253 +0.7%   95.1%   87.0% -8.5%

Table A5 – DAYSIM  / LICASO Simulation 5

Time
DAYSIM 176 minutes
LICASO 45 seconds

 

Parameter Value
-ab 6
-ad 5000
-ar 300
-as 500

 

Room Average Illuminance Difference sDA300/50% Difference
DAYSIM LICASO DAYSIM LICASO
West 1015.8    983 -3.2%   76.4%   71.0% -7.0%
South 2093.9 1991 -4.9% 100.0% 100.0%   0.0%
North    636.4    648 +1.8%   65.9%   65.0% -1.3%
East 1243.0 1253 +0.8%   95.1%   87.0% -8.5%

Table A6 – DAYSIM  / LICASO Simulation 6

Time
DAYSIM 182 minutes
LICASO 45 seconds

 

Parameter Value
-ab         7
-ad 5000
-ar   300
-as   500

 

Room Average Illuminance Difference sDA300/50% Difference
DAYSIM LICASO DAYSIM LICASO
West 1013.5    983 -3.0%   76.4%   71.0% -7.0%
South 2097.4 1991 -5.0% 100.0% 100.0%   0.0%
North    634.4    648 +2.1%   68.1%   65.0% -4.5%
East 1245.8 1253 +0.5%   95.1%   87.0% -8.5%

Table A7 – DAYSIM  / LICASO Simulation 7

Time
DAYSIM 221 minutes
LICASO 45 seconds

 

Parameter Value
-ab        7
-ad 6000
-ar    300
-as    600

 

Room Average Illuminance Difference sDA300/50% Difference
DAYSIM LICASO DAYSIM LICASO
West 1025.4    983 -4.1%   76.9%   71.0% -7.6%
South 2096.2 1991 -5.0% 100.0% 100.0%   0.0%
North    633.3    648 2.3%   68.1%   65.0% -4.5%
East 1249.3 1253 +0.2% 100.0%   87.0% -13.0%

NOTE: The DAYSIM sDA value for the East room appears to be an anomaly.

Table A8 – DAYSIM  / LICASO Simulation 8

Time
DAYSIM 261 minutes
LICASO 45 seconds

 

Parameter Value
-ab         7
-ad 7000
-ar    300
-as    700

 

Room Average Illuminance Difference sDA300/50% Difference
DAYSIM LICASO DAYSIM LICASO
West 1016.1 983 -3.2% 77.5% 71.0% -8.3%
South 2095.9 1991 -5.0% 100.0% 100.0% 0.0%
North 634.2 648 +2.1% 68.1% 65.0% -4.5%
East 1242.6 1253 +0.8% 94.5% 87.0% -7.9%

Climate-Based Daylight Modeling

From Theory to Practice

Ian Ashdown, P. Eng., FIES

Chief Scientist, Lighting Analysts Inc.

[ Please send all comments to allthingslighting@gmail.com ]

daylight, n. The light of day.

Apart from having a wonderfully circular definition in most English-language dictionaries, daylight really is just another form of illumination. As such, most people would expect lighting designers to be able to simulate daylight with the same ease that we simulate electric lighting … but ah, I see you blushing.

We have for the past one hundred and fifty years relied on daylight factors to predict the distribution of daylight in architectural spaces. The daylight factor metric is exceedingly simple to calculate, but it is not very useful in understanding how daylight illuminates an interior space. All it really tells us (and our clients) is whether there will be sufficient daylight to read a newspaper indoors on an overcast day. We have, in other words, good reason to blush.

We know better, of course. Given the architectural plans of a building, including both its geographical coordinates and orientation, we know we can use historical weather data to determine the typical distribution of daylight within the building for every hour of every day of the year. We even have a name for this: climate-based daylight modeling (CBDM), an expression introduced a decade ago at a CIBSE lighting conference appropriately called “Engineering the Future” (Mardaljevic 2006).

With CBDM, we can calculate daylight metrics such as spatial Daylight Availability and Annual Sunlight Exposure (IES 2012), Useful Daylight Illuminance (Nabil et al. 2006), Daylight Glare Probability (Wienold et al. 2006), and more.  The first two metrics are important in that they are necessary for earning all three LEED v4 daylighting credit points (USGBC 2013). As consultants to architectural firms, lighting design professionals have an obligation to provide these metrics. Glare metrics take this one step further, offering the ability to identify potential design problems with large expanses of glazing.

Going further still, we can design and validate daylight harvesting systems for energy savings, and address building energy modeling issues involving solar insolation. Most important, we can work with architects during the conceptual design phase to take full advantage of what daylighting has to offer. From a lighting design perspective, CBDM expands the consulting services we can provide.

… and yet we persist in using daylight factors for our lighting design work.

The problem for most lighting designers is that the practice of climate-based daylight modeling is anything but simple. Until recently, the only software capable of performing CBDM calculations has been Lawrence Berkeley National Laboratory’s Radiance … and here I must pause.

Radiance

Radiance is – there are no other words to describe it – an exceedingly powerful, and indeed wonderful, set of software tools for electric lighting and daylighting research. It is not a monolithic program, but rather a set of one hundred or so Unix programs that can be linked together using command-line scripts.

Radiance, however, is first and foremost a research tool. It is reasonable to assume that most architects and engineers will prefer not to learn Unix, with command-line scripts such as:

CBDM - Unix Command-line Script

Fortunately, there are a number of free and commercial architectural-engineering design applications that are available, and which provide graphical user interfaces to the Radiance toolset. Those that support climate-based daylight modeling include DAYSIM and DIVA for Rhino (www.daysim.ning.com).

This article is not, however, about Radiance and its derivatives; it is about climate-based daylight modeling. More particularly, it is about a unique radiosity-based approach to CBDM calculations that does not involve the Radiance daylight calculation engine. It is the culmination of over twelve years of research and development, beginning with the paper “Modeling Daylight for Interior Environments” (Ashdown 2004). The details are disclosed herein for those interested in understanding how it works.

Follow the Light

In order to fully understand the radiosity approach, it is necessary to follow the light from source to receiver. The source is the combination of direct sunlight and diffuse daylight; the receiver is a room or similar naturally illuminated space within a building.

Modeling direct sunlight is straightforward. The solar position in the sky can be readily calculated from the equations presented in Section 7.1.5, Solar Position, of the IESNA Lighting Handbook, Tenth Edition (IES 2011). The solar disk is only 0.5 degrees in apparent width, and so it can be reasonably modeled as an infinitely distant point source that produces a beam of light whose rays are parallel.

Modeling diffuse daylight is more challenging. Some of the extraterrestrial radiation from the sun is scattered by the Earth’s atmosphere, resulting in the hemispherical diffuse light source that is the sky. The sky luminance (colloquially, “brightness”) spatial distribution varies with geographic location, site altitude, time of day, time of year, and weather conditions, including clouds and aerosols such as smoke and airborne dust, and also the dew point temperature.

It is clearly not practical to deterministically model realistic weather conditions, especially partly cloudy weather where the sky luminance distribution may vary on a time scale of minutes. Modeling diffuse daylight therefore requires a simplifying mathematical model, which in turn requires measured weather data.

Typical Meteorological Year

There are over 2,100 weather stations around the world (including 1,100 in North America) that measure weather data on an hourly basis. Through a complex set of empirical rules (Wilcox et al. 2008), thirty years or more of weather station data is compared on a per-month basis, and the hourly weather records for the twelve “most typical” months are assembled into a Typical Meteorological Year (TMY3) data set for the station’s geographic location[1].

Of the 68 elements in each hourly weather record, two are of primary importance for climate-based daylight modeling:

Element Unit Description
Direct normal irradiance Watt-hour per square meter Amount of solar radiation received in a collimated beam on a surface normal to the sun during the 60-minute period ending at the timestamp.
Diffuse horizontal irradiance Watt-hour per square meter Amount of solar radiation received from the sky on a horizontal surface during the 60-minute period ending at the timestamp.

Table 1 – TMY3 solar radiation elements. (Source: Wilcox et al. 2008)

Direct normal irradiance can be measured with a pyrheliometer, an instrument that measures the solar irradiance (including visible light and ultraviolet and infrared radiation from 300 nm to 2800 nm) incident upon its thermopile sensor (Muneer 2004). The device is always aimed directly at the Sun, and it includes a narrow tube that limits its field of view to six degrees (e.g., Figure 1).

FIG. 1 – Pyrheliometer. (Source: www.hukseflux.com)

FIG. 1 – Pyrheliometer. (Source: www.hukseflux.com)

Diffuse horizontal irradiance is usually measured with a pyranometer, an instrument (such as that shown in Figure 2) that measures irradiance from the sky incident upon its horizontal thermopile sensor (Muneer 2004). A shadow band may be positioned over the sensor to obscure a six degree-wide band following the path of the Sun, although it must be moved on a regular basis throughout the year. Alternatively (and more accurately), the global horizontal irradiance can be measured without a shadow band, and the measured direct normal irradiance measurement subtracted from it to determine the diffuse horizontal irradiance (ibid).

FIG. 2 – Pyranometer. (Source: www.hukseflux.com)

FIG. 2 – Pyranometer. (Source: www.hukseflux.com)

Perhaps surprisingly, these two measurements are all that are needed to model diffuse daylight.

Perez Sky Model

Based on some 16,000 full-sky scans made from Berkeley, California, Perez et al. (1993) proposed an empirical “all weather” sky model that predicts the absolute sky luminance distribution for weather conditions ranging from clear skies to totally overcast. The only two measured input parameters are direct normal and diffuse horizontal irradiance. (The model also includes the dew point temperature as a measure of atmospheric moisture content, but this has only a minor effect on the predicted luminance distribution.)

Other all-weather sky models have been proposed, but various validations studies (e.g., Noorian et al. 2008) have shown that the Perez sky model is better than most. More important, it has been implemented in the Radiance tool gendaylit to generate a single sky luminance distribution, and in gendaymtx to generate a set of hourly sky luminance distributions for the year from a TMY3 weather data file.

It should also be noted that where TMY3 weather data include direct normal and diffuse illuminance values, they have likely been calculated from the corresponding irradiance measurements using the Perez sky model. When direct normal and diffuse horizontal illuminance values are submitted to gendaylit, it uses an undocumented iterative algorithm to estimate the original measured irradiance measurements that the Perez sky model requires.

Sky Luminance Distribution

Prior to the introduction of calibrated all-sky digital cameras with fisheye lenses (e.g., Figure 3), mechanical scanners were used to measure the sky luminance distribution. Still manufactured by EKO Instruments, these instruments consist of a luminance meter mounted on an alt-azimuth platform with stepper motors, and measure the sky luminance at 145 different directions in about 4-1/2 minutes. They were previously used to obtain data for, among other purposes, the validation of various all-weather sky models, including the Perez sky model.

FIG. 3 – Digital all-sky camera. (Source: www.eko-usa.com)

FIG. 3 – Digital all-sky camera. (Source: www.eko-usa.com)

One legacy of these scanners has been the Tregenza sky subdivision (Tregenza 1987), wherein the sky dome is subdivided into eight 12-degree horizontal bands with 145 sky patches (Figure 4).

FIG. 4 – Tregenza sky subdivision. (Source: Muneer 2004)

FIG. 4 – Tregenza sky subdivision. (Source: Muneer 2004)

A particular advantage of this subdivision is that the sky luminance distribution can conveniently be represented as 145 discrete luminance values. Apart from the circumsolar region within a few degrees of the solar disk, the luminance of the sky dome in any direction can be interpolated with reasonable accuracy from these values.

Daylight Coefficients

Another advantage of the Tregenza subdivision scheme comes from a paper published over three decades ago, simply titled “Daylight Coefficients” (Tregenza et al. 1983). The researchers observed that each sky patch can be thought of as a separate and independent area light source that potentially illuminates a room through a window or opening (Figure 5). Using radiative transfer theory (aka radiosity), they demonstrated that – in theory – the luminance distribution in the room due to the intereflection of diffuse daylight from the sky patch (which they called “sky zones”) could be calculated.

FIG. 5 – Daylight coefficients and sky patches. (Source: Tregenza et al. 1983)

FIG. 5 – Daylight coefficients and sky patches. (Source: Tregenza et al. 1983)

Suppose, then, that each sky patch is assigned a luminance of 1000 cd/m2. If the luminance distribution in the room due to each patch is calculated and the results summed, the resultant luminance distribution is that due to a uniform sky[2] with a luminance of 1000 cd/m2. For the lack of any previous terminology, we can call this the canonical solution for the distribution of diffuse daylight in the environment (e.g., Figure 6).

FIG. 6 – Canonical diffuse daylight solution.

FIG. 6 – Canonical diffuse daylight solution.

Given, however, that each sky patch represents an independent light source, its resultant luminance distribution can be scaled by the average luminance of the sky patch for any given Perez sky model solution and building orientation. Summing these scaled luminance distributions therefore provides the luminance distribution of daylight in the room for the Perez sky model solution (e.g., Figure 7).

FIG. 7 – Example Perez sky model solution (overcast sky).

FIG. 7 – Example Perez sky model solution (overcast sky).

This is the key to climate-based daylight modeling. Given the architectural plans for a building, it is possible to calculate the canonical solution for the distribution of daylight in its rooms and other interior spaces. Then, given a TMY3 or similar weather dataset for the building site, the sky luminance distribution can be calculated for each daylight hour of the year using the Perez sky model , and with this the interior luminance distributions.

The advantage, of course, is that scaling and summing the contributions of each sky patch is much simpler and faster than calculating the canonical daylight solution. Even for complex environments with hundreds of thousands of polygonal elements, this can be done in milliseconds on a commodity desktop computer.

Again, however, this is in theory … in practice, the devil is very much in the details.

Direct Sunlight

The same approach can be applied to direct sunlight. Following an approach proposed by Bourgeois et al. (2008),  interior luminance distributions can be calculated for a selected number of solar positions and included with the canonical solution (although they do not appear in the renderings). For a given TMY3 weather record, the solar position can then be calculated and the direct sunlight contribution bilinearly interpolated from the luminance distributions of the four closest precalculated solar positions.

Bourgeois et al. (2008) proposed that 65 solar positions chosen at hourly intervals on five selected days would be sufficient (Figure 8). However, choosing 120 solar positions at hourly intervals on nine selected days (Figure 9) is arguably a better choice. In particular, the average separation between solar positions is eight degrees, and the maximum separation is ten degrees. While a difference of four to five degrees may be significant in calculating the direct sunlight distribution in interior spaces for static scenes (and particularly so for photorealistic renderings), it is likely acceptable for climate-based daylight modeling where the sun traverses 15 degrees of the sky between hourly weather records. (It must also be remembered that each hourly weather record represents the average direct normal and diffuse horizontal irradiances measured over the previous hour.)

FIG. 8 –Annual solar path (65 positions) for Freiberg, Germany. (Source: Bourgeois et al. 2008)

FIG. 8 –Annual solar path (65 positions) for Freiberg, Germany. (Source: Bourgeois et al. 2008)

FIG. 9 – Annual solar path for Boulder, Colorado.

FIG. 9 – Annual solar path for Boulder, Colorado.

An example interior luminance distribution, including both diffuse daylight and direct sunlight, is shown in Figure 10. It should be noted that unlike ray-traced images, the shadow edges are not sharply defined. This is due to the mesh spacing of the floor (in this case 0.2 meters), and it is intentional. (This issue will be addressed in greater detail further on in this article.)

FIG. 10 – Example Perez sky model solution (clear sky).

FIG. 10 – Example Perez sky model solution (clear sky).

Ground Reflections

Direct sunlight and diffuse daylight reflected from the ground must also be taken into account, even if the architectural model does not include exterior surfaces. This can be accomplished by modeling a virtual ground plane as an inverted sky dome (a ground dome) with the same number of ground patches (Figure 11). Each ground patch serves the same purpose as its corresponding sky patch as an area light source. Unlike sky patches, however, all ground patches have the same luminance because they are diffusely reflecting light from the entire sky dome.

Assuming that the average reflectance of outdoor scenes is 18 percent (the same as a photographic gray card), the luminance of all ground patches is equal to 18 percent of the horizontal illuminance. Therefore, only a single interior luminance distribution needs to be determined for exterior ground reflections, which is included with the canonical solution and subsequently scaled according to the horizontal irradiance for a given TMY3 weather record.

FIG. 11 – Sky dome and inverted ground dome.

FIG. 11 – Sky dome and inverted ground dome.

Sky Dome Discretization

Yet another detail: most radiosity methods require all surfaces to be represented as meshes of triangular and quadrilateral elements. This is problematic in terms of the Tregenza sky subdivision, whose sky patches have curved edges. If the sky dome is represented as planar trapezoidal elements (e.g., Figure 12), the resultant gaps will result in errors of several percent or more when calculating luminance distributions due to diffuse daylight.

FIG. 12 – Tregenza sky subdivision errors

FIG. 12 – Tregenza sky subdivision errors.

The solution is to instead represent the sky dome with a hemispherical geodesic dome. The vertex coordinates of each planar sky patch can be determined by recursively subdividing half of an octahedron, as shown in Figure 13. Three subdivisions result in a geodesic dome with 256 triangular (and, of course, planar) patches with no gaps.

FIG. 13 – Octahedron subdivision

FIG. 13 – Octahedron subdivision.

With this, there is another important detail to consider. In their “Daylight Coefficients” paper, Tregenza et al. (1983) assumed that each sky patch would be projected onto a point on each surface element in the environment, as shown in Figure 5. This works in theory, but it is computationally inefficient in the extreme in that all 256 sky patches would need to be projected onto each surface element of both the interior and exterior environments. With today’s architectural models, this could involve hundreds of thousands of elements … and hours to days of computer time (e.g., Müller et al. 1995).

 Parallel Sky Patch Projection

A much more efficient approach is to model each sky patch like the solar disk, as an infinitely distant point source that produces a parallel beam of light (Ashdown 2004). With this, the sky patch illumination can be projected onto the entire environment at once, with all surface elements being considered in parallel (Figure 14). This is a standard computer graphics operation that can be performed either in software (e.g., Ashdown 1994) or in hardware using the computer’s graphics processing unit (Rushmeier et al. 1990).

FIG. 14 – Parallel sky patch projection.

FIG. 14 – Parallel sky patch projection.

This approach works well for most exterior environments because the envelopes of exterior objects (such as buildings) are typically convex. As a result, each exterior surface element is visible to multiple sky patches, with their parallel light beams being averaged and so not noticeable in the computer graphics renderings.

This assumption fails, however, for the windows and openings of interior environments. An example is presented in Figure 15, where the parallel light beams from the discrete sky patches are clearly visible as light “spokes” when they are projected through a narrow window onto the surfaces of interior environment.

Another problem (not illustrated here) is that the light levels in interior environments are often orders of magnitude less than the exterior horizontal illuminance. When the CAD models (such as the simple box shown in Figure 15) are specified, their surface edges may align exactly. However, when these models are rotated, translated, and possibly scaled in world space coordinates for lighting calculations, floating-point round-off of the vertex coordinates may result in very small gaps between surfaces such as the walls and floor.

(This is not a software problem – it is basically impossible to avoid this problem. Even the most high-end computer graphic displays may exhibit occasional single-pixel flickering at the surface edges when the objects are rotated or orbited for viewing.)

Even though the gaps may be fractions of a millimeter wide in world space, they allow direct sunlight or even diffuse daylight to enter the interior environment. While the amount of light is usually insignificant in terms of lighting calculations, the resultant “light leaks” can be quite obvious in computer graphics renderings.

To address these two problems, a radically different approach is needed when handling windows and openings in daylight calculations.

FIG. 15 – Light “spokes” due to parallel beams from discrete sky patches.

FIG. 15 – Light “spokes” due to parallel beams from discrete sky patches.

Windows and Openings

As previously noted, each exterior surface element is visible to multiple sky patches, with their parallel light beams being averaged. Windows are no different – they typically have a full hemispherical view of the exterior environment. Imagine then placing a camera with a fisheye lens facing outwards on the window and capturing a hemispherical image of the visible sky patches and exterior surface elements.

Comparing this to a ray tracing approach, each pixel represents a light ray that is incident upon the camera. It is slightly better than this, as the pixels each “see” a rectangular cone of light with no gaps between them. This being done in software, the camera resolution is arbitrary. With (say) 1.5 million pixels, a highly detailed high dynamic range image can be captured. The value of each pixel is the red-green-blue spectral radiance[3] in its field of view.

If we now reverse the camera orientation such that it faces inwards, we can project this HDR image onto the interior surface elements. This effectively transfers diffuse daylight and direct sunlight through windows and openings without the problems of light spokes, as shown in Figure 16.

FIG. 16 – Virtual cameras on windows eliminates light spokes.

FIG. 16 – Virtual cameras on windows eliminates light spokes.

The virtual camera measures the exterior spectral radiance distribution at its position on the window. At the same time, the projection of the parallel light beams from the sky patches yields the average spectral irradiance of the entire window surface. Knowing this value and the window area, the total amount of spectral radiant flux that is to be transferred through the window can be calculated.

There are, of course, further details – many of them – that need to be considered with this approach. Glass windows and transparent plastic (collectively, dielectric) surfaces exhibit complex reflectance and transmittance characteristics that vary with incidence angle (described by Fresnel equations), and may also exhibit spectrally selective absorptance (e.g., colored glass). There may also be multiple surfaces (e.g., triple-pane glazing) present. All of these issues, however, can be efficiently dealt with using radiosity-based methods.

This approach addresses the light spokes issue, but not light leaks due to floating-point round-off of the environment geometry. Solving this problem requires that the environment be separated into “exterior” and “interior” surfaces. The difference is that interior surfaces are illuminated only by daylight that is transferred through windows and openings (collectively, transition surfaces). With this approach, light leaks can be completely eliminated, no matter how imprecise the environment geometry.

Eliminating Hot Spots

There is yet another problem to consider. If we place the virtual camera at the center of a large window that is too close to an interior surface, it may result in a “hot spot” on the surface (e.g., Figure 17).

FIG. 17 – Single window patch results in “hot spot” on interior surface.

FIG. 17 – Single window patch results in “hot spot” on interior surface.

The reason for this spot is clear: the camera is concentrating all the light received by the window at a single point and projecting it onto the surface. The window is being modeled as a point source, and so the inverse square law applies.

The solution is equally clear, and in fact is required by IES LM-83-12: model large windows as an array of 0.3-meter square window “patches” (IES 2012). As long as the nearest significantly large interior surface is at least 0.6 meters distant (or twice the width of each window patch), the problem of hot spots will be eliminated (e.g., Figure 18).

FIG. 18 – Multiple window patches eliminate hot spots.

FIG. 18 – Multiple window patches eliminate hot spots.

Annual Daylight Simulations

Without delving into the details of how the radiosity method works – see Ashdown (1994) for some 500 pages of explanation – roughly 95 percent of the calculation time involves computing form factors between surface elements for each “step” of the iterative radiosity solution (e.g., Rushmeier et al. 1990). Most of the remaining time per step is spent calculating the amount of light transferred between elements with each “bounce” of light.

With the radiosity method, each surface element is assigned a parameter that represents how much light (technically, spectral radiant exitance) it has received at each step in the calculations. Given that there are 256 sky patches, 120 solar positions, and a virtual ground plane, all that needs to be done (while blithely ignoring the myriad details) is to assign an additional 377 parameters per surface element. These are used to store the 377 separate radiosity solutions that together represent the canonical daylight solution for the environment. This requires a considerable but still manageable amount of memory for even complex enivronments.

The same approach can be applied to electric lighting channels, with one additional parameter per channel. This enables, for instance, the ability to model daylight harvesting systems with switched or dimmable luminaires.

The obvious question is, how quickly can these calculations be performed, particularly for complex environments with tens to hundreds of thousands of surface elements? While these are clearly empirical results, numerous experiments to date have shown that it takes between two and three times as long to calculate a canonical solution as it does to calculate an equivalent static radiosity solution for a specific date and time.

… but the devil still remains in the details …

Transition Surfaces

Looking again at the virtual cameras used to capture images of the exterior environment and project it onto the interior environment, it can be seen that each pixel needs to represent 377 spectral radiance values per pixel. Assuming a 1.5 megapixel image, this represents at least several gigabytes of memory. With a multithreaded program simultaneously processing eight to sixteen cameras, it is clear that the virtual camera approach will fail spectacularly.

The solution to this problem, however, is simple: assign a unique identifier to each surface element, which is then assigned to the pixels that “see” it in the virtual image. When the image is projected onto the interior surface elements, the identifiers can be used to access the exterior surface elements and their canonical solution values.

Looking more closely at windows and openings, Figure 19 diagrammatically shows virtual cameras positioned at the centers of multiple window patches, where each camera captures a hemispherical image of the exterior environment and projects the image into the interior environment.

FIG. 19 – Virtual cameras positioned on multiple window patches.

FIG. 19 – Virtual cameras positioned on multiple window patches.

While this approach clearly works (e.g., Figure 18), it is computationally demanding. With complex architectural models such as buildings with hundreds of windows and potentially thousands of virtual cameras, the calculation times could extend into hours.

With this, it is instructive to consider again the observation that the image is equivalent to tracing a million or more rays through the camera position. With (say) a large window subdivided into one hundred 0.3-meter square window patches, this represents a dense array of 100 million light rays emanating from a single window. This may be necessary because the interior surfaces are relatively close to the window, and so the window patches are needed to avoid the formation of hot spots. For the exterior environment, however, there may not be any similarly close surfaces, and so a single virtual camera positioned in the center of the window is sufficient to capture the hemispherical image (e.g., Figure 20).

FIG. 20 – Single virtual camera positioned on window exterior.

FIG. 20 – Single virtual camera positioned on window exterior.

By itself, this does little to reduce the computational burden of tracing perhaps 100 million rays per window. In practice,however, as few as one million rays are likely sufficient (especially when it is considered that each ray is actually a finite width cone rather than an infinitesimally narrow ray).

The solution to this problem is to first capture the exterior image, and then randomly assign each of its pixels (i.e., rays) to one of the multiple cameras projecting the image into the interior environment. This has the effect of distributing the projected rays across the interior side of the window, which is needed to prevent the formation of hot spots. At the same time, the number of rays that need to be traced is reduced to those of a single virtual camera.

This process can be made adaptive, because each image pixel also provides the distance from the camera to the first intersected surface. By first capturing and projecting an image from the center of the window, this depth information can be used to identify the closest surfaces and so decide how best to subdivide the window on both sides (exterior and interior) for camera placement.

As complicated as this may look, it is computationally efficient. As an example, the environment shown in Figure 10 required 28 seconds of calculation time on a commodity desktop computer when calculated as a static environment for a given time and date. It is a simple environment with only 800 surface elements, but most of the calculation time is spent on transferring the direct sunlight and diffuse daylight through 84 window patches.

By comparison, the canonical solution shown in Figure 6 required only 42 seconds of calculation time. Once completed, any of the 4,380 hourly weather records in the TMY3 weather data can be calculated and displayed in milliseconds.

Bidirectional Transmittance Distribution

One of the criticisms of the radiosity approach for daylight simulation is that it can only model diffuse reflections. Quoting IES RP-5-13, Recommended Practice for Daylighting Buildings (IES 2013):

Radiosity methods assume that all surface materials have perfectly diffuse reflectance, i.e., they reflect light equally in all outgoing directions for all incident angles of incoming radiations.

This was a true statement when radiosity methods were first developed in the late 1980s, but it is certainly not true today. Radiosity methods are capable of accurately modeling the optical and spectral properties of opaque, transmissive, and translucent surfaces, including Fresnel reflectance and transmittance. They are also capable of accurately modeling both isotropic and anisotropic bidirectional reflectance and transmittance distribution functions (BRDFs and BTDFs) using analytic functions or measured data represented by the LBNL bidirectional scattering distributions function (BSDF) data format (once its specification has been finalized and published).

If commercial lighting design software has yet to support all of this functionality, it is only because there has been insufficient demand to date. IES LM-83-12 (IES 2012), however, specifies the modeling of window shades and blinds using their BSDF properties if available. Once this information become widely available from manufacturers, its simulation for CBDM purposes can be supported.

There is, however, no need to wait for measured BSDF data to become available in order to comply with the requirements of LM-83-12. Apart from advanced fenestration devices that redirect light, most diffusing glazings and fabric window shades behave as simple diffusers. Given a virtual camera, it is straightforward to include a virtual diffusion filter to model these glazings and shades.

Figure 21 illustrates several example window bidirectional transmittance distributions at incidence angles ranging from 20 to 60 degrees, and with three different degrees of diffusion. These would normally be applied to the direct sunlight entering an interior environment, but the same approach can be applied to light redirection devices such as semispecular light shelves.

FIG. 21 – Example window bidirectional transmittance distributions.

FIG. 21 – Example window bidirectional transmittance distributions.

It is equally straightforward – in principle – to define a BSDF “filter” in software for the virtual camera that will redirect incident rays such that they are either transmitted or reflected by the window or other designated surface. This, however, is a work in progress for the radiosity approach.

Virtual Photometers

There is no point in calculating the distribution of daylight within an interior environment if it cannot be measured for daylight metric calculations and other purposes. Perhaps surprisingly, this is where radiosity-based CBDM has a distinct advantage over ray tracing methods.

The original radiosity methods only provided illuminance and luminance measurements for selected points on opaque surfaces after the radiosity calculations had been completed. However, it possible to extend these methods such that transparent surfaces can be specified that receive but do not otherwise influence the flow of light within the environment. By subdividing these surfaces, each surface element becomes a virtual photometer (e.g., Figure 22).

FIG. 22 – Virtual photometers.

FIG. 22 – Virtual photometers.

The advantage of these photometers is that they record all the light passing through them in one direction – which is exactly what is needed when performing most daylight metric calculations. By comparison, an array of photometers used in ray-traced CBDM calculations will only measure the light incident at each meter position. Depending on the resolution of the direct sunlight shadows cast by Venetian blinds, for example, those meters that are shadowed may result in the spatial Daylight Availability (sDA) and Annual Sunlight Exposure (ASE) values being underreported.

Another advantage of these photometers is that they are computationally efficient. In the above example, there are 800 surface elements and 1,500 virtual photometers placed on an imaginary workplane. Where it took 42 seconds to calculate the canonical solution without the photometers, it took 57 seconds with them. Again, this is for the canonical solution – calculating the meter values for a given TMY3 weather record is a matter of a few more milliseconds.

Looking Forward

Much of the above has been developed specifically for climate-based daylight modeling using radiosity methods. It has been implemented in commercial lighting design software, but there will undoubtedly be minor changes as users gain experience with its features and capabilities.

There are also opportunities for further improvements and optimizations, but these will not be discussed until the research and development work has been completed – hopefully in less time than the twelve years it has taken to reach this point!

To date, climate-based daylight modeling has relied on ray tracing methods, specifically Radiance and its derivatives. There is nothing wrong with this, other than that it has possibly hindered research into alternative approaches.

This, however, is symptomatic of a larger issue. With software products such as Lightscape in the early 1990s, radiosity methods were favored for Hollywood’s computer graphics requirements. However, the introduction of photon mapping techniques (Jensen 2001) and vastly increased computing power through “rendering farms” (thousands of dedicated computers on a network), followed by techniques such as multidimensional light cuts and Metropolis light transport, soon eclipsed radiosity methods.

Research into radiosity methods peaked in 1994 and has been declining ever since, as evidenced by the yearly total of academic papers on the topic shown in Figure 23. It is, in the terminology of computer scientists, a “solved problem.”

Fig. 23 – Radiosity papers publication frequency.

Fig. 23 – Radiosity papers publication frequency.

As this article has attempted to show, however, there is still life to be had in radiosity methods. As the humorist Mark Twain famously never said, “The reports of my death have been greatly exaggerated.”

Acknowledgements

Thanks to Dawn De Grazio of Lighting Analysts Inc. for her review and comments on this article.

References

  1. Ashdown, I. 1994. Radiosity: A Programmer’s Perspective. New York, NY: John Wiley & Sons. Freely available from helios32.com/resources.htm.
  2. Ashdown, I. 2004. “Modeling Daylight for Interior Environments,” 2004 IESANZ Annual Conference Proceedings, Illuminating Engineering Society of Australia and New Zealand.
  3. Bourgeois, D., C. F. Reinhart, and G. Ward. 2008. “Standard Daylight Coefficient Model for Dynamic Daylighting Simulations,” Building Research & Information 36(1):68-82.
  4. 2003. Spatial Distribution of Daylight – CIE Standard General Sky, CIE Standard S 011/E:2003. Vienna, Austria: CIE Central Bureau.
  5. 2010. IES Lighting Handbook, Tenth Edition. New York, NY: Illuminating Engineering Society.
  6. 2012. LM-83-12, IES Spatial Daylight Autonomy (sDA) and Annual Sunlight Exposure (ASE). New York, NY: Illuminating Engineering Society.
  7. 2013. RP-5-13, Recommended Practice for Daylighting Buildings. New York, NY: Illuminating Engineering Society.
  8. Jensen, H. W. 2001. Realistic Image Synthesis Using Photon Mapping. Natick, MA: A. K. Peters.
  9. Mardaljevic, J. 2006. “Examples of Climate-Based Daylight Modeling,” Paper No. 67, CIBSE National Conference 2006: Engineering the Future.
  10. Müller, S., W. Kresse, N. Gatenby, and F. Schöffel. 1995. “A Radiosity Approach for the Simulation of Daylight,” Rendering Techniques ’95 (Proceedings of the Sixth Eurographics Workshop on Rendering), pp. 137-146. New York, NY: Springer-Verlag.
  11. Noorian, A. M, I. Moradi, and G. Kamali. 2008. “Evaluation of 12 Models to Estimate Hourly Diffuse Irradiation on Inclined Surfaces,” Renewable Energy 33:1406-1412.
  12. Muneer, T. 2004. Solar Radiation and Daylight Models, Second Edition. Oxford, UK: Elsevier Butterworth-Heinemann.
  13. Nabil, A., and J. Mardaljevic. 2006. “Useful Daylight Illuminances: A Replacement for Daylight Factors,” Energy and Buildings 38(7):905-913.
  14. Perez, R., R. Seals, and J. Michalsky, 1993. “All-Weather Model for Sky Luminance Distribution – Preliminary Configuration and Validation,” Solar Energy 50(3):235-245.
  15. Rushmeier, H. E., D. R. Baum, and D. E. Hall. 1990. “Accelerating the Hemi-Cube Algorithm for Calculating Radiation Form Factors,” ASME Journal of Heat Transfer 113:1044-1047.
  16. Tregenza, P. R., and I. M. Waters. 1983. “Daylight Coefficients,” Lighting Research and Technology 15(2):65-71.
  17. Tregenza, P. R. 1987. “Subdivision of the Sky Hemisphere for Luminance Measurements,” Lighting Research and Technology 19:13-14.
  18. 2013. LEED v4 BD+C: Schools – Daylight. Washington, DC: U.S. Green Building Council (www.usgbc.org).
  19. Wienold, J., and J. Christoffersen. 2006. “Evaluation Methods and Development of a New Glare Prediction Model for Daylight Environments with the Use of CCD Cameras,” Energy and Buildings 38:743–757.
  20. Wilcox, S., and W. Marion. 2008. Users Manual for TMY3 Data Sets. Technical Report NREL/TP-581-43156, Revised May 2008.Golden, CO: National Renewable Energy Laboratory.

[1] The World Meteorological Organization (www.wmo.int) defines “climate” as the “average weather” over a period of thirty years. The weather at a given location will of course vary on a per-year basis – sometimes drastically – from that of the Typical Meteorological Year weather data for that location.

[2] A uniform sky luminance distribution is equivalent to CIE Standard General Sky Type 5 (CIE 2003).

[3] Spectral radiance in this context refers to representing visible light as a combination of red, green, and blue (RGB) light. Assuming a 6500K light source, the equivalent luminance value L is given by L = 0.2125 * R + 0.7154 * G + 0.0721 * B.

 

 

 

 

Solar Illumination

The Inconstancy of Constants

Ian Ashdown, FIES

Chief Scientist, Lighting Analysts Inc.

January 21, 2015

[ Please send comments to allthingslighting@gmail.com ]

Lighting design is based in part on the reasonable assumption that photometric units have precise definitions. The candela, for example, has a precise mathematical definition:

The candela is the luminous intensity, in a given direction, of a source that emits monochromatic radiation of frequency 540×1012 hertz and that has a radiant intensity in that direction of 1/683 watt per steradian.

and from this all other photometric units are derived, including luminous flux, luminous intensity, luminance, illuminance, and so forth.

You might expect then that the same applies for daylighting design … but you would be wrong.

Solar Constant

We begin with the solar constant, which is defined as the solar irradiance (measured in radiant watts per square meter) incident on a plane perpendicular to the Sun at a distance of one astronomical unit (AU), which is roughly the mean distance from the Sun to the Earth.

Defining the solar constant is easy; measuring it is not. Although it was first roughly measured in 1838, accurate measurements must be performed using a satellite above the Earth’s atmosphere. You then have to correct for the distance of the Earth from the Sun, which varies due to the ellipticity of the Earth’s orbit. (The actual direct solar irradiance at the top of the atmosphere varies by some 6.9% over the course of a year.)

The technology of precision spectroradiometers has of course improved over the years, and with it increasingly precise measurements of the solar spectrum and the solar constant. By 1993, satellite measurements had a precision of 0.01 percent (Kittler and Darula 1996).

Today, the extraterrestrial solar spectrum is formally defined by ASTM International Standard E490-00a (2014), “Standard Solar Constant and Zero Air Mass Solar Spectral Irradiance Tables.” This standard defines the solar spectrum (FIG. 1) over a range of wavelengths from 200 nanometers (far ultraviolet) to one millimeter (far infrared).

Solar Illumination - FIG 1FIG. 1 – Solar spectrum (Source: ASTM 2014)

Integrated over the entire solar spectrum, the solar constant is 1366.1 ± 0.58 watts per square meter. This is of course something of a fiction, as the total radiation output of the Sun varies by approximately 0.1 percent due to the 11-year sunspot cycle and other longer-term effects. In general however, the output has varied by less than 0.2 percent over the past two millennia.

Solar Illumination Constant

The solar constant is useful for things like satellite design and climate studies. For lighting design purposes, however, only visible light in the range of 380 nm to 780 nm is important. By multiplying each wavelength by the CIE 1931 luminous efficiency function V(l), according to:

Solar Illumination - EQN 1

the solar illumination constant can be calculated.

Solar Illumination - FIG 2FIG. 2 – CIE 1931 luminous efficiency function V(l)

Based on ASTM (2014), the calculated solar illumination constant Esc is 133.1 kilolux (IES 2010).

There is nothing controversial here … right?

Daylighting Design

The lighting design and simulation program Radiance has long been a standard tool within the daylighting research community. The name “program” is a bit of a misnomer, as Radiance is a set of hundreds of Unix-based utility programs that are executed via command-line scripts. One of these programs is gendaylit. This is an essential tool for daylighting research, as it calculates the absolute sky luminance distribution in accordance with the Perez sky model (Perez et al. 1993).

The issue is that gendaylit was written in 1994. At that time, the standard solar spectrum was Wehrli (1985), which resulted in a solar constant of 1367 watts per square meter and a solar illumination constant of 127.5 kilolux. These values were hard-coded as constants in gendaylit.

When the IES Lighting Handbook, Tenth Edition, was being prepared, the solar illumination constant was recalculated based on ASTM E-490-00a(2006). This document was first published in 2000, and has since been reapproved and published as ASTM E490-00a(2014).

So herein lies the problem: which solar illumination constant should be used? To be compatible with the IES Lighting Handbook, Tenth Edition, it should be 133.1 kilolux. However, to be compatible with the “gold standard” Radiance used by the daylighting research community, it should be 127.5 kilolux.

This is not exactly an academic question, as the difference is 4.4 percent. This difference permeates all absolute daylighting values calculated using the Perez sky model. For annual daylight metric calculations such as spatial Daylight Autonomy (sDA) and Annual Sunlight Exposure (ASE), this difference can be significant (IES 2012).

It is an open question.

References

ASTM. 2014. Standard Solar Constant and Zero Air Mass Solar Spectral Irradiance Tables, ASTM E490-00a(2014). West Conshohocken, PA : ASTM International.

IES. 2010. IES Lighting Handbook, Tenth Edition. New York, NY: Illuminating Engineering Society.

IES. 2013. Approved Method: IES Spatial Daylight Autonomy (sDA) and Annual Sunlight Exposure (ASE). New York, NY: Illuminating Engineering Society.

Kittler, R., and S. Darula. Standardization of Daylight Phenomena and Algorithms: Luminous Solar Constant and its Variations. Bratislava, Slovakia: Institute of Construction and Architecture, Slovak Academy of Sciences.

Perez, R., R. Seals, and J. Michalsky, 1993. “All-Weather Model for Sky Luminance Distribution – Preliminary Configuration and Validation,” Solar Energy 50(3):235-245

Wehrli, C. 1985. Extraterrestrial Solar Spectrum, Publication No. 615. Davos, Switzerland: World Radiation Center.

 

Daylight Factors

Ian Ashdown, P. Eng., FIES

Chief Scientist, Lighting Analysts Inc.

[ Please send comments to allthingslighting@gmail.com. ]

Nine out of ten daylight simulation programs agree … and therein lies a story worth retelling.

Daylight in History

The story begins in the sixth century with the publication of Corpus Juris Civilis (“Body of Civil Law”) by order of the eastern Roman emperor Justinian I [Scott 1932]. Written in four volumes, it included the Digest, being extracts from the writings of earlier Roman jurists. Book VIII, Title 2, “Concerning Servitudes of Urban Estates,” includes this legal distinction between daylight and views:

Light is the power of seeing the sky, and a difference exists between light and view; for a view of lower places may be had, but light cannot be obtained from a place which is lower.

In a time when artificial lighting consisted of oil lamps, access to daylight was a critical issue. The Digest therefore had some 40 legal rulings on the rights of property owners concerning daylight. In some cases, a property owner whose newly-constructed building blocked a neighbor’s access to daylight could be legally compelled to tear the building down.

Somewhat surprisingly, these rulings survived over the centuries to become what is referred to as the “ancient lights” law in European legal traditions. Modern use of this concept dates back to the British Prescription Act of 1832, which reads in part:

When the access and use of light to and for any dwelling house, workshop, or other building shall have been actually enjoyed therewith for the full period of twenty years without interruption, the right thereto shall be deemed absolute and indefeasible, any local usage or custom to the contrary notwithstanding, unless it shall appear that the same was enjoyed by some consent or agreement expressly made or given for that purpose by deed or writing.

This enactment led of course to a new profession: Rights to Light surveyors. These chartered professionals served as expert witnesses in legal disputes, and offered advice to architects. The surveyor Robert Kerr wrote a book on the topic in 1865, in which he took 55 pages to explain the practice of surveying access to daylight [Kerr 1865].

Sky Factors

Complexity begets uncertainty however, and so it was that a chartered surveyor and lighting engineer named Percy J. Waldram proposed a much simpler way of determining adequate access to daylight [Waldram 1909]. It is today known in the United Kingdom as a “sky factor,” and is defined as:

Sky factor: the ratio of the illuminance Eindoor of a horizontal plane at a given point inside a building due to the light received directly from an overcast sky of uniform luminance sky, to the illuminance Eoutdoor of the point due to the unobstructed sky.

DF = ( Eindoor / Eoutdoor ) * 100 %

FIG 1

Waldram and his son later suggested that a sky factor of 0.2 percent was sufficient “for ordinary purposes, comparable to clerical work” [Waldram and Waldram 1923]. Their seemingly offhand comment that this was the level at which “average reasonable persons would consistently grumble” became what is now referred to by British light surveyors in all seriousness as the “grumble point” [Chynoweth 2004].

The sky factor metric was accepted by the Commission Internationale de L’Eclairage in 1932 [CIE 1932], and is still widely used in the United Kingdom by chartered surveyors. (It is however no longer recognized as a unit of light measurement by the British Standards Institution.)

Waldram’s choice of 0.2 percent was simply a rule-of-thumb guess [Waldram and Waldram 1923], with no supporting research [Chynoweth 2005]. The Royal Institution of Chartered Surveyors recommends at least 0.5 percent (which on average is about 25 lux) [RICS 2010], but the British legal system still works on the assumption of 0.2 percent, or 10 lux [Chynoweth 2009]. This leads to the curious situation of a complex legal system that deliberately encourages and enforces poor daylighting practices.

Daylight Factors

The sky factor metric is important in the United Kingdom because “rights to light” is a still-valid legal concept as an “easement right” that has descended from the rulings of ancient Roman jurists. This is not however the situation in American civil law, where access to daylight is considered a right only in exceptional circumstances [Unger 2005].

Regardless, the closely related “daylight factor” metric is widely used [CIE 1970]. This has the same basic definition as the sky factor above, with the exception that the CIE Standard Overcast Sky is used instead of a uniform luminance sky [IES 2013]. It also takes into account ground reflections, window transmittance, and interreflections from room surfaces.

This distinction is important. Both sky types are defined in CIE Standard S 011 [CIE 2003], where the CIE Standard Overcast Sky (Standard Sky Type 1) is described as, “steep luminance gradation towards zenith, azimuthal uniformity.” The sky factor metric however assumes CIE Standard Sky Type 5, “sky of uniform luminance.” (The “traditional” CIE Standard General Sky is now referred to as CIE Standard Sky Type 16; its luminance distribution near the horizon varies slightly from Sky Type 1.)

The advantage of a uniform luminance sky is that determining the sky factor for a given room is simply a matter of geometry. That is, the value of the sky factor is completely independent of the sky luminance distribution. It can in theory be calculated with the aid of Waldram diagrams using photographs or hand drawings [e.g., RICS 2010].

Waldram Diagram

The same is not true of course for the calculation of daylight factors. The luminance distribution of the CIE Standard Overcast Sky varies continuously from the horizon to zenith, and so computer calculations with a 3D CAD model of the room are essential. (Accurate calculation of interreflections from room surfaces also requires computer calculations.)

Limitations

As a lighting design tool, the daylight factor metric has numerous limitations. As noted by IES RP-5-13, Recommended Practices for Daylighting Buildings [IES 2013], it addresses only a single sky condition that is prevalent in its country of origin. It does not consider such daylighting design issues as direct and reflected sunlight, latitude, building orientation, time and date, or climatic conditions.

These are serious limitations in that satisfying a daylight factor requirement may result in excess daylight under clear sky conditions. IES RP-5-13 therefore recommends modern climate-based annual daylight performance metrics such as spatial Daylight Autonomy (sDA) and Annual Sunlight Exposure (ASE) in combination with Typical Meteorological Year (TMY) weather data files specific to a given geographical location [IES 2012].

As an aside, it is a myth that the sky factor metric was developed specifically for northern European climates. Prior to the widespread availability of electric lighting, access to diffuse daylight was preferred to that of direct sunlight. (The studios of artists and photographers for instance had north-facing windows and skylight wherever possible.) As the Victorian-era Keller wrote [Keller 1865]:

In fact, it is this diffused daylight which is constituted, by the express intent of Nature, the standard medium of human vision; for where there is one purpose of sight specially served by the direct and unobscured light of the sun, there must be a thousand for which the eye prefers the more genial agency of the diffused light of the atmosphere.

(Is it any wonder that Keller needed 55 pages to explain the practice of rights to light surveying?)

Even by the time of Waldram a half-century later, instruments for measuring daylight were not commonly available. Assuming a uniform luminance sky that allowed purely geometric calculations was therefore a matter of practical necessity.

Limitations aside, reports of the death of the daylight factor metric have been greatly exaggerated – it is still a useful if primitive tool, especially for students learning the basics of daylighting design. Climate-based annual daylight performance metrics such as sDA and ASE may be recommended for detailed analysis of architectural designs, but the daylight factor metric provides the necessary sanity checks.

Nine Out Of Ten

As noted previously, calculating daylight factors requires a 3D CAD model and daylight simulation software, if only to accurately model the spatial luminance distribution of the CIE Standard Overcast Sky and room surface interreflections. The obvious question is, how accurate are daylight simulation programs?

This question was investigated some sixteen years ago, and the results were not encouraging. Point-by-point errors for clear sky conditions were as much as 18 times, and for overcast sky conditions as much as 10 times the measured values. Much worse, the average error for clear sky conditions was about ten times for two of the four programs investigated.

Today, the situation is markedly different. A recent study [Iverson et al. 2013] investigated the ability of nine daylight simulation programs to calculate the daylight factor metric in five typical rooms. Going well beyond the basic requirements of CIE 171, “Test Cases to Assess the Accuracy of Lighting Computer Programs” [CIE 2006], this study is exemplary of how the accuracy of lighting design and simulation software should be assessed.

The study included nine daylight simulation programs:

  • Radiance
  • Daysim
  • Desktop Radiance
  • IESve
  • DIALux
  • Relux
  • Ecotect
  • VELUX
  • LightCalc

Given that these programs use a wide variety of radiosity and ray tracing techniques, you might hope to see at least reasonable agreement among their daylight factor predictions.

What the study revealed however was stunning: all but one of the programs agreed to within a few percent of each other.

Lighting Analysts’ AGi32 was not included in the study, but the lead author of the report kindly provided the CAD files for the test rooms so that Lighting Analysts could perform its own tests. The results were the same: agreement to within a few percent.

Test Room Example

Nine out of ten lighting programs agree … this is important news for lighting designers and architects involved in daylighting design.

The test models used in the study took into consideration room dimensions, surface reflectances, glass transmittance, and exterior obstructions. What this study in effect says is that whatever daylight simulation program (with one important exception – and it was not LightCalc) is chosen, daylight factor calculations will be within the ±10 percent accuracy range expected for such programs [Reinhart and Andersen 2006].

The full “Daylight Calculations in Practice” study is available here, and the Lighting Analysts follow-up study is available here.

References

Chynoweth, P. 2004. “Progressing the Rights to Light Debate – Part 1: A Review of Current Practice,” Structural Survey 22(3):131-137.

Chynoweth, P. 2005. “Progressing the Rights to Light Debate – Part 2: The Grumble Point Revisited,” Structural Survey 23(4):251-264.

Chynoweth, P. 2009. “Progressing the Rights to Light Debate – Part 3: Judicial Attitudes to Current Practice,” Structural Survey 27(1):7-19.

CIE. 1932. Recueil des Travaux et Compte Rendu des Seances, Huiteme Session (Collection of Works and Minutes of Sittings, Eighth Session). Cambridge, UK: Cambridge University Press.

CIE. 1970. Daylight. CIE 016-1970. Vienna, Austria: CIE Central Bureau.

CIE. 2003. Spatial Distribution of Daylight – CIE Standard General Sky. CIE S011/E:2003. Vienna, Austria: CIE Central Bureau.

CIE. 2006. Test Cases to Assess the Accuracy of Lighting Computer Programs. CIE 171:2006. Vienna, Austria: CIE Central Bureau.

IES. 2012. IES Spatial Daylight Autonomy (sDA) and Annual Sunlight Exposure (ASE). IES LM-83-12. New York, NY: Illuminating Engineering Society of North America.

IES. 2013. Recommended Practice for Daylighting Buildings, IES RP-5-13. New York, NY: Illuminating Engineering Society of North America.

Iversen, A., N. Roy, M. Hvass, M. Jørgensen, J. Christoffersen, W. Osterhaus, and K. Johnsen. 2013. Daylight Calculations in Practice. Technical Report SBi 2013:26. Copenhagen, Denmark: Danish Building Research Institute, Aalborg University.

Kerr, R. M. 1865. On Ancient Lights and the Evidence of Surveyors Thereon. (Available as https://archive.org/details/onancientlights00kerrgoog.)

Reinhart, C. F., and M. Andersen. 2006 “Development and Validation of a Radiance Model for a Translucent Panel,” Energy and Buildings 38(7):890-904.

RICS. 2010. Rights of Light – Practical Guidance for Chartered Surveyors in England and Wales, First Edition. Conventry, UK: Royal Institution of Chartered Surveyors. (Available as www.msasurvey.com/map/RICS Rights of light guidance note 2010.pdf.)

Scott, S. P. The Civil Law, Including the Twelve Tables, the Institutes of Gaius, the Rules of Ulpian, the Opinions of Paulus, the Enactments of Justinian, and the Constitutions of Leo. English Translation. Cincinnati, OH: The Central Trust Company. (Available from http://droitromain.upmf-grenoble.fr/Anglica/digest_Scott.htm.)

Ubbeholde, M. S., and C. Humann. 1998. “A Comparative Evaluation of Daylighting Software: Superlite, Lumen Micro, Lightscape and Radiance,” Proceedings of the International Daylighting Conference (Daylighting ’98), pp. 97-104.

Unger, S. C. 2005. “Ancient Lights in Wrigleyville: An Argument for the Unobstructed View of a National Pastime,” Indiana Law Review Vol. 38, pp. 533-564.

Waldram, P. J. 1909. “The Measurement of Illumination, Daylight and Artificial, with Special Reference to Ancient Light Disputes,” Journal of the Society of Architects 3:131-140.

Waldram, P. J., and J. M. Waldram. 1923. “Window Design and the Measurement and Predetermination of Daylight Illumination,” The Illuminating Engineer, Vol. XVI, pp. 90-122.

 

Thoughts on Color Rendering

Ian Ashdown, P. Eng., FIES

Chief Scientist, Lighting Analysts Inc.

[ Please send comments to allthingslighting@gmail.com. ]

UPDATE 14/10/06 – LightingEurope, the “Voice of the Lighting Industry,” has just published their LightingEurope Position Paper on Color Quality. To summarize:

1. LightingEurope supports to continue the use of the existing Color Fidelity metric CRI including eight reference colors.

2. LightingEurope supports to keep legal minimum requirements on CRI on the current level as defined in the EU Eco-design Regulation.


Lamp with a CRI of 90 or above, good. Lamp with a CRI of less than 80, bad. Is there anything else that lighting designers need to know about the Color Rendering Index (CRI) metric?

To be brutally honest … no. Despite all that has been written on the topic over the past decade, the importance of CRI to everyday lighting design today is minimal at best.

What is (or perhaps was) important is the history of color rendering metrics and the influence they had on fluorescent and LED lamp design. We can mostly ignore the issues of color rendering today precisely because of the CRI metric.

If you want to make senior citizens shudder, ask them what it was like to work in an office in the 1950s with fluorescent lighting. The linear fluorescent lamps of the time used calcium halophosphate phosphors that had a well-deserved reputation for making skin tones appear a sickly gray-green. Brightly-colored fabrics also looked disagreeably different from being viewed under daylight or incandescent lighting conditions.

Figure 1

FIG. 1 – Halophosphate lamp spectrum

Fluorescent lamp manufacturers could address this problem by varying the phosphor composition of their lamps. However, phosphors are expensive, and it is always difficult to convince consumers to pay more for products on the promise that they will “look better.” What was needed was an industry-standard metric.

Beginning in 1948, the Commission Internationale de l’Eclairage (CIE) began the quarter-century process of developing what is now the CIE General Colour Rendering Index, commonly referred to as the CRI metric (CIE 1995). The first version was published in 1965, and it was revised in 1974 to include the psychophysiological effects of chromatic adaptation.

An excellent description of the metric is available from Wikipedia (“color rendering index”), so there is no reason to repeat it here. All that needs to be said is its definition:

Color rendering: effect of an illuminant on the color appearance of objects by conscious or subconscious comparison with their color appearance under a reference illuminant.

and a reminder that the two illuminants (i.e., light sources) must have the same correlated color temperature (CCT).

This metric worked reasonably well for ranking linear fluorescent lamps from the era. Quartz halogen lamps had CRIs of nearly 100, while warm white fluorescent lamps typically had CRIs between 50 and 60. These lamps were particularly deficient in the red region of the spectrum (see Fig. 1), with warm white lamps having CRI R9 values as low as -111. (No, that is not a misprint; CRI values for specific test colors can be negative.)

Fluorescent lamp manufacturers could increase the red emission by mixing strontium and calcium halophosphates to create so-called “deluxe” phosphors. Lamps using these phosphors could achieve CRIs of approximately 90, but at the cost of roughly a one-third decrease in luminous efficacy (lumens output per electrical watt input).

This situation changed in the 1970s with two important discoveries:

  1. Lamps with improved luminous efficacy and very good color rendering properties could theoretically be achieved with three narrow-band (red, green and blue) lamp spectra (Thornton 1971); and
  2. The development of rare-earth phosphors for Thornton’s “triphosphor” lamps with CRIs of approximately 85 (Verstegen et al. 1974).

Figure 2

FIG. 2 – Rare-earth triphosphor lamp spectrum

By themselves, these two discoveries may or may not have had a significant impact on the manufacture of fluorescent lamps. Looking back, the lighting industry at the time had little to no interest in color rendering issues.

What the lamp manufacturers did have however was a metric to compare products with, and with this the opportunity for an effective marketing campaign. As a result, the otherwise-obscure CRI metric appeared in every manufacturer’s catalogs and sales literature, and it sold lamps.

With CRIs in the range of 85, triphosphor lamps have very good but not excellent color rendering properties. Lamp manufacturers therefore developed so-called “broadband” phosphors with four or five emission bands. With these, CRIs of 90 or so could be achieved.

Figure 3

FIG. 3 – Rare-earth broadband lamp spectrum

Marketing aside, did this really matter? For most commercial applications, the answer was probably no. In 1986, the CIE Guide on Indoor Lighting (CIE 1986) offered this helpful table:

CRI Ra Examples of Usage
> 90 Color matching, art galleries
80 – 90 Homes, restaurants, textile industry
60 – 80 Offices, schools, light industry
40 – 60 Heavy industry
20 – 40 Outdoors

Table 1 – CRI examples of usage

True, this list was likely influenced by the availability of halophosphate fluorescent lamps for indoor use, clear mercury vapor HID lamps for high-bay factory luminaires, and low-pressure sodium (LPS) lamps for roadway and area lighting. Still, it indicated how the CIE viewed its own metric at the time. CRI values were meant to be used as design guides rather than as precise numbers.

More tellingly, van Trigt (1999) presented a scholarly review of the CRI metric, in which he stated that “only a difference of some five points in the index is considered meaningful.”

The problem is that while the CRI metric Ra can be calculated from the measured lamp spectral power distribution (SPD) with arbitrary precision, it is nevertheless based on a mathematical model (the von Kries transformation) of psychophysiological behavior. Given this, it makes sense that the difference between CRI values of for example 88 and 90 is essentially meaningless.

In a sense, the CRI metric has served its purpose in promoting the development and commercialization of rare-earth lamp phosphors. With fluorescent lamp CRIs typically being in the range of 85 to 95 these days, lighting designers and consumers have little need to know anything other than “90 and above good, less than 80 bad.”

But then along came solid-state lighting …

The first commercially-produced SSL luminaire designed expressly for architectural applications was the TIR Systems Lexel, introduced at LightFair in April 2005 (Whitaker 2005). Based on a red-green-blue, high-flux LED die design, it generated white light whose color temperature could be varied from 3000 to 6500 kelvins.

Compared to the cool white LEDs with YAG phosphors and CRIs of approximately 75 available at the time, the white light produced by the Lexel was widely acclaimed by the trade show attendees, particularly for its color rendering properties.

Figure 4

FIG. 4 – LightFair 2005 – TIR Systems Lexel™

What the attendees did not know was that (ahem) the CRI Ra value varied from 25 at 3000K to 40 at 6500K (Speier and Salsbury 2006). By the standards of CIE 29.2, this was barely good enough for outdoor lighting only, along with clear MV and LPS lamps.

This seems odd, particularly when you compare the RGB LED lamp spectrum (FIG. 5) with that of a typical rare-earth triphosphor fluorescent lamp spectrum (FIG. 2). Based on this, you might guess that the CRI should be 85 to 90, not 25 to 40.

Figure 5

FIG. 5 – RGB LED lamp module spectrum

The answer lies in the dominant wavelengths of the emission peaks. Thornton (1974) calculated that the ideal dominant wavelengths (what he called “prime colors”) for triphosphor lamps were 450 nm (blue), 545 nm (green), and 610 nm (red). The dominant wavelengths of the Lexel – and indeed most color-changing RGB LED luminaires on the market today – were 465 nm, 525 nm, and 615 nm. If you were to change the green wavelength from 525 nm to 545 nm, the CRI would be 85 or so.

Unfortunately, the maximum dominant wavelength of reasonably efficient green InGaN LEDs is approximately 530 nm. The white light produced by color-changing RGB LED luminaires looks wonderful, but it is doomed have unacceptably low CRI values.

If anything, this is an example of the abject failure of the CRI metric to predict the color rendering properties of RGB LED luminaires. CIE Technical Committee 1-62 acknowledged this problem (CIE 2007), and recommended the development of a new and improved color rendering metric for all white light sources.

CIE Technical Committee 1-69 was therefore established in 2008 to “investigate new methods for assessing the colour rendition properties of white-light sources used for illumination, including solid-state light sources, with the goal of recommending new assessment procedures.” The committee investigated over a dozen proposals, but could only agree to bitterly disagree on any new metrics. As of this writing, the committee has yet to release its final report.

In response, the CIE recently established two new committees to further study the issue:

TC 1-90: Colour Fidelity Index. To evaluate available indices based on colour fidelity for assessing the colour quality of white- light sources with a goal of recommending a single colour fidelity index for industrial use.

TC 1-91: New Methods for Evaluating the Colour Quality of White-Light Sources. To evaluate available new methods for evaluating the colour quality of white-light sources with a goal of recommending methods for industrial use.

with reports due no earlier than 2015.

Despite numerous calls from the lighting industry for some clarity on color rendering metrics (e.g., Whitaker 2010, Colombo 2013), it is unlikely that the CIE will respond for at least a few more years.

Again however, does any of this really matter? The solid-state lighting industry effectively gave up waiting a number of years ago and began using CIE Special Colour Rendering Index R9 in addition to Ra to quantify the color rendering properties of white light sources (including both semiconductor and organic LEDs) for saturated red colors. While this combination of Ra and R9 is not a perfect solution, it is nonetheless a recognized industry standard, and it generally works (especially for marketing literature).

More to the point however is that the solid-state lighting industry has, like the fluorescent lamp industry before it, mostly outgrown the need for color rendering metrics. As long as the lamp module CRI is 80 or above – which is the case for most commercial products these days – there is little need to worry about CRI except for applications requiring critical color judgment.

With this, it is interesting to look at another industry that relies on the CRI metric: architectural glass. If you think about it, daylight illuminating interior spaces is spectrally filtered by glass windows and curtain walls. What is the CRI of daylight inside the building? If it is less than 80, it is in danger of being banned altogether by the US Environmental Protection Agency in accordance with the minimum CRI requirements of its ENERGY STAR program (EPA 2013).  (Just kidding … I think.)

It may surprise lighting designers to know that there is a European standard (BSI 2011) that specifies the calculation of indoor daylight CRI values, assuming a 6500 K daylight (CIE D65) illuminant. Architectural glass manufacturers publish CRI values for different thicknesses of their glass and window assemblies, and it can be calculated using the International Glazing Database (windowoptics.lbl.gov/data/igdb) and the freeware Optics 6 program from LBNL (windows.lbl.gov/software/Optics/optics.html).

You might think that bronze glass for example would significantly affect the color of indoor daylight, but this is not the case. Taking one manufacturer (Pilkington Glass) as an example, most of their products have CRIs in the mid-90s, with the lowest (Ra of 77) being for their Solar-E Arctic Blue low-emissivity glass with a blue body tint.

More interesting perhaps is that most architectural glass products have similar spectral transmittance spectra (Gombos et al. 2008). They are so similar in fact that the CIE has defined two “indoor daylight illuminants” (ID50 and ID65) that correspond to CIE daylight illuminants D50 and D65 as seen through generic architectural glass (CIE 2009).

As you might expect from architectural glass, the slight greenish tint is due to absorption of red light. This increases the effective color temperature of the incident D50 and D65 daylight to 5100 and 6600 kelvins, respectively.

Figure 6

FIG. 6 – Indoor Daylight ID50 spectrum (5100K)

Figure 7

FIG. 7 – Indoor Daylight ID65 spectrum (6600K)

What makes this interesting for lighting designers is that architects and clients may obsess over the need for “high-CRI lighting” in their buildings. If the design involves both electric lighting and daylighting, one response could be to ask about the CRI of the building glass. If it is less than 90, there may be little point in worrying about the electric lighting.

Lamp with a CRI of 90 or above, good. Lamp with a CRI of less than 80, bad. This, plus the knowledge that “only a difference of some five points in the index is considered meaningful” is likely all you need to know (or talk) about CRI for most lighting design projects.

References

BSI. 2011. BS EN 410:2011, Glass in Building. Determination of Luminous and Solar Characteristics of Glazing. London, United Kingdom: British Standards Institution.

CIE. 1986. Guide on Interior Lighting. CIE 29:2-1986. Vienna, Austria: CIE Central Bureau.

CIE. 1995. Method of Measuring and Specifying Colour Rendering Properties of Light Sources. CIE 13.3-1995. Vienna, Austria: CIE Central Bureau.

CIE. 2007. Colour Rendering of White LED Light Sources, CIE 177:2007. Vienna, Austria: CIE Central Bureau.

CIE. 2009. Indoor Daylight Illuminants. CIE 184:2009. Vienna, Austria: CIE Central Bureau.

Colombo, E. 2013. “Opinion: Choose a Colour Rendering Metric Now,” Lighting Research and Technology 45:520.

EPA. 2013. ENERGY STAR Program Requirements for Integral LED Lamps. Washington, DC: Environmental Protection Agency.

Gombos, K., et al. 2008. “Proposal for an Indoor Daylight Illuminant,” Color Research & Application 34(1):18-25.

Speier, I., and M. Salsbury. 2006. “Color Temperature Tunable White Light LED System,” Proc.  6th International Conference on Solid State Lighting, SPIE Vol. 6337, 63371F.

Thornton, W. A. 1971. “Luminosity and Color Rendering Capability of White Light,” Journal of the Optical Society of America 61:1155-1163.

van Trigt, C. 1999. Color Rendering, a Reassessment. Color Research & Application 24(3):197-206.

Verstegen, J. M. P. J. 1974. “New Class of Phosphors for ‘Deluxe’ Fluorescent Lamps,” Lighting Research and Technology 6:31-32.

Whitaker, T. 2005. “TIR’s LEXEL Platform Provides Cost-Effective LED Lighting,” LEDs Magazine, July, pp. 19-21.

Whitaker, T. 2010. “Search Continues for Replacement for Color Rendering Index,” LEDs Magazine, July/August, p. 59.