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

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

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

3、杂结构，计算量大（矩阵运算）。（频域分析）适用范围：计算光波导的模式求解。现状：适用于较简单结构的分析。但有限差分（偏导 差分）法广泛应用于数值方法中（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 ,用差分的方法将横截面上的场离散化。这种方法已被成功地应用于分析 Y型波导及S型弯曲波导中的光波传输， 且对损耗的计算也得到了准确的结果；FD-BPM还被用于分析条形波导、三维 弯曲波导、二阶非线性效应以及有源器件。频域分析方面，同样可采用光束传播法进行分析：可采用 相关函数法获得， 还发展了一种称为虚轴光束传播法的方法，用于分析波导中的模式。其实，BPM与FDTD有不少相似的地方。其不同在于，FDTD每次都要同 时计算

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

9、Index Profile at Y=0-20-1001020X(nm)1.4545Contour Map of Index Profile at Y=0 1200“1I I- 1 11.45-1 | | i I i I | II | i | I i I i | i i ii | i I i i| ii r i | i i i 1 | i i ii i i i-i ii i i i i-20-1001020X (pm)BPM理论来源于波动方程，波动方程是建立在Maxwell方程基础Maxwell方程的一般形式为E(t)B(t)(1a)H(t)J(t)D(t)t(1b)(1c)(1d)D(t)B

10、(t)式中，E为电场强度，H为磁场强度,D为电矢量位移，B为磁感应强度，J为电流密度矢量，为体电荷密度，t为时间。对于各向同性、非磁性、电中J(t) E(t),B(t)0H(t),D(t) E(t)式中,为电导率，0为真空磁化率常数，为介电常数。将式(2)代入式(1)E(t) 0H(t) E(t)幽0 tE(t)t(3a)(3b)考虑到场对时间的依赖E(t) Eexp(i t),H(t)H exp(i t)式中，E Exi Evj Ezk,H HxiX y ZXHyjHzk为复振幅，为角频率，i为单位虚数把式(4)代入式，有0H)EE)0(H)(6)(5a) (5b)式(5a)可进一步写为将式

