Gyrokinetic particle simulation of microturbulence for general magnetic geometry and experimental profiles
Abstract
Developments in gyrokinetic particle simulation enable the gyrokinetic toroidal code (GTC) to simulate turbulent transport in tokamaks with realistic equilibrium profiles and plasma geometry, which is a critical step in the code–experiment validation process. These new developments include numerical equilibrium representation using Bsplines, a new Poisson solver based on finite difference using fieldaligned mesh and magnetic flux coordinates, a new zonal flow solver for general geometry, and improvements on the conventional fourpoint gyroaverage with nonuniform background marker loading. The gyrokinetic Poisson equation is solved in the perpendicular plane instead of the poloidal plane. Exploiting these new features, GTC is able to simulate a typical DIIID discharge with experimental magnetic geometry and profiles. The simulated turbulent heat diffusivity and its radial profile show good agreement with other gyrokinetic codes. The newly developed nonuniform loading method provides a modified radial transport profile to that of the conventional uniform loading method.
 Authors:
 Institute for Fusion Theory and Simulation, Zhejiang University, Hangzhou 310027 (China)
 (United States)
 Department of Physics and Astronomy, University of California, Irvine, California 92697 (United States)
 (China)
 Publication Date:
 OSTI Identifier:
 22408111
 Resource Type:
 Journal Article
 Resource Relation:
 Journal Name: Physics of Plasmas; Journal Volume: 22; Journal Issue: 2; Other Information: (c) 2015 AIP Publishing LLC; Country of input: International Atomic Energy Agency (IAEA)
 Country of Publication:
 United States
 Language:
 English
 Subject:
 70 PLASMA PHYSICS AND FUSION TECHNOLOGY; COMPUTERIZED SIMULATION; DOUBLET3 DEVICE; EXPERIMENTAL DATA; MAGNETIC FLUX COORDINATES; PARTICLES; PLASMA; POISSON EQUATION; TOROIDAL CONFIGURATION
Citation Formats
Xiao, Yong, Email: yxiao@zju.edu.cn, Department of Physics and Astronomy, University of California, Irvine, California 92697, Holod, Ihor, Wang, Zhixuan, Lin, Zhihong, Fusion Simulation Center, Peking University, Beijing 100871, and Zhang, Taige. Gyrokinetic particle simulation of microturbulence for general magnetic geometry and experimental profiles. United States: N. p., 2015.
Web. doi:10.1063/1.4908275.
Xiao, Yong, Email: yxiao@zju.edu.cn, Department of Physics and Astronomy, University of California, Irvine, California 92697, Holod, Ihor, Wang, Zhixuan, Lin, Zhihong, Fusion Simulation Center, Peking University, Beijing 100871, & Zhang, Taige. Gyrokinetic particle simulation of microturbulence for general magnetic geometry and experimental profiles. United States. doi:10.1063/1.4908275.
Xiao, Yong, Email: yxiao@zju.edu.cn, Department of Physics and Astronomy, University of California, Irvine, California 92697, Holod, Ihor, Wang, Zhixuan, Lin, Zhihong, Fusion Simulation Center, Peking University, Beijing 100871, and Zhang, Taige. Sun .
"Gyrokinetic particle simulation of microturbulence for general magnetic geometry and experimental profiles". United States.
doi:10.1063/1.4908275.
@article{osti_22408111,
title = {Gyrokinetic particle simulation of microturbulence for general magnetic geometry and experimental profiles},
author = {Xiao, Yong, Email: yxiao@zju.edu.cn and Department of Physics and Astronomy, University of California, Irvine, California 92697 and Holod, Ihor and Wang, Zhixuan and Lin, Zhihong and Fusion Simulation Center, Peking University, Beijing 100871 and Zhang, Taige},
abstractNote = {Developments in gyrokinetic particle simulation enable the gyrokinetic toroidal code (GTC) to simulate turbulent transport in tokamaks with realistic equilibrium profiles and plasma geometry, which is a critical step in the code–experiment validation process. These new developments include numerical equilibrium representation using Bsplines, a new Poisson solver based on finite difference using fieldaligned mesh and magnetic flux coordinates, a new zonal flow solver for general geometry, and improvements on the conventional fourpoint gyroaverage with nonuniform background marker loading. The gyrokinetic Poisson equation is solved in the perpendicular plane instead of the poloidal plane. Exploiting these new features, GTC is able to simulate a typical DIIID discharge with experimental magnetic geometry and profiles. The simulated turbulent heat diffusivity and its radial profile show good agreement with other gyrokinetic codes. The newly developed nonuniform loading method provides a modified radial transport profile to that of the conventional uniform loading method.},
doi = {10.1063/1.4908275},
journal = {Physics of Plasmas},
number = 2,
volume = 22,
place = {United States},
year = {Sun Feb 15 00:00:00 EST 2015},
month = {Sun Feb 15 00:00:00 EST 2015}
}

