Variational formulation of particle algorithms for kinetic plasma simulations

DOI10.1016/J.JCP.2013.03.006zbMATH Open1349.82078OpenAlexW2074709256MaRDI QIDQ347741FDOQ347741

Authors e.g. Evstatiev, B. A. Shadwick

Publication date 5 December 2016

Published in Journal of Computational Physics (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/1210.3743

zbMATH Keywords

plasmaenergy conservingkineticnumericalparticle-in-cellmomentum conservingvariational

Mathematics Subject Classification ID

Probabilistic methods, particle methods, etc. for initial value and initial-boundary value problems involving PDEs (65M75) Statistical mechanics of plasmas (82D10)

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.

Recommendations

• An energy- and charge-conserving, implicit, electrostatic particle-in-cell algorithm
• Application of the phase space action principle to finite-size particle plasma simulations in the drift-kinetic approximation
• A Particle-in-cell Method with Adaptive Phase-space Remapping for Kinetic Plasmas
• scientific article; zbMATH DE number 3922303
• The energy conserving particle-in-cell method

Cites work

• scientific article; zbMATH DE number 3945609 (Why is no real title available?)
• scientific article; zbMATH DE number 3720310 (Why is no real title available?)
• scientific article; zbMATH DE number 41859 (Why is no real title available?)
• scientific article; zbMATH DE number 3560976 (Why is no real title available?)
• A Lagrangian formulation of the Boltzmann-Vlasov equation for plasmas
• An energy- and charge-conserving, implicit, electrostatic particle-in-cell algorithm
• Classical Dynamics
• Discrete mechanics and variational integrators
• Energy-conserving numerical approximations for Vlasov plasmas
• HIGH CONVERGENCE ORDER FINITE ELEMENTS WITH LUMPED MASS MATRIX
• Hamiltonian reductions for modeling relativistic laser-plasma interactions
• Magnetohydrodynamics with implicit plasma simulation
• Mechanical integrators derived from a discrete variational principle
• Optimized Forest-Ruth- and Suzuki-like algorithms for integration of motion in many-body systems
• Plasma Oscillations of a Large Number of Electron Beams
• The Liouville theorem and accurate plasma simulation
• The energy conserving particle-in-cell method

Cited in

(37)

• A Particle-in-cell Method with Adaptive Phase-space Remapping for Kinetic Plasmas
• Energy and momentum conservation theorems for electrostatic simulations
• The plasma simulation code: a modern particle-in-cell code with patch-based load-balancing
• A fully implicit, conservative, non-linear, electromagnetic hybrid particle-ion/fluid-electron algorithm
• Variational framework for structure-preserving electromagnetic particle-in-cell methods
• Multi-dimensional, fully-implicit, spectral method for the Vlasov-Maxwell equations with exact conservation laws in discrete form
• Structure-preserving vs. standard particle-in-cell methods: the case of an electron hybrid model
• Simplex-in-cell technique for collisionless plasma simulations
• A pseudospectral implicit particle-in-cell method with exact energy and charge conservation
• High order maximum principle preserving semi-Lagrangian finite difference WENO schemes for the Vlasov equation
• Hamiltonian particle-in-cell methods for Vlasov-Poisson equations
• Adaptive symplectic model order reduction of parametric particle-based Vlasov–Poisson equation
• Partially linearized algorithms in gyrokinetic particle simulation
• Detection and prediction of equilibrium states in kinetic plasma simulations via mode tracking using reduced-order dynamic mode decomposition
• Highly efficient energy-conserving moment method for the multi-dimensional Vlasov-Maxwell system
• On the particle method for the Vlasov kinetic equation
• Efficient Fourier basis particle simulation
• Numerical simulations of one laser-plasma model based on Poisson structure
• A conservative semi-Lagrangian HWENO method for the Vlasov equation
• Variational formulation of classical and quantum models for intense laser pulse propagation
• On the velocity space discretization for the Vlasov-Poisson system: comparison between implicit Hermite spectral and particle-in-cell methods
• A particle-in-Fourier method with semi-discrete energy conservation for non-periodic boundary conditions
• On variational Fourier particle methods
• The development of an implicit full \(f\) method for electromagnetic particle simulations of Alfvén waves and energetic particle physics
• On the structure and statistical theory of turbulence of extended magnetohydrodynamics
• Charge-and-energy conserving moment-based accelerator for a multi-species Vlasov-Fokker-Planck-Ampère system, part I: Collisionless aspects
• Noise and error analysis and optimization in particle-based kinetic plasma simulations
• Hamiltonian reductions of the one-dimensional Vlasov equation using phase-space moments
• Variance reduced particle solution of the Fokker-Planck equation with application to rarefied gas and plasma dynamics
• A charge-momentum-energy-conserving 1D3V hybrid Lagrangian-Eulerian method for Vlasov-Maxwell system
• Application of the phase space action principle to finite-size particle plasma simulations in the drift-kinetic approximation
• Finite grid instability and spectral fidelity of the electrostatic particle-in-cell algorithm
• Quasi-Helmholtz decomposition, Gauss' laws and charge conservation for finite element particle-in-cell
• MHD-kinetic hybrid code based on structure-preserving finite elements with particles-in-cell
• Vlasov and Drift Kinetic Simulation Methods Based on the Symplectic Integrator
• Solving the Vlasov-Maxwell equations using Hamiltonian splitting
• An asymptotic-preserving and energy-conserving particle-in-cell method for Vlasov–Maxwell equations

This page was built for publication: Variational formulation of particle algorithms for kinetic plasma simulations

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q347741)