Efficient Fourier basis particle simulation

DOI10.1016/J.JCP.2019.07.023MaRDI QIDQ2222459zbMATH OpenOpenAlexWikidataFDO

Authors Matthew S. Mitchell, Matthew T. Miecnikowski, Gregory Beylkin, Scott E. Parker

Publication date 27 January 2021

Published in Journal of Computational Physics (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/1808.03742

zbMATH Keywords

energy conservingFourier transformmomentum conservingnumericalparticle-in-cellplasma

Mathematics Subject Classification ID

Vlasov equations (35Q83) Probabilistic methods, particle methods, etc. for initial value and initial-boundary value problems involving PDEs (65M75) Complexity and performance of numerical algorithms (65Y20)

Abstract: 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). 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 and $N_m$ is the number of Fourier modes, and $D$ is 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.

Recommendations

• A fast algorithm for particle simulations
• A fast algorithm for particle simulations
• Fast multipole methods for particle dynamics
• scientific article; zbMATH DE number 16957
• Fast exact evaluation of particle interaction vectors in the finite volume particle method
• Parallel three-dimensional nonequispaced fast Fourier transforms and their application to particle simulation
• A Fast Adaptive Multipole Algorithm for Particle Simulations
• On a practical implementation of particle methods
• An Effective Improved Algorithm for Finite Particle Method
• Particle methods for multi-scale and multi-physics

Cites work

• A Lagrangian formulation of the Boltzmann-Vlasov equation for plasmas
• A method to suppress the finite-grid instability in plasma simulations
• An energy- and charge-conserving, implicit, electrostatic particle-in-cell algorithm
• Energy-conserving numerical approximations for Vlasov plasmas
• Fast Fourier Transforms for Nonequispaced Data
• Finite grid instability and spectral fidelity of the electrostatic particle-in-cell algorithm
• High-order nodal discontinuous Galerkin particle-in-cell method on unstructured grids
• scientific article; zbMATH DE number 41859 (Why is no real title available?)
• On the fast Fourier transform of functions with singularities
• The energy conserving particle-in-cell method
• Variational formulation of particle algorithms for kinetic plasma simulations

Cited in

(10)

• The development of an implicit full \(f\) method for electromagnetic particle simulations of Alfvén waves and energetic particle physics
• Parallel three-dimensional nonequispaced fast Fourier transforms and their application to particle simulation
• Fast exact evaluation of particle interaction vectors in the finite volume particle method
• A pseudospectral implicit particle-in-cell method with exact energy and charge conservation
• Grid-free weighted particle method applied to the Vlasov-Poisson equation
• A particle-in-Fourier method with semi-discrete energy conservation for non-periodic boundary conditions
• On variational Fourier particle methods
• A Fourier spectral immersed boundary method with exact translation invariance, improved boundary resolution, and a divergence-free velocity field
• Structure preservation using discrete gradients in the Vlasov-Poisson-Landau system
• A fast spectral particle method for the Landau equation

This page was built for publication: Efficient Fourier basis particle simulation

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2222459)