Rsoft软件(BPM-光束传播法模拟)中文课件(第四章-光束传播法基础-)江苏大学陈明阳教授编写-

江苏大学陈明阳光通信器件设计第四章光束传播法基础江苏大学光电信息科学与工程系陈明阳2018.9江苏大学陈明阳第一节电磁场数值计算方法简介1．电磁场数值计算电磁场数值计算是一种基于麦克斯韦方程组，通过建立逼近实际工程电磁场问题的连续型的数学模型，然后采用相应的数值计算方法，经离散化处理，把连续型数学模型转化为等价的离散型数学模型，计算出待求离散数学模型的离散解(数值解)，从而获得相应结果的一种方法。江苏大学陈明阳2.常用电磁场数值方法简介（1）有限差分法(频域有限差分法)有限差分法是利用划分网格的方法将定解区域离散化为网格离散节点的集合，然后基于差分原理，以各离散点上函数的差商来近似替代该点上的偏导数，这样待求的偏微分方程定解问题可转化为一组相应的差分方程的问题。根据差分方程组，解出各离散点上的待求函数值，即为所求定解问题的离散解。江苏大学陈明阳方法特点原理:偏导差分方法特点：原理简单、通用性好；对复杂结构，计算量大(矩阵运算)。(频域分析)适用范围：计算光波导的模式求解。现状：适用于较简单结构的分析。但有限差分(偏导差分)原理被广泛应用于各种数值计算江苏大学陈明阳(2)有限元法有限元法是以变分原理为基础，把所要求解的微分方程转化为相应的变分问题，即泛函求极值问题。将分析的区域划分为很多三角形(每个三角形成为一个基元)，每个基元内的场用多项式来表达，然后加入不同基元间场的连续条件，就可得到整个横截面的场分布。江苏大学陈明阳网格化单元示例：江苏大学陈明阳方法特点较复杂---需要前处理(三角化，剖分)；后处理：（场分布，伪解剔除）（通用性强，精度高）根据该方法对于各种各样的电磁计算问题具有较强的适应能力性，所形成的代数方程矩阵求解容易、收敛性好。对于形状和分布复杂的三维问题，由于其变量多和剖分要求细，往往因计算机内存而受到限制。程序设计复杂、计算量较大。适用范围：求解光波导的模式（有效折射率、色散、双折射、传输损耗等）。江苏大学陈明阳现状：功能昀强大的数值方法之一。特别是上世纪90年代出现的矢量有限元方法，完全解决了有限元方法出现的伪解问题，大大降低了有限元法的后处理过程。发展：有限元光束传播法。软件：COMSOLMULTIPHYSICS江苏大学陈明阳(3)时域有限差分法时域有限差分法是近年来开始流行的一种数值模拟方法，它通过将麦克斯韦方程在时间空间上离散化的方法实现对电磁波传播的模拟。它能够得到电磁波传输的瞬态（即时域）信息，通过傅里叶变换即可得到相应的频域信息。江苏大学陈明阳时域有限差分法由K.S.Yee于1966年首先提出，此后经过众多学者的努力，使之不断完善，现已比较成熟。计算过程为：设置初始场，然后依时间步推进计算，并在每一时间步交替地计算每一离散点的电场和磁场。江苏大学陈明阳特点：不需要矩阵运算，只需简单的加减乘除运算由前一时刻的场来获得下一时刻场的值。而且，非常适合于并行计算。适用范围：计算光波导的模场分布、有效折射率；研究波导之间的连接、耦合问题。注：主要用于一维和二维光波导的分析。三维波导分析计算量稍大。现状：ADIFDTD，可应用于各向异性介质，非线性介质，PML吸收边界江苏大学陈明阳(4)光束传播法光束传播法(BeamPropagationMethod，简写BPM)是目前光波导器件研究与设计领域昀流行的方法之一，其基本思想是在给定初始场的前提下，一步一步地计算出各个传播截面上的场。江苏大学陈明阳应用范围这种方法已被成功地应用于分析Y型波导及S型弯曲波导中的光波传输，且对损耗的计算也得到了准确的结果；FD-BPM还被用于分析条形波导、三维弯曲波导、二阶非线性效应以及有源器件。频域分析方面，同样可采用光束传播法进行分析：即虚轴光束传播法或基于相关函数法的光束传播法。江苏大学陈明阳方法特点BPM与FDTD的区别：FDTD每次都要同时计算整个波导的模场，而BPM只算一个面。计算量较小，应用范围非常广泛适用范围：计算光波导的模式、色散、双折射、传输损耗等；分析波导传输、连接、耦合，光栅的传输特性等。江苏大学陈明阳第二节有限差分光束传播法基本原理光束传播法(BPM)的基本思想就是把波导沿着传播方向剖分成若干个截面，根据前一个或几个截面上的已知场分布得到下一个截面上的场分布。江苏大学陈明阳1剖分示例锥形波导结构图江苏大学陈明阳计算参数(z方向步长取为10微米)江苏大学陈明阳计算时波导的实际折射率分布江苏大学陈明阳计算参数(z方向步长取为100微米)江苏大学陈明阳计算时波导的实际折射率分布江苏大学陈明阳2光束传播法公式推导BPM理论来源于波动方程，波动方程是建立在Maxwell方程基础上的。推导目的：了解BPM计算要求和限制江苏大学陈明阳Maxwell方程的一般形式()()0tttBE()()()ttttDHJ()tD()0tB(1a)(1b)(1c)(1d)江苏大学陈明阳(1)减少矢量数0()(),()(),()()ttttttJEBHDE(2)0()()0tttHE()()()ttttEHE(3a)(3b)将式(2)代入式(1)有：江苏大学陈明阳(2)去掉时间依赖项考虑到场对时间的依赖()exp(),()exp()tittitEEHH(4)00iEH()iHE把式(4)代入式(3),有(5a)(5b)江苏大学陈明阳(3)减少单个方程的变量数利用矢量公式及方程(5)，可得结果：20()()0iEE(6)00()riii22220022()c再定义：江苏大学陈明阳得到:电场的矢量波方程采用同样的过程，可以得到关于磁场矢量波方程2()EE(7a)21()()()HHH(7b)江苏大学陈明阳利用矢量公式进行化简，得到22rr1()0EEE221()0HHH(8a)(8b)江苏大学陈明阳(4)去除矢量的z方向标量0zE0zH221()0rrEEE2211()()0zzHHkkHH考虑准TE模()和准TM()模，有:(9a)(9b)江苏大学陈明阳这里：xyijxyEEEijxyHHHij江苏大学陈明阳展开得到标量方程形式2210rrxxyxrEEEExxy2210rryxyyrEEEEyxy22110yxxxxHHHHHyxyzz22110yyxyyHHHHHxyxzz(10a)(10b)(10c)(10d)江苏大学陈明阳222222110xxrxxryryrEEEExxyzEExyxy(10a)进一步可以写为(11a)江苏大学陈明阳(5)缓变包络近似把Ex、Ey、Hx、Hy随z的快速周期变化部分分离(,,)exp()xxExyziz(,,)exp()yyExyziz(,,)exp()xxHxyziz(,,)exp()yyHxyziz(12a)(12b)(12c)(12d)江苏大学陈明阳缓变包络近似原理示意图xyxyE,E,H,Hxyxy,,,江苏大学陈明阳222212xxxrxrixxyzz22210yxryrxyxy将式(12)代入式(11)可得包络函数的矢量波方程:(13a)江苏大学陈明阳整理后得到:222xxxxxxyyiAAzz(14a)此即推导得到的昀终公式江苏大学陈明阳三维半矢量形式222xxxxxiAzz忽略x，y场之间的耦合，则有(15a)江苏大学陈明阳222iAzz三维标量形式(忽略场的方向性)相应的还有：二维半矢量形式与三维半矢量类似，但是折射率分布更简单二维标量形式(16a)江苏大学陈明阳(6)方程离散化处理220z2222iizzz在纵向(即沿着光的传播方向z）的近似处理有缓变包络近似(SVEA)、广角近似等。如果包络函数随z的变化足够缓慢，使得则有(17a)(一)纵向数值处理采用有限差分方法，其核心就是把导数写成差分的形式。(Φ代表包络函数。)江苏大学陈明阳因此(14a)式可简化为:222xxxxxxyyiAAzz2xxxxxyyiAAz(18a)江苏大学陈明阳l1lzZ设相邻的两个截面分别用l和l+1标志，两截面的间距为Δz相邻截面差分格式示意图江苏大学陈明阳1llzz1lllfz1lllzf式(18)中右边的偏微分项可差分为式18可表示为(20)若右边值已知，则由l面上的场分布，就可以获得l+1面上的场分量。(19)即：(21)江苏大学陈明阳这样，式(18)可进一步化为11lllxxxxxxxyyCACAA(22a)，江苏大学陈明阳(二)横向数值处理目的:解决横向(x、y方向)偏导问题xyy(m-1,n+1)(m,n+1)(m+1,n+1)(m,n)(m+1,n)(m-1,n)(m-1,n-1)(m,n-1)(m+1,n-1)x有限差分网格结构江苏大学陈明阳差分公式：(,)(1,)(1,)2mnmnmnxx(,)(,1)(,1)2mnmnmnyy222(,)(1,)(1,)2(,)mnmnmnmnxx(23)(24)(25)江苏大学陈明阳222(,)(,1)(,1)2(,)mnmnmnmnyy2(,)11,11,14mnmnmnxyxy1,11,1mnmn2(,)11,11,14mnmnmnyxxy1,11,1mnmn(26)(27)(28)江苏大学陈明阳将以上各式应用到式(22a)，就可以得到有限差分光束传播法的基本方程111111,,1,1,,1,1xxxllllExExllExPmnmnPmnmnPmnmn,,1,1,,1,1xxxllllExExllExPmnmnPmnmnPmnmn(29a)江苏大学陈明阳光束传播法计算步骤设置折射率分布设置波源(初始场分布):高斯型，模式场等步进求解（三角化，递推法）后处理[有效折射率，传输损耗，色散，双折射，模式耦合]江苏大学陈明阳总结1.每次只能获得一个波长的结果2.计算得到的场分布是场的强度分布3.不考虑光的反射问题4.要求光近似为直线传播（Ez=0,Hz=0;缓变包络近似)5.步长影响计算结果6.矢量类型影响计算结果（弱导：标量，非弱导：半矢量，全矢量（场耦合））江苏大学陈明阳参考资料陈明阳，新型光子晶体光纤的数值模拟，燕山大学博士学位论文，2004江苏大学陈明阳江苏大学陈明阳第三节基本操作学习内容：1监视器（及路径）设置2初始场设置3BeamPROP计算