Variational Framework for Structure-Preserving Electromagnetic Particle-in-Cell Methods

In this article we apply a discrete action principle for the Vlasov-Maxwell equations in a structure-preserving particle-field discretization framework. In this framework the finite-dimensional electromagnetic potentials and fields are represented in a discrete de Rham sequence involving general finite element spaces, and the particle-field coupling is represented by a set of projection operators that commute with the differential operators. With a minimal number of assumptions which allow for a variety of finite elements and shape functions for the particles, we show that the resulting variational scheme has a general discrete Poisson structure and thus leads to a semi-discrete Hamiltonian system. By introducing discrete interior products we derive a second type of space discretization which is momentum preserving, based on the same finite elements and shape functions. We illustrate our method by applying it to spline finite elements, and to a new spectral discretization where the particle-field coupling relies on discrete Fourier transforms

Hierarchy of second order gyrokinetic Hamiltonian models for particle-in-cell codes

International audienceThe reduced-particle model is the central element for the systematic derivation of the gyrokinetic Vlasov-Maxwell equations from first principles. Coupled to the fields inside the gyrokinetic field-particle Lagrangian, the reduced-particle model defines polarization and magnetization effects appearing in the gyrokinetic Maxwell equations. It is also used for the reconstruction of the gyrokinetic Vlasov equation from the particle characteristics. Various representations of reduced-particle models are available according to the choice of the gyrokinetic phase space coordinates. In this paper, the Hamiltonian representation of the reduced particle dynamics at an order suitable for the implementation in particle-in-cell simulations is explicitly derived from the general reduction procedure. The second-order (with respect to the fluctuating electromagnetic fields), full Finite Larmor Radius (FLR) Hamiltonian gyrokinetic particle model as well as the second-order model suitable specifically for the long-wavelength approximation (i.e., containing up to the second-order FLR corrections), are derived and compared to the model recently implemented in the particle-in-cell code ORB5. We show that the same long-wavelength approximate equations can also be derived by taking the proper limit of the full FLR model