第二章光束传播法基本原理.docx

1、第二章光束传播法基本原理第四章光束传播法基础第一节数值计算方法1.电磁场数值计算它是一种基于麦克斯韦方程组，建立逼近实际工程电磁场问题的连续型的 数学模型,（合理的假设）然后采用相应的数值计算方法，经离散化处理， （合适的方法，使离散化的模型既能反映连续型模型的特性，又便于计算机分析 ）把连续型数学模型转化为等价的离散型数学模型，计算出待求离散数学模型的离散 解（数值解），从而获得相应结果的一种方法。2数值方法分类：时域分析、频域分析。时域分析：模拟光在波导中的传播过程频域分析：求解波导模式时域分析 逼真：把原来因为速度太快、结构太小、不可见的现象模拟出来， 能够直观地展示。求解：波导连接、耦

2、合、非线性特性、波导模式。频域分析：光场分布、给定具体结构波导的模式的有效折射率 （色散、偏振）、损 耗（材料吸收、结构本身导致）等。问题：频域结果能否推得时域信息？反之？3. 常用数值方法简介（1）有限差分法（频域有限差分法）（20世纪50年代出现）利用划分网格的方法将定解区域离散化为网格离散 节点的集合，然后基于差分原理，以各离散点上函数的差商来近似替代该点上 的偏导数，这样待求的偏微分方程定解问题可转化为一组相应的差分方程的问 题。根据差分方程组，解出各离散点上的待求函数值，即为所求定解问题的离 散解，再应用插值方法便可从离散解得到定解问题在整个场域上的近似解。原理:偏导差分方法特点：原

3、理简单、通用性好；对复杂结构，计算量大 （矩阵运算）0 （频域分析）适用范围：计算光波导的模式求解。现状：适用于较简单结构的分析。但有限差分（偏导 差分）法广泛应用于数值方 法中（2）有限元法20世纪40年代提出，其在电磁问题方面的应用有约 40多年历史。以变分原理为基础，把所要求解的微分方程转化为相应的变分问题，即泛函求极值问题。常见方法为把要分析的区域划分为很多三角形 （每个三角形成为一个基元），每个基元内的场用多项式来表达，然后加入不同基元间场的连续条 件，就可得到整个横截面的场分布。特点：较复杂-需要前处理（三角化，剖分）；后处理：（场分布，伪解剔除） （通用性强，精度高）根据该方法对

4、于各种各样的电磁计算问题具有较强 的适应能力性，所形成的代数方程矩阵求解容易、收敛性好。主要缺点：对于形状和分布复杂的三维问题，由于其变量多和剖分要求细，往往因计算 机内存而受到限制。程序设计复杂、计算量较大。适用范围：求解光波导的模式（有效折射率、色散、双折射、传输损耗等） 。现状：功能最强大的数值方法之一。特别是上世纪 90年代出现的矢量有限元方法，完全解决了有限元方法出现的伪解问题，大大降低了有限元法的后处理 过程。有限元光束传播法。（3）时域有限差分法时域有限差分法是近年来开始流行的一种数值模拟方法， 它通过将麦克斯韦方程在时间空间上离散化的方法实现对电磁波传播的模拟 。它能够得到电磁

5、 波传输的瞬态（即时域）信息，通过傅里叶变换即可得到相应的频域信息。时域有限差分法由K.S.Yee于1966年首先提出，此后经过众多学者的努力， 使之不断完善，现已比较成熟。但是，在许多方面它仍在继续发展，解决问题 的能力和应用范围仍在不断地提高和扩大。计算过程为：设置初始场，然后依时间步推进计算，并在每一时间步交替地 计算每一离散点的电场和磁场。特点：不需要矩阵运算，只需简单的加减乘除运算由前一时刻的场来获得下 一时刻场的值。而且，它还非常适合于并行计算，这正好与当今计算机的发展 趋势相吻合，这就更加提高了时域有限差分法解决实际复杂问题的能力。 适用范围：计算光波导的模场分布、有效折射率；研

