高频腔本征值问题的三维数值模拟Ξ-硬件和射频工程师.PDFVIP

　第 15 卷 　第 2 期 强 激 光 与 粒 子 束 Vol. 15 ,No. 2 　 　2003 年 2 月 HIGH POWER LASER AND PARTICL E BEAMS Feb. ,2003 　 文章编号 : (2003) 高频腔本征值问题的三维数值模拟 雷文强, 　杨中海 , 　肖 　礼 ( 电子科技大学 物理电子学院, 四川 成都 610054) 　　摘 　要 : 　运用有限积分法软件 MAFIA 和有限元软件 ANSYS 对两类高频腔的本征值问题进行数值模 拟 ,将模拟结果与部分实验数据进行验证 ,得到了比较一致的结果 ,其平均计算精度在 5 %以内,证实了两种电 磁软件都可用于高频腔本征值问题的数值模拟 ,在计算精度上为满足工程设计需要做了一些研究。另外将两 种模拟软件在计算时间上作了对比。 　　关键词 : 　有限积分法 ; 　MAFIA ; 　有限元法 ; 　ANSYS ; 　数值模拟 　　中图分类号 : 　TN124 　　　　文献标识码 : 　A 　　高频谐振腔在微波电子学以及加速器领域得到了广泛的运用 ,谐振频率、品质因数是高频谐振腔的重要参 ( ) 数。对于形状复杂、填充介质的高频腔难以得到解析解 ,通常由数值模拟方法求解 ,有限元法 FEM 是其中常 用的计算方法之一 ,它将场矢量离散处理并求解波动方程或泊松方程 ,在求解形状复杂或非均匀媒质问题时可 满足计算精度 , 成为 ANSYS 软件模拟计算的理论基础。然而在场矢量不连续的媒质分界面处无法满足 [1 ] ( ) [2 ,3 ] Maxwell 微分方程 ,而 Maxwell 积分方程却在函数不连续的情况下仍然适用 。因此有限积分法 FIT 由 Maxwell 积分方程出发 , 以独特的网格单元形式直接将场矢量离散在三维网格中并进行差分求解,成为 MAFI A 软件模拟计算的理论基础。 　　本文叙述了两种方法求解本征值问题的理论公式 ,并运用 MAFIA 和 ANSYS 软件计算高频加速腔[4 ]和休 斯行波管耦合腔[5 ] 的本征频率和 Q 值 ,将模拟结果与部分实验数据进行对比 ,得到了基本一致的结果 ;并对两 种软件在计算时间上的优劣进行了比较。 1 　有限积分法 　　有限积分法将实体模型分为若干小立方体的三维网格单元 G , 电场矢量离散分配到 G 单元 ,称之为电场 单元。电场分布在单元 G 的每条边的中点处 ,磁场分布在相应网格单元 G 的每个面的中点处 ,各个磁场方向 之间构成与 G 相互正交的另一个立方体网格单元 G ,称之为磁场单元。图 1 表示电场和磁场的网格单元 G , ε μ κ G ,其形式类似于 Yee 氏网格单元。同样, 电流密度 J 和各个媒质材料的介电常数 , 磁导率 , 电导率 都可 如此分配到相应的网格单元中。各个相关物理量在离散分配后形成若干矩阵向量 e , d , b , h , j , Dε, Dμ, Dκ,其 εμκ 中D , D , D 分别是离散化的 , , 矩阵向量 ; d , h分别为离散化的电位移矢量D 和磁场强度H 的矩阵向量。e = ε μ κ Fig.