Monte Carlo particle-in-cell methods for the simulation of the Vlasov–Maxwell gyrokinetic equations
Open Access
- 13 July 2015
- journal article
- research article
- Published by Cambridge University Press (CUP) in Journal of Plasma Physics
- Vol. 81 (5)
- https://doi.org/10.1017/s0022377815000574
Abstract
The particle-in-cell (PIC) algorithm is the most popular method for the discretisation of the general 6D Vlasov–Maxwell problem and it is widely used also for the simulation of the 5D gyrokinetic equations. The method consists of coupling a particle-based algorithm for the Vlasov equation with a grid-based method for the computation of the self-consistent electromagnetic fields. In this review we derive a Monte Carlo PIC finite-element model starting from a gyrokinetic discrete Lagrangian. The variations of the Lagrangian are used to obtain the time-continuous equations of motion for the particles and the finite-element approximation of the field equations. The Noether theorem for the semi-discretised system implies a certain number of conservation properties for the final set of equations. Moreover, the PIC method can be interpreted as a probabilistic Monte Carlo like method, consisting of calculating integrals of the continuous distribution function using a finite set of discrete markers. The nonlinear interactions along with numerical errors introduce random effects after some time. Therefore, the same tools for error analysis and error reduction used in Monte Carlo numerical methods can be applied to PIC simulations.This publication has 42 references indexed in Scilit:
- Variational formulation of particle algorithms for kinetic plasma simulationsJournal of Computational Physics, 2013
- Synergy between ion temperature gradient turbulence and neoclassical processes in global gyrokinetic particle-in-cell simulationsPhysics of Plasmas, 2012
- Derivation via free energy conservation constraints of gyrofluid equations with finite-gyroradius electromagnetic nonlinearitiesPhysics of Plasmas, 2010
- Nonequilibrium gyrokinetic fluctuation theory and sampling noise in gyrokinetic particle-in-cell simulationsPhysics of Plasmas, 2007
- A global collisionless PIC code in magnetic coordinatesComputer Physics Communications, 2007
- Nonlinear low noise particle-in-cell simulations of electron temperature gradient driven turbulencePhysics of Plasmas, 2007
- Discrete particle noise in particle-in-cell simulations of plasma microturbulencePhysics of Plasmas, 2005
- A revised δf algorithm for nonlinear PIC simulationComputer Physics Communications, 2003
- A unified Monte Carlo interpretation of particle simulations and applications to non-neutral plasmasPhysics of Plasmas, 1994
- Theory of plasma transport in toroidal confinement systemsReviews of Modern Physics, 1976