Lumerical官方案例、FDTD时域有限差分法仿真学习（十一）——米氏散射 （Mie scattering 2D）

最新推荐文章于 2024-04-25 18:42:53 发布
最新推荐文章于 2024-04-25 18:42:53 发布
阅读量7.5k 收藏 55
点赞数 13
分类专栏： FDTD学习笔记 文章标签： FDTD Lumerical 时域有限差分法 算法 米氏散射
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_45719141/article/details/124109718
版权

官网链接

我们模拟二维介电圆柱体的 Mie 散射。 将 FDTD 模拟的角散射结果与解析解进行比较。

在这里插入图片描述

一、结构设置

1.1 光源

全场散射场（total-field scattered-field,TFSF）源注入有限跨度的平面波，并将计算区域划分为两个区域。在源区内部同时存在入射场和散射场，而在源区外部只存在散射场。全场散射场源的工作原理是在其边界处，减去直接透射通过源区域或从平坦基板反射的任何光。因此，只有被完全包含在源区域内部的特征散射的光，才能穿过边界并传播到散射场区域。

1.2 分析组

总场监视器
在这里插入图片描述
Scirpt：

####################################################
# Total field
# This object can be used to calculate nano-particle
# cross sections. 
# # 这个对象可以用来计算纳米粒子的横截面。

# Output properties
# sigma: cross section 横截面
# lam: vector of wavelengths 波长向量
# f: vector of frequencies 频率向量
#
# Tags: FDTD total field monitor
#
# Copyright 2012 Lumerical Solutions Inc
################################################################

f=getdata("y2","f");   # get frequency data
lam=c/f*1e9;              # convert to wavelength in nm.

# Calculate power flowing outwards in each direction in 
# scattered field region. 
# Note minus sign for bottom and left monitors.
# 计算散射场区域每个方向向外流动的功率。
# 注意底部和左侧显示器的减号。
Ptop    =  transmission("y2");
if(havedata("y1")) { Pbottom = -transmission("y1"); }
else { Pbottom = Ptop; }
Pright  =  transmission("x2");
if(havedata("x1")) { Pleft   = -transmission("x1"); }
else { Pleft = Pright; }

# Calculate total power flowing outwards.
# power scattered (in Watts)
# 计算向外流动的总功率。
# 功率分散（以瓦特为单位）
Pscat = (Ptop + Pbottom + Pleft + Pright)*sourcepower(f);

# Calculate scattering cross-section
# 计算散射截面,也即面积
sigma=Pscat/sourceintensity(f);

havedata
在这里插入图片描述
sourcepower

在这里插入图片描述
在这里插入图片描述
sourceintensity

在这里插入图片描述
散射场监视器
在这里插入图片描述
Scirpt：

####################################################
# Scattered field
# This script calculates scattering cross-section
# and far field projection.
# # 这个脚本计算散射截面和远场投影。
# Input properties
# do far field: 0 for no, 1 for yes
#
# Output properties
# theta: vector of angles for far field projection 远场投影的角度向量
# E2_far: vector of far field intensities at angles theta 1 mm from cylinder
# 距离圆柱 1 毫米角度的远场强度矢量
# sigma: scattering cross section 距离圆柱 1 毫米角度的远场强度矢量
# lam: vector of wavelengths 波长向量
# f: vector of frequencies 频率向量
#  
#
# Tags: FDTD scattered field monitor
#
# Copyright 2012 Lumerical Solutions Inc
################################################################

f=getdata("y2","f");   # get frequency data
lam=c/f*1e9;              # convert to wavelength in nm.

# Calculate power flowing outwards in each direction in 
# scattered field region. 
# Note minus sign for bottom and left monitors.
Ptop    =  transmission("y2");
if(havedata("y1")) { Pbottom = -transmission("y1"); }
else { Pbottom = Ptop; }
Pright  =  transmission("x2");
if(havedata("x1")) { Pleft   = -transmission("x1"); }
else { Pleft = Pright; }

