Adaptive time-stepping Monte Carlo integration of Coulomb collisions
From MaRDI portal
Publication:6155750
DOI10.1016/J.CPC.2017.09.025arXiv1701.05043MaRDI QIDQ6155750FDOQ6155750
Authors: K. Särkimäki, E. Hirvijoki, J. Terävä 
Publication date: 7 June 2023
Published in: Computer Physics Communications (Search for Journal in Brave)
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"uttner 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. 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 self-contained Fortran 95 module and can be included into existing orbit-following tools that trace either the full Larmor motion or the guiding center dynamics. The adaptive time-stepping algorithm is expected to be useful in situations where the collision frequencies vary greatly over the course of a simulation. Examples include the slowing-down 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.
Full work available at URL: https://arxiv.org/abs/1701.05043
Cites Work
- Variable Step Size Control in the Numerical Solution of Stochastic Differential Equations
- Monte Carlo evaluation of transport coefficients
- Title not available (Why is that?)
- Higher-order time integration of Coulomb collisions in a plasma using Langevin equations
- On the simulation of iterated Itô integrals.
- An adaptive timestepping algorithm for stochastic differential equations.
- The Langevin equation. With applications to stochastic problems in physics, chemistry and electrical engineering.
- Multilevel Monte Carlo simulation of Coulomb collisions
Cited In (4)
- The Adaptive Monte Carlo Toolbox for Phase Space Integration and Generation
- Adjoint Monte Carlo Calculations of Pulse-Height-Spectrum
- Convergence issues in derivatives of Monte Carlo null-collision integral formulations: a solution
- Collisional \(N\)-body numerical integrator with applications to charged particle dynamics
This page was built for publication: Adaptive time-stepping Monte Carlo integration of Coulomb collisions
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6155750)
