lumerical FDTD、MODE自学日记

稻草月亮

已于 2023-09-19 17:30:30 修改

阅读量6.8k

点赞数 23

文章标签： python 1024程序员节

于 2022-09-07 09:15:36 首次发布

本文链接：https://blog.csdn.net/weixin_43413070/article/details/126698924

版权

寻学习路径和资料

官网问题搜索：

https://optics.ansys.com/hc/en-us

2022/9/07 3

看b站的lumerical视频

Lumerical 基础入门型教程（中文）6—— MODE Solutions EME操作实例演示_哔哩哔哩_bilibili

看到第6个。第78个讲dievce。很多人安利官网，但是这时候官网上不去两三次我就放弃了。应该找不同时间段尝试。多登几次，就能上去。

积极学习优化教程，学会优化与扫描。官网的教程如下链接：学完官网，还不会就找百度知乎csdn

https://optics.ansys.com/hc/en-us/sections/360005085013-Optimizations-and-Sweeps

2023 9 19

调mode看模式，原始设置是看5个模式，覆盖更多模式可以调监视器：

2023 9 1 5

为了搞明白基横模和高阶模怎么样弄出图来，问了下y老师，现在满脑子都是知识非常快乐。基横模是要看有效折射率的，调有效折射率最大的模式当作基模。论文里对基模的定义是q量子数为1的时候是基横模，也就是一圈的时候。mode可以看模场分布，捣鼓一下，尝试找到不同的光场图。

2023 6 30 5

解决了基本的Q因子怎么测，看到官网有讲解，我看了官网Q因子分析三四遍，以为要用脚本解决，其实用器件库的窗口，直接引入就可以了。问师姐很有帮助，自己看不明白那句话是什么意思，因为没用过那个窗口。给师姐买水果感谢一下！仔细看了一遍讲解，大致看懂了，要一句一句弄明白。

2023 6 21 3

明天端午节，我之前在捣鼓Q因子怎么测，今天得知2种Q因子测量方法.除了从object 小窗里添加Q分析组，另一种是用波峰所在波长除以半高宽，一个师姐和我说是激射光谱的实验一般这样处理。但是我的图是谐振的纵波波谱，可能这么算不准确。是一个突破吧，解决了

2023 6 9 5

光源调整的时候看到了添加光源代码的中文教程，还挺好，我需要再找找详细介绍不同光源的区别的文章。

FDTD快速入门之Lumerical脚本语言Pickup（五）光源（Source）_fdtd光源_ArcSight的博客-CSDN博客

想知道电偶集资和磁偶几次两种光源有什么区别。

2023/3/7 er

耦合器比较好的理想版本命名为plan1133，从监视器的view，T中可以看到耦合效率。产生这种稳定性的原因是通过增加一个谐振腔对波产生了筛选，三棱柱的槽产生了谐振的效果，使光波传输更加稳定。

查报错的链接网址：

Overlapping farfield profile with first gaussian mode

脚本控制学习

官网有脚本语言的查询文档，可以对照着找相应的圆环建立语言，https://optics.ansys.com/hc/en-us/articles/360037228834-Lumerical-scripting-language-By-category

设置材料就set（“metrial”，“GaAs”）；

deleteall;
r = 0.5e-6;
R = 3.5e-6;
addring;
set("name","GaAs buffer");
set("x",0);
set("y",0);
set("inner radius",r);
set("outer radius",R);
set("z",11.75e-6);
set("z span",1.5e-6);
set("theta start",0);
set("theta stop",0);
set("material","GaAs");
#first ceng
#ALGaAs ceng
addring;
set("name","Al0.4Ga0.6As");
set("x",0);
set("y",0);
set("inner radius",r);
set("outer radius",R);
set("z",13.375e-6);
set("z span",1.75e-6);
set("theta start",0);
set("theta stop",0);
set("index",3.25);
#
addring;
set("name","GaAs");
set("x",0);
set("y",0);
set("inner radius",r);
set("outer radius",R);
set("z",14.4775e-6);
set("z span",0.455e-6);
set("theta start",0);
set("theta stop",0);
set("material","GaAs");
#
addring;
set("name","Al0.4Ga0.6As");
set("x",0);
set("y",0);
set("inner radius",r);
set("outer radius",R);
set("z",15.58e-6);
set("z span",1.75e-6);
set("theta start",0);
set("theta stop",0);
set("index",3.25);
#
addring;
set("name","GaAs");
set("x",0);
set("y",0);
set("inner radius",r);
set("outer radius",R);
set("z",16.605e-6);
set("z span",0.3e-6);
set("theta start",0);
set("theta stop",0);
set("material","GaAs");

