Projekteja vuodessa
Abstrakti
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.
Alkuperäiskieli  Englanti 

Sivut  374383 
Sivumäärä  10 
Julkaisu  Computer Physics Communications 
Vuosikerta  222 
DOI  pysyväislinkit  
Tila  Julkaistu  1 tammik. 2018 
OKMjulkaisutyyppi  A1 Alkuperäisartikkeli tieteellisessä aikakauslehdessä 
Sormenjälki
Sukella tutkimusaiheisiin 'Adaptive timestepping Monte Carlo integration of Coulomb collisions'. Ne muodostavat yhdessä ainutlaatuisen sormenjäljen.Tietoaineistot

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

also: Fuusiota vauhtikierteellä: nopeita hiukkasia jäljittämässä jopa W7X stellaraattorissa
Särkimäki, K., KurkiSuonio, T., Kontula, J., Ollus, P., Varje, J., Sanchis, L. & Ikäheimo, A.
01/09/2016 → 31/08/2020
Projekti: Academy of Finland: Other research funding