PIC方法在RFQ加速器z—code和t-code模拟程序中的应用研究.pdfVIP

维普资讯 第28卷 第5期 核 技 术 Vo1．28，No．5 2005年 5月 NUCLEARTECHNIQUES Mav2005 PIC方法在 I Q加速器z-code和t-code 模拟程序中的应用研究 邢庆子 林郁正 傅世年 1(清华大学工程物理系 北京 100084) 2(中国科学院高能物理研究所 北京 100049) 摘要 本文介绍了二维PIC(Particleincel1)方法，这种方法常用于粒子动力学模拟中空间电荷作用的计算。 并比较了以时间为自变量 (t-code)和以纵向位置为 自变量 (z．code)的两种动力学模拟程序；针对 “国家重点 基础研究发展规划”洁净核能项目中的射频四极 (RFQ)加速器结构参数，给出了单束加速和正、负离子束同 时加速两种情况下，t-code和z．code模拟得出的传输效率。结果表明，当束团的相位宽度大或能散大时，z．code 在计算空间电荷作用时会引入相对较大的误差，从而应该使用 t-code来进行动力学模拟，以获得更准确的结 果。 关键词 PIC(Particleincel1)方法，空间电荷作用，t-code，z．code，射频四极 (RFQ)加速器 中图分类号 TL53 强流是粒子加速器的一个重要发展方向。随着 目成正比I】。因此，国际上许多粒子模拟程序，如 束流流强的增大，空间电荷效应也随之增强，从而 PARMILA 、 PARTRAN 、 LrNAC 、 IMPACT 、 导致束流发射度的增长和束晕的形成。正、负离子 PARMELA等，均使用 PIC方法计算空间电荷作 束同时加速，是人们为了抑制质子直线加速器低能 用I。 端的空间电荷效应，同时提高加速器的流强而提出 PIC方法包括有2D和3D计算方法。在上面提 的一种新方案I1J。从离子源引出的等流强的正、负 到的5个程序中，除PARMILA能够使用2D和3D 离子连续束，经混束磁铁混合成为中性束，经匹配 计算外，其余程序均使用 3DPIC方法计算。国际 后，进入RFQ(Radiorfequencyquadrupole)加速器。 通用的 RFQ设计和动力学模拟程序 PARMTEQJ 正、负离子束团在RFQ的成形段中逐渐分离，在纵 使用2DPIC方法计算，其中用到K．R．Crandall编 向相距 180。高频相位，然后同时被加速和聚束。为 写的空间电荷计算子程序 SCHEFF(SpaceCHarge 了研究清楚这个物理过程，需要进行必要的粒子动 EFFect)，它的局限性是束团在横向必须为角向对 力学模拟，其中空间电荷作用的计算是非常重要的 称分布 (或椭圆对称分布 ) J。3DPIC子程序有 一 部分。计算空间电荷作用比较常用的方法有 PICNIC(ParticleInCellNumericalIntegrationbe- PPI(Particletoparticleinteraction)和PIC(Particlein tweenCube)I9】，它的基本方法与SCHEFF相同，但 cel1)两种。PPI方法中，为避免两个宏粒子距离过 对束团分布没有任何要求。文献 1【0】中使用3DPPI 近从而产生不切实际的很大的库仑排斥力，一般将 方法改写了PARMTEQ程序，但完全的3DPIC方 每个宏粒子近似为具有一定半径，均匀带电的球体， 法还没有应用到PARMTEQ中。尽管3D模拟程序 计算每个宏粒子受到的空间电荷力，必须要求出其 更精确，但它需要的计算时间长，在加速结构的设 它所有宏粒子对其产生的库仑力，因此PPI方法的