Exportar Publicação

Li, F., Yu, P. C., Xu, X., Fiuza, F., Decyk, V. K., Dalichaouch, T....Mori, W. B. (2017). Controlling the numerical Cerenkov instability in PIC simulations using a customized finite difference Maxwell solver and a local FFT based current correction. Computer Physics Communications. 214, 6-17

F. Li et al.,  "Controlling the numerical Cerenkov instability in PIC simulations using a customized finite difference Maxwell solver and a local FFT based current correction", in Computer Physics Communications, vol. 214, pp. 6-17, 2017

@article{li2017_1715118543860,
	author = "Li, F. and Yu, P. C. and Xu, X. and Fiuza, F. and Decyk, V. K. and Dalichaouch, T. and Davidson, A. and Tableman, A. and An, W. and Tsung, F. S. and Fonseca, R. A. and Lu, W. and Mori, W. B.",
	title = "Controlling the numerical Cerenkov instability in PIC simulations using a customized finite difference Maxwell solver and a local FFT based current correction",
	journal = "Computer Physics Communications",
	year = "2017",
	volume = "214",
	number = "",
	doi = "10.1016/j.cpc.2017.01.001",
	pages = "6-17",
	url = "http://www.sciencedirect.com/science/article/pii/S0010465517300012"
}

TY  - JOUR
TI  - Controlling the numerical Cerenkov instability in PIC simulations using a customized finite difference Maxwell solver and a local FFT based current correction
T2  - Computer Physics Communications
VL  - 214
AU  - Li, F.
AU  - Yu, P. C.
AU  - Xu, X.
AU  - Fiuza, F.
AU  - Decyk, V. K.
AU  - Dalichaouch, T.
AU  - Davidson, A.
AU  - Tableman, A.
AU  - An, W.
AU  - Tsung, F. S.
AU  - Fonseca, R. A.
AU  - Lu, W.
AU  - Mori, W. B.
PY  - 2017
SP  - 6-17
SN  - 0010-4655
DO  - 10.1016/j.cpc.2017.01.001
UR  - http://www.sciencedirect.com/science/article/pii/S0010465517300012
AB  - In this paper we present a customized finite-difference-time-domain (FDTD) Maxwell solver for the particle-in-cell (PIC) algorithm. The solver is customized to effectively eliminate the numerical Cerenkov instability (NCI) which arises when a plasma (neutral or non-neutral) relativistically drifts on a grid when using the PIC algorithm. We control the EM dispersion curve in the direction of the plasma drift of a FDTD Maxwell solver by using a customized higher order finite difference operator for the spatial derivative along the direction of the drift (1 direction). We show that this eliminates the main NCI modes with moderate broken vertical bar k(1)broken vertical bar, while keeps additional main NCI modes well outside the range of physical interest with higher broken vertical bar k(1)broken vertical bar. These main NCI modes can be easily filtered out along with first spatial aliasing NCI modes which are also at the edge of the fundamental Brillouin zone. The customized solver has the possible advantage of improved parallel scalability because it can be easily partitioned along (1) over bar which typically has many more cells than other directions for the problems of interest. We show that FFTs can be performed locally to current on each partition to filter out the main and first spatial aliasing NCI modes, and to correct the current so that it satisfies the continuity equation for the customized spatial derivative. This ensures that Gauss' Law is satisfied. We present simulation examples of one relativistically drifting plasma, of two colliding relativistically drifting plasmas, and of nonlinear laser wakefield acceleration (LWFA) in a Lorentz boosted frame that show no evidence of the NCI can be observed when using this customized Maxwell solver together with its NCI elimination scheme.
ER  -