-
本节分别介绍束流生成算法、束流传输模型、电源故障模型以及真空管道孔径模型。
-
生成横向相空间匹配的初始粒子参考了文献[14]中的方法,其中水平方向粒子生成的详细步骤如下。
1) 确定束流初始分布的参数:Twiss参数[4](
$ \alpha_x^{} $ ,$ \beta_x^{} $ ,$ D_x^{} $ ,$ D_x^\prime $ ,$ X_{co.}^{} $ ,$ X_{co.}^\prime $ ),分布类型$ m $ ,粒子个数$ N $ ,束流发射度$ \epsilon $ 及截断倍数$ t $ ,束流动量分散$ \delta $ 等,得到如下参数:2) 生成两组[0, 1]之间均匀分布的随机数
$ S_1^{} $ 和$ S_2^{} $ ,应用如下变化:得到的
$x,x^\prime$ 即为满足条件的初始粒子分布。对于垂直方向而言,其步骤完全一致。值得注意的是,若生成粒子的位置存在色散或非零的闭合轨道,则需要考虑其影响,否则在进行多粒子跟踪时会造成错误的束流包络和束流中心较大的振荡。参数$ m $ 控制着分布类型,其对应关系如表1所列:参数m 束流分布类型 0.0 Hollow shell 0.5 Flat profile 1.0 Uniform 1.5 Elliptical 2.0 Parabolic $\infty$ Gaussian 生成高斯分布粒子分布的算法除可将
$ m $ 设为一个较大的参数外(如$ 1\times10^3 $ ),也可采用著名的“Box-Muller”算法[15],相比其他算法,该算法使用简单,同时具有更好的性能[16]。上述算法也可用以生成一维分布的粒子分布。 -
生成的粒子沿着全环进行逐元件传输,基础元件的传输模型如下所示。
1) 漂移节、四极磁铁和二极磁铁线性元件采用一阶传输矩阵进行束流传输,传输矩阵如下所示:
其中
$ (c_x^{}\,,s_x^{},\,c_y^{},\,s_y^{}) $ 为对希尔方程求解得到的类余弦、正弦项;$ J_1^{} $ 为$ s_x^{} $ 关于元件长度$ L $ 的积分项;$h = 1/\rho$ 为二极磁铁偏转半径的倒数;$ \beta_s^{}、\gamma_s^{} $ 为洛伦兹因子。二极磁铁边缘角对束的流影响可以等效为长度为0的冲击元件,其传输矩阵如下所示:
其中
$ \phi_i^{} $ 是二极铁边缘角;$ \overline{\phi}_i^{} $ 是边缘场的垂直方向修正项[17];g是二极磁铁全高宽;$ I_1^{} $ 是一阶边缘场积分常数,在MAD-X软件中它是我们所熟知的“FINT”值,其由二极磁铁的边缘场分布决定,具体模型这里不再赘述。为了考虑色品效应,对于具有动量偏差的粒子,需要对其磁场聚焦强度进行修正:
其中
$ K_n^{} $ 为n阶磁场聚焦强度;$ \delta $ 为当前粒子和参考粒子的动量偏差,这种修正方法可以实现对具有非零动量偏差粒子的精确传输,其中MAD-X和Sixtrack均采用这种修正。对于二极磁铁仅对边缘聚焦项进行修正而对磁铁主体并未作修正处理,其粒子传输精度为关于动量分散的一阶精度。2) 场误差采用
$ \frac{1}{2}R + {\rm{Kicker}} + \frac{1}{2}R $ 传输模型,其中$ R $ 是漂移节或对应元件的一阶传输矩阵,位于元件中心的Kicker对束流产生的冲击如下:其中
$ L $ 是元件长度;$K_0^{} = 1/\rho$ ;$ K_n^{} $ 为$ n $ 阶磁场聚焦强度;$ \beta_s^{} $ 是洛伦兹因子;Re和Im为对P取其实部或虚部。 -
BRing的二极磁铁为常温磁铁,其磁铁-电源的拓扑结构如图1所示。在正常情况下,磁铁由电源驱动,而一旦电源出现故障,保护开关
$ S $ 会立即导通,储存在磁铁线圈中的能量会通过简单的$ R-L $ 回路放电[13, 18]。磁铁线圈中的电流演化如式(7)所示:图中
$ L $ 为磁铁线圈电感,$ R $ 为磁铁线圈电阻,$ R_{\rm{BD}}^{} $ 为电路的剩余电阻和某些特殊二极磁铁的吸能电阻[2, 19],$\tau = 4L/\left(R_{\rm{BD}}^{}+4R\right)$ 为该放电电路的时间常数。对于四极磁铁和其他常温磁铁,放电规律和二极磁铁类似。 -
图2列出了几种典型的真空室孔径形状,如常见的圆形、跑道形、椭圆形以及在BRing注入、引出段的异形形状等。在程序中每个加速器元件都有相对应的孔径,当粒子跟踪至元件位置处时,通过比较粒子的坐标和孔径尺寸来判断粒子是否超出损失,一旦粒子坐标超出孔径边界,则该粒子标记为“已损失”而在后续的模拟中跳过该粒子而不再进行传输,同时记录该粒子的信息如损失的位置和时刻以便结果分析。
对于磁铁准直误差的影响以及Twiss函数等其他光学参数的计算方法在文献[20-21]中已有详细描述,由于其形式较为冗长,这里不再赘述。
-
在相同的条件下,我们就Twiss参数(包括
$ \beta $ 函数、色散函数)、误差下的闭合轨道、多粒子跟踪统计的束流发射度稳定性以及betatron工作点和色品这几个方面和MAD-X的模拟结果作了检验。Twiss参数如图4(a)所示,我们选取了BRing三个超周期中的一个,在不考虑各种误差的条件下,本程序计算的
$ \beta_x^{} $ 、$ \beta_y^{} $ 以及色散函数$ D_x^{} $ 均与MAD-X得到的结果完全一致。当考虑了BRing二极磁铁和四极磁铁的安装误差后,我们计算了其中一组误差下的闭合轨道,得到了如图4(b)的结果,其水平方向的闭合轨道$ X_{co}^{} $ 和垂直方向的闭合轨道$ Y_{co}^{} $ 和MAD-X的结果相比亦无区别。我们跟踪10 000个初始宏粒子在环内循环10 000圈的情况,束流的初始分布参数如表2所列,在模拟中我们不考虑安装误差和高阶场误差并逐圈记录所有粒子的坐标。图5统计了束流发射度随圈数变化的情况。结果表明,在10 000圈内,水平和垂直方向的束流发射度均保持稳定,没有增长或收缩情况;而对于更长圈数(
$ 1\times10^7 $ 圈)情况亦是如此。这说明本程序可以用来对束流进行长期稳定跟踪。参数 值 束流种类 ${^{238}{\rm{U}}^{35+}}$ 束流能量/(MeV·u−1) 17 束流分布 Gaussian 截断倍数 6 最大发射度/(πmm·mrad) 200/100 统计6倍rms.发射度/(πmm·mrad) 167.2/84.6 基于上述跟踪得到的粒子坐标,我们计算了每个粒子在归一化相空间中的逐圈相移,得到了工作点的小数部分。同时,我们采用多项式拟合的方法得到了统计的工作点和一阶色品,并和MAD-X计算的结果进行了比较(见表3)。如图6所示,除水平方向色品
$ Q_x^\prime $ 外,其余的结果完全一致。$ Q_x^\prime $ 有细微差别的原因是我们跟踪的所用的二极磁铁模型关于动量分散$ \delta p/p $ 是一阶精度的。该差别足够小,对本文要模拟的束流动力学过程几乎没有任何影响。项目 $Q_x^{}$ $Q_y^{}$ $Q_x^\prime$ $Q_y^\prime$ BTracker 0.47 0.43 −11.53 −11.59 MAD-X 9.47 9.43 −11.47 −11.59 -
我们模拟BRing中两台典型的二极磁铁和四极磁铁电源故障后的束流损失,并和MAD-X模拟的结果作比较。图7所示为二极磁铁电源PS41D01和四极磁铁电源PS41QD01故障下的束流损失沿全环的分布。可见PS41D01故障后,束流主要损失在D170漂移节和其附近的四极磁铁QF24上,束流损失相对集中,而PS41QD01这台散焦四极磁铁电源故障后,束流损失相对分散,主要损失在D91和D181附近的漂移节上。MAD-X得到的结果和本程序相比,无论是束流损失的主要位置还是每个位置处的束流占比,均有很好的一致性。
图8所示为上述模拟下的束流损失时间分布。结果可见,二极磁铁故障后,束流损失过程较快,束流在250 μs时刻开始损失,在1.5 ms内完全损失掉,而四极磁铁对应的束流损失相对较慢,大约在15 ms的时间内全部损失。以上结果也表明二极磁铁故障相比四极磁铁故障通常会造成更严重的后果。MAD-X软件得到的结果和本程序的模拟结果在时间分布上也是基本符合的。
-
在相同模拟条件下(参数如表4所列)利用Btracker和MAD-X进行了多粒子跟踪,得到了如图9所示的运行时间结果比较,结果表明本程序相较MAD-X在计算性能上有明显的提升。
参数 值 CPU Intel Core i9-9900K@3.6 GHz 平台 Windows 10 64 Bit 宏粒子数 10 000 跟踪圈数 10 000 通过外部脚本逐圈修改MAD-X软件的输入文件以改变磁铁的磁场强度,可以实现电源故障下的束流损失模拟,但其存在着使用繁琐和计算速度慢的问题。MAD-X采用Fortran和C++语言编写,由于设计之初并未针对多粒子模拟,相关算法并未做相关优化,并且只能使用单个线程,因此对于多粒子跟踪这种计算密集型任务而言,其速度较慢。对于BTracker,由于其粒子跟踪模块基于C++语言编写,跟踪算法也作了大量的优化以保证基础性能,同时利用CPU缓存和OpenMP并行技术进一步提高运算能力,可以充分、高效地利用计算机资源以实现高性能计算,在速度上较MAD-X有将近80倍的提升。利用该软件我们可以非常迅速地在分钟量级甚至秒量级时间完成模拟,得出结果;而使用MAD-X软件,同样的单次模拟则会需要若干小时的时间。基于BTracker计算性能高的特点,还可以通过多粒子跟踪进行加速器参数优化,这类优化过程通常需要进行非常可观的模拟次数,若使用MAD-X软件,其所需要的耗时将是非常冗长的。在HIAF机器保护模拟研究中本程序已得到应用[19]。
The Development of BTracker, a Simulation Code for Beam Losses From Magnet Power Failure
doi: 10.11804/NuclPhysRev.40.2022100
- Received Date: 2022-09-13
- Rev Recd Date: 2022-10-25
- Available Online: 2024-02-04
- Publish Date: 2023-12-20
Abstract: The booster ring(BRing) of High Intensity heavy-ion Accelerator Facility(HIAF) is characterized by high intensity and high energy. Failures of key equipment such as magnets and RF system will lead to significant beam losses and damage to the accelerator. A beam dynamics simulation code BTracker is developed to study beam losses under equipment failures and beam dynamics during beam dumping for the HIAF machine protrection. With this code, beam loss distribution under power failure of dipole and quadrupole magnets is simulated, and the results are benchmarked with MAD-X code. In addition, BTracker provides user-friendliness, interface flexibility and high performance to meet the requirements of machine protection simulation for HIAF-BRing.
Citation: | Hang REN, Jiancheng YANG, Geng WANG, Guodong SHEN, Liping YAO, Jie LIU, Shuang RUAN, Lina SHENG. The Development of BTracker, a Simulation Code for Beam Losses From Magnet Power Failure[J]. Nuclear Physics Review, 2023, 40(4): 547-554. doi: 10.11804/NuclPhysRev.40.2022100 |