Lumerical——模拟（Simulation）详解

模拟实体对象用来定义模拟参数如边界条件和网格尺寸。

模拟区域（Simulation region）

        模拟区域定义模拟计算的大多数参数，包括模拟区域和网格尺寸。

网格覆盖区域（Mesh override region）

        网格覆盖区域用来覆盖模拟区域某些部分的缺省网格。通常网格化参数是在模拟区域设定的。但是，如果模拟区域的某些部分需要一些特定的网格化条件，就可以指定网格覆盖区域。例如，某些包含材料的模拟计算对网格化效果非常敏感，金属表面可能需要网格覆盖区。

注意：每个模拟计算仅有一个模拟区域，但可有多个网格覆盖区域。

           任何包含后者部分包含在在模拟区域的实体对象均会包括在模拟计算中，而完全在模拟区域之外的实体对象不参与模拟计算。

一、基本属性选项卡（General tab）

        基本属性选项卡中选项取决于被编辑项是模拟区域还是网格细化区域。

1.1 模拟区域（Simulation region）

① 维数（DIMENSION）: 模拟区域的维（2D 或者 3D）。

② 模拟时间（SIMULATION TIME）: 模拟计算需要的最长时间。如果在设定的最大模拟时间前已达到自动结束的标准，实际模拟时间可能会短些。

③ 背景反射率（BACKGROUND INDEX）: 模拟区域周围、背景媒介的反射率。

1.2 网格覆盖区域（Mesh override region）

        对于网格覆盖区域，基本属性选项卡部分包含的选项有：覆盖的不同方向，以及输入所要求的最大网格步长。由于网格步长通常由模拟计算材料的反射率所决定，也就可以设定等效反射率

（EQUIVALENT INDEX） 用来确定网格步长。反射率较大时网格会更细。

二、几何形状选项卡（Geometry tab）

        几何形状选项卡包含的选项用来设定模型的大小和位置，以及网格细化区域。

三、网格设定选项卡（Mesh settings tab）

        可以设定选项卡中的属性可用来控制网格几何形状，以及相关参数。

3.1 网格类型（Mesh type）

        有三种网格生成算法。

（1）自动非均匀（缺省）网格(auto non-uniform)

        根据网格精度滑动条自动生成非均匀网格。强烈推荐初步模拟时网格精度为 1~2，返样模拟计算会快些。如果必要，收敛性测试时可使用高精度的网格。

        “网格精度（MESH ACCURACY）”参数是一个 1~8 的整数，其中 1 是低精度，8 是高精度（较小的网格）。网格化算法中考虑了许多因素，包括光源波长、材料性质、以及结构的几何形状。单位波长范围内网格点数（ppw）是网格化算法中的一个重要的考虑因素。精度 1 对应于 6 ppw，精度2 →10 ppw，精度 3 →14 ppw，依次类推，在滑动条上每点的增量是 4 ppw。重要的一点应当记住，反射率与波长是反比关系。在高反射率的材料中，有效波长较小，也就意味着对于高反射率材料网格化算法将选用比较小的网格。

（2）定制非均匀（Custom non-uniform）网格

        允许用户定制如何创建非均匀网格。如果采用波长设定网格格子，通常采用缺省的设定值 10 就可以了，在进行不严格的模拟计算也可以减少到 6-8。

        等级因子（grading factor）规定了修改网格的最大速度。例如，如果 dx(i+1) = α*dx(i)，返

是，1/（GRADING FACTOR） ≤ α ≤ GRADING FACTOR。等级因子应当在 1~2 之间，缺省

为 sqrt(2)。

（3）均匀（uniform）网格

        均匀网格应用在全体模拟空间体，不考虑任何材料的特性。如果网格覆盖区域同本选项一起使用，则在覆盖区域内会对网格大小加以限制，覆盖区域外部不受影响（网格是均匀的）

① 时间步长选项（Time Step Options）

        时间步长稳定因子（DT STABILITY FACTOR）: 给出模拟运算期间时间步长的大小，定义为 库朗（Courant）数字稳定极限。数值较大时模拟快，数值较小时模拟时间长。根据库朗稳定条件， 如果要保证 FDTD 算法稳定，该设定值应小于 1。

        时间步长（DT）: FDTD/传播模拟算法的时间步长。由空间格子的值所决定，以保证数值的 稳定性，用户不能直接设定该值。

