Time-step considerations in particle simulation algorithms for coulomb collisions in plasmas

Bruce I. Cohen, Andris M. Dimits, Alex Friedman, Russel Caflisch

Research output: Contribution to journalArticle

Abstract

The accuracy of first-order Euler and higher-order time-integration algorithms for grid-based Langevin equations collision models in a specific relaxation test problem is assessed. We show that statistical noise errors can overshadow time-step errors and argue that statistical noise errors can be conflated with time-step effects. Using a higher-order integration scheme may not achieve any benefit in accuracy for examples of practical interest. We also investigate the collisional relaxation of an initial electron-ion relative drift and the collisional relaxation to a resistive steady-state in which a quasi-steady current is driven by a constant applied electric field, as functions of the time step used to resolve the collision processes using binary and grid-based, test-particle Langevin equations models. We compare results from two grid-based Langevin equations collision algorithms to results from a binary collision algorithm for modeling electron-ion collisions. Some guidance is provided on how large a time step can be used compared to the inverse of the characteristic collision frequency for specific relaxation processes.

Original languageEnglish (US)
Article number5475281
Pages (from-to)2394-2406
Number of pages13
JournalIEEE Transactions on Plasma Science
Volume38
Issue number9 PART 1
DOIs
StatePublished - Sep 2010

Fingerprint

Coulomb collisions
collisions
simulation
grids
ions
electrons
electric fields

Keywords

  • Algorithms
  • collision processes
  • computer applications
  • numerical analysis
  • particle collisions
  • plasmas

ASJC Scopus subject areas

  • Condensed Matter Physics
  • Nuclear and High Energy Physics

Cite this

Time-step considerations in particle simulation algorithms for coulomb collisions in plasmas. / Cohen, Bruce I.; Dimits, Andris M.; Friedman, Alex; Caflisch, Russel.

In: IEEE Transactions on Plasma Science, Vol. 38, No. 9 PART 1, 5475281, 09.2010, p. 2394-2406.

Research output: Contribution to journalArticle

Cohen, Bruce I. ; Dimits, Andris M. ; Friedman, Alex ; Caflisch, Russel. / Time-step considerations in particle simulation algorithms for coulomb collisions in plasmas. In: IEEE Transactions on Plasma Science. 2010 ; Vol. 38, No. 9 PART 1. pp. 2394-2406.
@article{55666b2da02246fa8ba0092c8e805052,
title = "Time-step considerations in particle simulation algorithms for coulomb collisions in plasmas",
abstract = "The accuracy of first-order Euler and higher-order time-integration algorithms for grid-based Langevin equations collision models in a specific relaxation test problem is assessed. We show that statistical noise errors can overshadow time-step errors and argue that statistical noise errors can be conflated with time-step effects. Using a higher-order integration scheme may not achieve any benefit in accuracy for examples of practical interest. We also investigate the collisional relaxation of an initial electron-ion relative drift and the collisional relaxation to a resistive steady-state in which a quasi-steady current is driven by a constant applied electric field, as functions of the time step used to resolve the collision processes using binary and grid-based, test-particle Langevin equations models. We compare results from two grid-based Langevin equations collision algorithms to results from a binary collision algorithm for modeling electron-ion collisions. Some guidance is provided on how large a time step can be used compared to the inverse of the characteristic collision frequency for specific relaxation processes.",
keywords = "Algorithms, collision processes, computer applications, numerical analysis, particle collisions, plasmas",
author = "Cohen, {Bruce I.} and Dimits, {Andris M.} and Alex Friedman and Russel Caflisch",
year = "2010",
month = "9",
doi = "10.1109/TPS.2010.2049589",
language = "English (US)",
volume = "38",
pages = "2394--2406",
journal = "IEEE Transactions on Plasma Science",
issn = "0093-3813",
publisher = "Institute of Electrical and Electronics Engineers Inc.",
number = "9 PART 1",

}

TY - JOUR

T1 - Time-step considerations in particle simulation algorithms for coulomb collisions in plasmas

AU - Cohen, Bruce I.

AU - Dimits, Andris M.

AU - Friedman, Alex

AU - Caflisch, Russel

PY - 2010/9

Y1 - 2010/9

N2 - The accuracy of first-order Euler and higher-order time-integration algorithms for grid-based Langevin equations collision models in a specific relaxation test problem is assessed. We show that statistical noise errors can overshadow time-step errors and argue that statistical noise errors can be conflated with time-step effects. Using a higher-order integration scheme may not achieve any benefit in accuracy for examples of practical interest. We also investigate the collisional relaxation of an initial electron-ion relative drift and the collisional relaxation to a resistive steady-state in which a quasi-steady current is driven by a constant applied electric field, as functions of the time step used to resolve the collision processes using binary and grid-based, test-particle Langevin equations models. We compare results from two grid-based Langevin equations collision algorithms to results from a binary collision algorithm for modeling electron-ion collisions. Some guidance is provided on how large a time step can be used compared to the inverse of the characteristic collision frequency for specific relaxation processes.

AB - The accuracy of first-order Euler and higher-order time-integration algorithms for grid-based Langevin equations collision models in a specific relaxation test problem is assessed. We show that statistical noise errors can overshadow time-step errors and argue that statistical noise errors can be conflated with time-step effects. Using a higher-order integration scheme may not achieve any benefit in accuracy for examples of practical interest. We also investigate the collisional relaxation of an initial electron-ion relative drift and the collisional relaxation to a resistive steady-state in which a quasi-steady current is driven by a constant applied electric field, as functions of the time step used to resolve the collision processes using binary and grid-based, test-particle Langevin equations models. We compare results from two grid-based Langevin equations collision algorithms to results from a binary collision algorithm for modeling electron-ion collisions. Some guidance is provided on how large a time step can be used compared to the inverse of the characteristic collision frequency for specific relaxation processes.

KW - Algorithms

KW - collision processes

KW - computer applications

KW - numerical analysis

KW - particle collisions

KW - plasmas

UR - http://www.scopus.com/inward/record.url?scp=77956617263&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=77956617263&partnerID=8YFLogxK

U2 - 10.1109/TPS.2010.2049589

DO - 10.1109/TPS.2010.2049589

M3 - Article

AN - SCOPUS:77956617263

VL - 38

SP - 2394

EP - 2406

JO - IEEE Transactions on Plasma Science

JF - IEEE Transactions on Plasma Science

SN - 0093-3813

IS - 9 PART 1

M1 - 5475281

ER -