Progress on the Kinetic MHD model

We have developed a Lorentz force ion, fluid electron kinetic MHD hybrid model.\footnote{D. Barnes, \emph{et al}, Phys. Plasmas {\bf{15}}, 055702 (2008)} Here we implement the algorithm in the GEM gyrokinetic turbulence code. Numerically, for some applications, relatively large ion Larmor radius and gyrofrequency motivates the use of Lorentz force ions. An implicit algorithm,\footnote{Yang Chen, Scott E. Parker, to be submitted (2008).} which extends the GEM code to use Lorentz force ions and drift kinetic electrons has already been developed and work is already underway to include gyrokinetic electrons. Based on this ion model but with isothermal fluid electrons, we will use GEM to reproduce our previous results on shearless Alfven waves, ion sound waves and the Whistler waves.\footnote{D. Barnes, \emph{et al}, Phys. Plasmas {\bf{15}}, 055702 (2008)} This model will be used to investigate the linear stability of 1-D Harris model with a guide field. Eventually we will add gyrokinetic electrons. Work is under way to implement the Harris equilibrium in GEM.\footnote{X. Y. Wang, \emph{et al}, Phys. Plasmas {\bf{15}}, 072103 (2008).}$^,$\footnote{E. G. Harris, Nuovo Cimento, {\bf{23}} 115 (1962).}