Energy-momentum-conserving stochastic differential equations and algorithms for nonlinear Landau-Fokker-Planck equation
arXiv:2410.12079 · doi:10.1103/PhysRevE.111.025211
Abstract
Coulomb collision is a fundamental diffusion process in plasmas that can be described by the Landau-Fokker-Planck (LFP) equation or the stochastic differential equation (SDE). While energy and momentum are conserved exactly in the LFP equation, they are conserved only on average by the conventional corresponding SDEs, suggesting that the underlying stochastic process may not be well-defined by such SDEs. In this study, we derive new SDEs with exact energy-momentum conservation for the Coulomb collision by factorizing the collective effect of field particles into individual particles and enforcing Newton's third law. These SDEs, when interpreted in the Stratonovich sense, have a particularly simple form that represents pure diffusion between particles without drag. Numerical algorithms that preserve discrete conservation laws are developed and benchmarked in various relaxation processes. Techniques to reduce computational complexity are also discussed.
References in corpus (7)
- Structure and structure-preserving algorithms for plasma physics
- A guiding-center Fokker-Planck collision operator for nonuniform magnetic fields
- Multilevel Monte Carlo simulation of Coulomb collisions
- Monte Carlo implementation of a guiding-center Fokker-Planck kinetic equation
- Structure-preserving marker-particle discretizations of Coulomb collisions for particle-in-cell codes
- Simulating Pitch Angle Scattering Using An Explicitly Solvable Energy Conserving Algorithm
- Multispecies structure-preserving particle discretization of the Landau collision operator