《2.5维圆波导开放边界截断的修正完全匹配层方法》.pdfVIP

! 第# 卷! 第$ 期 强 激 光 与 粒 子 束 %’( # ，)( $! ! *++, 年$ 月 -./- 01234 56734 6)8 0649.:53 ;367 =( ，*++,! 文章编号：! ++?@$**（*++, ）+$?+@A,?+B ! # 维圆波导开放边界截断的修正完全匹配层方法! 李小泽，! 王建国，! 王! 玥 （西北核技术研究所，西安, ++*@ ） ! ! 摘! 要：! 在时域有限差分算法中，完全匹配层是一种高效的吸收边界条件。通过在标准完全匹配层 （05 ）方程中添加具有物理意义的电导率和导磁率附加项，提出了修正完全匹配层（05 ）方法，给出了圆柱 坐标下$ 维05 方程，并应用于* ( B 维圆波导开放边界的截断计算。计算结果表明，在不同频率和模式的激 励下，加入平滑电导率和导磁率附加项后，与标准05 方法相比，05 方法的相对误差降低C D + E; 。 ! ! 关键词：! 修正完全匹配层；! 圆波导；! 开放边界；! 截断；! 时域有限差分 ! ! 中图分类号：! 9)C** ( C ；9)C$! ! ! ! 文献标识码：! 6 ! ! 利用时域有限差分方法（F898 ）模拟开放区域的电磁问题时，吸收边界是必需的。;GGHIG 于##@ 年首 ［］ 先提出完全匹配层（05 ）的概念 ，这是一种建立在分裂场基础上的媒质吸收边界条件。理论上可以设置 05 吸收媒质，使其对于任意角度、任意频率的入射波反射无限小。由于数值离散的原因，05 仍会产生一定 的数值反射。文献［*?@ ］中报道了05 产生数值反射的原因。 ! ! 许多高功率微波（-0 ）器件，如返波管、行波管、虚阴极振荡器等，采用轴对称的圆柱波导结构，所有物理 量关于角向 对称。在这些-0 器件的模拟中，大部分采用* ( B 维模拟，即空间* 维，电场、磁场及带电粒子 ! 速度是$ 维分布。文献［@ ?B ］分析了通用粒子模拟软件中吸收边界的种类，比较了吸收边界的优缺点，并提出 用卷积完全匹配层（:05 ）截断直角坐标系中$ 维开放波导和圆柱坐标系* ( B 维圆波导开放边界的新方法， 模拟结果表明这种:05 对色散电磁波具有很好的吸收效果，但在实现粒子通过截断面时很复杂。文献［, ］ 采用分裂形式的05 研究直角坐标系下粒子模拟开放边界的截断问题，提出了非对称完全匹配层（605 ）方 法。605 可以考虑粒子通过截断面的情形，但有许多缺点：（ ）形式复杂，难以实现；（* ）中间存储量多，所需 ［@ ］ 内存量较大；（$ ）满足弱因果关系，在长时间的模拟中会出现后期反射和不稳定性 ；（@ ）吸收凋落波的性能较 差。本文在综合605 和:05 优点的基础上，提出了修正完全匹配层（05 ）方法，给出圆柱坐标系下的 05 方程，并计算了* ( B 维圆柱波导边界截断的吸收效果。 $% 修正’( 公式 ! ! 为便于推导，对=JKG’’ 方程做如下变换 $ $ $ ! # ! # ，! ! # ! ，! # （ ） * 式中：# 为真空中光速；L （ ） ，为介质波阻抗；和 分别为介电常数和磁导率。通过上述变换，得到 # $ $ # %! %