2023 5 9

418nm谐振的微盘激光器

https://optics.ansys.com/hc/en-us/articles/360042088353-Whispering-gallery-modes-of-a-microdisk

测谐振频率，可以使用偶极子光源和时间监视器来测量谐振。偶极子光源设置在400nm

到450nm波段之间，激励谐振。光源放在边缘上。

2023 5 7

要找到对称边界的设置方法，把环形区域设置成四分之一的状态，减少运行时间，细化网格。

2023、5、6

测量Q的官网教程https://optics.ansys.com/hc/en-us/articles/360041611774-Quality-factor-calculations-for-a-resonant-cavity

2023 4 26 3

single-mode

学习单模激光器的设计，有一篇文章：

https://opg.optica.org/oe/fulltext.cfm?uri=oe-29-9-13193&id=450156

2023 4 25 2

看官网的激光器仿真的过程，模仿它设定监视器和阅读出图的方法。

耳语廊微盘谐振频率：resonant frequencies

https://optics.ansys.com/hc/en-us/articles/360042088353-Whispering-gallery-modes-of-a-microdisk

apodization

为了找到结构的谐振频率，我们可以从时间监视器中获取场随时间变化的傅里叶变换，然而，对于高Q值的设备，使用具有端点化的频域监视器是更好的，因为场在模拟结束时不会完全衰减。更多信息可以在振镜页面找到。在仿真文件中，使用了一个点监测器，也是放在靠近结构边缘的地方。

在找到共振频率的位置后，我们可以通过编辑 "常规 "选项卡，设置频域剖面监测器来记录共振频率下的场剖面。

（To find the resonant frequencies of the structure, we can take the Fourier transform of the field over time from a time monitor, however, for high Q devices, using a frequency domain monitor with end apodization is preferable since the fields do not decay fully by the end of the simulation. More information can be found on the apodization page. In the simulation file, a point monitor is used and is also placed close to the edge of the structure.

After finding the location of the resonant frequencies, we can set the frequency domain profile monitors to record the field profile at the resonant frequencies by editing the General tab.）

下面的图显示了来自频域功率监视器的场谱。使用findpeaks脚本显示，谐振频率为x和y，对应的波长约为418.2纳米和428.6纳米。使用更细的网格可以得到更准确的结果。

（The following plot shows the field spectrum from the frequency domain power monitor. Using the findpeaks script shows that the resonant frequencies are x and y which corresponds to wavelengths of approximately 418.2 nm and 428.6 nm. Using a finer mesh can give more accurate results.）

看了一篇关于振镜的光网教程，没看懂。https://optics.ansys.com/hc/en-us/articles/360034902473-Understanding-time-apodization-in-frequency-domain-monitors

2023 4 24 一

被安排了任务，看看fdtd和device的联合使用方法，准备瞅瞅lumerical怎么看激射。

2023.4.20 4谷雨

不行得看看监视器的文档，不了解监视器是干嘛用的，找到才能复现文章。光源也得看看，有改动的地方才能灵活。

2023.4.18 2

阅读微盘激光器的耳语模式例子：Whispering gallery modes of a microdisk

https://optics.ansys.com/hc/en-us/articles/360042088353-Whispering-gallery-modes-of-a-microdisk

教程中通过使用偶极子源和时间监视器运行初始模拟来定位结构的谐振频率。将偶极子源注入的波长范围设定为400和450纳米之间，以激发这个范围内的共振频率。还将源放置在结构的边缘附近，因为这是预计模态场会很强的地方。频域功率监视器的场谱。

2023、4.17 1

找到一个Ring Vernier laser，浏览一下。感觉不能用

找到一个像微盘激光器的例子：Microdisk laser

https://optics.ansys.com/hc/en-us/articles/360041721234

再找一个ring的结构教学。环形可以设置绘制角度，绘制缺角圆环：https://optics.ansys.com/hc/en-us/articles/360034382194-Ring-Simulation-Object

关于激光器的案例目录，里面只有一个长得有点儿像微环激光器，链接没打开：

