Proposing a Method Reducing UGR Calculations for LED Luminaires with Diffusers in Indoor Applications

Featured Application: Proposing a method reducing UGR calculations for LED luminaires with diffuser (uniform luminous intensity distribution). Abstract: This paper studies discomfort glare in indoor LED lighting systems and proposes a method reducing UGR calculations. There are several methods for expressing discomfort glare, among which the uniﬁed glare rating (UGR) method is one of the most popular. It is safe to indicate that the UGR method was introduced before LED lighting became popular in indoor applications. Thus, the UGR method is not versatile currently. In this paper, key limitations of the UGR method are discussed. A method for reducing UGR calculations based on UGR tables for both LED luminaires with diffusers (uniform luminous intensity distribution) and with lenses is proposed. The method is mathematically discussed and its applicability for luminaires is discussed. Also, the luminaires are put through a test and the results from the tests are presented and the feasibility of the method separately from mathematical proofs is examined. The proof and case studies both indicate that the method can be used accurately for luminaires with diffusers. They also illustrate that the accuracy of applying the proposed method to the LED luminaires with lenses is not satisfactory. Simulations are carried out using DIALux lighting design software.


Introduction 1.LEDs as Low-Energy Indoor Light Sources
In 2015, the International Energy Agency (IEA) reported that 14% of residential electrical energy consumption in Organisation for Economic Co-operation and Development (OECD) countries was related to lighting.For non-OECD countries, the consumption share is usually higher [1].However, in 2021, the lighting share in building energy consumption became 5% and continues to decrease as a result of measures such as using LEDs instead of less energy-efficient bulbs [2] and electrification of heating systems [3].The use of LEDs in lighting systems is the result of the energy performance of buildings directive (EPBD)'s nearly zero-energy concept for European countries, which other countries have also adopted [1].Despite the mentioned decrease in energy consumption share, the absolute amount of energy consumption by lighting continues to increase [3].This indicates the increasing use of lighting around the world.
During the last decades, various energy initiatives have been undertaken.The use of incandescent lamps is discouraged.Energy-efficient light sources such as light-emitting diodes (LED) and organic light-emitting diodes (OLED) have been introduced.The EU has published "Lighting the Future-Accelerating the deployment of Innovative Lighting Technologies", proposing to apply policy innovations in Europe for the pervasive use of solid-state lighting products.The mentioned products will play a significant role in energy saving in the near future [1,4].The efficacy of LEDs can reach 200 lm/W.Thus, these light sources are very attractive and are regarded as energy-efficient light sources [5].However, LEDs are relatively small and compact light sources, which can deliver very high luminous intensity, resulting in significant discomfort glare.Dealing with these high intensities is an ongoing challenge faced by many manufacturers of luminaires.
Discomfort glare can be expressed numerically by the unified glare rating UGR, which was first adopted to lighting standards in 1995.Over the years, the UGR has been applied as a primary and common metric for quantifying discomfort glare [6,7].
This paper discusses challenges related to the use of the UGR evaluation method, especially in modern LED lighting systems, and proposes a novel method for reducing UGR calculations.

Glare Evaluation by UGR
UGR is a glare quantification tool which is based on contrast noticeable by an observer.Figure 1 shows an example of a space with highly contrasting LED luminaires.
The efficacy of LEDs can reach 200 lm/W.Thus, these light sources are very attractive and are regarded as energy-efficient light sources [5].However, LEDs are relatively small and compact light sources, which can deliver very high luminous intensity, resulting in significant discomfort glare.Dealing with these high intensities is an ongoing challenge faced by many manufacturers of luminaires.
Discomfort glare can be expressed numerically by the unified glare rating UGR, which was first adopted to lighting standards in 1995.Over the years, the UGR has been applied as a primary and common metric for quantifying discomfort glare [6,7].
This paper discusses challenges related to the use of the UGR evaluation method, especially in modern LED lighting systems, and proposes a novel method for reducing UGR calculations.

Glare Evaluation by UGR
UGR is a glare quantification tool which is based on contrast noticeable by an observer.Figure 1 shows an example of a space with highly contrasting LED luminaires.
For an indoor observer, the UGR formula is given as [8]: where   is the background luminance in cd/m 2 ,   is the luminance of the luminaire in cd/m 2 ,   is the solid angle of the luminaire in sr,   is the Guth position index of the luminaire, and n is the number of luminaires in the observer's field of view.Assumptions and limitations of the UGR method: • UGR cannot be applied for completely indirect lighting systems or luminous ceilings [8].

•
It is assumed the room walls are uniformly illuminated [8].
Alongside the mentioned limitations, discomfort glare evaluation is a challenging effort.Quek et al. said that the two categories of effects which lead to discomfort glare are saturation and contrast.UGR is only based on contrast effect.The performance of metrics using contrast effect is poor in brighter scenes.Quek et al. also stated that applicability of metrics using two or more effects simultaneously is limited by their dataset development For an indoor observer, the UGR formula is given as [8]: where L b is the background luminance in cd/m 2 , L i is the luminance of the luminaire in cd/m 2 , ω i is the solid angle of the luminaire in sr, p i is the Guth position index of the luminaire, and n is the number of luminaires in the observer's field of view.Assumptions and limitations of the UGR method: • UGR cannot be applied for completely indirect lighting systems or luminous ceilings [8].

