COMSOL光学波导传输仿真 光纤等波导的三维弯曲 模场分布 波束包络方法 FDTD计算模式弯曲损耗，模场分析

三维弯曲的光纤仿真总是让人头大——传输模式扭曲得像麻花，损耗突然飙升，模场分布更是难以捉摸。最近在COMSOL里折腾波导弯曲仿真时，发现波束包络法（Beam Envelopes）比传统有限元更适合处理这类问题。

先看个典型场景：半径5mm的90度光纤弯曲结构。用COMSOL建模时，几何部分可以用参数化螺旋线实现：

% 螺旋线参数化方程
radius = 5e-3;  // 弯曲半径
pitch = 0;      // 直波导设为0
theta = linspace(0, pi/2, 50);  // 90度弯曲
x = radius*cos(theta);
y = radius*sin(theta);
z = pitch*theta/(2*pi);

这里有个坑——直接画圆弧会导致端口模式失配。我的解决方案是先建立直波导段，再用"变形几何"功能实现平滑弯曲过渡，这样初始模场更容易收敛。

材料设置需要特别注意各向异性。光纤纤芯通常用二氧化硅，但弯曲会导致等效折射率变化。通过自定义坐标系实现各向异性材料：

n_eff = neff + (x/R)*n2;  // R为弯曲半径，n2为非线性系数
material = mphcreate('mat1','Geometry','comp1','CoordinateSystem','sys1');

波束包络法的核心在于用缓慢变化的包络函数代替快速振荡的场分量。在物理场设置中启用"Wave Envelope"模块，选择"User defined"波矢方向。这里有个实用技巧——用"mphinterp"函数提取初始模式的波矢量作为输入：

k0 = 2*pi/lambda;
k_vector = [0, 0, k0*n_core];  // 初始波矢方向

求解时建议开启"渐进扫描"功能，先计算直波导模式，再逐步增加弯曲角度。遇到不收敛的情况，试试调整"数值耗散"参数，0.1-0.3之间往往有奇效。

后处理阶段重点关注模场畸变和功率损耗。COMSOL的线积分功能可以直接计算弯曲损耗：

P_in = integrate(emw.PoavZ,2);  // 输入功率
P_out = integrate(emw.PoavZ,5); // 输出功率
loss = -10*log10(P_out/P_in);   // 单位dB

对比FDTD方法时发现，虽然FDTD能更好处理大角度弯曲的辐射损耗，但计算三维结构需要恐怖的网格量。波束包络法在内存占用上优势明显——同样的光纤模型，FDTD需要32GB内存，而COMSOL只用8GB就能跑。

最后分享个实用脚本，自动提取模场分布并生成GIF动画：

for ii = 1:num_steps
    mphloadframe(model,ii);
    F = mpheval(model,'emw.E');
    contourf(F.x,F.y,abs(F.E));
    writeAnimation(gcf,'mode_evolution.gif');
end

折腾了半个月的教训：弯曲半径小于3mm时，传统LP模式近似完全失效，必须用全矢量仿真；渐变折射率光纤比阶跃型更难收敛，建议用"人工阻尼层"吸收反射波。