FDTD_mie散射_项目研究（1）
mie散射
项目翻译
计算由平面波激发的纳米粒子的散射和吸收截面、局部场增强以及远场散射分布（米氏散射）。将截面和远场结果与解析解进行比较，以验证模拟的准确性。
纳米粒子的散射特性一般从场增强、横截面和远场分布等方面来描述。这个例子展示了如何从单个FDTD模拟中获得这些结果。
运行和结果
局部场增强 电磁场与纳米粒子的相互作用可在粒子表面产生强烈的场增强。频域场监视器可直接测量局部场增强。以下各图分别展示了在最接近脚本中指定的“目标波长”的波长点处，通过粒子中心的 XY、XZ 和 YZ 平面上的 |E|²。可以注意到，TFSF 源的边缘在图像中表现为颜色的突然变化。源内部的场为“总”场（即入射场 + 散射场），而源外部仅可见“散射”场。
吸收和散射截面
吸收截面（总吸收功率除以入射光束单位面积的功率）由位于 TFSF 源内部的分析组计算得出。该分析组测量流入粒子的净功率，并将其归一化为源强度，从而得出吸收截面。同样，散射截面由位于 TFSF 源外部的分析组计算得出。
根据定义，对于三维模拟，截面的单位为平方米；对于二维模拟，单位为米。截面测量值通常会根据散射体的大小进行归一化处理，如以下图示所示。米氏效率被定义为截面与几何面积（对于球体为πr²，对于圆为2r）的比值，并通常以尺寸参数（2πn₁/λ）为横坐标绘制，其中n₁为 FDTD 区域的背景折射率，在空气中为 1。
将 FDTD 结果与通过 mie3d 脚本获得的解析解进行了比较。两者结果之间的差异较为明显，需要对模拟设置进行一些改进。这将是接下来关于收敛性测试部分的主题。
远场角散射
在大多数散射实验中，对散射场（辐射模式）的测量是在远离散射体的位置进行的，相对于所考虑的波长尺度而言。 "scat_ff" 监测器返回远场中的散射场分布。以下极坐标图展示了 X-Y、X-Z 和 Y-Z 平面中远场的散射场。每个图都包含两种颜色的线：蓝色表示 FDTD 模拟结果，绿色表示来自 mie3ds12 脚本命令的解析结果。第一张图说明了每个平面中极角的定义方式。
重要的模型设置
模型设置脚本
模型对象中的设置脚本用于设置网格大小、模拟范围和粒子位置。该脚本是一种便捷的方式，可确保模拟区域、网格覆盖区域、源、散射体和吸收监视器的位置正确。例如，重要的是 TFSF 源位于散射体和吸收监视器之间，且这些对象之间至少有两个网格单元的距离。这些对象的位置必须通过设置脚本进行设置。其他属性，例如模拟时间，可以直接在对象中修改。
TFSF 源
TFSF 源专门针对此类情况而设计，即非周期性物体被平面波照射。它通过将散射场与入射场分离，使纳米粒子的散射分析变得简单直接。为了使散射分析正常工作，确保散射体完全位于 TFSF 源内至关重要。
使用 TFSF 源进行功率归一化
使用 TFSF 源进行功率归一化可能会令人感到困惑。与其将结果归一化到源功率（对于理想的平面波，由于其范围无限，源功率也是无限的），不如通过源强度进行归一化。这会导致功率测量值以横截面类型的单位返回。有关更多信息，请参阅 TFSF 源页面的功率归一化部分。
“abs”和“scat”分析组
在三维模拟中，横截面分析组由六个二维监视器组成，形成一个封闭的盒子，用于测量进出该盒子的净功率。这些分析组的位置非常重要。“abs”分析组用于测量吸收功率，必须完全位于 TFSF 源内部但粒子外部。“scat”监视器必须完全位于 TFSF 源外部。
网格覆盖区域
在涉及金属的模拟中，网格覆盖区域通常用于更精确地确定金属界面的位置，尤其是在曲面的情况下。在此模拟中，网格覆盖区域设置得足够大，不仅覆盖了金球，还覆盖了整个 TFSF 区域。这样做是有意为之，因为 TFSF 源在均匀网格区域中效果最佳。另外，请注意网格尺寸会影响总和散射监视器与源的距离。最好在源和监视器之间保持至少两个网格单元的距离，以避免放置在灰色阴影的源注入区域内的监视器返回不物理的结果。请注意，这些条件由“模型”设置脚本强制执行。有关更多信息，请参阅有关监视器与其它对象距离的页面。
对称性
此模拟在 X 和 Z 两个维度上都具有对称面。为了将模拟时间和内存减少四分之三，将 X 方向的最小边界条件设置为对称，将 Z 方向的最小边界条件设置为反对称。请注意，只有当粒子和源都具有必要的对称性时，才能使用对称性。
基于您的设备参数更新模型的说明
该模拟文件已参数化，以便更轻松地设置模拟。当前模板使用球形粒子，但也可用于任意形状的粒子或多个粒子。在“模型”中指定参数后，其余模拟对象的大小将自动调整。
设置源波长范围和偏振。设置纳米粒子的材料或折射率。在“模型”中设置纳米粒子的跨度和位置、网格覆盖的网格大小以及模拟跨度。源和“abs”/“scat”分析组将自动间隔最多两个网格单元，纳米粒子将完全被“abs”分析组包围。当模拟非球形粒子或多个粒子时，可能需要更新边界条件以匹配新结构的对称性。相关的脚本文件也需要修改，以校正几何区域和散射体的尺寸参数。