31 December 2013 New approach to construct freeform surface by numerically differential formulation
Author Affiliations +
Abstract
We proposed a new design method for single freeform reflective (or refractive) surface tailored to redistribute the radiant flux onto a prescribed illumination pattern. Unlike the conventional optimization approaches based on the grid mapping, in this study we estimated each segmental freeform surface by locally solving a second-order differential equation, which formulates the energy transportation between each domain cell. With finite element method via Hermite element, we validated a series of smooth reflective/refractive surfaces to reallocate the radiant flux from a point source toward a target plane with specific patterns. The proposed technique offers a large flexibility by varying the vectors of each ray with multiple refraction (or reflection), which imposes no restriction on the target distribution, collective solid angle, or even target topography.

1.

Introduction

Construction 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 η=N·(×N)=0 is necessary in order to make a smooth freeform surface, where N 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).56.7.8.9 The MA PDE also appears in other research fields, such as differential geometry, optimal control and mass transportation, meteorology, and geostrophic fluid.1011.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 PDE

To 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 I (Watts/sr) of a point source (in polar domain Ωθ) and prescribed illuminance distribution E (Watts/m2) of the target region (in Cartesian domain Ωx), in one dimension, can be related as

(1)

ΩθI(θ)dθ=ΩxE(x)dx.

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 Ωx into many local subdomains Ωθi and Ωxi, so that each corresponding pair satisfies the equal flux i=1,2,,n, 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 []0 with basis vectors {e10,e20} and origin point O. The target plane T is expressed by a plane equation C0·[X]0d0=0, where C0=[c1c2]T, d0, and [X]0 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 i’th local coordinate of any vector V is denoted by [V]i with corresponding orthogonal basis {e1i,e2i} and local origin Oi. Coordinate transformation from i’th to j’th coordinate system can be conducted via a sequence of operators, ij=TjRij, where Rij and Tj represent the rotation of basis vectors and translation: Tj(V)=[V]j+[OiOj]j, respectively. By means of subsequent intermediate local coordinates, the light source O and the target plane T can be linked as [O]i=0i[O]0 and C0·i0[X]i=d0, respectively.

The scenario of our freeform construction is schematically in Fig. 1. The segmental freeform surface to be constructed on i’th local coordinate system can be analytically presented as a function [A]i=[u,f(u)] with surface normal as

(2)

[N]i=(fu,1)fu2+1.

Fig. 1

(a) Global coordinate system: origin (O) locates at light source with basis vectors {e10,e20}. Local i’th coordinate system: origin (Oi) locates at freeform surface with local basis vectors {e1i,e2i}.The angles θ and ϕ0i are the angle of OOi with OA and e2i, respectively. A 1-D freeform model assumes that the global coordinate and local coordinate are coplanar. (b) Partitioning Ωθ and Ωx into n subdomain Ωθi and Ωxi for i=1,2,,n is based on equal flux theory. Computation domain Ωui that unifying variable θ and x into local variable u is mapped from subdomain Ωθi and Ωxi.

OE_53_3_031307_f001.png

As shown in Fig. 1(a), while a ray from point source O is hitting on the point A at the freeform reflector, the ray is guaranteed to be transported toward the target plane T. The refraction law in vector form is

(3)

n2n1[XA]i[XA]i=[OA]i[OA]i[([OA]i[OA]i·[N]i)n22n12+n12([OA]i[OA]i·[N]i)2][N]i,
where n1 and n2 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 n1=1.0 and n2=1.0. Here, we size each cell of subdomain pair (Ωθi and Ωxi) to allow equal flux transfer. In order to derive the freeform surface via the preceding equations, we must introduce an intermediate variable u in the computational domain (Ωui) to link the computation in both source and target domains (Ωθi and Ωxi). By changing the variable, Eq. (3) can be rewritten as

(4)

ΩuiI(u)dθidudu=ΩuiE(u)dxidudu.

Obviously, the integral form of energy conservation, Eq. (1) can be represented by a differential equation in summation:7

