Exact charge-conserving scatter-gather algorithm for particle-in-cell simulations on unstructured grids: A geometric perspective

Abstract We describe a charge-conserving scatter–gather algorithm for particle-in-cell simulations on unstructured grids. Charge conservation is obtained from first principles, i.e., without the need for any post-processing or correction steps. This algorithm recovers, at a fundamental level, the scatter–gather algorithms presented recently by Campos-Pinto et al. (2014) (to first-order) and by Squire et al. (2012), but it is derived here in a streamlined fashion from a geometric viewpoint. Some ingredients reflecting this viewpoint are (1) the use of (discrete) differential forms of various degrees to represent fields, currents, and charged particles and provide localization rules for the degrees of freedom thereof on the various grid elements (nodes, edges, facets), (2) use of Whitney forms as basic interpolants from discrete differential forms to continuum space, and (3) use of a Galerkin formula for the discrete Hodge star operators (i.e., “mass matrices” incorporating the metric datum of the grid) applicable to generally irregular, unstructured grids. The expressions obtained for the scatter charges and scatter currents are very concise and do not involve numerical quadrature rules. Appropriate fractional areas within each grid element are identified that represent scatter charges and scatter currents within the element, and a simple geometric representation for the (exact) charge conservation mechanism is obtained by such identification. The field update is based on the coupled first-order Maxwell’s curl equations to avoid spurious modes with secular growth (otherwise present in formulations that discretize the second-order wave equation). Examples are provided to verify preservation of discrete Gauss’ law for all times.

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

[2]  M. Clemens,et al.  Discrete Electromagnetism With the Finite Integration Technique - Abstract , 2001 .

[3]  T. Esirkepov,et al.  Exact charge conservation scheme for Particle-in-Cell simulation with an arbitrary form-factor , 2001 .

[4]  Lie-Quan Lee,et al.  PARALLEL HIGHER-ORDER FINITE ELEMENT METHOD FOR ACCURATE FIELD COMPUTATIONS IN WAKEFIELD AND PIC SIMULATIONS , 2009 .

[5]  F. Teixeira,et al.  On the Degrees of Freedom of Lattice Electrodynamics , 2008 .

[6]  Carretera de Valencia,et al.  The finite element method in electromagnetics , 2000 .

[7]  Michael C. Huang,et al.  Particle-in-cell simulations with charge-conserving current deposition on graphic processing units , 2010, J. Comput. Phys..

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

[9]  Hiroshi Matsumoto,et al.  A new charge conservation method in electromagnetic particle-in-cell simulations , 2003 .

[10]  A. Bruce Langdon,et al.  On enforcing Gauss' law in electromagnetic particle-in-cell codes , 1992 .

[11]  B. M. Marder,et al.  A method for incorporating Gauss' lasw into electromagnetic pic codes , 1987 .

[12]  J. W. Eastwood,et al.  The virtual particle electromagnetic particle-mesh method , 1991 .

[13]  Fernando L. Teixeira,et al.  Lattice Maxwell's Equations , 2014 .

[14]  Brian Osserman.  DIFFERENTIAL FORMS , 2005 .

[15]  H. Flanders.  Differential Forms with Applications to the Physical Sciences , 1964 .

[16]  A. Bossavit.  Computational Electromagnetism: Variational Formulations, Complementarity, Edge Elements , 1997 .

[17]  Rolf Schuhmann,et al.  Conservation of Discrete Energy and Related Laws in the Finite Integration Technique , 2001 .

[18]  Z. J. Cendes,et al.  A generalized tree-cotree gauge for magnetic field computation , 1995 .

[19]  Fernando L. Teixeira,et al.  Lattice Maxwell's Equations (Invited Paper) , 2014 .

[20]  John P. Verboncoeur,et al.  Charge conservation in electromagnetic PIC codes; spectral comparison of Boris/DADI and Langdon-Marder methods , 1997 .

[21]  F. Teixeira,et al.  Parallel and Explicit Finite-Element Time-Domain Method for Maxwell's Equations , 2011, IEEE Transactions on Antennas and Propagation.

[22]  F. Teixeira,et al.  Geometric finite element discretization of Maxwell equations in primal and dual spaces , 2005, physics/0503013.

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

[24]  Peter P. Silvester,et al.  Finite elements for electrical engineers: Finite Elements for Electrical Engineers , 1996 .

[25]  F. Teixeira,et al.  Differential Forms, Galerkin Duality, and Sparse Inverse Approximations in Finite Element Solutions of Maxwell Equations , 2007, IEEE Transactions on Antennas and Propagation.

[26]  R. Lee.  A note on mass lumping in the finite element time domain method , 2006, IEEE Transactions on Antennas and Propagation.

[27]  H. Whitney.  Geometric Integration Theory , 1957 .

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

[29]  G. Deschamps.  Electromagnetics and differential forms , 1981, Proceedings of the IEEE.

[30]  W. Chew,et al.  Lattice electromagnetic theory from a topological viewpoint , 1999 .

[31]  Ruey-Beei Wu,et al.  Treating late-time instability of hybrid finite-element/finite-difference time-domain method , 1999 .

[32]  Alain Bossavit,et al.  Generating Whitney forms of polynomial degree one and higher , 2002 .

[33]  R. Albanese,et al.  Integral formulation for 3D eddy-current computation using edge elements , 1988 .

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

[35]  John P. Verboncoeur,et al.  Particle-in-cell simulations of plasma accelerators and electron-neutral collisions , 2001 .

[36]  Fernando L. Teixeira,et al.  Trade-Offs for Unconditional Stability in the Finite-Element Time-Domain Method , 2014, IEEE Microwave and Wireless Components Letters.

[37]  Karl F. Warnick,et al.  Teaching Electromagnetic Field Theory Using Differential Forms , 1997, Teaching Electromagnetics.

[38]  A. Candel,et al.  Parallel finite element particle-in-cell code for simulations of space-charge dominated beam-cavity interactions , 2007, 2007 IEEE Particle Accelerator Conference (PAC).

[39]  Igor V. Sokolov.  Alternating-order interpolation in a charge-conserving scheme for particle-in-cell simulations , 2013, Comput. Phys. Commun..

[40]  John H. Booske,et al.  Plasma physics and related challenges of millimeter-wave-to-terahertz and high power microwave generationa) , 2008 .

[41]  A. Bossavit.  Whitney forms: a class of finite elements for three-dimensional computations in electromagnetism , 1988 .

[42]  D. Arnold,et al.  Finite element exterior calculus, homological techniques, and applications , 2006, Acta Numerica.

[43]  Fernando L. Teixeira,et al.  Differential Forms in Lattice Field Theories: An Overview , 2013 .

[44]  R W Hockney,et al.  Computer Simulation Using Particles , 1966 .

[45]  D. J. Larson,et al.  Fast-ignition transport studies: Realistic electron source, integrated particle-in-cell and hydrodynamic modeling, imposed magnetic fields , 2012, 1205.1594.

[46]  J. S. Shang,et al.  Computational electromagnetics , 1996, CSUR.

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