Helmholtz方程在三層平板光波導中的特征計算.pdf

1、本文主要提出了一種在上下無界的平板光波導區(qū)域中計算Helmholtz方程的特征問題的有效數(shù)值解法。 對于求解無界多層均勻的光波導特征問題，以往的方法首先是簡單地設立假的邊界條件將無界的求解區(qū)域有界化，如設這些邊界條件為第一、第二或第三類邊界條件。這種做法只是將原方程特征問題粗糙地近似化，所以無論是對此近似問題進行精確求解還是離散方程后進行數(shù)值計算，解的偏差會變得很大，導致解的逼近效果變差。產(chǎn)生偏差的主要原因在于這樣的邊界條件實際

2、上并不符合無界區(qū)域中光波傳播的特性。因此，人們在邊界處改進化引入讓波向外傳播的條件，即所謂“OutgoingCondition"，從而將無界的求解區(qū)域有界化，并根據(jù)包層與芯層之間的界面條件，建立關于特征值的非線性方程。由于求解非線性方程一般采用迭代法，如牛頓法，所以其收斂性嚴重依賴于迭代初始值；然而，選取合適的初始值使迭代法收斂在數(shù)值計算上是一件非常困難之事。 因此，本文采取從Helmholtz方程出發(fā)，引入完美的匹配層(PML

3、)將無界的求解區(qū)域有界化，將方程離散化(包括層與層之間的界面條件及上下邊界條件)，得到一復矩陣，利用非對稱廣義Rayleigh商迭代具有良好的局部收斂性和快速收斂的特點，構造多重廣義Rayleigh商迭代的方法來計算矩陣特征值逼近方程特征值，并對計算誤差做了一定的分析。數(shù)值模擬計算結(jié)果表明，利用這種方法，不僅可以近似求出Helmholtz方程在無界光波導中的傳播模的特征分布，而且可以大大改善數(shù)值求解的精確度，又具有易于數(shù)值計算的優(yōu)點，如