Numerical instability due to relativistic plasma drift in EM-PIC simulations
arXiv:1211.0953 · doi:10.1016/j.cpc.2013.07.003
Abstract
The numerical instability observed in the Electromagnetic-Particle-in-cell (EM-PIC) simulations with a plasma drifting with relativistic velocities is studied using both theory and computer simulations. We derive the numerical dispersion relation for a cold plasma drifting with a relativistic velocity and find an instability attributed to the coupling between the beam modes of the drifting plasma and the electromagnetic modes in the system. The characteristic pattern of the instability in Fourier space for various simulation setups and Maxwell Equation solvers are explored by solving the corresponding numerical dispersion relations. Furthermore, based upon these characteristic patterns we derive an asymptotic expression for the instability growth rate. The asymptotic expression greatly speeds up the calculation of instability growth rate and makes the parameter scan for minimal growth rate feasible even for full three dimensions. The results are compared against simulation results and good agreement is found. These results can be used as a guide to develop possible approaches to mitigate the instability. We examine the use of a spectral solver and show that such a solver when combined with a low pass filter with a cutoff value of essentially eliminates the instability while not modifying modes of physical interest. The use of spectral solver also provides minimal errors to electromagnetic modes in the lowest Brillouin zones.
Version 2 differs from 1 in that: in expressions of version 1 we inadvertently omitted a term (-1)^{sum_iν_i}, which was discovered upon seeing [ B. Godfrey and J. -L. Vay, arXiv: 1211.0232v1]; we rewrote the asymptotic expressions in 3D, and discussed the use of asymptotic expression to study the relationship between the instability growth rate, time step and grid sizes in 2D and 3D scenarios
References in corpus (5)
- Weibel-instability-mediated collisionless shocks in laboratory with ultraintense lasers
- Numerical stability of relativistic beam multidimensional PIC simulations employing the Esirkepov algorithm
- Suppressing the Numerical Cherenkov Instability in FDTD PIC Codes
- Elimination of the numerical Cerenkov instability for spectral EM-PIC codes
- Modeling of Laser wakefield acceleration in Lorentz boosted frame using EM-PIC code with spectral solver
Cited by in corpus (32)
- A spectral, quasi-cylindrical and dispersion-free Particle-In-Cell algorithm
- Implementation of a hybrid particle code with a PIC description in r-z and a gridless description in into OSIRIS
- High quality electron bunch generation using a longitudinal density-tailored plasma-based accelerator in the three-dimensional blowout regime
- Suppressing the Numerical Cherenkov Instability in FDTD PIC Codes
- Persistence of precursor waves in two-dimensional relativistic shocks
- Numerical stability analysis of the Pseudo-Spectral Analytical Time-Domain PIC algorithm
- Controlling the Numerical Cerenkov Instability in PIC simulations using a customized finite difference Maxwell solver and a local FFT based current correction
- Elimination of Numerical Cherenkov Instability in flowing-plasma Particle-In-Cell simulations by using Galilean coordinates
- Stable discrete representation of relativistically drifting plasmas
- Improved Numerical Cherenkov Instability Suppression in the Generalized PSTD PIC Algorithm
- A quasi-static particle-in-cell algorithm based on an azimuthal Fourier decomposition for highly efficient simulations of plasma-based acceleration: QPAD
- Elimination of the numerical Cerenkov instability for spectral EM-PIC codes
- Modeling of Laser wakefield acceleration in Lorentz boosted frame using EM-PIC code with spectral solver
- Mitigation of numerical Cerenkov radiation and instability using a hybrid finite difference-FFT Maxwell solver and a local charge conserving current deposit
- Finite grid instability and spectral fidelity of the electrostatic Particle-In-Cell algorithm
- Accurately simulating nine-dimensional phase space of relativistic particles in strong fields
- Stability Property of Numerical Cherenkov Radiation and its Application to Relativistic Shock Simulations
- On numerical errors to the fields surrounding a relativistically moving particle in PIC codes
- Ultra-Bright Electron Bunch Injection in a Plasma Wakefield Driven by a Superluminal Flying Focus Electron Beam
- On the Numerical Dispersion of Electromagnetic Particle-In-Cell Code : Finite Grid Instability
- Lorentz boosted frame simulation of Laser wakefield acceleration in quasi-3D geometry
- Time-step dependent force interpolation scheme for suppressing numerical Cherenkov instability in relativistic particle-in-cell simulations
- On the generation of ultra-bright and low energy spread electron beams in laser wakefield acceleration in a uniform plasma
- Accurate simulation of direct laser acceleration in a laser wakefield accelerator
- Overcoming timestep limitations in boosted-frame Particle-In-Cell simulations of plasma-based acceleration
- Integrating a ponderomotive guiding center algorithm into a quasi-static particle-in-cell code based on azimuthal mode decomposition
- Stopping Power Enhancement From Discrete Particle-Wake Correlations in High Energy Density Plasmas
- A new field solver for modeling of relativistic particle-laser interactions using the particle-in-cell algorithm
- An \mathcal{O}(N) Maxwell solver with improved numerical dispersion properties
- Modeling of Advanced Accelerator Concepts
- On Geometric Fourier Particle In Cell Methods
- A Hybrid Nodal-Staggered Pseudo-Spectral Electromagnetic Particle-In-Cell Method with Finite-Order Centering