Mobile QR Code QR CODE : Journal of the Korean Institute of Illuminating and Electrical Installation Engineers

Journal of the Korean Institute of Illuminating and Electrical Installation Engineers

ISO Journal TitleJ Korean Inst. IIIum. Electr. Install. Eng.

  1. (Ph.D. Student, Department of Electrical Engineering, Pukyong National University, Korea)



3-dimensional space vector, Inverse reflector design, Lambertian distribution, Simulation

1. Introduction

Fabricating an optical system is an inherently time-consuming and costly process, posing significant limitations on conducting experiments under diverse conditions. Consequently, computer-supported simulations have become increasingly critical in optical design[1]. Simulations enable performance verification and design optimization without the need for physical prototyping. However, optimization-focused simulations rely on iterative trial-and-error across extensive experimental ranges, resulting in substantial computational demands[2].

To address this challenge, various methods have been developed to improve simulation efficiency, including statistical approaches such as Design of Experiments (DOE)[3]. More recently, machine learning-based techniques have also been explored[4]. While these iterative simulations are frequently automated using simulator-provided application programming interfaces (APIs), the computational burden can become prohibitive especially when employing machine learning algorithms or handling complex systems with multiple components. In such cases, standard laboratory PCs are often inadequate, necessitating high-performance computational systems.

Providing appropriate initial values for optimization can significantly reduce computational costs and streamline the simulation process, particularly in resource-constrained environments. In this study, the basic shape of the optical system was derived using a straightforward approach based on three- dimensional spatial vectors. The initial design was then used as the starting point for optimization within the simulation. The efficiency of this method was validated by comparing the total number of simulations required when using the proposed approach against conventional trial-and-error and DOE methods.

A Study on Combining Simulation and 3D Space Vector for Efficient UV-C LED Reflector Design

2. Theoretical Background

2.1 Lambertian Distribution

In this study, a reflector was designed to redirect light emitted from an LED with a Lambertian distribution onto a rectangular target area(35㎝×35㎝) was examined, as shown in Fig. 1.

Fig. 1. Basic configuration of a reflector and LED source

../../Resources/kiiee/JIEIE.2025.39.3.167/fig1.png

A UV-C LED with a Lambertian distribution follows Lambert’s the cosine law[5], as given in Equation (1).

(1)
$I=\dfrac{\phi}{D^{2}}\cos\theta$

where $I$ is light intensity, $D$ is distance between

light source and the object, and $\phi$ is radiant flux.

As shown in Equation (1), the intensity of emitted light varies with the angle, and this variation is graphically illustrated in Fig. 2, where approximately 95% of the total light is confined within a range of ±60° (i.e., 120°). Therefore, under the assumption of a Lambertian distribution, it can be reasonably assumed that the LED emits light only within this ±60° range without significantly impacting the overall simulation results[6-8].

Fig. 2. Light distribution as a function of angle $\theta$

../../Resources/kiiee/JIEIE.2025.39.3.167/fig2.png

2.2 3D Vector Calculation

When considering a ray from the LED as a single line, the light reflection process between the source, reflector, and target can be represented using two straight vectors[9, 10], as illustrated in Fig. 3.

Fig. 3. Reflection process between source-reflector-target

../../Resources/kiiee/JIEIE.2025.39.3.167/fig3.png

Assume $P_{1}$ is a position of the light source, $P_{2}$ is the reflection point, and $P_{3}$ is a point in the target area. The reflection process can be expressed using two vectors going from $P_{1}$ → $P_{2}$ and $P_{3}$ → $P_{2}$ , respectively. Since the distance between a reflector and a target area (≥200mm) is significantly larger than the light source size (4mm), the light source was approximated as a point. If the former is vector $\vec{V_{1}}$ and the latter is vector $\vec{V_{2}}$ , using the vector expression method using parameters[11], $\vec{V_{1}}$ and $\vec{V_{2}}$ can be expressed as Equation (2) and (3).

(2)
$\vec{V_{1}}=\vec{r_{1}}+t\vec{u_{1}}$
(3)
$\vec{V_{2}}=\vec{r_{2}}+s\vec{u_{2}}$

$\vec{r_{1}}$ and $\vec{r_{2}}$ denote the vectors extending from the origin to point $P_{1}$ and $P_{2}$, respectively. Similarly, $\vec{v_{1}}$ and $\vec{v_{2}}$ represent the unit vectors corresponding to $\vec{V_{1}}$ and $\vec{V_{2}}$. The parameters $s$ and $t$ can be determined using the cross product of the vectors, as shown in Equation (4) and (5)[11].

