Implementation of 2D Domain Decomposition in the UCAN Gyrokinetic Particle-in-Cell Code and Resulting Performance of UCAN2
Year: 2016
Communications in Computational Physics, Vol. 19 (2016), Iss. 1 : pp. 205–225
Abstract
The massively parallel, nonlinear, three-dimensional (3D), toroidal, electrostatic, gyrokinetic, particle-in-cell (PIC), Cartesian geometry UCAN code, with particle ions and adiabatic electrons, has been successfully exercised to identify non-diffusive transport characteristics in present day tokamak discharges. The limitation in applying UCAN to larger scale discharges is the 1D domain decomposition in the toroidal (or z-) direction for massively parallel implementation using MPI which has restricted the calculations to a few hundred ion Larmor radii or gyroradii per plasma minor radius. To exceed these sizes, we have implemented 2D domain decomposition in UCAN with the addition of the y-direction to the processor mix. This has been facilitated by use of relevant components in the P2LIB library of field and particle management routines developed for UCLA's UPIC Framework of conventional PIC codes. The gyro-averaging specific to gyrokinetic codes is simplified by the use of replicated arrays for efficient charge accumulation and force deposition. The 2D domain-decomposed UCAN2 code reproduces the original 1D domain nonlinear results within round-off. Benchmarks of UCAN2 on the Cray XC30 Edison at NERSC demonstrate ideal scaling when problem size is increased along with processor number up to the largest power of 2 available, namely 131,072 processors. These particle weak scaling benchmarks also indicate that the 1 nanosecond per particle per time step and 1 TFlops barriers are easily broken by UCAN2 with 1 billion particles or more and 2000 or more processors.
You do not have full access to this article.
Already a Subscriber? Sign in as an individual or via your institution
Journal Article Details
Publisher Name: Global Science Press
Language: English
DOI: https://doi.org/10.4208/cicp.070115.030715a
Communications in Computational Physics, Vol. 19 (2016), Iss. 1 : pp. 205–225
Published online: 2016-01
AMS Subject Headings: Global Science Press
Copyright: COPYRIGHT: © Global Science Press
Pages: 21
-
A Primer on Complex Systems
Laboratory Plasmas: Dynamics of Transport Across Sheared Flows
Sánchez, Raúl | Newman, David2018
https://doi.org/10.1007/978-94-024-1229-1_9 [Citations: 0] -
Parallelization and implementation of multi-spin Monte Carlo simulation of 2D square Ising model using MPI and C++
Hassani, Dariush | Rafibakhsh, ShahnooshJournal of Theoretical and Applied Physics, Vol. 12 (2018), Iss. 3 P.199
https://doi.org/10.1007/s40094-018-0301-4 [Citations: 2]