A hybrid method for accelerated simulation of coulomb collisions in a plasma

Russel Caflisch, Chiaming Wang, Giacomo Dimarco, Bruce Cohen, Andris Dimits

Research output: Contribution to journalArticle

Abstract

If the collisional time scale for Coulomb collisions is comparable to the characteristic time scales for a plasma, then simulation of Coulombcollisions may be important for computation of kinetic plasma dynamics. This can be a computational bottleneck because of the large number of simulated particles and collisions (or phase-space resolution requirements i continuum algorithms), as well as the wide range of collision rates over the velocity distribution function. This paper considers Monte Carlo simulation of Coulomb collisions using the binary collision models of Takizuka and Abe and of Nanbu. It presents a hybrid method for accelerating the computation of Coulomb collisions. The hybrid method represents the velocity distribution function as a combination of a thermal component (a Maxwellian distribution) and a kinetic component (a set of discrete particles). Collisions between particles from the thermal component preserve the Maxwellian; collisions between particles from the kinetic component are performed using the method of Takizuka and Abe or of Nanbu. Collisions between the kinetic and thermal components are performed by sampling a particle from the thermal component and selecting a particle from the kinetic component. Particles are also transferred between the two components according to thermalization and dethermalization probabilities, which are functions of phase space.

Original languageEnglish (US)
Pages (from-to)865-887
Number of pages23
JournalMultiscale Modeling and Simulation
Volume7
Issue number2
DOIs
StatePublished - 2008

Fingerprint

Coulomb collisions
Hybrid Method
Plasma
Collision
collision
Plasmas
plasma
Kinetics
simulation
Simulation
Velocity distribution
collisions
Distribution functions
kinetics
Plasma simulation
Velocity Distribution
velocity distribution
distribution functions
Phase Space
Distribution Function

Keywords

  • Bump on tail
  • Coulomb collisions
  • Hybrid method
  • Plasma
  • Simulation
  • Thermalization

ASJC Scopus subject areas

  • Chemistry(all)
  • Modeling and Simulation
  • Ecological Modeling
  • Physics and Astronomy(all)
  • Computer Science Applications

Cite this

A hybrid method for accelerated simulation of coulomb collisions in a plasma. / Caflisch, Russel; Wang, Chiaming; Dimarco, Giacomo; Cohen, Bruce; Dimits, Andris.

In: Multiscale Modeling and Simulation, Vol. 7, No. 2, 2008, p. 865-887.

Research output: Contribution to journalArticle

Caflisch, Russel ; Wang, Chiaming ; Dimarco, Giacomo ; Cohen, Bruce ; Dimits, Andris. / A hybrid method for accelerated simulation of coulomb collisions in a plasma. In: Multiscale Modeling and Simulation. 2008 ; Vol. 7, No. 2. pp. 865-887.
@article{d13b869e83cd4f0594c3158a6b33c40f,
title = "A hybrid method for accelerated simulation of coulomb collisions in a plasma",
abstract = "If the collisional time scale for Coulomb collisions is comparable to the characteristic time scales for a plasma, then simulation of Coulombcollisions may be important for computation of kinetic plasma dynamics. This can be a computational bottleneck because of the large number of simulated particles and collisions (or phase-space resolution requirements i continuum algorithms), as well as the wide range of collision rates over the velocity distribution function. This paper considers Monte Carlo simulation of Coulomb collisions using the binary collision models of Takizuka and Abe and of Nanbu. It presents a hybrid method for accelerating the computation of Coulomb collisions. The hybrid method represents the velocity distribution function as a combination of a thermal component (a Maxwellian distribution) and a kinetic component (a set of discrete particles). Collisions between particles from the thermal component preserve the Maxwellian; collisions between particles from the kinetic component are performed using the method of Takizuka and Abe or of Nanbu. Collisions between the kinetic and thermal components are performed by sampling a particle from the thermal component and selecting a particle from the kinetic component. Particles are also transferred between the two components according to thermalization and dethermalization probabilities, which are functions of phase space.",
keywords = "Bump on tail, Coulomb collisions, Hybrid method, Plasma, Simulation, Thermalization",
author = "Russel Caflisch and Chiaming Wang and Giacomo Dimarco and Bruce Cohen and Andris Dimits",
year = "2008",
doi = "10.1137/070704939",
language = "English (US)",
volume = "7",
pages = "865--887",
journal = "Multiscale Modeling and Simulation",
issn = "1540-3459",
publisher = "Society for Industrial and Applied Mathematics Publications",
number = "2",

}

TY - JOUR

T1 - A hybrid method for accelerated simulation of coulomb collisions in a plasma

AU - Caflisch, Russel

AU - Wang, Chiaming

AU - Dimarco, Giacomo

AU - Cohen, Bruce

AU - Dimits, Andris

PY - 2008

Y1 - 2008

N2 - If the collisional time scale for Coulomb collisions is comparable to the characteristic time scales for a plasma, then simulation of Coulombcollisions may be important for computation of kinetic plasma dynamics. This can be a computational bottleneck because of the large number of simulated particles and collisions (or phase-space resolution requirements i continuum algorithms), as well as the wide range of collision rates over the velocity distribution function. This paper considers Monte Carlo simulation of Coulomb collisions using the binary collision models of Takizuka and Abe and of Nanbu. It presents a hybrid method for accelerating the computation of Coulomb collisions. The hybrid method represents the velocity distribution function as a combination of a thermal component (a Maxwellian distribution) and a kinetic component (a set of discrete particles). Collisions between particles from the thermal component preserve the Maxwellian; collisions between particles from the kinetic component are performed using the method of Takizuka and Abe or of Nanbu. Collisions between the kinetic and thermal components are performed by sampling a particle from the thermal component and selecting a particle from the kinetic component. Particles are also transferred between the two components according to thermalization and dethermalization probabilities, which are functions of phase space.

AB - If the collisional time scale for Coulomb collisions is comparable to the characteristic time scales for a plasma, then simulation of Coulombcollisions may be important for computation of kinetic plasma dynamics. This can be a computational bottleneck because of the large number of simulated particles and collisions (or phase-space resolution requirements i continuum algorithms), as well as the wide range of collision rates over the velocity distribution function. This paper considers Monte Carlo simulation of Coulomb collisions using the binary collision models of Takizuka and Abe and of Nanbu. It presents a hybrid method for accelerating the computation of Coulomb collisions. The hybrid method represents the velocity distribution function as a combination of a thermal component (a Maxwellian distribution) and a kinetic component (a set of discrete particles). Collisions between particles from the thermal component preserve the Maxwellian; collisions between particles from the kinetic component are performed using the method of Takizuka and Abe or of Nanbu. Collisions between the kinetic and thermal components are performed by sampling a particle from the thermal component and selecting a particle from the kinetic component. Particles are also transferred between the two components according to thermalization and dethermalization probabilities, which are functions of phase space.

KW - Bump on tail

KW - Coulomb collisions

KW - Hybrid method

KW - Plasma

KW - Simulation

KW - Thermalization

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

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

U2 - 10.1137/070704939

DO - 10.1137/070704939

M3 - Article

VL - 7

SP - 865

EP - 887

JO - Multiscale Modeling and Simulation

JF - Multiscale Modeling and Simulation

SN - 1540-3459

IS - 2

ER -