(4)
$t=\dfrac{\det(\vec{r_{1}}-\vec{r_{2}},\: \vec{v_{1}},\: \vec{v_{1}}\times\vec{v_{2}})}{\Vert\vec{v_{1}}\times\vec{v_{2}}\Vert}$
(5)
$s=\dfrac{\det(\vec{r_{1}}-\vec{r_{2}},\: \vec{v_{2}},\: \vec{v_{1}}\times\vec{v_{2}})}{\Vert\vec{v_{1}}\times\vec{v_{2}}\Vert}$

Additionally, the vector that bisects the angle between two intersecting vectors can be obtained by summing their respective unit vectors, as presented in Equation (6)[11].

(6)
$\vec{v_{1}}+\vec{v_{2}}=\dfrac{\vec{V_{1}}}{\vert \vec{V_{1}}\vert}+\dfrac{\vec{V_{2}}}{\vert \vec{V_{2}}\vert}$

A unit vector has a magnitude of 1, and the sum of two unit vectors aligns with the diagonal of the rhombus formed by these vectors. Since the diagonal of a rhombus inherently bisects its angles, the sum of these unit vectors also acts as the angle bisector[12]. Consequently, the angle bisector can be computed as shown in Equation (6).

2.3 Material for the UV-C Band

The design method presented in Section 2.2 is mathematically derived from the general characteristics of visible LEDs. Since the reflector designed in this study is intended for use with LEDs operating in the UV-C band, the material choice for the reflector is critical.

The UV-C band, characterized by wavelengths ranging from 200 to 280nm[13], is known for the low reflectance with most materials[14]. As illustrated in Fig. 4, aluminum is the only material that exhibits high reflectance within the UV-C band. Therefore, aluminum is recommended as the optimal material for reflectors operating in this wavelength range. However, prior research on aluminum surfaces has shown that they lack smoothness, leading to diffuse reflection. So at points where the surface is uneven, the angle of incidence and the angle of reflection are not equal[15]. Thus, Snell’s Law of Reflection cannot be directly applied in such cases.

To design an accurate reflector, it is essential to account for light scattering on the aluminum surface. Mathematical models used to describe light scattering on rough surfaces include Kirchhoff Diffraction Theory and the Rayleigh Approach[16]. However, these methods are often difficult to intuitively understand and have significant limitations when applied to practical design scenarios. Alternatively, the Bidirectional Reflectance Distribution Function (BRDF) can be used to represent surface scattering. The BRDF characterizes the distribution of reflected light relative to incident light and is mathematically defined as a function of parameters such as the incident angle vector, the reflection point, and the reflection angle vector[17, 18]. In recent industrial applications, the BRDF is commonly measured directly, using specialized equipment.

The scattering characteristics of the reflector surface were modeled based on the BRDF properties of aluminum, as determined from previous research [19]. The detailed scattering characteristics are shown in Fig. 5.

Fig. 4. Comparison of reflectance by wavelength among each material[14]

../../Resources/kiiee/JIEIE.2025.39.3.167/fig4.png

Fig. 5. Measured BRDF of Al-on-steel[19]

../../Resources/kiiee/JIEIE.2025.39.3.167/fig5.png

3. Experiments and Results

3.1 Designing a Reflector

In Fig. 3, $P_{1}$ and $P_{3}$ are predetermined values, while $P_{2}$ must be determined through calculation. The coordinates of $P_{2}$ can be obtained using Equation (2) to (6), but multiple possible solutions for $P_{2}$ may exist. Hence, to calculate the value, several reference points were set, as shown in Fig. 6. Therefore, among the multiple possible values of $P_{2}$, it is necessary to determine the appropriate coordinates that enable the continuous construction of the reflector’s curve. Assuming a radiation angle of 120° for the LED, the light originating from $P_{1}$ reaches x-Start and x-End respectively, as illustrated in Fig. 6.

Points labeled $P_{2}'$ and $P_{2}''$ in Fig. 6. represent some of the possible solutions for $P_{2}$ near x-Start. Similar to Fig. 3, the bisectors $n_{1}$ and $n_{2}$ as well as their corresponding tangential lines can be drawn at these points. Since the tangential line at $P_{2}''$ intersects x-Start, $P_{2}''$ can be selected as the value for $P_{2}$.

By repeating this process to x-End, an entire line of points for constructing the reflector can be obtained. Subsequently, by iterating this procedure along the y-axis, the full surface coordinates of the reflector can be determined. These sequential steps were automated using Python, and the resulting reflector shape was visualized using the 3D modeling software Blender, as shown in Fig. 7.

Fig. 6. A schematic representation of the method for determining the value of $P_{2}$

../../Resources/kiiee/JIEIE.2025.39.3.167/fig6.png

