3D electromagnetic monte carlo particle-in-cell simulations on MIMD parallel computers

摘要

A three-dimensional electromagnetic particle-in-cell code with Monte Carlo collision (PIC-MCC) has been developed for MIMD parallel supercomputers. This code uses a standard relativistic leapfrog scheme incorporated with Monte Carlo calculations to push plasma particles and to calculate collisional effects on particle orbits. A local finite-difference time-domain method is used to update the self-consistent electromagnetic fields. The code is implemented using the General Concurrent PIC (GCPIC) algorithm which uses a domain decom-position to divide the computation among the processors. Particles must be exchanged between processors as they move among subdomains. We evaluate the performance of this code using a 512-processor Intel Touch-stone Delta, a 512-processor Intel Paragon, and a 256-processor Cray T3D. It is shown that our code runs with a high parallel efficiency of > 95