CN101695214B - 等时性回旋加速器提高轴向聚焦力的磁极非对称镶条方法 - Google Patents
等时性回旋加速器提高轴向聚焦力的磁极非对称镶条方法 Download PDFInfo
- Publication number
- CN101695214B CN101695214B CN2009102111563A CN200910211156A CN101695214B CN 101695214 B CN101695214 B CN 101695214B CN 2009102111563 A CN2009102111563 A CN 2009102111563A CN 200910211156 A CN200910211156 A CN 200910211156A CN 101695214 B CN101695214 B CN 101695214B
- Authority
- CN
- China
- Prior art keywords
- isochronous
- magnetic field
- panel
- magnetic pole
- cyclotron
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 230000005291 magnetic effect Effects 0.000 title claims abstract description 74
- 238000000034 method Methods 0.000 title claims abstract description 21
- 230000010355 oscillation Effects 0.000 claims description 24
- 239000002245 particle Substances 0.000 claims description 19
- 235000011888 snacks Nutrition 0.000 claims description 19
- 238000005516 engineering process Methods 0.000 abstract description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 241000065675 Cyclops Species 0.000 description 2
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 229920000535 Tan II Polymers 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005266 casting Methods 0.000 description 1
- 230000019643 circumnutation Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005684 electric field Effects 0.000 description 1
- 239000003302 ferromagnetic material Substances 0.000 description 1
- 238000005242 forging Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 229910052742 iron Inorganic materials 0.000 description 1
- 238000003754 machining Methods 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 238000004804 winding Methods 0.000 description 1
Images
Landscapes
- Particle Accelerators (AREA)
Abstract
本发明涉及回旋加速器的磁铁设计技术,具体涉及一种等时性回旋加速器提高轴向聚焦力的磁极非对称镶条方法。该方法通过有限元模型中镶条任意半径处需要增加的螺旋角,得到主磁铁有限元模型所采用的极坐标系下镶条的轮廓坐标,根据镶条的轮廓坐标加工镶条,进行等时性磁场的垫补。这种磁极非对称镶条方法可以提高等时性回旋加速器在大半径区域的轴向聚焦力,有效防止非线性共振产生。
Description
技术领域
本发明涉及回旋加速器的磁铁设计技术,具体涉及一种等时性回旋加速器提高轴向聚焦力的磁极非对称镶条方法。
背景技术
回旋加速器是一种利用磁场使带电粒子作回旋运动,在运动中经高频电场反复加速的装置,其基本结构包括磁铁系统、离子源、注入系统和引出系统。其中,磁铁系统包括磁铁和调节线圈,使回旋加速器能够获得等时性磁场。磁场的垫补是使得回旋加速器中心平面磁场满足设计要求的重要手段,特别是对于回旋加速器来说,磁场的垫补是无法避免的,因为磁铁材料的不均匀性,磁铁加工和安装等因素将引起中心平面的磁场偏离设计的磁场分布。
在径向扇型等时性回旋加速器中,离子在轴向(即垂直于离子旋转平面的方向)总的受力方程为: 式中m为离子质量,ωc为离子回旋频率,为沿方位角平均磁场强度,f是磁场沿方位角的调变度,ωz表示轴向振荡角频率,r、z分别为粒子径向与轴向的坐标。由于ωz 2>0时,公式(1)才有振荡解,否则解发散,束流会损失,因此,必须有而且ωz 2越大,轴向聚焦力越大。轴向自由振荡频率vz=ωz/ωc,类似的,由离子在径向总的受力方程可以得到径向自由振荡频率vr=ωr/ωc。若轴向聚焦力过小,当vz减小至vr/2附近就会出现Walkinshaw共振。
对于螺旋扇型等时性回旋加速器,其中ξ为螺旋角。可见由于螺旋角的作用,螺旋扇型等时性回旋加速器的轴向聚焦力增大了,这就可以很好的避免Walkinshaw共振的出现。但采用螺旋型磁极,不仅增加了主磁铁本身的设计和建造难度,也使得高频谐振腔、径向靶等束流诊断设备、剥离引出系统等的工程复杂程度明显增大。
在使用径向扇形磁极的等时性回旋加速器主磁铁的工程建设过程中,用于磁极和盖板、磁轭等主体部件的大型铸锻件的磁性能,或加工过程中的公差控制出现严重超差,均有可能出现轴向聚焦力不够而导致非线性共振,尤其是靠近引出区的大半径区出现Walkinshaw等非线性共振会导致束流包络长大而损失。对于使用径向扇形磁极的回旋加速器,镶条垫补是主要的等时性磁场垫补方式之一。镶条一般是由与主磁铁磁极材料相同的材料制成的配铁,通常以可拆卸的方式设置于主磁铁磁极的两边,镶条垫补方法主要是利用镶条比主磁铁磁极易于加工的特点,通过改变镶条的外形尺寸来垫补获得等时性磁场。
发明内容
本发明的目的在于针对目前径向扇型等时性回旋加速器的缺陷,结合螺旋扇型等时性回旋加速器中螺旋型磁极的特点,提供一种等时性回旋加速器提高轴向聚焦力的磁极非对称镶条方法,从而有效的防止径向扇型等时性回旋加速器非线性共振的产生。
本发明的技术方案如下:一种等时性回旋加速器提高轴向聚焦力的磁极非对称镶条方法,包括如下步骤:
(1)根据回旋加速器当前磁场的数据,利用粒子跟踪程序输出的轴向自由振荡频率vz0,设计一条平滑且自由振荡频率值增大了的粒子轨道半径的曲线vzideal;
(2)利用轴向自由振荡频率公式得到磁极任意半径处的如下关系式:vzidieal 2-vz0 2=F2(tan2(ξ0+ξreal)-tan2ξ0),通过该关系式解出有限元模型中磁极任意半径处需要增加的螺旋角ξreal;其中,ε为粒子总能量,ε0为粒子静止能量,ξ为螺旋角,F2为对称垫补至等时性磁场时,粒子跟踪程序计算得到的调变度,ξ0为等时性磁场时粒子跟踪程序计算出轴向振荡频率vz后,根据轴向自由振荡频率公式得到的等效螺旋角;
(3)根据步骤(2)中得到的有限元模型中磁极任意半径处需要增加的螺旋角ξreal,得到主磁铁有限元模型所采用的极坐标系下镶条的轮廓坐标;
(4)根据镶条的轮廓坐标加工镶条,进行等时性磁场的垫补。
进一步,如上所述的等时性回旋加速器提高轴向聚焦力的磁极非对称镶条方法,步骤(1)中所述的回旋加速器当前磁场为设计阶段电磁场计算程序计算得到的磁场,或者实际测磁测得的待垫补磁场。
进一步,如上所述的等时性回旋加速器提高轴向聚焦力的磁极非对称镶条方法,步骤(3)中得到主磁铁有限元模型所采用的极坐标系下镶条的轮廓坐标的具体方式如下:
设有限元模型中磁极任意半径处需要增加的螺旋角其中代表幅角为θn的单位复数,(rn,θn)、(rn+1,θn+1)为镶条相邻两个控制点在极坐标下的坐标,根据等时性磁场对应的轴向和径向自由振荡频率曲线图,以及镶条与高频边界来选定rn和rn+1,
则,
得到
选定θn+1,解出xn=cos(θn+1-θn)后即可定出θn,然后再以θn与rn-1处的螺旋角作为已知条件解出θn-1,依次递推解出主磁铁有限元模型所采用的极坐标系下整个镶条的轮廓坐标。
本发明的有益效果如下:本发明在径向扇型磁极设计的基础上,融入了螺旋扇型磁极的特点,利用螺旋扇型磁极螺旋角的技术特征,导出径向扇型磁极两侧镶条的轮廓坐标,从而在径向扇型磁极设计的基础上增加额外的轴向聚焦力,特别是提高了等时性回旋加速器在大半径区域的轴向聚焦力,有效防止非线性共振产生。
附图说明
图1为本发明的方法流程图;
图2为主磁铁有限元模型所采用的极坐标系下镶条的结构示意图;
图3为对称垫补及不对称垫补至等时性磁场对应的自由振荡频率比较示意图。
具体实施方式
下面结合附图和实施例对本发明进行详细的描述。
为了在直边扇型磁极设计的基础上增加额外的轴向聚焦力,本发明提出了如下的非对称镶条方法,具体步骤如图1所示。
首先,从现有的磁场(可以是设计阶段电磁场计算程序计算得到的磁场,也可以是实际测磁测得的待垫补磁场)出发,利用粒子跟踪程序输出的轴向自由振荡频率vz0,设计出一条平滑且自由振荡频率值增大了的目标vzideal曲线(vzideal-粒子轨道半径的曲线)。在规划vzideal曲线时,应当使得vzideal与vr曲线避免回旋加速器中广知的有害共振。
然后,利用下面的轴向自由振荡频率公式:
得到磁极任意半径处的如下关系式:
vzidieal 2-vz0 2=F2(tan2(ξ0+ξreal)-tan2ξ0) (2-2)
式中,ε为粒子总能量,ε0为粒子静止能量,ξ为螺旋角,F2为对称垫补至等时性磁场时,粒子跟踪程序计算得到的调变度;ξ0为等时性磁场时粒子跟踪程序计算出轴向振荡频率vz后,根据式(2-1)计算而输出的等效螺旋角,利用(2-2)式解出实际有限元模型中需要增加的螺旋角ξreal。
接下来,由ξreal导出主磁铁极坐标系下的镶条的轮廓坐标。考察粒子转出磁极的一边的镶条轮廓,由螺旋角定义提供聚焦力的分析,该侧为散焦区,镶条角度随半径增大。由图2中所示的主磁铁有限元模型所采用的极坐标系下镶条的结构,可知这里的vn即为有限元模型中不同半径处的螺旋角ξreal,其中代表幅角为θn的单位复数,(rn,θn)、(rn+1,θn+1)为镶条相邻两个控制点在极坐标下的坐标。
故,
于是,选定θn+1,解出xn=cos(θn+1-θn)后即可定出θn,然后再以θn与rn-1处的螺旋角作为已知条件解出θn-1,依次递推解出主磁铁有限元模型所采用的极坐标系下镶条的轮廓坐标
根据上述分析,采用非对称镶条垫补方法得到的轴向振荡频率的增量为:Δvz=vzideal-vz0 (2-5)
下面介绍一个具体的实例。某中能回旋加速器的为紧凑型回旋加速器,有上下各四个磁极、四个磁轭和上下盖板组成,四个磁极对称分布在盖板上,磁极的半径均为2000mm,直边扇形结构,磁极的扇角为47°;四个谷区的深度为785mm,谷区的角宽度为43°,在其主磁铁的90度模型中,利用粒子跟踪程序CYCLOP得到对称垫补至等时性磁场对应的vr/2、vz,如图3中虚线与点划线所示。由图3可以看出从半径R=155cm附近开始,vz与vr/2的差值开始变小。由于一般有vr∝γ,故为了避免vz=vr/2的Walkinshaw共振,vz应当增大,即需要加大轴向聚焦力,因此不对称镶条垫补从R=183cm开始到R=155cm。利用非对称镶条方法,并根据镶条与高频边界限制,初始输入条件为:R=183cm前一个镶条控制点极坐标(190,24.94953°),R=183cm对应螺旋角ξ=1.6°设根据公式(2-3)、(2-4),设计了一组镶条不同半径对应的螺旋角,如表1所示。
垫补至等时性磁场后,利用Cyclop程序计算对应的径向与轴向自由振荡频率如图3双点划线与实线所示。由图3中对称镶条与非对称镶条均垫补至等时场时的振荡频率比较可以看出,非对称镶条方法可以明显提高轴向振荡频率,提高轴向聚焦,避免共振。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其同等技术的范围之内,则本发明也意图包含这些改动和变型在内。
Claims (3)
1.一种等时性回旋加速器提高轴向聚焦力的磁极非对称镶条方法,包括如下步骤:
(1)根据回旋加速器当前磁场的数据,利用粒子跟踪程序输出的轴向自由振荡频率vz0,设计一条平滑且自由振荡频率值增大了的粒子轨道半径的曲线vzideal;
(2)利用轴向自由振荡频率公式 得到磁极任意半径处的如下关系式:vzidieal 2-vz0 2=F2(tan2(ξ0+ξreal)-tan2ξ0),通过该关系式解出有限元模型中镶条任意半径处需要增加的螺旋角ξreal;其中,ε为粒子总能量,ε0为粒子静止能量,ξ为螺旋角,F2为对称垫补至等时性磁场时,粒子跟踪程序计算得到的调变度,ξ0为等时性磁场时粒子跟踪程序计算出轴向振荡频率vz后,根据轴向自由振荡频率公式得到的等效螺旋角;
(3)根据步骤(2)中得到的有限元模型中磁极任意半径处需要增加的螺旋角ξreal,得到主磁铁有限元模型所采用的极坐标系下镶条的轮廓坐标;
(4)根据镶条的轮廓坐标加工镶条,进行等时性磁场的垫补。
2.如权利要求1所述的等时性回旋加速器提高轴向聚焦力的磁极非对称镶条方法,其特征在于:步骤(1)中所述的回旋加速器当前磁场为设计阶段电磁场计算程序计算得到的磁场,或者实际测磁测得的待垫补磁场。
3.如权利要求1或2所述的等时性回旋加速器提高轴向聚焦力的磁极非对称镶条方法,其特征在于:步骤(3)中得到主磁铁有限元模型所采用的极坐标系下镶条的轮廓坐标的具体方式如下:
设有限元模型中镶条任意半径处需要增加的螺旋角 其中 代表幅角为θn的单位复数,(rn,θn)、(rn+1,θn+1)为镶条相邻两个控制点在极坐标下的坐标,根据等时性磁场对应的轴向和 径向自由振荡频率曲线图,以及镶条与高频边界来选定rn和rn+1,
选定θn+1,解出xn=cos(θn+1-θn)后即可定出θn,然后再以θn与rn-1处的螺旋角作为已知条件解出θn-1,依次递推解出主磁铁有限元模型所采用的极坐标系下整个镶条的轮廓坐标。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009102111563A CN101695214B (zh) | 2009-11-06 | 2009-11-06 | 等时性回旋加速器提高轴向聚焦力的磁极非对称镶条方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009102111563A CN101695214B (zh) | 2009-11-06 | 2009-11-06 | 等时性回旋加速器提高轴向聚焦力的磁极非对称镶条方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101695214A CN101695214A (zh) | 2010-04-14 |
CN101695214B true CN101695214B (zh) | 2011-12-07 |
Family
ID=42094135
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2009102111563A Active CN101695214B (zh) | 2009-11-06 | 2009-11-06 | 等时性回旋加速器提高轴向聚焦力的磁极非对称镶条方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101695214B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106211540B (zh) * | 2016-07-29 | 2018-10-09 | 中国原子能科学研究院 | 230MeV超导回旋加速器防止引出区有害共振的机械结构 |
CN106804091B (zh) * | 2017-03-21 | 2018-04-10 | 华中科技大学 | 一种回旋加速器等时性磁场垫补方法及系统 |
CN107466152B (zh) * | 2017-07-25 | 2019-09-20 | 中国科学院高能物理研究所 | 多极磁铁及其磁场谐波垫补方法 |
CN110740561B (zh) * | 2018-07-21 | 2021-07-30 | 中国原子能科学研究院 | 超导回旋加速器中引出磁通道元件的边缘场垫补方法 |
CN109526132B (zh) * | 2018-11-27 | 2019-10-15 | 中国原子能科学研究院 | 一种异型大尺寸高温超导磁铁 |
-
2009
- 2009-11-06 CN CN2009102111563A patent/CN101695214B/zh active Active
Non-Patent Citations (5)
Title |
---|
张天爵等.100MeV回旋加速器主磁铁设计与工艺技术研究.《Chinese Physics C(HEP & NP)》.2008,第32卷(第增刊I期),1-3. * |
张天爵等.紧凑型回旋加速器主磁铁一次谐波垫补方法的研究与实践.《中国原子能科学研究院年报2006》.2007,44-46. * |
王川等.CYCIAE-100 的非对称镶条对提高轴向聚焦力的影响.《中国原子能科学研究院年报2007》.2008,105-106. * |
钟俊晴.100MeV回旋加速器中心区实验台架的磁场测量和垫补.《中国优秀硕士学位论文全文数据库》.2008,(第4期),1-72. * |
钟俊晴等.强流回旋加速器综合试验装置的磁场测量与垫补.《Chinese Physics C(HEP & NP)》.2008,第32卷(第增刊I期),56-58. * |
Also Published As
Publication number | Publication date |
---|---|
CN101695214A (zh) | 2010-04-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101695214B (zh) | 等时性回旋加速器提高轴向聚焦力的磁极非对称镶条方法 | |
CN102819033B (zh) | 一种离子束流自动测量系统及测量方法 | |
CN105245146B (zh) | 一种实现永磁同步电动机转速单向混沌运动的方法 | |
CN106411213A (zh) | 用于提供电机定子电流值的方法及装置 | |
Li et al. | Microscopic analysis of order parameters in nuclear quantum phase transitions | |
CN101697659A (zh) | 回旋加速器采用连续曲面的等时性磁场精密垫补方法 | |
CN101697658A (zh) | 中等能量回旋加速器强聚焦直边扇形磁铁 | |
CN101827489B (zh) | 用于加速负氢、h2+的紧凑型回旋加速器 | |
Shao et al. | Design and optimisation of energy‐efficient PM‐assisted synchronous reluctance machines for electric vehicles | |
CN105978434B (zh) | 一种永磁同步电机最大转矩电流比控制方法及控制系统 | |
CN103983927A (zh) | 根据霍尔推力器中耦合振荡伴生的动态磁场确定线圈安匝变化百分率范围的方法 | |
CN104362917A (zh) | 车用交流发电机漏磁问题的优化设计方法 | |
Tanigaki et al. | Present status of FFAG accelerators in KURRI for ADS study | |
Wang | Magnetic field distribution of injection chicane dipoles in Spallation Neutron Source accumulator ring | |
Karamysheva et al. | Magnetic system for SC200 superconducting cyclotron for proton therapy | |
Shin et al. | Measurement and analysis of a 13 MeV cyclotron magnetic field | |
CN112541154A (zh) | 一种磁路功率的计算方法 | |
Li et al. | Benchmark simulations of gyro-kinetic electron and fully-kinetic ion model for lower hybrid waves in linear region | |
CN113791345B (zh) | 一种交流铜耗计算方法 | |
Namgoong et al. | Design of cyclotron magnet with tenth harmonics for accelerator mass spectrometry | |
CN103338579B (zh) | 一种条带型横向四极分量激励器 | |
CN103440432B (zh) | 一种利用瞬变等离子体板新增频率点的方法 | |
Xing et al. | RF-induced frequency-shift resistant design of the resonant cavity of the radio frequency quadrupole with the high average-power operation | |
Harrison | Resonant extraction at the tevatron | |
Park et al. | Design Features Of K= 100 Cyclotron Magnet For ISOL RIB Production |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |