Particle-in-cell simulations with charge-conserving current deposition on graphic processing units

  • Authors:
  • Xianglong Kong;Michael C. Huang;Chuang Ren;Viktor K. Decyk

  • Affiliations:
  • Department of Mechanical Engineering, University of Rochester, Rochester, NY 14627, USA and Laboratory for Laser Energetics, University of Rochester, Rochester, NY 14627, USA;Department of Electrical and Computer Engineering, University of Rochester, Rochester, NY 14627, USA;Department of Mechanical Engineering, University of Rochester, Rochester, NY 14627, USA and Laboratory for Laser Energetics, University of Rochester, Rochester, NY 14627, USA and Department of Phy ...;Department of Physics and Astronomy, University of California Los Angeles, Los Angeles, CA 90095, USA

  • Venue:
  • Journal of Computational Physics
  • Year:
  • 2011

Quantified Score

Hi-index 31.45

Visualization

Abstract

We present an implementation of a 2D fully relativistic, electromagnetic particle-in-cell code, with charge-conserving current deposition, on parallel graphics processors (GPU) with CUDA. The GPU implementation achieved a one particle-step process time of 2.52ns for cold plasma runs and 9.15ns for extremely relativistic plasma runs, which are respectively 81 and 27 times faster than a single threaded state-of-art CPU code. A particle-based computation thread assignment was used in the current deposition scheme and write conflicts among the threads were resolved by a thread racing technique. A parallel particle sorting scheme was also developed and used. The implementation took advantage of fast on-chip shared memory, and can in principle be extended to 3D.