LED diffused transmission freeform surface design for uniform illumination

LED diffuse transmission has been widely applied in backlight display area. However, the application of diffuse transmission is mostly based on experimental method and lack of theoretical analysis. Based on the above problems, an approximate numerical method is proposed to design diffuse transmission freeform surfaces. First, based on the transmission surface and the property of the LED, a mathematical model is established according to the property of LED, and the nonlinear equation of the illuminance distribution on the diffuse transmission inner surface is deduced. Then, by solving the equation, the profiles of diffuse transmission freeform surface for an indirect illumination system are obtained. Finally, the non-sequential ray tracing is used to detect the performance of the illuminance distribution on the target plane. Experimental results show that the diffuse transmission freeform surface has a better lighting performance than traditional diffuse transmission lighting systems.


Introduction
LED has many advantages such as high efficiency, energy saving, environmental protection, long life and small volume [1][2][3]. It is rapidly replacing traditional light sources in many applications. However, LED has the disadvantages of too large divergence angle and uneven distribution of illuminance, making it difficult to directly apply to lighting [2]. In order to improve the lighting performance in the effective area of the target plane, it is usually necessary to perform a secondary optical design to redistribute LED illuminance pattern [4,5].
In recent years, many scholars have conducted secondary optical design of LEDs. Moreno had put forward an optical design method of direct lighting [6,7]. Compared with direct lighting methods, redistribution of emitted light by catadioptric is an effective way to improve illuminance performance [8][9][10][11][12][13]. However, the traditional reflector or refractive lens has a high requirement for the processing precision and the material. In addition, in many cases, we need to use diffuse transmission to improve the uniformity of illumination system such as industrial inspections and/or real life [14,15]. Ripoll J proposed the expression of the reflection and transmission coefficients at the diffusiondiffusion interface [16]. M. Pang proposed a novel measurement method for obtaining a large laser beam intensity distribution based on CCD diffuse transmission imaging [17]. Serra et al. [18] investigate the CW transmitted light profiles when an inhomogeneity is immersed in highly scattering media slabs. Wu et al. [19] through the analysis of single-layer and multilayer diffuse transmission plane influence on the luminous intensity distribution concluded that increasing the distance between the first diffuse screen and the target plane can improve the spot uniformity. Bressel and Reich [20] applied the Mie theory to calculate the scattering properties of the material. These scholars studied the characteristics of the illuminance distribution through the diffuse transmission plane. In the aspect of diffuse applications, Kim et al. [21] presented an integrated prism with a bridge pattern of quantum dot (QD) film to improve the efficiency of the LED backlight unit. Qin [22] enhanced the brightness of backlight by solving the related characteristics of prism film angle. Wang et al. [23] compared the distance-height ratio of different LED lighting systems and proposed a new reversing design method, and a kind of ultrathin direct-lit LED backlight system is designed. The above research on diffuse transmission has achieved desired results, but their design is only based on experiments or direct applications of diffuse transmission, with no direct design of diffuse surface to enhance performance of the whole system.
A mathematical model for diffuse transmission is presented in this paper, and the diffuse transmission surface is designed, which improved the illumination uniformity in the specific area of the target plane. Firstly, according to the illumination requirements of the transmission surface and the target plane, a series of nonlinear functions on profiles of the diffuse transmission freeform surface are established. Then, solve the equation numerically to get a set of point coordinates that describe the contour of the surface. Using the MATLAB data interaction and optical software TracePro, a solid model of the light source is obtained. Finally, the illumination performance of several lighting systems is compared through experiments, and the experimental results testify the validity of the design method.

