做COMSOL仿真的人多少都会走到“复现论文”这一步。尤其像“光子晶体光纤”这种微结构光波导想完全靠自己去构想结构参数基本不可能只能从已发表的文献里找结构图、找折射率分布、找模式场图。我这次接到的复现任务目标很明确复现论文里左侧的仿真图一个是SPR传感器一个是三芯分束器。这两个器件的物理机制完全不同一个靠表面等离子体共振的损耗谱变化做传感一个靠多芯间的模式耦合做分束但它们在COMSOL里的实现路径高度重合几何建模、材料赋值、模式分析、后处理出图。这篇博文把我从拆解题目到最后出图的全过程整理出来重点讲清楚每一步为什么要这么做以及哪些地方容易踩坑。想动手复现同类论文的可以直接照着参数抄。1. 复现之前先把这篇论文的“图”拆明白1.1 标题里的信息量到底有多大“COMSOL光学仿真 光子晶体光纤 论文复现图是仿的一个SPR传感器和一个三芯分束器图左原文”这句话看起来只是零散的笔记但拆开来看其实包含了三个明确的信息第一仿真工具是COMSOL不是RSoft也不是Lumerical第二仿真对象属于光子晶体光纤范畴换句话说几何结构里有周期性排列的空气孔功能依赖微结构引入的等效折射率差第三需要复现的具体器件有两个SPR传感器和三芯分束器而且目标非常明确——把论文左侧的图“仿”出来。在光学仿真领域SPR传感器和三芯分束器是两个很有代表性的研究热点。SPR传感器通常是在光子晶体光纤的某个空气孔内壁镀一层金属膜利用光纤传导模与金属表面等离子体波之间的共振耦合实现对周围环境折射率的灵敏检测。三芯分束器则是在同一根光纤里设计三个纤芯利用相邻纤芯间的倏逝场耦合实现光功率在两个或多个通道之间的分配。两者的共同点是物理图像复杂、参数敏感、对网格精度要求高非常适合用COMSOL的波动光学模块来求解。这里我也直接说一句如果只是想把论文里的示意图照着画出来那用CAD或者PPT就够了。但如果要复现的是“仿真结果图”——比如模式场分布图、有效折射率曲线、损耗谱曲线——那就要把整个仿真流程跑通拿到真实的计算数据后再出图。题目里说的“图左原文”意思就是原论文的仿真结果在页面上靠左复现时要做到和它在数值、颜色分布、模式形态上对得上。1.2 复现论文图本质是复现什么数据很多人觉得复现论文图就是把那张图重新画一遍。但做了一次就知道真正的复现其实分三个层次。第一层是复现几何结构。这个相对容易从论文里提取空气孔直径、孔间距、层数、纤芯位置就可以了。第二层是复现物理结果比如SPR传感器在某个波长下的模式损耗峰值或者分束器在一定传播距离下的功率分配比这需要把材料参数、边界条件、求解设置全部对齐。第三层才是复现视觉效果包括模式场图的光强分布、颜色映射、坐标轴范围。这三个层次里第一层是基础第二层是灵魂第三层则是论文评审和读者最容易感知的部分在实际操作中最容易被忽视。以SPR传感器为例论文里通常会给出“限制损耗随波长变化”的曲线以及共振波长处的模式电场分布图。要复现这个结果你不仅要搭建出带金属膜的光子晶体光纤截面结构还要在波长扫描中逐步求解模式的有效折射率提取模式损耗再找到损耗最大的波长点在这个波长下重新计算并导出电场模分布。整个过程本质上是在复现一组“数据映射关系”结构参数 → 有效折射率 → 损耗 → 场分布。三芯分束器也是类似只是把“损耗”换成了“耦合长度”把“共振波长”换成了“拍长周期”。所以决定复现成败的关键不是你会不会按按钮而是你有没有把论文里隐藏的参数条件推断出来。比如SPR传感器的金属膜厚度论文正文里可能只写了一个范围真正出图用的却是某个确定值你只能靠反复试算去反推。1.3 COMSOL版本与模块的前置准备COMSOL做这类光学仿真需要用到的是Wave Optics Module波动光学模块核心物理场接口是“电磁波频域ewfd”。如果你只想算本征模式可以直接用“模式分析”研究步骤如果想看传播方向上的功率变化也可以用“边界模式分析”加“频域传播”的方式。这里我习惯的做法是截面结构用二维模型边界模式分析求解有效折射率和模式场分布如果需要沿光纤长度方向的演化再用二维模型得到各个模式的耦合系数手动计算或者用耦合模方程处理。版本方面我目前常用的是COMSOL 6.2但复现过程中也遇到用户在5.6版本上操作。不同版本在材料库、求解器默认设置上会有一些细微差别尤其是“模式分析”的搜索基准折射率这个参数老版本如果没有设置好很容易搜出完全不同的模式。建议在开始复现前先确认自己用的版本翻阅一下论文的提交时间尽量选择与论文结果发布时间接近的版本否则有的材料库数据会对不上。另外一个容易被忽略的地方是COMSOL的“案例库”功能。软件自带里有一个“Surface Plasmon Resonance Sensor”的案例还有光纤模式分析的案例这些都可以作为起点。我第一次做SPR传感器复现时就是从案例库里拉了一个金属-介质界面的模型然后修改成光子晶体光纤结构这样材料库和边界条件的默认设置都有保障减少了很多低级错误。2. 光子晶体光纤SPR传感器从几何结构到损耗谱复现2.1 光子晶体光纤的几何建模光子晶体光纤PCF的核心几何特征是周期性排列的空气孔。传统结构里空气孔以正六边形晶格排列中心位置去掉一个或多个空气孔形成纤芯。对于SPR传感器通常会在靠近纤芯的某个空气孔内壁沉积一层金属膜金、银或氧化铟锡ITO然后在金属膜内填充待测分析物。建模时第一步是在COMSOL中定义一个二维平面用“参数”节点设置空气孔直径d、孔间距Λ也叫节距、包层空气孔层数N、金属膜厚度t_metal等关键参数。比如常见的结构是Λ2.0μmd1.0μm金属膜厚度t40nm。这里千万注意单位换算COMSOL默认单位是米m你要是直接输入2它会被当成2米后面网格剖分和求解会直接崩溃。我习惯把所有尺寸参数先换算成微米为单位并在参数表里注明单位。第二步是画空气孔阵列。这里我用的是“正六边形阵列”的方式先在工作平面上画出中心的六个最近邻空气孔位置再按层数一圈一圈往外铺。不建议手动一个个画圆而是用COMSOL的“阵列”功能或者直接用脚本传入坐标数组。手画十几个孔还能忍几十个孔绝对会想哭而且手画还容易导致几何对象彼此重叠后续布洛赫周期边界或网格剖分都会出问题。第三步是挖掉中心的一个或多个空气孔形成纤芯区域。然后选择SPR检测孔在它周围画一个环形区域用来表示金属膜层。这个环形区域在COMSOL里可以用“圆环”对象实现内半径是检测孔半径外半径是内半径加金属膜厚度。完成金属环形域后再画一个稍大的圆表示待测分析物区域最外层还要再画一个圆作为完美匹配层PML的基础域。在实际操作中很多人画完几何后直接设置物理场结果模式分析时找不到想要的模式回头一看才发现空气孔之间有很多微小的重叠或者缝隙导致模式被限制在了奇怪的位置。我的建议是几何画完后使用“几何”工具栏里的“检查”功能把所有对象做一次布尔并集或者差集操作再查看最终域数目是否与预期一致。2.2 材料参数与金属膜层建模材料参数是非线性光学仿真里最容易“差之毫厘谬以千里”的部分。PCF的基体材料一般是石英玻璃二氧化硅SiO₂其折射率可以用Sellmeier公式计算COMSOL材料库里有内置数据。但要注意论文里可能用的是纯石英的折射率而你手头材料库里的“Silica glass”在可见光和近红外波段的数据可能有偏差最好还是自己用Sellmeier公式做一个插值函数输进去。对SPR传感来说金属膜层的介电常数是最关键的参数。金、银在光波段的介电常数通常用Drude模型或Johnson与Christy的实测数据描述。COMSOL的“材料库-光学-金属”里提供了金和银的光学常数但在某些波长范围内数据不全或者插值后介电常数虚部偏大导致SPR损耗峰位置漂移。建议的做法是从文献里找到所用的金属介电常数表达式比如金用ε(ω) ε∞ - ω_p² / (ω(ω iω_c))在COMSOL中这个表达式可以直接写成随角频率变化的解析函数赋给金属域。很多论文会明确给出ε∞、ω_p、ω_c这几个参数直接抄过来用就可以。如果论文没给就用材料库数据但不修改并在笔记里记下“此处依赖材料库结果可能与原图有偏差”。注意膜层的建模方式还有一层讲究。金属膜厚度通常是纳米量级例如40nm而整个光纤截面是几十微米几何尺度差了上千倍。如果用自由三角形网格去剖分这层膜网格会非常密集计算量爆炸。我通常会在金属环形域上用“映射”网格径向划分2-4层周向按角分辨率划分再将外部石英区域接上自由三角形网格。这样既保证了膜层内的高精度又避免了全局网格过密。2.3 边界条件与PML的合理配置COMSOL仿真中边界条件决定了你计算的是“理想孤立结构”还是“带有反射边界干扰的结构”。对于光子晶体光纤的模式分析最外层不能直接使用默认的边界条件因为光场虽然被限制在纤芯附近但包层里总有倏逝场泄漏直接截断会产生不物理的反射。这里的标准做法是在最外层添加完美匹配层并要求PML区域足够厚通常取最大波长的1/4到1/2。我一般将PML设置为一个环形域宽度约为工作波长的1/3材料属性用内置的PML特征自动设置。在“电磁波频域”界面中需要把PML域选上并在“完美匹配层”节点中确认厚度方向设置有误。如果你把PML设置错了方向求解结果通常会“不收敛”或者在边界处出现异常的高电场。对于周期性结构的光子晶体光纤如果只取了中心的一部分空气孔来分析还可以在包层外边界设置周期性边界条件但这只适用于完整的周期结构。SPR传感器因为引入了非周期的金属膜层和分析物通道整个结构不再严格周期这种情况下用PML代替周期边界反而更符合实际也更稳定。这里我必须强调一个很多人踩过的坑模式分析中如果PML设置不当搜索出来的模式可能是“PML模”一种被边界捕获的伪模式而不是光纤导模。判断方法很简单看模式有效折射率实部是否在包层有效折射率附近或之下或者看电场分布是否集中在PML区域。如果电场图里大部分能量集中在边界赶紧回到PML设置里去检查参数。2.4 模式分析与结果提取损耗谱怎么算模式分析是SPR传感器仿真的核心求解步骤。在COMSOL中你可以用“研究-模式分析”求解一个特征值问题得到该截面上的本征模式及其有效折射率n_eff。SPR传感器的检测原理是当某一波段的表面等离子体波与光纤基模满足相位匹配条件时基模能量会大量耦合到金属-分析物界面上的等离子体模式表现为该模式的有效折射率虚部迅速增大。具体数值关系是限制损耗 α 8.686 × (2π/λ) × Im(n_eff) dB/m。这里的λ是波长单位需要统一。在COMSOL中模式分析完成后可以直接在“全局计算”里输出实部和虚部然后在后处理中按公式计算损耗。实际操作时需要扫描波长。比如论文里给出了波长从0.6μm到1.1μm范围内的损耗谱你就需要建立一个参数化扫描依次改变波长并重新运行模式分析。这里要特别留心“模式搜索基准”的设置。因为波长变化后有效折射率也会变化如果搜索基准始终不变可能会锁定到高阶模。我的做法是在波长扫描中把“搜索频率附近的模式数”设置为6-8个然后把每轮的基模数据手动挑出来。也可以利用COMSOL的“模式跟踪”技巧在扫描过程中保持同一个模式序号但前提是模式间没有发生交叉。提取完损耗谱后要和论文曲线进行对比。首次复现时曲线形状对得上但共振波长偏移了可能有三个原因金属膜厚度偏差、金属介电常数表达式不同、分析物折射率设置不同。逐一排查后我那次偏移的原因是膜厚度用了30nm而论文实际用的是45nm调整后峰位立刻对齐。做这类复现耐心反推参数是常态第一次就对上的概率基本为零。3. 三芯分束器耦合模式与超模分析3.1 三芯结构的几何构建与参数值三芯分束器和SPR传感器的最大区别是它不需要金属膜整体就是纯介质波导依赖的是多个芯子之间的耦合。结构上中心芯和两个对称分布的外芯之间的距离、芯子半径、折射率差共同决定了耦合强度和分束比。采用COMSOL建模时第一步同样是创建空气孔阵列只是这次不去掉中心一个孔而是去掉三个孔——通常是一个中心孔加对称的两侧孔形成一字型或三角形排列的三芯结构。这三个纤芯周围的空气孔构成包层将光限制在三个纤芯里。常见参数Λ2.0μmd0.8μm三个缺陷孔中心的间距可以通过去掉的孔的数量控制比如中心芯与旁边芯相隔两个孔距。在几何实现上我习惯把所有空气孔用“圆形”对象参数化给定孔的编号和坐标。这里有一个很实用的小技巧用COMSOL的“参数化曲线”功能把孔的位置坐标用一个矩阵或循环生成避免反复手动输入坐标。如果你使用的是新版COMSOL还可以借助外部参数扫描和阵列功能一次生成整个包层。但需要注意COMSOL的“阵列”生成的是对象副本如果每个圆还要做布尔并集网格剖分前必须先合并几何。3.2 超模分析三芯波导为什么有四个特征模式这个环节是三芯分束器复现的核心知识点。三芯波导本质上是一个耦合波导系统当三个芯子靠得很近时单个芯子的模式不再是独立本征态系统的本征模式变成了“超模”supermodes。对于三个完全相同的芯子超模的组合方式是一个对称模symmetric mode所有芯同相位两个反对称模anti-symmetric modes还有一个在特定排列下会退化的模式。不同的超模具有不同的有效折射率传播常数也有差异这才导致了光束在传输过程中出现周期性功率转移的现象——也就是分束器的物理基础。在COMSOL中做超模分析需要把模式分析的“所需模式数”设置为足够多。比如你期望看到基模和第一组高阶模至少要设置8-12个目标模式。最重要的一个操作是在模式分析设置里把“搜索有效折射率基准”设置为略低于纤芯材料折射率的值。因为超模的有效折射率通常介于包层材料折射率和芯层折射率之间如果基准值设置太高或太低求解器会漏掉关键模式。实际复现时我还发现一个规律三芯分束器的超模电场分布图与单芯光波导完全不同。普通单模光纤的基模是一个圆对称光斑而三芯超模会出现明显的多峰结构。例如对称模是三个等强度光斑同相位反对称模则会出现相位反转中间芯的电场与两侧芯的电场方向相反。论文图里那些颜色分明的多芯光斑图就是在后处理中选择了“电场模normE”并隐藏了空气孔轮廓得到的。3.3 耦合长度与分束行为的验证手段分束器参数里最重要的一项是耦合长度L_c。对两个相同波导L_c的定义是当入射光从一根芯子输入经过长度L_c后功率完全转移到另一根芯子的长度。对三芯系统耦合长度与超模间的传播常数差直接相关L_c π / |Δβ|其中Δβ是两个超模传播常数的差值。在COMSOL中模式分析得到的是有效折射率n_eff而传播常数β k₀ × n_eff因此Δβ k₀ × Δn_effk₀是真空中的波数。所以复现的关键步骤是算出对称模和反对称模的有效折射率做差值然后除以波长相关的k₀得到耦合长度。把这个长度与论文中的设计长度对比就能验证你的结构参数是否拿对了。我复现时发现三芯系统的耦合长度对芯间距极其敏感。芯间距从2个孔距增加到3个孔距耦合长度可能增大到原来的5倍以上。论文里如果只给了截面结构图没有明确标注对称轴你就得自己试着调整间距直到计算出的耦合长度和论文给出的传输图特征周期一致。这个反推过程比较耗时但一旦数值吻合后面一切就都顺了。另外要提醒一句二维模式分析得到的是横截面的本征解不等于沿长度方向的分束行为。如果论文里有光在光纤中传播的纵向强度演化图那你还需要借助耦合模理论或者建一个三维模型把传输问题显式求解出来。不过大多数论文复现需求停留在截面模式场图层面先算二维超模就够了。4. 实操踩坑记录复现中容易翻车的七个细节4.1 模式搜索范围不当导致老是找错模式模式分析里最折磨人的问题就是求解结果不是你看论文时想要的模式。要么漏了基模要么一大堆模式都是包层模真正关心的SPR模式和超模却找不到。这个问题我几乎每次复现都会遇到核心原因是“搜索时使用的模式数”和“搜索基准折射率”没设置好。我在SPR传感器复现中把模式数从默认的6提高到10并且把基准折射率设置为1.44附近石英折射率附近才稳定锁定了表面等离子体模式。三芯分束器中基准折射率设置为1.44-1.45之间则能优先找到芯层限制模式。还有一个容易忽视的选项“模式分析”默认寻找最接近基准值的有效折射率模式。如果你把基准值设成1.0求解器会优先给你一堆包层模甚至辐射模所以你一定要计算一下包层有效折射率的大致范围。可以直接先用一个简化的等效介质模型估算也可以先跑一次粗网格的模式分析观察结果集中出现在哪个区间再细化搜索。4.2 网格粗导致SPR共振峰偏移网格密度对SPR传感器的损耗峰位置影响非常显著。SPR传感器中金属膜厚度只有几十纳米表面等离子体波的能量集中在金属-分析物界面附近几纳米到几十纳米的范围内普通网格无法分辨这么薄的区域共振位置会偏向不真实的波长。我的建议是金属膜径向至少划分3层网格每层厚度不超过金属膜厚度的1/3分析物通道与金属界面附近的网格最大尺寸不要超过10nm。如果一个COMSOL模型计算量太大可以先用粗网格试跑整个波长范围找到大致共振波长再在共振波长附近做一次加密网格的精细计算。这种“两步法”在复现论文时非常实用。4.3 导入STEP几何文件时警告不断很多人喜欢在SolidWorks里建好几何另存为STEP文件后导入COMSOL这是效率最高的方式之一但往往会遇到一堆警告比如“实体退化”、“面方向不一致”、“边重叠”等。这些警告在大多数情况下不影响最终仿真但如果出现“几何对象未完全闭合”之类的错误就可能导致后续物理场无法选择域。遇到这种情况我有几条处理经验导入后先运行“几何-修复”功能再尝试用“并集”操作把所有实体合并如果修复后仍有问题检查CAD模型里是否有极小圆角或倒角特征这些特征在导入后容易造成短边和退化对象移除后通常能正常使用。另外文件的单位一定要在导入时确认SolidWorks里如果用的是毫米COMSOL默认是米选错单位会让几何尺寸放大一千倍那种情况下所有网格都会失败。4.4 后处理出图与论文配图对齐论文复现的“最后一公里”是把结果图像做得和原图一致。这里不只是审美问题更重要的是信息传达一致坐标轴范围、颜色映射、物理量单位、光斑的强度分布比例都会影响读者对结果是否吻合的判断。COMSOL的后处理中我一般使用“二维绘图组”里的“表面”图绘制“电场模”或者“有效折射率相关的物理量”。首先要把颜色表设置成论文常用的“彩虹”或“灰度”方案然后固定颜色范围。这里有一个很重要的细节COMSOL默认的颜色范围是自动适配当前帧的最大最小值这会导致不同波长下的图颜色标准不一致和论文图对比时误以为强度分布变了。所以一定要在“颜色图例”里取消“自动”范围手动设定同样的上下限。其次是视角问题。截面模式场图通常用正视视角把绘图组中“视图”的X和Y方向固定好取消“自动缩放”让图形居中且比例不变。还可以打开“网格/坐标轴”设置把空气孔的轮廓用细线叠加上去便于和论文结构图相对照。4.5 内存不足与求解时间过长COMSOL仿真对内存的消耗主要在网格求解阶段二维问题还好但如果做了参数化扫描每个波长点都重新剖分和求解时间成本和内存占用都会上来。我复现SPR传感器时波长扫描点几十个每个点都需要求解模式数总耗时一个多小时。优化方法主要有四种一是在保证精度的前提下把PML用映射网格划分每层径向网格数不要太大二是采用“辅助扫描”先把一个波长下算好的模式结果作为下一个波长的初始猜测加速收敛三是关闭多余的后处理默认计算比如不计算所有域中的损耗密度只保留模式分析结果四是在多核机器上开启“分布式求解”COMSOL在集群上可以按波长点分配计算任务省时间很可观。4.6 版本差异导致结果与论文不一致最后说一个很多人没想到的问题COMSOL版本不同材料库插值数据甚至求解器数值精度都有差异。老版本里默认的电磁波频域接口在金属域中可能会自动把“完美电导体”或“阻抗边界”作为默认选项而新版本对金属域的处理更准确。使用材料库时不同版本给出的金介电常数在某一波段也可能差几个百分点这足以影响SPR峰位的精度。复现时如果发现损耗峰与论文偏差较大不要执着地调结构参数先检查材料库版本和介电常数来源。可以自己输入论文里明确给出的折射率或介电常数数据避免材料库版本差异干扰。如果你使用的是老版本COMSOL建议在“模型开发器”里手动检查每个域的物理设置看是否有默认边界条件被错误添加到内部边界上。结语一次复现成功后的几点实在体会跑完这两个复现案例我最深的感受是论文复现最大的价值不是“把图仿出来”而是逼着你把一连串原本模糊的物理概念变成一个个具体参数和操作。拿SPR传感器来说只有亲手计算过损耗谱你才会真正理解为什么金属膜厚度、分析物折射率、环境温度都是传感性能的决定因素。拿三芯分束器来说只有亲手搜出那几个超模你才会明白多芯系统里的“本征模式”和直觉里的“单个纤芯模式”完全是两码事。最后给准备动手复现的朋友一点建议开工之前把所有论文里出现的参数列成一张表格标注清楚是“原文明确给出”还是“需要反推”再决定每次修改哪个变量。同时养成保存多个版本模型的习惯每改变一个关键参数就存一个.mph文件避免调参数调乱了想回退却找不到原始状态。过程会很磨人但当你跑出来的曲线最终和原论文叠在一起几乎重合时那种感觉确实挺值。