Higher-order particle representation for particle-in-cell simulations

Abstract In this paper we present an alternative approach to the representation of simulation particles for unstructured electrostatic and electromagnetic PIC simulations. In our modified PIC algorithm we represent particles as having a smooth shape function limited by some specified finite radius, r 0 . A unique feature of our approach is the representation of this shape by surrounding simulation particles with a set of virtual particles with delta shape, with fixed offsets and weights derived from Gaussian quadrature rules and the value of r 0 . As the virtual particles are purely computational, they provide the additional benefit of increasing the arithmetic intensity of traditionally memory bound particle kernels. The modified algorithm is implemented within Sandia National Laboratories' unstructured EMPIRE-PIC code, for electrostatic and electromagnetic simulations, using periodic boundary conditions. We show results for a representative set of benchmark problems, including electron orbit, a transverse electromagnetic wave propagating through a plasma, numerical heating, and a plasma slab expansion. Good error reduction across all of the chosen problems is achieved as the particles are made progressively smoother, with the optimal particle radius appearing to be problem-dependent.

[1]  A. B. Langdon,et al.  Performance and optimization of direct implicit particle simulation , 1989 .

[2]  Liang Wang,et al.  The Plasma Simulation Code: A modern particle-in-cell code with patch-based load-balancing , 2016, J. Comput. Phys..

[3]  Hong Qin,et al.  Geometric integration of the Vlasov-Maxwell system with a variational particle-in-cell scheme , 2012, 1401.6723.

[4]  T. H. Dupree.  KINETIC THEORY OF PLASMA AND THE ELECTROMAGNETIC FIELD , 1963 .

[5]  E. Quataert,et al.  PARTICLE-IN-CELL SIMULATIONS OF CONTINUOUSLY DRIVEN MIRROR AND ION CYCLOTRON INSTABILITIES IN HIGH BETA ASTROPHYSICAL AND HELIOSPHERIC PLASMAS , 2014, 1402.0014.

[6]  M. Berger,et al.  Adaptive mesh refinement for hyperbolic partial differential equations , 1982 .

[7]  Jan S. Hesthaven,et al.  High-order nodal discontinuous Galerkin particle-in-cell method on unstructured grids , 2006, J. Comput. Phys..

[8]  R. Hockney,et al.  Quiet high resolution computer models of a plasma , 1974 .

[9]  Francis F. Chen,et al.  Introduction to Plasma Physics and Controlled Fusion , 2015 .

[10]  J. Dawson.  Particle simulation of plasmas , 1983 .

[11]  Wei Lu,et al.  OSIRIS: A Three-Dimensional, Fully Relativistic Particle in Cell Code for Modeling Plasma Based Accelerators , 2002, International Conference on Computational Science.

[12]  W. Zhang,et al.  Warp-X: A new exascale computing platform for beam–plasma simulations , 2017, Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment.

[13]  Rony Keppens,et al.  A Comprehensive Comparison of Relativistic Particle Integrators , 2017, 1710.09164.

[14]  Qing Zhang,et al.  The Gyrokinetic Particle Simulation of Fusion Plasmas on Tianhe-2 Supercomputer , 2016, 2016 7th Workshop on Latest Advances in Scalable Algorithms for Large-Scale Systems (ScalA).

[15]  Benjamin Bergen,et al.  0.374 Pflop/s trillion-particle kinetic modeling of laser plasma interaction on roadrunner , 2008, 2008 SC - International Conference for High Performance Computing, Networking, Storage and Analysis.

[16]  Phillip Colella,et al.  Application of adaptive mesh refinement to particle-in-cell simulations of plasmas and beams , 2003 .

[17]  A. Bruce Langdon,et al.  THEORY OF PLASMA SIMULATION USING FINITE-SIZE PARTICLES. , 1970 .

[18]  P. W. Rambo,et al.  Numerical Heating in Hybrid Plasma Simulations , 1997 .

