Tech-X, Boulder, Colorado
A 4D version of BNL's thin element spin tracking code SPINK* with limited elements has been successfully ported to a C++/GPU platform using GPULib**. This prototype used only quadrupoles, simple snakes, dipoles and drifts. It achieved an 80 fold speed up over serial CPU version of the code when pushing 100,000 particles. We present the approach used to track spin and orbit degrees of freedom of polarized proton beams simultaneously. We also present recent results of prototyping a general-purpose particle tracking on GPUs, discussing our CUDA implementation of maps for single-particle dynamics in the Argonne National Lab's accelerator code Elegant***.
*A.Luccio,"Spin tracking in RHIC code SPINK," Proc. of Adriatico Res Conf (1995)
**http://GPULib.txcorp.com
***M.Borland,"elegant:A Flexible SDDS-compliant Code for Accel. Sim.", APS LS-287 (2000)
GSI, Darmstadt
An effective analytical and semi-analytical method for internal electrical field calculations was proposed for ellipsoidal shaped beam as well as for a beam of arbitrary longitudinal shape with an elliptical transverse cross section. This method combines acceptable accuracy with a high speed of computation. The existing version of the DYNAMION code uses the particle-particle method to calculate the electrical field, which needs a significant time for computation. A semi-analytical algorithm of electrical field calculation was introduced into DYNAMION code. It allows much faster beam dynamics simulations than the old one (above 5·103 particles). The DYNAMION parameter "macroparticle size" was investigated in connection with the new space charge algorithm. The beam dynamics simulations were performed through the 1st Alvarez tank of the UNILAC using the original and the new method. The RMS emittance growth as a benchmark parameter shows sufficient agreement between both solvers.
Tsinghua University, Beijing
PSI, Villigen
TUB, Beijing
CIAE, Beijing
The 1.3 MW of beam power delivered by the PSI 590 MeV Ring Cyclotron together with stringent requirements regarding the controlled and uncontrolled beam losses poses great challenges with respect to predictive simulations. A new particle matter interaction model in OPAL is taking into account energy loss, multiple Coulomb scattering and large angle Rutherford scattering. This model together with the 3D space charge will significantly increase the predictive capabilities of OPAL. We describe a large scale simulation effort, which leads to a better quantitative understanding of the existing PSI high power proton cyclotron facility. The initial condition for the PSI Ring simulations is obtained from a new time structure measurements and the many profile monitors available in the 72 MeV injection line. A large turn separation and narrow beam size at the extraction turn is obtained. We show that OPAL can precise predict the radial beam pattern at extraction with large dynamic range (3-4 orders of magnitude). The described capabilities are mandatory in the design and operation of the next generation high power proton drivers.