(5)

EΩui·dxiduIΩui·dθidu=0foruΩui.

Let (uhi,0) be the coordinate of the hitting point where the incident ray OA with light source [O]i=(Lxi,Lzi) hits on i’th local x-axis. The x-coordinate of the hitting point uhi can be expressed as two separate functions of u and θi by trigonometric formula:

(6)

uhi(θi)=ρ0sin(θi)cos(θi+ϕ0i),

(7)

uhi(u)=f(u)(uLxi)f(u)Lzi+u,
where θi and ϕ0i are the angle of OOi with OA and e2i, respectively, and ρ0 is the distance between O and Oi. The transformation between variables u and θi can be linked by the joint point uh and the Jacobian chain rule becomes

(8)

dθidu=dθiduhiduhidu.

Similarly, suppose the incident ray OA is reflected onto the point [X]i=(xi,zi) at the target plane T. The variables xi and u related by the law of reflection can be expressed as

(9)

[OA]i[XA]i(uxi)=2([OA]i·[N]i)[f(u)](uLxi),

(10)

[OA]i[XA]i[f(u)zi]=2([OA]i·[N]i)[f(u)][f(u)Lzi].

Dividing Eqs. (9) by (10), the equation that defines a local freeform surface can be obtained:

(11)

uxif(u)zi=2([OA]i·[N]i)(uLxi)2([OA]i·[N]i)[f(u)][f(u)Lzi]i.

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 x-coordinate of the i’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 zi by xi through the equation C0·i0[X]i=d0:

(12)

xi(u)=C0if(u)+C1if(u)+C2iu+C3iC4if(u)+C5if(u)+C6iu+C7i+R(u,f),
where R(u,f)=O(u2,uf,uf,f2,ff,f2) is the high order nonlinear term that can be neglected under the small local planar approximation, and the coefficients, Cji, j=1,2,,7, are constants with respect to the i’th local coordinate system. For the refraction case, the variables xi and u can be related by the law of refraction following a similar derivation. After differentiating Eq. (12) to obtain Jacobian dxi/du and inserting Eqs. (8) and (12) into Eq. (5), we can devise the local reflection (refraction) freeform surface to a nonlinear differential equation:

(13)

α1if+α2if+α3if+α4i=0,
where α1i, α2i, α3i, and α4i are coefficient functions subject to low order derivatives of f. This equation describes the ray-tracing map and the energy redistribution. The total energy is conserved as ΩxiE(u)x(u)du in each local domain. Therefore, the energy deviation is mainly contributed by the neglected term R(u,f) and can be controlled when the integral ΩxiE(u)R(u,f)du 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 (Ωθi and Ωxi) 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 Ri1 and Ri in the polar cell toward the corresponding boundary [bi1,bi] in the Cartesian cell, so the radiant flux is fully transferred without stray loss. The edge ray,Ri, reflected by the freeform surface is propagating toward the corresponding point bi. We shall set a reasonable right end point of the intermediate cell Ωui, 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:

Algorithm

Freeform surface via local reconstruction (FS-LRA).

Input:Ωθ: Domain of light source: Ωx: Domain of target; O: Light source; T; Target plane equation; Global coordinate frame []0; O1: initial point of freeform surface:
Output: f(x): freeform surface function;
1. Partition the source and target domain into Ωθi and Ωxi for i=1,2,…,n;
2. Setup position and orientation of the freeform (i.e., determine global and the first local coordinate frame []1):
3. fori=1 to ndo
4.  Determine Ωui;
5.  Solve equation (13) on Ωui;
6.  Setup the i+1-th coordinate frame []i+1;
7.  Apply transformation operator Hii+1 to the light source and target plane;
8. end for

Fig. 2

