On the elimination of numerical Cerenkov radiation in PIC simulations

Greenwood, Andrew; Cartwright, Keith; Luginsland, J.W.; Baca, Ernest A.

doi:10.1016/j.jcp.2004.06.021

Cited by 104 publications

(83 citation statements)

References 5 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…For each species Np = 2 × 10 8 computational particles and tsim = 200 Tp with ∆t = 0.0325 Tp. For the P -mode, simulations performed using the standard Finite Difference Time Domain (FDTD) Maxwell solver algorithm of PICCANTE were strongly affected by numerical Čherenkov radiation (NCR; for details see Greenwood et al (2004)) due to high-frequency waves which propagate slower than high-energy particles. Thus, for the P -mode simulations we set up another PIC code (PICcolino) implementing a spectral Maxwell solver based on the Fast Fourier Transform, which is free from NCR.…”

Section: Numerical Set-upmentioning

confidence: 99%

Particle acceleration and radiation friction effects in the filamentation instability of pair plasmas

D'Angelo¹,

Fedeli²,

Sgattoni³

et al. 2015

Mon. Not. R. Astron. Soc.

View full text Add to dashboard Cite

The evolution of the filamentation instability produced by two counter-streaming, ultrarelativistic pair plasmas is studied with particle-in-cell simulations. Radiation friction effects are taken into account. Two dimensional simulations are performed for both cases of the initial momenta being perpendicular (T -mode) or parallel (P -mode) to the simulation plane. In the initial stage the instability is purely transverse for both modes and generates small-scale filaments which later merge into larger structures. Particle acceleration leads to a strong broadening of the energy spectrum with the formation of a peak at twice the initial energy for the T -mode. In the nonlinear stage significant differences between T -and P -modes in the evolution of the fields and in the spectra of accelerated particles are apparent. The presence of radiative losses does not change the dynamics of the instability but strongly affects the structure of the particle spectra in the ultra-relativistic regime (particle energy > 100 MeV) and for high plasma densities (> 10 21 cm −3 ).

show abstract

Section: Numerical Set-upmentioning

confidence: 99%

Particle acceleration and radiation friction effects in the filamentation instability of pair plasmas

D'Angelo¹,

Fedeli²,

Sgattoni³

et al. 2015

Mon. Not. R. Astron. Soc.

View full text Add to dashboard Cite

show abstract

“…It is also possible to alter the computational stencil used to calculate the spatial derivates when solving Maxwell's equations. 33 This allows a better control of the numerical dispersion introduced by the finite difference time domain method, and is relevant to avoid numerical Cherenkov radiation and to more precisely model the propagation of the laser over long distances.…”

Section: Numerical Algorithmsmentioning

confidence: 99%

Modeling laser wakefield accelerator experiments with ultrafast particle-in-cell simulations in boosted frames

Martins¹,

Fonseca²,

Vieira³

et al. 2010

Physics of Plasmas

View full text Add to dashboard Cite

The development of new laser systems at the 10 Petawatt range will push laser wakefield accelerators to novel regimes, for which theoretical scalings predict the possibility to accelerate electron bunches up to tens of GeVs in meter-scale plasmas. Numerical simulations will play a crucial role in testing, probing, and optimizing the physical parameters and the setup of future experiments. Fully kinetic simulations are computationally very demanding, pushing the limits of today’s supercomputers. In this paper, the recent developments in the OSIRIS framework [R. A. Fonseca et al., Lect. Notes Comput. Sci. 2331, 342 (2002)] are described, in particular the boosted frame scheme, which leads to a dramatic change in the computational resources required to model laser wakefield accelerators. Results from one-to-one modeling of the next generation of laser systems are discussed, including the confirmation of electron bunch acceleration to the energy frontier.

show abstract

“…Hence, the high energy particles can travel in vacuum faster than their own radiation. This effect is commonly referred to as numerical Cherenkov radiation [4], which (due to its accumulative character) corrupts the simulation. Hence, the electromagnetic field computation for short relativistic bunches in long structures remains a challenging problem even with the fastest computers available.…”

Section: Introductionmentioning

confidence: 99%

TE/TM field solver for particle beam simulations without numerical Cherenkov radiation

Zagorodnov

Weiland²

2005

Phys. Rev. ST Accel. Beams

View full text Add to dashboard Cite

The Yee finite-difference time domain method (FDTD) is commonly used in wake field and particle-incell simulations. However, in accelerator modeling the high energy particles can travel in vacuum faster than their own radiation. This effect is commonly referred to as numerical Cherenkov radiation and is a consequence of numerical grid dispersion. Several numerical approaches are proposed to reduce the dispersion for all angles and for a given frequency range, that justifies itself for domains big in all three directions. On the contrary, in accelerator modeling the transverse dimensions and transverse beam velocity are small, but it is extremely important to eliminate the dispersion error in the well-defined direction of the beam motion for all frequencies. In this paper we propose a new two-level economical conservative scheme for electromagnetic field calculations in three dimensions. The scheme does not have dispersion in the longitudinal direction and is staircase-free (second order convergent). Unlike the FDTD method, it is based on a ''transversal-electric/transversal-magnetic'' (TE/TM)-like splitting of the field components in time. The scheme assures energy and charge conservation. Additionally, the usage of damping terms allows suppressing high frequency noise generated due to the transverse dispersion and the current fluctuations. The dispersion relation of the damping scheme is analyzed. As numerical examples show, the new scheme is much more accurate on the long-time scale than the conventional FDTD approach.

show abstract

On the elimination of numerical Cerenkov radiation in PIC simulations

Cited by 104 publications

References 5 publications

Particle acceleration and radiation friction effects in the filamentation instability of pair plasmas

Particle acceleration and radiation friction effects in the filamentation instability of pair plasmas

Modeling laser wakefield accelerator experiments with ultrafast particle-in-cell simulations in boosted frames

TE/TM field solver for particle beam simulations without numerical Cherenkov radiation

Contact Info

Product

Resources

About