基于区间遍历法的纤芯及包层模有效折射率的计算.docxVIP

基于区间遍历法的纤芯及包层模有效折射率的计算 1 包层模有效转化率计算 目前，光绪以其良好的传输特性已广泛应用于光纤通信、光纤传感器、光纤存储等领域。需要精确地计算纤芯模型和包层模型的有效折射，以便计算间隔色散模型和不同模型之间的耦合。为了准确计算纤芯模型和包层模型的有效折射，需要更高的计算精度、更快的计算速度和更精确的计算方法。这在计算光纤纤芯模型和包层模型的有效折射方面发挥着非常重要的作用。 目前,对于光纤纤芯及包层模有效折射率的计算及仿真研究的文献主要分为两大类,第一类是通过求解在弱导近似条件下得到的光纤本征方程来获得有效折射率值的方法,这由于本征方程简单而可以化简很多计算过程,但其计算结果精度较低.第二类,使用精确特征方程计算纤芯及包层模的有效折射率,通常采用二分法求解窄线宽光源下纤芯模的有效折射率,但当光源线宽较大时,二分法的收敛速度很难满足计算要求;由于包层模的本征方程是带有奇点的超越方程(类tan函数),二分法并不适用于此类问题的求解,并且至今未见文章详细阐述求解多值超越方程得到包层模有效折射率的算法. 本文利用三层结构光纤模型,对光纤纤芯模及包层模有效折射率进行了计算及仿真研究;应用截弦法快速精确求解了纤芯单值超越本征方程,得到了纤芯模的有效折射率;并通过COMSOL模场仿真软件验证了仿真结果的正确性;通过使用区间遍历及截弦法求解包层本征方程这一多值超越方程,得到了各次包层模的有效折射率;采用Mathematica软件对求解过程进行仿真研究,绘制了纤芯模、包层模有效折射率与波长关系曲线.本文详细论述了具体求解过程及方法,其求解方法适用于光源带宽较宽、要求精度较高情况下光纤中纤芯及包层模有效折射率的求解. 2 理论分析与模拟 2.1 模拟结果与分析 光纤三层结构模型如图1所示,将光纤看作纤芯、包层和环境层的三层结构.最内层为纤芯,半径为r1,折射率为n1;其外层为包层,半径为r2,折射率为n2;最外面环境层的折射率为n3,半径可视为无穷大.单模光纤纤芯模的本征方程为 式中,归一化横向相位参数U=ur1,u2=k02(n12-n2effco);归一化横向衰减参数W=wr1,w2=k02(n2effco-n22);Jl和Kl分别为第l阶第一类贝塞尔函数和第l阶第二类虚宗量贝塞尔函数;Jl′(U)和Kl′(W)是分别对U,W的一阶微分;βco为纤芯模的纵向传播常数,且βco=k0neffco,neffco为所要求解的纤芯模的有效折射率,k0=2π/λ为光在真空中的波数,λ为光在真空中的波长. 对(1)式中Jl′(U),Kl′(W)做如下变换,有 令 选用参数为n1=1.4681,n2=1.4628,r1=8.3μm,r2=62.5μm的Coning SMF-28光纤为模型,使用Mathematica绘制的λ为1550 nm时neffco与f(neffco)的关系曲线如图2,从中可以看出,曲线与横轴有且只有一个交点,曲线与横轴的交点为neffco的值,该图中neffco的取值为1.46698. 在单模条件下,函数f(neffco)=0为单值超越方程,无法求得准确的解析解,使用f(neffco)=0的近似数值解来得到纤芯模的有效折射率neffco.每一个输入波长λ都有惟一的neffco与之相对应,但如果研究涉及的光源带宽较宽、计算精度要求较高,那么计算量就会变得很大,该条件下二分法难以满足计算需求.本文使用截弦法来求解超越方程(4)的根,该方法收敛速度快,计算量相对较小.首先取[na,nb]?[n2,n1],且na,nb无限接近n2,n1,计算(na,f(na))与(nb,f(nb))连成的直线与横轴的交点记为nc,若f(na)f(nc)0,则取(na,f(na))与(nc,f(nc))与横轴的交点nd,否则取(nb,f(nb))与(nc,f(nc))与横轴的交点为nd.重复上述计算,直到|nc-nd|ε,ε为计算精度(本文精度为10-4). 利用COMSOL场分析软件对计算结果进行验证,当选取的建模参数与上述模拟相同时,也就是当设定入射波长为1550 nm时,设置的纤芯的有效折射率为1.46698,获得的光纤中模场分布图如图3(a)所示.从图中可以看出,在纤芯的有效折射率为1.46698时,光场模被限制在纤芯中,也就是实现了该光场模式的有效传输.而设置纤芯的有效折射率为另外值时(如1.465333),光场的分布如图3(b)所示,这时光场的传输并不集中在光纤纤芯部分,该有效折射率不是纤芯模式对应下的有效折射率,这进一步验证了我们上述模拟求解纤芯有效折射率方法的准确性. 利用此方法,使用Mathematica软件计算了波长λ在1400—1600 nm,波长取样间隔为2 nm时所对应的neffco.光纤参数为n1=1.468