skip to main content

SciTech ConnectSciTech Connect

Title: Numerical stability of relativistic beam multidimensional PIC simulations employing the Esirkepov algorithm

Rapidly growing numerical instabilities routinely occur in multidimensional particle-in-cell computer simulations of plasma-based particle accelerators, astrophysical phenomena, and relativistic charged particle beams. Reducing instability growth to acceptable levels has necessitated higher resolution grids, high-order field solvers, current filtering, etc. except for certain ratios of the time step to the axial cell size, for which numerical growth rates and saturation levels are reduced substantially. This paper derives and solves the cold beam dispersion relation for numerical instabilities in multidimensional, relativistic, electromagnetic particle-in-cell programs employing either the standard or the Cole–Karkkainnen finite difference field solver on a staggered mesh and the common Esirkepov current-gathering algorithm. Good overall agreement is achieved with previously reported results of the WARP code. In particular, the existence of select time steps for which instabilities are minimized is explained. Additionally, an alternative field interpolation algorithm is proposed for which instabilities are almost completely eliminated for a particular time step in ultra-relativistic simulations.
Authors:
 [1] ;  [2]
  1. University of Maryland, College Park, MD 20742 (United States)
  2. Lawrence Berkeley National Laboratory, Berkeley, CA 94720 (United States)
Publication Date:
OSTI Identifier:
22230782
Resource Type:
Journal Article
Resource Relation:
Journal Name: Journal of Computational Physics; Journal Volume: 248; Other Information: Copyright (c) 2013 Elsevier Science B.V., Amsterdam, The Netherlands, All rights reserved.; Country of input: International Atomic Energy Agency (IAEA)
Country of Publication:
United States
Language:
English
Subject:
43 PARTICLE ACCELERATORS; ACCELERATORS; ALGORITHMS; ASTROPHYSICS; CHARGED PARTICLES; COMPUTERIZED SIMULATION; DISPERSION RELATIONS; INSTABILITY; INTERPOLATION; PARTICLE BEAMS; RELATIVISTIC RANGE