② 网格最小网格步长（Min Mesh Step）

        网格最小步长（MIN MESH STEP）设定所有计算区域的绝对最小网格尺寸。会覆盖对网格的 其它所有设定，包括网格覆盖区域。

四、边界条件选项卡（Boundary conditions tab）

4.1 完美匹配层（PML）边界条件

        许多模拟算法采用吸收边界条件，在这种条件下传播出模拟区域的的光辐射同内部场无作用。 完美匹配层（PML）吸收入射的电磁波。PML 对垂直入射光可有效吸收，对临界入射光会产生较大的反射。

        在特定边界，层数取决于网格尺寸大小（注意：层数越多占用的内存会越大）。更进一步，当所围结构完全扩展到边界区域时PML边界的效果最好。不论是否结构绘制在PML区域内部还是外部，这是其缺省行为。

4.2 金属（Metal）边界条件

        金属边界条件用在完全反射的特定边界，在模拟空间的边界无能量逃逸。

4.3  周期（Periodic）边界条件

        当结构体和电磁场均呈周期时，采用周期边界条件。周期结构条件可用来在一个或者多个方向模拟结构体，结构体的在一个方向而不必在其它方向呈周期性。

4.4 布洛赫（Bloch （FDTD Solutions/Propagator））边界条件

        当结构体和电磁场呈周期性、但每个周期间存在相位移时可采用布洛赫边界条件。布

洛赫边界条件用在下述的两个 FDTD Solutions 和传播模拟算法中：

① 在周期结构体的某个角度已激励一个平面波--在这种情况下，对于一个给定的模拟，可以得

到一个频率点的精确的反射和传播数据。

② 计算周期实体对象的带结构--在这种情况下。偶极子光源的带宽脉冲入射到周期结构体。

4.5 对称（Symmetric ）边界条件

        当用户处理某些具有一个或多个对称面的问题时，可采用对称边界条件。结构体和光源都必项是对称的。对称边界是电场的镜像，是磁场的反镜像。下图可视化的解释了对称性边界条件。

必项认真的考虑所需要的是对称性边界条件还是反对称性边界条件，给出要解决问题的向量对称图。 为了保证获得有意义的结果，光源和边界的对称性必项是同样的。

4.6 反对称（Asymmetric ）边界条件

        当用户处理某些具有一个或多个对称面的问题时，可采用反对称边界条件。对称边界是电场的反镜像，是磁场的镜像。

4.7 边界条件选项卡选项（Boundary condition tab options）

① XMIN, XMAX, YMIN, YMAX, ZMIN, ZMAX BOUNDARIES：给出沿模拟区域周界的边界条件。对称和反对称条件应用于最小（x min, y min 和 z min）边界条件。

② ALLOW SYMMETRY ON ALL BOUNDARIES: 缺省情况下，对称和反对称条件仅应用于最小（x min, y min 和 z min）边界条件。 逅择本选项后，运行将对称和反对称条件仅应用于最大边界条件，以便模拟计算具有对称性的周期结构体。

③ SET BASED ON SOURCE ANGLE: 布洛赫 边界条件通常用来将平面波以一定角度入射到周期结构体中。根据当前模拟计算的光源设定 kx、 ky 和 kz 值。注意：如果定义了多个光源，所有光源布洛赫 边界条件设置并行相同。缺省情况下，本项是选择的。如果不选择，应直接设定 kx、 ky 和 kz 值。

④ BLOCH UNITS: 有两种单位制用来规定 kx、ky 和 kz 值：

        bandstructure: 在这种单位制下，在 kx、ky、kz 的单位为（2π/），其中 是 FDTD 模拟区域的 x，y 或 z 跨距。对于带结构计算，采用返些单位方程方便。  

      SI: 在国际单位中， 的单位定义为 。对于按一定角度入射的平面波，采用这种单位 非常方便。

⑤ KX, KY, KZ: 在 x, y 和 z 方向设定的、单位如上所述的布洛赫边界波矢。

五、高级选项选项卡（Advanced options tab）

(1)SET SIMULATION BANDWIDTH: 缺省情况下，模拟的带宽继承自光源的带宽。选择本选项后可以直接设置模拟带宽。模拟带宽对模拟计算的影响是多方面的，包括：网格生成、材料适应、以 及模拟器的频率范围等。

