A quasi-static particle-in-cell algorithm based on an azimuthal Fourier decomposition for highly efficient simulations of plasma-based acceleration: QPAD
Journal Article
·
· Computer Physics Communications
- Univ. of California, Los Angeles, CA (United States); SLAC National Accelerator Lab., Menlo Park, CA (United States)
- Beijing Normal Univ. (China)
- Univ. of California, Los Angeles, CA (United States)
- SLAC National Accelerator Lab., Menlo Park, CA (United States)
The three-dimensional (3D) quasi-static particle-in-cell (PIC) algorithm is a very efficient method for modeling short-pulse laser or relativistic charged particle beam–plasma interactions. In this algorithm, the plasma response, i.e., plasma wave wake, to a non-evolving laser or particle beam is calculated using a set of Maxwell’s equations based on the quasi-static approximate equations that exclude radiation. The plasma fields are then used to advance the laser or beam forward using a large time step. The algorithm is many orders of magnitude faster than a 3D fully explicit relativistic electromagnetic PIC algorithm. It has been shown to be capable to accurately model the evolution of lasers and particle beams in a variety of scenarios. Additionally, at the same time, an algorithm in which the fields, currents and Maxwell equations are decomposed into azimuthal harmonics has been shown to reduce the algorithmic complexity of a 3D explicit PIC algorithm to that of a 2D algorithm when the expansion is truncated while maintaining accuracy for problems with near azimuthal symmetry. This hybrid algorithm uses a PIC description in r–z and a gridless description in . We describe a novel method that combines the quasi-static and hybrid PIC methods. This algorithm expands the fields, charge and current density into azimuthal harmonics. A set of the quasi-static field equations is derived for each harmonic. The complex amplitudes of the fields are then solved using the finite difference method. The beam and plasma particles are advanced in Cartesian coordinates using the total fields. Details on how this algorithm was implemented using a similar workflow to an existing quasi-static code, QuickPIC, are presented. The new code is called QPAD for QuickPIC with Azimuthal Decomposition. Benchmarks and comparisons between a fully 3D explicit PIC code (OSIRIS), a full 3D quasi-static code (QuickPIC), and the new quasi-static PIC code with azimuthal decomposition (QPAD) are also presented.
- Research Organization:
- SLAC National Accelerator Laboratory (SLAC), Menlo Park, CA (United States)
- Sponsoring Organization:
- National Science Foundation (NSF); USDOE; USDOE Office of Science (SC), Basic Energy Sciences (BES)
- Grant/Contract Number:
- AC02-76SF00515; SC0010064
- OSTI ID:
- 1774944
- Alternate ID(s):
- OSTI ID: 1775682
- Journal Information:
- Computer Physics Communications, Journal Name: Computer Physics Communications Vol. 261; ISSN 0010-4655
- Publisher:
- ElsevierCopyright Statement
- Country of Publication:
- United States
- Language:
- English
Similar Records
Integrating a ponderomotive guiding center algorithm into a quasi-static particle-in-cell code based on azimuthal mode decomposition
Implementation of a Mesh refinement algorithm into the quasi-static PIC code QuickPIC
An improved iteration loop for the three dimensional quasi-static particle-in-cell algorithm: QuickPIC
Journal Article
·
Tue Sep 06 20:00:00 EDT 2022
· Journal of Computational Physics
·
OSTI ID:2419983
Implementation of a Mesh refinement algorithm into the quasi-static PIC code QuickPIC
Journal Article
·
Tue Sep 30 20:00:00 EDT 2025
· Journal of Computational Physics
·
OSTI ID:3012400
An improved iteration loop for the three dimensional quasi-static particle-in-cell algorithm: QuickPIC
Journal Article
·
Tue Oct 01 00:00:00 EDT 2013
· Journal of Computational Physics
·
OSTI ID:22230797