标量FDTD法分析渐变折射率光波导模场分布.docx

标量FDTD法分析渐变折射率光波导模场分布.docx

《标量FDTD法分析渐变折射率光波导模场分布.docx》由会员分享，可在线阅读，更多相关《标量FDTD法分析渐变折射率光波导模场分布.docx（11页珍藏版）》请在冰点文库上搜索。

标量FDTD法分析渐变折射率光波导模场分布

收稿日期:

2001208215.

作者简介:

刘　靖（19692,男,硕士研究生;武汉,华中科技大学光电子工程系（430074.基金项目:

教育部高等学校骨干教师资助计划（GG214021051121007.

标量FDTD法分析渐变折射率光波导模场分布

刘　元秀华　华中科技大学

光电子工程系

黄重庆

岳阳师范学院物理系

摘要:

用标量时域有限差分法（Finite2differeceTime2domainMethod,FDTD法分析了渐变折射率平面光波导中基模的模场分布情况.采用平面光波、球面光波、高斯光波等5种不同形态的光波激励同一波导,借助计算机进行数值求解,得到的光场分布图在光场传播达到稳定后完全相同,说明波导的模式与光波形态无关,只与波导结构和光波波长有关.结果表明该方法直观、精确、快速,并与解析法分析的结果一致.关　键　词:

渐变折射率光波导;基模;模场分布;时域有限差分法

中图分类号:

TN252　　文献标识码:

A　　文章编号:

167124512（20020220063203

　　随着集成光学和光通信的发展,光波导的应用范围越来越广,

精确地分析和模拟光波导器件的性能对于发展光电子学是十分重要的.对光波导性能的分析方法虽多,但一般都有各种各样的局限性.随着计算机技术的发展和相应CAD软件的研制,时域有限差分（FDTD法[1]注目,它具有直观、快捷、和计算时间、系列优点,对各种物体的作用,而用于光波导研究的报导还很少见.本课题将这种方法应用于光波导模场的研究,分析了渐变折射率平面光波导中基模的模场分布情况.

1　基本原理

渐变折射率光波导薄膜层折射率为抛物线型

分布,波导截面上的折射率变化可用图1表示.中线x=x1处的折射率最大,其值为n1;x0为曲线与x轴的交点到x1之间的距离;覆盖层及衬底

图1　渐变折射率波导横截面结构

的折射率为n2,且是均匀的.折射率分布具体可

写为:

n2（x=n21[1-（x-x12/x20]

（|x-x1|≤h/2;

|x-x1|>h/2.

（1

h为薄膜层的厚度;x方向为折射率变化方向.在|x-x1|比较小时,式（1可写成如下平方律近似公式:

n（x=

n1[1-（x-x12/（2x20]

　　　　　（|x-x1|≤h/2;

n2　　　　（|x-x1|>h/2.

光波的传播方向为z方向,波导结构与y无关.由麦克斯韦方程可知,平面光波导TE模的场分量只有Ey,Hz,Hx,其中Hx和Hz均可用Ey表示,所以电场强度Ey的分布情况就是总场的分布情况,对于电磁场矢量的每个坐标分量均满足标量波动方程.如果从标量波动方程入手,用二维标量FDTD法求解更为简单.Ey满足如下标量波动方程:

52Ey/5x2+52Ey/5z2-μ

ε（52Ey/5t2=0.（2

　　根据FDTD法原理[2],用符号EN（i,k=E（iΔx,kΔz,NΔt代表场分量E（x,z,t,其中Δx和Δz分别为沿x方向和z方向的空间网

格步长,Δt为时间步长,i和k为空间步长个数;

N为网格步长个数.对Ey采用中心差分近似并

取正方形网格,即令Δx=Δz=Δs,则式（2变

第30卷第2期　　　　华　中　科　技　大　学　学　报（自然科学版　　　　　Vol.30　No.22002年　2月　　　J.HuazhongUniv.ofSci.&Tech.（NatureScienceEdition　　　　Feb.　2002

为差分方程

　EN+1

y（i,k=p（i[ENy（i-1,k+

ENy（i+1,k+ENy（i,k-1+EN

y（i,k+1-4ENy（i,k]+2ENy（i,k-EN-1

y

（i,k,式中p（i=Δt2/（μ0ε0Δs2n2i.

由式（1得:

n2

i=

n2

1

1-4

（i-iΔs

h

2

1-

n2

n21

　　　　（|（i-i1Δs|≤h/2;

n2

2　　　（|（i-i1Δs|>h/2,

式中i1=X1/Δs.

为保证FDTD法的数值稳定性[3],时间步长与空间步长的选取应满足

Δt≤[1/（Δx2+1/（Δz2]1/2/vmax.　　若取Δs=2vmaxΔt=0.05λmin,为了消除数值色散[4],要求空间步长Δs与波长λ必须满足关

系:

　p（i=0.25

n22

n21

1-

（Δh

2

1-

2

2

1

　　　　　　（|i-i1（20.25　　　　|1|（2

Δs.,网

格在x和z方向必须加以限制,这可通过把吸收边界条件[5]应用于网格边界来实现.例如若x=Δx为下边界（i=1,其二阶近似吸收边界条件为

