Designing optical free-form surfaces for extended ... - OSA Publishing

Designing optical free-form surfaces for extended sources ∗

2 A. V¨ ¨ R. Wester,1, G. Muller, oll,2 M. Berens,2 J. Stollenwerk,1,2 and P. Loosen1,2 2

1 Fraunhofer Institute for Laser Technology ILT, Aachen, Germany Chair for Technology of Optical Systems, RWTH Aachen University, Germany ∗ [email protected]

Abstract: LED lighting has been a strongly growing field for the last decade. The outstanding features of LED, like compactness and low operating temperature take the control of light distributions to a new level. Key for this is the development of sophisticated optical elements that distribute the light as intended. The optics design method known as tailoring relies on the point source assumption. This assumption holds as long as the optical element is large compared to the LED chip. With chip sizes of 1 mm2 this is of no concern if each chip is endowed with its own optic. To increase the power of a luminaire, LED chips are arranged to form light engines that reach several cm in diameter. In order to save costs and space it is often desirable to use a single optical element for the light engine. At the same time the scale of the optics must not be increased in order to trivially keep the point source assumption valid. For such design tasks point source algorithms are of limited usefulness. New methods that take into account the extent of the light source have to be developed. We present two such extended source methods. The first method iteratively adapts the target light distribution that is fed into a points source method while the second method employs a full phase space description of the optical system. © 2014 Optical Society of America OCIS codes: (080.4298) Nonimaging optics; (220.2740) Geometric optical design.

References and links 1. S. J. Schruben, “Formulation of a Reflector-Design Problem for a Lighting Fixture,” J. Opt. Soc. Am. 62, 1498– 1501 (1972). 2. H. Ries and J. Muschaweck, “Tailored freeform optical surfaces,” J. Opt. Soc. Am. A 19, 590–595 (2002). 3. B. Parkyn and D. Pelka, “Free-form illumination lenses designed by a pseudo-rectangular lawnmower algorithm,” Proc. SPIE 6338, 633808–633808–7 (2006). 4. A. Bruneton, A. Bäuerle, M. Traub, R. Wester, and P. Loosen, “Irradiance tailoring with two-sided, Fresnel-type freeform optics,” Proc. SPIE 8485, 84850H–84850H–7 (2012). 5. A. Bäuerle, A. Bruneton, R. Wester, J. Stollenwerk, and P. Loosen, “Algorithm for irradiance tailoring using multiple freeform optical surfaces,” Opt. Express 20, 14477–14485 (2012). 6. F. R. Fournier, W. J. Cassarly, and J. P. Rolland, “Optimization of single reflectors for extended sources,” Proc. of SPIE 7103, 71030I–71030I–10 (2008). 7. P. Ben´ıtez, J.-C. Miñano, J. Blen, R. Mohedano, J. Chaves, O. Dross, M. Hernández, and W. Falicoff, “Simultaneous multiple surface optical design method in three dimensions,” SPIE Proc. Optical Engineering 43, 1489–1502 (2004). 8. A. Rabl and J. Gordon, “Reflector design for illumination with extended sources: the basic solutions,” Applied Optics 33, 6012–6021 (1994). 9. T. Glimm and V. Oliker, “Optical design of single reflector systems and the mongekantorovich mass transfer problem,” Journal of Mathematical Sciences 117, 4096–4108 (2003).

#202374 - $15.00 USD Received 24 Dec 2013; revised 27 Feb 2014; accepted 27 Feb 2014; published 7 Mar 2014 (C) 2014 OSA 10 March 2014 | Vol. 22, No. S2 | DOI:10.1364/OE.22.00A552 | OPTICS EXPRESS A552

10. S. Seroka and S. Sertl, “Modeling of refractive freeform surfaces by a nonlinear PDE for the generation of a given target light distribution,” International Light Simulation Symposium (2012). 11. R. Wu, L. Xu, P. Liu, Y. Zhang, Z. Zheng, H. Li, and X. Liu, “Freeform illumination design: a nonlinear boundary problem for the elliptic Monge-Ampére equation.” Optics letters 38, 229–231 (2013). 12. S. A. Kochengin and V. I. Oliker, “Determination of reflector surfaces from near-field scattering dataII. Numerical solution,” Numer. Math. 79, 553–568 (1998). 13. F. R. Fournier, W. J. Cassarly, and J. P. Roland, “Designing freeform reflectors for extended sources,” Proc. of SPIE 7423, 743202–743202–12 (2009). 14. S. Haker, L. Zhu, A. Tannenbaum, and S. Angenent, “Optimal Mass Transport for Registration and Warping,” Int. J. of Comput. Vision 60, 225–240 (2004). 15. A. Bruneton, A. Bäuerle, R. Wester, J. Stollenwerk, and P. Loosen, “High resolution irradiance tailoring using multiple freeform surfaces,” Opt. Express 21, 10563–10571 (2013). 16. J. Bortz and N. Shatz, “Iterative generalized functional method of nonimaging optical design,” Proc. SPIE 6670, 66700A–66700A–15 (2007). 17. R. Wester, A. Bruneton, A. Bäuerle, J. Stollenwerk, and P. Loosen, “Irradiance tailoring for extended sources using a point-source freeform design algorithm,” Proc. of SPIE, Optical Systems Design 8550, 85502S–85502S– 8 (2012).

