运用全波离散镜象理论计算二维空域格林函数(二)%3a数值结果.pdfVIP

运用全波离散镜象理论计算 二维空域格林函数(二)：数值结果 周桂蓉冯宁宁方大纲 南京理工大学电光学院毫米波技术研究室，南京，210094 【摘要】：本文以电光调制器的分析为背景，用全波离散镜象(FWDI)的两层 近似方法结合表面波极点贡献的抽取技术计算了一个典型电光调制器两层介质结 构矢量位与标量位空域格林函数，其结果与数值积分的结果吻合很好，而且时间 大大节省。 1． 引言 电光调制器是高速集成光通讯系统中的关键部件。目前．具有SiO：缓冲层的 TI扩散的LiNbO，电光调制器得到了越米越广泛的应用【1】【2】。由于电光调制器 两层介质厚薄悬殊，使得有限差分(FD)等方法的应用不太方便。混合位积分 时，需要重复计算电流源或磁流源激励的矢量位和标量位格林函数【3】，这些格 林函数的谱域解析形式可以用谱域导抗法(SDlA)等方法很方便地得到【4】【5】。 全波离敌复镜象(FWDI)技术的提出使空域计算的慢收敛和快振荡特性带来的 计算困难得以大大缓解【6]。本文在Aksun的两层近似方法的基础上结合表面波 极点抽取技术计算了典型电光调制器结构的矢量位和标量位二维空域格林函数， 其结果与用平均值方法处理的数值积分的结果吻合很好。 2． 有关计算公式 (z，Y，z) 如图l的坐标所示，为一个典型的电 ‘，(，*)＼(z·，，·，z-)光调制器的结构：一层很薄的SiO，缓冲 层和单轴各向异性的LiNbO，介质层。下 ‘：、：占j×一：‘：s姆：：：×·：：、：·： 面以G裂分量为例，说明计算的方法： r}r+r÷r}，+r}f●’÷’+，’’■’■’■’■’■’{1j—d1■1* 在本文的计算中，我们采用如下的 燃鬻黪舞谚蝌 Fourier逆变换形式： 拳麓参麓鲁∥墨享》； } ZⅡ一” 图1一个典型的电光调制器结构：SiO：其中G=是谱域格林函数，它的具体形 缓冲层(厚度为d)以及单轴各向异性 LiNbO，晶体的两层介质结构。 式可参见[5】的推导过程，同样在[5]中还 给出了，它与谱域反射系数的关系。根据如 下的二维Sommerfeld恒等式【81， (2) L筹e-zGo-x,)耻弘弋㈣ 我们可以将G?表示成如p四部分： 1)直射项： 一69一 G∥’,2Uzo (3) Z2C／tf。{2协) 其中，，=、瓜j了i两。 2)表面波项： Ree√咿J’l≈。 (4) Gn“s一_磊Po瓦1 其中t Res，代表谱函数在每个极点的留数和。 (5) 胁，2磊，[筹L。plim吖t，^慨 式中的RⅡ为TE波的谱域反射系数，同时也是被取样拟合的谱函数。在谱函数 中表面波极点部分的贡献为 驴篙等脚∥止小圳 (6) 3)第一层近似的复镜象项： G∥=尝去缸雕协。)