A circular equilibrium model for local gyrokinetic simulations of ion temperature gradient fluctuations in reversed field pinches

A simple large-aspect-ratio (R0/r) circular equilibrium model is developed for low-beta reversed field pinch (RFP) geometry. The model is suitable for treating small scale instability and turbulent transport driven by ion temperature gradient (ITG) and related electron drift modes in gyrokinetic simulations. The equilibrium model is an RFP generalization of the common tokamak s-α model to small safety factor (q), where the poloidal field dominates the toroidal field. The model accommodates the RFP toroidal field reversal (where q vanishes) by generalizing the cylindrical force-free Bessel function model (BFM) [J. B. Taylor, Phys. Rev. Lett. 33, 1139 (1974)] to toroidal geometry. The global equilibrium can be described in terms of the RFP field reversal and pinch parameters [F,Θ]. This new toroidal Bessel function model (TBFM) has been incorporated into the gyrokinetic code GYRO [J. Candy and R. E. Waltz, J.Comput. Phys. 186, 545 (2003)] and used here to explore local electrostatic ITG adiabatic electron instability rates for typical low-q RFP parameters.