Numerical stability analysis of the pseudo-spectral analytical time-domain PIC algorithm

DOI10.1016/J.JCP.2013.10.053MaRDI QIDQ348631zbMATH OpenOpenAlexFDO

Authors Brendan B. Godfrey, Jean-Luc Vay, Irving Haber

Publication date 5 December 2016

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

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

zbMATH Keywords

numerical stabilityparticle-in-cellpseudo-spectralrelativistic beam

Mathematics Subject Classification ID

Spectral, collocation and related methods for initial value and initial-boundary value problems involving PDEs (65M70) Probabilistic methods, particle methods, etc. for initial value and initial-boundary value problems involving PDEs (65M75) Statistical mechanics of plasmas (82D10)

Abstract: The pseudo-spectral analytical time-domain (PSATD) particle-in-cell (PIC) algorithm solves the vacuum Maxwell's equations exactly, has no Courant time-step limit (as conventionally defined), and offers substantial flexibility in plasma and particle beam simulations. It is, however, not free of the usual numerical instabilities, including the numerical Cherenkov instability, when applied to relativistic beam simulations. This paper derives and solves the numerical dispersion relation for the PSATD algorithm and compares the results with corresponding behavior of the more conventional pseudo-spectral time-domain (PSTD) and finite difference time-domain (FDTD) algorithms. In general, PSATD offers superior stability properties over a reasonable range of time steps. More importantly, one version of the PSATD algorithm, when combined with digital filtering, is almost completely free of the numerical Cherenkov instability for time steps (scaled to the speed of light) comparable to or smaller than the axial cell size.

Recommendations

• Numerical stability of relativistic beam multidimensional PIC simulations employing the Esirkepov algorithm
• On the numerical dispersion of electromagnetic particle-in-cell code: finite grid instability
• A spectral, quasi-cylindrical and dispersion-free particle-in-cell algorithm
• Numerical instability due to relativistic plasma drift in EM-PIC simulations
• Controlling the numerical Cerenkov instability in PIC simulations using a customized finite difference Maxwell solver and a local FFT based current correction

Cites work

• A domain decomposition method for pseudo-spectral electromagnetic simulations of plasmas
• Exact charge conservation scheme for particle-in-cell simulation with an arbitrary form-factor
• High-accuracy Yee algorithm based on nonstandard finite differences: new developments and verifications
• scientific article; zbMATH DE number 486467 (Why is no real title available?)
• Numerical instability due to relativistic plasma drift in EM-PIC simulations
• Numerical solution of initial boundary value problems involving maxwell's equations in isotropic media
• Numerical stability of relativistic beam multidimensional PIC simulations employing the Esirkepov algorithm
• On the elimination of numerical Čerenkov radiation in PIC simulations

Cited in

(23)

• A particle-in-cell method for the simulation of plasmas based on an unconditionally stable field solver
• A spectral, quasi-cylindrical and dispersion-free particle-in-cell algorithm
• Finite grid instability and spectral fidelity of the electrostatic particle-in-cell algorithm
• Controlling the numerical Cerenkov instability in PIC simulations using a customized finite difference Maxwell solver and a local FFT based current correction
• Integrating a ponderomotive guiding center algorithm into a quasi-static particle-in-cell code based on azimuthal mode decomposition
• Diagnosing numerical Cherenkov instabilities in relativistic plasma simulations based on general meshes
• Detailed analysis of the effects of stencil spatial variations with arbitrary high-order finite-difference Maxwell solver
• Mitigation of numerical Cerenkov radiation and instability using a hybrid finite difference-FFT Maxwell solver and a local charge conserving current deposit
• On the numerical dispersion of electromagnetic particle-in-cell code: finite grid instability
• Enabling Lorentz boosted frame particle-in-cell simulations of laser wakefield acceleration in quasi-3D geometry
• A hybrid nodal-staggered pseudo-spectral electromagnetic particle-in-cell method with finite-order centering
• Efficiency of the perfectly matched layer with high-order finite difference and pseudo-spectral Maxwell solvers
• Numerical stability of relativistic beam multidimensional PIC simulations employing the Esirkepov algorithm
• Suppressing the numerical cherenkov instability in FDTD PIC codes
• Long-Time Instability Analysis of Pseudospectral Time-Domain Method
• A new field solver for modeling of relativistic particle-laser interactions using the particle-in-cell algorithm
• A quasi-static particle-in-cell algorithm based on an azimuthal Fourier decomposition for highly efficient simulations of plasma-based acceleration: QPAD
• A pseudospectral implicit particle-in-cell method with exact energy and charge conservation
• On variational Fourier particle methods
• High-order exponential solver method for particle-in-cell simulations
• A HIP-Compatible Accelerator Backend for Fourier-Bessel Particle-in-Cell Simulations on CPU/DCU Heterogeneous Clusters
• Time-step dependent force interpolation scheme for suppressing numerical Cherenkov instability in relativistic particle-in-cell simulations
• On numerical errors to the fields surrounding a relativistically moving particle in PIC codes

Describes a project that uses

Uses Software

• WARP

This page was built for publication: Numerical stability analysis of the pseudo-spectral analytical time-domain PIC algorithm

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