1.

Introduction

Driven by the rapid development of high power light emitting diodes (LED) as light sources, the need for lenses and mirrors that can distribute the light from LEDs in a predetermined manner has strongly increased. The small size of single chip LEDs enables the realization of very general light distributions at least as long as the size of the optics is large compared to the LED size so that the point source approximation remains valid. Design algorithms that rely on the point source assumption can be considered to be state of the art ( [1–5]), although the numerical realization still is a challenging task. The success of LED light sources triggered the demand for ever higher luminous fluxes in order to enter markets that today are still dominated by conventional light sources. The power output per chip of LED is limited so that multi-chip modules are being developed that now are entering the market with a luminous flux of up to 5000 Lumen. The size of such modules amounts to a few cm2 . In most cases customers prefer small-sized lighting systems. This forbids scaling the light-forming elements along with the new LED module sizes. Optical free-form design algorithms thus have to account for the extent of light sources. There have been some enhancements of point source algorithms to account for extended sources. A distribution that is fed into a point source algorithm is interactively adapted so that the distribution that is generated by the resulting optical surface matches the desired distribution as closely as possible. Besides the adaptation of the target distribution, Fournier et.al. [6] optimize the position of a point source in order to get the best results for a stationary extended source. The SMS method developed by Ben´ıtez, Miñano et.al. [7] allows for the simultaneous construction of two optical surfaces that exactly transform two input congruences (zero e´ tendue wave fields) into two prescribed output congruences. With more surfaces more congruences can be transformed. SMS can thus best be described as a multi-étendue-zero-source method. The inclusion of source non-uniformities is not easily possible. The method was originally designed for the 2D case but was extended to approximately handle 3D cases as well. Rabl on the other hand [8] employs the edge ray principle to construct 2D freeform reflector surfaces. Source non-uniformities cannot be taken into account. In the following section we supply some definitions. In section 2 we give a short overview over existing point source algorithms and in section 3 we outline a new approach to improve optical design for extended sources utilizing phase space concepts.


1.1.

Parametrization of optical surfaces

Optical surfaces have to be described or parametrized in some way. Table (1) shows a list of surface parametrization schemes used in the field of optical design methods. The design task consists in determining the parameters of the optical surfaces in an optical system. In the case of free-form surfaces the number of parameters can exceed 105 . Table 1. Optical surface parametrizations. surface representation conic surfaces polynomials radial basis functions Bezier-patches B-Splines triangulated surfaces

1.2.

parameters curvature, conic coefficient polynomial coefficients coefficients, centers control points control points nodes

Light flux conservation

In lossless optical systems the differential fluxes dΦ are equal on the source and target side, i.e.: ns Ls (~us ,~ss )dΩs cos(θs )dAs = nt Lt (~ut ,~st )dΩt cos(θt )dAt

(1)

Index s designates the source side and index t the target side respectively. L(~u,~s) is the radiance (refer to Table (2)). With the assumption of e´ tendue conservation Table 2. Table of symbols used in this article. ~u ~s cos(θ ) ~N dA ns nt

-

position vector on the surface ray unit direction vector direction cosine (~s · ~N) surface normal area element source-side refractive index target-side refractive index

ns dΩs cos(θs )dAs =1 nt dΩt cos(θt )dAt

(2)

Ls (~us ,~ss ) = Lt (~ut ,~st )

(3)

Equation (1) yields

i.e. the radiance is constant along a ray path. 2.

Point source algorithms

A lot of design algorithms have been published that exploit the point source assumption or assumption of zero source extent. Table 3 gives an overview over the principal methods. Direct


Table 3. Principal point source algorithms. Method Directly optimizing a merit function employing Monte-Carlo ray tracing. Ray targeting, i.e. ray mapping with subsequent surface reconstruction. Method of ellipsoids. Determine surface such that ∆E = E(prescribed) − E(computed) = 0. Explicit formulation and solution of a Monge-Ampère type equation.

Authors

Parkyn [3], Bruneton et.al. [4] Oliker and Glimm [9] Ries and Muschaweck [2], Bäuerle et.al. [5] Seroka and Sertl [10], Wu et.al. [11]

