折腾光纤模型的手记

comsol仿真-W型光子晶体光纤色散与损耗分析效果展示

最近在实验室被导师催着搞光子晶体光纤的仿真，W型结构这种带双包层设计的玩意儿确实有点意思。作为COMSOL萌新，边啃说明书边试错，折腾一周终于把色散曲线和损耗谱给整明白了。

先说建模这个重头戏，画六角晶格的时候手滑设置错了对称性，结果跑出来的模式场跟抽象画似的。老老实实重新建了个参数化模型，空气孔直径d=2.4μm，晶格常数Λ=4μm。关键代码段长这样：

model.param.set('d', '2.4[um]'); model.param.set('Lambda', '4[um]'); geometry = model.geom.create('geom1', 3); hexagon = geometry.create('hex1', 'Hexagon'); hexagon.set('type', 'circumradius'); hexagon.set('r', 'Lambda');

这里设置晶格参数时容易踩的坑是单位换算，COMSOL默认用国际单位制，但咱们工程上习惯用微米作单位。后面那个Circumradius模式比Edge length模式更方便控制孔间距，亲测有效。

comsol仿真-W型光子晶体光纤色散与损耗分析效果展示

跑模式分析时发现边界条件设置是门玄学。原本随便选了个散射边界，结果有效折射率虚部异常得离谱。后来换成完美匹配层(PML)才稳定下来，具体厚度得试个三五次——太薄了吸收不够，太厚了计算量爆炸。建议用参数扫描功能批量测试，类似这种操作：

study = model.study.create('param_sweep'); study.set('param', 'pml_thickness'); study.set('vals', [0.5,1,1.5,2]);

仿真跑通后处理数据才是真战场。用内置的ewfd.neff提取有效折射率时，注意要取基模的实部计算色散。色散系数D的计算公式得自己写脚本处理：

lambda = np.linspace(1.3,1.6,100) # 波长范围 D = - (lambda/c) * np.gradient(neff_real, lambda) # 色散系数公式 plt.plot(lambda, D*1e6) # 单位转ps/(nm·km)

这个梯度计算容易翻车，建议用中心差分法代替简单差分。有次偷懒没处理噪声，出来的曲线跟心电图似的，被组会同门笑了半小时。

损耗方面重点看材料吸收和泄露损耗的竞争关系。在1.55μm窗口，W型结构的双包层设计确实给力。对比传统单包层结构，弯曲损耗降低了约40%。不过当空气填充比超过75%时，模式约束变弱导致损耗反升，这个转折点需要重点关注。

最后秀个成果图（此处脑补色散曲线和损耗谱对比图）。发现当第二包层空气孔直径控制在1.8μm时，在C+L波段实现了近零色散平坦区，同时保持损耗低于0.5dB/km。不过导师看完幽幽来了句："这参数加工难度堪比绣花..." 得，看来下周得接着调工艺兼容性参数了。