Mitigation of numerical Cerenkov radiation and instability using a hybrid finite difference-FFT Maxwell solver and a local charge conserving current deposit

DOI10.1016/J.CPC.2015.08.026MaRDI QIDQ2374070zbMATH OpenOpenAlexDBLPWikidataFDO

Authors Peicheng Yu, Xinlu Xu, Adam Tableman, Viktor K. Decyk, Frank S. Tsung, Frederico Fiuza, Asher Davidson, Jorge Vieira, Ricardo A. Fonseca, Wei Lü, Luis O. Silva, Warren B. Mori

Publication date 14 December 2016

Published in Computer Physics Communications (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/1502.01376

zbMATH Keywords

relativistic plasma driftPIC simulationnumerical Cerenkov instabilityhybrid Maxwell solverquasi-3D algorithm

Mathematics Subject Classification ID

Probabilistic methods, particle methods, etc. for initial value and initial-boundary value problems involving PDEs (65M75) Statistical mechanics of plasmas (82D10)

Abstract: 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 correctly account for the EM calculation in hybrid Yee-FFT solver. 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.

Recommendations

• Controlling the numerical Cerenkov instability in PIC simulations using a customized finite difference Maxwell solver and a local FFT based current correction
• Suppressing the numerical Cherenkov radiation in the Yee numerical scheme
• On numerical errors to the fields surrounding a relativistically moving particle in PIC codes
• Numerical instability due to relativistic plasma drift in EM-PIC simulations
• On the elimination of numerical Čerenkov radiation in PIC simulations

Cites work

• A domain decomposition method for pseudo-spectral electromagnetic simulations of plasmas
• Exact charge conservation scheme for particle-in-cell simulation with an arbitrary form-factor
• Numerical instability due to relativistic plasma drift in EM-PIC simulations
• Numerical simulations of laser wakefield accelerators in optimal Lorentz frames
• Numerical solution of initial boundary value problems involving maxwell's equations in isotropic media
• Numerical stability analysis of the pseudo-spectral analytical time-domain PIC algorithm
• Numerical stability of relativistic beam multidimensional PIC simulations employing the Esirkepov algorithm
• Suppressing the numerical cherenkov instability in FDTD PIC codes

Cited in

(10)

• Accurately simulating nine-dimensional phase space of relativistic particles in strong fields
• Suppressing the numerical Cherenkov radiation in the Yee numerical scheme
• A pseudospectral implicit particle-in-cell method with exact energy and charge conservation
• Numerical instability due to relativistic plasma drift in EM-PIC simulations
• A new field solver for modeling of relativistic particle-laser interactions using the particle-in-cell algorithm
• On numerical errors to the fields surrounding a relativistically moving particle in PIC codes
• A systematic approach to numerical dispersion in Maxwell solvers
• A pseudo energy-invariant method for relativistic wave equations with Riesz space-fractional derivatives
• Controlling the numerical Cerenkov instability in PIC simulations using a customized finite difference Maxwell solver and a local FFT based current correction
• On the elimination of numerical Čerenkov radiation in PIC simulations

Describes a project that uses

Uses Software

• UPIC-EMMA
• OSIRIS

This page was built for publication: Mitigation of numerical Cerenkov radiation and instability using a hybrid finite difference-FFT Maxwell solver and a local charge conserving current deposit

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2374070)