　[52/（5x5t-（1/v（52/5t2+

（v/2（52/5z2]Ey=0.

其差分格式为

　EN+1y（i,k=-EN+1

y

（2,k/3+7[ENy（1,k+ENy（2,k]/6+[EN

y（1,k-1+ENy（i,k+1+ENy（2,k-1+EN

y（2,k+1]/12-EN-1y（1,k/3-EN-1y

（2,k.其他各边界也有类似的差分公式.

2　数值计算结果与讨论

本文所分析的渐变折射率光波导参数为:

薄

膜层的厚度h=1μm,λ0=1.3μm

波导中线处折射率为n1=

1.563,覆盖层和衬底的折射率均为n

2=1.

550,Δs=0.04159μs.最后由计算机模拟,得到5种不同光源激励下的渐变折射率光波导基模光场分布如图2~6所示.由于平面光波、

球面光波、高斯光波、sin（axsin（bt光波在i方

向上的场源分布情况与基模的稳态模场分布比较接近,因此光场只需传播较短的距离和时间便达到了稳定.而cos（axsin（bt光波与基模的稳态模场分布趋势刚好相反,光场需要传播较远的距离和较长的时间才能趋于稳定;图中可以看出,用

图2　平面光波导中基模的光场分布

图3　球面光波导中基模的光场分布

图4　高斯光波导中基模的光场分布

图5　sin（axsin（bt光波导中基模的光场分布

4

6　　　　　　　　　　　　　　华　中　科　技　大　学　学　报（自然科学版　　　　　　　　　　　第30卷

图6　cos（axsin（bt光波导中基模的光场分布

cos（axsin（bt光波作为激励时,光场分布的起

始部分（k=1~100与基模的分布情形相差较远,还没有形成模式;随着传播距离的增加,在光场分布的较后面部分（k=110~250,光场传播渐渐达到稳定而形成基模模式（图6所示.以上结果说明:

场源形态与稳态场分布越接近,则光场达到稳定所需要的时间越短,反之需要较长时间才能形成模式.对同一波导,无论用哪种光源激励,达到稳定后i方向的二维稳态模场分布都是相同的（图7所示,的结构和波长有关,

图7　渐变折射率光波导基模的二维模场分布

形显示技术简单易行,显示结果直观、快捷,数据

精确、可信度高等特点,而且与其他方法所得结果一致,通过设置5种不同的激励对基模进行分析,借助计算机进行数值求解,发现场源的形态和稳态模场越接近,则达到稳定所需要的时间和传播的距离越短.但对于同一波导,无论用哪种激励,经过一定的时间后都可以形成相同的稳定模式,且与解析方法讨论的结果一致[6~8],说明波导的模式与激励源的形态无关,只由波导本身的结构参数和场源的波长决定.

参

考

文

献

[1]YeeKS.Numericalsolutionofinitialboundaryvalue

problemsinvolvingMaxwell’sequationsinisotripicmedia.

IEEETrans.

AntennasPropagat.,1966,

14（1:

302~307

[2]高本庆.时域有限差分法（FD2TDMethod.北京:

国

防工业出版社[A,solutionofsteady2

scatteringproblemsusingthetime2Maxwell’sequations.

IEEETrans.on

MicrowaveTheoryTech.,1975,23:

623~630[4]王长清,祝西里.电磁场计算中的时域有限差分法.北

京:

北京大学出版社,1994.

[5]葛俊祥,李玉星.广义正交坐标系下FD2TD算法的吸

收边界条件.电子学报,1998,26（3:

115~116

[6]徐保强,杨秀峰,夏秀兰.光纤通信及网络技术.北京:

北京航空航天大学出版社,1999.

[7]明　海,张国平,谢建平.光电子技术.合肥:

中国科学

技术大学出版社,1998.

[8]彭江得.光电子技术基础.北京:

清华大学出版社,

1999.

Theanalysisofmodefielddistributionsofthegraded

indexopticalwaveguidebythescalarFDTD

LiuJing　YuanXiuhua　HuangDexiu　HuangChongqing

Abstract:

ThispaperanalyzesthemodefielddistributionsTE0modeofthegradedindexplanarwaveguidesbythescalarfinite2differencetime2domain（FDTD.Fivekindsofopticalsourcesareusedasstimuliwhenanalyzingthesamewaveguide,namelyplaneopticalwave,Guassopticalwaveandsphericalwaveetc.Thenumericalsolutionsareresolvedbycomputer,thefinalfiguresofmodefielddistributionsarethesameafterthesteadilyspreadingopticalwavesthatindicatethatthedistributionsofthestableopticalfieldshavenothingtodowiththestimulatedopticalwaves.Thedistributionsarejustconnectedwiththestructuralparametersoftheopticalwaveguideandtheopticalwavelength.Theresultisdirect2viewing,precise,rapidandidenticaltothediscussionbytheanalyticmethod.

Keywords:

gradedindexopticalwaveguide;TE0mode;modefielddistributions;FDTDmethodLiuJing　Postgraduate;Dept.ofOptoelectronicEng.,HuazhongUniversityofScienceandTechnology,

Wuhan430074,China.

5

6第2期　　　　　　　　　刘　靖等:

标量FDTD法分析渐变折射率光波导模场分布