相关论文: Numerical stability of relativistic beam multidime…
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…
The Particle-In-Cell (PIC) method is widely used in relativistic particle beam and laser plasma modeling. However, the PIC method exhibits numerical instabilities that can render unphysical simulation results or even destroy the simulation.…
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…
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…
Numerical heating in particle-in-cell (PIC) codes currently precludes the accurate simulation of cold, relativistic plasma over long periods, severely limiting their applications in astrophysical environments. We present a spatially…
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…
Many astrophysical plasmas are prone to beam-plasma instabilities. For relativistic and dilute beams, the {\it spectral} support of the beam-plasma instabilities is narrow, i.e., the linearly unstable modes that grow with rates comparable…
Numerical Cherenkov radiation (NCR) or instability is a detrimental effect frequently found in electromagnetic particle-in-cell (EM-PIC) simulations involving relativistic plasma beams. NCR is caused by spurious coupling between…
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…
We present an algorithm for solving the linear dispersion relation in an inhomogeneous, magnetised, relativistic plasma. The method is a generalisation of a previously reported algorithm that was limited to the homogeneous case. The…
In clusters of galaxies, the specific entropy of intracluster plasma increases outwards. Nevertheless, a number of recent studies have shown that the intracluster medium is subject to buoyancy instabilities due to the effects of cosmic rays…
The WT scheme, a piecewise polynomial force interpolation scheme with time-step dependency, is proposed in this paper for relativistic particle-in-cell (PIC) simulations. The WT scheme removes the lowest order numerical Cherenkov…
The interaction of TeV blazars emitted gamma-rays with the extragalactic background photons gives rise to a relativistic beam of electron-positron ($e^- e^+$) pairs propagating through the intergalactic medium, producing a cascade through…
We discuss the irreversibility, nonlocality, and fluctuations, as well as the Lyapunov and hydrodynamic instabilities characterizing atomistic, smooth-particle, and finite-difference solutions of the two-dimensional Rayleigh-B\'enard…
We explore how inhomogeneity in the background plasma number density alters the growth of electrostatic unstable wavemodes of beam plasma systems. This is particularly interesting for blazar-driven beam-plasma instabilities, which may be…
A Particle-in-Cell (PIC) numerical simulation of the electron Weibel instability is applied in a frame of Darwin (radiationless) approximation of the self-consistent fields of sparse plasma. As a result, we were able to supplement the…
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…
We describe an algorithm that computes the linear dispersion relation of waves and instabilities in relativistic plasmas within a Vlasov-Maxwell description. The method used is fully relativistic and involves explicit integration of…
We studied the stability property of numerical Cherenkov radiation in relativistic plasma flows employing particle-in-cell simulations. Using the implicit finite-difference time-domain method to solve Maxwell equations, we found that…
The filamentation instability of counterpropagating symmetric beams of electrons is examined with 1D and 2D particle-in-cell (PIC) simulations, which are oriented orthogonally to the beam velocity vector. The beams are uniform, warm and…