Overview of the freeform construction process. Step 1: The source domain and target domain are divided into n subdomain pairs (Ωθi,Ωxi), i=1,2,,n. Step 2: The global information including the location of light source and the target plane is transformed into i’th coordinate. Step 3: The computational domain Ωui is determined by the corresponding subdomain pair (Ωθi,Ωxi), the refraction or reflection laws and the convexity of the freeform. Step 4: The freeform surface is calculated by solving Eq. (13) on Ωui using Hermite finite element method. Step 5: The (i+1)’th coordinate system is defined by the rightmost boundary point Oi+1 and the normal and tangent vectors at Oi+1. If i is not equal to n, go to step 2. Step 6: The freeform surface is assembled from local surfaces.

OE_53_3_031307_f002.png

3.

Freeform Surface via Local Reconstruction Algorithm: Error Evaluation

In 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 (x,y)=(0,10) of an ideal parabolic reflector y=x2/12+13. Here, we limit the emitted angle from the point source within the range of Ωθ=[119deg,119deg], and the radiant flux reflected by the ideal parabolic reflector was directed toward the given target plane Ωx=[10,10] 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 es=maxxΩxSFF(x)SP(x), maximum deviation between the constructed freeform surface (SFF) and the ideal parabolic curve (SP), is on the order of 104 to 106, 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.1617.18

Fig. 3

(a) S: point light source. The ideal parabolic curve is y=x2/12+13. (b) The maximum surface difference between ideal parabolic curve and numerical surface. The surface error with respect to subdomain numerical approximates a function S(x)1/x2. The amount of error reduces quadratically as we increase the number of partition cells.

OE_53_3_031307_f003.png

4.

Construction Algorithm: Two Design Types

4.1.

Uniform Illuminance

After 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 Ωθ=[30deg,30deg]. 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.

Fig. 4

(a) and (b) Layout of the freeform surface reflector and refractor is designed for uniform illuminance, respectively. (c) Illuminance distribution on the target plane. Inset is the magnification of illuminance plateau pattern. The nonuniformity of reflection case was well controlled within 1.05% RMS and 1.89% peak-to-valley. (d) The nonuniformity of refraction case was well controlled within 2.10% RMS and 4.01% peak-to-valley.

OE_53_3_031307_f004.png

4.2.

Linear Varying Illuminance

In 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 (x,y) along the x-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 x-direction. There are only about 5.58% and 1.52% RMS nonuniformity with the perfect chromatic line for reflection and refraction cases, respectively.

Fig. 5

(a) and (b) Layout of the freeform surface reflectors and refractor is designed for chromatic varying illuminance. The CIE chromatic of two point sources is red (x,y)=(0.73,0.27) and green (x,y)=(0.38,0.62). (c) and (d) Comparison of illuminance pattern by the constructed freeform surface (dash line) and that with ideally linear shape (solid line). The nonuniformity with the ideal linear case was well controlled within 5.58% RMS and 1.52% RMS in reflection case and refraction case, respectively. (e) and (f) Five sampling points along x-direction (10, 5, 0, 5, and 10 m) exhibit a linearly chromatic shift in CIE 1931 x, y chromaticity coordinate in both case.

OE_53_3_031307_f005.png

5.

Conclusion

In 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 x) in both domains (Ωθi and Ωxi) into a common intermediate variable u in the computational domain (Ωui). 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.

Acknowledgments

This 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.

References

1. W. A. Parkyn, “Segmented illumination lenses for step lighting and wall-washing,” Proc. SPIE 3779, 363–370 (1999).PSISDG0277-786X http://dx.doi.org/10.1117/12.368228 Google Scholar

2. L. WangK. QianY. Luo, “Discontinuous free-form lens design for prescribed irradiance,” Appl. Opt. 46(18), 3716–3723 (2007).APOPAI0003-6935 http://dx.doi.org/10.1364/AO.46.003716 Google Scholar

3. F. R. Fournieret al., “Fast freeform reflector generation using source-target maps,” Opt. Express 18(5), 5295–5304 (2010).OPEXFF1094-4087 http://dx.doi.org/10.1364/OE.18.005295 Google Scholar

