Reuse & Permissions

It is not necessary to obtain permission to reuse this article or its components as it is available under the terms of the Creative Commons Attribution 3.0 License. This license permits unrestricted use, distribution, and reproduction in any medium, provided attribution to the author(s) and the published article's title, journal citation, and DOI are maintained. Please note that some figures may have been included with permission from other third parties. It is your responsibility to obtain the proper permission from the rights holder directly for these figures.

Export citation

Export citation

Choose format for download:

Download Citation
  • Open Access

Particle-in-cell beam dynamics simulations with a wavelet-based Poisson solver

Balša Terzić1, Ilya V. Pogorelov2, and Courtlandt L. Bohn1

  • 1Beam Physics and Astrophysics Group, Northern Illinois University, DeKalb, Illinois 60115, USA
  • 2Accelerator and Fusion Research Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA

Phys. Rev. ST Accel. Beams 10, 034201 – Published 2 March, 2007

DOI: https://doi.org/10.1103/PhysRevSTAB.10.034201

Abstract

We report on a successful implementation of a three-dimensional wavelet-based solver for the Poisson equation with Dirichlet boundary conditions, optimized for use in particle-in-cell (PIC) simulations. The solver is based on the operator formulation of the conjugate gradient algorithm, for which effectively diagonal preconditioners are available in wavelet bases. Because of the recursive nature of PIC simulations, a good initial approximation to the iterative solution is always readily available, which we demonstrate to be a key advantage in terms of overall computational speed. While the Laplacian remains sparse in a wavelet representation, the wavelet-decomposed potential and density can be rendered sparse through a procedure that amounts to simultaneous compression and denoising of the data. We explain how this procedure can be carried out in a controlled and near-optimal way, and show the effect it has on the overall solver performance. After testing the solver in a stand-alone mode, we integrated it into the IMPACT-T beam dynamics particle-in-cell code and extensively benchmarked it against the IMPACT-T with the native FFT-based Poisson solver. We present and discuss these benchmarking results, as well as the results of modeling the Fermi/NICADD photoinjector using IMPACT-T with the wavelet-based solver.

View figure in article

Article Text