# Calculate total power flowing outwards.
# power scattered (in Watts)
Pscat = (Ptop + Pbottom + Pleft + Pright)*sourcepower(f);

# Calculate scattering cross-section
sigma=Pscat/sourceintensity(f);

###########################################################
# far field calculation
if(%do far field%)
{

  #define a vector of angles for far field calculation # 定义远场计算的角度向量
  theta=linspace(0,180,200);

  #calculate the X and Y vertices at a radius of 1mm from the cylinder
  #these are the locations where the far field calculation is done
  #计算距离圆柱1mm半径的X和Y顶点
   #这些是远场计算完成的位置
  r0=1e-3;
  x=r0*cos(theta*pi/180);
  y=r0*sin(theta*pi/180);

  #read the monitor frequencies #读取监视器频率
  far_fj=matrix(2); # find indices of wavelengths for far field proj. # 找到远场项目的波长指数,1um和2um。
  far_fj(1)=find(c/f,1e-6);  
  far_fj(2)=find(c/f,2e-6); 
  #allocate a matrix for the result 为结果分配一个矩阵
  E2_far = matrix(2,length(theta));

  # loop over all angles and far field positions
  # 循环所有角度和远场位置
  for(j=1:2){ # the frequency points

      tempE = farfieldexact("x2",x,y,far_fj(j))-farfieldexact("x1",x,y,far_fj(j))
             +farfieldexact("y2",x,y,far_fj(j))-farfieldexact("y1",x,y,far_fj(j));
	# tempE大小为200×3，即L×3
      # sum over all field components E2 = |Ex|^2 + |Ey|^2 + |Ez|^2
      # 对所有场分量求和 E2 = |Ex|^2 + |Ey|^2 + |Ez|^2
      E2_far(j,1:length(theta)) = sum(abs(tempE)^2,2); #第二个维度求和，即各场分量
 }
}

farfieldexact

在这里插入图片描述
在这里插入图片描述
数据mie_cylinder
在这里插入图片描述
数据以 3（行）X 1000 列的行/列格式组织
第一行是以度为单位的散射角
第二行是 lambda = 1 微米的散射幅度
第三行是 lambda = 2 微米的散射幅度

