skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: On the Numerical Dispersion of Electromagnetic Particle-In-Cell Code : Finite Grid Instability

Abstract

The Particle-In-Cell (PIC) method is widely used in relativistic particle beam and laser plasma modeling. However, the PIC method exhibits numerical instabilities that can render unphysical simulation results or even destroy the simulation. For electromagnetic relativistic beam and plasma modeling, the most relevant numerical instabilities are the finite grid instability and the numerical Cherenkov instability. We review the numerical dispersion relation of the electromagnetic PIC algorithm to analyze the origin of these instabilities. We rigorously derive the faithful 3D numerical dispersion of the PIC algorithm, and then specialize to the Yee FDTD scheme. In particular, we account for the manner in which the PIC algorithm updates and samples the fields and distribution function. Temporal and spatial phase factors from solving Maxwell's equations on the Yee grid with the leapfrog scheme are also explicitly accounted for. Numerical solutions to the electrostatic-like modes in the 1D dispersion relation for a cold drifting plasma are obtained for parameters of interest. In the succeeding analysis, we investigate how the finite grid instability arises from the interaction of the numerical 1D modes admitted in the system and their aliases. The most significant interaction is due critically to the correct representation of the operators in themore » dispersion relation. We obtain a simple analytic expression for the peak growth rate due to this interaction.« less

Authors:
 [1];  [2];  [1];  [1];  [1];  [1]
  1. Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
  2. (United States) Dept. of Physics and Astronomy
Publication Date:
Research Org.:
Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
Sponsoring Org.:
USDOE
OSTI Identifier:
1143957
Report Number(s):
LA-UR-14-25341
DOE Contract Number:  
AC52-06NA25396
Resource Type:
Technical Report
Country of Publication:
United States
Language:
English
Subject:
70 PLASMA PHYSICS AND FUSION TECHNOLOGY

Citation Formats

Meyers, Michael David, Univ. of California, Los Angeles, CA, Huang, Chengkun, Zeng, Yong, Yi, Sunghwan, and Albright, Brian James. On the Numerical Dispersion of Electromagnetic Particle-In-Cell Code : Finite Grid Instability. United States: N. p., 2014. Web. doi:10.2172/1143957.
Meyers, Michael David, Univ. of California, Los Angeles, CA, Huang, Chengkun, Zeng, Yong, Yi, Sunghwan, & Albright, Brian James. On the Numerical Dispersion of Electromagnetic Particle-In-Cell Code : Finite Grid Instability. United States. doi:10.2172/1143957.
Meyers, Michael David, Univ. of California, Los Angeles, CA, Huang, Chengkun, Zeng, Yong, Yi, Sunghwan, and Albright, Brian James. Tue . "On the Numerical Dispersion of Electromagnetic Particle-In-Cell Code : Finite Grid Instability". United States. doi:10.2172/1143957. https://www.osti.gov/servlets/purl/1143957.
@article{osti_1143957,
title = {On the Numerical Dispersion of Electromagnetic Particle-In-Cell Code : Finite Grid Instability},
author = {Meyers, Michael David and Univ. of California, Los Angeles, CA and Huang, Chengkun and Zeng, Yong and Yi, Sunghwan and Albright, Brian James},
abstractNote = {The Particle-In-Cell (PIC) method is widely used in relativistic particle beam and laser plasma modeling. However, the PIC method exhibits numerical instabilities that can render unphysical simulation results or even destroy the simulation. For electromagnetic relativistic beam and plasma modeling, the most relevant numerical instabilities are the finite grid instability and the numerical Cherenkov instability. We review the numerical dispersion relation of the electromagnetic PIC algorithm to analyze the origin of these instabilities. We rigorously derive the faithful 3D numerical dispersion of the PIC algorithm, and then specialize to the Yee FDTD scheme. In particular, we account for the manner in which the PIC algorithm updates and samples the fields and distribution function. Temporal and spatial phase factors from solving Maxwell's equations on the Yee grid with the leapfrog scheme are also explicitly accounted for. Numerical solutions to the electrostatic-like modes in the 1D dispersion relation for a cold drifting plasma are obtained for parameters of interest. In the succeeding analysis, we investigate how the finite grid instability arises from the interaction of the numerical 1D modes admitted in the system and their aliases. The most significant interaction is due critically to the correct representation of the operators in the dispersion relation. We obtain a simple analytic expression for the peak growth rate due to this interaction.},
doi = {10.2172/1143957},
journal = {},
number = ,
volume = ,
place = {United States},
year = {Tue Jul 15 00:00:00 EDT 2014},
month = {Tue Jul 15 00:00:00 EDT 2014}
}

Technical Report:

Save / Share: