Importance of gyrokinetic exact Fokker-Planck collisions in fusion plasma turbulence.

Gyrokinetic simulations of turbulence are fundamental to understanding and predicting particle and energy loss in magnetic fusion devices. Previous works have used model collision operators with approximate field-particle terms of unknown accuracy and/or have neglected collisional finite Larmor radius effects. This work moves beyond models to demonstrate important corrections using a gyrokinetic Fokker-Planck collision operator with the exact field-particle terms, in realistic simulations of turbulence in magnetically confined fusion plasmas. The exact operator shows significant corrections for temperature-gradient-driven trapped electron mode turbulence and zonal flow damping, and for microtearing modes in a Joint European Torus pedestal under ITER-like wall conditions. Analysis of the corrections using parameter scans motivates an accurate model which closely reproduces the exact results while reducing computational demands.

Gyrokinetic Landau collision operator in conservative form.

A gyrokinetic linearized exact (not model) Landau collision operator is derived by transforming the symmetric and conservative Landau form. The formulation obtains the velocity-space flux density and preserves the operator's conservative form as the divergence of this flux density. The operator contains both test-particle and field-particle contributions, and finite Larmor radius effects are evaluated in either Bessel function series or gyrophase integrals. While equivalent to the gyrokinetic Fokker-Planck form with Rosenbluth potentials [B. Li and D. R. Ernst, Phys. Rev. Lett. 106, 195002 (2011)10.1103/PhysRevLett.106.195002], the gyrokinetic conservative Landau form explicitly preserves the symmetry between test-particle and field-particle contributions, which underlies the conservation laws and the H theorem, and enables discretization with a finite-volume or spectral method to preserve the conservation properties numerically, independent of resolution. The form of the exact linearized field-particle terms differs from those of widely used model operators. We show the finite Larmor radius corrections to the field-particle terms in the exact linearized operator involve Bessel functions of all orders, while present model field-particle terms involve only the first two Bessel functions. This new symmetric and conservative formulation enables the gyrokinetic exact linearized Landau operator to be implemented in gyrokinetic turbulence codes for comparison with present model operators using similar numerical methods.