6、究波导之间的连接、耦合 问题。注：主要用于一维和二维光波导的分析。三维波导分析计算量稍大。 现状：ADI FDTD，可应用于各向异性介质，非线性介质， PML吸收边界光束传播法（Beam Propagation Method，简写 BPM）光束传播法是目前光波导器件研究与设计领域最流行的方法之一， 其基本思想是在给定初始场的前提下，一步一步地计算出各个传播截面上的场 。光束传播法最早是由M.D.Feit等人于1978年研究光场及大气激光束传播时提出的。 最早的BPM是以快速傅里叶变换（Fast-Fourier Transform，称FFT）为数学手段 实现的，称为FFT-BPM。FFT-BPM

7、源于标量波方程，只能得到标量场（即只能 处理一个偏振分量），不能分辨出场的不同偏振（TE模或TM模）以及场之间的耦 合。由于上述缺点，D.Yevick等人于1989年提出了一种新方法一有限差分光束 传播法FD-BPM，用差分的方法将横截面上的场离散化。这种方法已被成功地应用于分析 丫型波导及S型弯曲波导中的光波传输， 且对损耗的计算也得到了准确的结果；FD-BPM还被用于分析条形波导、三维 弯曲波导、二阶非线性效应以及有源器件。频域分析方面，同样可采用光束传播法进行分析：可采用 相关函数法获得， 还发展了一种称为虚轴光束传播法的方法，用于分析波导中的模式。其实，BPM与FDTD有不少相似的地方

8、。其不同在于，FDTD每次都要同 时计算整个波导的模场，而BPM只算一个面 特点：计算量较小，应用范围非常广泛适用范围：计算光波导的模式、色散、双折射、传输损耗等；分析波导传输、 连接、耦合，光栅的传输特性等。4数值方法发展趋势：方法融合现象明显(有限元法与光束传播法的结合形成了另外一种方法 一有限元光束传播法(FE-BPM )。)、相互推动(PML FDTD , BPM, FEM)。第二节有限差分光束传播法基本原理光束传播法(BPM)的基本思想就是把波导沿着传播方向剖分成若干个截面,根据前一个或几个截面上的已知场分布得到下一个截面上的场分布BeaBPROP Sisulation Paraae

9、tersDomain Min;Domain Max:Compute Step:Slice Step:Monitor Step:Current Value 28Default UseValue Defs0.20.2Current Value 24 2iDefault Use Vilue DefsCurrentValue|5-U2QDefault Use Value Defs r 1112010100lio-|100 LEstimated Time:0.088 minSave SettingsOKCancelContour Map of Index Profile at Y=01000 -800

10、-600 -400 -200 -20 -100 10 20X(nm)I ISIbuIat ion P ara*eterE Compute Index ProTileXCu rentValueDefaultValueUse DeFsDomw Mir:17Domain Max：|27.2PCompute Step:02-j 0.2审Slice 5切 Moriit&r Step:Current Default Ute Value Value DeFsCurrent Default Use Value Value Defs| 卫2 |-2：2 170 0 |7f23.2 . 2 破1100|7El-

11、PC2 pn Pioo|100|1QOio r17100 |7S/mbals.Diplav M&deOutput Fie Piefix:Estimated Timer0.006 minIContauiMap)亍Save SettingsDisplay. utpul._5k1! kkb CancelContour Map of Index Profile at Y=01200woo -8006004002000-201020 =LN0X (pm)BPM理论来源于波动方程，波动方程是建立在Maxwell方程的一般形式为E(t)畀1.4545Maxwell方程基础上的。(1a)H(t)J(t)D(t

12、)D(t)B(t)(lb)(lc)(ld)式中，E为电场强度，H为磁场强度，D为电矢量位移，B为磁感应强度，J 为电流密度矢量， 为体电荷密度，t为时间。对于各向同性、非磁性、电中性介质，有J(t)E(t),B(t) oH(t),D(t) E(t)式中，为电导率,0为真空磁化率常数,为介电常数。将式(2)代入式有H(t)门E(t) o0 t(3a)E(t)H(t) E(t)t(3b)考虑到场对时间的依赖E(t)Eexp(i t),H(t)H exp(i t)式中，E ExiEyj Ezk,H Hxi H yj H zk为复振幅，为角频率，i为单位虚数。把式(4)代入式，有E i 0HH ( i

13、0)E(5a)(5b)式(5a)可进一步写为(E) i 0(H) 0(6)将式(5b)代入式，有(E) 2 0(i )E 0定义复相对介电常数式中， 由下式表示量波方程对于任何矢量G，有从而可以进一步得到(15)把式(12)和式(16)代入到式(9)和(11)，可以得到2 12E (% % E)2 1 22H %( H) 2H 0 (18)%考虑准TE模(Ez 0)和准TM( Hz 0)模，有2 1 22E ( % E ) 2E 0 (19a)%2 1 1 % 2H %( H) k (：k H ) H 0(19b)式中i x j y(20)EExiEyj(21)HHxiHyj(22)将式(19

14、)写成分量形式如下 (23a)(23b)% HyHx%Hx2Hx2Hy进一步可以写为1 %Hy% zz2 .2 .笔XEx22yz2Ey0x y2Ey2Ey22xzx0%EyEyHy%Ey2Ex2Ex%Ex(24a)(24b)%y式(24)中，Ex，2Hx2Hx1%Hx2 x2 z%zz2Hx -2Hy% y1Hy x% x化1%Hy2yz2%zzH -2Hx%x1Hx y% yy xy(24c)y，z的函数,H y是空间坐标x，%Ex1 Hy% x1 Hx% y把Ex, Ey, Hx, H y随z的快速周期变化部分分离，令Exx(x,y ,z)exp( iz)(25a)Eyy(x ,y,z)

