第二章光束传播法基本原理.docx
《第二章光束传播法基本原理.docx》由会员分享,可在线阅读,更多相关《第二章光束传播法基本原理.docx(34页珍藏版)》请在冰豆网上搜索。
![第二章光束传播法基本原理.docx](https://file1.bdocx.com/fileroot1/2023-3/31/ca7ec9a6-b8f6-488b-870d-91a104595a95/ca7ec9a6-b8f6-488b-870d-91a104595a951.gif)
第二章光束传播法基本原理
第四章光束传播法基础
第一节数值计算方法
1.电磁场数值计算
它是一种基于麦克斯韦方程组,建立逼近实际工程电磁场问题的连续型的数学模型,(合理的假设)然后采用相应的数值计算方法,经离散化处理,(合适
的方法,使离散化的模型既能反映连续型模型的特性,又便于计算机分析)把连
续型数学模型转化为等价的离散型数学模型,计算出待求离散数学模型的离散—解(数值解),从而获得相应结果的一种方法。
2•数值方法分类:
时域分析、频域分析。
时域分析:
模拟光在波导中的传播过程
频域分析:
求解波导模式
时域分析逼真:
把原来因为速度太快、结构太小、不可见的现象模拟出来,能够直观地展示。
求解:
波导连接、耦合、非线性特性、波导模式。
频域分析:
光场分布、给定具体结构波导的模式的有效折射率(色散、偏振)、损耗(材料吸收、结构本身导致)等。
问题:
频域结果能否推得时域信息?
反之?
3.常用数值方法简介
(1)有限差分法(频域有限差分法)
(20世纪50年代出现)利用划分网格的方法将定解区域离散化为网格离散节点的集合,然后基于差分原理,以各离散点上函数的差商来近似替代该点上的偏导数,这样待求的偏微分方程定解问题可转化为一组相应的差分方程的问题。
根据差分方程组,解出各离散点上的待求函数值,即为所求定解问题的离散解,再应用插值方法便可从离散解得到定解问题在整个场域上的近似解。
原理:
偏导差分
方法特点:
原理简单、通用性好;对复杂结构,计算量大(矩阵运算)0(频域分
析)
适用范围:
计算光波导的模式求解。
现状:
适用于较简单结构的分析。
但有限差分(偏导差分)法广泛应用于数值方法中
(2)有限元法
20世纪40年代提出,其在电磁问题方面的应用有约40多年历史。
以变分原理为基础,把所要求解的微分方程转化为相应的变分问题,即泛
函求极值问题。
常见方法为把要分析的区域划分为很多三角形(每个三角形成为
一个基元),每个基元内的场用多项式来表达,然后加入不同基元间场的连续条件,就可得到整个横截面的场分布。
特点:
较复杂---需要前处理(三角化,剖分);后处理:
(场分布,伪解剔除)(通用性强,精度高)根据该方法对于各种各样的电磁计算问题具有较强的适应能力性,所形成的代数方程矩阵求解容易、收敛性好。
主要缺点:
对于形状和分布复杂的三维问题,由于其变量多和剖分要求细,往往因计算机内存而受到限制。
程序设计复杂、计算量较大。
适用范围:
求解光波导的模式(有效折射率、色散、双折射、传输损耗等)。
现状:
功能最强大的数值方法之一。
特别是上世纪90年代出现的矢量有限元方
法,完全解决了有限元方法出现的伪解问题,大大降低了有限元法的后处理过程。
有限元光束传播法。
(3)时域有限差分法
时域有限差分法是近年来开始流行的一种数值模拟方法,它通过将麦克斯
韦方程在时间空间上离散化的方法实现对电磁波传播的模拟。
它能够得到电磁波传输的瞬态(即时域)信息,通过傅里叶变换即可得到相应的频域信息。
时域有限差分法由K.S.Yee于1966年首先提出,此后经过众多学者的努力,使之不断完善,现已比较成熟。
但是,在许多方面它仍在继续发展,解决问题的能力和应用范围仍在不断地提高和扩大。
计算过程为:
设置初始场,然后依时间步推进计算,并在每一时间步交替地计算每一离散点的电场和磁场。
特点:
不需要矩阵运算,只需简单的加减乘除运算由前一时刻的场来获得下一时刻场的值。
而且,它还非常适合于并行计算,这正好与当今计算机的发展趋势相吻合,这就更加提高了时域有限差分法解决实际复杂问题的能力。
适用范围:
计算光波导的模场分布、有效折射率;研究波导之间的连接、耦合问题。
注:
主要用于一维和二维光波导的分析。
三维波导分析计算量稍大。
现状:
ADIFDTD,可应用于各向异性介质,非线性介质,PML吸收边界
⑷光束传播法(BeamPropagationMethod,简写BPM)
光束传播法是目前光波导器件研究与设计领域最流行的方法之一,其基本
思想是在给定初始场的前提下,一步一步地计算出各个传播截面上的场。
光束
传播法最早是由M.D.Feit等人于1978年研究光场及大气激光束传播时提出的。
最早的BPM是以快速傅里叶变换(Fast-FourierTransform,称FFT)为数学手段实现的,称为FFT-BPM。
FFT-BPM源于标量波方程,只能得到标量场(即只能处理一个偏振分量),不能分辨出场的不同偏振(TE模或TM模)以及场之间的耦合。
由于上述缺点,D.Yevick等人于1989年提出了一种新方法一有限差分光束传播法FD-BPM,用差分的方法将横截面上的场离散化。
这种方法已被成功地应用于分析丫型波导及S型弯曲波导中的光波传输,且对损耗的计算也得到了准确的结果;FD-BPM还被用于分析条形波导、三维弯曲波导、二阶非线性效应以及有源器件。
频域分析方面,同样可采用光束传播法进行分析:
可采用相关函数法获得,还发展了一种称为虚轴光束传播法的方法,用于分析波导中的模式。
其实,BPM与FDTD有不少相似的地方。
其不同在于,FDTD每次都要同时计算整个波导的模场,而BPM只算一个面特点:
计算量较小,应用范围非常广泛
适用范围:
计算光波导的模式、色散、双折射、传输损耗等;分析波导传输、连接、耦合,光栅的传输特性等。
4•数值方法发展趋势:
方法融合现象明显(有限元法与光束传播法的结合形成了另外一种方法一有
限元光束传播法(FE-BPM)。
)、相互推动(PMLFDTD,BPM,FEM)。
第二节有限差分光束传播法基本原理
光束传播法(BPM)的基本思想就是把波导沿着传播方向剖分成若干个截面,
根据前一个或几个截面上的已知场分布得到下一个截面上的场分布
BeaBPROPSisulationParaaeters
DomainMin;
DomainMax:
ComputeStep:
SliceStep:
MonitorStep:
CurrentValue^28~
DefaultUse
ValueDefs
0.2
0.2
CurrentValue^24~[2i
DefaultUseVilueDefs
Current
Value
|5-
[U2Q~
DefaultUseValueDefsr~11120
10
100
lio-
|100L
EstimatedTime:
0.088min
SaveSettings
OK
Cancel
ContourMapofIndexProfileatY=0
1000-
800-
600-
■
■
400-
■
■
200-
■
■
-20-10
01020
X(nm)
IISIbuIationPara*eterE—ComputeIndexProTile
X
Curent
Value
Default
Value
UseDeFs
DomwMir:
17
DomainMax:
|27.2
P
ComputeStep:
[02-
j0.2
审
Slice5切Moriit&rStep:
CurrentDefaultUteValueValueDeFs
CurrentDefaultUseValueValueDefs
|卫2|-2':
217
[0[0|7
f23.2[.2破
[1100
|7
[El-P
[C2~pn~P
[ioo
|100
|1QO
ior
17
100|7
S/mbals...
Di^plavM&de
OutputFiePiefix:
EstimatedTimer
0.006min
IContauiMap^)亍
SaveSettings
Display...
□utpul...
_"5k……"'1
!
■■■»■■■■■kkb■■■■■■■■■■■■
Cancel
ContourMapofIndexProfileatY=0
1200
woo-
800
600
400
200
0
-20
10
20
£=L
N
0
X(pm)
BPM
理论来源于波动方程,波动方程是建立在
Maxwell方程的一般形式为
E(t)畀
1.4545
Maxwell方程基础上的。
(1a)
H(t)
J(t)
D(t)
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)o
0t
(3a)
E(t)
H(t)E(t)
t
(3b)
考虑到场对时间的依赖
E(t)
Eexp(it),H(t)
Hexp(it)
⑷
式中,
EExi
EyjEzk,HHx
iHyjHzk为复振幅,
为
角频率,i为单位虚数。
把式(4)代入式⑶,有
Ei0H
H(i
0
)E
(5a)
(5b)
式(5a)可进一步写为
(
E)i0(
H)0
(6)
将式(5b)代入式⑹,有
(
E)20(
i)E0
⑺
定义复相对介电常数
式中,由下式表示
量波方程
对于任何矢量G,有
从而可以进一步得到
(15)
把式(12)和式(16)代入到式(9)和(11),可以得到
21
2E(%%E)
212
2H%(H)2H0(18)
%
考虑准TE模(Ez0)和准TM(Hz0)模,有
212
2E(—%E)2E0(19a)
%
211%2
H%(H)^k(:
kH)H0
(19b)
式中
■
ix
■jy
(20)
E
Exi
Eyj
(21)
H
Hxi
Hyj
(22)
将式(19)写成分量形式如下
(23a)
(23b)
%Hy
Hx
%Hx
2Hx
2Hy
进一步可以写为
1%Hy
%z
z
2.
2.
笔X
Ex
2
2
y
z
2Ey
0
xy
2Ey
2Ey
2
2
x
z
x
0
%Ey
'Ey
Hy
%Ey
2Ex
2Ex
%Ex
(24a)
(24b)
%—
y
式(24)中,Ex,
2Hx
2Hx
1
%Hx
2x
2z
%
z
z
2Hx-
2Hy
%—y
1
H
yx
%x
化
1
%
Hy
2
y
z2
%
z
z
H-
2Hx
%—
x
1
H
xy
%y
yx
y
(24c)
y,z的函数,
Hy是空间坐标x,
%Ex
1Hy
%x
1Hx
%y
把Ex,Ey,Hx,Hy随z的快速周期变化部分分离,令
Ex
x(x,y,z)exp(i
z)
(25a)
Ey
y(x,y,z)exp(i
z)
(25b)
Hx
x(x,y,z)exp(i
z)
(25c)
Hy
y(x,y,z)exp(i
z)
(25d)
n。
(26)
c
式中,n°为参考折射率,选择时应尽量接近导模的有效折射率,否则会影响计算精度。
x,y,x,y为包络函数,如图i所示。
EEHH
xyxy
图2-1包络函数示意图
将式(25)代入式(24)可得包络函数的矢量波方程
x%x°%x
22
xx2i
22
yz
2
y
""2
z
2i
%—
y
%—
x
2
x
2
z
进一步整理得
2i
y
u
x
(27b)
1
%
x
%
z
z
1
%
y
x
0
(27c)
1
%
y
%
z
z
1
%
x
y
0
(27d)
Axx
x
Axyy
x
2i
%
y
2i
%
x
x
2
x
式中
x
Ayy
z2
2i
(28b)
2
x
~2~"
z
2i
Bxx
Bxy
(28c)
Ayx
2i
Byy
B
yxx
(28d)
1__
%x
2
x
2
y
1
2
x
2
x
1%
■
i
%z
Bxy
Byyy%
三维半矢量形式
Byx
忽略x,y场之间的耦合,则有
2
x
2
z
2i
Axx
2i
2
x
2
z
2i
B
XXX
2
y
2
z
2i
B
yyy
3.三维标量形式
(忽略场的方向性)
2i
2i
4.二维半矢量形式
与三维半矢量类似,但是折射率分布更简单
5.二维标量形式
2.1.2方程离散数值处理
(1)纵向数值处理直接求解方程式(27)是非常困难的,因此需要对它进行离散化处理,通过数值方法来求解。
BPM的数值离散化处理方法很多,这里采用有限差分方法来实现,有限差分法的核心就是把导数写成差分的形式。
为了便于方程的求解,有必要对方程进行近似处理,在纵向(即沿着光的传播方向Z)的近似处理有缓变包络近似(SVEA)、广角近似等。
在这里,我们采用缓变包络近似方法。
如果包络函数随Z的变化足够缓慢,使得
2
齐0
z
则有
2
2i—2i—㈣
ZZZ
其中,代表上述各个包络函数。
纵向处理主要是解决传播方向上相邻两个截面上场的关系问题,从式(29)可
以看出,纵向处理就转化为对Z的一阶偏微分处理。
设相邻的两个截面分别用I
和I1标志,第I个截面上的场为已知,第I1截面上的场待求,两截面的间
距为Z,见图
图2相邻截面差分格式示意图
式(29)中右边的偏微分项可差分为
(30)
式28可表示为
(31)
类似的方法还有:
稳定性:
是指计算过程中积累误差是无限增加还是可以控制,对于沿z方向
数值损耗:
由数值计算引起的沿传输方向上的能量损失,是非物理损耗。
研究表明,0.5时,数值损耗最小,1时,数值损耗最大。
因此,综合考虑到数值计算的稳定性及数值损耗,在计算中,要仔细选择合
适的值,使得稳定性和数值损耗都可以接受。
这样,式(28)可进
C
步化为
Axx
l1
x
C
1
Axx
xAxyy
(32a)
C
Ayy
l1y
C
1
Ayy
yAyxx
(32b)
D
Bxx
l1
x
D
1
Bxx
xAxyy
(32c)
D
Byy
l1y
D
1
Byy
lBlyyxx
(32d)
其中
C
2i
D
2i
1
%
z
z
z%
z
缓变包络近似特点:
最早提出,方法最简单
广角近似(Pade近似特点:
近似更少,更高精度
图3有限差分网格结构
在有限差分光束传播法中,横向处理可采用上述的九点差分格式,设要差分的
变量为(x,y),则在点(m,n)上,的一阶和二阶导数可差分为
(m1,n)(m1,n)
X(m,n)2X
(33)
_(m,n1)(m,n1)
y(m,n)2y
(34)
2
(m1,n)
(m1,n)
2
(m,n)
2
2
X
(m,n)
X
(35)
2
(m,n1)
(m,n1)
2
(m,n)
2
2
y
(m,n)
y
(36)
2
1
1
m1,n1
m1,n1
Xy|(m,n)
4xy
(m,n)
(38)
m1,n1m1,n1
另外,对于变量K,有
m,n
m1/2,n
m1/2,n
1
Km
1,n
m
1,n
21III
2x
Km,nm,n
T
m
1,n
K
m,n
m,n
Km
1,n
m1,n
(39)
1
K
1
1
K
1
K
xK
y
m,n
x
K
y
m1,n
K
y
m1,n
P
m1,n
1
m
1,n
1-
P
m1,n
1
m
1,n
1+
式中
P
m
1,n
1
m1,n
1
(40)
Tm
1,n
K
m
1,nK
m,n
2Km
1,nK
J
m,n
Km1,n1Pm1,n1
Km1,n
将以上各式应用到式(32),进行整理,就可以得到有限差分光束传播法的基本方程
11111111
FExm,nxm,nPexm1,nxm1,n
i1,i1/
FExm,n1xm,n1
PExm,n;m,nPExm1,n;m1,n
PExm,n1xm,n1
QExm1,n1ym1,n1
(41a)
l1
PEm,n
匚y
l1y
m,n
pE
Ey
1
m
1,n
i1
ym1,n
l1
PEm,n
y
1
i1y
m,n
1
ll
PEym,ny
m,n
pE
m
y
1,n
i
ym
1,n
pEm,n1
匚y1
iy
m,n
1
QEy
m
1,n1
xm1,n1
(41b)
PHx1m,n
l1
x
m,nPHx
1m1,n
l1
xm
1,n
PHx1m,n
1
i1
xm,n
1
PHxm,n
l
x
m,n
P;m
1,n:
m
1,n
PHxm,n
1
l
x
m,n1
qH
nX
m1,n1
l
ym
1,n
1
(41c)
PHy1m,n
l1y
m,nPHy1m
i1
i1,n
m1,n
PHy1m,n
1
i1y
m,n1
PHym,n
i
y
m,n
pHm
Hy
1,nym
1,n
PHym,n
1
iy
m,n1
qH
ny
m1,n1
i
xm
1,n
1
(41d)
式中
i1
PEm,n
C
l1
Tem
1,nTE1
2
m1,n
4
x
m,n,l
PE1m1,n
x
TeE1m1,n
2
i1/
Pexm,n1
PExm,nC1
T;m1,nT;m1,n4
2
X
m,n,l
P;m1,n
X1
1T;m1,n
2
x
P;xm,n1
Q;m1,n1
X
m1mn1n
%m1,n1,l
%m1,n,l
T;m1,n
2%m1,n
%m1,n%m,n
P;1
匚y
m,n
T;1m,n
1T;1m,n14
2m,n,l
i1/
电m,n1
T;1m,n1
2
PE
匚y
m,n
qE
pHx1
pE1
匚y
m1,n
tEm,n
P;
Ey
1,n1
m,n
m,n1
p;
Ey
T;m,n14
1T;m,n
1,n-2
x
%m1,n
m,n,l
1,l
1
%m,n1,l
Te
D1
m,n
1
m,n
2%m,n
%m,n
TH1
1
1%m,n
m,n1tH1m,n1
m,n,l1
%m,n,l
PHxm,n
QHm1,n
x
pH
D1
pHx1
m,n
m,n
%m,n,l
pH
m,n
THm,n1
1,n
tH1
m,n1
2
y
tHm,n1THm,n1
m,n,l
%m,n,l1%m,n,l1
1,n
tHm,n1
2%m,n
%m,n,l
%m,n1,l
%m,n%m,n1
Dim,n丄2i1力口①11
z2z%m,n,l
町1m,nD11m,n
TH1m1,nTH1m1,n
2
X
m,n,l1
m,n
m,n,l
%m,n,l1%m,n,l1
%m,n,l
pHm1,n
x
1THm1,n
2
x
lm1
mn1n
%m,n,l
QHym1,n1
1
y
4xy
%m1,n,l
2
PHym,n1
2.2边界条件的选取
光束传播法是在计算机的数据存储空间中对连续的实际电磁波传播过程在空间上进行数值模拟。
在电磁场的辐射、散射等问题中,边界总是开放的,即电