BNCT02加速器机器保护系统设计

何泳成; 吴煊; 张玉亮; 朱鹏; 薛康佳; 王林; 李明涛; 傅世年

doi:10.11884/HPLPB202537.240153

BNCT02加速器机器保护系统设计

doi: 10.11884/HPLPB202537.240153

何泳成^{1, 2, 3,},
吴煊^{1, 2, 3},
张玉亮^{1, 2, 3},
朱鹏^{1, 2, 3},
薛康佳^{1, 2},
王林^{1, 2, 3},
李明涛^{1, 2, 3},
傅世年^{1, 2, 3}

1.
中国科学院高能物理研究所，北京 100049
2.
散裂中子源科学中心，广东东莞 523803
3.
中国科学院大学，北京 100049

基金项目: 广东省引进创新创业团队项目(2017ZT07S225)

详细信息

作者简介:
何泳成，heyongcheng@ihep.ac.cn

中图分类号: TL503.6
计量
- 文章访问数: 282
- HTML全文浏览量: 113
- PDF下载量: 40
- 被引次数: 2
出版历程
- 收稿日期: 2024-05-10
- 修回日期: 2024-10-10
- 录用日期: 2024-10-18
- 网络出版日期: 2024-11-19
- 刊出日期: 2025-12-13

Design of machine protection system for BNCT02 accelerator

He Yongcheng^{1, 2, 3
,},
Wu Xuan^{1, 2, 3},
Zhang Yuliang^{1, 2, 3},
Zhu Peng^{1, 2, 3},
Xue Kangjia^{1, 2},
Wang Lin^{1, 2, 3},
Li Mingtao^{1, 2, 3},
Fu Shinian^{1, 2, 3}

1.
Institute of High Energy Physics, Chinese Academy of Sciences (CAS), Beijing 100049, China
2.
Spallation Neutron Source Science Center, Dongguan 523803, China
3.
University of Chinese Academy of Sciences, Beijing 100049, China

摘要

摘要: 中国科学院高能物理研究所第二台硼中子俘获治疗（BNCT02）加速器主要由1台离子源、1条低能束流传输线、1台射频四极加速器和3条高能束流传输线组成。为了保障BNCT02加速器的安全运行，设计了基于横河PLC和实验物理及工业控制系统（EPICS）软件工具包的机器保护系统。为了增强安全性，该系统采用了冗余设计，由两套完全独立且主要输入、输出信号一致的子系统构成。测试结果表明，BNCT02加速器机器保护系统的响应时间小于1.6 ms，且具有稳定可靠性高的特点，满足BNCT02加速器运行的需要。
- 硼中子俘获治疗 /
- 加速器 /
- 机器保护系统 /
- 实验物理及工业控制系统
Abstract: The RF linac dedicated to boron neutron capture therapy (BNCT02) in our institute is mainly composed of an ion source, a low energy beam transport line, a radio frequency quadrupole accelerator and three high energy beam transport lines. To ensure the safe operation of the BNCT02 accelerator, a machine protection system (MPS) was designed based on Yokogawa PLC and Experimental Physics and Industrial Control System (EPICS) software toolkit. The MPS adopts a redundant design, consisting of two completely independent subsystems with consistent main input and output signals. The test results show that the response time of the MPS is less than 1.6 ms, and it has the characteristics of high stability and reliability, which meets the operational requirements of the BNCT02 accelerator.
- boron neutron capture therapy /
- accelerator /
- machine protection system /
- Experimental Physics and Industrial Control System

HTML全文