https://optics.ansys.com/hc/en-us/articles/360043683333-Lasers-and-gain-list-of-examples

2023/4/13 4

最近在复现一个微环激光器，需要恶补一下微环激光器的知识。看看论文吧，微环激光器的论文，微盘的也行，之后可能要用它的优化微环的宽度。

2023/4/6 4

找到一些mode模型，结果模型里没有带扫描，所i只能自己查查相关交叠积分和耦合效率以及对准容差的方法。

2023/3/31 5

半懂不懂之间，做了一个新的扫描。扫出来看看什么情况。

果然，报错了。

https://optics.ansys.com/hc/en-us/articles/360034405254-overlap-Script-command

再看看已有扫描的代码，看看官网扫描部分的案例。然后是交叠积分的案例。

2023/3/29 3

捣鼓一下已有的扫描，尽量看懂

捣鼓自己需要做的扫描

2023/3/28 2

学会选择光源模式.

从光源的设置里，用用户选择模式，可以点select mode，

从里面选好，点确认select就可以，能看得到带结构黑色轮廓的光场图片。

然后直接从图片上读数，确定激光光斑直径。

从光源的监视器里找E，可以得到出射光斑直径大小。

2023/3/27 1

交叠积分overlap

https://optics.ansys.com/hc/en-us/articles/360034405254-overlap-Script-command

2023/3/20 1

一个锥形耦合器的FDTD官网教程没看明白的部分：

Change taper length

To recalculate results for a longer taper, go to the CELL GROUP DEFINITION section of the EME Analysis window and change the GROUP SPAN of the second cell group region. Change this from 10um to 100um, which will correspond to the same device geometry but now with a 100um long taper region.
Press the EME PROPAGATE button to re-calculate the results, then visualize the new results for the longer taper. Note that this calculation is almost instant since it does not require any additional mode calculations.

Scan taper length

To scan the taper length over a range of values, the propagation sweep widget in the EME analysis window can be used.

In the EME Analysis window, select the PROPAGATION SWEEP checkbox, and set the settings of the propagation sweep according to the following table.

setting	value
parameter	group span 2
start	10
stop	200
number of points	191

Press the EME SWEEP button to run the sweep over taper lengths.

Plot sweep results

After the sweep is complete, press the VISUALIZE EME SWEEP button. The S-parameters will be plotted in a new visualizer.
The transmission corresponds to |S12|^2. To plot this, remove the other results from the list of attributes in the Visualizer, and set the SCALAR OPERATION of the result to Abs^2.

Discussion and results

The EME solver in MODE is a fully vectorial bi-directional Maxwell's equations solver. The solver relies on modal decomposition of electromagnetic fields into a basis set of eigenmodes, which are computed by dividing the geometry into multiple cells and solving for the modes at the interface between adjacent cells. This method accounts for multiple-reflection events, and only one simulation is needed for all input/output modes and polarization so it is ideal for simulating tapers and performing length scanning.

In the EME solver setup, we define the cross sections where the modes are solved by defining cell groups. For uniform regions where the cross section of the structure does not change in the propagation direction (ex. cell group 1 and 3, or the input/output waveguide regions), only one cell is necessary in the cell group since using additional cells will not affect the results.

For regions such as a taper where the cross section of the device varies, you can specify a number of cells within the cell group where the modes of the structure will be calculated, and in these regions we want to set the subcell method to CVCS which reduce the staircasing effect due to the discrete changes in cross section of the structure between each adjacent cell. The number of basis modes to use for the calculation can also be set in the EME solver object. It is recommended to start with a small number of modes for the initial calculation. Once everything is working as expected, one can increase the number of modes until the result converges.

The cell boundaries of the structure can be seen in the CAD view below.

The modes at the center each cell are calculated on a finite mesh transverse mesh. Mesh override regions can be added to force a finer mesh where necessary. For this spot size converter, we add a mesh override region over the tapered silicon waveguide to better resolve the geometry. The view mesh button displays the transverse mesh in the CAD view as shown below.

We can select the mode (or a set of modes by multi-selecting) of interest by editing the ports and choosing the desired modes. The user s-matrix result that is calculated by the EME solver will return the results for the selected mode(s) only. For this device we are interested in the fraction of power transmitted from the fundamental mode of the silicon waveguide at port 1 to the fundamental polymer mode at port 2 which is given by |S21|^2 with port 1 at the input side, and port 2 at the output. However, since the device behaves symmetrically, we can get the same result by looking at |S12|^2. For more information about the S-matrix index mapping see EME solver analysis.