(注：列数太多了，所以没有完整显示为一行）

mie_theory_2d.lsf

###########################################################################
# Scriptfile: mie_theory_2d.lsf
#
# Description:
# This script analyses the results of mie_theory_2d.fsp
#
# Copyright 2009, Lumerical Solutions, Inc.
###########################################################################

do_far_field = 1;         # set 0 to not perform this test, 1 to perform the test
# 设置 0 不执行此测试，设置 1 执行测试
do_cross_sections = 0;    # set 0 to not perform this test, 1 to perform the test

select("scat");
set("do far field",do_far_field);
runanalysis;

###########################################################################
# Far field projections
if (do_far_field) {

  E2_far = getdata("scat","E2_far");
  theta = getdata("scat","theta");
  
  for(j=1:2){ # the frequency points
      E2_far(j,1:length(theta)) = E2_far(j,1:length(theta))/E2_far(j,1);
  }

  #load analytic result from file, there are 1000 points
  #从文件中加载解析结果，有1000个点
  mie_theory_data = readdata("mie_cylinder.txt");
  theta_analytic = mie_theory_data(1,1:1000);

  # plot data for the 2 wavelengths
  E2_analytic_1=interp(mie_theory_data(2,1:1000),theta_analytic,theta);
  E2_analytic_2=interp(mie_theory_data(3,1:1000),theta_analytic,theta);
  plot(theta,E2_far(1,1:length(theta)),E2_analytic_1,
             E2_far(2,1:length(theta)),E2_analytic_2,
             "angle (degrees)","scattered intensity (log scale)",
             "scattered intensity (log scale)","logplot");#"logplot"表示绘制对数图
  legend("FDTD 1 micron","Analytic 1 micron","FDTD 2 microns","Analytic 2 microns");
  # polar(theta*pi/180,E2_far(1,1:length(theta)),E2_analytic_1, E2_far(2,1:length(theta)),E2_analytic_2);
  # legend("FDTD 1 micron","Analytic 1 micron","FDTD 2 microns","Analytic 2 microns");
}




###########################################################################
# Cross sections
if (do_cross_sections) {
  ?"  Calculating scattering and absorption cross sections";

  # collect data on wavelength range
  f=getdata("scat","f");      # get freqency data
  lambda=c/f;                 # convert to wavelength
  r = getnamed("circle","radius"); # the radius of the mie particle

  # Calculate scattering cross-section
  sigmascat = getdata("scat","sigma");
  Qscat = sigmascat/(2*r);

  # Calculate absorption cross-section
  sigmaabs = -getdata("total","sigma");
  Qabs=sigmaabs/(2*r);

  # calculate the size parameter
  size_parameter = (2*r)/lambda;

  # Plot results
  plot(size_parameter,Qscat, Qabs ,"size parameter (d/lambda)","Mie efficiency","cross section");
  legend("Scattering","Absorption");
}

readdata

在这里插入图片描述
interp
在这里插入图片描述

二、结果

打开并运行名为 mie_theory_2d.fsp 的仿真文件。 模拟在下面创建 mpeg 电影。 全场散射场 (TFSF) 源用于产生入射到介质圆柱体上的平面波。 只有被圆柱体散射的光传播到 TFSF 边界之外。

在这里插入图片描述
mie_theory_2d_movie_TM

运行模拟后，运行名为 mie_theory_2d.lsf 的脚本文件。 脚本文件计算 2 个波长的一系列远场投影，并将结果与存储在文件 cylinder.txt 中的电介质圆柱体散射的分析数据进行比较。 请注意，监视器记录了 25 个波长，但分析组中的脚本选择了 1 和 2 微米的监视器数据。 结果图（对数刻度）如下所示。 分析和模拟结果在所有角度都非常吻合，超过 5 个数量级。

在这里插入图片描述
该脚本还可以计算散射和吸收截面。 吸收将为零，因为圆柱体是由没有吸收的电介质组成的。 要计算横截面，请编辑脚本并将 do_cross_sections 变量设置为 1。

子非鱼icon
关注
  • 13
    点赞
  • 踩
  • 55
    收藏
    觉得还不错? 一键收藏
  • 6
    评论
专栏目录
技术美术个人笔记（八）——瑞利散射与米氏散射
小谁的博客
08-10 1万+
散射是指光子与空气中的分子或是悬浮粒子弹性碰撞之后，改变行进方向。生活中常见的散射可概分为两种：瑞利散射（Rayleigh scattering）和米氏散射（Mie scattering），示意图如图。 图解：图中的横轴是粒子直径与波长长度的比值，纵轴是经过数据标准化后的散射强度。蓝色区域是粒子直径比波长小很多的瑞利散射，波长越短，散射越强，而且散射角度较大；黄色区域是粒子较大的米氏散射，散射强度和波长之间没有固定关系，且大部分散射沿着阳光原本行进的方向前进。 瑞利散射是指当光照到直径比光波长小很多的
利用FDTD软件仿真拓扑光子（五）-抗散射仿真与软件设置
yangzhao2014的博客
11-29 3392
本系列主要讲解如何利用lumerical公司的FDTD软件仿真拓扑光子绝缘体的能带结构。主要包括以下几方面的内容： 1）前言 2）光子晶体结构分析 3）能带结构仿真与软件设置 4）边缘态仿真与软件设置 5）抗散射仿真与软件设置 6）单向仿真与软件设置 7）其它仿真 1、 仿真目标 本系列的仿真工作是基于文献[1]，采用石墨烯型光子晶体构成拓扑界面，下图左红色和蓝色界面处所示...
金纳米-散射-画图代码_米散射_fdtd_FDTD消光_nanogold_金纳米颗粒_
10-01
纳米金颗粒散射，在FDTD环境下运行，包括一个画吸收，散射，消光的代码文件。
Lumerical官方案例、FDTD时域有限差分法仿真学习（六）——等离子体超材料吸收器（Plasmonic metamaterial absorber）
weixin_45719141的博客
03-20 1万+
在本例中，将根据参考文献 [6] 计算等离子体超材料吸收体的吸收特性。 这些器件的吸收机制主要受局域电磁共振的激发支配，并且高度依赖于顶部金属层的几何形状和介电层的厚度。
Comsol金属贴片建模与多极子展开分析：透反射计算与极子分析详细教程
最新发布
XKJQWJnJa的博客
04-25 388
本文介绍了使用Comsol进行金属贴片建模和多极子展开分析的方法，并结合透反射计算和极子分析对金属贴片的电磁性能进行评估和优化。通过本文的介绍，读者可以学习到金属贴片建模和多极子展开分析的基本方法，并了解到Comsol在金属贴片设计和优化中的应用。通过对金属贴片进行多极子展开分析，可以得到金属贴片的透射率、反射率等重要参数，对于优化金属贴片的性能具有重要意义。通过分析金属贴片的极子分布，可以了解金属贴片的电磁响应机制，为优化金属贴片的设计提供依据。Comsol金属贴片建模，金属体系多极子展开分析。
Ansys Lumerical | 米氏散射 FDTD
ueotek的博客
08-10 1003
如果颗粒在基板上，则必须修改分析的远场部分。Mie 效率定义为横截面与几何面积的比， πr2对于球体（3D）和对2r于圆（2D），并且通常相对于尺寸参数 （2πn1/λ)，其中n1是 FDTD 区域的背景指数，对于空气为1。下图显示|E|的平方在XY，XZ和YZ平面中，在最接近脚本中指定的“目标波长”的波长点穿过粒子中心。与其将结果归一化为源功率（对于理想平面波来说，这是无限的，因为它具有无限的范围），不如按源强度进行归一化。此外，较小的网格会产生更高分辨率的场轮廓，从而更好地解析金属界面附近的场。
[个人笔记]FDTD100
Temmie1024的博客
01-23 1万+
边界条件 名称 含义 PML 吸收边界条件 periodic 周期性边界条件
FDTD散射复现
weixin_63709084的博客
03-24 2247
本文用FDTD Solutions复现了粒子散射问题
基于FDTD的雷达散射截面积（RCS）分析
congcongnanian1221的博客
09-10 4020
FDTD作为一种数值计算方法，需要将计算空间离散成适用于FDTD计算的Yee网格。根据远区场的定义可知，远场区的距离是近场区的几十甚至上百倍。由于计算机的内存有限，这就给FDTD计算造成了很大的困难。目前主要是利用等效原理结合近远场外推技术来计算远场。通过在计算目标周围形成一个包络面，通过包络面上的等效电磁流Js和Ms来计算出散射远场的信息。 电磁流计算如下表： 外推长方体...
电磁场与仿真软件(38)
weixin_46649493的博客
03-22 764
FDFD
FDTD solutions——TFSF（全场散射场）光源及斜入射
热门推荐
Temmie1024的博客
09-24 3万+
初始建立仿真文件： 1.建立物理结构  几何结构  材料特性 2.定义仿真区域  边界条件 3.设置激励源（光源） 4.设置监视器  至少使用一个时间监视器  分析组  只记录需要的数据 运行与仿真： 1.运行仿真文件  分布式并行运算  并发式并行运算 2.参数扫描与优化 3.结果分析  视觉化器visualize  使用script文本程式 光源使用原则 Plane wave：平面波，周期性结构（或非周期结构，但PML阶段产生衍射时）；如各类光栅、超导材料等。 Gaussian/NA source：高斯
光学教学中的FDTD方法和模型
02-25
本文提出了时域有限差分法（FDTD）作为光学教育的教学方法。 同时，提出了一种基于FDTD算法的仿真软件FDTD解决方案，作为一种新工具，可以帮助初学者建立光学模型并分析光学问题。 FDTD算法的核心是将与时间相关的麦克斯韦方程组离散为空间和时间偏导数，然后模拟电子脉冲与理想导体或半导体之间相互作用的响应。 由于电磁场的求解是时域的，因此减少了存储器的使用，可以很容易地获得宽带的仿真结果。 因此，在光学教学中推广FDTD算法是有效的。 FDTD使我们能够设计，分析和测试用于波传播，散射，反射，衍射，极化和非线​​性现象的现代无源和非线性光子组件（例如生物粒子，纳米粒子等）。 不同的FDTD模型可以帮助教师和学生解决光学教育中的几乎所有光学问题。此外，FDTD解决方案的GUI对初学者非常友好，学习者可以很快地掌握它。
Mie散射计算
09-10
matlab下 用mie散射 计算球形粒子 散射效率 消光效率 或者散射截面 消光截面
FDTD方法在柱面波及半空间有耗介质散射问题中的应用
04-07
FDTD方法在柱面波及半空间有耗介质散射问题中的应用
米氏理论经典版--包括计算方法
12-17
经典的米氏理论介绍学习，强烈推荐 介绍了米氏理论的公式获得和计算方法。 如果你想计算散射截面，吸收界面，可以用这个教材。
FDTD.rar_FDTD radar_fdtd MATLAB_radar cross section _三维 FDTD 散射
09-19
用c++语言，编写了基于FDTD的三维散射体雷达散射截面的计算程序。程序注释清楚。稍微改动必要的数据就可以直接使用。已经多次验证了程序的正确性。
专业FDTD软件lumerical破解版,时域有限差分法
01-20
利用时域有限差分法的专业光学建模工具----lumerical破解版，32位\64为系统具有
关于有限时域差分法，Lumerical-FDTD Solutions的中文指导手册
09-02
关于有限时域差分法，Lumerical-FDTD Solutions的中文指导手册，非常详细，每一个图标的意义和用法都有解释说明，从软件各个图标含义，到使用步骤和建模流程，以及背后理论解释，都有详解。从事光学建模仿真相关学习...
时域有限差分法对平面TE波的MATLAB仿真.zip_时域有限差分 fdtd 吸收边界
07-15
基于时域有限差分TE电磁波仿真,含MATLAB代码
时域有限差分法信道建模matlab仿真
08-18
时域有限差分法（Finite Difference Time Domain, FDTD）是一种数值求解电磁波传播问题的方法，可以用于信道建模的MATLAB仿真。以下是一个简单的示例代码，用于信道建模的FDTD仿真： ```matlab % 信道参数 c = 3e8; % 光速 f0 = 2.4e9; % 中心频率 lambda = c/f0; % 波长 dx = lambda/20; % 空间步长 dy = lambda/20; % 空间步长 dt = dx/(2*c); % 时间步长 % 仿真区域尺寸 x_size = 100*lambda; y_size = 100*lambda; t_end = 200*dt; % 创建网格 nx = round(x_size/dx); ny = round(y_size/dy); nt = round(t_end/dt); % 初始化电场和磁场 ez = zeros(nx,ny); hy = zeros(nx,ny); % 更新步骤 for t = 1:nt % 更新磁场 for i = 1:nx-1 for j = 1:ny-1 hy(i,j) = hy(i,j) + (ez(i+1,j) - ez(i,j))/dx; hx(i,j) = hx(i,j) + (ez(i,j+1) - ez(i,j))/dy; end end % 更新电场 for i = 2:nx-1 for j = 2:ny-1 ez(i,j) = ez(i,j) + (hy(i,j) - hy(i-1,j))/dx; end end % 激励源 ez(nx/2,ny/2) = sin(2*pi*f0*t*dt); % 绘制结果 imagesc(ez'); title('FDTD信道建模仿真'); xlabel('X轴'); ylabel('Y轴'); colorbar; drawnow; end ``` 这是一个简单的二维FDTD仿真，使用了时域的差分方法来求解电磁场的传播。你可以根据需要修改参数和边界条件，来进行不同场景下的信道建模仿真。注意，在实际应用中可能需要考虑更多的物理效应和精细调节，这里只是一个简单的示例。 希望对你有所帮助！

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 6
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值