A three-dimensional, nonlinear gyrokinetic simulation of plasma turbulence resolving scales from the ion to electron gyroradius with a realistic mass ratio is presented, where all damping is provided by resolved physical mechanisms. The resulting energy spectra are quantitatively consistent with a magnetic power spectrum scaling of k(-2.8) as observed in in situ spacecraft measurements of the "dissipation range" of solar wind turbulence. Despite the strongly nonlinear nature of the turbulence, the linear kinetic Alfvén wave mode quantitatively describes the polarization of the turbulent fluctuations. The collisional ion heating is measured at subion-Larmor radius scales, which provides evidence of the ion entropy cascade in an electromagnetic turbulence simulation.

The Hasegawa-Wakatani equations, coupling plasma density and electrostatic potential through an approximation to the physics of parallel electron motions, are a simple model that describes resistive drift wave turbulence. We present numerical analyses of bifurcation phenomena in the model that provide new insights into the interactions between turbulence and zonal flows in the tokamak plasma edge region. The simulation results show a regime where, after an initial transient, drift wave turbulence is suppressed through zonal flow generation. As a parameter controlling the strength of the turbulence is tuned, this zonal flow dominated state is rapidly destroyed and a turbulence-dominated state re-emerges. The transition is explained in terms of the KelvinHelmholtz stability of zonal flows. This is the first observation of an upshift of turbulence onset in the resistive drift wave system, which is analogous to the well-known Dimits shift in turbulence driven by ion temperature gradients. * Electronic address: ryusuke.numata@anu.edu.au

a b s t r a c tThe gyrokinetic simulation code AstroGK is developed to study fundamental aspects of kinetic plasmas and for applications mainly to astrophysical problems. AstroGK is an Eulerian slab code that solves the electromagnetic gyrokinetic-Maxwell equations in five-dimensional phase space, and is derived from the existing gyrokinetics code GS2 by removing magnetic geometry effects. Algorithms used in the code are described. The code is benchmarked using linear and nonlinear problems. Serial and parallel performance scalings are also presented.