Analysis and Results

Pressing the run button will calculate the modes at each cell. You can visualize the calculated modes by expanding the EME solver and cell group in the Objects Tree, then right-clicking the individual cell and selecting the result to visualize.

To see the final field profile of the device as well as the S-matrix results, press the EME PROPAGATE button in the EME analysis window. Once the propagation is complete, profile monitor results and S-matrix results will be available, and can be visualized by right-clicking on the objects in the Objects Tree. The results for different propagation lengths can also be changed without having to re-calculate any modes. The field profile for a tapered region of length 10 um and 100 um are shown below.

10 um taper (xz plot)

100 um taper (xz plot)

Scattering parameters relate the transmission and reflection coefficients for each port and input/output modes of the device. This is automatically calculated by the EME solver, and returned as the result of an EME solver region. The internal s-matrix includes all of the s-parameters for all the modes of all the ports, whereas the user s-matrix will contain only the s-parameters for the modes selected in the ports. Since we have 2 ports, and we are only interested in the fundamental mode at each port, the user s-matrix will be a 2 by 2 matrix, with elements S11, S12, S21 and S22.

Length scanning

The propagation sweep widget allows you to scan the length of any cell group and calculate s-matrix results automatically. The S-matrix index mapping table allows you to quickly identify which s-matrix components correspond to which port and mode.

Below, the transmission through the taper is plotted over taper lengths from 10 um to 200 um.

The length scanning can also be done by running the script spot_size_converter.lsf.

EME vs 3D FDTD

We also compare the EME results with 3D FDTD. The results between two solvers agree reasonably well, however they are done with completely different time scale. The EME simulation takes 3 minutes to simulate 101 different taper lengths (blue squares), whereas 3D FDTD takes 6 hours to simulate 11 different taper lengths (green squares).

None vs CVCS subcell method

To see the effect of staircasing, change the subcell method for group span 2 from "CVCS" to "none" and re-run the eme sweep.

One can see that when the CVCS subcell method is not used for the tapered portion of the structure in cell group 2, the staircasing effect will result in a transmission curve that is much rougher than before.

References

[1] T. Tsuchizawa et al, “Microphotonics devices based on silicon microfabrication technology”, IEEE J. Select. Topics Quantum Electron., 11, 2005, 232-240

#########################################################
# file: optimization_AR_coating_example_script.lsf
#
# Description: This script file
#   is used to generate and run the thickness
#   optimization as well as plot the results
#
# Copyright 2015, Lumerical Solutions, Inc.
#########################################################

closeall;
clear;

# add a new optimization and set basic properties
addsweep(1);
setsweep("optimization", "name", "thickness_optimization_script");
setsweep("thickness_optimization_script", "Type", "Minimize");
setsweep("thickness_optimization_script", "algorithm", "Particle Swarm");
setsweep("thickness_optimization_script", "maximum generations", 20);
setsweep("thickness_optimization_script", "generation size", 10);
setsweep("thickness_optimization_script", "tolerance", 0);

# define the parameter thickness
para = struct;
para.Parameter = "::model::AR structure::thickness";
para.Type = "Length";
para.Min = 0.05e-6;
para.Max = 0.15e-6;
para.Units = "microns";

# add the parameter thickness to the optimization
addsweepparameter("thickness_optimization_script", para);

# define figure of merit
result_1 = struct;
result_1.Name = "R";
result_1.Result = "::model::R::T";
result_1.Optimize = true;

result_2 = struct;
result_2.Name = "T";
result_2.Result = "::model::T::T";
result_2.Optimize = false;


# add the figure of merits R & T to the optimization
addsweepresult("thickness_optimization_script", result_1);
addsweepresult("thickness_optimization_script", result_2);

# run optimization
runsweep("thickness_optimization_script");

# get & view the results - parameter value
R = getsweepresult("thickness_optimization_script", "parameter trend");
value = R.getattribute("parameter value");
gen = R.getparameter("generation");

plot(gen, value / 1e-9, "generation", "parameter value (nm)", "thickness");
legend("");

阅读Creating parameter sweeps using a script，今天学习的两个都是纯脚本命令。阅读了一遍加强熟悉。

第四课全脚本生成优化

