Publications Details
New field solver for the buckshot code
The Buckshot code is a gridless particle simulation code which is used extensively at Sandia to study relativistic electron beam propagation in the ion-focused regime. The present version of the code evaluates the force on a particle by summing over all the other particles, thus the execution time is proportional to N{sup 2} where N is the number of simulation particles. A new gridless field solver has been developed with time scaling of N log N and the breakeven point between the old and new code is around N = 64 on the Cray X-MP computer. For N = 1000, the new solver is about nine times faster than the old version. The new solver is based on a solution of the two-dimensional Poisson equation which if Fourier decomposed in the azimuthal direction and the solution of the radial equation is represented by integrals over the charge density. These integrals are then replaced by sums over the simulation particles which are assumed to be point particles. The near-field singularity is removed by the Fourier decomposition so long as the number of Fourier modes is much less than the number of simulation particles. The algorithm is written in such a way that the field due to a given species is Fourier decomposed with respect to the center of mass of that species, thus it is possible to study linear and nonlinear ion hose physics with a very small number of azimuthal modes. Typically M = 2 is found to be sufficient for most IFR problems. The old and new solvers have been compared and the agreement is excellent. 2 refs., 3 figs.