Mathematical modeling

 

 

Mathematical simulation and analysis of the cyclotron U-120M

 

Main modelling systems:


Matrix model of cyclotron

  • computation and optimization of cyclotron parameters for arbitrary selected accelerating regime set-up 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:

  • 2-D preliminary magnetic field topography measurement and processing

  • computation of 3-D 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.6x4mm2

  • 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 U-120M