skip to main content

DOE PAGESDOE PAGES

Title: Adaptive time-stepping Monte Carlo integration of Coulomb collisions

Here, we report an accessible and robust tool for evaluating the effects of Coulomb collisions on a test particle in a plasma that obeys Maxwell–Jüttner statistics. The implementation is based on the Beliaev–Budker collision integral which allows both the test particle and the background plasma to be relativistic. The integration method supports adaptive time stepping, which is shown to greatly improve the computational efficiency. The Monte Carlo method is implemented for both the three-dimensional particle momentum space and the five-dimensional guiding center phase space.
Authors:
ORCiD logo [1] ;  [2] ;  [1]
  1. Aalto Univ., Espoo (Finland)
  2. Princeton Plasma Physics Lab. (PPPL), Princeton, NJ (United States)
Publication Date:
Grant/Contract Number:
259675; 298126; AC02–09CH11466
Type:
Accepted Manuscript
Journal Name:
Computer Physics Communications
Additional Journal Information:
Journal Volume: 222; Journal Issue: C; Journal ID: ISSN 0010-4655
Publisher:
Elsevier
Research Org:
Princeton Plasma Physics Lab. (PPPL), Princeton, NJ (United States)
Sponsoring Org:
USDOE
Country of Publication:
United States
Language:
English
Subject:
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS; 97 MATHEMATICS AND COMPUTING; Coulomb collision; Monte Carlo; Fokker–Planck equation; Milstein method
OSTI Identifier:
1441004

Sarkimaki, Konsta, Hirvijoki, E., and Terava, J.. Adaptive time-stepping Monte Carlo integration of Coulomb collisions. United States: N. p., Web. doi:10.1016/j.cpc.2017.09.025.
Sarkimaki, Konsta, Hirvijoki, E., & Terava, J.. Adaptive time-stepping Monte Carlo integration of Coulomb collisions. United States. doi:10.1016/j.cpc.2017.09.025.
Sarkimaki, Konsta, Hirvijoki, E., and Terava, J.. 2017. "Adaptive time-stepping Monte Carlo integration of Coulomb collisions". United States. doi:10.1016/j.cpc.2017.09.025. https://www.osti.gov/servlets/purl/1441004.
@article{osti_1441004,
title = {Adaptive time-stepping Monte Carlo integration of Coulomb collisions},
author = {Sarkimaki, Konsta and Hirvijoki, E. and Terava, J.},
abstractNote = {Here, we report an accessible and robust tool for evaluating the effects of Coulomb collisions on a test particle in a plasma that obeys Maxwell–Jüttner statistics. The implementation is based on the Beliaev–Budker collision integral which allows both the test particle and the background plasma to be relativistic. The integration method supports adaptive time stepping, which is shown to greatly improve the computational efficiency. The Monte Carlo method is implemented for both the three-dimensional particle momentum space and the five-dimensional guiding center phase space.},
doi = {10.1016/j.cpc.2017.09.025},
journal = {Computer Physics Communications},
number = C,
volume = 222,
place = {United States},
year = {2017},
month = {10}
}