4. H. RiesJ. A. Muschaweck, “Tailored freeform optical surfaces,” J. Opt. Soc. Am. A 19(3), 590–595 (2002).JOAOD60740-3232 http://dx.doi.org/10.1364/JOSAA.19.000590 Google Scholar

5. P. GuanX. J. Wang, “On a Monge-Ampere equation arising in geometric optics,” J. Differ. Geom. 48(2), 205–223 (1998).JDGEAS0022-040X Google Scholar

6. V. I. Oliker, “Mathematical aspects of design of beam shaping surfaces in geometrical optics,” in Trends in Nonlinear Analysis, pp. 191–222, Springer, Berlin Heidelberg (2002). Google Scholar

7. J. S. Schruben, “Formulation of a reflector-design problem for a lighting fixture,” J. Opt. Soc. Am. 62(12), 1498–1501 (1972).JOSAAH0030-3941 http://dx.doi.org/10.1364/JOSA.62.001498 Google Scholar

8. 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).INPEEY0266-5611 http://dx.doi.org/10.1088/0266-5611/14/3/016 Google Scholar

9. P. Benítezet al., “Simultaneous multiple surface optical design method in three dimensions,” Opt. Eng. 43(7), 1489–1502 (2004).OPEGAR0091-3286 http://dx.doi.org/10.1117/1.1752918 Google Scholar

10. S. Y. ChengS. T. Yau, “On the regularity of the Monge-Ampere equation det(2u/xixj)=f(x,u),” Commun. Pur. Appl. Math. 30(1), 41–68 (1977).CPMAMV0010-3640 http://dx.doi.org/10.1002/cpa.3160300104 Google Scholar

11. 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

12. R. J. McCannA. M. Oberman, “Exact semi-geostrophic flows in an elliptical ocean basin,” Nonlinearity 17(5), 1891–1922 (2004).NONLE50951-7715 http://dx.doi.org/10.1088/0951-7715/17/5/017 Google Scholar

13. 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

14. K. Wanget al., “Design of compact freeform lens for application specific light-emitting diode packaging,” Opt. Express 18(2), 413–425 (2010).OPEXFF1094-4087 http://dx.doi.org/10.1364/OE.18.000413 Google Scholar

15. Z. ZhenrongH. XiangL. Xu, “Freeform surface lens for LED uniform illumination,” Appl. Opt. 48(35), 6627–6634 (2009).APOPAI0003-6935 http://dx.doi.org/10.1364/AO.48.006627 Google Scholar

16. J. PeteraJ. F. T. Pittman, “Isoparametric Hermite elements,” Int. J. Numer. Methods Eng. 37(20), 3489–3519 (1994).IJNMBH0029-5981 http://dx.doi.org/10.1002/(ISSN)1097-0207 Google Scholar

17. 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).CMMECC0045-7825 http://dx.doi.org/10.1016/S0045-7825(98)00295-3 Google Scholar

18. 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).JSCOEB0885-7474 http://dx.doi.org/10.1007/s10915-008-9221-9 Google Scholar

Biography

Yu-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.

© The Authors. Published by SPIE under a Creative Commons Attribution 3.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.
Yu-Lin Tsai, Yu-Lin Tsai, Ming-Chen Chiang, Ming-Chen Chiang, Ray Chang, Ray Chang, Chung-Hao Tien, Chung-Hao Tien, Chin-Tien Wu, Chin-Tien Wu, } "New approach to construct freeform surface by numerically differential formulation," Optical Engineering 53(3), 031307 (31 December 2013). https://doi.org/10.1117/1.OE.53.3.031307 . Submission:
JOURNAL ARTICLE
6 PAGES


SHARE
RELATED CONTENT

Tailored Fresnel optics: I. Global shape
Proceedings of SPIE (November 11 2003)
Diffuse reflection from smooth dielectric surfaces
Proceedings of SPIE (November 30 1993)
Three-Dimensional Robot Vision
Proceedings of SPIE (February 15 1984)
Maximum gap of ideal concentrators for cylindrical absorbers
Proceedings of SPIE (October 05 1999)

Back to Top