FDTD光学仿真
在这里十分感谢b站,up主Super-郭的幸苦搬运
FDTD(给兄弟们挂个🔗)
FDTD边看边更…
兄弟们,还有个事。就是这个视频,很尴尬,前两节课是中文的,后面实例就是英文的…所以就看了2节内容2333333(没字幕,很难受啊>_<…)
大家还可以看看这个,up主Nestor呐等的字幕翻译…的官网视频
FDTD100(给兄弟们挂个🔗)
一、交互界面初识别:
1.基本物件:
·结构Structures
·仿真设置Simulation
·光源Sources
·监视器Monitors
1.1结构
(1)结构与结构组:
(2)自定义物件:
GDSll可导入cad文件、image可导入照片…
(3)物件库 = 结构库 + 分析库:
(4)高级特性:
可以用于矩阵、介电常数的旋转等…
1.2 仿真区域的设置:
仿真选项 | 设置 |
---|---|
仿真区域结构 | 可以是2D、3D;3D较为常见,2D仿真区域Z轴是常量,(xy横截面的结构),则何种情况下,z方向的结构会被认为是(沿着z方向均匀和无限远的尺寸) |
仿真时间 | 仿真多需要的最长时间段,确保完全衰减;(大的仿真体积/强共振系统)需要长时间衰减,大的仿真时间,否则频域监视器里会有“伪影”(例如:T power 监视器里看见“人造波纹”) |
仿真温度 | 除非添加了“折射率扰动材质”(受温度变化的折射率材料),否则不需要设置模拟温度 |
背景折射率 | 没有明确设置“结构”的任何区域,中的“折射率”(没有设置结构 = 空气) |
设置求解器的:位置和跨度;仿真区域外的结构,不会进行仿真,不会记录数据点
_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3dlaXhpbl80NDk0MTM4NQ==,size_16,color_FFFFFF,t_70)
“自动非均匀”,网格精度为2,这种默认设置;这将建立一个每个波长十个网格点的“非均匀渐变网格”
网格单元,有两种不同的材质:
1.默认的“保形变体0”,大多数情况下都可以用
mesh:用于分辨率更高的"局部"区域仿真…例如某种更细节的特征…
imporant ! ! !这是仿真区域的边界条件:
1.PML:完美匹配层,当光波传播到仿真区域边界(电磁波给予完全吸收,没有反射,相当于传送到无限远…)
2.Metal:金属反射边界,相当于光波到达仿真区域边界,被金属反射回来了…(即,完全反射光场)
3.Periodic:周期性边界条件,适用于周期结构 + 平面波照明情况下(沿轴入射…);若与坐标轴有一定的角度,要有角度的平面使用Bloch边界条件。
4.Bloch:是Periodic的一般形式,Periodic是Bloch的特殊形式,即光波正入射的时候使用…
5.Symmetric/anti-symmetric:对称/反对称边界条件…要结构对称性、光源的偏振也要对称…
时间点:根据“稳定性判据”自动计算,一般不要动
1.3光源:
红色为传播方向 蓝色为偏正方向 灰色地方不能放置监视器
1.4监视器:
折射率监视器: 记录材料特性…折射率…
电影监视器: 记录场随时间的变化过程的动画…
时间监视器: 记录场随时间的变化衰减…记录矢量E(t)/H(t)…确认仿真时间足…一般用于监视点 线/面/体慎用!!!
频率监视器:profile傅里叶变换下…任何单色波长下的稳态结果(稳态传输频谱 和 空间场轮廓);power计算透射率/反射率
模式扩展监视器:用于分析波导或光纤的"特定模式"下的传播功率量…
1.5 设计 结构组/分析组:
anlaies group:
1.传输盒:返回流经监视器的净功率
2.功率吸收分析组:返回空间功率吸收曲线…
3.Q分析组:返回谐振腔的品质因子:
1.6 扫描/优化:
1.7 文本程序编辑:
二、实例(1):二氧化硅镀膜
1.设计:
(1)结构的设计:
基底:设置成sio2,长宽9微米,9微米的厚度
上面镀上50纳米si的层:
将其添加在一个结构组里:
(2)仿真区域的设计:
x,y方向是无限大…设计成周期性perioddic的边界条件
(3)光源的导入:
1.加入一个平面波:
2.波长的设置:
3.计算网格/观察网格:
(4)监视器的导入:
1.放置"折射率"监视器:让其在XZ平面上…电影监视器也是一样…
2."时间"监视器:一般设置的点监视器…勾掉snap to grid…可以自由移动…
3.profile 监视器…勾掉output power…
4.power 监视器…频率点设置成50个…
(5)Run 之前的两项的检测1.材料检测2.内存占用检测
(6)Run,一般的progress不要超过90%-95%…达到100%,说明时间太短了…需要修改…
2.仿真结果:
(1)"折射率"监视器:
(2)"时间"监视器:监视 E(t)/H(t) 的 变化…
观察"选取点"的光谱变化…
下面则是:其是"2"个不同监视器,对Re(振幅)的捕捉…
(3)"profile"监视器:…傅里叶变换下…任何单色波长达到稳态下的结果…
(4)"power"监视器:
R T 的分布…
3.script editor:
f=getdata("R","f");
T=transmission("T");
R=-transmission("R");
A=1-R-T;
plot(c/f*1e6,R,T,A,"wavelength μm","R T");
Lenged('R','T','A')
其中A为其"吸收曲线":
4.扫描/优化:
result:
1.si的厚度进行 sweep,扫描结果暂时没有,用分析组去设置
2.分析组进行设置:
3.sweep里再添加"扫描结果":
4.运行结果: