X-dispersionless Maxwell solver for plasma-based particle acceleration

JOURNAL OF COMPUTATIONAL PHYSICS

卷:418

Article

Pukhov, Alexander¹ 

[1] Univ Duesseldorf, Inst Theoret Phys 1, D-40225 Dusseldorf, Germany

关键词: Maxwell solver;    Dispersionless;    Numerical Cerenkov instability;   

DOI  :  10.1016/j.jcp.2020.109622

来源: Elsevier

【 摘 要 】

A semi-implicit finite difference time domain (FDTD) numerical Maxwell solver is developed for full electromagnetic Particle-in-Cell (PIC) codes for the simulations of plasma-based acceleration. The solver projects the volumetric Yee lattice into planes transverse to a selected axis (the particle acceleration direction). The scheme - by design - removes the numerical dispersion of electromagnetic waves running parallel the selected axis. The fields locations in the transverse plane are selected so that the scheme is Lorentz-invariant for relativistic transformations along the selected axis. The solver results in Galilean shift of transverse fields by exactly one cell per time step. This eases greatly the problem of numerical Cerenkov instability (NCI). The fields positions build rhombi in plane (RIP) patterns. The RIP scheme uses a compact local stencil that makes it perfectly suitable for massively parallel processing via domain decomposition along all three dimensions. No global/local spectral methods are involved. (C) 2020 Published by Elsevier Inc.

【 授权许可】

Free