arrow
Return

Efficient Fourier basis particle simulation

delete2019-11-01
delete9
delete
OA
AI
M
Matthew S. Mitchell *
M
Matthew Miecnikowski
G
Gregory Beylkin
S
Scott Parker
DOI:10.1016/j.jcp.2019.07.023delete
deleteOriginal
deleteShare
deleteSave
View PDF
Abstract

Abstract

En 中文
The standard particle-in-cell algorithm suffers from grid heating. There exists a gridless alternative which bypasses the deposition step and calculates each Fourier mode of the charge density directly from the particle positions. We show that a gridless method can be computed efficiently through the use of an Unequally Spaced Fast Fourier Transform (USFFT) algorithm. After a spectral field solve, the forces on the particles are calculated via the inverse USFFT (a rapid solution of an approximate linear system) [1,2]. We provide one and two dimensional implementations of this algorithm with an asymptotic runtime of O(N-p + N-m(D) log N-m(D)) for each iteration, identical to the standard PIC algorithm (where N-p is the number of particles, N-m is the number of Fourier modes, and Dis the spatial dimensionality of the problem). We demonstrate superior energy conservation and reduced noise, as well as convergence of the energy conservation at small time steps. (C) 2019 Elsevier Inc. All rights reserved.
Keywords:
Numerical
Plasma
Particle-in-cell
Energy conserving
Momentum conserving
Fourier transform
AI Summary

AI Summary

Key information extracted from the uploaded paper, including a brief overview, abstract, background, key highlights, visual analysis, and future outlook.

Journal

Journal of Computational Physics cover
Journal of Computational Physics
IF:
3.8
Papers:
1.5W
Citations:
7.4W

Organization

University of Colorado System cover
University of Colorado System
Scholars:
6.3W
Papers: 5.5W
Citations: 1.8K
U
university of colorado boulder
Scholars:
1.9W
Papers: 1.5W
Citations: 33