课程设计---锥形过渡波导的损耗特性.doc

课程设计锥形过渡波导的损耗特性所属课程：集成光学姓 名： 学 号： 班 级： 指导教师： 摘 要在集成光电子器件中，通常需要把两彼此分离的两波导结构连接在一起。此两分离的波导结构在结构参数上可能是相同的，也可能是不同的。同弯曲波导一样，锥形波导也是一种动态的波导器件，不能对其进行本征分析。近些年，已提出多种数值方法来优化设计过渡波导，其中光束传播法是人们最常用的方法。研究表明，光束传播法非常适宜用来研究象锥形过渡波导这样的沿光传播方向变化的动态器件。对wiwo的锥形过渡波导,利用有限差分光束传播法对几种形状侧面边界下的过渡波导进行了分析。本文就是基于该方法主要研究介绍指数型锥形过渡波导（凹形）在两种不同的折射率情况下:1) 窄端口为输入端，研究输出波导输出功率随宽端口宽度的变化规律；2） 宽端口为输入端，研究输出波导输出功率随宽端口宽度的变化规律；关键词：光束传播法（bpm）；指数型锥形过渡波导；波导输出功率i山东理工大学课程设计说明书 目录 集成光学目 录摘 要i目 录ii第一章 绪 论11.1 集成光学的发展11.2 optibpm软件简介11.3 本课程设计的目的2第二章 bpm基本原理及常用的边界条件22.1 bpm基本原理22.2 bpm中常用的边界条件3第三章 输入端口相对宽度对锥形过渡波导传输特性的影响53.1 几种常用形状光波导介绍63.2仿真软件应用简介93.3仿真结果12致 谢17参 考 文 献18ii山东理工大学课程设计说明书 正文 集成光学第一章 绪 论1.1 集成光学的发展集成光学的概念，起源于二十世纪60年代末、70年代初，其主要思想是在共同衬底上建立各种光学器件，然后用薄膜波导将他们连接起来，从而形成一个能完成特定功能的功能芯片。然而由于当时各种器件对材料、工艺、结构的要求差异很大，技术要求十分复杂，实现所有器件的全面集成仍存在困难，集成光学的目标因而发生了一定的变化，从全光集成到可集成部分电路的部分光集成，从注重集成光学技术到研究探索各种光波导光学器件。目前集成光学研究热点主要集中在理论与器件研究。理论研究热点主要集中在以下两个方面：（1）围绕新型集成光学器件的结构设计、功能模拟与特性、参数的计算等；（2）基础理论研究。集成光学基础理论研究主要集中在两类：一是基于固体物理学的基本理论和方法，研究和探讨制作微观集成光学器件的可能性。二是基于波动光学的基本理论和方法，从波导光学的角度研究集成光学器件。集成光学器件的研究主要是：在光通信领域的集成光学器件研究、集成光学传感器和其他集成光学器件，如mens光开关阵列等。1.2 optibpm软件简介optibpm 是一套功能强大、使用者界面友善且可利用计算机辅助设计的设计仿真软件，并可设计及解决不同的积体及光纤导波问题。光束传播法，或称为 bpm是optibpm的核心，是一种一步接着一步来仿真光通过任何波导物质的行为，在积体光学及光纤光学中，当光传播经过一可传导的结构时，其光场可以在任一点被追踪出来，bpm可以允许观察任一点被仿真出的光场分布，而且可以容许同时检查辐射光及被传播的光场。 光学波导是光组件中的重要组件，它可以在光讯号中扮演传导、耦合、开关、分光、多任务及解多任务的角色，被动波导、电光组件、发射器、接收器及电子部分装置被集成于一个芯片上，使用的技术为平面技术，其就好像微电子的技术。 optibpm是一套使用者界面非常友善的软件，它可以在二维及三维的波导组件上仿真光的传播，且optibpm三维仿真提供了任何所需要的步阶折射率(step index)的波导设计。应用范围：(1): y型波导、波导分波装置(2): 波导耦合器、各类波导元件(3): 非线性波导(4): 光纤设计(5): awg设计(6): 混合光源的运算1.3 本课程设计的目的1）了解光束传播法（bpm）的原理；2）了解指数型锥形过渡波导及其在集成光路的应用；3）掌握并熟练应用波导分析模拟分析软件optibpm；4）窄端口为输入端，研究输出波导输出功率随宽端口宽度的变化规律；宽端口为输入端，研究输出波导输出功率随宽端口宽度的变化规律；第二章 bpm基本原理及常用的边界条件2.1 bpm基本原理bpm理论来源于波动方程，波动方程是建立在helmhotze方程基础上的。helmhotze方程的一般形式为： （2-1）式中，波导的几何结构完全由折射率分布来确定。考虑到在典型的波导问题中，场的迅变部分是沿波导轴向传输引起的位相变化，假设波导轴只要是沿z向，那么，可以引入一个所谓的慢变场，设： （2-2）这里是参考波数，用于表示场的平均相位变化。参考波数通常通过以参考折射率的形式表示。是任意常数，要求它的选择使是z的慢变函数。把（2-2）式代入到helmholtz方程便可得到慢变场所满足的方程： （2-3）方程（2-3）同确切的helmhotze方程是完全等效的。若随z的变化足够慢，方程（2-3）的第一项同第二项相比就可以忽略不计，这就是大家熟悉的慢变包络近似，也称为傍轴或抛物近似。经运算，可得到： （2-4）这是基本的三维标量形式的bpm方程，若忽略与y有关的项就可得到二维标量形式的bpm方程。给定一个输入场，上述方程决定了在空间内的场分布。对许多具体的问题，迅变位相因子的引入，可使得在数值计算中剖分网格可在纵向（即z方向）大于光波长，这很大程度上提高了数值技术的效率。此外，和z有关的二次微分项的忽略，使二阶边值问题，转变为一阶初始值问题。这一点也同样提高了bpm方法的计算效率。但慢变包络近似使我们只能考虑波导中沿z轴附近的传播场，这也对参考折射率的选择有了限制。而对于如多模干涉波导器件中所存在的有复杂的位相变化的场，利用该式就不能进行精确地模拟。2.2 bpm中常用的边界条件边界的处理是人们应用bpm时所必须面对的一个重要的问题，不好的边界条件，会使入射到边界处的光反射回计算窗口，引起计算误差，因此选择合适的边界条件是很关键的。bpm中常用边界条件除了传统的dirichlet条件和neumann条件外，最常用的还有吸收边界条件（absorbing boundary conditions，abc）、透明边界条件（transparent boundary conditions，tbc）、完美匹配层条件（perfectly matched layer，pml）和补偿算子条件（complementary operators method，com）等。吸收边界条件是人为在计算区域外加一层损耗层，通过选择合适的吸收系数和吸收层厚度使出射光波在到达吸收层边界时衰减至零。最早的傅立叶变换bpm 方法对吸收边界处理是在边界处乘上一个合适的复介电常数, 以吸收边界的反射，较早的差分bpm 采用的也是类似的方法，需要选择合适的吸收常数是这种边界条件的关键。这种方法的难点在于如何选取合适的吸收系数和吸收层厚度，并且也没有一个普遍的选取方法。为了解决这些问题，人们又提出了许多新的吸收边界条件。这里主要介绍透明边界条件和完美匹配层条件。一、tbc边界条件hadely首先提出了透明边界条件(transparent boundary condition ，tbc)，它是在边界附近把光场近似看成是平面波，并以平面波的形式在边界处向外透射出去。假定光在z向传播过程中，有一沿x向右向无穷远处辐射的辐射模，其中e0、kx是复常数，e的实部表示边界处相位的变化，虚部则表示场幅度的微分。只要kx的实部是正数，就能够保证辐射光能量只会流出计算窗口，因此tbc算法的关键在于kx的优化选取。以下只考虑二维情况下的与z向平行的右边界，利用上述假定，当光沿z向传播了的距离时，光能流相应地向右流动了x的距离，对于差分格点m和m-1(临近计算窗口右侧)，此时有关系式： （2-5）当第n步计算完成后，可以得到n+1步计算所要的kx （2-6）但是在应用上式之前，必须保证kx的实部是负的以确保只有向外的辐射而没有反射。以上是简单tbc，完全tbc的算法相对要复杂些。kx的是通过下式来确定的： （2-7）上式中m是边界附近某一格点，m的值为 ： （2-8）式中z，x是计算区域的离散间隔，k=k0nr，是真空中的波数，nr为参考折射率，kx是上一步计算中的近似平面波矢的横向分量。当mwo的锥形过渡波导,利用有限差分光束传播法对几种形状侧面边界下的过渡波导进行了分析，并在入口宽度固定的情况下对出口宽度进行了优化。但在许多波导器件中，也会出现入口宽度大于出口宽度的情况，即wiwo时，对凸起型指数锥形过渡波导g0，对凹型指数锥形过渡波导g0。当 wiwo时，情况则相反。图3.3给出了窄口输入时指数型过渡波导在g取不同值时的形状。对凹型波导，形状因子g越大,过渡波导在入口端宽度的变化越缓慢，但在出口处波导宽度的变化却较大。对凸型波导则相反。宽口输入时指数型过渡波导的形状随形状因子的变化与此相类似。 （a） （b）图3.3 形状因子g对指数型锥形过渡波导形状的影响(a) 凸型 (b) 凹型 3.2仿真软件应用简介下面计算机模拟过程中软件的截图，以此说明软件的基本用法。图3.4软件界面图3.5 软件界面及波导端口宽度的定义 图3.6 仿真操作界面1图3.7 仿真操作界面2图3.8仿真操作界面3图3.9数据生成及保存界面图3.10 指数型光波导的光场分布图3.11 指数型光波导的折射 率分布3.3仿真结果（一）包层折射率为1.442情况下:(1) 表3-12为在窄端口为输入端，输出波导输出功率随宽端口宽度的变化规律宽度（m）输出功率d=4d=5d=6d=7d=8d=9d=10d=11d=120.999973121291 0.9996349464420.9994618930970.9993686769970.9990904632080.9989548903970.9986870255930.9983636747340.997532558282 表3-12下图是相对应的输出功率随输出端口宽度变化波形图（x轴为输出端口宽度，y轴为输出功率）图3-12由图可知, 输出功率随输出端口宽度增宽而减小。(2) 表3-13为宽端口为输入端，输出波导输出功率随宽端口宽度的变化规律宽度（m）输出功率d=4d=5d=6d=7d=8d=9d=10d=11d=120.999550011282 0.9992783377990.998957189293 0.998943541341 0.9991547389690.998974212158 0.9990653027350.998725072534 0.999392470029 表3-13下图是相应的输出功率随输入端口宽度变化波形图（x轴为输入端口宽度，y轴为输出功率）图3.13由图可知,输出功率随输入端口增大而波动变化,先减小再增大再减小。（二）包层折射率为1.300情况下:(1)表3-14为在窄端口为输入端，输出波导输出功率随宽端口宽度的变化规律宽度（m）输出功率d=4d=5d=6d=7d=8d=9d=10d=11d=12d=13d=140.999194334834 0.9992804246470.9991654048470.9994880305410.999239390010.9992509951280.9992259478650.9992686191250.9992826039490.9993246665620.9993571102表 3-14下图是相应的输出功率随输入端口宽度变化波形图：（x轴为输入端口宽度，y轴为输出功率）图3.14由图可得:输出功率随输出端口增宽先减小再增大,再减小的波动变化.(2) 表3-15为宽端口为输入端，输出波导输出功率随宽端口宽度的变化规律宽度（m）输出功率d=4d=5d=6d=7d=8d=9d=10d=11d=12d=13d=140.999194334834 0.9996635107110.9986908489790.9842964565520.9720702386090.9994881295680.9995312057170.9994549551990.9992826039490.9996999179660.850665369741表3-15下图是相应的输出功率随输入端口宽度变化波形图：（x轴为输入端口宽度，y轴为输出功率）图3.15由图可得:输出功率随输入端口增大先围绕一个值缓慢波动,然后到达某一特殊值时下降。16山东理工大学课程设计说明书 致谢 集成光学致 谢本次课程设计是在我的导师武继江老师的悉心指导下完成的，武老师严谨的治学态度给了我极大的帮助和影响。感谢武老师在这一期间对我的关心和指导。武老师不仅在课程设计期间给我悉心的指导，在学习和生活上都给予我很大的关心和帮助，在此我再次向武老师表示衷心的感谢!17 参 考 文 献1 shi y c, dai d x, he s l. improved performance of a silicon-on-insulator-based multimode interference coupler by using taper structures. opt. commun.,2005,253(4-6): 276-2822 mizuno t, kitoh t, itoh m, et al. optical spot size converter using narrow laterally tapered waveguide for planar lightwave circuits. j. lightwave technol., 2004, 22 (3):833-8393 hettrick s j, wang j, li c, et al. an experimental com