[19]  John D. Villasenor,et al.  Rigorous charge conservation for local electromagnetic field solvers , 1992 .

[20]  R. Hockney,et al.  Measurements of collision and heating times in a two-dimensional thermal computer plasma , 1971 .

[21]  Claus-Dieter Munz,et al.  Comparison of coupling techniques in a high-order discontinuous Galerkin-based particle-in-cell solver , 2011 .

[22]  Robert Bridson,et al.  A high‐order accurate particle‐in‐cell method , 2012 .

[23]  C. Birdsall,et al.  Plasma Physics via Computer Simulation , 2018 .

[24]  Martin Campos Pinto,et al.  Charge-conserving FEM–PIC schemes on general grids☆ , 2014 .

[25]  Rolf Schuhmann,et al.  Conformal FDTD-methods to avoid time step reduction with and without cell enlargement , 2007, J. Comput. Phys..

[26]  Stephen A. Jarvis,et al.  Performance of a Second Order Electrostatic Particle-in-Cell Algorithm on Modern Many-Core Architectures , 2018, UKPEW.

[27]  K. Bowers,et al.  Ultrahigh performance three-dimensional electromagnetic relativistic kinetic plasma simulationa) , 2008 .

[28]  R. Marchand,et al.  PTetra, a Tool to Simulate Low Orbit Satellite–Plasma Interaction , 2012, IEEE Transactions on Plasma Science.

[29]  Samuel Williams,et al.  Kinetic turbulence simulations at extreme scale on leadership-class systems , 2013, 2013 SC - International Conference for High Performance Computing, Networking, Storage and Analysis (SC).

[30]  D. Payan,et al.  SPIS Open-Source Code: Methods, Capabilities, Achievements, and Prospects , 2008, IEEE Transactions on Plasma Science.

[31]  R. G. Evans,et al.  Contemporary particle-in-cell approach to laser-plasma modelling , 2015 .

[32]  Derek M. Causon,et al.  A cartesian cut cell method for compressible flows Part A: static body problems , 1997, The Aeronautical Journal (1968).

[33]  J. Bittencourt.  Fundamentals of plasma physics , 1986 .

[34]  R. Mittra,et al.  A locally conformal finite-difference time-domain (FDTD) algorithm for modeling three-dimensional perfectly conducting objects , 1997 .

[35]  Gregory Fridman,et al.  Applied Plasma Medicine , 2008 .

[36]  Iu. L. Klimontovich,et al.  The statistical theory of non-equilibrium processes in a plasma , 1967 .

[37]  Fernando L. Teixeira,et al.  Exact charge-conserving scatter-gather algorithm for particle-in-cell simulations on unstructured grids: A geometric perspective , 2014, Comput. Phys. Commun..

[38]  Charles K. Birdsall,et al.  Clouds-in-Clouds, Clouds-in-Cells Physics for Many-Body Plasma Simulation , 1997 .

[39]  K. Yee.  Numerical solution of initial boundary value problems involving maxwell's equations in isotropic media , 1966 .

[40]  J. Hesthaven,et al.  Nodal high-order methods on unstructured grids , 2002 .

[41]  J. Nédélec.  Mixed finite elements in ℝ3 , 1980 .

[42]  A. Lamberts,et al.  SHARP: A Spatially Higher-order, Relativistic Particle-in-cell Code , 2017, 1702.04732.

[43]  Leonid Oliker,et al.  Extreme Scale Plasma Turbulence Simulations on Top Supercomputers Worldwide , 2016, SC16: International Conference for High Performance Computing, Networking, Storage and Analysis.

[44]  A. Bruce Langdon,et al.  EFFECTS OF THE SPATIAL GRID IN SIMULATION PLASMAS. , 1970 .

[45]  T. D. Pointon.  Second-order, exact charge conservation for electromagnetic particle-in-cell simulation in complex geometry , 2008, Comput. Phys. Commun..

[46]  Jian-Ming Jin,et al.  The Finite Element Method in Electromagnetics , 1993 .