The authors propose and illustrate a general numerical method to follow the probability distribution in phase space as a function of time. It applies to any multiparticle system governed by Liouville, Vlasov or Vlasov-Fokker-Planck dynamics. The technique, based on discretization of the local Perron-Frobenius operator, is simple in concept, easy to implement, and numerically stable in examples studied to date. The authors illustrate by treating longitudinal dynamics in electron storage rings with realistic wake field. Applied to the SLC damping rings, the method gives the observed current threshold for bunch lengthening, and several aspects of observed behavior above threshold, including the presence of a bursting or sawtooth mode. In contrast to previous particle-in-cell simulations, the authors have very low numerical noise and the ability to follow the motion over several damping times. The method has also been applied to the coherent beam-beam interaction. It appears likely that this approach will be of interest for some of the central problems of this workshop, for instance matching of space-charge dominated beams to a focusing channel, and coherent synchrotron radiation with self-consistent charge/current density.
[1]
Alexander W. Chao,et al.
Physics Of Collective Beam Instabilities In High Energy Accelerators
,
1993
.
[2]
K. Hepp,et al.
Nonlinear dynamics aspects of particle accelerators
,
1986
.
[3]
MATTHEW SANDS,et al.
PHYSICS OF ELECTRON STORAGE RINGS. AN INTRODUCTION.
,
2018
.
[4]
M. Tigner,et al.
The physics of particle accelerators
,
1984
.
[5]
Andrew G. Glen,et al.
APPL
,
2001
.
[6]
S. Zacks,et al.
Introduction to stochastic differential equations
,
1988
.
[7]
Christian Soize,et al.
The Fokker-Planck Equation for Stochastic Dynamical Systems and Its Explicit Steady State Solutions
,
1994,
Series on Advances in Mathematics for Applied Sciences.