3、数值仿真方法概述：FDTD、FEM、BPM、EIM
做集成光学仿真，说白了就是跟麦克斯韦方程组打交道。但解析解只存在于教科书里，真实波导结构——弯曲、锥形、耦合区——都得靠数值方法硬算。我个人习惯把这四种方法看作一套工具箱：什么场景用哪把扳手，心里得有数。
3.1 时域有限差分法（FDTD）
FDTD 的思路很暴力——直接把时空离散化，用差分代替微分，一步步推进电磁场。我在做 MMI 分束器 优化时，FDTD 是首选，因为它能捕捉所有后向反射和倏逝波细节。
| 特性 | 说明 |
|---|---|
| 核心原理 | Yee 网格 + 蛙跳迭代，时间步进求解 |
| 优势 | 宽频带、全波、支持非线性/色散材料 |
| 代价 | 网格尺寸需 < λ/10，高 Q 结构收敛极慢 |
| 典型场景 | 光栅耦合器、环形谐振器、超表面 |
3.2 有限元法（FEM）
FEM 把计算域切成三角形或四面体单元，用基函数逼近场分布。它的网格可以局部加密——在模场变化剧烈的区域（比如脊波导的尖角）自动细化，平坦区域用大网格。我习惯用 FEM 做 模式分析 和 弯曲波导 损耗计算。
为什么 FEM 适合弯曲结构？因为它的非结构化网格能完美贴合弯曲边界，而 FDTD 的矩形网格会产生阶梯近似误差。嗯，这里要注意：FEM 的矩阵组装和求解很吃内存，3D 模型动辄几十 GB。
3.3 光束传播法（BPM）
BPM 假设光主要沿一个方向传播（比如 z 轴），忽略后向反射。它把 3D 问题降维成 2D 切片迭代，速度极快。你想想看，一个 1mm 长的锥形波导，FDTD 要算几小时，BPM 几分钟就搞定了。
但 BPM 有硬伤——它不能处理 大角度弯曲（>10°）和 强耦合 场景。我记得有一次用 BPM 仿真 Y 分支，结果跟实验对不上，后来才发现后向反射被忽略了，改用 FDTD 才复现出干涉条纹。
| BPM 变体 | 适用场景 | 限制 |
|---|---|---|
| 标量 BPM | 弱导波导、长距离传输 | 忽略偏振 |
| 矢量 BPM | 高折射率差波导（如 SOI） | 计算量翻倍 |
| 广角 BPM | 弯曲波导、锥形结构 | 角度仍受限 < 30° |
3.4 模式求解法（EIM）
EIM 也叫有效折射率法，是这四种方法里最「轻量」的。它把 2D 波导截面分解成两个 1D 问题，先算垂直方向的有效折射率，再算水平方向的模式。说白了，就是降维打击。
我在做 阵列波导光栅（AWG） 的初步设计时，EIM 是首选。它能快速给出波导宽度、间距对有效折射率的影响趋势，帮我把设计空间从几百个参数缩小到十几个。但 EIM 的精度有限，尤其对高折射率差波导（如 Si 波导），误差可能达到 10%。
- 初步设计阶段，需要快速扫参
- 弱导波导（如 SiO₂ 波导），精度足够
- 与 FDTD/FEM 配合，先粗后精
3.5 方法对比与选择策略
四种方法各有脾气。我整理了一张对比表，方便你快速决策：
| 方法 | 精度 | 速度 | 内存 | 适用维度 | 典型误差源 |
|---|---|---|---|---|---|
| FDTD | ★★★★★ | ★☆☆☆☆ | 高 | 2D/3D | 阶梯近似、数值色散 |
| FEM | ★★★★★ | ★★☆☆☆ | 高 | 2D/3D | 基函数阶数不足 |
| BPM | ★★★☆☆ | ★★★★★ | 低 | 2D/3D | 忽略后向反射 |
| EIM | ★★☆☆☆ | ★★★★★ | 极低 | 1D/2D | 降维近似误差 |
我个人习惯的流程是：先用 EIM 快速扫参，锁定几个候选结构；再用 BPM 或 FEM 做精细优化；最后用 FDTD 验证关键性能指标。这样既保证了效率，又不会漏掉物理细节。