2024/06/29 by Sriramkrishnan Muralikrishnan, Muralikrishnan, Sriramkrishnan, Robert Speck +1
Engineering · Environmental Science · Materials Science · #35Q83 #65M75 #82D10 #Atmospheric aerosols and clouds #Computational Physics (physics.comp-ph) #Electron and X-Ray Spectroscopy Techniques #FOS: Mathematics #FOS: Physical sciences #Numerical Analysis (math.NA) #Plasma Physics (physics.plasm-ph) #Surface Roughness and Optical Measurements
paper · pdf · doi:10.48550/arxiv.2407.00485
openalex publication_date 2024/06/29 · openalex created_date 2024/07/04 · openalex updated_date 2026/07/28
We propose a parareal based time parallelization scheme in the phase-space for the particle-in-Fourier (PIF) discretization of the Vlasov-Poisson system used in kinetic plasma simulations. We use PIF with a coarse tolerance for the nonuniform fast Fourier transforms, or the standard particle-in-cell scheme, combined with temporal coarsening, as coarse propagators. This is different from the typical spatial coarsening of particles and/or Fourier modes for parareal, which are not possible or effective for PIF schemes. We perform an error analysis of the algorithm and verify the results numerically with Landau damping, two-stream instability, and Penning trap test cases in 3D-3V. We also implement the space-time parallelization of the PIF schemes in the open-source, performance-portable library IPPL and conduct scaling studies up to 1536 A100 GPUs on the JUWELS booster supercomputer. The space-time parallelization utilizing the parareal algorithm for the time parallelization provides up to 4-6 times speedup compared to spatial parallelization alone and achieves a push rate of around 1 billion particles per second for the benchmark plasma mini-apps considered.