A gyrokinetic ion/mass-less fluid electron hybrid model as implemented in the GEM code [Y. Chen and S. E. Parker, J. Comput. Phys. 220, 837 (2007)] is used to study the reverse shear Alfvén eigenmodes (RSAE) observed in DIII-D, discharge #142111. This is a well diagnosed case with measurement of the core-localized RSAE mode structures and the mode frequency, which can be used to compare with simulations. Simulations reproduce many features of the observation, including the mode frequency up-sweeping in time and the sweeping range. A new algorithmic feature is added to the GEM code for this study. Instead of the gyrokinetic Poisson equation itself, its time derivative, or the vorticity equation, is solved to obtain the electric potential. This permits a numerical scheme that ensures the E × B convection of the equilibrium density profiles of each species cancel each other in the absence of any finite-Larmor-radius effects. These nonlinear simulations generally result in an electron temperature fluctuation level that is comparable to measurements, and a mode frequency spectrum broader than the experimental spectrum. The spectral width from simulations can be reduced if less steep beam density profiles are used, but then the experimental fluctuation level can be reproduced only if a collision rate above the classical level is assumed. © 2013 American Institute of Physics.