衍射光栅数值模拟的模展开方法及Chebyshev配置Dirichlet-to-Neumann映射方法

来源 :中国科学技术大学 | 被引量 : 0次 | 上传用户:hfx285306638
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
衍射光栅是一种剧期结构,具有很多实际的应用,例如单色仪,光谱仪,激光器,波分复用器,光脉冲压缩设备及其他光学仪器等方面。数值方法在衍射光栅的设计,分析和结构优化等方面是必不可少的。理论上,当光栅结构在一个方向具有周期性的时候,可以用标准的数值方法,如有限元方法来解决问题。然而运用这些方法会产生一些很大的,复值的,不正定的线性方程组,求解这些方程组要用很多的计算时间。其他一些利用光栅结构现有的几何特性的方法反而更为有效。在衍射光栅方面现有的方法包括解析模展开法,Fourier模展开法,有限差分模展开法,差分法和积分方程法等等。所有的模展开法都要求结构具有均匀层状,这样波场才能在每一层用特征根和特征向量的展开形式来表示。在应用这些方法的过程中,求解每一层中的特征值问题是花费计算时间最多的一部分。   在本文中,首先介绍解析模展开法和Fourier模展开法。Fourier模展开法是用Fourier序列展开法来解特征根问题的,由于较容易实现,所以Fourier模展开法是很常用的一种方法。然后给出一个四阶有限差分的模展开法。所有的这些模展开方法都需要求解特征根问题。由于解特征根问题相对的用时较多,所以对于均匀层状的衍射光栅,我们给出一种Dirichlet-to-Neumann(DtN)算子法。为了避免在每一层解特征根问题,我们构造一个算子,这个算子把区域内的波场映射到它在每一层边界的法向导数上去。在实际应用中,这个算子被称为DtN算子,用一个矩阵来估计它。接下来在均匀层状方向用一个Chebyshev配置方法离散,在周期方向用一个四阶有限差分格式来离散,最后就可以高效地解出这个算子。在以前的工作中,DtN算子法被用于周期层状的圆柱体结构和分段均匀排列韵咣敲导。对于圆柱体结构,DtN算子用圆柱波来构造。对于分段均匀排列的光波导,在横截面方向,用Chebyshev配置法和二阶有限差分法来离散,从而构造DtN算子。在四阶有限差分的模展开法中,我们用了同样的四阶有限差分格式来离散周期方向。通过数值算例可以看出,DtN算子法比Fourier模展开法效率更高,因为这种方法避免了求解特征根方程问题,解DtN算子用时更少。最后,我们将Chebyshev配置DtN算子法拓展到衍射光栅的斜入射波问题。
其他文献
工程设计,最优控制,信息技术以及经济均衡等领域的许多实际问题的数学模型均为半无限规划模型,半无限规划已成为求解实际问题的强有力的工具,关于半无限规划问题的求解方法倍
本文主要研究了一类二阶微分方程组边值问题和一类奇异p-Laplacian方程及n维p-Laplacian方程组边值问题正解的存在性.本文共分为四章:   第一章,简述了问题产生的历史背景
经验似然方法是一类非常重要的构造非参数置信区间和检验的方法,Owen对此方法的一般性质进行了系统的研究.许多研究成果表明,它有类似bootstrap的抽样特性.与传统或者是现代
本文研究了一类非线性Petrovsky方程(组)初边值问题解的定性行为:解的局部存在性、整体存在性、渐近行为和爆破性质。第一章介绍了研究工作的应用背景和发展概况,同时概述了本
半定规划是近20年以来发展起来的一个新的数学规划分支.随着半定规划的发展,实际生活很多领域中的问题可以转化为半定规划模型来求解,例如:经济、金融、工程、管理等.由于它的约
连通图G中的任意两点u和v,一条u-v测地线是指u,v两点间的最短路。令Lu,v表示位于u-v测地线上所有点的集合。对于子集S,令I(S)=Uu,v∈I(u,v),如果I(S)=V(G),则我们称S是G的测地集。把
本文研究涉及带Hardy—Sobolev—Mazya项的奇异半线性椭圆型方程Neumann边界条件下正解的存在性,即方程—△u—λu/|y|2=|u|pt—1u/|y|t有光滑边界的有界区域Ω()RN上Neumann边界条