Enhanced Collisional Losses from a Magnetic Mirror Using the Lenard-Bernstein Collision Operator

Collisions play a crucial role in governing particle and energy transport in plasmas confined in a magnetic mirror trap. Modern gyrokinetic codes are used to model transport in magnetic mirrors, but some of these codes utilize approximate model collision operators. This study focuses on a Pastukhov-style method of images calculation of particle and energy confinement times using a Lenard-Bernstein model collision operator. Prior work on parallel particle and energy balances used a different Fokker-Planck plasma collision operator and the method needs to be extended in non-trivial ways to study the Lenard-Bernstein operator. To assess the effectiveness of our approach, we compare our results with a modern finite element solver. Our findings reveal that the particle confinement time scales like $a \exp(a^2)$ using the Lenard-Bernstein operator, in contrast to the more accurate scaling that the Coulomb collision operator would yield $a^2 \exp(a^2)$, where $a^2$ is approximately proportional to the ambipolar potential. We propose that codes modeling collisional losses in a magnetic mirrors utilizing the Lenard-Bernstein or Dougherty collision operator scale their collision frequency of any electrostatically confined species. This study illuminates the intricate role the collision operator plays in the Pastukhov-style method of images calculation of collisional confinement.
Comment: 21 pages, 4 figures, 3 tables, submitted to journal of plasma physics