1 INTRODUCTION
The SC200 proton therapy superconducting cyclotron will be able to accelerate protons to an energy of 200 MeV with a maximum beam intensity of 1 μA. The average magnetic field of the cyclotron is up to 3.6 T [1-2]. As a key component of the cyclotron, the superconducting magnet provides the isochronous magnetic field with a strong focusing force to restrict beam movement. However, the asymmetry of sectors and the superconducting coils (SCs) generates a radial component, Br, of the magnetic field in the median plane, directly affecting the vertical deviation of the beam orbit from the median plane [3-6].
A number of measurement methods have been proposed to measure the average radial component of the magnetic field (Brav), including nuclear magnetic resonance (NMR) [7], hall probes [8-10], and flux measurements with search coils [3][11-14]. To design such a measurement system, the field strength, homogeneity, variation in time, and required accuracy should be considered8. Search coils are highly sensitive and less strongly influenced by the vertical component of the magnetic field (Bz), which is more accurate than the hall probes for magnetic field measurement. The use of a variety of search coils is the most-often used tool for characterizing field qualities in cyclotrons [15-18].
For SC200, the Brav measurement system was designed to measure the Brav of the magnetic field by using search coils with a high testing accuracy. The Brav measurement system mainly includes a drive system, control system, and data acquisition system. In particular, based on the electromagnetic induction principle, there are 15 search coils moving in the magnetic field to measure the Brav. The minimum and maximum measurement radii were 30 mm and 620 mm, respectively. The search coils are made of copper and covered with an insulating varnish. Owing to the quality differences in the winding process, it is necessary to perform the commissioning of the search coil outside of the SC200 cyclotron.
In this study, we present develop a detailed design and investigate the required commissioning work to improve Brav measurement accuracy. This paper describes the Brav tolerance analysis, the commissioning platform design, and research. The analysis and results are of great significance for the Br measurement of the SC200 proton therapy superconducting cyclotron.
2 B_rav TOLERANCES FOR SC200 SUPERCONDUCTING CYCLOTRON
Owing to the asymmetry of sectors and SCs, the Brav existing in the Azimuth-Varied Field cyclotron, which causes the effective median plane (EMP) of the magnetic system, does not coincide with the median plane of the cyclotron. The EMP was formulated by J. I. M. Botman and H. L. Hagedorn for the cyclotron work region [3]. Thus, the tolerances for the Brav of the magnetic field should be estimated.
When the magnetic system symmetry is broken, stable vertical oscillations exist near the EMP. The horizontal field component exists in the median plane. To separately identify the tolerances for the horizontal field components from the permitted vertical beam offset, the Brav can be written as
where Zeff is the vertical position of the EMP, Bzav is the average vertical component of the magnetic field at radius R, Qz represents the vertical betatron frequency, and R is the magnetic field radius.
The SC200 superconducting cyclotron has been proposed and designed [17,19-21]. The isochronous field Bzav distribution is shown in Fig.1(a), and the maximum magnetic field was approximately 3.54 T [20,22]. From the computed vertical betatron frequency Qz in Fig.1(b), Eq. (1) shows the tolerance distribution of Brav for the acceptable vertical offset of the beam of 1 mm shown in FIG.2.
-201906/1001-8042-30-06-007/alternativeImage/1001-8042-30-06-007-F001.jpg)
-201906/1001-8042-30-06-007/alternativeImage/1001-8042-30-06-007-F002.jpg)
With the frequency Qz close to 0.32 in the extraction region, the isochronous field Bzav was approximately 3.54 T. The Brav tolerance in the median plane was approximately 3–7 G, shown in Fig.2. Consequently, the testing resolution of the Brav measurement system should be smaller than 1 G.
3 COMMISSIONING PLATFORM DESIGN
3.1 TESTING PRINCIPLE
The magnetic flux distribution produced by the test coil had an excitation current of 2 A, as shown in Fig.3. It generates a magnetic field of 3–7 G, the same as the magnetic field Brav tolerance of SC200 superconducting cyclotron.
-201906/1001-8042-30-06-007/alternativeImage/1001-8042-30-06-007-F003.jpg)
A commissioning diagram comprising a test coil and 15 concentric search coils is shown in Fig.4. According to the electromagnetic induction principle, when the search coils move in the magnetic field, the induced voltage U is generated, which can be integrated by the fluxmeter. The testing flux change △Φt can be written as
-201906/1001-8042-30-06-007/alternativeImage/1001-8042-30-06-007-F004.jpg)
The radial component change △Br can be calculated using
The effective sensitive area of the search coil A is
where N is the number of search coil turns, R is the coil radius, and ΔZ is the distance shift from the median plane.
The search coil was first connected to the fluxmeter, which was shifted from its initial position of –Δz to a terminal position of +Δz. The fluxmeter acquired the magnetic flux generated by the coil. Subsequently, the other coils were changed to connect to the fluxmeter, and the previous steps were repeated.
3.2 PLATFORM DESIGN
Based on the testing principle, the search coils must be a priority during design. The search coils comprised copper wire, with a diameter of 0.3 mm, covered with insulating varnish. The search coils were numbered from C1 to C15. Each coil could be connected to the LINKJOIN LZ840 fluxmeter by a coil selector. Then, the coil parameters, including the radius, the number of turns, and the coil resistance, were sent to the fluxmeter. The distribution of the search coils on the disk is shown in FIG.5.
-201906/1001-8042-30-06-007/alternativeImage/1001-8042-30-06-007-F005.jpg)
The input coil parameters are shown in Table I.
Coil No. | C1 | C2 | C3 | C4 | C5 | C6 | C7 | C8 | C9 | C10 | C11 | C12 | C13 | C14 | C15 |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
Radius R (mm) | 30 | 55 | 80 | 105 | 130 | 155 | 180 | 205 | 230 | 255 | 280 | 305 | 330 | 355 | 380 |
Number of turns N | 34 | 26 | 20 | 16 | 10 | 10 | 8 | 8 | 8 | 6 | 6 | 6 | 6 | 6 | 6 |
Coil resistance(Ω) | 7.705 | 7.822 | 8.038 | 7.847 | 7.914 | 8.15 | 8.051 | 8.236 | 8.547 | 7.91 | 8.108 | 8.237 | 8.53 | 8.622 | 8.855 |
The commissioning platform of the Brav measurement system was designed with a drive system, control system, and data acquisition system, as shown in Fig.6. The self-aligning system includes two groups of cross screws and linear guide structures. These institutions can adjust the disk center coinciding with the geometric center of the superconducting cyclotron. The drive system can drive the search coils to move along the vertical axis of the cyclotron. The control system comprises the parameter settings area, the data and status display area, and the logic control area. The data acquisition system is responsible for collecting the measurement signals. The LINKJOIN LZ840 fluxmeter collects the flux signal with high sensitivity, which should be warmed for five minutes and done zero setting and drift setting before using.
-201906/1001-8042-30-06-007/alternativeImage/1001-8042-30-06-007-F006.jpg)
To avoid the magnetic field influence from Bz, it is necessary to adjust and confirm the relative position between the test coil and search coil. An angular error of 1 mrad will produce an approximately tenfold error in the Br value. Therefore, the parallelism of the search coil and test coil must be less than 0.4 mm, and that of the concentricity must be less than 0.5 mm. In this way, the Br value will not be influenced by the Bz value.
4 COMMISSIONING AND DISCUSSION
To confirm that the Brav measurement system operates well, based on the testing principle above, system commissioning should be conducted outside the cyclotron. The drive system made the search coils move from their initial position of -10 mm to the terminal position of +10 mm. The fluxmeter acquired the first flux Φ1 generated by the first search coil C1. Subsequently, the second flux Φ'1 was generated when the disk moved in reverse. Following this, C2 was switched to connect to the fluxmeter, and the previous steps were repeated. This was performed sequentially until coil C15, and all the measurement data were collected for further analysis.
We found that the coil speed affects the Brav results. Thus, the testing was conducted with five different speeds in the range of ±10 mm. The collected flux was fitted by the least square method. The Brav was calculated using Eq.3 in FIG.7, and the trend of the Brav curves are consistent. The theoretical value (Brth)was provided by the TOSCA software. The Brav was closer to the Brth as the coil speed increased from 1 mm/s to 2.5 mm/s. There were no obvious relative error changes between the 2.5 mm/s and 3 mm/s speeds.
-201906/1001-8042-30-06-007/alternativeImage/1001-8042-30-06-007-F007.jpg)
Compared with Brth, the relative error of Brav was calculated and is shown in FIG.8. The relative error of Brav with a speed of 2.5 mm/s was the smallest. The relative error of the curves is greater in two edge radius ranges. When increasing the coil speed, the flux variation within a unit time will increase correspondingly. Thus, the signal-to-noise ratio will be promoted. For the fluxmeter, the disturbance caused by noise will comprise a smaller proportion and the measurement results will be more accurate at the points where the value of the signal-to-noise ratio is larger. Otherwise, owing to the structure parameters of the fluxmeter, the original right signal completely becomes the interfering signal. In the case of a severe signal drift, the effective signal is often submerged. Therefore, 2.5 mm/s is considered to be the optimal speed for the search coil.
-201906/1001-8042-30-06-007/alternativeImage/1001-8042-30-06-007-F008.jpg)
Following this, the search coils were shifted four times with a speed of 2.5 mm/s. The Brav distributions were deduced and are shown in Fig.9. The Brth and Brav curves are consistent, as shown in Fig.10. The Brav denotes the average measured value from four Brav values.
-201906/1001-8042-30-06-007/alternativeImage/1001-8042-30-06-007-F009.jpg)
-201906/1001-8042-30-06-007/alternativeImage/1001-8042-30-06-007-F010.jpg)
As shown in Fig.11, the relative error curve first increases then decreases with the radius increase. The minimum error occurred at R = 205 mm, corresponding to C8. The minimum relative error was 0.09513% with a Brav of 7 G, which continued increasing to approximately 12.37% for Brav of 1 G at the edge point. In addition, the smallest absolute error and biggest absolute error of △Brav were 0.006857 G and 0.05402 G, respectively, both less than 1 G.
-201906/1001-8042-30-06-007/alternativeImage/1001-8042-30-06-007-F011.jpg)
The measurement accuracy of the Brav changed with the radial location. For the radial range of 0-205 mm, the magnetic field intensity increased together with the measuring accuracy. For the radial range of 205-380 mm, the measuring accuracy decreased. There is a low intensity at both ends of the Br distribution curve as the measurement accuracy is susceptible to interference. The measurement accuracy was high under the high intensity field in the middle radial position.
5 CONCLUSION
In summary, this paper presents the detailed design and commissioning for a Brav measurement system. The commissioning platform was designed to verify the measurement precision based on the electromagnetic induction principle. Under the testing magnetic field, the coil speed of 2.5 mm/s is considered to be the optimal speed for the search coil. The relative error was approximately 0.1% under the maximum Br of 7 G. The measurement precision reached 1.0×10-3, which satisfies the measurement requirements.
The analysis and results are of great significance for the Br measurement of the SC200 proton therapy superconducting cyclotron. However, the measurement error of the system exceeds 10% in edge radius ranges; thus, in future, we must improve the measurement precision by optimizing the system structure or researching a reasonable optimization method.
Magnetic system for SC200 superconducting cyclotron for proton therapy
,Conceptual Design of Scanning Magnets for Superconducting Proton Therapy Facility
[J]. IEEE Transactions on Applied Superconductivity, 27, 1-8(2017). doi: 10.1109/TASC.2017.2698201Median plane effects in the Eindhoven AVF cyclotron
. IEEE Trans. On Nucl. Sci. 28, 2128-2130(1981). doi: 10.1109/TNS.1981.4331613Vertical beam motion in the AGOR cyclotron
.The axial injection system of the isochronous cyclotron
[J]. 1513-1515(1992). http://accelconf.web.cern.ch/AccelConf/e92/PDF/EPAC1992_1513.PDFHarmonic coil application for shaping of the magnetic field of isochronous cyclotron AIC-144
[J], (2000). http://www1.jinr.ru/Preprints/2000/P9200039.pdfMapping of the New IBA Superconducting SYNCHROCYCLOTRON (S2C2) for Proton Therapy
,Field mapping system for cyclotron magnet
[J]. Nuclear Inst & Methods in Physics Research A, 545, 533-541(2005). doi: 10.1016/j.nima.2005.02.009Precise integration of undulator Hall probe scans
[J]. IEEE Transactions on Magnetics Mag, 30, 2435-2438(1994). doi: 10.1109/20.305769The Magnetic Field Mapping System for the IBA C70 Cyclotron
[J]. Cyclotrons and Their Applications (2007). http://accelconf.web.cern.ch/accelconf/c07/PAPERS/78.PDFOverview of magnetic measurement techniques, CERN Academic Training Program
(2003). http://cds.cern.ch/record/610089Field Mapping System Design for the Superconducting Cyclotron CYCIAE-230
[J]. IEEE Transactions on Applied Superconductivity, 28, 1-1(2018). doi: 10.1109/TASC.2018.2789447Magnetic field measurements and mapping techniques
[J]. CERN Academic Training Program (2003). http://cds.cern.ch/record/610089?ln=zh_CN.Classification of magnetic measurement methods
[J], CERN, Geneva, Switzerland, 16, 70-83(1992). doi: 10.5170/CERN-1992-005.70Characteristics and performance of the system developed for magnetic mapping of the NSCL superconducting K800 cyclotron magnet
[J], IEEE Trans. On Nucl. Sci. 32, 3734 (1985). doi: 10.1109/tns.1985.4334485Development of a method for measuring the radial component of the magnetic field in AVF cyclotrons
,Median plane effects and measurement method for radial component of magnetic field in AVF cyclotrons, Magnet and Vacuum
:Median plane magnetic field mapping for superconducting cyclotron (SCC) in VECC
[J].Computer modeling of magnet for SC200 superconducting cyclotron
,Study of the beam extraction from superconducting cyclotron SC200
[C],Compact superconducting cyclotron SC200 for proton therapy
,