Projects per year
Abstract
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 threedimensional particle momentum space and the fivedimensional guiding center phase space. Detailed description is provided for both the physics and implementation of the operator. The focus is in adaptive integration of stochastic differential equations, which is an overlooked aspect among existing Monte Carlo implementations of Coulomb collision operators. We verify that our operator converges to known analytical results and demonstrate that careless implementation of the adaptive time step can lead to severely erroneous results. The operator is provided as a selfcontained Fortran 95 module and can be included into existing orbitfollowing tools that trace either the full Larmor motion or the guiding center dynamics. The adaptive timestepping algorithm is expected to be useful in situations where the collision frequencies vary greatly over the course of a simulation. Examples include the slowingdown of fusion products or other fast ions, and the Dreicer generation of runaway electrons as well as the generation of fast ions or electrons with ion or electron cyclotron resonance heating. Program summary Program Title: AMCC –(A)daptive (M)onte(C)arlo (C)oulomb collisions Program Files doi: http://dx.doi.org/10.17632/wd3t3t6dy7.1 Licensing provisions: MIT licence Programming language: Fortran 95 Nature of problem: Testparticle tracing is a common feat within existing fusion applications. While efficient adaptive methods exist for integrating the incompressible Hamiltonian flow, the effects of Coulomb collisions are commonly implemented with far less sophisticated algorithms. Solution method: The relativistic Fokker–Planck equation for testparticles in Maxwell–Jüttner background plasmas is converted into a stochastic differential equation. The stochastic differential equation is solved using adaptive Monte Carlo techniques. Methods to evaluate the effect of Coulomb collisions for both the threedimensional particle momentum space and the fivedimensional reduced guiding center phase space are included. Additional comments including restrictions and unusual features: The package includes optionality to evaluate the relativistic Fokker–Planck coefficients, a feature useful for constructing accurate orbit averaged collision operators. The package also provides explicit onestep symplectic integrator for the relativistic Lorentz force that can be used for tracing testparticles in given electromagnetic backgrounds.
Original language  English 

Pages (fromto)  374383 
Number of pages  10 
Journal  Computer Physics Communications 
Volume  222 
DOIs  
Publication status  Published  1 Jan 2018 
MoE publication type  A1 Journal articlerefereed 
Keywords
 Coulomb collision
 Fokker–Planck equation
 Milstein method
 Monte Carlo
Fingerprint Dive into the research topics of 'Adaptive timestepping Monte Carlo integration of Coulomb collisions'. Together they form a unique fingerprint.
Datasets

Adaptive timestepping Monte Carlo integration of Coulomb collisions
Särkimäki, K. (Creator), Hirvijoki, E. (Contributor) & Terävä, J. (Contributor), Mendeley Data, 1 Jan 2017
DOI: 10.17632/wd3t3t6dy7.1, https://data.mendeley.com/datasets/wd3t3t6dy7/1
Dataset
Projects
 1 Finished

also: Fusion with a twist: Racing particles (also) in W7X stellarator
Kontula, J., Varje, J., Särkimäki, K., KurkiSuonio, T., Ollus, P. & Sanchis Sanchez, S.
01/09/2016 → 31/08/2020
Project: Academy of Finland: Other research funding