Fig. 7. The shape of the reflector implemented through 3D modeling

../../Resources/kiiee/JIEIE.2025.39.3.167/fig7.png

3.2 Simulation Results

A simulation was conducted using LightTools software to evaluate the performance of the reflector designed in Section 3.1. The 3D-designed geometry was implemented in the simulation using the Freeform Entity function in LightTools. The simulation results are shown in Fig. 8.

The reflector material was specified as aluminum, and the BRDF characteristics shown in Fig. 5. were applied to account for surface scattering. Under these conditions, a simulation was conducted, and the resulting irradiance [W/m²] distribution in the target area is shown in Fig. 9(a).

The irradiance chart in Fig. 9(a) demonstrates that, while some light is concentrated at the center, the overall irradiance distribution is highly uniform. The target area was divided into a 10 × 10 grid, and the number of incident rays in each section was measured in the simulations. The uniformity $U$, as calculated using Equation (7), was found to be approximately 70%. Fig. 9(b) presents the measured results from a previous study[20], showing a uniformity of approximately 43%.

(7)
$U(\%)=\dfrac{\min . number \;of\; rays}{\max . number\; of \;rays}$

This study confirms that the designed reflector effectively fulfills its intended purpose. Without prior design considerations, directly fabricating the reflector’s shape would require numerous iterative experiments. For instance, the "Trough Parabola," a basic reflector included in LightTools, requires the specification of 14 design parameters. Table 1 summarizes these parameters, categorizing them by data type, range, and level.

For categorical data, the number of levels was set equal to the number of selectable items available in the simulation, while for numerical data, the levels were fixed at three. Under these conditions, performing a full factorial design of DOE would require a total of 8,748 runs[21]. Assuming each simulation takes 10seconds, the total simulation time would be approximately 24hours, with additional time potentially required, depending on the specifications of the computing resources. In contrast, the proposed method, which includes automated reflector design followed by minor adjustments in the simulation, required approximately 30minutes. This demonstrates that reflector design can be achieved more efficiently compared to traditional trial-and-error methods or statistical DOE approaches.

Fig. 8. The shape of the reflector implemented through simulation

../../Resources/kiiee/JIEIE.2025.39.3.167/fig8-1.png../../Resources/kiiee/JIEIE.2025.39.3.167/fig8-2.png

Fig. 9. (a) Measurement results from this study (b) Measurement results from previous study [20]

../../Resources/kiiee/JIEIE.2025.39.3.167/fig9.png

Table 1. Summary of experimental factors for the trough parabola

Factor

Data type

Range

Level

Half width

Numerical

data

0-50mm

3

Half height

0-50mm

x Eccentricity

0-1

y Eccentricity

0-1

Rotate

0-90˚

Radius

0-50mm

Offset

0-10mm

Shape

Categorical

data

Circle, Oval,

Square,

Rectangle

4

4. Conclusion

In this study, a reflector was designed using an inverse reflector design approach based on three-dimensional vector calculations to achieve uniform illumination of a target area from a UV-C LED with a Lambertian distribution. The performance of the designed reflector was evaluated through simulations, and the resulting irradiance measurements confirmed an overall uniform distribution, with minor deviations at the center.

The proposed method enables the derivation of a reflector’s basic shape through relatively simple calculations, and the design process can be efficiently automated using tools such as Python. Compared to traditional trial-and-error methods or full-factorial DOE designs without initial values, this approach offers significant time savings and improved efficiency.

Future research will focus on experimentally validating a reflecter fabricatedusing aluminum. Additionally, precise BRDF measurements tailored to the material properties will be conducted and incorporated into both the design and fabrication processes to enhance accuracy and effectiveness.

Acknowledgement

This work was supported by a Research Grant of Pukyong National University(2024).

References

