折腾光纤模型的手记

comsol仿真-W型光子晶体光纤色散与损耗分析效果展示最近在实验室被导师催着搞光子晶体光纤的仿真W型结构这种带双包层设计的玩意儿确实有点意思。作为COMSOL萌新边啃说明书边试错折腾一周终于把色散曲线和损耗谱给整明白了。先说建模这个重头戏画六角晶格的时候手滑设置错了对称性结果跑出来的模式场跟抽象画似的。老老实实重新建了个参数化模型空气孔直径d2.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时在CL波段实现了近零色散平坦区同时保持损耗低于0.5dB/km。不过导师看完幽幽来了句这参数加工难度堪比绣花... 得看来下周得接着调工艺兼容性参数了。