The design of diffuse transmission of freeform surface
In practice, this dependence turns out to be a power law that mostly depends on the encapsulant and semiconductor region shapes. The irradiance distribution is given by [6,7] E is the illuminance distribution on the light surface. E 0 is the distribution of illuminance on the freeform surface of LED light source when the \; is 0°. \; is the angle between the light from the LED and the direction of the vertical freeform surface. m, as Eq. (2), is generally provided by the LED manufacturer.
where / 1/2 is the angle between the direction of the normal direction and the 1/2 direction of the maximum luminous intensity, that is, the angle corresponding to half of the maximum luminous intensity value. In this article, the m value is 1 because of the most common Lambert body. The luminous intensity of the light emitted by the LED source in each direction of the receiver is I ; ð Þ as Eq. (3).
I 0 is the luminous intensity of the LED light source on the freeform surface when the \; is 0°. According to the inverse square law of illuminance, as shown in Fig. 1, the illuminance produced by a LED on the receiving surface is like Eq. (4).
The R is the vertical distance between the LED and the target plane. From Eq. (4), we can get the illuminance distribution at any point on the freeform surface, and its illuminance is obtained from Eq. (5).
E represents the illuminance distribution of any point in the target plane. Figure 2a, b shows the Cartesian rectangular coordinate system. Substitute Eq. (1) into Eq. (6).
E i x; y; z ð Þ represents the illuminance of arbitrary coordinate positions of diffuse transmission freeform surface.
Þ represents the coordinates on the freeform surface, x s ; y s ð Þ represents the position of the LED light source, and n represents the number of discrete points on the freeform surface. The illuminance generated by all LED array light sources is superimposed on the diffuse transmission freeform surface.
To simplify the calculation, a number of sample points (x 1 , 0, z 1 ), (x 2 , 0, z 2 ), …, (x u , 0, z u ) which can be regarded as Lambertian secondary point sources are selected to approximate the outline of the freeform surface. Suppose that an outgoing light L i of one LED x i ; y i ; 0 ð Þis L i x i ; 0; z i ð Þ on the transmission surface. As the Lambertian secondary point sources on the freeform surface, the illuminance of L i on the target plane is obtained from Eq. (8). where z j is the distance between LED and the target plane, ; j is the viewing angle and E 0j is the illuminance distribution at a z j distance and \; j ¼ 0 . As shown in Fig. 3, we establish the mathematical model of the surface profile of the diffuse transmission freeform surface. Considering that the diffuse transmission freeform surface model we have established is a rotational symmetry model, it is only necessary to design a curve in the freeform in this paper. Therefore, we simplify a three-dimensional problem to a two-dimensional problem; that is, only one of the contour lines in the freeform surface can be obtained through rotation to obtain the whole freeform surface model. Equation (8) is transformed into Cartesian coordinate system and then converted to Eq. (9).
u is the number of sampling points of the freeform surface in the x-axis direction. Through ray tracing, the vector from point x i ; 0; z i ð Þ to point x j ; 0; R À Á is emergent ray x j is the coordinates of the target plane along the x-axis, and R is a constant. The normal vector on the diffuse transmission freeform surface is Here dz and dx represent the differential of z and x, respectively. Assume that the included angle between OUT ! and NOR ! is \o, and Eq. (10) is obtained.
The ideal diffuse transmission freeform surface with the cosine characteristic can be regarded as a continuous distributed point light source, and its illuminance on the target plane is obtained from Eq. (11). The transmission light emitted from this freeform inner surface is scattered equally at random direction.
where I 0j is the central luminous intensity corresponding to \o ¼ 0 , the whole surface illuminance distribution as in Eq. (12).
where ''a'' is the length of the outline of the freeform surface in the x-axis direction and r s is the distance between the freeform surface to the target plane. da represents the area element on the freeform surface. BTDF is a bidirectional transmittance distribution function, which is used to represent the transmission characteristics of the transmission surface.
The Target Surface where M is the radiation emission, Q represents the radiant flux, L represents the radiance, / represents the solid angle and k represents the transmittance of the diffuse transmission material. It is not difficult to derive from Eq. (13) that the expression of BTDF is  The Eq. (12) is integrated on the target plane, and the expression of illuminance on the whole target plane is obtained, such as Eq. (15).
where r s is the length of the target plane in the x-axis direction. As in Fig. 3, a number of feature points are selected on the target lighting surface to calculate the structure parameters that can make the d 2 E s =dx 2 ¼ 0 get minimal value at each feature point. In this way, the coordinates of every point q i on the two-dimensional curve on the freeform surface can be calculated by numerical integration iteration. According to the discrete data obtained from MATLAB programming, it is fitted to a closed curve. As shown in Fig. 4, in order to prove the validity of this design, this paper selects four kinds of LED arrays, namely linear array of 2 LEDs, square LED array, hexagonal LED array and circular array. The coordinate position of the LED in the 2D (x, y) and the unit of axis x and axis y were both in millimeters:  (14.14, 14.14), (-14.14, 14.14), (14.14, -14.14), (-14.14, -14.14). Figure 5 shows the diffuse transmission freeform surface 2D contours required for the four LED arrays. Then the complete freeform surfaces are obtained by rotating 360°around the z-axis as shown in Fig. 6. Here (a), (b), (c) and (d) represent four different diffuse transmission freeform surfaces of 2-LED array, square array, hexagonal array and circular array.

Results
As shown in Fig. 5, based on the above design method of diffuse transmission freeform surface, we draw the surface profile of the diffuse transmission freeform surface and use the Monte Carlo ray tracing principle of TracePro software to simulate the experimental results. Select the LED with a power of 1 W. For the convenience of experiment, we set each LED as point source, and the emergent ray is 250,000. In this paper, the distance between LED and target plane is 600-1500 mm, with 100 mm to 500 mm as the effective radius of illumination area.
In order to prove the validity of the design of the diffuse transmission freeform surface, the experimental results of the diffuse transmission freeform surface are compared with the common diffuse transmission hemispherical surface, the diffuse transmission sheet and the direct illumination. Figure 7 shows the spot of light under a diffuse transmission freeform surface. Table 1 shows the illumination uniformity of the four lighting systems. Figure 8 shows the uniformity and efficiency of the illumination of a circular LED array under different conditions. It is obvious that with the increase in the distance between the light source and the target plane, the lighting conditions of each system have changed. In this paper, the uniformity of the designed diffuse transmission freeform surface is higher than that of other lighting systems. Table 1 shows that the uniformity of the designed freeform illumination system is about 3% higher than other illumination systems, which is an improvement for the diffuse transmission illumination system. Table 2 shows the lighting efficiency of the four lighting systems in the near field. Compared with other diffuse transmission surfaces, the efficiency of the transmission hemispherical surface is significantly lower than that of the other three lighting systems. In addition, the diffuse freeform surface is slightly better than the direct exposure but is slightly worse than the transmission sheet. In this regard, our explanation is that two times of the transmission of the LED rays from the freeform surface to the target plane have taken place so that the efficiency of the lighting is slightly reduced.

Conclusion
In this paper, a method of transmission freeform surface design for LED is presented. Based on several LED arrays and diffuse transmission characteristics, a set of nonlinear equations for solving the contour of freeform surface is derived by establishing an approximate mathematical model in a Cartesian coordinate system. The illuminance distribution of the diffuse transmission freeform surfaces are compared with direct illumination, transmission sheet and transmission hemispherical surface on the effective illumination area. Experimental results show that the illumination performance of the diffuse transmission freeform surface lighting system is superior to others. The rationality and effectiveness of this design are further verified. This method proposed in this paper also expanded the prospect of the indirect lightings in practical.