粒子模拟中波导激励源的设计与实现.pdfVIP

! 第# 卷! 第$ 期 强 激 光 与 粒 子 束 %’( # ，)( $! ! $**+ 年$ 月 ,-., /0123 45623 5)7 /538-942 :25;6 7=( ，$**+! 文章编号：! **?@$$（$**+ ）$$*A*+ 粒子模拟中波导激励源的设计与实现! 周! 俊，! 祝大军，! 刘大刚，! 刘盛纲 （电子科技大学 物理电子学院，成都+**B? ） ! ! 摘! 要：! 讨论了粒子模拟中用时域有限差分（C787 ）法设计波导激励源的常见问题，并根据波的传播特 性，提出一种模拟波导激励源的新方法，给出通用的计算公式，讨论了几种特殊情况下的算法修正。给出了$ 维柱坐标系下磁绝缘线振荡器（;-40 ）计算实例，结果表明：该方法能模拟较为复杂的波导激励源，并有效消除 了激励源边界处的虚假反射，与传统波导激励源模拟方法比较，验证了该方法的正确性与优越性。 ! ! 关键词：! C787 法；! 粒子模拟；! 电磁场；! 激励源 ! ! 中图分类号：! 0?+$ ( ? ；8)*B! ! ! ! 文献标识码：! 5 ［］ ! ! 时域有限差分（C787 ）法 近年来被广泛用于各类电磁场问题的计算，由于C787 采用时间和空间的差分 对;DEF’’ 方程组直接离散化，因此在计算复杂问题时具有非常大的灵活性。C787 用于模拟波导激励源时， 要尽可能保证源的特性与实际物理模型性质一致，这就要求必须合理考虑激励源边界处的散射场，所以目前除 了一些简单的应用采用强迫激励源外，大部分激励源引入方式都是建立在总场G 散射场体系上的，这类方法通 常将计算区域划分为总场区和散射场区，分区计算之后再在两区的连接边界上引入连接条件，从而严格计算出 场分布，该过程不含任何近似处理，形式紧凑，计算精确，但需专门划分出散射场区用以辅助计算总场，并且在 迭代过程中要对连接边界两侧的场信息进行多次传递，处理较为繁琐。 ! ! 本文提出的方法也建立在总场G 散射场体系上，但并没有对计算区域进行划分，而是从波的传播特性出发， 利用波传播方向上相邻网格点上场的权重关系，简便地计算出散射场，从而与已知入射场共同求得总场。文中 给出详细的推导过程和必要的算法分析，并将该方法设计的波导激励源引入到自主开发的粒子模拟软件9,- /-9 中。计算了$ 维柱坐标系下的;-40 器件，结果表明该方法能模拟较为复杂的波导激励源，并有效消除激 励源边界处的虚假反射，与强迫激励源相比，各方面都更为优越；与传统总场G 散射场体系波导激励源模拟方法 相比，计算精度相当，但更为简便灵活，并能节省内存使用和计算时间。 ! 波导激励源的传统设置方法 !( ! 强迫激励源 ! ! 通过直接对特定的电场或磁场分量强行赋予所需的时间变化形式，可以简便地建立起强迫激励源。强迫 激励源形式简单、使用方便，可用来处理许多工程问题；但它没有考虑散射场，容易在激励源边界处引起虚假反 射，造成数据误差。例如，常用的最简单的高斯脉冲激励源和正弦波激励源可表示为 $ ! （ ）# ! EI ｛$ ［（ $ ）% ］｝ （ ） H * * ! （ ）# ! JHK（ ） （$ ） H * ! ! （ ）表示引入的电场激励源形式，是时间 的函数；! 表示电场的最大值。对于式（ ）所示的高斯脉冲 式中：H * 源， 表示脉冲的中心点时间； 是一个基准时间，场值在 L M