Differential algebraic methods for space charge modeling and applications to the University of Maryland Electron Ring
The future of particle accelerators is moving towards the intensity frontier; the need to place more particles into a smaller space is a common requirement of nearly all applications of particle accelerators. Putting large numbers of particles in a small space means that the mutual repulsion of these charged particles becomes a significant factor, this effect is called space charge. In this work we develop a series of differential algebra based methods to simulate the effects of space charge in particle accelerators. These methods were used to model the University of Maryland Electron Ring, a small 3.8 meter diameter 10 KeV electron storage ring designed to observe the effects of space charge in a safe, cost effective manner. The methods developed here are designed to not only simulate the effects of space charge on the motions of the test particles in the system but to add their effects to the transfer map of the system. Once they have been added useful information about the beam, such as tune shifts and chromaticities, can be extracted directly from the map. In order to make the simulation self consistent, the statistical moments of the distribution are used to create a self consistent Taylor series representing the distribution function, which is combined with pre-stored integrals solved using a Duffy transformation to find the potential. This method can not only find the map of the system, but also advance the particles under most conditions. For conditions where it cannot be used to accurately advance the particles a differential algebra based fast multipole method is implemented. By using differential algebras to create local expansions, noticeable time savings are found.