optimization with Monte-Carlo ray tracing has only limited applicability in the field of freeform optics design. The ray mapping methods suffer from the problem of finding ray mappings that meet the integrability condition of the optical surface. In the case of one-dimensional flux density distributions this is of no concern but in the general two dimensional case this can be a severe limitation. The method of ellipsoids is mathematically sound 2some drawbacks. 4 but has K K 0 log τ [12], with The resulting surface is only C and the computation time is O γ γ K the number of ellipsoids and γ and τ being two characteristic numbers of the problem. This implies a steep increase of the computation time with the number of ellipsoids, so this method is best suited for small numbers of ellipsoids (e.g. [13] employ 100 ellipsoids). The last two methods in the above list are the most advanced point source algorithms that satisfy the integrability condition automatically and can be implemented efficiently. A drawback of these methods is that highly nonlinear equations have to be solved which requires starting configurations that are close enough to the solution in order to converge. Fig. 1 shows the general set up of a point source, an optical element with free-form surfaces and a target plane.

Fig. 1. The rays emitted by a point source travel through the optical system to the target plane. The optical system is designed so that ray paths don’t cross.


Flux conservation yields: Is (sx , sy )dsx dsy = Et (ux , uy )dux duy

(4)

In Eq. (4) sx , sy are the direction vector components parallel to the source surface, ux , uy are the Cartesian coordinates in the target plane, Is (sx , sy ) is the direction dependent intensity of the source and Et (ux , uy ) the irradiance on the target plane. 2.1.

Monge-Ampère transport problem

The surface design task is to find a bijective coordinate mapping: ux = ux (sx , sx ) uy = ux (sx , sy )

(5) (6)

that realizes the prescribed irradiance distribution on the target. Inserting into Eq. (4) yields: Is (sx , sy )dsx dsy = Et (ux (sx , sy ), uy (sx , sy ))det(J)dsx dsy with the determinant of the Jacobian matrix J: ∂ ux ∂ uy ∂ ux ∂ uy − det(J) = ∂ sx ∂ sy ∂ sy ∂ sx

(7)

(8)

Under certain conditions the mapping function is given by the gradient of a convex scalar function [14]: (ux , uy )> = ∇φ

Is (sx , sy ) = Et (φ , ∇φ )

2 2 ! ∂ φ ∂ 2φ ∂ 2φ − ∂ s2x ∂ s2y ∂ sx ∂ sy

(9)

(10)

This second order nonlinear PDE is called Monge-Ampère equation. In general, the mapping function is more complicated due to the nonlinearities of the optical system and the resulting equation is of General Monge-Ampère type. Because of the strong nonlinearities these equations are difficult to solve. Several solution methods have been described in literature: 1. Explicit solution: Express the mapping in terms of the parametrization of the optical surface, discretize the equation and solve it [10], [11]. 2. Implicit solution: Determine the free-form surface such that the difference between actual and prescribed irradiance vanishes. • The actual irradiance is computed using curvatures of the light fields and the optical surface [2]. • Flux instead of irradiance is employed [15]. 3.

Algorithms for extended sources

Besides the SMS method [7] extended source algorithms mostly rely on point source algorithms wrapped in an additional optimizing loop. #202374 - $15.00 USD Received 24 Dec 2013; revised 27 Feb 2014; accepted 27 Feb 2014; published 7 Mar 2014 (C) 2014 OSA 10 March 2014 | Vol. 22, No. S2 | DOI:10.1364/OE.22.00A552 | OPTICS EXPRESS A556

3.1.

Irradiance adaptation

One method to deal with extended sources is to adapt the target irradiance distribution in such a way that the resulting optical surface creates the desired distribution for an extended source. There are different adaption strategies described in literature. The simplest method is to adapt the irradiance distribution point-wise [16]. Fig. 2 schematically shows a method that employs a Monge-Ampère equation to compute the adapted distribution [17].

Fig. 2. Target irradiance distribution adaptation using the solution of a MongeAmpère equation. The desired target distribution a) is fed into a point source method. The distribution b) is computed with the lens of step a) using a point source. Distribution c) shows the case with extended source. Now a mapping is computed via a solution of a Monge-Ampère equation using the method outlined in [14] to transform the distribution c) into distribution b). This mapping is then applied to the original distribution a) which results in distribution d). Distribution d) is now fed into the point source method which results in the distribution e) for the extended source case.

3.2.

A phase space method

A more rigorous approach has to deal with flux conservation formulated in the optical phase space (see Eq. (1)). Multiplying Eq. (3) by cos(θt ) and integrating both sides over target side solid angles yields: Z

Z

Lt (~ut ,~st ) cos(θt )dΩt = Ωt

Ls (~us ,~ss ) cos(θt )dΩt

(11)

Ωt


The integral on the left hand side is the irradiance on the target plane: Z

