Second-order, exact charge conservation for electromagnetic particle-in-cell simulation in complex geometry

Abstract A second-order, exact charge-conserving algorithm for accumulating charge and current on the spatial grid for electromagnetic particle-in-cell (EM-PIC) simulation in bounded geometry is presented. The algorithm supports standard EM-PIC exterior boundary conditions and complex internal conductors on non-uniform grids. Boundary surfaces are handled by smoothly transitioning from second to first-order weighting within half a cell of the boundary. When a particle is exactly on the boundary surface (either about to be killed, or just created), the weighting is fully first-order. This means that particle creation and particle/surface interaction models developed for first-order weighting do not need to be modified. An additional feature is the use of an energy-conserving interpolation scheme from the electric field on the grid to the particles. Results show that high-density, cold plasmas with ω p e Δ t ∼ 1 , and Δ x / λ D ≫ 1 , can be modeled with reasonable accuracy and good energy conservation. This opens up a significant new capability for explicit simulation of high-density plasmas in high-power devices.

[1]  R. Morse,et al.  NUMERICAL SIMULATION OF THE WEIBEL INSTABILITY IN ONE AND TWO DIMENSIONS. , 1971 .

[2]  T. D. Pointon,et al.  Slanted conducting boundaries and field emission of particles in an electromagnetic particle simulation code , 1991 .

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

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

[5]  W. Stygar,et al.  Particle-in-cell simulations of electron flow in the post-hole convolute of the Z accelerator , 2001 .

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

[7]  M. Cuneo.  The effect of electrode contamination, cleaning and conditioning on high-energy pulsed-power device performance , 1999 .

[8]  J. U. Brackbill,et al.  CELEST1D: an implicit, fully kinetic model for low-frequency, electromagnetic plasma simulation , 1992 .

[9]  R. Commisso,et al.  Evaluation of self-magnetically pinched diodes up to 10 MV as high-resolution flash X-ray sources , 2004, IEEE Transactions on Plasma Science.

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

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

[12]  N. A. Krall,et al.  Simulation codes for light-ion diode modeling , 1994 .

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

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

[15]  J. Denavit,et al.  Fluid and field algorithms for time-implicit plasma simulation , 1991 .

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

[17]  D. Hewett,et al.  Electromagnetic direct implicit plasma simulation , 1987 .

[18]  Ryohei Itatani,et al.  High Order Spline Interpolations in the Particle Simulation , 1985 .

[19]  John W. Luginsland,et al.  A virtual prototyping environment for directed-energy concepts , 2002, Comput. Sci. Eng..

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