#########################################################
# file: sweep_AR_coating_example_script.lsf
#
# Description: This script file
#   is used to generate and run the thickness
#   sweep as well as plot the results
#
# Copyright 2015, Lumerical Solutions, Inc.
#########################################################

closeall;
clear;

# add a new sweep and set basic properties
addsweep;
setsweep("sweep", "name", "thickness_sweep_script");
setsweep("thickness_sweep_script", "type", "Ranges");

# define the parameter thickness
para = struct;
para.Name = "thickness";
para.Parameter = "::model::AR structure::thickness";
para.Type = "Length";
para.Start = 0.05e-6;
para.Stop = 0.15e-6;
para.Units = "microns";

# add the parameter thickness to the sweep
addsweepparameter("thickness_sweep_script", para);

# define results 
result_1 = struct;
result_1.Name = "R";
result_1.Result = "::model::R::T";

result_2 = struct;
result_2.Name = "T";
result_2.Result = "::model::T::T";

# add the results R & T to the sweep
addsweepresult("thickness_sweep_script", result_1);
addsweepresult("thickness_sweep_script", result_2);
setsweep("thickness_sweep_script", "number of points", 10);  

# run the sweep
runsweep("thickness_sweep_script");

# view the results
R = getsweepresult("thickness_sweep_script", "R");
T = getsweepresult("thickness_sweep_script", "T");

plot(R.thickness*1e9, R.T, "AR thickness (nm)","R");
plot(T.thickness*1e9, T.T, "AR thickness (nm)","T");

#visualize(R);
#visualize(T);

2022/10/24 1 九周

早9：37官网打不开，看看matlab编程关于优化的部分代码。FDTD是麦克斯韦方程旋度方程的离散空间网格取值求解，是广泛使用的数值方法。

reflection = getsweepresult("wavelength", "R"); 
R =     -reflection.T;
lambda = reflection.lambda_sweep*1e9;

plot(reflection.source_angle, pinch(R,2,1), pinch(R,2,2), pinch(R,2,3),
     "angle of incidence (degrees)","Reflection","Reflection vs angle");
legend('lambda = ' + num2str(lambda(1)),
       'lambda = ' + num2str(lambda(2)),
       'lambda = ' + num2str(lambda(3)));

下午3：51官网可以打开。

2022/10/21 5

学习扫描与优化，第二课官网案例里有一个扫描单个参数的例子，硅上材料层厚度优化：

硅上抗反射层厚度ARhttps://optics.ansys.com/hc/en-us/articles/360034922873，有获取结果的脚本。（获取结果之前，先分析模式，先运行相应设置好的扫描块，再运行完之后再获取结果）

第二课是波长与入射角两个参数，对反射情况的变化显示。中午11点左右网不好，登不上官网，改日再战。下午没有赶上听nlp，又把官网第二课案例学了一遍，感觉有所收获。

2022/10/20 4

上午网络不给力，官网链接打不开。看fdtd原理的书，巩固fdtd的麦克斯韦方程，麦克斯韦方程离散化，用一个方程设定离散步长。

2022/10/19 3

应对如何验证算法正确，是否使用的办法？一，原理是这样的，可以通过原理推导。二，用我的办法算过例子，结果一致。三，用其他计算方法计算，结果相似。四，科学理论提出的结果都有怀疑余地，不能因为有怀疑余地停滞不前。五，做出来实验。

找到一个很炫酷的博客FDTD学习笔记 - 随笔分类 - lihao_Q - 博客园

有一篇号称一篇文章入门fdtdhttp://www.360doc.com/content/22/0520/20/32196507_1032382435.shtml

2022/10/18 2

早上网不太好，调了一下，反复重登，15分钟以后就好了。可能是我妈保佑我。

官网上第一课:

https://optics.ansys.com/hc/en-us/articles/360043102714-Corner-analysis-utility

扫描流程还有相应的脚本，脚本流程很清晰，但是一些命令看不懂：

clear;
switchtolayout;
project_name = "corner_analysis.icp";
installdesignkit("lumfoundry_template.cml", "./", true);
load(project_name);

analysis_name = "corner_analysis";
try{deletesweep(analysis_name);}

# define corner analysis
addsweep(4);
setsweep("Corner sweep", "name", analysis_name);
setsweep(analysis_name, "resave files after analysis", 0);

