Computational Methods in the Warp Code Framework for Kinetic Simulations of Particle Beams and Plasmas

• Published in: IEEE transactions on plasma science Vol. 42; no. 5; pp. 1321 - 1334
• Main Authors: Friedman, Alex, Cohen, Ronald H., Grote, David P., Lund, Steven M., Sharp, William M., Vay, Jean-Luc, Haber, Irving, Kishek, Rami A.
• Format: Journal Article
• Language: English
• Published: New York IEEE 01.05.2014; The Institute of Electrical and Electronics Engineers, Inc. (IEEE)
• Subjects: Acceleration; Algorithms; computer; Electromagnetics; Electrostatics; Geometry; Laboratories; laser; Laser beams; Magnetic fields; Maxwell; Ned Birdsall; numerical simulation; particle beam; Particle beams; particle-in-cell; plasma; Plasma physics; Plasmas; Simulation; Structural beams; computer; particle-in-cell; Algorithms; laser; Maxwell; particle beam; plasma; Ned Birdsall; numerical simulation
• ISSN: 0093-3813; 1939-9375
• DOI: 10.1109/TPS.2014.2308546

The Warp code (and its framework of associated tools) was initially developed for particle-in-cell simulations of space-charge-dominated ion beams in accelerators, for heavy-ion-driven inertial fusion energy, and related experiments. It has found a broad range of applications, including nonneutral plasmas in traps, stray electron clouds in accelerators, laser-based acceleration, and the focusing of ion beams produced when short-pulse lasers irradiate foil targets. We summarize novel methods used in Warp, including: time-stepping conducive to diagnosis and particle injection; an interactive Python-Fortran-C structure that enables scripted and interactive user steering of runs; a variety of geometries (3-D x, y, z; 2-D r, z; 2-D x, y); electrostatic and electromagnetic field solvers; a cut-cell representation for internal boundaries; the use of warped coordinates for bent beam lines; adaptive mesh refinement, including a capability for time-dependent space-charge-limited flow from curved surfaces; models for accelerator lattice elements (magnetic or electrostatic quadrupole lenses, accelerating gaps, etc.) at user-selectable levels of detail; models for particle interactions with gas and walls; moment/envelope models that support sophisticated particle loading; a drift-Lorentz mover for rapid tracking through regions of strong and weak magnetic field; a Lorentz-boosted frame formulation with a Lorentz-invariant modification of the Boris mover; an electromagnetic solver with tunable dispersion and stride-based digital filtering; and a pseudospectral electromagnetic solver. Warp has proven useful for a wide range of applications, described very briefly herein. It is available as an open-source code under a BSD license. This paper describes material presented during the Prof. Charles K. (Ned) Birdsall Memorial Session of the 2013 IEEE Pulsed Power and Plasma Science Conference. In addition to our overview of the computational methods used in Warp, we summarize a few aspects of Ned's contributions to plasma simulation and to the careers of those he mentored.