11、(5b)代入式(6)，有E)0()E定义复相对介电常数(8)% % i% (一二)00将式(8)代入式(7),就可以得到关于电场的矢量波方程(E)2E式中，由下式表示2o(10)22 0 0 % % ()2%C其中，C为真空光速,为真空波长。采用同样的过程，可以得到关于磁场矢(11)(12)(13)(14)(15)(16)量波方程12(H)获 ( H) 2H对于任何矢量G ,有(G) ( G) 2G又B ( °H) 0D ( E) EE从而可以进一步得到H 011E -(E) - ( %E)%把式(12)和式(16)代入到式(9)和(11),可以得到o12E (-% E) 2E 0(

12、17)2H考虑准TE模(Ez2E2H式中1 2% ( H) 2H 0 %0)和准 TM( Hz0)模，有(18)2E 0(19a)H)(一k H2H 0(19b)将式(19)写成分量形式如下2Ex2Ey2Hx% Hy(20)E*iEyj(21)HxiXHyj(22)%Ey y2EX0X(23a)%E XHx% Hx2Ey 0(23b)2Hx 0(23c)2Hy 0(23d)2.1 % HxHy1 %HyHyx % xyx% z z进一步可以写为(24d)2Ex2Ex2Ex%Ex2Ey(24a)2Ey2 _2 _2e2ej-y-2"2-xz%一y%Ex2Ex(24b)1 Hx% y2H

13、x2 x2Hx 1Hx2Hx2Hy1 Hy % x(24c)% x1 Hy2H y2y2h2 z% Hy式(24)中，Ex, xEyHx2Hy2Hx% x % yHy是空间坐标x,y, z的函数,把Ex, Ey, H x, H y随z的快速周期变化部分分离，令Exx(x,y,z)exp(iz)(25a)Eyy(x,y,z)exp( iz)(25b)Hxx(x,y,z)exp( iz)(25c)Hyy(x,y,z)exp(iz)(25d)c式中，n0为参考折射率，选择时应尽量接近导模的有效折射率，否则会影响计算精度。 x , y, x，y为包络函数，如图1所示。x, y, x, yEx，Ey，H

14、图2-1包络函数示意图将式(25)代入式(24)可得包络函数的矢量波方程22x x 2i22y z12(28a)(27a)2y-2x2y""2z2i(27b)222ix x2-2-x z(27c)1» y% z z22% -yT T 2ix % x yz2221%212 i -x% -% z y x y x % y进一步整理得(27d)2T 2i xAxx xA<yyzz22i(28b)2x-2"z2i1 7o% z zBxxBxy(28c)2iB B yy y yx x(28d)式中x2x 222xyA yy y1_%、Ayx2x2x1 % i%

15、 z2BxyByy三维半矢量形式Byx忽略x, y场之间的耦合，则有2 x2z2ix2y ""2 z2x2z2i2i2iBxx xByy3.二维标量形式（忽略场的方向性）2-V 2i - A z z22i z4.二维半矢量形式与三维半矢量类似，但是折射率分布更简单5.二维标量形式2.1.2方程离散数值处理（1）纵向数值处理直接求解方程式（27）是非常困难的，因此需要对它进行离 散化处理，通过数值方法来求解。BPM的数值离散化处理方法很多，这里采用有限差分方法来实现，有限差分法的核心就是把导数写成差分的形式。为了便于方程的求解，有必要对方程进行近似处理，在纵向 （即沿着光的传

16、 播方向Z）的近似处理有缓变包络近似（SVEA）、广角近似等。在这里，我们采用 缓变包络近似方法。如果包络函数随Z的变化足够缓慢，使得20z则有22- 2i 2i （29）Z Z Z其中，代表上述各个包络函数。纵向处理主要是解决传播方向上相邻两个截面上场的关系问题，从式 （29）可 以看出，纵向处理就转化为对Z的一阶偏微分处理。设相邻的两个截面分别用l 和l 1标志，第l个截面上的场为已知，第l 1截面上的场待求，两截面的问 距为 Z,见图。Z图2相邻截面差分格式示意图式(29)中右边的偏微分项可差分为(30)式28可表小为(31)即：11 zf若右边值已知，则由1面上的场场分布,就可以获得1

17、+1面上的场分量1 1类似的方法还有：稳定性：是指计算过程中积累误差是无限增加还是可以控制， 对于沿z方向折射率缓变的情况,0.5时，上述差分格式是稳定的。数值损耗：由数值计算引起的沿传输方向上的能量损失，是非物理损耗。研究表明,0.5时，数值损耗最小,1时，数值损耗最大。因此，综合考虑到数值计算的稳定性及数值损耗， 在计算中，要仔细选择合适的 值，使得稳定性和数值损耗都可以接受。这样，式(28)可进步化为AxxAyyi y(32a)ixBxxBxxByyByy(32b) Ayy(32c)B iyx x (32d)2i其中2iz%缓变包络近似特点:最早提出，方法最简单 广角近似 (Pad e近

18、例特点：近似更少，更高精度(2)横向数值处理图3所示的是九点差分网格，横向数值处理就是图3有限差分网格结构在有限差分光束传播法中，横向处理可采用上述的九点差分格式，设要差分的变量为 (x, y),则在点(m,n)上，的一阶和二阶导数可差分为x (m,n)y (m,n)(m 1,n) (m 1,n)2 x(33)(m,n 1) (m,n 1) 2 y(m,n)(34)(m 1,n) (m 1,n) 2 (m,n)2x2y (m,n)(m,n)(35)(m,n 1) (m,n 1) 2 (m,n)2y(36)1 ， m 1,n 1 m 1,n 14 x ym 1,n 1m 1,n 1(37)(m,

19、n)m 1,n1,n 1m 1,n1,n 1(38)另外，对于变量K ,有m,nm 1/ 2,nm 1/2,n1,n1,n1,nK m,nm, n1,nm, nm,n1,n1,n(39)m,nm 1,nm 1,n1,nm 1,n 1Pm 1,n 1 m 1,n 1 -P m 1,n 1m 1,n 1 +式中P m 1,n 1 m 1,n 1(40)Km 1,n K m,nT m 1,n2K m 1,n K m,nK m 1,n 1P m 1,n 1K m 1,n将以上各式应用到式（32）,进行整理，就可以得到有限差分光束传播法的基本方11PEx m,ni 1x m,npE 1xi 11,n x

20、 m 1,ni 1/i 1/PEx m,n 1x m,n 1i - iiPex m,n x m,nPex m 1,nm 1,nPEx m,n 1 x m,n 1QE m 1,n 11V m 1,n 1xyi 1i 1i 1PE m,n v m,nPE mEyyEyi 1i 1PEm,n 1 v m,n 1yyiiiPE m,n v m,n PE m 1,ny Jyy/pE m,n 1' m,n 1yy(41a)i 11,n y m 1,nm 1,nQEy m 1,n 1 x m 1,n 1(41b)l 1Phx m,nxm,nPHx1,n1,nl 1PHx m,nlPhx m,nPHx

21、 m,n xm, nlx1m, nlPHx m1,n1,nm,n 1qH x1,n1,n(41c)l 1PHy m,nl 1PHy m,nm, nl 1PHym1,n1,nm, nlPHy m,nPHy m,nm, nl yPHym,n式中l 1PEx m,nPE1xqH1,nl 1TE mm 1,nl 1PEx m,n1,n1,n1,n1,n(41d)m,n,ltE1TE121,n1,npExm,n CTE m 1,nTE m 1,n4PExqEx1,nTepE 1ym,n11Te m, n1,nPEx m,n“ I2m, n,1i1Te m 1,nm 1 m n 1 n% m 1,n 1,

22、l1,n,l2% m 1,n1,n% m 1,n% m,ni 11Te m,n 142m, n,ll 1PEy m,n 1i 1Te m, npEym,nqEPEyi TE m, npE匚y1,n 1Tem 1,nm, n 1pEEym,nm,nD11m, niTE m, n 142m, n,11 tE m,n1,n-2x2% m,n% m,ntH1% m 1,n% m,n 1,l1% m,ni 1m, n 1 Thm,n 12m,n,1 1%m,n,l 2 %m,n,l%m,n,l 12 zi iPhx m 1,n2xxl 11 1 d Th m,n 1Phx m,n 12yPHX m,ni

23、D m, n 1iiTh m,n 1Th m,n 12y22122m,n,1 xi %m,n,1 1%m,n,1 1%m,n,12 zPH m 1,n xPHx m,n 111 Th m, n 12yQHX m 1,n 1 xm 1 m n 1 n% m,n,1% m, n 1,12 %m,nTh m,n 1%m, n%m,n 1D1 m, n12i z%m,n,l 1 %m,n,l 12 z%m,n,li iPHymEd11m, ni 1iThm 1,n Th1 m 1,npHlPHy m,nm,n,l 12i%m,n,lPHy1m,n1,niD m, n%m,n,l 2%m, n,l2 zi 1TH m 1,ntH m 1,nTH m 1,nm, n,l%m,n,l 1 %m,n,l 1%m,n,lpH m,n 1yP： m 1,nxlTH m 1,n2Xlm 1 m n 1 n % m, n,lQh m 1,n 1 1 y4 x y% m 1,n,lTh m 1,n2 %m,n%m, n%m 1,n2.2边界条件的选取光束传播法是在计算机的数据存储空间中对连续的实际电磁波传播过程在 空