# define corner
corners_L = struct;
corners_L.name = "MMI";
corners_L.Library = "./lumfoundry_template/lumfoundry_template.lib.x";
corners_L.Corners = "fast; slow; nominal";
addsweepparameter(analysis_name, corners_L);

# define result
results_L = struct;
results_L.name = "transmission";
results_L.result = "::Root Element::ONA_1::input 1/mode 1/transmission";
addsweepresult(analysis_name,results_L);

# run corner analysis
runsweep(analysis_name);

# visualize results
result = getsweepresult(analysis_name, results_L.name );
visualize(result);

今天官网很给力（10：00），服务器登不上去但是自己的电脑可以登录。第二课创建嵌套参数扫描，嵌套扫描即扫描多个参数，单个参数扫描单元是第二课多参数扫描的内置课程，是一个官网调厚度的例子：https://optics.ansys.com/hc/en-us/articles/360034922873

下午网不太行（4：00），看完第二课收工。

2022/10/17 1 10

学优化与扫描，（9：30开官网连不上。

百度相关优化，看到一篇文档，文字叙述新建一个结构组，在其中使用脚本语言优化参数，paraxx是优化对象的命令。百度优化透镜参数，看到一些命令行，看了一遍留个映像。

看到csdn上一个超透镜的优化案例，是官网的案例翻译，我没太看懂，看了一遍，感觉比百度那个文档还不好理解，但是有案例下载链接，码下来。Lumerical官方案例、FDTD时域有限差分法仿真学习（十四）——超透镜（Metalens）_子非鱼icon的博客-CSDN博客

这篇有配套的个人讲解理解，从代码角度阅读官网文章。image可以绘制相位随高度和半径变化的二维图像。【FDTD Lumerical sweep扫描参数数据提取及可视化【以S参数为例】】_Pale blue dot的博客-CSDN博客_fdtd扫描脚本

下午4：49 官网连上了，可以看超透镜的链接https://optics.ansys.com/hc/en-us/articles/360042097313

也可以看扫描与优化的第一个链接课

2022/10/13 4

I have learned to optimization（优化）

找到了可以切入的点，优化模斑转换器的长度宽度来反向设计一个两层的模斑转换器。如果学会粒子群算法，就加上用粒子群算法做优化的部分，设计一个双层的模斑转换器。

积极学习优化教程，学会优化与扫描。官网的教程如下链接：学完官网，还不会就找百度知乎csdn

https://optics.ansys.com/hc/en-us/sections/360005085013-Optimizations-and-Sweeps

2022/10/11 2

搜索官网有原理教程，但是视频需要翻墙才能看，想知道有什么办法可以看到这些视频。文字暂时能看到，先看文字教程吧。

In this course, you will learn about the underlying solver physics and numerics of Ansys Lumerical FDTD, the types of problems it can solve, and how it can best be used for parallel computation. By the end of this section, you will be able to:

Briefly explain what FDTD is and when it should be used
Identify some applications where FDTD can be used
Understand the consequences of using a finite-sized mesh
Know when to use 2D vs 3D simulations
Understand how frequency domain results are obtained from time-domain simulations

一片参考文献

. Dennis M. Sullivan, Electromagnetic simulation using the FDTD method. New York: IEEE Press Series, (2000).

Understanding Mesh Refinement and Conformal Mesh in FDTD – Ansys Optics

zai知乎看到一篇自学记录，可以参考下自学Lumerical FDTFD——FDTD Algorithm - 知乎 (zhihu.com)

2022/10/10 1

我有一个更加清晰的目标，通过这个星期的学习实现复现论文仿真的任务。

2022/10/07 5

今天是学习lumerical fdtd solution 的一个月，这个月从找到lumerical 软件的b站视频开始看，到现在对官网的使用刚入门，b站视频看了以后，对模型的初步搭建有了概念，接下来跟从官网教程，熟练一些步骤。官网不一定能顺利登陆。对fdtd原理的理解，麦克斯韦旋度方程，看了一遍讲解，感觉当下是理解了，过两天会不会忘就不一定。对官网的挖宝。或者对杨老师的模型的理解挖掘。

2022/10/5 3

find a 用户手册

https://optics.ansys.com/hc/en-us/articles/360033154434

2022/9/30 5

找到官网答疑的论坛，可以从上面提问或搜索。

2022/9/22 4

