感应加速腔横向阻抗模拟计算.pdfVIP

维普资讯 第 12卷 第 6期 强 激 光 与 粒 子 束 Vo1．12，No．6 2000年 11月 HIGH POW ER IASER AND PARTICIE BEAM S Nov．，2000 文章编号 ： 1001--4322(2OOO)O6一O763一O4 感应加速腔横 向阻抗模拟计算 ‘ 张开志h ， 林郁正 (1．中物 院流体物理研究所 ，绵阳 919信箱 lO6分箱 ，621900； 2．清华大学加速器实验室 ，l00084) 摘 要 ： 束腔相互作用导致的束流崩溃不稳定性 (BBU)将严重地恶化束品质 ，而加速腔横 向阻抗 反映了这种相互作用的强弱。用MAFIA程序对三种不同的感应加速的腔横 向阻抗进行 了模拟计算 。计 算结果与采用脉 冲双线法测量 的结果反映了相 同的物理性质 。 关键词 ： 感应加速腔 ； 束流崩溃不稳定性 ； 横 向阻抗 ； MAFIA 中图分类号 ： TL50 文献标识码 ： A 直线感应加速器 中能传输 的最大束流受到束流崩溃不稳定性 (BBU)的限制。束流在非理想导电或 不连续 的束管道 (如含加速腔的束管道)中运动时，会在其 中激励起尾场 ；而尾场又将反作用于束流 。在 一 定的条件下 ，横 向尾场与束流 的相互作用将形成正反馈 ，导致束流崩溃不稳定性 。横 向阻抗反映了这 种相互作用 的大小 ，是研究束流崩溃不稳定性 的重要参量 。因此 ，加速腔的横 向阻抗模拟计算是研究束 流崩溃不稳定性必不可少的过程 。 l 计算原理 横 向阻抗模拟计算 ，即求解给定束流在加速器组件中激励的电磁场这样 的电动力学 问题 ，可 以在频 域进行 ，也可 以在时域进行 。 MAFIA 的E模块可用于在频域求解阻抗 。具体应用到直线感应加速腔时 ，却存在诸多的局 限。首 先 ，高阶模计算量大 ，而且结果也不准确 。其次 ，感应加速腔 中的铁氧体是一种色散材料 ，其介 电常数和 磁导率 随频率变化 ，但 E模块却几乎无法处理色散材料 。 MAFIA 的T3模块可在时域计算横 向尾场 ，通过 Fourier变换 ，由尾 函数得到横向阻抗 。更重要 的 是 ，利用 T3模块可以处理色散材料 。但是 ，T3模块不适于计算长程尾场 。不过 ，感应加速腔 Q值很低 ， 导致 BBU 的是短程尾场 。于是最终选用 了T3模块计算横 向阻抗 。 MAFIA在 时域计算横 向尾场的理论基础是有 限积分原理 (FIT)，即在积 分形式下求解 Maxwell 方程 。分别用场 的平均值乘 以路径或 面积 ，取代 安培定律和法拉第定律 中的线或面积分，就得 到 Maxwell方程 的一级近似 。应该注意的是 ，电场定义于矩形网格边的中点，而磁场定义于 网格面 的中心 ， 将 电场和磁场连接起来就形成两个对偶 网格口]。 在 FIT方法 中，法拉第定律表示成积分形式 ，CDse=D b，其 中e和 b都是 3N 维的矢量 ，N是 网格 的节点数 。c是一个 3N×3N 的矩阵，对应于求旋度 。 是 3N×3N 的对角矩阵，包含 网格单元边 的长 度 ，D 也是 3N×3N 的矩阵，包含 网格单元的表面积 。同样 ，安培定律表示为 硒 ^ 一 + )，上标 “ 一 ”对应对偶网格 。考虑到材料分布时 ，MAFIA采用 了关于时间变量 的蛙跳技术 ，用递归方法来计算 时变 电磁场 b 一 b ——zSt(02CD )e 。 (1) e 。 一 e 一 ~ (D-／D2 _|sD )6 ——zStD~-J+ (2) 其 中上标代表时间步长 (=nAt)，D。和 是描述 网格 中填充介质的参数 e和 的对角矩阵。 ·收稿 日期 ：2000—01—25； 修订 日期 ；2000—04—2O 基金项 目：国防科技基础研 究基金 资助课题 作者简介 ；张开志 (19