A code for calculating the time evolution of beam parameters in high intensity circular accelerators

Beam emittances in a circular accelerator with a high beam intensity are strongly affected by the small angle intrabeam Coulomb scattering. In the computer simulation model we present here we used three coupled nonlinear differential equations to describe the evolution of the emittances in the transverse and the longitudinal planes. These equations include terms which take into account the intra-beam scattering, adiabatic damping, microwave instabilities, synchrotron damping, and quantum excitations. A code is generated to solve the equations numerically and incorporated into a FORTRAN code library. Circular high intensity physics routines are included in the library such as intrabeam scattering, Touschek scattering, and the bunch lengthening effect of higher harmonic cavities. The code runs presently in the PC environment. Description of the code and some examples are presented.