Loading and injection of Maxwellian distributions in particle simulations

Existing and new particle loading and injection algorithms for particle simulations are analyzed to determine numerical accuracy and computational efficiency. Emphasis has been placed on loading and emission of Maxwellian, drifting Maxwellian, and cutoff Maxwellian velocity distributions. Once a velocity distribution has been inverted for loading or injection, time-centering of the position and velocity is necessary in order to maintain second-order accuracy. Here, the accuracy of these methods is determined and compared to three analytic test cases with spatially varying, time-dependent, and time-independent electric fields in a homogeneous magnetic field and a self-consistent crossed-field diode. The initial push is shown to be important in calculating the correct electric field at the boundary where particles are injected, in relaxing constraints on the time step, and in providing reliable field fluctuations due to particle statistics.

[1]  J. Byers,et al.  PERPENDICULARLY PROPAGATING PLASMA CYCLOTRON INSTABILITIES SIMULATED WITH A ONE-DIMENSIONAL COMPUTER MODEL. , 1970 .

[2]  V. Vahedi,et al.  Simultaneous Potential and Circuit Solution for Two-Dimensional Bounded Plasma Simulation Codes , 1997 .

[3]  C. Birdsall,et al.  Rapid Current Transition in a Crossed-Field Diode , 1996 .

[4]  George Marsaglia,et al.  Simulating interpolation search , 1993 .

[5]  Paul Bratley,et al.  Algorithm 659: Implementing Sobol's quasirandom sequence generator , 1988, TOMS.

[6]  W. Lawson.  Particle simulation of bounded 1D Plasma Systems , 1989 .

[7]  O. Buneman,et al.  Time-Reversible Difference Procedures , 1967 .

[8]  A. W. Hull.  The Effect of A Uniform Magnetic Field on the Motion of Electrons Between Coaxial Cylinders. , 1921 .

[9]  M. Lieberman,et al.  Verification of frequency scaling laws for capacitive radio‐frequency discharges using two‐dimensional simulations* , 1993 .

[10]  S. J. Gitomer,et al.  Multibeam instability in a Maxwellian simulation plasma , 1976 .

[11]  J. Denavit.  Discrete particle effects in whistler simulation , 1974 .

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

[13]  Jacques Denavit,et al.  Numerical simulation of plasmas with periodic smoothing in phase space , 1972 .

[14]  P. C. Gray,et al.  Heat flow between species in one-dimensional particle plasma simulations , 1991 .

[15]  S. J. Gitomer.  Comments on “Numerical Simulation of the Weibel Instability in One and Two Dimensions” , 1971 .

[16]  L. A. Schwager.  Effects of low reflection on the collector and source sheaths of a finite ion temperature plasma , 1988 .

[17]  C. Birdsall,et al.  Similarity of stability characteristics of planar and coaxial crossed‐field diodes , 1996 .

[18]  H. Okuda,et al.  Collisions in multi‐dimensional plasma simulations , 1975 .

[19]  John P. Verboncoeur,et al.  Multipactor electron discharge physics using an improved secondary emission model , 1998 .

[20]  Simulation of transmitted current in a cylindrical cross-field diode , 1995 .

[21]  John P. Verboncoeur,et al.  Simultaneous potential and circuit solution for 1D bounded plasma particle simulation codes , 1990 .

[22]  Y. Lau,et al.  Limiting current in a crossed‐field gap , 1993 .

[23]  Artificial cooling due to quiet injection in bounded plasma particle simulations , 1987 .

[24]  John M. Dawson,et al.  PLASMA OSCILLATIONS OF A LARGE NUMBER OF ELECTRON BEAMS , 1960 .

[25]  John P. Verboncoeur,et al.  An object-oriented electromagnetic PIC code , 1995 .

[26]  Y. Lau,et al.  Transition to turbulence in a crossed‐field gap , 1994 .