Mathematical simulation and analysis of the cyclotron U120M
Main modelling systems:
Matrix model of cyclotron
 computation and optimization of cyclotron parameters for arbitrary selected accelerating regime setup via the computer control system
 numerical and graphical representation of calculated cyclotron parameters and main beam properties
 saving the data files for dynamic simulation of any accelerating regime.
Dynamic simulation of acceleration and extraction process
 equilibrium orbit analyses
 computation of ion trajectories in 3D space ( injection, acceleration and extraction )
 numerical and graphical representation of calculated beam properties

Auxiliary systems:
 2D preliminary magnetic field topography measurement and processing
 computation of 3D electric field of the actual accelerating structure
 interpolation of magnetic field topography for any accelerating regime

Matrix model of cyclotron
 Main goal: calculations of currents in the main and trim magnetic coils for any accelerated regime.
 Acquirement of a radial dependence of an average magnetic field as near as possible to the isochronous one. Optimization of beam phase radial dependencies as well as of the radial and azimutal free oscillations and the trim coils power consumption.
 The magnet currents are optimized by an automatically working algorithm. It is possible to minimize either the magnetic field deviations from isochronism or the beam phase history.
 A fine manual tuning with respect to the beam phase history and the frequency of radial and axial free oscillations is also possible.
 Graphical and numerical representation of acquired data. Saving the files with description of the corresponding magnetic field for an use in dynamic simulation programs.

Simplified scheme of matrix model of cyclotron
Input and otput windows of matrix model
Processing of Cyclotron magnetic field
Electrostatic field
Equilibrium orbits analyses
Ion beam dynamics computation
Simplified scheme of beam dynamics mathematic simulation
Setting up the accelerated beam parameters and program modification
Notes:
 Ions injected from an ion source are distributed regularly at the slit 1.6 x 4mm^{2}
 The contribution of each injected ion to total “beam current” depends on its starting position and the value of Dee voltage at the moment of injection (i.e. on the value of RF phase at that moment)
 It is also possible to inject ions at any radius with parameters obtained for example from equilibrium orbits
 All components of magnetic and electric field are calculated at each integration step at all region of acceleration. The loss of ion is then indicated by zero or maximal potential at actual position, corresponding either to the high voltage or to the grounded electrode
 The coordinate axes or the beam probes can be chosen as the beam stop lines. It is also possible to prepare the conditions for an automatic beam extraction to required point or for a fine manual tuning of extracted beam

Accelerated beam properties
Centering of accelerated beam by harmonic coils
Influence of centering on the beam (Rstop = 48.00 cm)
Numerical output of accelerated beam properties
*** ACCELERATED BEAM: *** Accelerated 215 from 255 ions in regime Hm_34_6b.bhp :
Phase from 40o to 40o with step 5.0o Zet from 2.0 to +2.0 with step 1.00mm drIZ from 0.80 to 0.80 with step 0.40mm Starting energy of ions 0.000050 MeV No dFi0
Rstop= 48.20 cm, fiStop= 134.00o Rorb= 48.34 cm, Xc= 0.19 cm, Yc= 0.09 cm. Current at 1. orbit= 1712.4, at RStop= 1693.0. Beam transfer from first orbit to RMax is 98.9 %.
Average values of accelerated ions properties : Rmean= 48.24cm, radial. angle= 2.01o. Zmean= 0.00 cm, axial. angle= 0.00o. Emean= 32.240 MeV.

Analyses of the accelerated ions energy distribution
Accelerated beam emittance at the stripping foil
Stripping method of H^{}/D^{} extraction
Setting up the extraction system parameters
The magnetic field of the correction magnet
Extraction of the 32.25MeV proton beam to the ion line
Numerical output of extracted beam properties
Distribution of ions extracted to the triplet entry
Extracted beam energy: calculation & foil activation
Extracted beam emittance at the triplet entry
Acceleration and extraction of H^{} and D^{} ions at cyclotron U120M