Abstract:High fidelity modeling of plasma based acceleration (PBA) requires the use of three dimensional, fully nonlinear, and kinetic descriptions based on the particle-in-cell (PIC) method. In PBA an intense particle beam or laser (driver) propagates through a tenuous plasma whereby it excites a plasma wave wake. Three-dimensional PIC algorithms based on the quasi-static approximation (QSA) have been successfully applied to efficiently model the interaction between relativistic charged particle beams and plasma. In a QSA PIC algorithm, the plasma response to a charged particle beam or laser driver is calculated based on forces from the driver and self-consistent forces from the QSA form of Maxwell's equations. These fields are then used to advance the charged particle beam or laser forward by a large time step. Since the time step is not limited by the regular Courant-Friedrichs-Lewy (CFL) condition that constrains a standard 3D fully electromagnetic PIC code, a 3D QSA PIC code can achieve orders of magnitude speedup in performance. Recently, a new hybrid QSA PIC algorithm that combines another speedup technique known as an azimuthal Fourier decomposition has been proposed and implemented. This hybrid algorithm decomposes the electromagnetic fields, charge and current density into azimuthal harmonics and only the Fourier coefficients need to be updated, which can reduce the algorithmic complexity of a 3D code to that of a 2D code. Modeling the laser-plasma interaction in a full 3D electromagnetic PIC algorithm is very computationally expensive due the enormous disparity of physical scales to be resolved. In the QSA the laser is modeled using the ponderomotive guiding center (PGC) approach. We describe how to implement a PGC algorithm compatible for the QSA PIC algorithms based on the azimuthal mode expansion. This algorithm permits time steps orders of magnitude larger than the cell size and it can be asynchronously parallelized. Details on how this is implemented into the QSA PIC code that utilizes an azimuthal mode expansion, QPAD, are also described. Benchmarks and comparisons between a fully 3D explicit PIC code (OSIRIS), as well as a few examples related to laser wakefield acceleration, are presented.

Concept of the generalized reduced-order particle-in-cell scheme and verification in an axial-azimuthal Hall thruster configuration

Generalized reduced-order particle-in-cell scheme for Hall thruster modeling: concept and in-depth verification in the axial-azimuthal configuration

Dimensionality reduction of Poisson's equation with application to particle-in-cell simulations of Hall thrusters

Enhancing one-dimensional particle-in-cell simulations to self-consistently resolve instability-induced electron transport in Hall thrusters

Reduced-order particle-in-cell simulations of a high-power magnetically shielded Hall thruster

Latest progress on the reduced-order particle-in-cell scheme: II. Quasi-3D implementation and verification

Latest progress on the reduced-order particle-in-cell scheme: I. refining the underlying formulation

Effects of the neutral dynamics model on the particle-in-cell simulations of a Hall thruster plasma discharge

A quasi-static particle-in-cell algorithm based on an azimuthal Fourier decomposition for highly efficient simulations of plasma-based acceleration: QPAD

Non-Maxwellian electron effects on the macroscopic response of a Hall thruster discharge from an axial-radial kinetic model

Physics, simulation and diagnostics of Hall effect thrusters

Integrating a ponderomotive guiding center algorithm into a quasi-static particle-in-cell code based on azimuthal mode decomposition

A two-dimensional hybrid model of the Hall thruster discharge

Coupling plasma physics and chemistry in the PIC model of electric propulsion: Application to an air-breathing, low-power Hall thruster

Ion angular distribution simulation of the Highly Efficient Multistage Plasma Thruster

Kinetic, 3-D, PIC-DSMC Simulations of Ion Thruster Plumes and the Backflow Region

High-order Implicit Particle-in-cell Method for Plasma Simulations at Solid Densities.

Geometric Particle-In-Cell discretizations of a plasma hybrid model with kinetic ions and mass-less fluid electrons

Study of scaling law for particle-in-cell/Monte Carlo simulation of low-temperature magnetized plasma for electric propulsion

Latest progress in Hall thrusters plasma modelling

Two-dimensional Modeling of the Hall Thruster Discharge with Non-uniform Propellant Supply in Azimuth