•
It is assumed the room walls are uniformly illuminated [8].
Alongside the mentioned limitations, discomfort glare evaluation is a challenging effort.Quek et al. said that the two categories of effects which lead to discomfort glare are saturation and contrast.UGR is only based on contrast effect.The performance of metrics using contrast effect is poor in brighter scenes.Quek et al. also stated that applicability of metrics using two or more effects simultaneously is limited by their dataset development scope [12].Fotios and Kent argue that in experiments, biases can affect reliability and validity of data.They indicate that among four psychological procedures used for explicit quantitative measurement, adjustment and category are the most prevalent procedures in discomfort glare evaluation.The two mentioned procedures are subjected to biases such as stimulus range, anchor, and order effect biases [13].Consequently, UGR as a metric which developed using adjustment procedure [14] is subjected to the mentioned biases.Pierson et al. suggested that there is possibility of the existence of some unknown psychological and physiological factors influencing discomfort glare perception.They also noted that the intensity of each influencing factor cannot be understood and explained accurately.Therefore, an exact correlation between discomfort glare as assessed by the metrics and as perceived by subjects cannot be expected.They also, based on previous research, reported that even age and culture have effects on discomfort glare [15].Wienold et al. based on investigating 22 glare prediction metrics, stated that still there is not a reliable metric that considers all scenarios [16].Słomi ński obtained UGR for bare and diffusercovered LED luminaires.UGR was calculated using both luminance measured by an image luminance measuring device (ILMD) and luminous intensity numbers in the luminaire IES file.The results showed for both luminaires that the UGRs calculated by the two methods have substantial differences.He recommended correct luminance distribution measurement according to human eye geometry and adding a factor to the UGR formula related to luminaire construction, particularly for luminaires with non-uniform luminance distribution on the whole surface [9].
According to limitations of the existing UGR method as pointed out above, alternatives such as the model based on receptive fields introduced in [17] are proposed.The model better interprets visual discomfort, especially in the case of luminaires with non-uniform luminous intensity distribution.However, the model still must be improved for indirect viewing angles and different background luminance levels.As mentioned previously, UGR is a widely accepted standard visual discomfort rating index, but according to limitations, researchers are working to enhance its applicability [10].T. Porsch et al. [18] introduced software for automated UGR calculation using imaging luminance measuring devices (ILMDs).In [19], evaluations were made by comparing UGRs obtained with the original method.In this research, the luminaires with non-uniform luminous intensity distribution were divided into smaller subdivisions.The results showed that UGR was not robust for luminaires with non-uniform luminous intensity distribution glare evaluation.Therefore, it is necessary to perform further research to find a robust method for glare evaluation for both luminaires with uniform and non-uniform luminous intensity distribution.Scheir et al. [20] defined an improved UGR standard method which took into account pixels of a luminaire with luminance above a threshold as the luminating part.The other pixels were considered as background luminance.The method can better describe the glare felt by observers than the standard CIE-proposed UGR calculation method.In [21], UGR for conventional fluorescent luminaires was obtained by three methods using DIALux software, LMK software, and a Python-based software.UGR was calculated by the mentioned methods for the same spaces.The results' similarity showed the applicability of UGR for luminaires with uniform luminous intensity distribution such as the fluorescent luminaires.Tyukhova et al. [22] used the UGR method for small sources to address glare produced by some types of LED sources via combining measurement and software analysis.Their study did not address LED luminaires with multiple light sources (LED matrix) or when several light sources were present in the field of view simultaneously.Sharma et al. [23] demonstrated the importance of light source positioning for controlling UGR.Additionally, they showed the correlation between UGR and observers' visual discomfort perception by simultaneously using a qualitative method.Wolska and Sawicki suggested a procedure for using high dynamic range imaging (HDRI) techniques for UGR assessment by determining the best exposure time.Their method uses ILMD.Consequently, their proposed procedure requires the use of technical devices [6].
The mentioned studies showed that, despite limitations of the UGR method in some conditions as explained previously, research to improve UGR to cope with different scenarios and conditions, especially modern luminaires such as LEDs, is plentiful.Proposed methods for UGR calculation improvement did not address reducing the number of calculations for obtaining UGR.
The reduction and simplification of the calculations and models have received attention in science fields as can be seen in [24][25][26][27].Due to the prevalence of UGR application, it is useful to devise methods for decreasing the number of UGR calculations.
The proposed tabular method makes it possible to decrease UGR calculations by using the UGR table and shows the UGR table's applicability for points inside a room which the table was not originally designed for.

Paper Structure
In this paper, we investigate whether approaches and methods with fewer calculations can be proposed for UGR calculation for different LED luminaire types.In Section 2, the proposed tabular method for UGR calculation for LED luminaires is introduced.Then, the feasibility of the method is proved.Afterwards, the method is verified by comparing the method's results with those obtained from the simulation.Finally, in Section 3, the UGR calculation method proposed in this paper and its validity are summarized and presented in brief, while future research areas in the field are suggested.

Introduction to the Method
UGR calculation for points inside an indoor space in lighting design is essential.Since, as stated in AS/NZS 1680, not only criteria such as minimum needed work-plane illuminance for a particular task should be met, but also UGR criterion observation for users' visual comfort is mandatory [28].By referring to Equation (1), it can be deduced that UGR calculation for a large number of luminaires in a room requires a high number of calculations.The reason is the number of summands in front of the sigma operator, which is equivalent to the number of luminaires.An approach that drastically reduces the calculation number can proportionally reduce the time and energy consumption for UGR calculation.
According to CIE 117-1995 [8], UGR tables are provided by luminaires' manufacturers.This table is designed as a very useful and efficient tool for fast estimation and calculation of glare only at the midpoints of the room's walls.A sample of a UGR table is illustrated in Table S1 in the Supplementary File.By finding the desired X and Y dimensions, the corresponding luminaire direction, and the surface's reflection factor, the UGR can be obtained.
Where intended number exists in the table, it can be used directly, and in case of not existence, simple interpolation (extrapolation) using numbers in the UGR table can be used as a tool to estimate UGR.(All interpolations (extrapolations) used in this research are well-known linear (1D) [29] (for 2 known points) or bilinear (2D) [30] (for 4 known points) interpolations (extrapolations).) Through theoretical and experimental investigations in this paper, it is proved that in addition to the points in the room for which the UGR table is designed, this table can be used for all points on the room's midlines to obtain the UGR (Figure 2).
A method for using numbers in the UGR table for obtaining UGR on the room's midlines is proposed.In the proposed method, by simply using numbers in the UGR table issued by luminaires' manufacturer, we attempt to obtain UGR for every point on the room's midlines.The proposed method is fast because, as mentioned previously, for a room with dimension XH × YH, for the very first points of the observer's movement path, the number of summands in front of the sigma operator in Equation ( 1) is n = X × Y (number of luminaires in front of the observer).However, in the proposed method, only one 2D interpolation or extrapolation between related adjacent numbers of the table is required.A method for using numbers in the UGR table for obtaining UGR on the room's midlines is proposed.In the proposed method, by simply using numbers in the UGR table issued by luminaires' manufacturer, we attempt to obtain UGR for every point on the room's midlines.The proposed method is fast because, as mentioned previously, for a room with dimension XH × YH, for the very first points of the observer's movement path, the number of summands in front of the sigma operator in Equation ( 1) is  = X × Y(number of luminaires in front of the observer).However, in the proposed method, only one 2D interpolation or extrapolation between related adjacent numbers of the table is required.

