|
1.IntroductionConstruction of a freeform surface to reallocate the radiant flux is a substantial challenge in terms of nonimaging optics. The basic concept is to convert the radiant flux emitted from a point source onto a prescribed illumination, usually in Cartesian coordinate. Parkyn1 firstly pioneered a tessellation method to map the cell area between a point source in polar grid with a particular rectangular illumination in Cartesian grid. The cells of both grid topologies were sized so that they share common energy. The methodology was based on the assumption that when the polar cell on the light source was exactly mapped to the corresponding rectangular target cell, the flux of the angular cell would be fully transferred to the target cell accordingly. Each corresponding pair of cells on these grids has an associated surface normal vector, and then governs the surface shape using extrinsic differential geometry. The major concern about equal flux grid method lies in a fact that the constructed surface, either reflective or refractive, is unlikely to guarantee the integrability condition in the overall surface, and thus leads to a discontinuous boundary. The cumulative deviations from the ideal surface normal become more pronounced at the far corners of the output grid. The discontinuous boundary could lead to a problem. The problem is that a small amount of radiant flux will be distorted by the sharp cliff. Hence, sufficient partition cells shall be numbered to confirm the estimated accuracy for surface construction. To alleviate the discontinuous issue of freeform, Wang et al.2 introduced the variable separation mapping strategy to take into account the normal deviation of the curve below a present value. The surface curve was regenerated with tangent vectors along the curve perpendicular to the normal vector, thus eliminating the accumulated normal deviation. The proposed scheme constructs a smooth surface by interpolating the discontinuity surface with non-uniform rational basis spline (NURBS).3 On the other hand, the integrability condition is necessary in order to make a smooth freeform surface, where is the surface normal.4 To fulfill the integrability condition, Fournier3 interpolated the ellipsoids, obtained from Oliker’s algorithm with appropriate grid partition, by using the NURBS where parameters of the NURBS were determined by minimizing the residual curl . Although the surface can be made smooth and the resulted radiant flux can be less distorted via the above mentioned approaches, the amount of energy deviation is still difficult to be estimated before the freeform surface is constructed. Furthermore, the computational costs of the constrained optimization for the NURBS’ parameters may rise significantly and the convergence is not guaranteed as the number of grid partition increases. In additional to the tessellation method, the freeform surface design problem is also related to the solution of the well-known Monge-Ampere (MA) equation. In 1972, Schruben first derived the MA equation obtained from reflector design. Beginning in 1980s until the present, Oliker et al. formulated the freeform problem into a highly nonlinear MA type of partial differential equation (PDE).5–9 The MA PDE also appears in other research fields, such as differential geometry, optimal control and mass transportation, meteorology, and geostrophic fluid.10–12 Although many works about existence and regularity of MA PDE have been conducted, numerical computation of freeform surface reconstruction by solving the MA PDE remains a challenge in practice. We are curious if there exists a systematic way to construct the freeform in which not only the surface is guaranteed to be smooth, but also the correspondence between light source and target plane is ensured under a controllable energy deviation. In this article, we formulate the ray forward and inverse propagation with negligible energy deviation within the MA PDE framework. By following the works of Caffarelli and Oliker13, the freeform considered here is assumed to be convex. The model is based on local energy conservation in which the flux correspondence associated with each local domain is prescribed. We take advantage of mathematical analysis of the MA PDE and successfully formulate the freeform design problem into the popular finite element paradigm. The formulation could maintain the ability to control the light beam and guarantee a smooth freeform surface. Moreover, this article describes this in detail by going through the different modeling steps of the model. All the main restrictions influencing the target distribution, collective solid angle, or source distribution are discussed, alongside the viability of realizing an optical surface for different specific requests. 2.Freeform Model in Localized PDETo develop a method for freeform construction subject to the source-target information, we must deduce an adequate, physically realistic model of ray propagation in the free space. The model can then be employed to understand the reflective (refractive) forward process of rendering illuminance according to a known optical configuration. Also, the model is applicable to the inverse algorithm of estimating the freeform surface from a given source-target corresponding condition. First of all, the ray transportation must satisfy the energy conservation among coordinate transformation.14,15 The radiant intensity () of a point source (in polar domain ) and prescribed illuminance distribution () of the target region (in Cartesian domain ), in one dimension, can be related as The goal is to design a freeform surface by analyzing the potential and limitation of using the PDE. In this article, the freeform surface is assumed to be continuous and convex. By partitioning the global domain and into many local subdomains and , so that each corresponding pair satisfies the equal flux , we are able to ensure the convergence of the integration. As long as the number of partitions is sufficiently large, the freeform can be considered to be a flat surface. As a result, a reduced but simplified representation of a freeform can be derived on each local domain. Hereafter, a flat surface criterion to the freeform in each local domain is called a small local planar approximation. The representation of the local freeform is indeed a nonlinear differential equation and shall be derived in the following. First, let us introduce some notations. Global coordinate frame is denoted by with basis vectors and origin point . The target plane is expressed by a plane equation , where , , and are the normal vector, shift and spatial variables, respectively. Although the laws of reflection (or refraction) are invariant with coordinate system, solution of freeform problem requires that the relations derived from these laws can be expressed in a coordinate system appropriate to the geometry of the given problem. Hence, we define a local coordinate system, by which the ’th local coordinate of any vector is denoted by with corresponding orthogonal basis and local origin . Coordinate transformation from ’th to ’th coordinate system can be conducted via a sequence of operators, , where and represent the rotation of basis vectors and translation: , respectively. By means of subsequent intermediate local coordinates, the light source and the target plane can be linked as and , respectively. The scenario of our freeform construction is schematically in Fig. 1. The segmental freeform surface to be constructed on ’th local coordinate system can be analytically presented as a function with surface normal as As shown in Fig. 1(a), while a ray from point source is hitting on the point at the freeform reflector, the ray is guaranteed to be transported toward the target plane . The refraction law in vector form is where and represent the refraction indexes near the source side and target plane side, respectively. For the reflection law,7 the refraction index can be defined by and . Here, we size each cell of subdomain pair ( and ) to allow equal flux transfer. In order to derive the freeform surface via the preceding equations, we must introduce an intermediate variable in the computational domain () to link the computation in both source and target domains ( and ). By changing the variable, Eq. (3) can be rewritten asObviously, the integral form of energy conservation, Eq. (1) can be represented by a differential equation in summation:7 Let be the coordinate of the hitting point where the incident ray with light source hits on ’th local -axis. The -coordinate of the hitting point can be expressed as two separate functions of and by trigonometric formula: where and are the angle of with and , respectively, and is the distance between and . The transformation between variables and can be linked by the joint point and the Jacobian chain rule becomesSimilarly, suppose the incident ray is reflected onto the point at the target plane . The variables and related by the law of reflection can be expressed as Dividing Eqs. (9) by (10), the equation that defines a local freeform surface can be obtained: If we fix the source and target plane position, the freeform surface defined in Eq. (11) can also be considered as a parameterized surface where the parameter is the -coordinate of the ’th coordinate system. The freeform parameters depend on the external geometry at which each optical route O-A-X constitutes an epipolar plane. In order to transform the underdetermined Eq. (11) to a determined problem, we replace by through the equation : where is the high order nonlinear term that can be neglected under the small local planar approximation, and the coefficients, , , are constants with respect to the ’th local coordinate system. For the refraction case, the variables and can be related by the law of refraction following a similar derivation. After differentiating Eq. (12) to obtain Jacobian and inserting Eqs. (8) and (12) into Eq. (5), we can devise the local reflection (refraction) freeform surface to a nonlinear differential equation: where , , , and are coefficient functions subject to low order derivatives of . This equation describes the ray-tracing map and the energy redistribution. The total energy is conserved as in each local domain. Therefore, the energy deviation is mainly contributed by the neglected term and can be controlled when the integral can be estimated.The construction of freeform surface by numerically differential formulation has an advantage in that we can regularize the smoothness tolerance by setting the partition number ( and ) in a tradeoff between computational efficiency and surface smoothness. In order to ensure the numerically estimated accuracy of freeform surface, sufficient partition cells shall be numbered in both domains. Here, we impose the edge rays and in the polar cell toward the corresponding boundary in the Cartesian cell, so the radiant flux is fully transferred without stray loss. The edge ray,, reflected by the freeform surface is propagating toward the corresponding point . We shall set a reasonable right end point of the intermediate cell , by which the freeform possesses convexity in the local domain. To solve the nonlinear PDE in Eq. (13), we employ the vanished moment method proposed by Feng and Neilan and use Newton’s method to find the adequate root function. The linearized equation at each iteration step can be solved by the finite element method. Here, we skip some straightforward but tedious expansion and rearrangement. Figure 2 demonstrates the key steps of the freeform construction process. The algorithm of freeform construction in this session can be summarized as following: AlgorithmFreeform surface via local reconstruction (FS-LRA).
3.Freeform Surface via Local Reconstruction Algorithm: Error EvaluationIn order to validate the proposed freeform scheme, we first constructed a surface by the previous algorithm and compared it with a known parabolic reflector; both are aimed to deliver a collimation beam as shown in Fig. 3. A point source was placed at the focal point of an ideal parabolic reflector . Here, we limit the emitted angle from the point source within the range of , and the radiant flux reflected by the ideal parabolic reflector was directed toward the given target plane with a uniform and parallel propagation. Figure 3(a) exhibits the surface comparison between the curve constructed via proposed algorithm and ideal parabolic curve. The surface error , maximum deviation between the constructed freeform surface () and the ideal parabolic curve (), is on the order of to , which is extremely small and quadratically inverse to the number of partition cells. The results are in agreement with the theoretical prediction that the error comes from the finite element method with Hermite element.16–18 4.Construction Algorithm: Two Design Types4.1.Uniform IlluminanceAfter validating our freeform algorithm with a comparative parabolic surface by forward testing, we now aim to achieve a uniform illuminance with a geometric correspondence by the freeform surface via local reconstruction algorithm (FS-LRA). As shown in Fig. 4, a 20-m-wide target plane was placed beneath a point source with 6-m distance. The emitted angle from the point source was confined for reflection and refraction within . With a mere 80 partition cells and 960 partition cells in both source-target subdomains, a freeform reflector and refractor can be generated in about 1.1 and 10 s, with 2.3 GHz Intel® dual core computer. The study employed a ray tracing simulation on the constructed freeform surface with a commercial tool LightTools®, in which light rays from the source were incidental to the target plane. The root-mean-square (RMS) nonuniformity is merely 1.05% (reflector) and 2.10% (refractor), respectively. It is noted that the accuracy or partition number of FS-LRA is highly dependent on the tolerance of the freeform surface. For uniform illuminance request in both reflect/refractive cases, a monotonic curve is accessible with the prescribed simple target surface. Because the mathematical form in the refraction case is more complex than that in reflection case, 960 partitions are required for convergence to achieve the design uniform illumination pattern, but only 80 partitions are required for the constructed reflector to achieve the same uniformity. 4.2.Linear Varying IlluminanceIn this example, we extended our technique to create a prescribed color illumination by a pair of independent sources with different colors. Figure 5 shows the schematic diagrams of the system layout for both reflect/refractive cases. With a 20-m target plane for both cases, we aimed to create a linearly varying international commission on illumination (CIE) chromatic along the -direction associated with two color sources. The distance between source and target plane is 10 and 70 m, respectively. The light sources leave the center at 10 m and 15 m, respectively. The different configuration conditions of reflection and refraction design are to achieve desired illuminance performance. As shown in Figs. 5(c) and 5(d), the illuminances feature a high linearity along the -direction. There are only about 5.58% and 1.52% RMS nonuniformity with the perfect chromatic line for reflection and refraction cases, respectively. 5.ConclusionIn this article, we proposed a freeform model based on local energy conservation associated with a series of coordinate transformation. The major advantage of this method is to exploit the well-proven existence and continuity of an elliptic nonlinear differential equation of the Monge-Ampere type. We successfully formulated the freeform design into a finite element paradigm with Newton’s iteration. The key maneuver is to unify the variables ( and ) in both domains ( and ) into a common intermediate variable in the computational domain (). After imposing the small local planar approximation, which is absolutely valid with a sufficient number of partition cells, a complete freeform reconstruction algorithm can be developed. In addition, the segmental differential formulation is capable of taking numerical error into account in each step. The proposed surface parameterization guarantees the surface smoothness with no restriction on the desired target distribution, collective solid angle, or source distribution, respectively. This technique still leaves many opportunities open and clearly more research must be carried out to explore its potential in full. First, at this moment, the issue was tackled in a one-dimensional case; complete treatment in consideration of twist deformation in freeform surface is underway. Second, point source approximation is another subject to be addressed when the freeform structure was placed in proximity of the light source such as light emitting diodes applications. The preliminary results presented in this article, however, indicate that parameterizing the freeform surface via local freeform PDE may create another route in the field of freeform optics, as it has the potential to take freeform design into many nonimaging applications. AcknowledgmentsThis work was financial supported by the National Science Council of Taiwan government under contract NSC101-2622-E-009-001-CC2 and 100-2115-M-009-001. The first author and second author contributed to this work equally. ReferencesW. A. Parkyn,
“Segmented illumination lenses for step lighting and wall-washing,”
Proc. SPIE, 3779 363
–370
(1999). http://dx.doi.org/10.1117/12.368228 PSISDG 0277-786X Google Scholar
L. WangK. QianY. Luo,
“Discontinuous free-form lens design for prescribed irradiance,”
Appl. Opt., 46
(18), 3716
–3723
(2007). http://dx.doi.org/10.1364/AO.46.003716 APOPAI 0003-6935 Google Scholar
F. R. Fournieret al.,
“Fast freeform reflector generation using source-target maps,”
Opt. Express, 18
(5), 5295
–5304
(2010). http://dx.doi.org/10.1364/OE.18.005295 OPEXFF 1094-4087 Google Scholar
H. RiesJ. A. Muschaweck,
“Tailored freeform optical surfaces,”
J. Opt. Soc. Am. A, 19
(3), 590
–595
(2002). http://dx.doi.org/10.1364/JOSAA.19.000590 JOAOD6 0740-3232 Google Scholar
P. GuanX. J. Wang,
“On a Monge-Ampere equation arising in geometric optics,”
J. Differ. Geom., 48
(2), 205
–223
(1998). JDGEAS 0022-040X Google Scholar
V. I. Oliker,
“Mathematical aspects of design of beam shaping surfaces in geometrical optics,”
Trends in Nonlinear Analysis, 191
–222 Springer, Berlin Heidelberg
(2002). Google Scholar
J. S. Schruben,
“Formulation of a reflector-design problem for a lighting fixture,”
J. Opt. Soc. Am., 62
(12), 1498
–1501
(1972). http://dx.doi.org/10.1364/JOSA.62.001498 JOSAAH 0030-3941 Google Scholar
S. A. KochenginV. I. OlikerO. von Tempski,
“On the design of reflectors with prespecified distribution of virtual sources and intensities,”
Inverse Probl., 14
(3), 661
–677
(1998). http://dx.doi.org/10.1088/0266-5611/14/3/016 INPEEY 0266-5611 Google Scholar
P. Benítezet al.,
“Simultaneous multiple surface optical design method in three dimensions,”
Opt. Eng., 43
(7), 1489
–1502
(2004). http://dx.doi.org/10.1117/1.1752918 OPEGAR 0091-3286 Google Scholar
S. Y. ChengS. T. Yau,
“On the regularity of the Monge-Ampere equation ,”
Commun. Pur. Appl. Math., 30
(1), 41
–68
(1977). http://dx.doi.org/10.1002/cpa.3160300104 CPMAMV 0010-3640 Google Scholar
T. GlimmV. I. Oliker,
“Optical design of single reflector systems and the Monge–Kantorovich mass transfer problem,”
J. Math. Sci., 117
(3), 4096
–4108
(2003). http://dx.doi.org/10.1023/A:1024856201493 Google Scholar
R. J. McCannA. M. Oberman,
“Exact semi-geostrophic flows in an elliptical ocean basin,”
Nonlinearity, 17
(5), 1891
–1922
(2004). http://dx.doi.org/10.1088/0951-7715/17/5/017 NONLE5 0951-7715 Google Scholar
L. A. CaffarelliV. I. Oliker,
“Weak solutions of one inverse problem in geometric optics,”
J. Math. Sci., 154
(1), 39
–49
(2008). http://dx.doi.org/10.1007/s10958-008-9152-x Google Scholar
K. Wanget al.,
“Design of compact freeform lens for application specific light-emitting diode packaging,”
Opt. Express, 18
(2), 413
–425
(2010). http://dx.doi.org/10.1364/OE.18.000413 OPEXFF 1094-4087 Google Scholar
Z. ZhenrongH. XiangL. Xu,
“Freeform surface lens for LED uniform illumination,”
Appl. Opt., 48
(35), 6627
–6634
(2009). http://dx.doi.org/10.1364/AO.48.006627 APOPAI 0003-6935 Google Scholar
J. PeteraJ. F. T. Pittman,
“Isoparametric Hermite elements,”
Int. J. Numer. Methods Eng., 37
(20), 3489
–3519
(1994). http://dx.doi.org/10.1002/(ISSN)1097-0207 IJNMBH 0029-5981 Google Scholar
K. H. Muci-KuchlerJ. C. Miranda-ValenzuelaT. J. Rudolphi,
“Formulation of a new error indicator for adaptive meshing with Hermite boundary elements,”
Comput. Methods Appl. Mech. Eng., 173
(3–4), 419
–431
(1999). http://dx.doi.org/10.1016/S0045-7825(98)00295-3 CMMECC 0045-7825 Google Scholar
X. FengM. Neilan,
“Vanishing moment method and moment solutions for fully nonlinear second order partial dierential equations,”
J. Sci. Comput., 38
(1), 74
–98
(2009). http://dx.doi.org/10.1007/s10915-008-9221-9 JSCOEB 0885-7474 Google Scholar
BiographyYu-Lin Tsai is a PhD candidate in the Department of Photonics at National Chiao Tung University. He received his MA degree in the Institute of Mathematical Modeling and Scientific Computing from National Chiao Tung University, Hsinchu, Taiwan, in 2010. His research interests include nonimaging optical design and freeform design problems. Ming-Chen Chiang is an engineer in the Department of Wireless Product of Universal Scientific Industrial. He received his MA degrees both in IMMSC (Institute of Mathematical Modeling and Scientific Computing) and Institute of Communication Engineering from National Chiao Tung University, Hsinchu, Taiwan, in 2011. His interests include the freeform design problems and communication channel analysis. Ray Chang is a master student in the Institute of Display at National Chiao Tung University. He received his MA degree from the Institute of Display from National Chiao Tung University, Hsinchu, Taiwan, in 2013. His research includes nonimaging freeform design problems. Chung-Hao Tien is an associate professor in the Department of Photonics at the National Chiao Tung University (NCTU) where he has been a faculty member since 2004. He received his PhD degree in the Institute of Electro-Optical Engineering from NCTU in 2004. His research interests include computational imaging, free form nonimaging, display and\ lighting optics. Chin-Tien Wu is an associate professor in the Department of Applied Mathematics at National Chiao Tung University. He received his PhD degree in Scientific Computation and Mathematical Modeling from the University of Maryland in 2003. His research interests include scientific computing and freeform design problems. |