Turbulent transport of toroidal angular momentum in low flow gyrokinetics

We derive a self-consistent equation for the turbulent transport of toroidal angular momentum in tokamaks in the low flow ordering that only requires solving gyrokinetic Fokker-Planck and quasineutrality equations correct to second order in an expansion on the gyroradius over scale length. We also show that according to our orderings the long wavelength toroidal rotation and the long wavelength radial electric field satisfy the neoclassical relation that gives the toroidal rotation as a function of the radial electric field and the radial gradients of pressure and temperature. Thus, the radial electric field can be solved for once the toroidal rotation is calculated from the transport of toroidal angular momentum. Unfortunately, even though this methodology only requires a gyrokinetic model correct to second order in gyroradius over scale length, current gyrokinetic simulations are only valid to first order. To overcome this difficulty, we exploit the smallish ratio B-p/B, where B is the total magnetic field and B-p is its poloidal component. When B-p/B is small, the usual first order gyrokinetic equation provides solutions that are accurate enough to employ for our expression for the transport of toroidal angular momentum. We show that current delta f and full f simulations only need small corrections to achieve this accuracy. Full f simulations, however, are still unable to determine the long wavelength, radial electric field from the quasineutrality equation.