Ultrahigh performance three-dimensional electromagnetic relativistic kinetic plasma simulationa)

The algorithms, implementation details, and applications of VPIC, a state-of-the-art first principles 3D electromagnetic relativistic kinetic particle-in-cell code, are discussed. Unlike most codes, VPIC is designed to minimize data motion, as, due to physical limitations (including the speed of light!), moving data between and even within modern microprocessors is more time consuming than performing computations. As a result, VPIC has achieved unprecedented levels of performance. For example, VPIC can perform ∼0.17 billion cold particles pushed and charge conserving accumulated per second per processor on IBM’s Cell microprocessor—equivalent to sustaining Los Alamos’s planned Roadrunner supercomputer at ∼0.56 petaflop (quadrillion floating point operations per second). VPIC has enabled previously intractable simulations in numerous areas of plasma physics, including magnetic reconnection and laser plasma interactions; next generation supercomputers like Roadrunner will enable further advances.

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

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

[3]  Jay P. Boris,et al.  Proceedings of the Conference on the Numerical Simulation of Plasmas (4th) Held at the Naval Research Laboratory, Washington, D.C. on 2, 3 November 1970 , 1971 .

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

[5]  R. Higdon.  Absorbing boundary conditions for difference approximations to the multi-dimensional wave equation , 1986 .

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

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

[8]  R. W. Hockney,et al.  Body-fitted electromagnetic PIC software for use on parallel computers , 1995 .

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

[10]  Dan Winske,et al.  Modeling ion interpenetration, stagnation, and thermalization in colliding plasmas , 1996 .

[11]  G. E. Sasser,et al.  3-D ICEPIC simulations of the relativistic klystron oscillator , 2000 .

[12]  Kevin J. Bowers,et al.  Accelerating a paricle -in-cell simulation using a hybrid counting sort , 2001 .

[13]  J. Cary,et al.  VORPAL: a versatile plasma simulation code , 2004 .

[14]  H. Peter Hofstee,et al.  Introduction to the Cell multiprocessor , 2005, IBM J. Res. Dev..

[15]  K. Bowers,et al.  Saturation of backward stimulated scattering of a laser beam in the kinetic regime. , 2007, Physical review letters.

[16]  Michael P Eastwood,et al.  A common, avoidable source of error in molecular dynamics integrators. , 2007, The Journal of chemical physics.

[17]  J. Vay,et al.  Noninvariance of space- and time-scale ranges under a Lorentz Transformation and the implications for the study of relativistic interactions. , 2007, Physical review letters.

[18]  K. Bowers,et al.  Spectral energy transfer and dissipation of magnetic energy from fluid to kinetic scales. , 2007, Physical review letters.

[19]  William Daughton,et al.  Fully Kinetic 3D Simulations of Collisionless Reconnection in Large-Scale Pair Plasmas , 2007 .

[20]  Brian James Albright,et al.  Monoenergetic and GeV ion acceleration from the laser breakout afterburner using ultrathin targets , 2007 .

[21]  Brian James Albright,et al.  Relativistic Buneman instability in the laser breakout afterburner , 2007 .