间断伽辽金有限元(DG)方法是近年来计算电磁学领域发展最为迅速的方法之一^[1-2]。它具有精度高^[3]、不规则网格易处理^[4]、时间步长选取灵活^[5]、并行计算易实现^[6]等优点，在雷达、通讯、微波、集成光学等领域有重要的应用^[7-9]。为了提高间断元方法的计算精度，需要加密模型划分的单元格(h方法)或者提高基函数插值阶数(p方法)。一般来说, p方法要比h方法收敛更快。基函数阶数的增加必然会带来质量矩阵的增大，从而使存储开销急剧增长(增长速度与阶数的三次方成正比)。为了解决这一困难，人们提出了质量集中算法。其基本思想是，将基函数的插值点与控制方程弱形式中的数值积分点重合，从而达到增加质量矩阵中零元素的目的。对标量基函数而言，可以使质量矩阵对角化，而对矢量基函数，质量矩阵可化为3×3的块对角矩阵。显然，对两种类型的基函数，质量集中都极大地降低了存储和计算开销。文献[10]介绍了一种新的质量集中算法，并将其成功应用到波动问题上。关于该方法在电磁学领域中的应用，以及新算法与传统质量集中算法之间的差异性需要进一步探讨。本文通过数值算例对这两种算法进行了比较，重点考察了两种方法的精度和稳定特征，为质量集中算法在DG方法中的应用提供参考。

1. 间断有限元算法

1.1 二维问题的质量矩阵

本文仅讨论二维问题。TE波模式下(TM模式类似)，无源区间的二维麦克斯韦方程组可以表示为

