做光纤仿真的人都知道单看一个波长下的模式分布远不能判断这根光纤能不能用。我上次调一个色散补偿光纤模型折射率剖面只动了百分之零点几1550nm附近的色散值却变了近一倍。那种情况下脑子里如果没有一张有效折射率随波长变化的完整曲线根本连问题出在哪都说不清。所以“光纤模式色散图”这个题目看起来基础实际上是整个光纤设计的标尺——从模式的有效折射率到群速度、群速度色散、材料色散与波导色散的叠加再到器件的工作带宽全都从这张图里读出来。这篇文章把这些环节完整过一遍从物理原理、COMSOL建模、模式分析、参数扫描到出图解读顺手把容易踩的坑也交代清楚。适合刚入门光学仿真、准备用有限元方法算光纤模式的同学参考也适合已经在算但被模式跳变和伪模式折腾过的老手。我要强调一点色散图不是画出来好看的它是设计判断的依据。同样的光纤参数在1200nm和1650nm下表现可能完全不同而这些差异全都能从模式色散图上提前看出来。1. 先搞清楚色散图到底在算什么东西1.1 模式是什么为什么要关心模式光纤里的光并不是像小时候玩手电筒那样直直地一束光往前走。它实际上是以一种特定的电磁场分布形式在光纤里传播这种分布由光纤的折射率剖面和边界条件共同决定。一组满足麦克斯韦方程组、在传播过程中形态保持不变的电磁场解就叫做一个模式。不同模式的区别在于横向场分布不同纵向传播速度也不同。从数学上看求解光纤模式就是在解一个本征值问题。假设光纤沿z轴方向结构不变电场可以写成E(x, y, z) E(x, y) · exp(-jβz)其中β是传播常数它决定了光在这个模式下的纵向传播行为。常用有效折射率来替代β定义为neff β / k0k0是真空中的波数。neff直观地告诉我们这个模式在光纤里感受到的“平均折射率”是多少。典型的单模光纤中芯层折射率在1.450左右包层在1.444左右基模的有效折射率大概落在1.444到1.450之间。为什么要关心模式因为通信、传感、激光器领域几乎所有的性能参数——模式色散、弯曲损耗、非线性系数、模场直径、截止波长——都建立在模式这个概念之上。多模光纤里不同模式传播速度不同会带来模式色散单模光纤里虽然只有一个模式在传播但它的neff随波长变化一样会产生色散。想要定量地理解这些现象就必须先把模式随波长的变化规律算清楚。1.2 色散图的三张面孔“色散图”这个名字在不同场合下有不同所指最常看到的是这三种第一种neff随波长λ变化的曲线。这是COMSOL模式分析最直接能算出来的东西横坐标是波长纵坐标是有效折射率。对单模光纤来说就是一条单调下降的曲线对多模光纤来说是一簇曲线每个模式一条。第二种传播常数β随角频率ω变化的曲线。这个在物理上更本质因为曲线的斜率直接对应群速度的倒数dω/dβ而曲线的二阶导数对应群速度色散。不过工程上直接看β-ω图的人少多数人只是把它当作中间结果。第三种色散参数D随波长变化的曲线。这是工程设计最常用的形式D的单位是ps/(nm·km)直接告诉你在1nm的波长间隔和1km的长度下脉冲会拉伸多少皮秒。需要说明的是D本质上就是neff对波长的二阶导数换算出来的并不是一个新的独立量。这三种图的换算关系很简单β ω·neff / c。麻烦不在于换算而在于数值精度。neff的绝对数值通常只在小数点后第四位上有变化二阶求导之后这几位有效数字就是全部家底。很多人的仿真结果一求导就全是噪声就是因为最初算neff的精度不够或者扫描波长的步长选得不合适。这一点后面还会专门说。2. 建模型前的物理准备与边界条件设计2.1 从麦克斯韦方程组到二维特征值问题有人第一次用COMSOL算光纤模式时会困惑为什么不用建一段几十微米长的光纤实体然后让光在里面跑其实完全不需要。光纤沿传播方向是均匀的既然场沿z方向只有相位变化、没有形态变化那我们只需要在光纤的横截面上求解本征值问题就够了。这是一个非常漂亮的降维处理从三维问题变成二维问题计算量大幅下降物理图像也更清晰。具体做法是把电场形式代入频域波动方程得到横截面上的亥姆霍兹方程。求解这个方程时传播常数β是未知数被当成特征值来寻找。在COMSOL的波动光学模块里这对应“电磁波频域”接口下的“模式分析”研究步骤。软件会在横截面上求解找出满足边界条件的若干个β再转换为我们熟悉的neff。要注意的是模式分析求出的neff通常是一个复数。实部对应模式的有效折射率虚部对应损耗或增益。如果算的是理想的、无损的光纤虚部基本为零但如果在材料里加了吸收项或用了PML边界虚部就有实际意义了。后面我们判断伪模式时虚部的大小也是一个很管用的筛子。2.2 几何、材料与边界条件的取舍建二维几何时我们要画的是光纤的横截面。以一个标准的阶跃折射率单模光纤为例我通常用这些参数芯层半径4.15μm包层半径62.5μm包层外面再画一层PML区域厚度取8~10μm左右。芯层和包层的折射率用Sellmeier公式定义。Sellmeier公式把折射率表达为波长的函数里面有几组系数不同材料的系数不同。如果你在COMSOL材料库里找不到光纤玻璃的对应项也可以手动输入SiO2或掺杂SiO2的Sellmeier系数这对材料色散的准确计算至关重要。边界条件是模式分析里最关键的一环。很多人直接在最外面用默认的边界条件结果算出来的特征值又密又杂乱。我推荐的做法是在包层外画一个圆环作为PML区域。PML的作用是吸收向外辐射的场让真实的导模不受外部反射的干扰。只关心导模时使用PML比简单地用完美磁导体更稳妥因为实际的求解域边界总会反射掉一部分泄漏的场把结果搞脏。还有一个小技巧如果光纤结构本身是圆对称的有人会考虑用四分之一截面加对称边界来减少计算量。但我要提醒一句在模式分析中不建议这么做。圆对称光纤里的很多模式是简并的例如HE21的两个正交偏振它们的场分布呈旋转对称而不是镜像对称。如果强行用对称边界会把真正的简并模式人为地袭掉一批后面模式识别的麻烦会很大。老老实实算全截面反而最省心。2.3 网格怎么剖才能保住近简并模式网格是模式分析里最容易被低估的环节。弱导光纤中芯层和包层折射率差往往只有0.3%~0.5%不同模式的neff差异可以小到10^-4甚至在10^-5量级。网格如果太粗有限元方法的数值误差完全可能掩盖掉真实模式之间的区别表现为该分开的模式没分开或者出现了很多不存在的假模式。我的经验是芯层区域沿径向至少保证每个波长10~20个单元的密度具体来说在1550nm附近芯层内的最大网格尺寸不要超过0.3μm。芯层和包层交界处必须加密因为折射率突变的位置是所有模式场分布的敏感区域在这里我会把最大单元尺寸压到0.1μm左右。包层可以用较粗的网格但PML内部要能保证至少两个单元层让吸收效果正常工作。网格划分完别急着求解。先看一下单元质量分布确保最差单元质量不低于0.3否则求解时局部误差会非常大。我见过一个案例某个模式在纤芯边缘出现了明显的“锯齿”场分布查到最后就是网格拉伸过度在那个区域单元质量只有0.1。把网格修正之后模式场形变得光滑特征值也立刻合理了。3. COMSOL实操从模式分析到多波长扫描3.1 模式分析研究步骤的配置细节现在进入实操环节。我以COMSOL的波动光学模块为例步骤是通用的不同版本界面上可能会略有差别。首先新建模型选择二维空间维度添加“电磁波频域”物理场接口。研究类型直接选“模式分析”它会自动创建一个带模式分析步骤的研究节点。模式分析的关键设置有两处。第一处是“模式数”一般在初期不确定光纤里有多少模式时我会先填8到10个再根据求解结果收敛。第二处是“搜索基准值”这个值表示你希望求得传播常数附近的模式。对石英光纤基准值设在1.45附近比较合适。基准值设得太偏有可能把模式漏掉设得太靠近包层折射率则会搜出一堆辐射模式来。求解器在后台做的事是寻找满足特征值方程的β并把结果转换成neff。默认情况下COMSOL给出的neff是复数且特征值列表一般按虚部大小排列。这意味着如果一次找到了10个模式列表里最靠前的可能是那些损耗很大的伪模式真正的导模反而排在后面。所以求解完成后不要只看前几个结果要按neff的实部在1.44~1.45范围内筛选。3.2 参数扫描与数据收集的完整流程单波长下的模式分析只是起点。要得到色散图我们需要对整个波长范围做扫描。我的做法是先在模型里定义一个参数lambda_w初始值1.55μm然后把材料定义里所有跟波长有关的表达式都改成这个参数。Sellmeier公式里的波长项也要引用它这一步如果漏了扫出来的曲线就是一条没有物理意义的假直线。扫描范围通常结合应用场景来定。单模通信波段我会覆盖1260nm到1675nm如果是在设计少模光纤或特种光纤范围可能从800nm一路扫到2000nm。步长的选择原则是先粗扫步长取10nm把趋势和模式数量摸清楚精细扫描再取1nm到2nm用来看色散零点的精确位置。如果你直接上1nm步长在全波段扫描很可能会因为某个波长的模式阶次判断出错导致整条曲线在后期处理时很难收拾。参数扫描在COMSOL中可以直接用“辅助扫描”挂到模式分析研究下面。开启辅助扫描后选定lambda_w为扫描参数输入扫描列表。每一次扫描求解都会保留下来并把结果存成“解1”、“解2”等形式。导出数据时我习惯使用“全局计算”节点定义表达式计算各模式的real(neff)然后导出为表格。记住给每个模式也加一个编号变量否则多模式情况下你根本分不清哪条曲线对应哪个模式。3.3 简并模式与模式识别的现场处理等你真的开始扫描会遇到一个让所有新手头疼的问题模式排序在扫描过程中会变。比如在短波长范围某个模式可能是第3个被解出来的到了长波长范围它变成了第5个。这不是软件出错了而是特征值求解器对特征值排序的方式不是按物理模式身份排序而是按数值大小排序。如果你直接按“第几个模式”去画neff随波长的曲线画出来的图就像一堆乱麻。处理这个问题我有两个办法。第一个是在模式分析设置里勾选“在频域中搜索相近特征值”之类的选项让求解器尝试追踪模式。具体名称因版本而异但原理是以上一个解为初值继续搜索这样同一个物理模式在相邻波长点之间的neff变化很小求解器容易延续到同一条分支上。第二个方法更稳妥扫描完成后不要急着画图把所有模式的neff结果导出来按波长升序排列再用脚本或电子表格重新分组。判断依据是同一个物理模式在相邻波长之间的neff差值应该很小一般小于0.001。如果某条曲线的neff在某点突然跳了0.01以上那几乎可以肯定是排序接错了。这时候看场分布是最后的仲裁手段把疑似跳变波长处的几个模式场形画出来逐个对号入座曲线就理清楚了。还有一个细节圆对称光纤里的LP模经常简并例如HE21模的两种正交偏振态在理想圆对称结构下neff完全一样。有限元网格本身不可能是绝对对称的所以你会看到两个neff值非常接近但不完全相等的解这很正常。画图的时候这两条线会几乎重合不要以为是重复模式而删掉其中一条。4. 色散图的绘制、后处理与结果解读4.1 三种出图方式与选择逻辑模式扫描的数据在手接下来就是出图。如果只是看模式数量和大体趋势neff-λ曲线最直观。横轴波长纵轴有效折射率。每条模式一条线大致呈下降趋势。单模光纤只有一条线它会随着波长增加平缓下降多模光纤是一簇线不同模式在短波长方向散开在长波长方向逐渐靠近截止状态。如果想对照理论公式我更喜欢画β-ω图。传播常数β 2π·neff/λ角频率ω 2πc/λ把扫描解中的neff和lambda_w带进去就能画。这条曲线的局部斜率意义非常明确dω/dβ是群速度也就是光脉冲能量在光纤里实际传播的速度。曲线的弯曲程度则直接对应群速度色散。第三种是参考用的D-λ图这是工程上最关心的。后面会专门讲计算方式。我的个人习惯是先用neff-λ图做模式识别和正确性检查等确认曲线干净、模式归属无误之后再做D-λ图用于设计讨论。不建议一上来就直接求二阶导否则噪声会把真实信号全吃掉你只会得到一条乱抖的线。4.2 从有效折射率计算色散参数D色散参数D的标准定义是群时延对波长的导数除以长度换算到neff之后可以写成D - (λ / c) · d²neff / dλ²这个式子看起来简单但实际计算时非常考验数据质量。d²是二阶导数意味着neff数据里哪怕是10^-6量级的随机误差求导之后都会放大成明显的锯齿。所以计算D时主要有三条经验。第一扫描步长不能太稀疏。做色散参数计算时我把精细扫描步长定为1nm甚至0.5nm。数据点越密差分格式引入的截断误差越小。第二不要用COMSOL默认的有限差分方式直接对原始曲线求导除非数据已经被平滑过。我通常先把neff-λ曲线导出在数据统计分析软件里先做平滑例如Savitzky-Golay滤波再求二阶导数。第三算出来的D要跟文献或量级直觉对比一下。普通标准单模光纤在1550nm附近D约17ps/(nm·km)正值且随着波长增大而增大。如果你的计算结果在正负几十之间乱跳先别怀疑公式回去看原始neff曲线是否够光滑。4.3 一张合格的色散图应该长什么样色散图画完怎么判断它算得对不对我的检查清单第一项是neff的落点范围。所有导模的neff应该严格介于包层折射率和芯层折射率之间。低于包层折射率的解基本是辐射模式或伪模式可以忽略。第二项是曲线形态。正常导模的neff随波长单调递减曲线光滑没有拐弯处的异常抖动。如果有模式在某个波长突然消失说明它在这个波长下已经截止。第三项是模式数是否符合预期。一个V参数在2.4以下的阶跃光纤应该只有基模V参数接近5时会出现几个高阶模。如果算出来的模式数远多于理论预测很可能是伪模式混进来了。还有一点我想单独强调看多模色散图时要学会识别截止波长。某个模式的neff曲线在长波长方向逐渐靠近包层折射率这个交点就对应截止波长。设计多模光纤或少模光纤时模式截止的位置决定了一大堆关键指标所以这张图上“曲线在哪里终结”与“曲线本身长什么样”同样重要。5. 常见问题与排查技巧实录5.1 模式少了或排序混乱怎么办模式数不够是最常遇到的第一类问题。明明理论上有5个导模结果模式分析只给出3个。这通常有几种原因搜索基准值离真实值太远模式数设置太少或者网格太粗导致近简并模式被数值误差吞掉了。我的排查顺序是首先把模式数加到10以上基准值分别试1.44、1.45、1.46几个点看看新增解是否稳定出现。其次检查网格在芯层做一次细化如果细化后模式数变多了就说明之前确实是网格背锅。另外在扫描长波长范围时很多高阶模式已经截止这时候求解出“模式少”是正常的不要硬凑。排序混乱的问题我的建议是始终以“neff值连续性”为分组依据。任何时候都不要相信模式编号的绝对顺序它在不同波长之间可能完全是乱的。写完扫描数据之后写一个小脚本按neff值排序分组比人工逐条核对高效得多。5.2 伪模式满天飞怎么过滤伪模式几乎每个人都会碰到。它们的特点很鲜明neff要么低于包层折射率要么虚部大得离谱场分布看起来在边界上乱跳或者主要能量滞留在PML内部。伪模式的来源主要是两个。一个是PML吸收效果不理想辐射模式在PML内部反射回来形成不真实的驻波解。这种情况把PML加厚并增加PML内部的网格密度效果立竿见影。另一个是求解域边界太靠近真实模式把本来向外扩散的场硬憋成了一个高损耗的假导模。解决办法是增大包层区域或者增强PML的衰减效果。在过滤伪模式时一个很实用的标准是看neff的虚部。真正的导模在无损材料里虚部应该非常小通常小于实部的千分之一。伪模式的虚部则往往跟实部同一数量级或更大。扫描过程中就可以设置筛选条件只保留虚部小于某个阈值的解能省去很多后续的清理工作。5.3 扫描结果曲线“毛刺”的根因与修复画出来的色散曲线如果出现毛刺不要急着觉得是物理上正常。毛刺通常只有两类原因。第一类是前面说过的模式排序跳变某条曲线在一个波长点后突然接到另一个模式上去了。修复办法就是把数据重新分组或者勾选模式追踪选项重新求解。第二类是网格误差导致的数值抖动特征是在很短的波长区间内neff有小幅波动波动幅度在10^-5到10^-4量级。网格误差型毛刺有一个隐蔽但常见的来源波长变化时如果网格划分是自动生成的每个波长下网格可能都会变导致相邻波长点的有限元误差不一样相当于给曲线叠加了一个随机噪声。解决办法是在扫描前固定网格。先选定一个代表波长比如工作波长划分好网格然后让物理场设置里的默认网格“不随波长变化”。这个操作让neff曲线的光滑程度提升非常明显再求D时你会直接感觉到好处。如果做完以上所有清理曲线还是有不规则的凹凸那就要回头看看材料折射率的设置。Sellmeier公式如果只在部分波段有效或者折射率表达式中用了不连续的插值表也会造成曲线出现非物理的弯折。这类问题在长波段尤其容易出现检查方法很简单把材料折射率随波长的单点曲线单独画出来如果材料曲线本身就不光滑那么模式色散图的毛刺就从这里来。最后再分享一个实操中的小技巧每次跑完模式扫描我会顺手把neff数据导出一份原始表格存下来。这个习惯救过我很多次——有时候后处理阶段发现曲线不对不用重跑十几个小时的扫描直接拿原始数据重新分组就能定位问题。色散图这个东西前期数据质量决定了后面所有分析的可靠性。你愿意多花一点时间在网格和扫描设置上后面画图和解图就会舒服很多。
