Exportar Publicação

Yu, P., Xu, X., Tableman, A., Decyk, V. K., Tsung, F. S., Fiuza, F....Mori, W. B. (2015). Mitigation of numerical Cerenkov radiation and instability using a hybrid finite difference-FFT Maxwell solver and a local charge conserving current deposit. Computer Physics Communications. 197, 144-152

P. Yu et al.,  "Mitigation of numerical Cerenkov radiation and instability using a hybrid finite difference-FFT Maxwell solver and a local charge conserving current deposit", in Computer Physics Communications, vol. 197, pp. 144-152, 2015

@article{yu2015_1714843143357,
	author = "Yu, P. and Xu, X. and Tableman, A. and Decyk, V. K. and Tsung, F. S. and Fiuza, F. and Vieira, J. and Davidson, A. and Fonseca, R. A. and Lu, W. and Silva, L. O. and Mori, W. B.",
	title = "Mitigation of numerical Cerenkov radiation and instability using a hybrid finite difference-FFT Maxwell solver and a local charge conserving current deposit",
	journal = "Computer Physics Communications",
	year = "2015",
	volume = "197",
	number = "",
	doi = "10.1016/j.cpc.2015.08.026",
	pages = "144-152",
	url = "http://www.sciencedirect.com/science/article/pii/S0010465515003252"
}

TY  - JOUR
TI  - Mitigation of numerical Cerenkov radiation and instability using a hybrid finite difference-FFT Maxwell solver and a local charge conserving current deposit
T2  - Computer Physics Communications
VL  - 197
AU  - Yu, P.
AU  - Xu, X.
AU  - Tableman, A.
AU  - Decyk, V. K.
AU  - Tsung, F. S.
AU  - Fiuza, F.
AU  - Vieira, J.
AU  - Davidson, A.
AU  - Fonseca, R. A.
AU  - Lu, W.
AU  - Silva, L. O.
AU  - Mori, W. B.
PY  - 2015
SP  - 144-152
SN  - 0010-4655
DO  - 10.1016/j.cpc.2015.08.026
UR  - http://www.sciencedirect.com/science/article/pii/S0010465515003252
AB  - A hybrid Maxwell solver for fully relativistic and electromagnetic (EM) particle-in-cell (PIC) codes is described. In this solver, the EM fields are solved in k space by performing an FFT in one direction, while using finite difference operators in the other direction(s). This solver eliminates the numerical Cerenkov radiation for particles moving in the preferred direction. Moreover, the numerical Cerenkov instability (NCI) induced by the relativistically drifting plasma and beam can be eliminated using this hybrid solver by applying strategies that are similar to those recently developed for pure FFT solvers. A current correction is applied for the charge conserving current deposit to ensure that Gauss's Law is satisfied. A theoretical analysis of the dispersion properties in vacuum and in a drifting plasma for the hybrid solver is presented, and compared with PIC simulations with good agreement obtained. This hybrid solver is applied to both 2D and 3D Cartesian and quasi-3D (in which the fields and current are decomposed into azimuthal harmonics) geometries. Illustrative results for laser wakefield accelerator simulation in a Lorentz boosted frame using the hybrid solver in the 2D Cartesian geometry are presented, and compared against results from 2D UPIC-EMMA simulation which uses a pure spectral Maxwell solver, and from OSIRIS 2D lab frame simulation using the standard Yee solver. Very good agreement is obtained which demonstrates the feasibility of using the hybrid solver for high fidelity simulation of relativistically drifting plasma with no evidence of the numerical Cerenkov instability. 
ER  -