The {delta}f particleincell method for gyrokinetic simulations with kinetic electrons and electromagnetic perturbations [Y. Chen, S. Parker, J. Comput. Phys. 189 (2003) 463] is extended to include arbitrary toroidal equilibrium profiles and fluxsurface shapes. The domain is an arbitrarily sized toroidal slice with periodicity assumed in toroidal direction. It is global radially and poloidally along the magnetic field. The differential operators and Jacobians are represented numerically which is a quite general approach with wide applicability. Discretization of the field equations is described. The issue of domain decomposition and particle load balancing is addressed. A derivation of the splitweight scheme ismore »

Electromagnetic formulation of global gyrokinetic particle simulation in toroidal geometry
The fluidkinetic hybrid electron model for global electromagnetic gyrokinetic particle simulations has been formulated in toroidal geometry using magnetic coordinates, providing the capabilities to describe low frequency processes in electromagnetic turbulence with electron dynamics. In the limit of long wavelength and no parallel electric field our equations reduce to the ideal magnetohydrodynamic equations. The formulation has been generalized to include equilibrium flows. The equations for zonal components of electrostatic and vector potentials have been derived, demonstrating the electron screening of the zonal vector potential. 
Gyrokinetic particle simulation of driftcompressional modes in dipole geometry
Gyrokinetic particle simulation of low frequency compressional modes has been developed using flux coordinates in the global magnetic dipole geometry. The compressional component is formulated in a scalar form of the parallel magnetic perturbation, and the gyroaveraging is performed explicitly in the configuration space. A reduced gyrokinetic model, in which the compressional perturbations are decoupled from the shear Alfven and electrostatic perturbations, has been implemented. Linear simulation results have been verified using a numerical Nyquist analysis of the dispersion relation in the slab limit. Global simulations of unstable driftcompressional modes in the dipole geometry with kinetic ions find that finitemore » 
Gyrokinetic simulations of mesoscale energetic particledriven Alfvenic turbulent transport embedded in microturbulence
Energetic particle (EP) transport from local highn toroidal Alfven eigenmodes (TAEs) and energetic particle modes (EPMs) is simulated with a gyrokinetic code. Linear and nonlinear simulations have identified a parameter range where the longwave TAE and EPM are unstable alongside the wellknown iontemperaturegradient (ITG) and trappedelectronmode (TEM) instabilities. A new eigenvalue solver in GYRO facilitates this mode identification. States of nonlinearly saturated local TAE/EPM turbulent intensity are identified, showing a 'soft' transport threshold for enhanced energetic particle transport against the TAE/EPM drive from the EP pressure gradient. The very longwavelength (mesoscale) TAE/EPM transport is saturated partially by nonlinear interaction withmore » 
Magnetic stochasticity in gyrokinetic simulations of plasma microturbulence
Analysis of the magnetic field structure from electromagnetic simulations of tokamak ion temperature gradient turbulence demonstrates that the magnetic field can be stochastic even at very low plasma pressure. The degree of magnetic stochasticity is quantified by evaluating the magnetic diffusion coefficient. We find that the magnetic stochasticity fails to produce a dramatic increase in the electron heat conductivity because the magnetic diffusion coefficient remains small.