Variational formulation of particle algorithms for kinetic plasma simulations

Abstract

Common time-explicit numerical methods for kinetic simulations of plasmas in the low-collisions limit fall into two classes of algorithms: momentum conserving and energy conserving. Each has certain drawbacks. The PIC algorithm does not conserve total energy, which may lead to spurious numerical heating (grid heating). Its overall accuracy is at most second due to the nature of the force interpolation between grid and particle position. Energy-conserving algorithms do not exhibit grid heating, but because their formulation uses potentials, computationally undesirable matrix inversions may be necessary. In addition, compared to PIC algorithms for the same accuracy, these algorithms have higher numerical noise due to the restricted choice of particle shapes. Here we formulate time-explicit, finite-size particle algorithms using particular reductions of the particle distribution function. These reductions are used in two variational principles, a Lagrangian-based and a Hamiltonian-based in conjunction with a non-canonical Poisson bracket. The Lagrangian formulations here generalize previous such formulations. The Hamiltonian formulation is presented here for the first time. Many drawbacks of the two classes of particle methods are mitigated. For example, restrictions on particle shapes are relaxed in energy conserving algorithms, which allows to decrease the numerical noise in these methods. The Hamiltonian formulation of particle algorithms is done in terms of fields instead of potentials, thus avoiding solving Poisson's equation. An algorithm that conserves both energy and momentum is presented. Other features of the algorithms include a natural way to perform coordinate transformations, the use of various time integrating methods, and the ability to increase the overall accuracy beyond second order, including all generalizations.

Cited by in corpus (61)

• GEMPIC: Geometric ElectroMagnetic Particle-In-Cell Methods
• Structure and structure-preserving algorithms for plasma physics
• Explicit high-order non-canonical symplectic particle-in-cell algorithms for Vlasov-Maxwell systems
• Canonical symplectic particle-in-cell method for long-term large-scale simulations of the Vlasov-Maxwell system
• Structure-Preserving Geometric Particle-in-Cell Methods for Vlasov-Maxwell Systems
• Explicit symplectic algorithms based on generating functions for charged particle dynamics
• Hamiltonian integration methods for Vlasov-Maxwell equations
• High Order Maximum Principle Preserving Semi-Lagrangian Finite Difference WENO schemes for the Vlasov Equation
• Hamiltonian particle-in-cell methods for Vlasov-Maxwell equations
• Metriplectic Integrators for the Landau Collision Operator
• Explicit high-order noncanonical symplectic algorithms for ideal two-fluid systems
• A fully implicit, conservative, non-linear, electromagnetic hybrid particle-ion/fluid-electron algorithm
• High order volume-preserving algorithms for relativistic charged particles in general electromagnetic fields
• Finite grid instability and spectral fidelity of the electrostatic Particle-In-Cell algorithm
• Variational principles for the guiding-center Vlasov-Maxwell equations
• 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
• Action Principles for Relativistic Extended Magnetohydrodynamics: A Unified Theory of Magnetofluid Models
• Degenerate Variational Integrators for Magnetic Field Line Flow and Guiding Center Trajectories
• Detection and Prediction of Equilibrium States in Kinetic Plasma Simulations via Mode Tracking using Reduced-Order Dynamic Mode Decomposition
• A Symplectic Multi-Particle Tracking Model for Self-Consistent Space-Charge Simulation
• A conservative semi-Lagrangian HWENO method for the Vlasov equation
• Lagrangian and Dirac constraints for the ideal incompressible fluid and magnetohydrodynamics
• Structure-preserving marker-particle discretizations of Coulomb collisions for particle-in-cell codes
• On the structure and statistical theory of turbulence of extended magnetohydrodynamics
• Efficient Fourier Basis Particle Simulation
• Field theory and structure-preserving geometric particle-in-cell algorithm for drift wave instability and turbulence
• Explicit High-Order Gauge-Independent Symplectic Algorithms for Relativistic Charged Particle Dynamics
• An Action Principle for Relativistic MHD
• A symplectic particle-in-cell model for space-charge beam dynamics simulation
• The Plasma Simulation Code: A modern particle-in-cell code with load-balancing and GPU support
• Classical molecular dynamics simulations of hydrogen plasmas and development of an analytical statistical model for computational validity assessment
• Coherent buildup of high harmonic radiation: the classical perspective
• Explicit symplectic algorithms based on generating functions for relativistic charged particle dynamics in time-dependent electromagnetic field
• Highly efficient energy-conserving moment method for the multi-dimensional Vlasov-Maxwell system
• Hamiltonian Particle-in-Cell methods for Vlasov-Poisson equations
• Local Energy Conservation Law for Spatially-Discretized Hamiltonian Vlasov-Maxwell System
• Quasi-Helmholtz Decomposition, Gauss' Laws and Charge Conservation for Finite Element Particle-in-Cell
• Application of the phase space action principle to finite-size particle plasma simulations in the drift-kinetic approximation
• Stochastic variational principles for the collisional Vlasov-Maxwell and Vlasov-Poisson equations
• An Asymptotic-Preserving and Energy-Conserving Particle-In-Cell Method for Vlasov-Maxwell Equations
• Subcycling of particle orbits in variational, geometric electromagnetic particle-in-cell methods
• Noise and error analysis and optimization in particle-based kinetic plasma simulations
• On the velocity space discretization for the Vlasov-Poisson system: comparison between Hermite spectral and Particle-in-Cell methods. Part 2: fully-implicit scheme
• Multispecies structure-preserving particle discretization of the Landau collision operator
• Metriplectic foundations of gyrokinetic Vlasov-Maxwell-Landau theory
• Accelerating Particle-in-Cell Kinetic Plasma Simulations via Reduced-Order Modeling of Space-Charge Dynamics using Dynamic Mode Decomposition
• Testing the conservative character of particle simulations: II. Spurious heating of guiding centers and full orbits subject to fluctuations expressed in terms of $\mathbf{E}$ and $\mathbf{B}$
• Hamiltonian structure of the gauge-free gyrokinetic Vlasov-Maxwell equations
• A gauge-compatible Hamiltonian splitting algorithm for particle-in-cell simulations using finite element exterior calculus
• Classical versus quantum views of intense laser pulse propagation in gases
• Time simulation of the nonlinear wave-particle interaction in meters long traveling-wave tubes
• Variational Formulation of Macro-Particle Models for Electromagnetic Plasma Simulations
• Variational formulation of classical and quantum models for intense laser pulse propagation
• Grid instability growth rates for explicit, electrostatic momentum- and energy-conserving particle-in-cell algorithms
• Simplex-in-Cell Technique for Collisionless Plasma Simulations
• Hamiltonian reductions of the one-dimensional Vlasov equation using phase-space moments
• Charge-conserving, variational particle-in-cell method for the drift-kinetic Vlasov-Maxwell system
• On Geometric Fourier Particle In Cell Methods
• A Vlasov Algorithm Derived from Phase Space Conservation
• Symplectic Modeling of Beam Loading in Electromagnetic Cavities
• Quasi Monte Carlo inverse transform sampling for phase space conserving Lagrangian particle methods and Eulerian-Lagrangian coupling