Et (~ut ) =

Ls (~us ,~ss ) cos(θt )dΩt

(12)

Ωt

The difficulty here consists in efficiently calculating the integral on the right hand side of Eq. (12). A pragmatic approach is to trace rays from the target plane back to the source, relying on the conservation of radiance of Eq. (3) to obtain each ray’s contribution to the resulting irradiance distribution on the target plane. Fig. 3 shows an optical system with an extended source, a spherical lens and a target plane. On the left and right side, two-dimensional slices of the optical phase space are shown. Whereas the phase space on the source side is densely occupied, this does not hold for the target side. In order to avoid tracing significantly more rays than necessary, the target phase space over which is integrated can be limited to the relevant parts. One possibility to constrain it is by coarsely forward ray tracing from the light source and assessing which parts of the target phase space are subsequently relevant for the integral. Another possibility is to limit the integral to those rays that, from the considered target point, intersect the optics. This does not assess which portion of the source phase space is relevant, but we will consider the amount of additional rays that are traced, while not contributing to the irradiance, to be negligible.

Fig. 3. Example setup of a source, a spherical lens and a target plane. Rays are traced from the target plane back to the source plane. Rays that hit the source contribute to the integral Eq. (12). The squares on the left and right show the points in phase space of those rays that hit the source. The phase space on the source is densely occupied whereas on the target side only a small part of phase space is occupied. In order to avoid unnecessary computations the target side phase space has to be restricted to the relevant parts. This is achieved by first tracing forward from the source to the target. The curve designated by ”integration” is the result of integral Eq. (12). The dots represent results obtained by using a Monte-Carlo ray tracer.

The following section will give an example of an approximate evaluation of the integral Eq. (12) and application to the optimization of one free-form surface.


3.3.

Example computation

In the following example, the integral Eq. (12) is evaluated over the portion of the target phase space that originates from the optics, for each considered target point. To achieve this, rays are traced backwards from the respective target points through the optics onto the source plane, where their radiance information is extracted. In the present example we treat a 2-dimensional case. Each of the two surfaces of a lens are described by polynomials of degree 8. The coefficients of the two polynomials are optimized simultaneously. The integral is discretized with respect to the angle θt : Et (~ut ) ≈

1 [Ls,i (~us ,~ss ) cos θt,i + Ls,i+1 (~us ,~ss ) cos θt,i+1 ] ∆Ωt,i,i+1 2∑ i

(13)

∆Ωt,i,i+1 is the solid angle between rays i and i + 1. Evaluating the summation numerically for a set of target points yields an approximation of the irradiance pattern across these target points. For the optimization of the free-form surfaces, a gradient method together with an appropriate merit function is employed. The merit function is defined as the sum over the squares of the differences between prescribed and calculated target irradiances. The power of the prescribed and achieved distributions on the target is normalized to 1 W , as in this example, the light flux passing through the optics is enclosed in the analyzed target area. Fig. 4 shows the set up and the shapes of the two free-form surfaces of the lens. The resulting target distribution is shown in Fig. 5.

50 mm target

25 mm

14 mm lens source 6 mm

Fig. 4. The source has an extend of 6 mm, the lens is 14 mm wide, the target has a width of 50 mm and the distance from source to target is 25 mm. The output angle of the source is restricted to 46.5◦ .

4.

Conclusion

With the advent of high power LED sources, designing general free-form optical elements to realize arbitrary irradiance distributions has attracted much attention within the last decade. Point source algorithms are quite mature today and have been used to design different optical systems for automotive and general illumination, image projection and other applications. #202374 - $15.00 USD Received 24 Dec 2013; revised 27 Feb 2014; accepted 27 Feb 2014; published 7 Mar 2014 (C) 2014 OSA 10 March 2014 | Vol. 22, No. S2 | DOI:10.1364/OE.22.00A552 | OPTICS EXPRESS A559

Fig. 5. Target irradiance distribution generated by the setup Fig. 4 computed using FRED (www.photonengr.com). The prescribed target distribution was uniform within ±25mm.

Most of these algorithms rely on the fact that ray directions are unique, i.e. the sources have zero e´ tendue. Real sources have finite e´ tendue, the point source assumption thus is only an approximation more or less justified. As a rule of thumb the dimensions of the optical elements have to be at least five times larger compared to the source extent in order for the point source assumption to be valid. For smaller ratios the source extent has to be taken into account. Besides the target distribution adaptation method, that is an approximation as well, phase space methods are the only approach to deal with extended sources without employing approximations and thus being limited to large optic to source ratios. The 2D phase space method presented in the present article demonstrates good performance even in case of a small optic to source extend ratio of about two. It can handle non-uniform sources and two surfaces simultaneously, its extension to 3D surfaces is straightforward and currently underway.
