×

A novel adaptive time stepping variant of the Boris-Buneman integrator for the simulation of particle accelerators with space charge. (English) Zbl 1351.81023

Summary: We show that adaptive time stepping in particle accelerator simulation is an enhancement for certain problems. The new algorithm has been implemented in the OPAL (Object Oriented Parallel Accelerator Library) framework. The idea is to adjust the frequency of costly self-field calculations, which are needed to model Coulomb interaction (space charge) effects. In analogy to a Kepler orbit simulation that requires a higher time step resolution at the close encounter, we propose to choose the time step based on the magnitude of the space charge forces. Inspired by geometric integration techniques, our algorithm chooses the time step proportional to a function of the current phase space state instead of calculating a local error estimate like a conventional adaptive procedure. Building on recent work, a more profound argument is given on how exactly the time step should be chosen. An intermediate algorithm, initially built to allow a clearer analysis by introducing separate time steps for external field and self-field integration, turned out to be useful by its own, for a large class of problems.

MSC:

81-08 Computational methods for problems pertaining to quantum theory
81Vxx Applications of quantum theory to specific physical systems

Software:

OPAL
PDFBibTeX XMLCite
Full Text: DOI arXiv

References:

[1] Adelmann, A., The OPAL (Object Oriented Parallel Accelerator Library) Framework (2014), Paul Scherrer Institut, Technical Report PSI-PR-08-02
[2] Agarwalla, S. K.; Conrad, J. M.; Shaevitz, M., Short-baseline neutrino oscillation waves in ultra-large liquid scintillator detectors (2011)
[3] Abs, M.; Adelmann, A., Multimegawatt DAE \(δ\) ALUS cyclotrons for Neutrino Physics (2012)
[4] Bungau, A.; Adelmann, A., An electron antineutrino disappearance search using high-rate 8Li production and decay, Phys. Rev. Lett., 109 (2012), 141802
[5] Toggweiler, M., An adaptive time integration method for more efficient simulation of particle accelerators (2011), Computer Science Department: Computer Science Department ETH Zurich, Master thesis
[6] Gear, C. W., Numerical Initial Value Problems in Ordinary Differential Equations (1971), Prentice-Hall: Prentice-Hall Englewood Cliffs, NJ · Zbl 1145.65316
[7] Hairer, E.; Lubich, C.; Wanner, G., Geometric Numerical Integration: Structure-Preserving Algorithms for Ordinary Differential Equations (2006), Springer: Springer Berlin · Zbl 1094.65125
[8] Landau, L. D.; Lifshitz, E. M., Electrodynamics of Continuous Media (1984), Pergamon: Pergamon Oxford
[9] Qiang, J.; Ryne, R. D., Parallel 3D Poisson solver for a charged beam in a conducting pipe, Comput. Phys. Commun., 138, 18-28 (2001) · Zbl 0982.78024
[10] Qiang, J.; Gluckstern, R. L., Three-dimensional Poisson solver for a charged beam with large aspect ratio in a conducting pipe, Comput. Phys. Commun., 160, 120-128 (2004) · Zbl 1196.74280
[11] Birdsall, C. K.; Langdon, A. B., Plasma Physics via Computer Simulation (1985), McGraw-Hill
[12] Sundman, K., Memoire sur le probleme des trois corps, Acta Math., 36, 105-179 (1912) · JFM 43.0826.01
[13] Huang, W.; Leimkuhler, B., The adaptive Verlet method, SIAM J. Sci. Comput., 18, 239-256 (1997) · Zbl 0877.65048
[14] Bond, S.; Leimkuhler, B., Time-transformations for reversible variable stepsize integration, Numer. Algorithms, 19, 55-71 (1998) · Zbl 0923.70002
[15] Merson, R. H., Numerical Integration of the Differential Equations of Celestial Mechanics (1975), European Space Research Organisation, Royal Aircraft Establishment: European Space Research Organisation, Royal Aircraft Establishment Great Britain
[16] Berry, M. M., A variable-step double-integration multi-step integrator (2004), Virginia Polytechnic Institute and State University, Ph.D. thesis
[17] Yang, J.; Adelmann, A.; Humbel, M.; Seidel, M.; Zhang, T., Beam dynamics in high intensity cyclotrons including neighboring bunch effects: model, implementation and application, Phys. Rev. Spec. Top., Accel. Beams, 13 (2010), 064201
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.