Gyrokinetic particle simulation of fusion plasmas: path to petascale computing

Gyrokinetic particle simulation of fusion plasmas for studying turbulent transport on state-of-theart computers has a long history of important scientific discoveries. The primary examples are: (i) the identification of ion temperature gradient (ITG) drift turbulence as the most plausible process responsible for the thermal transport observed in tokamak experiments; (ii) the reduction of such transport due to the presence of zonal flows; (iii) the confinement scaling trends associated with size of the plasma and also with the ionic isotope species. With the availability of terascale computers in recent years, we have also been able to carry out simulations with improved physics fidelity using experimentally relevant parameters. Computationally, we have demonstrated that our lead Particle-in- Cell (PIC) code, the Gyrokinetic Turbulence Code (GTC), is portable, efficient, and scalable on various MPP platforms. Convergence studies with unprecedented phase-space resolution have also been carried out. Since petascale resources are expected to be available in the near future, we have also engaged in developing better physics models and more efficient numerical algorithms to take advantage of this exciting opportunity. For the near term, we are interested in understanding some basic physics issues related to burning plasmas experiments in International Thermonuclear Experimental Reactor (ITER) - a multi-billion dollar device to be constructed over the next decade. Our long range goal is to carry out integrated simulations for ITER plasmas for a wide range of temporal and spatial scales, including high-frequency short-wavelength wave heating, low-frequency meso-scale transport, and low-frequency large scale magnetohydrodynamic (MHD) physics on these computers.