Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userElimination of Numerical Cherenkov Instability in flowing-plasma Particle-In-Cell simulations by using Galilean coordinates
87
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Particle-In-Cell (PIC) simulations of relativistic flowing plasmas are of key interest to several fields of physics (including e.g. laser-wakefield acceleration, when viewed in a Lorentz-boosted frame), but remain sometimes infeasible due to the well-known numerical Cherenkov instability (NCI). In this article, we show that, for a plasma drifting at a uniform relativistic velocity, the NCI can be eliminated by simply integrating the PIC equations in Galilean coordinates that follow the plasma (also sometimes known as comoving coordinates) within a spectral analytical framework. The elimination of the NCI is verified empirically and confirmed by a theoretical analysis of the instability. Moreover, it is shown that this method is applicable both to Cartesian geometry and to cylindrical geometry with azimuthal Fourier decomposition.
rate research
Read More
We design and develop a new Particle-in-Cell (PIC) method for plasma simulations using Deep-Learning (DL) to calculate the electric field from the electron phase space. We train a Multilayer Perceptron (MLP) and a Convolutional Neural Network (CNN) to solve the two-stream instability test. We verify that the DL-based MLP PIC method produces the correct results using the two-stream instability: the DL-based PIC provides the expected growth rate of the two-stream instability. The DL-based PIC does not conserve the total energy and momentum. However, the DL-based PIC method is stable against the cold-beam instability, affecting traditional PIC methods. This work shows that integrating DL technologies into traditional computational methods is a viable approach for developing next-generation PIC algorithms.
A procedure for largely suppressing the numerical Cherenkov instability in finite difference time-domain (FDTD) particle-in-cell (PIC) simulations of cold, relativistic beams is derived, and residual growth rates computed and compared with WARP code simulation results. Sample laser-plasma acceleration simulation output is provided to further validate the new procedure.
Several simulations of turbulence in the Large Plasma Device (LAPD) [W. Gekelman et al., Rev. Sci. Inst. 62, 2875 (1991)] are energetically analyzed and compared with each other and with the experiment. The simulations use the same model, but different axial boundary conditions. They employ either periodic, zero-value, zero-derivative, or sheath axial boundaries. The linear stability physics is different between the scenarios because the various boundary conditions allow the drift wave instability to access different axial structures, and the sheath boundary simulation contains a conducting wall mode instability which is just as unstable as the drift waves. Nevertheless, the turbulence in all the simulations is relatively similar because it is primarily driven by a robust nonlinear instability that is the same for all cases. The nonlinear instability preferentially drives $k_parallel = 0$ potential energy fluctuations, which then three-wave couple to $k_parallel e 0$ potential energy fluctuations in order to access the adiabatic response to transfer their energy to kinetic energy fluctuations. The turbulence self-organizes to drive this nonlinear instability, which destroys the linear eigenmode structures, making the linear instabilities ineffective.
When using an electromagnetic particle-in-cell (EM-PIC) code to simulate a relativistically drifting plasma, a violent numerical instability known as the numerical Cerenkov instability (NCI) occurs. The NCI is due to the unphysical coupling of electromagnetic waves on a grid to wave-particle resonances, including aliased resonances, i.e., $omega + 2pimu/Delta t=(k_1+ 2pi u_1/Delta x_1)v_0$, where $mu$ and $ u_1$ refer to the time and space aliases and the plasma is drifting relativistically at velocity $v_0$ in the $hat{1}$-direction. Recent studies have shown that an EM-PIC code which uses a spectral field solver and a low pass filter can eliminate the fastest growing modes of the NCI. Based on these studies a new spectral PIC code for studying laser wakefield acceleration (LWFA) in the Lorentz boosted frame was developed. However, we show that for parameters of relevance for LWFA simulations in the boosted frame, a relativistically drifting plasma is susceptible to a host of additional unstable modes with lower growth rates, and that these modes appear when the fastest growing unstable modes are filtered out. We show that these modes are most easily identified as the coupling between modes which are purely transverse (EM) and purely longitudinal (Langmuir) in the rest frame of the plasma for specific time and space aliases. We rewrite the dispersion relation of the drifting plasma for a general field solver and obtain analytic expressions for the location and growth rate for each unstable mode, i.e, for each time and space aliased resonances. We show for the spectral solver that when the fastest growing mode is eliminated a new mode at the fundamental resonance ($mu= u_1=0$) can be seen. (Please check the whole abstract in the paper).
The family of generalized Pseudo-Spectral Time Domain (including the Pseudo-Spectral Analytical Time Domain) Particle-in-Cell algorithms offers substantial versatility for simulating particle beams and plasmas, and well written codes using these algorithms run reasonably fast. When simulating relativistic beams and streaming plasmas in multiple dimensions, they are, however, subject to the numerical Cherenkov instability. Previous studies have shown that instability growth rates can be reduced substantially by modifying slightly the transverse fields as seen by the streaming particles . Here, we offer an approach which completely eliminates the fundamental mode of the numerical Cherenkov instability while minimizing the transverse field corrections. The procedure, numerically computed residual growth rates (from weaker, higher order instability aliases), and comparisons with WARP simulations are presented. In some instances, there are no numerical instabilities whatsoever, at least in the linear regime.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