15、exp( iz)(25b)Hxx(x,y,z)exp( iz)(25c)Hyy(x,y,z)exp( iz)(25d)n。(26)c式中，n为参考折射率，选择时应尽量接近导模的有效折射率，否则会影响计 算精度。 x， y， x， y为包络函数，如图i所示。E E H Hx y x y图2-1包络函数示意图将式(25)代入式(24)可得包络函数的矢量波方程x % x % x2 2x x 2i2 2y z2y2z2i%y%x2x2z进一步整理得2iyux(27b)1%x%zz1%yx0(27c)1%y%zz1%xy0(27d)AxxxAxy yx2i%y2i%xx2x式中xAyyz22i(28b)

16、2x2z2iBxxBxy(28c)Ayx2iByyByx x(28d)1_% x2x2y12x2x1 %i% zBxyByy y %三维半矢量形式Byx忽略x , y场之间的耦合，则有2x2z2iAxx2i2x2z2iBXX X2y2z2iByy y3.三维标量形式（忽略场的方向性）2i2i4.二维半矢量形式与三维半矢量类似，但是折射率分布更简单5 .二维标量形式2.1.2方程离散数值处理(1)纵向数值处理直接求解方程式(27)是非常困难的，因此需要对它进行离 散化处理，通过数值方法来求解。 BPM的数值离散化处理方法很多，这里采用 有限差分方法来实现，有限差分法的核心就是把导数写成差分的形式

17、。为了便于方程的求解，有必要对方程进行近似处理，在纵向 (即沿着光的传 播方向Z)的近似处理有缓变包络近似(SVEA)、广角近似等。在这里，我们采用 缓变包络近似方法。如果包络函数随Z的变化足够缓慢，使得2齐0z则有22i 2i Z Z Z其中， 代表上述各个包络函数。纵向处理主要是解决传播方向上相邻两个截面上场的关系问题，从式 (29)可以看出，纵向处理就转化为对Z的一阶偏微分处理。设相邻的两个截面分别用I和I 1标志，第I个截面上的场为已知，第I 1截面上的场待求，两截面的间距为 Z,见图图2相邻截面差分格式示意图式(29)中右边的偏微分项可差分为(30)式28可表示为(31)类似的方法还

18、有:稳定性：是指计算过程中积累误差是无限增加还是可以控制， 对于沿z方向数值损耗：由数值计算引起的沿传输方向上的能量损失，是非物理损耗。研 究表明， 0.5时，数值损耗最小， 1时，数值损耗最大。因此，综合考虑到数值计算的稳定性及数值损耗， 在计算中，要仔细选择合适的 值，使得稳定性和数值损耗都可以接受。这样，式(28)可进C步化为Axxl 1xC1Axxx Axy y(32a)CAyyl 1 yC1Ayyy Ayx x(32b)DBxxl 1xD1Bxxx Axy y(32c)DByyl 1 yD1Byyl B l y yx x(32d)其中C2iD2i1%zzz%z缓变包络近似 特点:最早

