直线感应加速器束流崩溃不稳定性数值模拟.pdfVIP

第 l5卷 第 4期 强 激 光 与 粒 子 束 Vo1．15，NO．4 2003年 4月 HIGH POW ER LASER AND PARTICLE BEAM S Apr．，2003 文章编号： 1001—4322(2003)04—0397—04 直线感应加速器束流崩溃不稳定性数值模拟 张开志 ，林郁正 ， 王华岑 (1．中国工程物理研究院 流体物理研究所 ，四川 绵阳621900， 2．清华大学 加速器实验室，北京 100084) 摘 要 ： 在理论分析 的基础上 ，开发 了直线感应加速器柬流崩溃不稳定性数值模拟程序 。描述 了利用该 程序开展的研究工作．这些研究揭示 了柬流崩溃不稳定性的一般规律，分析了相关参数对柬流崩溃不稳定性的 影响 。最后提出了直线感应加速器束流崩溃不稳定性抑制方法。 关键词： 直线感应加速器 ； 柬流崩溃不稳定性； 横向阻抗 ； 横向尾场； 质心横向位移 ； 数值模拟 中图分类号 ： TL50 文献标识码 ： A 直线感应加速器 中的束流崩溃不稳定性 (BBU)是束流脉冲与感应加速腔相互作用的结果 。束流脉冲经 过加速腔时将在其 中激励起横 向尾场 ，而横 向尾场又将使束流质心产生横 向高频振荡，这两种因素互相耦合， 最终导致束流崩溃不稳定性 。由于束流崩溃的过程非常复杂，采用解析分析方法得出的结论 比较粗略叫 ，必须 采用数值模拟以及试验才能进行 比较全面深入的研究 。 l 数值模拟原理 将束流脉冲等分成 m个束片，研究它们在 个加速腔 中的运动 。束流脉冲经过加速腔时 ，将受到加速 电 场 ，横 向尾场和螺线管线圈磁场 的作用 ，因而能量 、动量和质心的横 向位移都会发生变化 。将加速腔分成漂移 段和加速间隙两段，在漂移段只考虑轴 向磁场 的作用，我们假定在一个加速腔 中，螺线管磁场是均匀分布的。 从文献 [1]中的理论分析出发，可以导出第 i个加速腔漂移段 中第 ．『个束片的质心横 向位移r(f，．)『和横 向动量 ra(，_)『的递推关系 rd(，)一r(i—l，j)cos(k )+ ／(—l，)／[ksin(kf1)] (1a) rd(，．)『一一r(i—l，j)kasin(kf1)+r(—l，_)『cos(kf1) (1b) 式中：kp=qB／2m。 ，q是电荷 ，B是螺线管线圈磁场 ；z是加速腔长度 。 与漂移段相 比，加速间隙很窄，电子经过加速间隙后 ，可以近似认为质心横 向偏移没有变化 ，只是能量和动 量发生 了变化 ，于是 ，加速间隙出口的参量可 以表示为 r(i，．『)一 rd(i，)+ r (2a) rt(，)一 ／d(，．)『+K(i，)+r，，K(i，-)『=(q／m。c。)IW上(￡一￡)I(t)／d(￡)d￡ (2b) (2a)式 中r 及 (2b)式中 ／ 表示两个加速腔轴心之间的偏差 ，代表加速器的准直度 ，在程序 中，它 由一个随机 函数表示 。K(i，．)『为横 向冲量 ；W上(￡)代表横 向尾函数 。模拟计算 K(i，．)『时，将积分近似成求和 一 客exp[‘) ]sin[叫 At]J(d(f△ △￡越小，求和就越接近积分 。电子束能量的变化为 7(0一y(—1)+A7，A7即归一化的能量增益 。模拟计算 时需要 的加速腔参数包括横 向阻抗 、品质因数 、谐振频率和加速腔长度及加速腔数量 ，束流参数包括束流强度 、 平顶、上升前沿 以及注入束质心横 向位移。束流脉冲波形采用 (3)式来模拟时，与实际测量的波形符合得很好 。 J(f)一 J0[1一exp(一n1￡)]{l—exp[a2(￡一t )]) (3) 式中：n一3．5／3，nz一4／3，r代表上升前沿；￡ 为脉冲宽度 。 2 束流崩溃不稳定性发展 的一般规律 假定加速器 由8O个加速腔组成 ，加速腔长度为533mm。取加速腔横