我学fdtd达到了第一个20小时，我知道了模型怎么搭，照着一个流程我会搭起光栅，会算E T H P，照着指引会搭建模型了，大致了解搭建模型的流程。不会用脚本搭建模型和计算。下一个阶段的目标，学会用算法优化模型，预计20h彻底熟练。

fdtd快速入门博客教程看完，脚本可以完成从物理模型建立到fdtd区域建立，光源监视设置的全部流程。数据获取可以用get得到。感觉可以看看那本中文书和官网教程了，看几个实例，自己算一下。也可以看网盘里的两本教程，浏览，选择部分精读。

找到了官网计算损失的帖子：计算吸收的光功率 - 简单方法 – Ansys光学

论坛：Ansys Photonics中文论坛.
模型：Simulation objects.
脚本命令合集：Scripting Language.

用代码分析波段损耗：FDTD快速入门之Lumerical脚本语言Pickup（一）序言_ArcSight的博客-CSDN博客_fdtd脚本语言

在FDTD中创建锥体在FDTD中创建锥体_Temmie1024的博客-CSDN博客_fdtd圆台结构

2022/9/23/5

需要加强的，边界条件的使用情形效果和优缺点，网格精度的选择和dx的设置情形和优缺点。快速入门教程里仿真区域的搭建部分提及mesh的设置，但不详细，找到资料就在看一下。

设置光源的偏振有什么意义？需要查一下

快速入门博客教程有第七部分(15条消息) FDTD快速入门之Lumerical脚本语言Pickup（七）数据获取与分析（Analysis）_ArcSight的博客-CSDN博客_fdtd监视器结果怎么分析

用脚本搭建角锥、椭圆：

(15条消息) FDTD快速入门之Lumerical脚本语言Pickup（三）结构(Structure)搭建_ArcSight的博客-CSDN博客_fdtd脚本语言

阅读一篇博客[个人笔记]FDTD代码——边用边更_Temmie1024的博客-CSDN博客_fdtd代码

Lumerical的脚本语言允许你自动完成任务和分析，如操作仿真对象、启动仿真和分析结果。脚本命令可以直接输入到脚本提示符中，也可以从保存的脚本文件（.lsf）中运行。你也可以参加介绍性的脚本AIC学习课程以了解更多信息。该学习课程的视频也可以直接从脚本--视频列表中获得。脚本语言清单链接：（这个网页今天下午不能上，上午可）https://optics.ansys.com/hc/en-us/articles/360037228834-Lumerical-scripting-language-By-category

FDTD add语句部分详细内容_Temmie1024的博客-CSDN博客

可以将fdtd与matlab link：matlab link fdtd 边用边更_Temmie1024的博客-CSDN博客

2022/9/20 2

复现了b站视频第二个例子，流程学习。是我构建的第7个模型。

2022/9/16 5

找到一个例子，复现它的计算和搭建

理解搭建和计算重要自变量的选择原因

熟练化

照着祝一个微盘激光器搭建，在弄完物理模型之后，对分析器和区域的建立顺序产生了困惑。

解决方案：再看一下b站视频，在心中彻底熟悉器件结构，看一下参考手册的例子和说明

我用的是fdtd和mode

第一个视频讲三个部分。第二个部分讲界面交互。fdtd 有限差分时域法：为什么使用fdtd？波长等级的时候请使用fdtd。优点是可以查看宽光谱和主要模式分布。交互界面。文件视图设置，左侧橙色格子是仿真格子。设置边界条件和监视器。一个非常重要的点是RUN进度在不超过95%的时候可以准确预测，方法是预先设置autoshutoff min。参数扫描和优化。扫描优化和设置。

第二个视频是操作实例。建立物理结构，添加仿真区域和边界，添加光源和监视器。检查材料拟之后就可以运算分析。（视频不是高清的，凑乎看一下。）讲一个硅衬底上用etch刻蚀一个槽。加周期边界，加光源，加所有监视器。然后分配资源算一下。扫描设置参数就可以优化。

第三个mode solution。40：00以后有个例子光纤。最后还有一个例子是脊波导，没详细指南。搭建物理模型，选择仿真区域，设置边界条件偏好金属（但金属不能算弯曲损耗）。仿真光纤，光子晶体，还有脊形波导。设置时注意细节。

第四个2.5Dfdtd，设置一下绿色点和四个蓝色点的位置。有例子，设置细节。