文章来源：Bäuerle A, Bruneton A, Wester R, et al. Algorithm for irradiance tailoring using multiple freeform optical surfaces[J]. Optics express, 2012, 20(13): 14477-14485.
主要内容：本文基于最优化质量传递（optimal mass transport）多自由曲面光学的设计算法，结合照明设计软件FRED模拟了洗墙灯实例。使用该算法，可直接调整多面光学系统（在本文中，我们用了双边自由曲面）获得了近似于预先确定照度分布，而同时可捕获从光源发出的大部分光线。增加多曲面自由度可用来减小菲涅尔损失，包含加工上的约束以及减小零件尺寸。
Axel Bauerle, 1, 2 Adrien Bruneton, 1∗ Rolf Wester, 2
Jochen Stollenwerk, 1, 2 and Peter Loosen1, 2
1Chair for the Technology of Optical Systems, RWTH Aachen University, 52056 Aachen,
Germany
2Fraunhofer Institute for Laser Technology ILT, Steinbachstr. 15, 52074 Aachen, Germany
摘要
自由曲面透镜和反射镜的设计能够获得非径向对称的照度分布，且同时保持光学系统的紧凑性。对于点状光源来说，比如LED，为了增加发光效率往往需要捕获一个宽角度的光源。这样往往导致的结果是产生强弯曲光学元件，需要两个透镜面作用于总的光的折射，从而最小化菲涅尔损失。在本篇文章中，我们报道了一个基于最优化质量传递（optimal mass transport）多自由曲面光学的设计算法，并结合光学工程仿真软件FRED对一般照明问题给出了应用实例。
1.前言
在照明应用中，透镜和反射镜可以以一种预先确定的方式上分配光线，市场上对这种透镜及反射镜的需求越来越多。为了获得常用的光分布，光学设计中自由度数量必须高于传统的光学元件，这就引入了自由曲面光学的概念，对此有众多设计算法提出[1-5]。折射式光学元件(透镜)在材料-空气界面会遇到菲涅尔反射，控制其精准的路径来避免损失难度较大。因为随着光线角度倾斜菲涅尔反射在增加，因此理想情况是使用几个自由曲面来增加系统的光学效率，在这种情况下，一个单透镜面足以调整照度分布。
到目前为止文献中发布的设计算法，只有Minano、Benıtez[5]提出的SMS3D方法可以直接裁剪多平滑表面，且同时获得预定的照度分布，甚至在一定程度上可把扩展光源考虑进去。然而，据作者所知，SMS3D算法只是在他们团队内部使用。对于点光源的情况，Ries 和Muschaweck[2]得出了一组偏微分方程来描述单个光学面，但一般对于多光学表面目前还没有报道。
在本文中，基于传输理论的公式[6]，对两个自由曲面和一个点光源发射器的照度调整问题，我们提出了灵活的近似解算法。此外，我们利用FRED软件演示了对于一般照明设计任务的可行性。
2.光线映射：有关光学设计的传输理论
在一般的公式中，质量传输理论讲述了最优路径的计算，此计算允许从初始质量分布到目标质量分布的连续传递。就光学而言，光通量扮演了重要的角色，投射光源描述为在2维空间 上光通量密度 (图1)。类似的，目标空间 上的光通量密度为 。为了清楚起见，在3维空间中，假定 平行于2维平面，光通量密度 和 在各自的局部笛卡尔坐标（x,y）下被参数化（图1）。
图1.映射计算图，点光源投射到平面  及目标照度投射到  ，自由曲面位置在 和 之间。
光学系统（自由曲面）的设计任务相当于发现一个微分同胚映射（光线映射），以便于照度分布转换匹配目标分布：
 （1）
(tx,ty)代表了光源光线透过 上的（x,y）在 上的目标点。因为沿着无穷小灯管从光源到探测目标上光通量守恒，照度转换公式可以写为：
 （2）
