Simulating pitch angle scattering using an explicitly solvable energy-conserving algorithm
Particle distribution functions evolving under the Lorentz operator can be simulated with the Langevin equation for pitch angle scattering. This approach is frequently used in particle based Monte-Carlo simulations of plasma collisions, among others. However, most numerical treatments do not guarantee energy conservation, which may lead to unphysical artifacts such as numerical heating and spectra distortions. In this paper, we present a novel structure-preserving numerical algorithm for the Langevin equation for pitch angle scattering. Similar to the well-known Boris algorithm, the proposed numerical scheme takes advantage of the structure-preserving properties of the Cayley transform when calculating the velocity-space rotations. The resulting algorithm is explicitly solvable, while preserving the norm of velocities down to machine precision. We demonstrate that the method has the same order of numerical convergence as the traditional stochastic Euler-Maruyama method. The numerical scheme is benchmarked by simulating the pitch angle scattering of a particle beam, and comparing with the analytical solution. Benchmark results show excellent agreement with theoretical predictions, showcasing the remarkable long time accuracy of the proposed algorithm.
- Research Organization:
- Princeton Plasma Physics Laboratory (PPPL), Princeton, NJ (United States)
- Sponsoring Organization:
- USDOE
- Grant/Contract Number:
- AC02-09CH11466
- OSTI ID:
- 1660496
- Journal Information:
- Physical Review. E, Vol. 102, Issue 3; ISSN 2470-0045
- Publisher:
- American Physical Society (APS)Copyright Statement
- Country of Publication:
- United States
- Language:
- English
Similar Records
An explicitly solvable energy-conserving algorithm for pitch-angle scattering in magnetized plasmas
Higher-order time integration of Coulomb collisions in a plasma using Langevin equations