1 
Ingo Sieber, Horst Eggert, Helmut Guth, and Wilfried Jakob, “Design simulation and optimization of micro-optical components,” SPIE, vol. 3430, 1998.DOI
2 
S. Andradottir, “A review of simulation optimization techniques,” in Proc. of 1998 winter simulation conference (Cat. No. 98CH36274), vol. 1, pp. 151-158, 1998.DOI
3 
J.P.C. Kleijnen, “Experimental design for sensitivity analysis, optimization and validation of simulation models,” Operations research, pp. 1-25, 1997.URL
4 
M. Luo, et al., “Efficient lens design enabled by a multilayer perceptron-based machine learning scheme,” Optik, vol. 273, 2023.DOI
5 
Herlina, Randa Adriano, and Rudy Setiabudy, “The Mapping of Lighting Intensity from the Light Distribution on LED and CFL Lamps,” in Proc. of 2019 4th ICITISEE, pp. 429-434, 2019.DOI
6 
Ron Lenk and Carol Lenk, “Practical lighting design with LEDs,” John Wiley & Sons, 2nd Ed., 2017.URL
7 
J. H. Lee, K. B. Nahm, J. H. Ko, and J. H. Kim, “Optimization of optical structure of lightguide panel for uniformity improvement of edge-lit backlight,” Korean Journal of Optics and Photonics, vol. 21, no. 2, pp. 61-68, 2010.DOI
8 
Bo Mei and Han Long, “Research on chaotic characteristics of array LED light system,” in Proc. of International Conference on Frontiers of Electronics, Information and Computation Technologies, vol. 121, pp. 1-6, 2021.DOI
9 
H. Wu, et al., “Design of a freeform reflector for LED sources based on a feedback method,” Lighting Research & Technology, vol. 48, no. 3, pp. 365-375, 2016.DOI
10 
L. B. Romijn, Jan H.M. ten Thije Boonkkamp, and W. L. IJzerman, “Inverse reflector design for a point source and far-field target,” Journal of Computational Physics, vol. 408, 2020.DOI
11 
Kreyszig Erwin, “Advanced engineering mathematics,” John Wiley & Sons, 8th Ed., 1998.URL
12 
S. K. Lee and C. G. Lee, “A study on equations of bisector and trisectors of angle,” Communications of Mathematical Education, vol. 21, no. 3, pp. 515-525, 2007.URL
13 
C. Wang, et al., “Efficiency improvement of batch reactors for water sterilization using UV-C LED arrays,” Environmental Technology, vol. 42, no. 25, pp. 4038-4046, 2021.DOI
14 
Y. X. Cui, “Study and development of near-infrared reflective and absorptive materials for energy saving application,” Ph.D. Thesis, Chemistry, Carleton University, 2012.URL
15 
Theodore V. Vorburger, Egon Marx, and Thomas R. Lettieri, “Regimes of surface roughness measurable with light scattering,” Applied optics, vol. 32, no. 19, pp. 3401-3408, 1993.DOI
16 
John C. Stover, “Optical scattering: Measurement and analysis,” SPIE Press, 3rd Ed., 1999.URL
17 
Rosana Montes and Carlos Ureña, “An overview of BRDF models,” University of Grenada, Technical Report LSI-2012-001, 2012.URL
18 
H. Zhang, Z. Wu, Y. Cao, and G. Zhang, “Measurement and statistical modeling of BRDF of various samples,” Optica applicata, vol. 40, no. 1, pp. 197-208, 2010.URL
19 
Maria Brogren, et al., “Optical properties, durability, and system aspects of a new aluminium-polymer-laminated steel reflector for solar concentrators,” Solar Energy Materials and Solar Cells, vol. 82, no. 3, pp. 387-412, 2004.DOI
20 
J. W. Ahn, K. H. Choi, and U. C. Ryu, “Designing bottom reflector for UV-C LED for HEPA filter sterilization,” in Proc. of 2024 spring Conference of KIIEE, pp. 50, 2024.URL
21 
S. H. Park, “Design of Experiments,” MinYoungSa, 2nd Ed., pp. 305-338, 2012.URL

Biography

Jae-Wook Ahn
../../Resources/kiiee/JIEIE.2025.39.3.167/au1.png

He received a B.S. degree in electronic engineering from Korea Aviation University, Korea, in 2004, and his M.S. degree in electrical engineering from Pukyong National University, Busan, Korea, in 2021. He is now a Ph.D. student in the Department of Electrical Engineering at Pukyong National University, Busan, Korea, and a government worker of ROKAF (Republic of Korea, Air Force).

Byeong-Seok Yu
../../Resources/kiiee/JIEIE.2025.39.3.167/au2.png

He received B.S., M.S., and Ph.D. degrees in Naval Architecture and Ocean Engineering from Seoul National University, Seoul, Korea, in 1990, 1993, and 1998, respectively. Until 2019, he worked for Hanjin Heavy Industries and Construction. Since Sep. 2019, he has been a professor in the Major of Smart Mechanics and Mobility Engineering at Pukyong National University, Busan, Korea.

Uh-Chan Ryu
../../Resources/kiiee/JIEIE.2025.39.3.167/au3.png

He received B.S. and M.S. degrees in electrical engineering from Seoul National University, Seoul, Korea, in 1997 and 1999, respectively, and his Ph.D. degree in the School of Information and Communications from Gwangju Institute of Science and Technology, Gwangju, Korea. Since Sep. 2013, he has been a professor in the Major of Electrical, Electronic, and Software Engineering at Pukyong National University, Busan, Korea.