(2)FORCE SYMMETRIC X, Y, Z MESH: 这些选项会使网格在 x, y 或 z轴方向对称。选择这些选项后，网格算法仅考虑模拟区域的上半部分。负半部分是正半部分的拷贝。网格算法不会考虑所有物理结构体和网格覆盖区的负半部分。本选项强迫网格位于模拟区域的中心，强制了一个对称性网格后，可保证从对称模拟计算转换到非对称模拟计算时网格不变。

(3)OVERRIDE SIMULATION BANDWIDTH FOR MESH GENERATION: 选择本选项后，允许可根据定制的波长或频率范围生成网格，而不是根据带宽（缺省情况下继承自光源的带宽）。

(4)SNAP PEC TO YEE CELL BOUNDARY: 本选项强制任何理想电寻体（PEC）结构体的界面同Yee 胞元的边界对准。这样可以保证在理想电导体的界面电场分量同界面相切。如果在 PEC 界面将电场的法量设定为 0 时，这样设定会避免在某些情况出现异常。要达到这种目的，当生成模拟网格时， 规定的 PEC 界面可以便宜 dx/2 ，dx 是 Yee 胞元的大小。

(5)ALWAYS USE COMPLEX FIELDS: 选择本选项后，在模拟计算时强制算法中使用复场。这会导致的仿真时间变慢和增加对内存的需求，叧有必要时才选择。缺省情况下，叧有在布洛赫 边界条件才采用复场。

(6)USE EARLY SHUTOFF:选择本选项后，当大多数能量离开模拟空间时模拟计算自动结束。

(7)AUTO SHUTOFF MIN: 当模拟空间的总能量下降到入射最大能量的设定分数时，模拟计算将

结束。

(8)USE DIVERGENCE CHECKING: 选择本选项后，当模拟空间的总能量出现多次大于入射次数

等于设定值数时，模拟计算自动结束。

(9)AUTO SHUTOFF MAX: 当模拟空间的总能量达到最大能量的次数等于设定值数时，模拟计算

将结束。模拟数据会自动保存。

(10)DOWN SAMPLE TIME: 检查自动关闭条件（auto shutoff conditions）在 dT 时间步长中每个下采样（Down Sample）的时间次数。

(11)PML KAPPA: 在 PML 边界使用的归一化电导率和磁导率。

(12)PML SIGMA: 在 PML 边界使用的最大电导率和磁导率 。

(13)MINIMUM PML LAYERS: PML 边界条件下单元格的最小数目。增加 PML 局将减少来自边界

的背景反射，但会增加模拟时间并增大对内存的需求。缺省设定为 12。

(14)MAXIMUM PML LAYERS: PML 边界条件下单元格的最大数目。增加 PML 局将减少来自边

界的背景反射，但会增加模拟时间并增大对内存的需求。缺省设定为 64。

(15)PML POLYNOMIAL: 多项式的指数，当当辐射垂直入射进入 PML 时，该值确定电导率和磁导率是以多快的速度增加的。缺省设定为 3，表明随着进入到 PML 深度的增加，电导率和磁导率按3 次方增大。

(16)TYPE OF PML: 使用的 PML 类型为标准的（STANDARD）或 稳定的（STABILIZED）。在某些情况下，STABILIZED PML 数值计算比较稳定，但性能略有降低。

(17)EXTEND STRUCTURE THROUGH PML: 在大多数情况下，我们希望结构体扩展过 PML， 这样离开模拟区域的场就会被完全吸收，否则可能会有一些不必要的反射。本选项会自动地将结构体延伸过 PML，尽管在 CAD 中结构体不是这样绘制的。

(18)MAX SOURCE TIME SIGNAL LENGTH: 是光源监视器用来保存光源“时间”和“时间信号” 属性数据的最大长度。如果使用了大量的光源，以及模拟时间大约为 100ps（非常罕见），高级用户可能想要减少该数值以节省内存。然而，由于“时间”和“时间信号”对计算光源功率（sourcepower）、光源归一化（sourcenorm）以及传播函数的归一化非常重要，必项注意光源的归一化。特别是在大多数情况下，应当使用 nonorm 选项。

(19)SET PROCESS GRID: 将模拟空间划分为可由不同的进程用来处理子区域后，经人工设定以实现较高效率计算。区域为方形时计算效率较好，例如：对于有 20 个区域计算来说，采用 5×2×2 的计算效率比采用 20×1×1 高。

(20)NX, NY, NZ: 并行计算程序布置。

(21)SET DEFAULTS: 将设置恢复到原始缺省状态。