$\left\{\begin{array}{l} \varepsilon \frac{\partial E_x}{\partial t}=\frac{\partial H_z}{\partial y} \\ \varepsilon \frac{\partial E_y}{\partial t}=-\frac{\partial H_z}{\partial x} \\ \mu \frac{\partial H_z}{\partial t}=\frac{\partial E_y}{\partial x}-\frac{\partial E_x}{\partial y} \end{array}\right.$

(1)

式中: E_x和E_y分别是电场沿x轴和y轴的分量; H_z是磁场的分量; ε和μ分别是腔体内的电容率和磁导率。

将计算区域Ω划分为K个四边形网格单元，每个单元记作Ω^k。单元Ω^k内的电磁场量E和H可以近似为

$\begin{aligned} \boldsymbol{E}_x^k(x, y, t) &=\sum\limits_{i=0}^N e_{x, i}^k(t) \boldsymbol{\varPhi}_i(x, y) \\ \boldsymbol{E}_y^k(x, y, t) &=\sum\limits_{i=0}^N e_{y, i}^k(t) \boldsymbol{\varPhi}_i(x, y) \\ \boldsymbol{H}_z^k(x, y, t) &=\sum\limits_{i=0}^N h_{z, i}^k(t) \boldsymbol{\varPhi}_i(x, y) \end{aligned}$

(2)

式中: (e_{x, i}^k, e_{y, i}^k, h_{z, i}^k)为待求的电磁场量; Φ(x, y)是二维基矢量。将电磁场量的离散形式代入到方程组(1)中，用基矢量点乘公式两端，并在Ω^k上积分，可以得到

$\begin{gathered} \varepsilon \sum\limits_{i=0}^N \phi_{i, j}^k \frac{\partial e_{x, i}^k}{\partial t}=\sum\limits_{i=0}^N \xi_{i, j}^k h_{z, i}^k \\ \varepsilon \sum\limits_{i=0}^N \phi_{i, j}^k \frac{\partial e_{y, i}^k}{\partial t}=-\sum\limits_{i=0}^N \delta_{i, j}^k h_{z, i}^k \\ \mu \sum\limits_{i=0}^N \phi_{i, j}^k \frac{\partial h_{z, i}^k}{\partial t}=\sum\limits_{i=0}^N \delta_{i, j}^k e_{y, i}^k-\sum\limits_{i=0}^N \xi_{i, j}^k e_{x, i}^k \end{gathered}$

(3)

其中

$\begin{gathered} (\cdot , \cdot)_{\Omega^k}=\int_{\Omega^k} \cdot \mathrm{d} s \\ \phi_{i, j}^k=\left(\boldsymbol{\varPhi}_i^k, \boldsymbol{\varPhi}_j^k\right)_{\Omega^k} \\ \delta_{i, j}^k=\left(\boldsymbol{\varPhi}_i^k, \frac{\partial \boldsymbol{\varPhi}_j^k}{\partial x}\right)_{\Omega^k}, \quad \xi_{i, j}^k=\left(\boldsymbol{\varPhi}_i^k, \frac{\partial \boldsymbol{\varPhi}_j^k}{\partial y}\right)_{\Omega^k} \end{gathered}$

两次使用分部积分可以进一步得到间断有限元法的强形式

$\begin{gathered} \sum\limits_{i=0}^N\left(M_{i j}^k \cdot \frac{\partial e_{x, i}^k}{\partial t}-S_{i j}^{y, k} \cdot h_{z, i}^k\right)=\int_{\partial \Omega^k} \boldsymbol{\varPhi}_i^k(x, y) \cdot \hat{\boldsymbol{n}} \times\left(\boldsymbol{H}_z^{k *}-\boldsymbol{H}_z^k\right) \mathrm{d} l \\ \sum\limits_{i=0}^N\left(M_{i j}^k \cdot \frac{\partial e_{y, i}^k}{\partial t}+S_{i j}^{x, k} \cdot h_{z, i}^k\right)=\int_{\partial \Omega^k} \boldsymbol{\varPhi}_i^k(x, y) \cdot \hat{\boldsymbol{n}} \times\left(\boldsymbol{H}_z^{k *}-\boldsymbol{H}_z^k\right) \mathrm{d} l \\ \sum\limits_{i=0}^N\left(M_{i j}^k \cdot \frac{\partial h_{z, i}^k}{\partial t}-S_{i j}^{x, k} \cdot e_{y, i}^k+S_{i j}^{y, k} \cdot e_{x, i}^k\right)= \\ \int_{\partial \Omega^k} \boldsymbol{\varPhi}_i^k(x, y)\left[\hat{\boldsymbol{n}} \times\left(\boldsymbol{E}_x^{k *}-\boldsymbol{E}_x^k\right)+\hat{\boldsymbol{n}} \times\left(\boldsymbol{E}_y^{k *}-\boldsymbol{E}_y^k\right)\right] \mathrm{d} l \end{gathered}$

(4)

式中: $\hat{\boldsymbol{n}}$ 为∂Ω^k外向法向量，(H_z^k*, E_x^k*, E_y^k*)是边界∂Ω^k上的数值通量。

M^k表示该单元上的质量矩阵，S^{x, k}, S^{y, k}表示相应的刚度矩阵，具体表示为

$\begin{aligned} M_{i j}^k &=C \int_{\Omega^k} \boldsymbol{\varPhi}_i^k(x, y) \boldsymbol{\varPhi}_j^k(x, y) \mathrm{d} s \\ S_{i j}^{x, k} &=\int_{\Omega^k} \boldsymbol{\varPhi}_i^k(x, y) \frac{\partial \boldsymbol{\varPhi}_j^k(x, y)}{\partial x} \mathrm{~d} s \\ S_{i j}^{y, k} &=\int_{\Omega^k} \boldsymbol{\varPhi}_i^k(x, y) \frac{\partial \boldsymbol{\varPhi}_j^k(x, y)}{\partial y} \mathrm{~d} s \end{aligned}$

(5)

式中: C表示ε或者μ。

由于积分单元Ω^k并不是一个标准的四边形，所以我们使用坐标变换，通过雅克比行列式|J_k|将其转换到[-1, 1]²的正方形单元 $\hat{\Omega}$ 上积分，质量矩阵的元素表示为

$M_{i j}^k=\int_{\hat{\Omega}} C \boldsymbol{\varPhi}_i^k \boldsymbol{\varPhi}_j^k\left|\boldsymbol{J}_k\right| \mathrm{d} s=C \sum\limits_{\alpha=0}^{N^{\prime}} \omega_\beta \sum\limits_{\beta=0}^{N^{\prime}} \omega_\alpha \boldsymbol{\varPhi}_i\left(x_\alpha, y_\beta\right) \boldsymbol{\varPhi}_j\left(x_\alpha, y_\beta\right)\left|\boldsymbol{J}_k\right|$

(6)

式中: ω_α, ω_β是对应积分点的权重。

通过观察M_ij^k的形式可以发现，对于选用p=N阶插值多项式的单元，需要一个(N+1)²×(N+1)²的方阵来存储。由于系数C和雅克比行列式|J_k|不同，每个单元的M^k也不尽相同。因此，对于拥有n个单元的网格，程序中总共需要n个(N+1)²×(N+1)²的二维数组来存储质量矩阵。

1.2 两种方法质量矩阵存储量比较

为了减少质量矩阵带来的内存开销，传统的质量集中间断有限元算法是选用相同的插值点和积分点来近似计算，则二维基矢量Φ (x, y)表示为

$\varPhi_{i j}(x, y)=\phi_i^N(x) \phi_j^N(y)=\prod\limits_{k=0, k \neq i}^N \frac{x-x_k}{x_i-x_k} \cdot \prod\limits_{k=0, k \neq j}^N \frac{y-y_k}{y_j-y_k}$

(7)

当插值点与积分点重合时有

$\sum\limits_{\alpha=0}^N \omega_\beta \sum\limits_{\beta=0}^N \omega_\alpha \boldsymbol{\varPhi}_i\left(x_\alpha, y_\beta\right) \boldsymbol{\varPhi}_j\left(x_\alpha, y_\beta\right) \begin{cases}\neq 0, & i=j \\ =0, & i \neq j\end{cases}$

(8)

由此方法得到的质量矩阵M^k，是一个(N+1)²阶的对角矩阵，只需要用一个长度为(N+1)²的一维数组来存放其对角线上的元素，大大减少了存储空间的损耗。

文献[]中介绍的Weight-Adjust算法与之类似，在其基础之上给出了另外一种 $\hat{\Omega}$ 单元上逆质量矩阵(M^k)^-1的近似求法，可以表示为

$\left(\boldsymbol{M}^k\right)^{-1} \approx \boldsymbol{M}^{-1} \boldsymbol{M}_{1 / \varepsilon\left|J_k\right|} \boldsymbol{M}^{-1}$

(9)

其中

$M_{i j}=\int_{\hat{\Omega}} \boldsymbol{\varPhi}_i \boldsymbol{\varPhi}_j \mathrm{~d} s$

(10)

$\left(M_{1 / \varepsilon\left|J_k\right|}\right)_{i j}=\int_{\hat{\Omega} \varepsilon} \frac{\boldsymbol{\varPhi}_i \boldsymbol{\varPhi}_j}{\left|\boldsymbol{J}_k\right|} \mathrm{d} s=\sum\limits_{\beta=0}^N \omega_\beta \sum\limits_{\alpha=0}^N \omega_a \frac{\boldsymbol{\varPhi}_i\left(x_\alpha, y_\beta\right) \boldsymbol{\varPhi}_j\left(x_\alpha, y_\beta\right)}{\varepsilon\left|\boldsymbol{J}_k\right|}$

(11)

通过观察式(11)我们可以发现，M矩阵同样需要一个(N+1)²阶的方阵来存储，但是由于每个参考单元上的M矩阵均相同，拥有n个单元的网格可以共用同一个(N+1)²阶方阵，造成的开销可以忽略不计。当选用相同的插值点和积分点时，矩阵M_{1/ε|J_k|}是一个(N+1)²阶的对角矩阵，只需要用一个长度为(N+1)²的一维数组来存放其对角线上的元素。

通过上面的分析，两种近似方法所需要的质量矩阵存储量相同，下文的数值算例部分则从精度以及计算效率等方面比较这两种方法。

2. 数值算例和讨论

2.1 谐振腔

首先考虑计算一个有解析解的方形谐振腔问题，来分析两种方法的精度以及计算效率的差异。设计一个[0, 1]² m²的方形区域，划分为12个不规则的四边形单元。在区域内给定基模作为初始值后，两种方法均使用相同的时间步长(时间积分使用Runge-Kutta法^[4])和6阶插值多项式进行模拟，并记录测量点(0.5, 0.85)处磁场分量H_z随着时间变化的值。最后分别将两种方法模拟得到的数据与参考值进行比较，得到绝对误差。

在图 1(a)中，给出了测量点处的参考时域波型以及两种近似方法下得到的实际波型，图 1(b)中进一步得到测量点处数值解与参考解比较的绝对误差。可以明显地看出，传统DG算法模拟结果的误差振幅大约为WADG算法的两倍。

图 1 测量点处得到的模拟结果与参考值的误差比较

Figure 1. Comparison of simulation results with reference values at measuring point

下载: 全尺寸图片幻灯片

另外，通过简单地调整程序的时间步长，我们发现在保证模拟结果收敛于解析解的前提下，传统的质量集中算法可以使用的最大时间步长是Weight-Adjust算法的4倍。

2.2 滤波器

下面我们将用上述两种方法去模拟计算一个带通滤波器的散射参数，具体的截面模型如图 2所示，各项参数均参考文献[11]。整个截面模型包括左右两端的PML吸收层，并且被划分为3488个长方形单元。

滤波器的散射系数可以表示为

$S_{11}=\sqrt{\frac{p_{\text {ref }}}{p_{\text {inc }}}}, \quad S_{21}=\sqrt{\frac{p_{\text {tra }}}{p_{\mathrm{inc}}}}$

(12)

图 2 带通滤波器截面的几何模型

Figure 2. Geometrical model of bandpass filter's cross-section

下载: 全尺寸图片幻灯片

为了方便计算S参数所需的系数，我们添加了TF/SF界面，并在该界面上设置了TE₁₀的模式波。通过电磁场的能流密度公式，我们求得TF/SF面左侧面上的反射系数p_ref，右侧面上的发射系数p_inc以及滤波器与右侧PML交界面上的透射系数p_tra。

我们将最后模拟结果，即不同频率模式波下的散射系数表示在图 3中。两种方法在长方形网格下模拟得到的散射系数结果差别较小，并且与文献[11]中给出的模拟结果基本一致。

图 3 带通滤波器的散射系数

Figure 3. Scattering parameters of bandpass filter

下载: 全尺寸图片幻灯片

2.3 光子晶体

近年来，间断有限元法在光子晶体领域的应用比较广泛，为了进一步探讨两种算法处理更加复杂的问题时的性能，我们将其应用到模拟光子晶体功率分配器内的电磁场^[12]。如图 4所示，整个截面被划分为3488个四边形单元，每个单元均选用4阶插值多项式。模型的截面尺寸如图 5所示，包括一层厚度为1 μm的PML吸收层以及其最外侧的PEC边界。

图 4 光子晶体的截面网格

Figure 4. Meshes of photonic crystal's cross-section

下载: 全尺寸图片幻灯片

图 5 光子晶体功率分配器的截面模型以及电磁场分布

Figure 5. Cross-section geometry and electromagnetic field distribution of photonic crystol power divider

下载: 全尺寸图片幻灯片

光子晶体中圆柱体的半径为0.18 μm，其材料的相对电容率和相对磁导率分别为ε(r)=11.56和μ(r)=1。相邻圆柱体中心的间距为1 μm，其中的填充介质相对介电常数和相对磁导率分别为ε(r)=1和μ(r)=1。在图中标注的位置我们给定了一个高斯脉冲 $E(r, t)=\hat{z} G (t)$ ，其中G(t)=e^－t²/2τ²cos(ωt)，频率ω=2.513 3×10⁶c(c是光速)，时延τ=1.591 5。

分别用两种方法进行模拟，记录下每点处的电磁场大小，并进行绘制。这里给出了Weight-Adjusted算法下，光子晶体功率分配器截面电磁场模拟结果，如图 5中所示，本文模拟得到的电磁场分布与文献[12]给出的参考分布相同。当然，传统的质量集中算法也能得到相似的模拟结果，这里便不再使用图片说明。

为了进一步说明两种方法在模拟精度上的差异，我们选取分配器通道内的一点，分别记录两组模拟中该点处的时域波型，并与参考数据进行比较(参考数据由高阶间断有限元法模拟得到)。从图 6给出的两种模拟方法各自的绝对误差中可以看出，Weight-Adjusted算法获得的数据精确度更高。

图 6 通道内某测量点处两组数据的绝对误差

Figure 6. Absolute error of two sets of data at measuring point

下载: 全尺寸图片幻灯片

3. 结论

Weight-Adjusted算法是一种能够简化质量矩阵，从而节省内存开销的近似算法。本文首次尝试将其应用到处理四边形网格上的二维时域电磁场问题，并与该领域内一种已经成熟使用的质量集中算法在精度、内存、计算效率等方面进行了比较。谐振腔、带通滤波器以及光子晶体内电磁场的数值算例的结果表明，Weight-Adjusted算法的精度明显更高。

图 1 BNCT02结构示意图

Figure 1. Structural diagram of BNCT02

下载: 全尺寸图片幻灯片

图 2 机器保护系统与相关设备及系统间的关系

Figure 2. Relationship between machine protection system (MPS) and related equipment and systems

下载: 全尺寸图片幻灯片

图 3 机器保护系统结构图

Figure 3. Structural diagram of the machine protection system

下载: 全尺寸图片幻灯片

图 4 梯形图程序开发环境及部分程序

Figure 4. Development environment of a ladder diagram program and part of the program

下载: 全尺寸图片幻灯片

图 5 关键设备及汇总状态界面

Figure 5. Interface of key equipment and summary status

下载: 全尺寸图片幻灯片

图 6 响应时间测试

Figure 6. Test result of response time

下载: 全尺寸图片幻灯片

图 7 MPS1 EPICS IOC 状态

Figure 7. The status of MPS1 EPICS IOC

下载: 全尺寸图片幻灯片

表 1 机器保护系统的主要输入信号数量

Table 1. Numbers of main input signals of the machine protection system

subsystem	number of signals
subsystem	PS	VV	BS	RMS	PPS	TMS	TBS	MI	DDS	ITS	FPS	total
MPS1	47	8	2	10	5	4	2	4	2	2	2	88
MPS2	31	8	2	9	4	4	2	4	2	2	1	69

下载: 导出CSV

表 2 机器保护系统的主要输出信号数量

Table 2. Numbers of main output signals of the machine protection system

subsystem	number of signals
subsystem	ISC	RMS	MI	FPS	HEBTBPS01	total
MPS1	2	1	2	1	3	9
MPS2	2	1	2	1	0	6

下载: 导出CSV

表 3 机器模式与需就绪子区对应关系

Table 3. Relationship between machine mode and sub areas

machine mode	linac sub area	dump sub area	treatment room 1 sub area	treatment room 2 sub area
ion source mode	/	/	/	/
dump mode	√	√	/	/
treatment room 1 mode	√	/	√	/
treatment room 2 mode	√	/	/	√

下载: 导出CSV

表 4 机器模式与HEBTBPS01状态信号对应关系

Table 4. Relationship between machine mode and HEBTBPS01 status signal

machine mode	power on	power off	positive	negative	current raised	fault	error exceeding limit
ion source mode	/	/	/	/	/	/	/
dump mode	×	√	/	/	×	√	/
treatment room 1 mode	√	×	√	×	√	√	√
treatment room 2 mode	√	×	×	√	√	√	√

下载: 导出CSV

表 5 机器模式与PPS控制区对应关系

Table 5. Relationship between machine mode and PPS control areas

machine mode	accelerator control area	treatment room 1 control area	treatment room 2 control area	vertical beam transport line control area
ion source mode	√	/	/	√
dump mode	√	/	/	√
treatment room 1 mode	√	√	/	√
treatment room 2 mode	√	/	√	√

下载: 导出CSV

参考文献(19)

[1]	Nakagawa Y, Hatanaka H. Boron neutron capture therapy: clinical brain tumor studies[J]. Journal of Neuro-Oncology, 1997, 33(1): 105-115.
[2]	Nakagawa Y, Pooh K, Kobayashi T, et al. Clinical review of the Japanese experience with boron neutron capture therapy and a proposed strategy using epithermal neutron beams[J]. Journal of Neuro-Oncology, 2003, 62(1/2): 87-99.
[3]	Kato I, Ono K, Sakurai Y, et al. Effectiveness of BNCT for recurrent head and neck malignancies[J]. Applied Radiation and Isotopes, 2004, 61(5): 1069-1073. doi: 10.1016/j.apradiso.2004.05.059
[4]	傅世年, 梁天骄, 陈和生. BNCT中子源的研发现况与展望[J]. 科学通报, 2022, 67(14):1471-1478 doi: 10.1360/TB-2021-1254 Fu Shinian, Liang Tianjiao, Chen Hesheng. Status and outlook: Research and development on the neutron source for BNCT[J]. Chinese Science Bulletin, 2022, 67(14): 1471-1478 doi: 10.1360/TB-2021-1254
[5]	陈和生, 邢更妹. 硼中子俘获治疗的发展机遇与挑战: 含硼药物的研发[J]. 科学通报, 2022, 67(14):1469-1470 doi: 10.1360/TB-2022-0474 Chen Hesheng, Xing Gengmei. Opportunities and challenges in boron neutron capture therapy: development of boron-containing drugs[J]. Chinese Science Bulletin, 2022, 67(14): 1469-1470 doi: 10.1360/TB-2022-0474
[6]	欧阳华甫, 肖永川, 刘盛进, 等. BNCT加速器设计和调试[J]. 白城师范学院学报, 2020, 34(2):1-9 Ouyang Huafu, Xiao Yongchuan, Liu Shengjin, et al. Designing and commissioning of BNCT accelerator[J]. Journal of Baicheng Normal University, 2020, 34(2): 1-9
[7]	欧阳华甫, 肖永川, 曹秀霞, 等. BNCT02加速器设计及离子源调试[J]. 白城师范学院学报, 2022, 36(5):1-8 Ouyang Huafu, Xiao Yongchuan, Cao Xiuxia, et al. Designing of BNCT02 accelerator and commissioning of ECR ion source[J]. Journal of Baicheng Normal University, 2022, 36(5): 1-8
[8]	赵籍九, 尹兆升. 粒子加速器技术[M]. 北京: 高等教育出版社, 2006 Zhao Jijiu, Yin Zhaosheng. Particle accelerator technology[M]. Beijing: Higher Education Press, 2006
[9]	李明涛, 张玉亮, 袁月, 等. CSNS加速器一键开关机程序设计和实现[J]. 强激光与粒子束, 2023, 35:084001 doi: 10.11884/HPLPB202335.230014 Li Mingtao, Zhang Yuliang, Yuan Yue, et al. Design and implementation of one key on/off program at the accelerator of CSNS[J]. High Power Laser and Particle Beams, 2023, 35: 084001 doi: 10.11884/HPLPB202335.230014
[10]	张玉亮, 李明涛, 吴煊, 等. 加速器故障分析软件及其在CSNS应用[J]. 强激光与粒子束, 2021, 33:074004 doi: 10.11884/HPLPB202133.210001 Zhang Yuliang, Li Mingtao, Wu Xuan, et al. A software based post mortem system and its application in CSNS[J]. High Power Laser and Particle Beams, 2021, 33: 074004 doi: 10.11884/HPLPB202133.210001
[11]	朱文超, 林汉尚, 周泽然, 等. 基于嵌入式EPICS的中子监测仪研制[J]. 强激光与粒子束, 2021, 33:026001 doi: 10.11884/HPLPB202133.200168 Zhu Wenchao, Lin Hanshang, Zhou Zeran, et al. Development of a neutron monitor based on embedded EPICS[J]. High Power Laser and Particle Beams, 2021, 33: 026001 doi: 10.11884/HPLPB202133.200168
[12]	Zhang Yuliang, Jin Dapeng, Zhu Peng, et al. The accelerator control system of CSNS[J]. Radiation Detection Technology and Methods, 2020, 4(4): 478-491. doi: 10.1007/s41605-020-00203-y
[13]	Zhang Yuliang, Kang Mingtao, Jin Dapeng, et al. The run management system for CSNS[J]. Radiation Detection Technology and Methods, 2019, 3: 35. doi: 10.1007/s41605-019-0116-0
[14]	程司农, 张昭, 朱鹏, 等. 基于高精度时间数据的加速器故障分析样机设计[J]. 核技术, 2022, 45:100203 doi: 10.11889/j.0253-3219.2022.hjs.45.100203 Cheng Sinong, Zhang Zhao, Zhu Peng, et al. Design of accelerator failure analysis system prototype based on high precision time data[J]. Nuclear Techniques, 2022, 45: 100203 doi: 10.11889/j.0253-3219.2022.hjs.45.100203
[15]	洪春霞, 周平, 滑文强, 等. 上海光源小角X射线散射实验站自动校准系统研究[J]. 核技术, 2021, 44:010102 doi: 10.11889/j.0253-3219.2021.hjs.44.010102 Hong Chunxia, Zhou Ping, Hua Wenqiang, et al. Automatic calibration system of small angle X-ray scattering experiment station at SSRF[J]. Nuclear Techniques, 2021, 44: 010102 doi: 10.11889/j.0253-3219.2021.hjs.44.010102
[16]	曾孟麒, 尹亮, 尹聪聪, 等. 基于EPICS的SHINE束线站定时设备控制系统[J]. 核技术, 2023, 46:070102 doi: 10.11889/j.0253-3219.2023.hjs.46.070102 Zeng Mengqi, Yin Liang, Yin Congcong, et al. Development and implementation of an EPICS-based timing equipment control system for SHINE beamlines and endstations[J]. Nuclear Techniques, 2023, 46: 070102 doi: 10.11889/j.0253-3219.2023.hjs.46.070102
[17]	张根灿. 上海质子治疗装置注入器远程控制系统的构建[D]. 上海: 中国科学院大学(中国科学院上海应用物理研究所), 2018 Zhang Gencan. The construction of remote control system for proton injector of Shanghai proton therapy facility[D]. Shanghai: Shanghai Institute of Applied Physics, Chinese Academy of Science, 2018
[18]	温特. 基于EPICS的HUST-PTF束损探测研究[D]. 武汉: 华中科技大学, 2020 Wen Te. Design of an EPICS based beam loss detection for HUST-PTF[D]. Wuhan: Huazhong University of Science & Technology, 2020
[19]	何泳成, 王春红, 赵卓, 等. 基于F3RP61的嵌入式EPICS IOC应用研究与实现[J]. 核电子学与探测技术, 2012, 32(1):17-20 doi: 10.3969/j.issn.0258-0934.2012.01.004 HE Yongcheng, WANG Chunhong, ZHAO Zhuo, et al. Research and implementation of embedded EPICS IOC based on F3RP61[J]. Nuclear Electronics & Detection Technology, 2012, 32(1): 17-20 doi: 10.3969/j.issn.0258-0934.2012.01.004

施引文献

期刊类型引用(0)
其他类型引用(2)

资源附件(0)

访问统计

点击查看大图

图(7) / 表(5)

计量

文章访问数: 282
HTML全文浏览量: 113
PDF下载量: 40
被引次数: 2

1. 间断有限元算法
1.1 二维问题的质量矩阵
1.2 两种方法质量矩阵存储量比较
2. 数值算例和讨论
2.1 谐振腔
2.2 滤波器
2.3 光子晶体
3. 结论

姓名
邮箱
手机号码
标题
留言内容
验证码

留言板

BNCT02加速器机器保护系统设计

doi: 10.11884/HPLPB202537.240153

作者简介:
何泳成，heyongcheng@ihep.ac.cn

计量

Design of machine protection system for BNCT02 accelerator

1. 间断有限元算法

1.1 二维问题的质量矩阵

1.2 两种方法质量矩阵存储量比较

2. 数值算例和讨论

2.1 谐振腔

2.2 滤波器

2.3 光子晶体

3. 结论

期刊类型引用(0)

其他类型引用(2)

计量

目录

1. 间断有限元算法

1.1 二维问题的质量矩阵

1.2 两种方法质量矩阵存储量比较

2. 数值算例和讨论

2.1 谐振腔

2.2 滤波器

2.3 光子晶体

3. 结论

留言板

BNCT02加速器机器保护系统设计

doi: 10.11884/HPLPB202537.240153

作者简介: 何泳成，heyongcheng@ihep.ac.cn

计量

出版历程

Design of machine protection system for BNCT02 accelerator

1. 间断有限元算法

1.1 二维问题的质量矩阵

1.2 两种方法质量矩阵存储量比较

2. 数值算例和讨论

2.1 谐振腔

2.2 滤波器

2.3 光子晶体

3. 结论

期刊类型引用(0)

其他类型引用(2)

计量

出版历程

目录

1. 间断有限元算法

1.1 二维问题的质量矩阵

1.2 两种方法质量矩阵存储量比较

2. 数值算例和讨论

2.1 谐振腔

2.2 滤波器

2.3 光子晶体

3. 结论

作者简介:
何泳成，heyongcheng@ihep.ac.cn