Numerical stability analysis of the pseudo-spectral analytical time-domain PIC algorithm
From MaRDI portal
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
- scientific article; zbMATH DE number 486467 (Why is no real title available?)
- 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
- 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
(21)- Diagnosing numerical Cherenkov instabilities in relativistic plasma simulations based on general meshes
- Enabling Lorentz boosted frame particle-in-cell simulations of laser wakefield acceleration in quasi-3D geometry
- Mitigation of numerical Cerenkov radiation and instability using a hybrid finite difference-FFT Maxwell solver and a local charge conserving current deposit
- A particle-in-cell method for the simulation of plasmas based on an unconditionally stable field solver
- Detailed analysis of the effects of stencil spatial variations with arbitrary high-order finite-difference Maxwell solver
- A pseudospectral implicit particle-in-cell method with exact energy and charge conservation
- Efficiency of the perfectly matched layer with high-order finite difference and pseudo-spectral Maxwell solvers
- Long-Time Instability Analysis of Pseudospectral Time-Domain Method
- On variational Fourier particle methods
- A new field solver for modeling of relativistic particle-laser interactions using the particle-in-cell algorithm
- On numerical errors to the fields surrounding a relativistically moving particle in PIC codes
- Time-step dependent force interpolation scheme for suppressing numerical Cherenkov instability in relativistic particle-in-cell simulations
- Numerical stability of relativistic beam multidimensional PIC simulations employing the Esirkepov algorithm
- Suppressing the numerical cherenkov instability in FDTD PIC codes
- A quasi-static particle-in-cell algorithm based on an azimuthal Fourier decomposition for highly efficient simulations of plasma-based acceleration: QPAD
- 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
- A hybrid nodal-staggered pseudo-spectral electromagnetic particle-in-cell method with finite-order centering
- On the numerical dispersion of electromagnetic particle-in-cell code: finite grid instability
- Integrating a ponderomotive guiding center algorithm into a quasi-static particle-in-cell code based on azimuthal mode decomposition
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)