Method Explanation
When an observer is located at a point with a certain distance from the walls, it is rational to have an estimation for UGR at that point.Reducing the number of calculations of Equation (1) for that point would be desirable for reducing UGR calculation time and energy.
A UGR table is a table as illustrated in Supplementary Table S1.As mentioned before, it is designed for obtaining UGR in the middle of each wall in a room.In the proposed method, when a point on a room midline is considered, for obtaining UGR for observer at that point, it is assumed that the space behind the observer is removed and UGR is obtained from the UGR table for a new room with space equal to the remaining space of the main room (The new room dimensions in the observer's field of view are smaller, but in other directions, dimensions of the main and new rooms are equal) in front of the observer (Figure 3).If the exact corresponding dimensions do not exist in the table, the related UGR can be interpolated or extrapolated.Illustrating the points UGR table was originally designed for (black arrows), and points for which the UGR table can be used (red lines).

Method Explanation
When an observer is located at a point with a certain distance from the walls, it is rational to have an estimation for UGR at that point.Reducing the number of calculations of Equation (1) for that point would be desirable for reducing UGR calculation time and energy.
A UGR table is a table as illustrated in Supplementary Table S1.As mentioned before, it is designed for obtaining UGR in the middle of each wall in a room.In the proposed method, when a point on a room midline is considered, for obtaining UGR for observer at that point, it is assumed that the space behind the observer is removed and UGR is obtained from the UGR table for a new room with space equal to the remaining space of the main room (The new room dimensions in the observer's field of view are smaller, but in other directions, dimensions of the main and new rooms are equal) in front of the observer (Figure 3).If the exact corresponding dimensions do not exist in the table, the related UGR can be interpolated or extrapolated.
In this research, the Y direction is assumed to be the direction of the observer's movement, and the X direction is assumed to be perpendicular to the direction of the observer's movement.
As previously indicated, for a sample luminaire, the UGR table is illustrated in the Supplementary File as Table S1.If it is assumed the room has dimensions of 8 m × 8 m, then, based on a luminaire height of H = 2 m from the observers' eye, the length and width dimensions of the room will be 4H × 4H.Reflection factors for the ceiling, walls, and floors are 70, 50, and 20, respectively.Viewing direction is at right angles to the lamp axis.Therefore, for the observer position at the wall midpoint, UGR will be 21.2; for the observer position in the center of the room with the same viewing direction, UGR will be 18.9.The latter number is obtained based on the proposed method and corresponds to a room with dimensions of 8 m × 4 m (4H × 2H), while in reality, the dimensions of the room have not been changed.In this research, the Y direction is assumed to be the direction of the observer's movement, and the X direction is assumed to be perpendicular to the direction of the observer's movement.
As previously indicated, for a sample luminaire, the UGR table is illustrated in the Supplementary File as Table S1.If it is assumed the room has dimensions of 8 m × 8 m, then, based on a luminaire height of H = 2 m from the observers' eye, the length and width dimensions of the room will be 4H × 4H.Reflection factors for the ceiling, walls, and floors are 70, 50, and 20, respectively.Viewing direction is at right angles to the lamp axis.Therefore, for the observer position at the wall midpoint, UGR will be 21.2; for the observer position in the center of the room with the same viewing direction, UGR will be 18.9.The latter number is obtained based on the proposed method and corresponds to a room with dimensions of 8 m × 4 m (4H × 2H), while in reality, the dimensions of the room have not been changed.

Mathematical and Practical Proof of Feasibility of the Method
As explained previously, in the proposed method, UGR table numbers or their 2D interpolation or extrapolation are used as the UGR for the observer's position on the room's midline.Therefore, it is assumed that UGR for an observer at any point on a room's midline is equal to UGR for an observer in a room with size equal to the space in front of the first observer when the second observer is placed in the middle of one wall looking toward the opposite-side wall in a smaller room.For proving the validity of idea, it should be proved that UGRs obtained for both observers are adequately close to each other.
To start the proof, it is necessary to refer to Equation ( 1), the well-known CIE 117:1995 formula for UGR calculation [8].The formula takes into account all but one of the variables

Mathematical and Practical Proof of Feasibility of the Method
As explained previously, in the proposed method, UGR table numbers or their 2D interpolation or extrapolation are used as the UGR for the observer's position on the room's midline.Therefore, it is assumed that UGR for an observer at any point on a room's midline is equal to UGR for an observer in a room with size equal to the space in front of the first observer when the second observer is placed in the middle of one wall looking toward the opposite-side wall in a smaller room.For proving the validity of idea, it should be proved that UGRs obtained for both observers are adequately close to each other.
To start the proof, it is necessary to refer to Equation ( 1), the well-known CIE 117:1995 formula for UGR calculation [8].The formula takes into account all but one of the variables which are directly related to the observer's field of view.Therefore, except for one variable, all variables-L i as the luminance of the luminaire i, ω i as the solid angle of the luminaire i, and p i as the Guth position index of the luminaire i-are directly related to the luminaires which are in the observer's sight.The only parameter which considers the effect of luminaires out of the observer's sight is L b , the background luminance.Therefore, for two mentioned observers in two different rooms, it is enough to prove that the background luminances are close to each other such that when using them in Equation ( 1), the resulted UGRs are reasonably close to each other.
Consider two mentioned observers.One is in a room in a position which places some of the luminaires in front of them (in their field of view), and some of the luminaires behind them (Figure 4a).The other observer exists in a room containing only the luminaires which are in the first observer's field of view.The nearest wall behind the second observer is at a distance of H/2 from the first row of the luminaires in front of them (Figure 4b).The arrangement of the luminaires is as the standard arrangement explained in CIE 117:1995 [8].
of the luminaires in front of them (in their field of view), and some of the luminaires behind them (Figure 4a).The other observer exists in a room containing only the luminaires which are in the first observer's field of view.The nearest wall behind the second observer is at a distance of H/2 from the first row of the luminaires in front of them (Figure 4b).The arrangement of the luminaires is as the standard arrangement explained in CIE 117:1995 [8].For UGR calculation in the two rooms illustrated in Figure 4a,b, Equation ( 1) can be stated as: UGR main room ( f or observer position) UGR corresponding smaller room(Observer is located in the middle of one wall) = 8 log 0.25 where the parameters with indices a and b correspond to rooms a and b as illustrated in Figure 4a,b.
It should be proved that: UGR main room ( f or current observer s position) = UGR corresponded smaller room(Observer is located in the middle o f one wall) Therefore, from Equations ( 2)-( 4): Or: As a result, by proving: the equality in Equation ( 6) will be proved.
As can be deduced from Figure 4a,b, since there are same luminaires in both observers' fields of view in the mentioned figures, Equations ( 7)-( 9) are true.
However, L b in Equation ( 1) is the only variable which is not related only to the luminaires in the observer's field of view.Therefore, the aim is to prove that the values of L b for both of the rooms and observers' positions are reasonably close to each other.For this aim, the range of changes in background luminance will be obtained based on observers' position changes; we also examine whether background luminance change can be ignored.
In this research, the CIE 190:2010 [31] method for L b calculation will be used: where, E WID is an indirect component of the illuminance on the wall produced by luminaires.In turn, E WID is: where F UWID is the indirect utilization factor for walls, N is the number of luminaires, A W is the total area of walls (m 2 ) (W stands for wall) between the reference plane and the luminaires' installation plane, and Φ 0 equals 1000 lm (reference luminous flux) when the luminaire's luminous intensity distribution is given in cd per 1000 lm.It should be noted that in case of giving luminous intensity distribution based on luminaire luminous flux, the luminous intensity distribution numbers should be divided by Φ luminaire /1000 to change the luminous intensity unit to cd per 1000 lm, which makes the use of Equation ( 12) possible.
As the observer moves toward the opposite-side wall in the main (first) room, the size of the smaller (second) room corresponding to the observer's position in the first room (In the second room, the observer's distance from the opposite-side wall is the same as the first observer, but the nearest wall behind the second observer is exactly H/2 distance from the first row of the luminaires in front of them) changes as illustrated in Figure 4a,b, and the length of the second room in the first observer's path decreases.Meanwhile, the number of luminaires in both observers' fields of view decreases.
If there is a room with dimensions of XH × Y L H (the main room) and an observer moves in Y direction, it can be said when the observer in the main room is in the position (Y 0 + ∆Y)H (distance from opposite-side wall), and for H = 2 m, the number of luminaires effecting background luminance would be: and the area of wall which reflects background luminance is: where Y 0 is the smallest distance of the observer from the opposite-side wall, ∆Y is added to Y 0 for obtaining the Y corresponding to the observer's position, and Y L is the length of the main room (Y and X should be multiplied by H for obtaining dimensions in meters).
For the corresponding smaller room: and By putting aside A w , N, and Φ 0 , the remaining part of Equation ( 12) is related to the indirect utilization factor of the walls, which is the next part of the formula to be considered.
Indirect utilization factor for walls is [31]: In Equation ( 17), distribution factors F DF , F DW , and F DC are defined as [31]: F DF and F DW are derived from Φ ZL , which is the sum of the cumulative CIE zonal fluxes multiplied by the geometric factor values for luminaires in the standard array.Geometric factors are related to the room's dimensions.Therefore, they change in small rooms when the observer's position changes.However, F DC is derived from R ULO , which is only related to the luminaire's luminous intensity specifications.Then, it is independent of the observer's position in the room.Based on changing parts of F UWID as the observer's position changes, its range of change will be investigated as follows.
Φ ZL as an important part of F UWID is expressed as: Φ ZL1 , Φ ZL2 , Φ ZL3 , and Φ ZL4 are only related to each luminaire's luminous flux properties and are not related to the observer's position or the room's dimensions.Then, they are independent of the observer's position in the room.By rewriting F UWID formula (Equation ( 17)) considering the change in the observer's position: Following that, for finding the changing range for F UWID , the changing range for the right side of Equation ( 22) will be investigated in detail.
For F DF , on the right side of Equation ( 22), it can be written: As previously mentioned, in Equation ( 23), Φ ZLn s are related to the luminaires' specifications and do not change when the room's specifications change.Geometric factors (F GLn s) are related to the room's dimensions and ranges of their changing with room dimensions are illustrated in Table 1.From Table 1, it can be understood that when a room's Y dimension increases, F GL1 and F GL2 decrease and F GL3 and F GL4 increase.From each row of Table 1 and considering Equation ( 21), it can be deduced that for luminaires whose Φ ZLn values are close to each other, the increasing and decreasing of related geometric factors can cancel each other and CIE zonal flux (Φ ZL ) can change little when a room's dimensions change.Knowing that, and also since in Equation ( 18), F DF is CIE zonal flux (Φ ZL ) divided by 1000 lm, F DF change based on changes in room dimensions is not prominent.
From Equation (19) and keeping in mind that R DLO is not dependent on room dimensions: Consequently, F DW change which is equal to F DF change but with a negative sign is not also prominent.
Since luminaires' R ULO value is usually negligible, the F DC also can be assumed to be a negligible value.Therefore, F DF and F DW are the only distribution factors influencing the change of the indirect utilization factor of the walls.
Table 2 illustrates a range of transfer factors' change according to the change in a room's dimensions.Validity of using data illustrated in Tables 1 and 2 requires maintaining a trend of change for geometric and transfer factors linearly, both inside and outside of the studied range (X = 2H, 4H, 8H, and 12H, Y = 4H, 6H, and 8H).Therefore, linear regression was applied to trends of mentioned factors for inside and outside of the studied range where data were available, and R 2 -values of regression for all factors for a common range are illustrated in Tables S3 and S4.Based on data illustrated in Tables S2 and S3, the average of R 2 -value for linear regression of geometric factors is 0.8701 and for transfer factors is 0.9651.The mentioned R 2 -values show acceptable linear correlation between room dimension in the Y direction and geometric and transfer factors.It shows the possibility of generalization of the results obtained in the studied dimension range for outside the range.
From Equation (22), the average change of the indirect utilization factor for walls based on the change of distribution and transfer factors can be stated as Equation ( 25): From Equations ( 19), (20), and (25): Since, in many luminaires, luminous intensity in the upper hemisphere is near zero, R ULO also equals zero.Therefore, Equation ( 26) would be: By expanding Equation ( 27): Multiplication of two differential values makes very small values.Therefore, the terms including multiplication of two differential values are removed from Equation (28): By subtracting F UWID (Y 0 ) obtained by Equation (17) from F UWID (Y 0 +∆Y) in Equation ( 29): Previously, it was concluded that ∆F DF does not have a prominent value, hence by removing the terms including ∆F DF in Equation ( 30): In other words: By considering Table 2, it can be deduced that: and From Equations ( 18) and ( 21) and by investigating luminaires with their luminous intensity distribution in the lower hemisphere around the luminaire, F DF can be expressed as follows: Usually in luminaires with uniform luminous intensity distribution in bottom hemisphere: Also, by estimation it can be assumed that Equation (35) for F DF calculation can be expressed as: From Equation (36) it can be concluded that: Hence, from Equations (37) and (38): Also, R DLO and F DF are numbers smaller than 1.Hence, from Equations ( 32)-(34), and (39): According to Equation ( 12), the indirect component of the illuminance on the walls in terms of its change based on room dimension change would be: N A W (Y 0 ) for a sample dimension changing range (8 m to 16 m) in both directions of X and Y can be seen in Figure 5.  From Figure 5 it can be deduced that the average change of Therefore: For finding the change of the indirect component of the illuminance on the walls according to the change in the room's dimensions: From Figure 5 it can be deduced that the average change of N A W in the considered dimension change range is 0.125.Consequently, Equation (41) would be: Therefore: For finding the change of the indirect component of the illuminance on the walls according to the change in the room's dimensions: The right side of Equation (44) can be rewritten as: Therefore, average change in E WID for a ±4 m change in observer position would be: From Equation (17), transfer factors F T,FW , F T,WW-1 , and F T,CW and distribution factors F DF , F DW , and F DC are used for obtaining the indirect utilization factor for walls (F UWID ).The average of transfer factors for usual calculations as mentioned in Table 5 in CIE 190:2010 are 0.135, 0.251, and 0.371, for F T,FW , F T,WW-1 , and F T,CW , respectively.F DF and F DW are smaller than one and F DC for luminaires with their luminous intensity distribution in lower hemisphere is approximately zero.Therefore: Since, as shown in Equation ( 11), E W ID (Y 0 ) is proportional to background luminance average, from Equation (46) it is known that for finding average change in background luminance, F UW ID (Y 0 ) and N A W (Y 0 ) are needed.Also, as mentioned before, R DLO is smaller than one and an average value of 0.5 can be considered for it.Therefore, from Equation (39) and average value of R DLO , in Equation (39) F DF would be equal to 0.225.Hence, F UW ID (Y 0 ) in Equation ( 48) would be equal to 0.099.Then, ∆E W ID (Y 0 ), based on average values of N/A W and F UWID in the studied range, would be: Ultimately, the ratio of background luminance change to its absolute value is: By applying Equations ( 7)-( 9) to Equation (51): = 8 log 51.76−1.851.76 = 0.13 (52) Therefore, on average, while an observer moves from a distance of 16 m from the opposite-side wall toward it, at distance of 8 m, only an approximate change in UGR of 0.13 units occurs for them.This change is far from the 3 UGR units of change which is perceptible according to the Hopkinson criterion.Remembering that previously it was shown that geometric and transfer factors have linear correlation with dimension change, the obtained result can be generalized and extrapolated for dimensions greater or smaller than values studied in this proof.
Then, having distribution and transfer factors, N/A W , and Φ 0 = 1000 lm, values for E WID and in turn, L b (background luminance) for {X, Y | (X = 12, 8, 2) & (Y = 12, 2, 1)} can be obtained.The background luminances for mentioned Xs and Ys obtained by manual calculation and simulation for different types of luminaires are illustrated in Tables 3 and 4, respectively.In the mentioned tables, for each X as the observer's path width, deviation of background luminance for Y = 1 and 2 from background luminance at the start point of path (Y = 12) is also illustrated as a percentage of background luminance at the starting point.Due to what is mentioned in CIE 117:1995 about insensitivity of the UGR to the errors produced in L b calculation (in which +33% error in L b leads to only 1-unit error in UGR), by using data in Table 3 for manual calculation and Table 4 for simulation of L b , it can be deduced that starting from a distance of 24 m from the opposite-side wall, for the manual calculation example ending at a distance of at least 4 m and for the simulation ending at a distance of less than 2 m from the opposite-side wall, the UGR calculation error is limited to only 1 unit, which is far from the 3 units of UGR introduced in the Hopkins criterion for perceptible UGR difference.
Also, for the previously obtained range of UGR change for 4 ≤ X ≤ 8 and 4 ≤ Y ≤ 8 for an average value of R DLO equal to 0.5, the UGR change was obtained as 3.67%, which corresponds to 0.13 units of UGR error.
Hence, in addition to L i , ω i , and p i , which are totally independent of the luminaires and space out of an observer's sight, background luminance also has relatively small changes according to the number of luminaires and change in dimensions of the space both inside and outside the observer's sight.Therefore, in addition to Equations ( 7)-( 9), it was proved that Equation ( 10) is also true with approximation and it can be assumed that UGR obtained for the main room and corresponding smaller room illustrated in Figure 4a,b are close enough and UGR table can be used for points on the room's midline.
More practical samples for evaluating the feasibility of the method will be examined in Section 2.4.
The tabular method was originally designed for UGR calculation in the middle of one wall in a room.When using it at a point inside a room on the room's midline other than wall midpoints, it should be shown that for a smaller room size corresponding to the space in front of the observer in a larger room when observer's location is other than the room's margins, by using Equation ( 1), the obtained UGR for smaller room is close enough to the real UGR for the observer in the larger room.Therefore, the only difference between the parameters of the main and corresponding rooms exists in background luminance.Then, by showing acceptable difference between the two background luminances, it was shown that the proposed tabular method is applicable for an observer's location on the room's midline other than the wall midpoints.
All of these lead to the valid use of UGR table for changes in an observer's position in a room instead of its conventional use for observers in the middle of one wall of the room.Also, since the proposed method has been proved mathematically, repeatability of the method's results is obvious.
Since the proposed method uses the UGR table introduced to CIE standards [8], the limitations of the UGR calculation method in CIE 117:1995 such as limitations for indirect lighting, luminaires' dimensions, luminous intensity distribution uniformity [4], etc. are applied to the proposed method.Iacomussi et al. suggested that initiatives in formulations regarding investigation and evaluation of solid-state light source (SSL) discomfort glare is necessary.They also concluded that existing criteria for discomfort glare evaluation can be applied to LED luminaires provided that they are similar to conventional fluorescent luminaires in luminous intensity and luminance distribution.They recommended equipping SSLs with diffusers [7].Additionally, position index is used in UGR calculation and was developed by examining light sources with uniform luminous intensity; it appears that it is not suitable for luminaires without a set circular geometry and a constant spectrum.In addition, the results of studies on position index of luminaires with non-uniform luminous intensity distribution were contradictory.Furthermore, the form and structure of modern luminaires make it problematic to calculate the luminance needed for UGR calculation from luminous intensity distribution and apparent area of the luminaire.It is not usually easy to specify an apparent area for LED luminaires, particularly for those with complex retrofits [4].The aforementioned limitations make applicability of the UGR method for evaluation of discomfort glare for LED luminaires vague.
Here, by performing case studies, the proposed method will be examined.First, eight types of LED luminaire with diffusers are studied.The luminous intensity distribution for all the eight luminaire types can be seen in Figure 6a-h.
Although feasibility of the method is shown in Section 2.3 for luminaires with diffusers (uniform luminous intensity distribution), luminaires with variety in technical specification were selected for practical investigation of the method's efficacy.The technical specifications of selected luminaires are illustrated in Table S5.
For studying eight types of LED luminaires with diffusers, UGR at the eye height of 1.2 m is simulated in DIALux software.The room size is 5H × 5H and the observer moves on the room's midline from the left wall to the right wall parallel to the luminaires' axis.For meticulous calculation, UGR is simulated each 0.5 m of the observer's path.Detailed results for all the types of luminaire are integrated and illustrated in Figure 7.
As can be seen in Figure 7, UGR changes for all LED luminaires with diffusers along the midline of the room were more or less similar.They started with a value which decreased after passing each row of the luminaires perpendicular to the observer's viewing direction (moving path).After reaching the last two meters, the decrement became such that the UGR was under 10 units of UGR and was not given by the simulation software.
Due to the prevalence of LED luminaires with diffusers, it is helpful to find methods to calculate UGR for this type of luminaire in a reliable manner, with the smallest possible number of calculations, and with acceptable accuracy.In this research, a tabular approach is proposed for the estimation of the UGR of LED luminaires with diffusers, as well as conventional luminaire types with uniform luminous intensity distribution for UGR calculation reduction with acceptable accuracy.Although feasibility of the method is shown in Section 2.3 for luminaires with diffusers (uniform luminous intensity distribution), luminaires with variety in technical specification were selected for practical investigation of the method's efficacy.The technical specifications of selected luminaires are illustrated in Table S5.

Case Study and Illustration of Implementation
For studying eight types of LED luminaires with diffusers, UGR at the eye height of 1.2 m is simulated in DIALux software.The room size is 5H 5H and the observer moves on the room's midline from the left wall to the right wall parallel to the luminaires' axis.For meticulous calculation, UGR is simulated each 0.5 m of the observer's path.Detailed results for all the types of luminaire are integrated and illustrated in Figure 7.As can be seen in Figure 7, UGR changes for all LED luminaires with diffusers along the midline of the room were more or less similar.They started with a value which decreased after passing each row of the luminaires perpendicular to the observer's viewing direction (moving path).After reaching the last two meters, the decrement became such that the UGR was under 10 units of UGR and was not given by the simulation software.
Due to the prevalence of LED luminaires with diffusers, it is helpful to find methods to calculate UGR for this type of luminaire in a reliable manner, with the smallest possible number of calculations, and with acceptable accuracy.In this research, a tabular approach is proposed for the estimation of the UGR of LED luminaires with diffusers, as well as conventional luminaire types with uniform luminous intensity distribution for UGR calculation reduction with acceptable accuracy.
For luminaires, manufacturers usually publish UGR calculating tables for obtaining UGR in the worst case, which corresponds to the midpoints of the walls.An example of a UGR table (luminaire with diffuser type A's UGR table) is illustrated in Supplementary Table S1.
By accurate investigation of the curves in Figure 7 along with UGR tables published by manufacturers of eight types of LED luminaires with diffusers, a relation between them can be understood.
As explained before, the proposed method uses each number in the UGR table related to a room with dimensions corresponding to the observer's position on a room midline.The X dimensions for two rooms are equal, but the Y dimension of the second (corresponding smaller) room is equal to the observer's distance from the opposite-side wall in the observer's field of view in the main room.
As showed mathematically in Section 2.3, the method works properly.For further validation of the method, examples for luminaires with luminous intensity distributions as in Figure 6a-h will be examined.By examining eight luminaires A to H and using two sets of the table data for X = 2H and 4H as well as Y = 2H, 3H, 4H, and 6H (as within the red border in Supplementary Table S1 and illustrated in Supplementary Tables S6-S13) for 2D interpolation or extrapolation, the similarities between the results of simulation For luminaires, manufacturers usually publish UGR calculating tables for obtaining UGR in the worst case, which corresponds to the midpoints of the walls.An example of a UGR table (luminaire with diffuser type A's UGR table) is illustrated in Supplementary Table S1.
By accurate investigation of the curves in Figure 7 along with UGR tables published by manufacturers of eight types of LED luminaires with diffusers, a relation between them can be understood.
As explained before, the proposed method uses each number in the UGR table related to a room with dimensions corresponding to the observer's position on a room midline.The X dimensions for two rooms are equal, but the Y dimension of the second (corresponding smaller) room is equal to the observer's distance from the opposite-side wall in the observer's field of view in the main room.
As showed mathematically in Section 2.3, the method works properly.For further validation of the method, examples for luminaires with luminous intensity distributions as in Figure 6a-h will be examined.By examining eight luminaires A to H and using two sets of the table data for X = 2H and 4H as well as Y = 2H, 3H, 4H, and 6H (as within the red border in Supplementary Table S1 and illustrated in Supplementary Tables S6-S13) for 2D interpolation or extrapolation, the similarities between the results of simulation and the proposed method for the mentioned types of luminaires are checked and confirmed (Tables S14-S21).X dimensions (2H and 4H) are selected since the room dimension 5H is near these two values and the ranges of Y from 2H, 3H, 4H, and 6H correspond to 4 m, 6 m, 8 m, and 12 m, respectively, which properly cover the observer's distance from 10 to 0 m from the opposite-side wall.
In Tables S14-S21, data obtained from simulation and 2D interpolation or extrapolation of UGR table data for X = 10 m and Y = {y| y = n × 0.5 m , n ∈ Z & 0 ≤ n ≤ 19} for each type of luminaire with a diffuser are illustrated.Data used in the UGR table are for X = (2, 4) H and Y = (2, 3, 4, 6) H.The deviation of 2D-interpolated or extrapolated data average from the simulated UGR average for luminaires A to H is illustrated in Table 5 as a percentage of 30 units of UGR.The maximum of absolute deviation is approximately 4.73% of 30 units of UGR.A perceptible change in UGR according to the Hopkinson criteria is 3 units, equaling 10% of 30 [19].Therefore, deviation in the results of the proposed method is in the acceptable range.Additionally, from Supplementary Figures S1-S8, it is obvious that UGR obtained by 2D-interpolated or extrapolated table data follows the simulated UGR very closely.Therefore, for reducing UGR calculations in the indoor space, table data can be used.Hence, UGR at the calculation points, other than the wall midpoints, can be obtained from UGR tables.Since in most applications, LED luminous flux must be formed to meet lighting needs, using lenses is one of the approaches for this aim [32].Therefore, lenses as optical accessories enhance luminous intensity in LED luminaires in the required direction.An LED luminaire with lenses is illustrated in Figure 8.
each type of luminaire with a diffuser are illustrated.Data used in the UGR table are for X = (2,4) H and Y = (2,3,4,6) H.The deviation of 2D-interpolated or extrapolated data average from the simulated UGR average for luminaires A to H is illustrated in Table 5 as a percentage of 30 units of UGR.The maximum of absolute deviation is approximately 4.73% of 30 units of UGR.A perceptible change in UGR according to the Hopkinson criteria is 3 units, equaling 10% of 30 [19].Therefore, deviation in the results of the proposed method is in the acceptable range.Additionally, from Supplementary Figures S1-S8, it is obvious that UGR obtained by 2D-interpolated or extrapolated table data follows the simulated UGR very closely.Therefore, for reducing UGR calculations in the indoor space, table data can be used.Hence, UGR at the calculation points, other than the wall midpoints, can be obtained from UGR tables.Since in most applications, LED luminous flux must be formed to meet lighting needs, using lenses is one of the approaches for this aim [32].Therefore, lenses as optical accessories enhance luminous intensity in LED luminaires in the required direction.An LED luminaire with lenses is illustrated in Figure 8.
Due to the popularity of LED luminaires with lenses, the UGR tabular calculation method proposed in this paper, with demonstrated validity for luminaires with diffusers (uniform luminous intensity distribution) and limitations explained in Section 2.3, is investigated in this case study to estimate efficacy of the method for luminaires with lenses.Although the feasibility of the method was shown in Section 2.3 for luminaires with diffusers (uniform luminous intensity distribution), luminaires equipped with lenses with variety in technical specifications were selected for this case study.The technical specifications of selected luminaires are illustrated in Table S22.Due to the popularity of LED luminaires with lenses, the UGR tabular calculation method proposed in this paper, with demonstrated validity for luminaires with diffusers (uniform luminous intensity distribution) and limitations explained in Section 2.3, is investigated in this case study to estimate efficacy of the method for luminaires with lenses.
Although the feasibility of the method was shown in Section 2.3 for luminaires with diffusers (uniform luminous intensity distribution), luminaires equipped with lenses with variety in technical specifications were selected for this case study.The technical specifications of selected luminaires are illustrated in Table S22.
For studying eight types of LED luminaires with lenses (luminous intensity distributions of the luminaires are illustrated in Figure 9a-h), UGR at an eye height of 1.2 m from the floor was simulated in DIALux software.The room size is 5H × 5H.The observer moves on the room's midline from the left wall to the right wall of the room.For meticulous calculation, UGR is simulated in 0.5 m intervals and the results for all the types of luminaire are divided into luminaires with mild changes and those with drastic changes of UGR during the observer's movement on the room's midline, as illustrated in Figure 10a,b.Mild UGR changes correspond to UGR changes in a range of less than or equal to 6 units when the observer moves along the room's midline.For drastic UGR changes, the range of change is more than 6 units.
As can be seen in Figure 10a,b, the change in UGR values for all LED luminaires with lenses along the midline of the room is more or less similar.They have an alternating behavior with a similar period.Luminaires E and C alternate drastically, while B, G, and F have smoother changes.For each luminaire type, the magnitude and the changes of the values are more or less similar along the observer's path for each UGR curve.
For studying eight types of LED luminaires with lenses (luminous intensity distributions of the luminaires are illustrated in Figure 9a-h), UGR at an eye height of 1.2 m from the floor was simulated in DIALux software.The room size is 5H × 5H.The observer moves on the room's midline from the left wall to the right wall of the room.For meticulous calculation, UGR is simulated in 0.5 m intervals and the results for all the types of luminaire are divided into luminaires with mild changes and those with drastic changes of UGR during the observer's movement on the room's midline, as illustrated in Figure 10a,b.Mild UGR changes correspond to UGR changes in a range of less than or equal to 6 units when the observer moves along the room's midline.For drastic UGR changes, the range of change is more than 6 units.As can be seen in Figure 10a,b, the change in UGR values for all LED luminaires with lenses along the midline of the room is more or less similar.They have an alternating behavior with a similar period.Luminaires E and C alternate drastically, while B, G, and F have smoother changes.For each luminaire type, the magnitude and the changes of the values are more or less similar along the observer's path for each UGR curve.
The points which can be deduced from Figure 10a,b can be summarized as follows: • Continuous curve represents lower contrast as UGR is a ratio of luminance (A, B, D, F, G and H).The points which can be deduced from Figure 10a,b can be summarized as follows: • Continuous curve represents lower contrast as UGR is a ratio of luminance (A, B, D, F, G and H).

•
The average UGR matters (UGR of G < UGR of B).

•
Discontinuous curve introduces a "torch" effect to the design (C and E).

•
If UGR is not available, there is a chance that its value reaches 30 (E).
Due to the prevalence of LED luminaires with lenses, it is helpful if methods can be found to calculate UGR for this luminaire type in a reliable manner, with the lowest possible number of calculations, and with acceptable accuracy.As for LED luminaires with diffusers, the tabular method proposed in this research is also investigated for UGR estimation for LED luminaires with lenses.However, previously in Section 2.3, it was anticipated that the proposed method likely does not have enough efficacy for luminaires with non-uniform luminous intensity distribution.
An example of a UGR table (luminaire with lens type A's UGR table) is illustrated in Table S23.
Similar to studying luminaires with diffusers, by accurately investigating the curves in Figure 10a,b along with UGR tables published by the manufacturers for eight types of LED luminaires with lenses, a relation between them can be understood.By examining eight luminaires A to H and using two sets of the table data for X = 2H and 4H and Y = 2H, 3H, 4H, and 6H (in the red border in Table S23, for example, and illustrated in Tables S24-S31) for 2D interpolation or extrapolation, the similarities between the proposed method and simulation results for the mentioned luminaires are checked (Tables S32-S39).X dimensions (2H and 4H) are selected since the room's X dimension 5H is near these two values and the range of Y from 2H, 3H, 4H, and 6H corresponds to 4 m, 6 m, 8 m, and 12 m, respectively, which properly covers an observer's distance from 10 to 0 m from the opposite-side wall in the room.
Tables S32-S39 show data obtained from 2D interpolation or extrapolation of table data and simulation data for X = 10 m and Y = {y| y = n × 0.5 m , n ∈ Z & 0 ≤ n ≤ 19 } for each luminaire with lenses.Data used in the UGR table is for X = (2, 4) H and Y= (2,3,4,6) H.The deviation of 2D-interpolated or extrapolated data average from the simulated UGR average, for luminaires A to H, is illustrated in Table 6 as a percentage of 30 units of UGR.As mentioned before, according to Hopkinson criteria, below 10% deviation of full-scale UGR changing (0 to 30) is acceptable.Therefore, the proposed method results are not in the reliable range.Also, from Figures S9-S16, it is obvious that UGR obtained by 2D-interpolated or extrapolated table data does not properly follow the simulated UGR.In addition, for further study and finding any possibility for adapting the proposed tabular method for LED luminaires with lens (with non-uniform luminous intensity distribution), difference between UGR obtained by tabular method and simulation, according to supplementary Tables S32-S39 are illustrated in Table 7 for distances from opposite-side wall equal to 4, 6, and 8 m separately.For mentioned distances from opposite-side wall, corresponding UGRs exist in UGR table and no interpolation or extrapolation in y axis is needed.From Table 7, it can be deduced that even at points that the UGR table directly gives UGR, there are great differences (more than 10% Hopkinson criterion as perceivable UGR change) between the proposed method and simulation results.Consequently, since the proposed tabular method is based on a UGR standard table which is constructed according to the standard UGR calculation method [4], and the result of the proposed method even for points the UGR table directly gives, the UGR value has great difference with results from simulation based on the standard UGR calculation formula.Hence, the UGR calculation method applicability for LED luminaires with lenses could be a unclear as stated in former research (e.g., [4,7]).Therefore, further research and efforts to enhance the usual method are needed to find ways for accurate UGR calculation for luminaires with non-uniform luminous intensity distribution, including LED luminaires equipped with lenses.Consequently, as a UGR-calculation-reducing method in indoor spaces along room midlines for luminaires with lenses, the table data cannot be used.This corresponds exactly to what was anticipated for proposed method efficacy for LED luminaires with lenses in Section 2.3.* Difference is incalculable since simulation software according to UGR exceeding of 30 units limit does not calculate the UGR.

Conclusions
In this paper, a comprehensive investigation was performed on the efficacy of the proposed easy-to-use and calculation-reducing method for UGR calculation for LED luminaires.UGR as a metric for evaluation of indoor lighting user visual discomfort is introduced in related standards, and advised UGR values should be observed by indoor lighting designers.Indoor lighting designers must simulate their designation to find out whether it meets the standard limits or not.The proposed method makes it possible to find by referral to the UGR table given by the manufacturer, or by an interpolation or extrapolation of the table numbers, the UGR of points on the room's midline.Efficacy and validation limits of the proposed method were mathematically proved and discussed in detail.Since LED luminaires have attracted great attention, it is crucial that approaches be found to calculate UGR while consuming less time and energy.For luminaires with diffusers or uniform luminous intensity distribution, the tabular method was proved, calculated, and simulated.The results confirmed the validity of this method.
For luminaires with lenses, as in the method proof anticipated, the tabular method failed to model the UGR accurately.It can be based on non-uniform luminous intensity distribution of these types of luminaires.Further research is needed for modelling and finding a fast and easy-to-use way to calculate UGR for LED luminaires with lenses.

Figure 1 .
Figure 1.Visual discomfort glare of luminaires installed in the ceiling in a closed space.

Figure 1 .
Figure 1.Visual discomfort glare of luminaires installed in the ceiling in a closed space.

Figure 2 .
Figure 2.Illustrating the points UGR table was originally designed for (black arrows), and points for which the UGR table can be used (red lines).

Figure 2 .
Figure 2.Illustrating the points UGR table was originally designed for (black arrows), and points for which the UGR table can be used (red lines).

Figure 3 .
Figure 3. Illustration of the part of the main room (in grey color) which is used as a new room with different dimensions for UGR calculation in the proposed method.

Figure 3 .
Figure 3. Illustration of the part of the main room (in grey color) which is used as a new room with different dimensions for UGR calculation in the proposed method.

Figure 4 .
Figure 4. (a) Observer in a position in a room in which some of the luminaires are in front of them (in their field of view), (b) observer in a room which contains only the luminaires that are in the first observer's field of view.

Table 1 .
Geometric factors changing when room dimensions change for X = 2~12 and Y = 4~8.(The negative numbers show decreasing and the positive numbers show increasing of the geometric factors with increasing Y).
for a sample dimension changing range (8 m to 16 m) in both directions of X and Y can be seen in Figure5.

Figure 5 .
Figure 5.    as a function of X and Y (room dimensions).

Figure 5 .
Figure 5. N A W as a function of X and Y (room dimensions).
is concluded that for an observer position change (or room dimension change) from Y = 16 m to 8 m, background luminance change average is 3.67% for average values of R DLO = 0.5 and F DF = 0.225.This value for background luminance change can lead to UGR change as: ∆UGR = UGR a − UGR b = 8 log 0.25

Figure 6 .
Figure 6.(a-h) (Luminaires types A-H) Luminous intensity distribution of the eight types of luminaires with diffusers studied in this research.

Figure 7 .
Figure 7. Eight LED luminaires' UGR curves (luminaires described in Figure 6a-h) along the room's midline when observer's viewing direction is parallel to the luminaires' axis.

Figure 7 .
Figure 7. Eight LED luminaires' UGR curves (luminaires described in Figure 6a-h) along the room's midline when observer's viewing direction is parallel to the luminaires' axis.

Figure 8 .
Figure 8.A sample LED luminaire with lenses.

Figure 8 .
Figure 8.A sample LED luminaire with lenses.

Figure 9 .Figure 10 .
Figure 9. (a-h).(Luminaires types A-H) Luminous intensity distribution of eight luminaires with lenses studied in this research.Figure 9. (a-h).(Luminaires types A-H) Luminous intensity distribution of eight luminaires with lenses studied in this research.

Figure 10 .
Figure 10.Eight LED luminaires' UGR curves (a) for luminaires with UGR mild changes (∆ max−min UGR ≤ 6), and (b) for drastic changes of UGR (∆ max−min UGR > 6) along the room's midline.(For luminaires described in Figure 9a-h along the room's midline when the observer's viewing direction is parallel to the luminaires' axis).

Table 2 .
Transfer factors changing with the changes in a room's dimensions for X = 2~12 and Y = 4~8.(All changes are in form of decreasing with increasing Y).

Table 3 .
An example of manually calculated background luminance for different room widths while the length also is changing.

Table 4 .
An example of simulated background luminance for different room widths while the length also is changing.

Table 5 .
The deviation of 2D-interpolated or extrapolated UGR average from the simulated UGR average for luminaires A to H.

Table 5 .
The deviation of 2D-interpolated or extrapolated UGR average from the simulated UGR average for luminaires A to H.

Table 6 .
The deviation of 2D-interpolated or extrapolated data average from average of the simulated UGR as a percentage of 30 units of UGR for luminaires A to H.

Table 7 .
The deviation of 2D-interpolated or extrapolated UGR table data from the simulated UGR as a percentage of 30 units of UGR for luminaires A to H.