Design and implementation of vacuum control system of China Spallation Neutron Source
-
摘要: 中国散裂中子源(CSNS)加速器真空控制系统负责真空数据采集、设备监控和闸板阀控制与联锁,是设备运行和故障诊断以及超高真空保持的重要保障。本文介绍了加速器真空需求,基于实验物理及工业控制系统EPICS软件框架的真空控制系统设计与实现,使用横河可编程逻辑控制器PLC控制与联锁设备,摩莎MOXA工控机监测真空状态,EPICS PV数据直接进入声音报警系统和历史数据库系统,为工作人员及时发现和处理问题、进行后续数据分析和机器研究等提供了便捷途径和可靠保障。目前,该系统已完成现场安装和调试,并已正式投入运行。运行结果表明,该系统具有稳定性好、可靠性高、人机交互友好的特点,很好地满足了加速器真空控制系统运行的需要。
-
关键词:
- 中国散裂中子源 /
- 超高真空 /
- 真空控制系统 /
- 实验物理及工业控制系统
Abstract: The accelerator vacuum control system of China Spallation Neutron Source (CSNS) is designed to be responsible for vacuum data acquisition, equipment monitoring and gate valve control and interlock, which is an important guarantee for safe operation, fault diagnosis and ultra-high vacuum. This paper introduces the vacuum requirements of the accelerator, and the vacuum control system design and implementation. Based on the experimental physics and industrial control system (EPICS) software framework, the Yokogawa PLC is used to control and interlock equipment, the MOXA IPC monitors the vacuum status, the EPICS PV data directly enters the BEAST alarm server and the channel archiver, providing convenience and reliable guarantee for the staff to discover and solve problems in time, as well as carry out subsequent data analysis and machine research. At present, the vacuum control system has been installed and debugged on site, and has been put into operation. The results show that the system has the characteristics of high stability, high reliability and friendly operator interface, which satisfies the operational needs. -
激光等离子体加速器可以达到TV/m的加速梯度,在极短的距离便可将离子加速到较高能量[1-2],在医用等方向具有重大应用潜力[3],相关应用要求我们必须精准测量束流横向位置。激光驱动质子束流具有短脉冲、低电量等性质,常见于质子装置的BPM,如斜切型BPM[4],shoebox BPM[5]难以满足要求,而电子装置中的腔式BPM[6]具有精确测量低电量束流位置的能力,因此我们探索在激光驱动质子装置中使用腔式BPM的可行性,以期望其实现束流位置的精确测量并伴随激光等离子加速器广泛应用到医疗等领域。本文根据激光驱动质子束流的特点,推导了其通过腔式BPM的输出信号,基于该原理进行了腔式BPM的设计和优化,并讨论了该BPM对于不同激光驱动质子束流的适用性。
1. 腔式BPM对激光驱动质子束流的响应
1.1 腔式BPM的原理
腔式BPM是利用束流经过谐振腔后激励起的偶极模来监测束流的横向位置偏移。以矩形腔的偶极模之一TM210模为例,其电场形式为[7]
{Ex=Ey=0Ez=E0sin(2πax)sin(πby)eiωt (1) 式中: E0为常数,a、b分别为矩形腔的长、短边。考虑一束特征脉宽长度为σt的高斯分布的粒子束,从横向位置(x, y)入射谐振腔,应用近轴条件x=a/2+x′(x′≈0)、y=b/2,有
{RQ0=|V|2ωU≈8LT2ωε0ab(2πa)2x′2∝x′2Vout 0=√ZP=ωq2√ZQext RQ0exp(−ω2σ2t2)∝q|x′| (2) 式中: RQ0为归一化分路阻抗,L为腔体纵向长度,T=sin(ωL/2c)ωL/2c为渡越时间因子,Z为负载特征阻抗,P为输出到负载的功率,q为束团电荷量,Qext为外部品质因数。由此可见在近轴条件下,腔式BPM输出信号幅度Vout0与偏移轴线距离x′成正比,偏移方向则由参考腔提供相位信息进行判断。考虑衰减时间常数τ=QL/ω,则信号形式为
Vout =Vout 0exp(−t2τ)sin(ωt+φ) (3) 式中: QL为有载品质因数,φ为初始相位。
1.2 激光驱动质子束流产生的输出信号
1.2.1 激光驱动质子束流位置偏移的信号
激光加速产生的质子束团具有短脉冲、宽能谱和大散角的特点。因为初始激光打靶的靶点只有μm量级,通常其会具有比较小的发射度,例如内华达大学Cowan等人测量得到的发射度εn, rms~0.004 mm·mrad[8]。在经过束流收集系统收集后, 大散角束流变为小散角束流。因为质子束团有较大的横向面积,所以输出信号必须要考虑束团横向分布的影响。为了得到在较大束斑情况下的信号,我们利用激光驱动质子束团层流的性质[8],将其分割为许多微束(小到可以不计横向分布),叠加所有微束的信号即可得到整体束团的位置信号。此时每个微束由于位置偏移产生的输出信号可以表达为
Vp-bunch =Aqbunch |x|sin(ωt+φ) (4) 式中:A=ω2√ZQext8LT2ωε0ab(2πa)2exp(−ω2σ2t2)exp(−t2τ),qbunch为该微束的电荷量,x为该微束的偏移轴线位置。对所有微束的信号在横向分布条件下加权叠加,即可以得到实际质子束团的输出信号。由于此时所有微束是同时入射进腔式BPM,故在同一侧的束团输出信号同向(φ相等),而在相反两侧的束团输出信号反相(φ相差π),则可以借用偏移位置x的符号来表征相位关系。因此我们考虑相位关系后统一对
Vp-bunch =A1qbun ch x (5) 进行叠加,就可以得到大束斑束团的输出信号, 其中A1=Asin(ωt),并不失一般性取正向入射时φ=0。
为了对微束信号进行加权叠加,我们要得到整个束斑的横向分布。基于北京大学激光等离子体实验室(有目前世界首条且唯一在运行的激光驱动质子束线)的束流实验结果,如图 1所示,发射度约为εn~0.4 mm·mrad,我们做出质子束团的横向分布假设:总电量为Q的束团的横向分布f(x, y)是以位置(x0, y0)为对称中心的中心对称分布,对称中心(x0, y0)定义为大束斑束流的横向偏移位置。此时,叠加所有微束的信号有
Vp-total =∬DA1Qf(x,y)x ds,(x,y)∈D (6) 图 1 (a) RCF测得约25%能散束流横向分布,(b)闪烁体测得约4%能散束流横向分布,(c)MCP测得约1%能散束流横向分布Figure 1. Experimental results: (a) beams' transverse distribution with about 25% energy dispersion on RCF, (b) beams' transverse distribution with about 4% energy dispersion on scintillator, and (c) beams' transverse distribution with about 1% energy dispersion on MCP做变换x=x0+x′,y=y0+y′,则有
Vp-total =∬D′A1Qf(x0+x′,y0+y′)(x0+x′)ds′=A1Qx0∬D′f(x0+x′,y0+y′)ds′+A1Q∬D′f(x0+x′,y0+y′)x′ds′=A1Qx0(x′,y′)∈D′ (7) 其中∬D′f(x0+x′,y0+y′)ds′为分布函数平移后的积分,值仍为1,∬D′(f(x0+x′,y0+y′)x′ds′为二元偶函数与奇函数的乘积在对称定义域上的积分,值为0。因此可以得出,虽然束斑横向面积较大,如果其横向分布满足关于某点(x0, y0)中心对称时,其入射进腔式BPM产生的输出信号等价于相同电量的微束从点(x0, y0)入射进腔式BPM产生的输出信号。
1.2.2 激光驱动质子束流发散角的信号
为了得到质子束团发散角的输出信号,我们仍将其分割为许多微束(小到可以不计横向分布),则叠加所有微束倾斜入射产生的信号即为发散角的信号。激光驱动质子束团在经过收集系统后,散角变小,所有微束均满足较小的倾斜角条件(tanθ≈θ),则微束倾斜θ入射产生的信号可以表示为[7]
Va-bunch =A2qbunch θ (8) 式中: A2=ω2√ZQext8T2ωε0ab(2πa)2L2√L2sin(ωL4c)exp(−ω2σt22)exp(−t2τ)cos(ωt),其中已由θ的符号表征不同倾斜角入射后的信号之间相位关系(同相或反相),并且比较A1和A2的振荡项可知位置信号与倾斜信号相位正交(π/2相差)。
为了对各倾斜微束信号进行加权叠加,我们要得到整个束斑的倾斜角分布。束团倾斜入射进腔式BPM的轨迹在x-z平面的投影如图 2,其中s为聚焦时的某个固定距离。根据上小节中束团横向分布关于点(x0, y0)中心对称的假设有
f(x0−x)=f(x0+x) (9) 再依据激光驱动质子束团小倾斜角条件(tanθ≈θ)有
f(sθ0−sθ)=f(sθ0+sθ) (10) 即束团的倾斜角分布g(θ)=f(sθ) 为关于角θ0的对称分布。此时,叠加所有微束的倾斜信号有
Va-total =∫θmin (11) 做类似上小节的变换,最终可以得到Va-total=A2Qθ0。因此可以得出,得益于激光驱动质子束团较小的发射度,在束团横向分布中心对称假设下,其入射进腔式BPM后由于发散角产生的输出信号等价于相同电量的微束倾斜θ0入射时产生的输出信号。
综合位置信号和发散角信号,合输出信号为
V_{\text {total }}=V_{\mathrm{p} \text {-total }}+\mathrm{j} V_{\mathrm{a} \text {-total }} (12) 两信号之间正交(存在π/2相差),故通过相位探测可以分别获得对应信号。通过上述推导,我们证明了在腔式BPM中中心对称束流以小发散角入射与微束倾斜入射的等价性,继而实现了可以利用微束来模拟设计应用于激光驱动质子束流的腔式BPM。
2. 腔式BPM设计
2.1 腔式BPM的参数设计及仿真
为了防止信号串扰[9],我们采用矩形谐振腔方案[10]。考虑到质子束流横向尺寸大,需要在大谐振腔上开足够大的管道直径,我们将工作频率选定为f210=1.32 GHz和f120=1.42 GHz。根据理想矩形谐振腔本征模频率公式
f_{\text {m n } \mathrm{p}}=\frac{c}{2} \sqrt{\left(\frac{m}{a}\right)^2+\left(\frac{n}{b}\right)^2+\left(\frac{p}{l}\right)^2}, \quad m n \neq 0 (13) 得到腔体长边a=260.84 mm,短边b=230.86 mm。使用CST微波模块进行本征模仿真得到矩形腔两个偶极模的磁场如图 3。
接下来,我们不断调整腔体的长度L和管道直径D,来仿真式(2)中的归一化分路阻抗 \frac{R}{Q_0}。当D=80 mm时,改变L分别得到的两个偶极模在不同位置处的 \frac{R}{Q_0},如图 4。同时不断调整波导和同轴耦合件参数,来仿真优化式(2)中的外部品质因数Qext,使两个模式的Qext均为2000左右。
由于对谐振腔的开孔等因素,本征模频率会出现漂移。我们通过微调腔体及波导位置等将频率修正回设计值,最终腔式BPM优化尺寸如表 1所示,其中波导采用标准波导BJ-18型,其位置指波导边缘距离谐振腔中轴线距离。腔式BPM的理论参数如表 2所示,其中 \frac{R}{Q_0}为1 mm处归一化分路阻抗。
表 1 腔式BPM优化尺寸Table 1. Optimum geometry of cavity BPMa /mm b/mm L/mm D/mm waveguide length X/mm waveguide position X/mm waveguide length Y/mm waveguide position Y/mm 241.6 210 75 80 250 90.5 280 76 表 2 腔式BPM理论参数Table 2. Theoretical parameters of cavity BPMf110/GHz (R/Q0)110/Ω f210/GHz (Qext)210 (R/Q0)210/Ω f120/GHz (Qext)120 (R/Q0)120/Ω 0.9167 224.5 1.320 1955 0.1088 1.420 2103 0.1328 2.2 适用范围讨论及分辨率估计
激光驱动质子束流初始能谱宽,经过束线选能后不同能量的束团电量可达到几十pC到nC不等[11-12],其脉宽σt主要由中心能量E、能散Δ及传输距离s决定,可以估算为
\begin{aligned} &\sigma_t \approx \frac{1}{2.3556} \frac{s}{v_2-v_1} \approx \\ &\frac{1}{2.3556} \frac{s}{c}\left(1 / \sqrt{1-\left(\frac{\varepsilon_0}{\varepsilon_0+E(1-{\mathit{\Delta}} / 2)}\right)^2}-1 / \sqrt{1-\left(\frac{\varepsilon_0}{\varepsilon_0+E(1+{\mathit{\Delta}} / 2)}\right)^2}\right) \end{aligned} (14) 式中: 2.355 6为高斯分布半高宽与特征脉宽σt的比例常数,v2, v1为正负最大能散处速度,ε0为质子静止能量,取s=10 m处,代入激光驱动质子束流的能量及能散信息可获得脉宽分布如图 5(a)。进一步,将式(2)变化得到极限分辨率处归一化分路阻抗表达式
\frac{R}{Q_0}=\left(\frac{2 V_{\text{out} 0}}{\omega q \exp \left(-\frac{\omega^2 \sigma_{\mathrm{t}}^2}{2}\right)}\right)^2 \frac{Q_{\text {ext }}}{Z} (15) 由于存在暗电流等干扰,参考文献[10]中的测量值,在分辨率极限处至少应有Vout0=5 mV,同时取Z=50 Ω,代入激光驱动质子束流的电量和脉宽信息可得 \frac{R}{{{Q_0}}}的分布如图 5(b),与图 3中各位置处 \frac{R}{{{Q_0}}}值比较可以得到不同激光驱动质子束团条件下能达到的位置分辨率。
以目前PW级激光加速系统的典型质子加速结果[11]为例进行分辨率预估,得到表 3, 其中σt为束团按1%能散传输10 m时的脉宽,结果表明对于不同能量束团可实现好于1 mm的位置分辨率。
表 3 分辨率估计Table 3. Resolution estimationE/MeV σt/ps q/pC (R/Q0)210/Ω resolution/mm 30 273 500 0.039 2 0.6 50 208 150 0.050 7 0.8 100 142 50 0.093 2 1 3. 结论
本文推导了应用于激光驱动质子束流的腔式BPM的输出信号,证明了当束流横向分布满足关于点(x0, y0) 中心对称分布、束团发散角较小时,其产生的输出信号等价于同电量微束从点(x0, y0) 倾斜θ0入射进腔式BPM时产生的输出信号。基于这个原理,我们使用CST软件进行了应用在激光驱动质子装置上的腔式BPM的设计仿真,采用了矩形腔波导耦合的方案,且选择的工作频率较低、间隔较大,最大限度防止信号串扰。对本腔式BPM在不同激光驱动质子束团参数条件下可获得的分辨率进行了讨论,结果表明其应用在PW激光加速系统中可以实现较好的位置分辨率。
-
表 1 加速器真空环境需求表
Table 1. Requirements of vacuum system
accelerator vacuum degree/ Pa ion source(IS)low energy beam transport line(LEBT) 2.0×10−3 radio frequency quadrupole(RFQ) 1.0×10−5 medium energy beam transport line(MEBT) 1.0×10−5 drift tube linac(DTL) 1.0×10−5 linac to ring beam transport line(LRBT) 1.0×10−5 rapid cycling synchrotron(RCS) 5.0×10−6 ring to target beam transport line(RTBT) 1.0×10−5 表 2 真空系统主要设备统计
Table 2. The vacuum system devices
device molecular pump ion pump gauge valve quantity 7 147 30 34 表 3 真空数据电气接口表
Table 3. Serial ports of vacuum system
device model quantity analog interface COM TPG-300 54 not used RS-422 MKS 937B 9 used not used ion pump power supply 147 none RS-485 表 4 其他系统设备真空压力联锁信号
Table 4. Interlock signals of other systems
system ion source power supply linac RF DTL quantity 2 10 4 9 -
[1] Wei Jie, Chen Hesheng, Chen Yanwei, et al. China Spallation Neutron Source: Design, R&D, and outlook[J]. Nuclear Instruments & Methods in Physics Research Section A—Accelerators Spectrometers Detectors and Associated Equipment, 2009, 600(1): 10-13. [2] Wei Jie, Fu Shinian, Tang Jingyu, et al. China Spallation Neutron Source—An overview of application prospects[J]. Chinese Physics C, 2009, 33(11): 1033-1042. doi: 10.1088/1674-1137/33/11/021 [3] Wang Sheng, Fang Shouxian, Fu Shinian, et al. Introduction to the overall physics design of CSNS accelerators[J]. Chinese Physics C, 2009, 33(S2): 1-3. doi: 10.1088/1674-1137/33/S2/001 [4] 刘华昌, 彭军, 巩克云, 等. 中国散裂中子源漂移管直线加速器研制进展[J]. 原子能科学技术, 2015, 49(S2):556-559. (Liu Huachang, Peng Jun, Gong Keyun, et al. DTL construction status of China Spallation Neutron Source[J]. Atomic Energy Science and Technology, 2015, 49(S2): 556-559 [5] 董海义, 宋洪, 李琦, 等. 中国散裂中子源(CSNS)真空系统研制[J]. 真空, 2015(4):1-6. (Dong Haiyi, Song Hong, Li Qi, et al. Research and development of vacuum system for China Spallation Neutron Source[J]. Vacuum, 2015(4): 1-6 [6] 黄涛, 董海义, 王鹏程, 等. CSNS离子泵性能测试[J]. 真空, 2017(6):43-46. (Huang Tao, Dong Haiyi, Wang Pengcheng, et al. Performance test of CSNS ion pump[J]. Vacuum, 2017(6): 43-46 [7] 王鹏程, 黄涛, 刘佳明, 等. 中国散裂中子源(CSNS)LRBT输运线真空系统[J]. 真空, 2019(5):21-25. (Wang Pengcheng, Huang Tao, Liu Jiaming, et al. The vacuum system of LRBT at CSNS[J]. Vacuum, 2019(5): 21-25 [8] 李刚, 赵籍九, 陈伯飞, 等. BEPCⅡ储存环真空控制系统[J]. 原子能科学技术, 2010(4):504-507. (Li Gang, Zhao Jijiu, Chen Bofei, et al. Vacuum control system of BEPCⅡ storage ring[J]. Atomic Energy Science and Technology, 2010(4): 504-507 [9] 何泳成, 王春红, 王金灿, 等. 基于嵌入式IOC的CSNS真空控制系统样机研制[J]. 核电子学与探测技术, 2013(3):280-283. (He Yongcheng, Wang Chunhong, Wang Jincan, et al. Development of a prototype vacuum control system based on embedded IOC for CSNS[J]. Nuclear Electronics & Detection Technology, 2013(3): 280-283 doi: 10.3969/j.issn.0258-0934.2013.03.004 [10] 赵籍九, 尹兆升. 粒子加速器技术[M]. 北京: 高等教育出版社, 2006.Zhao Jijiu, Yin Zhaosheng. Particle accelerator technology[M]. Beijing: Higher Education Press, 2006 [11] 吴煊, 王春红, 王金灿, 等. 基于Modbus RS232的数字电源控制器的远控设计与实现[J]. 核电子学与探测技术, 2012(8):867-870. (Wu Xuan, Wang Chunhong, Wang Jincan, et al. The design and implementation of the remote control of a digital power supply control module based on Modbus/RS232[J]. Nuclear Electronics & Detection Technology, 2012(8): 867-870 doi: 10.3969/j.issn.0258-0934.2012.08.001 [12] 张玉亮. 基于StreamDevice的VME机箱监测及其应用研究[J]. 核电子学与探测技术, 2014:834-836. (Zhang Yuliang. Study of StreamDevice based VME chassis monitoring and its implementation[J]. Nuclear Electronics & Detection Technology, 2014: 834-836 doi: 10.3969/j.issn.0258-0934.2014.07.007 [13] 何泳成, 李刚, 康明涛, 等. CSNS反角白光中子源束线控制系统研制[J]. 强激光与粒子束, 2018, 30:114006. (He Yongcheng, Li Gang, Kang Mingtao, et al. Development of CSNS white neutron source beam line control system[J]. High Power Laser and Particle Beams, 2018, 30: 114006 doi: 10.11884/HPLPB201830.180193 -