An Energy- and Charge-conserving, Implicit, Electrostatic Particle-in-Cell Algorithm
arXiv:1101.3701 · doi:10.1016/j.jcp.2011.05.031
Abstract
This paper discusses a novel fully implicit formulation for a 1D electrostatic particle-in-cell (PIC) plasma simulation approach. Unlike earlier implicit electrostatic PIC approaches (which are based on a linearized Vlasov-Poisson formulation), ours is based on a nonlinearly converged Vlasov-Ampère (VA) model. By iterating particles and fields to a tight nonlinear convergence tolerance, the approach features superior stability and accuracy properties, avoiding most of the accuracy pitfalls in earlier implicit PIC implementations. In particular, the formulation is stable against temporal (CFL) and spatial (aliasing) instabilities. It is charge- and energy-conserving to numerical roundoff for arbitrary implicit time steps. While momentum is not exactly conserved, errors are kept small by an adaptive particle sub-stepping orbit integrator, which is instrumental to prevent particle tunneling. The VA model is orbit-averaged along particle orbits to enforce an energy conservation theorem with particle sub-stepping. As a result, very large time steps, constrained only by the dynamical time scale of interest, are possible without accuracy loss. Algorithmically, the approach features a Jacobian-free Newton-Krylov solver. A main development in this study is the nonlinear elimination of the new-time particle variables (positions and velocities). Such nonlinear elimination, which we term particle enslavement, results in a nonlinear formulation with memory requirements comparable to those of a fluid computation, and affords us substantial freedom in regards to the particle orbit integrator. Numerical examples are presented that demonstrate the advertised properties of the scheme. In particular, long-time ion acoustic wave simulations show that numerical accuracy does not degrade even with very large implicit time steps, and that significant CPU gains are possible.
29 pages, 8 figures, submitted to Journal of Computational Physics
References in corpus (1)
Cited by in corpus (73)
- Magnetic reconnection in the era of exascale computing and multiscale experiments
- Exactly Energy Conserving Semi-Implicit Particle in Cell Formulation
- Variational formulation of particle algorithms for kinetic plasma simulations
- A multi-dimensional, energy- and charge-conserving, nonlinearly implicit, electromagnetic Vlasov-Darwin particle-in-cell algorithm
- A Hamiltonian splitting for the Vlasov-Maxwell system
- Energy-conserving discontinuous Galerkin methods for the Vlasov-Ampère system
- An energy- and charge-conserving, nonlinearly implicit, electromagnetic 1D-3V Vlasov-Darwin particle-in-cell algorithm
- Structure-Preserving Geometric Particle-in-Cell Methods for Vlasov-Maxwell Systems
- CPIC: a Curvilinear Particle-In-Cell code for plasma-material interaction studies
- An efficient mixed-precision, hybrid CPU-GPU implementation of a fully implicit particle-in-cell algorithm
- Energy-conserving time propagation for a geometric particle-in-cell Vlasov--Maxwell solver
- Temporal resolution criterion for correctly simulating relativistic electron motion in a high-intensity laser field
- Energy-conserving Discontinuous Galerkin Methods for the Vlasov-Maxwell System
- Fluid preconditioning for Newton-Krylov-based, fully implicit, electrostatic particle-in-cell simulations
- The Fluid-Kinetic Particle-in-Cell Solver for Plasma Simulations
- SHARP: A Spatially Higher-order, Relativistic Particle-in-Cell Code
- Gauss's Law Satisfying Energy-Conserving Semi-Implicit Particle-in-Cell Method
- An energy-conserving and asymptotic-preserving charged-particle orbit implicit time integrator for arbitrary electromagnetic fields
- Metriplectic Integrators for the Landau Collision Operator
- On Energy and Momentum Conservation in Particle-in-Cell Simulation
- A semi-implicit, energy- and charge-conserving particle-in-cell algorithm for the relativistic Vlasov-Maxwell equations
- Quadratic conservative scheme for relativistic Vlasov--Maxwell system
- A fully implicit, conservative, non-linear, electromagnetic hybrid particle-ion/fluid-electron algorithm
- Enabling particle applications for exascale computing platforms
- Finite grid instability and spectral fidelity of the electrostatic Particle-In-Cell algorithm
- Arbitrarily high order Convected Scheme solution of the Vlasov-Poisson system
- Multiple scale kinetic simulations with the energy conserving semi implicit particle in cell (PIC) method
- Performance analysis and implementation details of the Energy Conserving Semi-Implicit Method code (ECsim)
- Explicit Structure-Preserving Geometric Particle-in-Cell Algorithm in Curvilinear Orthogonal Coordinate Systems and Its Applications to Whole-Device 6D Kinetic Simulations of Tokamak Physics
- Electron dynamics in planar radio frequency magnetron plasmas: I. The mechanism of Hall heating and the μ-mode
- Electron dynamics in planar radio frequency magnetron plasmas: II. Heating and energization mechanisms studied via a 2d3v particle-in-cell/Monte Carlo code
- ECsim-CYL: Energy Conserving Semi-Implicit particle in cell simulation in axially symmetric cylindrical coordinates
- The geometric theory of charge conservation in particle-in-cell simulations
- Modeling of very high frequency large-electrode capacitively coupled plasmas with a fully electromagnetic particle-in-cell code
- PolyPIC: the Polymorphic-Particle-in-Cell Method for Fluid-Kinetic Coupling
- Advances in the Implementation of the Exactly Energy Conserving Semi-Implicit (ECsim) Particle in Cell Method
- Electron dynamics in planar radio frequency magnetron plasmas: III. Comparison of experimental investigations of power absorption dynamics to simulation results
- Efficient Fourier Basis Particle Simulation
- On the Numerical Dispersion of Electromagnetic Particle-In-Cell Code : Finite Grid Instability
- 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
- Highly efficient energy-conserving moment method for the multi-dimensional Vlasov-Maxwell system
- An efficient, conservative, time-implicit solver for the fully kinetic arbitrary-species 1D-2V Vlasov-Ampère system
- Verification of a Fully Implicit Particle-in-Cell Method for the Formalism of Electromagnetic Gyrokinetics in the XGC Code
- A Maxwell-Ampère Nernst-Planck Framework for Modeling Charge Dynamics
- Time Integrator Agnostic Charge Conserving Finite Element PIC
- Hamiltonian Particle-in-Cell methods for Vlasov-Poisson equations
- An asymptotic preserving scheme for the relativistic Vlasov--Maxwell equations in the classical limit
- Application of the phase space action principle to finite-size particle plasma simulations in the drift-kinetic approximation
- Arbitrary-order time-accurate semi-Lagrangian spectral approximations of the Vlasov-Poisson system
- A multilevel Monte Carlo method for a class of McKean-Vlasov processes
- An Asymptotic-Preserving and Energy-Conserving Particle-In-Cell Method for Vlasov-Maxwell Equations
- Numerical Study of the Two-Species Vlasov-Ampère System: Energy-Conserving Schemes and the Current-Driven Ion-Acoustic Instability
- Subcycling of particle orbits in variational, geometric electromagnetic particle-in-cell methods
- A Semi-Lagrangian Spectral Method for the Vlasov-Poisson System based on Fourier, Legendre and Hermite Polynomials
- Noise and error analysis and optimization in particle-based kinetic plasma simulations
- Reducing Noise for PIC Simulations Using Kernel Density Estimation Algorithm
- Implicit temporal discretization and exact energy conservation for particle methods applied to the Poisson-Boltzmann equation
- An implicit, conservative electrostatic particle-in-cell algorithm for paraxial magnetic nozzles
- Conservative Projection Between Finite Element and Particle Bases
- Explicit energy-conserving modification of relativistic PIC method
- Variational Formulation of Macro-Particle Models for Electromagnetic Plasma Simulations
- Asymptotic-preserving Particle-In-Cell methods for the Vlasov-Maxwell system near quasi-neutrality
- Grid instability growth rates for explicit, electrostatic momentum- and energy-conserving particle-in-cell algorithms
- Comparison of Integrators for Electromagnetic Particle in Cell Methods: Algorithms and Applications
- Machine learning accelerated particle-in-cell plasma simulations
- Finite spatial-grid effects in energy-conserving particle-in-cell algorithms
- Kinetic Simulation of Collisional Magnetized Plasmas with Semi-Implicit Time Integration
- Simplex-in-Cell Technique for Collisionless Plasma Simulations
- Charge-conserving, variational particle-in-cell method for the drift-kinetic Vlasov-Maxwell system
- Bridging kinetic plasma descriptions and single fluid models
- Optimal Landau-type closure parameters for two-fluid simulations of plasma turbulence at kinetic scales
- Projective multiscale time-integration for electrostatic particle-in-cell methods