An energy- and charge-conserving, nonlinearly implicit, electromagnetic 1D-3V Vlasov-Darwin particle-in-cell algorithm

Abstract

A recent proof-of-principle study proposes a nonlinear electrostatic implicit particle-in-cell (PIC) algorithm in one dimension (Chen, Chacon, Barnes, J. Comput. Phys. 230 (2011) 7018). The algorithm employs a kinetically enslaved Jacobian-free Newton-Krylov (JFNK) method, and conserves energy and charge to numerical round-off. In this study, we generalize the method to electromagnetic simulations in 1D using the Darwin approximation of Maxwell's equations, which avoids radiative aliasing noise issues by ordering out the light wave. An implicit, orbit-averaged time-space-centered finite difference scheme is applied to both the 1D Darwin field equations (in potential form) and the 1D-3V particle orbit equations to produce a discrete system that remains exactly charge- and energy-conserving. Furthermore, enabled by the implicit Darwin equations, exact conservation of the canonical momentum per particle in any ignorable direction is enforced via a suitable scattering rule for the magnetic field. Several 1D numerical experiments demonstrate the accuracy and the conservation properties of the algorithm.

24 pages, 4 figures

References in corpus (1)

• Fluid preconditioning for Newton-Krylov-based, fully implicit, electrostatic particle-in-cell simulations

Cited by in corpus (20)

• A multi-dimensional, energy- and charge-conserving, nonlinearly implicit, electromagnetic Vlasov-Darwin particle-in-cell algorithm
• Temporal resolution criterion for correctly simulating relativistic electron motion in a high-intensity laser field
• Fluid preconditioning for Newton-Krylov-based, fully implicit, electrostatic particle-in-cell simulations
• An energy-conserving and asymptotic-preserving charged-particle orbit implicit time integrator for arbitrary electromagnetic fields
• Metriplectic Integrators for the Landau Collision Operator
• A semi-implicit, energy- and charge-conserving particle-in-cell algorithm for the relativistic Vlasov-Maxwell equations
• Neutral Vlasov kinetic theory of magnetized plasmas
• Quadratic conservative scheme for relativistic Vlasov--Maxwell system
• ViDA: a Vlasov-DArwin solver for plasma physics at electron scales
• Verification of long wavelength electromagnetic modes with a gyrokinetic-fluid hybrid model in the XGC code
• Modeling of very high frequency large-electrode capacitively coupled plasmas with a fully electromagnetic particle-in-cell code
• An implicit, conservative and asymptotic-preserving electrostatic particle-in-cell algorithm for arbitrarily magnetized plasmas in uniform magnetic fields
• Rubrics for Charge Conserving Current Mapping in Finite Element Particle in Cell Methods
• Verification of a Fully Implicit Particle-in-Cell Method for the $v_{\parallel }$ Formalism of Electromagnetic Gyrokinetics in the XGC Code
• Comparison of linear modes in kinetic plasma models
• Subcycling of particle orbits in variational, geometric electromagnetic particle-in-cell methods
• An implicit, conservative electrostatic particle-in-cell algorithm for paraxial magnetic nozzles
• A Particle-in-cell Method for Plasmas with a Generalized Momentum Formulation, Part II: Enforcing the Lorenz Gauge Condition
• Finite spatial-grid effects in energy-conserving particle-in-cell algorithms
• Charge-conserving, variational particle-in-cell method for the drift-kinetic Vlasov-Maxwell system