Du是u的雅克比行列式，表示沿着路径灯管横截面的压缩或放大， 是常用的算子符号。在整个平面 上积分此公式得到了总的能量守恒关系  。
因为映射u并不是独一无二的[7]，在光学设计任务中主要的限制是发现可引导光学表面连续可微的映射。这也就是所谓的表面法向矢量N可积条件[8]
N*curl(N)=0 (3)
计算强制满足方程(3)的光线映射并不是简单的事情，因为由雅克比判决式(Monge-Ampere-type 方程)可知方程(2)一般等效于非线性二阶偏微分方程。处理两个光学表面而不是一个光学表面使这个问题更具挑战性。
3.近似最优化光线映射
表面法线矢量场直接关系到映射信息（通过斯涅耳折射定律）。因此，即使现在还没有被证实，看起来似乎是可信的：如果光线映射的旋度自身减小，表面法线矢量场的旋度可大幅度减小。
使用最新的质量传递理论有助于实现这一目标。在处理预定的目标函数这方面大部分的工作集中在寻找最佳映射。每个元件的位移、质量权重是具有代表性的二次函数。几种不同的解决问题的方法已经被提出[9]。本文中，我们关注于在图像变换的背景下由Haker提出的一阶、无参数方案[7]。
点光源的光通量分布投射到方形的2维空间 ，因此获得了平坦的光通量密度μ。一个改进的立体投影通常用于获得这个通量密度，与此同时控制可捕获光锥角度。同样的目标上的光通量分布投射到一个平行的方形区域 。如图1所示。
Haker的步骤[7]是首先找到μ0和μ1初始映射关系（通常是在笛卡尔坐标轴上两个连续1维数值积分），初始映射结果用  表示。
Haker演示了μ0和μ1所有映射可用连续变量t表示，同样的，u可以看作t的函数，  。通量符合如下演化方程：
 （4）
Du表示映射的雅克比行列式，  代表在二维空间旋转90度，  是f的解，表示为泊松公式  。
这个方程的稳态解（当t--->∞ ）已经被Haker证明了是二次型性能指标最优化映射。本文中，此方程的解使用Hakert数值技术，使用 作为初始点。注意演化方程本身不使用μ1作为相关信息，因为它已经包含在初始映射 里面了。
用这种方法减小映射的旋度会得到一个好的积分近似条件，如下面第五部分所示。
4.由光线映射构建光学表面
使用映射信息，光学表面的计算可促使光源光线偏转到想要的目标位置。总的光的偏转分成几个部分（在如下的例子中描述了两种），对于一个给定的光学面，每个都可实现。在本文中，构建方法是使用标准的最小二乘优化法。表面使用三角形网格画法顶点为i=1,……N，他们的位置为由下列式子给出：
 （5）
r0(i)是光线i的起点，可以是点光源的位置或光线通过另一个光学面的位置。S(i)是光线的方向矢量，λ(i)是标量参数定义表面点i。
S(i)和N(i)给出了光线折射或反射后方向矢量S0(i)。法向矢量N(i)在三角形面顶点位置由顶点相邻面法向平均加权计算得到。
鉴于S0(i)、光线在光学面上的位置r(i)、光线交于目标面的点T(i)都可以计算。由下式(6)最小化作为目标函数
 （6）
Tx(i)和Ty(i)是给定的参数矢量 的光线局部目标面的坐标，tx(i)和ty(i)分别是映射算法中想得到的局部目标坐标。
表面构建是基于目标面上的光线坐标符合积分条件。算法与在目标面上的光线位置协调进行，因此可在多个光学表面偏转光线，引入额外的自由度到设计上而并非算法上，因为算法只考虑单个面。
使用三角网格算法的优势体现在可非常灵活的进一步处理结果表面。例如他们可以细化处理或快速切割，因此，引入了出色的透镜边界控制条件。结合多面设计允许包含加工限制条件（如注塑成型可避免凹陷）。
5.建筑照明：洗墙灯实例应用
图2给出了标准的、高光学效率，建筑照明上具有挑战性几何结构，吸顶灯距离均匀照度接受面80cm。第一步，建立一个简单的轴上演示设计，并对比了传统的非-优化的映射(1维逐次积分)设计。基于此，建立了更具有挑战性的倾斜设计。
图2.（a）自定义坐标轴的洗墙灯应用 (b)简单的轴上方案
5.1 轴上实例
依据上面列出的算法，一个初始的光线映射计算出来，并构建了相应的光学表面结构（假设材料是PMMA，n=1.49），最小化目标函数(公式6)。为完善设计，对如上的初始映射进行了优化，第二组表面随后被构建。图3给出了相应的映射，分别表示为规则网格的变形及z分量的旋度。当系数为200时总的旋度量急剧减小。
图3. 光线映射（表示为规则网格的变形）及相应的Z分量旋度映射，(a) 优化前 (b) 优化后
在两种表面处执行了蒙特卡洛光线追迹，图4给出了每种情况的最终照度分布。500万条光线、目标面上探测器分辨率点为71*71，所有的这些均由FRED[10]光学工程仿真软件来完成。光源是标准的朗伯体点发射器，半角分布为70度圆锥角，方向沿着z轴。
尽管透镜的形状极其相似，图4（d）中的光分布显示了与均匀照度分布相比显示了极小的变形。证实了旋度减少得到了令人满意的光线映射，例如与初始映射相比更好的与积分条件（公式3）的匹配，同时更接近想要的光通量变化。
图.4 FRED软件线框图示及用蒙特卡罗光线追迹后辐照度分布（任意单位）。(a)和(c)没有映射优化;(b)和(d)映射优化后 |