References (48)

  1. R. Hockney and J. Eastwood, Computer Simulation Using Particles (Institute of Physics Publishing, London, 1988).
  2. J. Dawson, Rev. Mod. Phys. 55, 403 (1983).
  3. J. Qiang, R. Ryne, S. Habib, and V. Decyk, J. Comput. Phys. 163, 434 (2000).
  4. F. Combes, F. Debbasch, D. Friedli, and D. Pfenniger, Astron. Astrophys. 233, 82 (1990).
  5. H. Couchman, Astrophys. J. 368, L23 (1991).
  6. A. Melott, R. Splinter, S. Shandarin, and Y. Suto, Astrophys. J. 479, L79 (1997).
  7. R. Splinter, A. Melott, S. Shandarin, and Y. Suto, Astrophys. J. 497, 38 (1998).
  8. A. Romeo, C. Horellou, and J. Bergh, Mon. Not. R. Astron. Soc. 342, 337 (2003).
  9. A. Romeo, C. Horellou, and J. Bergh, Mon. Not. R. Astron. Soc. 354, 1208 (2004).
  10. K. Heitmann, P. Ricker, M. Warren, and S. Habib, Astrophys. J. Suppl. Ser. 160, 28 (2005).
  11. C. Birdsall and A. Langdon, Plasma Physics via Computer Simulations (Institute of Physics Publishing, London, 1991).
  12. A. Friedman, D. Grote, and I. Haber, Phys. Fluids B 4, 2203 (1992).
  13. F. H. Harlow, Meth. Comput. Phys. 3, 319 (1964).
  14. D. Hamilton, F. Lindholm, and A. Marshak, Principles and Applications of Semiconductor Device Modeling (Rinehart and Winston, New York, 1971).
  15. G. Beylkin, in Wavelets: Mathematics and Applications, edited by J. Benedetto and M. Frazier (CRC Press LLC, Boca Raton, Florida, 1993).
  16. A. Averbuch, G. Beylkin, R. Coifman, P. Fischer, and M. Israeli, in Signal and Image Representation in Combined Spaces, edited by Y. Zeevi and R. Coifman (Academic Press, San Diego, 1995).
  17. Z. Huang, M. Borland, P. Emma, J. Wu, C. Limborg, G. Stupakov, and J. Welch, Phys. Rev. ST Accel. Beams 7, 074401 (2004).
  18. I. Pogorelov, J. Qiang, R. Ryne, M. Venturini, A. Zholents, and R. Warnock, Proceedings of the 9th International Computational Accelerator Physics Conference, 2006, p. 182 (WEPPP01).
  19. M. Berger and P. Colella, J. Comput. Phys. 82, 64 (1989).
  20. J. Bella, M. Berger, J. Saltzman, and M. Welcome, SIAM J. Sci. Comput. 15, 127 (1994).
  21. R. Coifman and M. Wickerhauser, IEEE Trans. Inf. Theory 38, 713 (1992).
  22. R. Coifman, Y. Meyer, and M. Wickerhauser, in Wavelets and Their Applications, edited by M. B. Ruskai (Jones and Bartlett, Boston, 1992).
  23. D. Donoho and I. Johnstone, Biometrika 81, 425 (1994).
  24. D. Donoho and I. Johnstone, Stanford University internal report, 1994, http://www-stat.stanford.edu/~donoho/Reports/1994/idealbasis.pdf.
  25. N. Saito, in Wavelets in Geophysics, edited by E. Foufoula-Georgiou and P. Kumar (Academic Press, San Diego, 1994).
  26. C. Bohn, AIP Conf. Proc. 647, 81 (2002).
  27. G. Bassi, T. Agoh, M. Dolhus, L. Giannessi, R. Hajima, A. Kabel, T. Limberg, and M. Quattromini, Nucl. Instrum. Methods Phys. Res., Sect. A 557, 189 (2006).
  28. A. Averbuch, G. Beylkin, R. Coifman, P. Fischer, and M. Israeli, Institut de Mathématiques de Bordeaux Internal Report, 2003.
  29. J. Qiang and R. Ryne, Comput. Phys. Commun. 138, 18 (2001).
  30. J. Qiang, S. Lidia, R. Ryne, and C. Limborg-Deprey, Phys. Rev. ST Accel. Beams 9, 044204 (2006).
  31. I. Daubechies, Ten Lectures on Wavelets (SIAM, Philadelphia, 1992).
  32. M. Wickerhauser, Adaptive Wavelet Analysis From Theory to Software (A K Peters, Wellesley, Massachusetts, 1994).
  33. M. Misiti, Y. Misiti, G. Oppenheim, and J. Poggi, Wavelet Toolbox for Use with MATLAB (The MathWorks, Inc., Natick, Massachusetts, 1997).
  34. S. Goedecker, Wavelets and Their Applications (Presses Polytechniques et Universitaires Romandes, Lausanne, 1998).
  35. S. Mallat, A Wavelet Tour of Signal Processing (Elsevier, San Diego, 1998).
  36. S. Goedecker and O. Ivanov, Comput. Phys. 12, 548 (1998).
  37. G. Golub and C. Van Loan, Matrix Computations (Johns Hopkins, Baltimore, 1996).
  38. M. R. Hestenes and E. Stiefel, J. Res. Nat. Bur. Stand. 49, 409 (1952).
  39. W. Press, S. Teukolsky, W. Vetterling, and B. Flannery, Numerical Recipes in Fortran: The Art of Scientific Computing, 2nd edition (University of Cambridge, Cambridge, England, 1992).
  40. R. Coifman and D. Donoho, in Wavelet and Statistics, edited by A. Antoniadis and G. Oppenheim (Springer, New York, 1995).
  41. F. Anscombe, Biometrika 35, 246 (1948).
  42. A. Stuart and J. Ord, Advanced Theory of Statistics (Oxford University Press, New York, 1991), Vol. 2.
  43. D. Donoho, Stanford University Internal Report, 1993, http://www-stat.stanford.edu/~donoho/Reports/1993/toulouse.ps.Z.
  44. F. Murtagh, J. Starck, and A. Bijaoui, Astron. Astrophys. Suppl. Ser. 112, 179 (1995).
  45. J. Starck, D. Donoho, and E. Candés, Astron. Astrophys. 398, 785 (2003).
  46. H. Plummer, Mon. Not. R. Astron. Soc. 71, 460 (1911).
  47. See http://www.nicadd.niu.edu/fnpl/
  48. R. Tikhoplav, A. Melissinos, J. Li, and P. Piot, Particle Accelerator Conference, 2005.

Outline

Information

Sign In to Your Journals Account

Filter

Filter

Article Lookup

Enter a citation