Scalable spectral solver in Galilean coordinates for eliminating the numerical Cherenkov instability in particle-in-cell simulations of streaming plasmas
Abstract
Discretizing Maxwell's equations in Galilean (comoving) coordinates allows the derivation of a pseudospectral solver that eliminates the numerical Cherenkov instability for electromagnetic particle-in-cell simulations of relativistic plasmas flowing at a uniform velocity. Here we generalize this solver by incorporating spatial derivatives of arbitrary order, thereby enabling efficient parallelization by domain decomposition. This allows scaling of the algorithm to many distributed compute units. We derive the numerical dispersion relation of the algorithm and present a comprehensive theoretical stability analysis. The method is applied to simulations of plasma acceleration in a Lorentz-boosted frame of reference.
- Authors:
- Publication Date:
- Research Org.:
- Lawrence Berkeley National Laboratory (LBNL), Berkeley, CA (United States)
- Sponsoring Org.:
- USDOE National Nuclear Security Administration (NNSA); Federal Ministry of Education and Research (BMBF); USDOE Office of Science (SC), High Energy Physics (HEP)
- OSTI Identifier:
- 1637566
- Alternate Identifier(s):
- OSTI ID: 1756341
- Grant/Contract Number:
- AC02-05CH11231; 05K19GUA
- Resource Type:
- Published Article
- Journal Name:
- Physical Review. E
- Additional Journal Information:
- Journal Name: Physical Review. E Journal Volume: 102 Journal Issue: 1; Journal ID: ISSN 2470-0045
- Publisher:
- American Physical Society (APS)
- Country of Publication:
- United States
- Language:
- English
- Subject:
- 70 PLASMA PHYSICS AND FUSION TECHNOLOGY; Laser wakefield acceleration; Plasma acceleration & new acceleration techniques; Relativistic plasmas; Particle-in-cell methods
Citation Formats
Kirchen, Manuel, Lehe, Remi, Jalas, Soeren, Shapoval, Olga, Vay, Jean-Luc, and Maier, Andreas R. Scalable spectral solver in Galilean coordinates for eliminating the numerical Cherenkov instability in particle-in-cell simulations of streaming plasmas. United States: N. p., 2020.
Web. doi:10.1103/PhysRevE.102.013202.
Kirchen, Manuel, Lehe, Remi, Jalas, Soeren, Shapoval, Olga, Vay, Jean-Luc, & Maier, Andreas R. Scalable spectral solver in Galilean coordinates for eliminating the numerical Cherenkov instability in particle-in-cell simulations of streaming plasmas. United States. https://doi.org/10.1103/PhysRevE.102.013202
Kirchen, Manuel, Lehe, Remi, Jalas, Soeren, Shapoval, Olga, Vay, Jean-Luc, and Maier, Andreas R. Tue .
"Scalable spectral solver in Galilean coordinates for eliminating the numerical Cherenkov instability in particle-in-cell simulations of streaming plasmas". United States. https://doi.org/10.1103/PhysRevE.102.013202.
@article{osti_1637566,
title = {Scalable spectral solver in Galilean coordinates for eliminating the numerical Cherenkov instability in particle-in-cell simulations of streaming plasmas},
author = {Kirchen, Manuel and Lehe, Remi and Jalas, Soeren and Shapoval, Olga and Vay, Jean-Luc and Maier, Andreas R.},
abstractNote = {Discretizing Maxwell's equations in Galilean (comoving) coordinates allows the derivation of a pseudospectral solver that eliminates the numerical Cherenkov instability for electromagnetic particle-in-cell simulations of relativistic plasmas flowing at a uniform velocity. Here we generalize this solver by incorporating spatial derivatives of arbitrary order, thereby enabling efficient parallelization by domain decomposition. This allows scaling of the algorithm to many distributed compute units. We derive the numerical dispersion relation of the algorithm and present a comprehensive theoretical stability analysis. The method is applied to simulations of plasma acceleration in a Lorentz-boosted frame of reference.},
doi = {10.1103/PhysRevE.102.013202},
journal = {Physical Review. E},
number = 1,
volume = 102,
place = {United States},
year = {Tue Jul 07 00:00:00 EDT 2020},
month = {Tue Jul 07 00:00:00 EDT 2020}
}
https://doi.org/10.1103/PhysRevE.102.013202
Web of Science
Works referenced in this record:
Noninvariance of Space- and Time-Scale Ranges under a Lorentz Transformation and the Implications for the Study of Relativistic Interactions
journal, March 2007
- Vay, J. -L.
- Physical Review Letters, Vol. 98, Issue 13
Generation of finite difference formulas on arbitrarily spaced grids
journal, January 1988
- Fornberg, Bengt
- Mathematics of Computation, Vol. 51, Issue 184
Elimination of the numerical Cerenkov instability for spectral EM-PIC codes
journal, July 2015
- Yu, Peicheng; Xu, Xinlu; Decyk, Viktor K.
- Computer Physics Communications, Vol. 192
A domain decomposition method for pseudo-spectral electromagnetic simulations of plasmas
journal, June 2013
- Vay, Jean-Luc; Haber, Irving; Godfrey, Brendan B.
- Journal of Computational Physics, Vol. 243
Numerical Stability Improvements for the Pseudospectral EM PIC Algorithm
journal, May 2014
- Godfrey, Brendan B.; Vay, Jean-Luc; Haber, Irving
- IEEE Transactions on Plasma Science, Vol. 42, Issue 5
Numerical simulations of laser wakefield accelerators in optimal Lorentz frames
journal, May 2010
- Martins, Samuel F.; Fonseca, Ricardo A.; Silva, Luís O.
- Computer Physics Communications, Vol. 181, Issue 5
Accurate modeling of plasma acceleration with arbitrary order pseudo-spectral particle-in-cell methods
journal, March 2017
- Jalas, S.; Dornmair, I.; Lehe, R.
- Physics of Plasmas, Vol. 24, Issue 3
The PSTD algorithm: A time-domain method requiring only two cells per wavelength
journal, June 1997
- Liu, Q. H.
- Microwave and Optical Technology Letters, Vol. 15, Issue 3
On numerical errors to the fields surrounding a relativistically moving particle in PIC codes
journal, July 2020
- Xu, Xinlu; Li, Fei; Tsung, Frank S.
- Journal of Computational Physics, Vol. 413
Novel methods in the Particle-In-Cell accelerator Code-Framework Warp
journal, January 2012
- Vay, J-L; Grote, D. P.; Cohen, R. H.
- Computational Science & Discovery, Vol. 5, Issue 1
Elimination of numerical Cherenkov instability in flowing-plasma particle-in-cell simulations by using Galilean coordinates
journal, November 2016
- Lehe, Remi; Kirchen, Manuel; Godfrey, Brendan B.
- Physical Review E, Vol. 94, Issue 5
Suppressing the numerical Cherenkov instability in FDTD PIC codes
journal, June 2014
- Godfrey, Brendan B.; Vay, Jean-Luc
- Journal of Computational Physics, Vol. 267
Numerical growth of emittance in simulations of laser-wakefield acceleration
journal, February 2013
- Lehe, R.; Lifschitz, A.; Thaury, C.
- Physical Review Special Topics - Accelerators and Beams, Vol. 16, Issue 2
Mitigation of numerical Cerenkov radiation and instability using a hybrid finite difference-FFT Maxwell solver and a local charge conserving current deposit
journal, December 2015
- Yu, Peicheng; Xu, Xinlu; Tableman, Adam
- Computer Physics Communications, Vol. 197
Controlling the numerical Cerenkov instability in PIC simulations using a customized finite difference Maxwell solver and a local FFT based current correction
journal, May 2017
- Li, Fei; Yu, Peicheng; Xu, Xinlu
- Computer Physics Communications, Vol. 214
Particle Acceleration in Relativistic Collisionless Shocks: Fermi Process at Last?
journal, July 2008
- Spitkovsky, Anatoly
- The Astrophysical Journal, Vol. 682, Issue 1
Generalized algorithm for control of numerical dispersion in explicit time-domain electromagnetic simulations
journal, April 2013
- Cowan, Benjamin M.; Bruhwiler, David L.; Cary, John R.
- Physical Review Special Topics - Accelerators and Beams, Vol. 16, Issue 4
Numerical stability of relativistic beam multidimensional PIC simulations employing the Esirkepov algorithm
journal, September 2013
- Godfrey, Brendan B.; Vay, Jean-Luc
- Journal of Computational Physics, Vol. 248
Numerical stability analysis of the pseudo-spectral analytical time-domain PIC algorithm
journal, February 2014
- Godfrey, Brendan B.; Vay, Jean-Luc; Haber, Irving
- Journal of Computational Physics, Vol. 258
Stable discrete representation of relativistically drifting plasmas
journal, October 2016
- Kirchen, M.; Lehe, R.; Godfrey, B. B.
- Physics of Plasmas, Vol. 23, Issue 10
Warp-X: A new exascale computing platform for beam–plasma simulations
journal, November 2018
- Vay, J. -L.; Almgren, A.; Bell, J.
- Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, Vol. 909
A spectral, quasi-cylindrical and dispersion-free Particle-In-Cell algorithm
journal, June 2016
- Lehe, Rémi; Kirchen, Manuel; Andriyash, Igor A.
- Computer Physics Communications, Vol. 203
Detailed analysis of the effects of stencil spatial variations with arbitrary high-order finite-difference Maxwell solver
journal, March 2016
- Vincenti, H.; Vay, J. -L.
- Computer Physics Communications, Vol. 200
Canonical momenta and numerical instabilities in particle codes
journal, September 1975
- Godfrey, Brendan B.
- Journal of Computational Physics, Vol. 19, Issue 1
Numerical solution of initial boundary value problems involving maxwell's equations in isotropic media
journal, May 1966
- Kane Yee,
- IEEE Transactions on Antennas and Propagation, Vol. 14, Issue 3
Physics of laser-driven plasma-based electron accelerators
journal, August 2009
- Esarey, E.; Schroeder, C. B.; Leemans, W. P.
- Reviews of Modern Physics, Vol. 81, Issue 3
Numerical instability due to relativistic plasma drift in EM-PIC simulations
journal, November 2013
- Xu, Xinlu; Yu, Peicheng; Martins, Samual F.
- Computer Physics Communications, Vol. 184, Issue 11
Exploring laser-wakefield-accelerator regimes for near-term lasers using particle-in-cell simulation in Lorentz-boosted frames
journal, March 2010
- Martins, S. F.; Fonseca, R. A.; Lu, W.
- Nature Physics, Vol. 6, Issue 4
Interpenetrating Plasma Shells: Near-Equipartition Magnetic Field Generation and Nonthermal Particle Acceleration
journal, September 2003
- Silva, L. O.; Fonseca, R. A.; Tonge, J. W.
- The Astrophysical Journal, Vol. 596, Issue 1
A high-accuracy realization of the Yee algorithm using non-standard finite differences
journal, June 1997
- Cole, J. B.
- IEEE Transactions on Microwave Theory and Techniques, Vol. 45, Issue 6
Modeling of 10 GeV-1 TeV laser-plasma accelerators using Lorentz boosted simulations
journal, December 2011
- Vay, J. -L.; Geddes, C. G. R.; Esarey, E.
- Physics of Plasmas, Vol. 18, Issue 12
JUWELS: Modular Tier-0/1 Supercomputer at Jülich Supercomputing Centre
journal, February 2019
- Support, Supercomputing
- Journal of large-scale research facilities JLSRF, Vol. 5
Ultrahigh-order Maxwell solver with extreme scalability for electromagnetic PIC simulations of plasmas
journal, July 2018
- Vincenti, Henri; Vay, Jean-Luc
- Computer Physics Communications, Vol. 228
Numerical Cherenkov instabilities in electromagnetic particle codes
journal, August 1974
- Godfrey, Brendan B.
- Journal of Computational Physics, Vol. 15, Issue 4
Improved numerical Cherenkov instability suppression in the generalized PSTD PIC algorithm
journal, November 2015
- Godfrey, Brendan B.; Vay, Jean-Luc
- Computer Physics Communications, Vol. 196
Generation of a Small-Scale Quasi-Static Magnetic Field and Fast Particles during the Collision of Electron-Positron Plasma Clouds
journal, May 1998
- Kazimura, Y.; Sakai, J. I.; Neubert, T.
- The Astrophysical Journal, Vol. 498, Issue 2
Particle simulation of plasmas
journal, April 1983
- Dawson, John M.
- Reviews of Modern Physics, Vol. 55, Issue 2
Pseudospectral Maxwell solvers for an accurate modeling of Doppler harmonic generation on plasma mirrors with particle-in-cell codes
journal, September 2017
- Blaclard, G.; Vincenti, H.; Lehe, R.
- Physical Review E, Vol. 96, Issue 3