19、提出，方法最简单广角近似 (Pad e近似特点:近似更少，更高精度图3有限差分网格结构在有限差分光束传播法中，横向处理可采用上述的九点差分格式，设要差分的变量为 (x, y)，则在点(m, n)上, 的一阶和二阶导数可差分为 (m 1,n) (m 1, n)X (m, n) 2 X(33)_ (m,n 1) (m,n 1)y(m,n) 2 y(34)2(m 1,n)(m 1,n)2(m, n)22X(m,n)X(35)2(m, n 1)(m, n 1)2(m,n)22y(m,n)y(36)211m 1,n 1m 1,n 1X y|(m,n)4 x y(m,n)(38)m 1, n 1 m 1,

20、n 1另外，对于变量K，有m,nm 1/ 2,nm 1/2, n1K m1,nm1,n2 1 III2 xK m,n m, nTm1,nKm, nm,nK m1,nm 1,n(39)1K11K1Kx Kym,nxKym 1,nKym 1,nPm 1,n1m1,n1 -Pm 1,n1m1,n1 +式中Pm1,n1m 1,n1(40)T m1,nKm1,n Km,n2K m1,n KJm,nKm 1, n 1 Pm 1, n 1K m 1,n将以上各式应用到式(32)，进行整理，就可以得到有限差分光束传播法的基本方 程11 11 11 11FEx m,n x m,n Pex m 1, n x m

21、1, ni1 ,i1 /FEx m, n 1 x m, n 1PEx m,n ； m,n PEx m 1,n ； m 1,nPEx m, n 1 x m, n 1QEx m 1, n 1 y m 1, n 1(41a)l 1PE m,n匚yl 1 ym,npEEy1m1,ni 1y m 1,nl 1PE m,ny1i 1 ym,n1l lPEy m,n ym,npEmy1,niy m1,npE m, n 1匚y 1i ym,n1QEym1,n 1x m 1, n 1(41b)PHx1 m,nl 1xm,n PHx1 m 1,nl 1x m1,nPHx1 m,n1i 1x m,n1PHx m,n

22、lxm, nP； m1,n ： m1,nPHx m,n1lxm,n 1qHn Xm 1,n 1ly m1,n1(41c)PHy1 m,nl 1 ym,n PHy1 mi 1i 1,nm 1,nPHy1 m,n1i 1 ym, n 1PHy m,niym, npH mH y1,n y m1,nPHy m,n1i ym,n 1qHn ym 1, n 1ix m1,n1(41d)式中i 1PE m,nCl 1Te m1,n TE12m 1,n4xm,n,lPE 1 m 1,nxTeE 1 m 1,n2i 1 /Pex m, n 1PEx m,n C 1T； m 1, n T； m 1, n 42Xm

23、, n,lP； m 1,nX 11 T； m 1,n2xP；x m,n 1Q； m 1,n 1Xm 1 m n 1 n% m 1,n 1,l% m 1,n,lT； m 1,n2% m 1,n% m 1,n % m,nP； 1匚ym,nT； 1 m, n1 T； 1 m, n 1 42 m,n,li 1 /电 m,n 1T； 1 m,n 12PE匚ym,nqEpHx1pE 1匚ym 1,ntE m, nP；Ey1,n 1m,nm, n 1p；EyT； m, n 1 41 T； m,n1,n - 2x% m 1,nm, n,l1,l1% m,n 1,lTeD1m, n1m, n2% m,n% m,

24、 nTH111 % m,nm,n 1 tH 1 m, n 1m,n,l 1%m, n,lPHx m,nQH m 1,nxpHD1pHx1m, nm, n%m, n,lpHm,nTH m,n 11,ntH1m, n 12ytH m, n 1 TH m,n 1m,n,l%m, n,l 1 %m, n,l 11,ntH m, n 12 %m,n% m, n,l% m, n 1,l%m, n %m,n 1Di m,n 丄 2i 1 力口1 1z 2 z%m,n,l町1 m,n D1 1 m, nTH1 m 1,n TH 1 m 1,n2Xm, n,l 1m, nm, n,l%m, n,l 1 %m, n,l 1%m, n,lpH m 1,nx1 TH m 1,n2xl m 1m n 1 n% m, n,lQHy m 1,n 11y4 x y% m 1,n,l2PHy m,n 12.2边界条件的选取光束传播法是在计算机的数据存储空间中对连续的实际电磁波传播过程在 空间上进行数值模拟。在电磁场的辐射、散射等问题中，边界总是开放的，即 电