Efficient Numerical and Parallel Methods for Beam Dynamics Simulations

TitleEfficient Numerical and Parallel Methods for Beam Dynamics Simulations
Publication TypeJournal Article
Year of Publication2010
AuthorsXu, J, Mustapha, B, Min, MS
Date Published10/2010
Other NumbersANL/MCS-P1803-1010

This paper describes various numerical and parallel methods used at Argonne National Laboratory during the past five years to develop scalable software packages for beam dynamics simulations. It is based on the paper \"Developing petascalable algorithms for beam dynamics simulations,\" which appeared in the proceedings of the first International Particle Accelerator Conference in 2010 [Jin Xu et al. Proceedings of IPAC10, Kyoto, Japan, May 23-28 (2010)]. Our focus here is on the standard particle-in-cell (PIC) method, direct Vlasov solvers, and scalable Poisson solvers. Among these methods the most challenging are scalable Poisson solvers, designed in three dimensions to account for space charge effects. Several approaches have been used to solve Poisson\'s equation efficiently in different situations. High-order numerical methods have been adopted to increase the accuracy. Domain decomposition methods have been used for parallelizing the solvers, and good scaling has been achieved. Preliminary results for the direct Vlasov solvers have been obtained for up to four dimensions. The parallel beam dynamics code PTRACK, which uses the PIC method to solve Poisson\'s equation, has been used as a workhorse for end-to-end simulations and large-scale design optimization for linear accelerators. We have successfully run the parallel methods on tens of thousands of processors of the IBM Blue Gene/P system at the Argonne Leadership Computing Facility, which provides an excellent environment for conducting large-scale beam dynamics simulations.