Neoclassical Simulation of Tokamak Plasmas using Continuum Gyrokinetc Code TEMPEST

Journal Article · Fri Nov 09 04:00:00 EST 2007 · Physical Review E, vol. 78, N/A, July 29, 2008, pp. 016406

OSTI ID:941379

Xu, X Q

We present gyrokinetic neoclassical simulations of tokamak plasmas with self-consistent electric field for the first time using a fully nonlinear (full-f) continuum code TEMPEST in a circular geometry. A set of gyrokinetic equations are discretized on a five dimensional computational grid in phase space. The present implementation is a Method of Lines approach where the phase-space derivatives are discretized with finite differences and implicit backwards differencing formulas are used to advance the system in time. The fully nonlinear Boltzmann model is used for electrons. The neoclassical electric field is obtained by solving gyrokinetic Poisson equation with self-consistent poloidal variation. With our 4D ({psi}, {theta}, {epsilon}, {mu}) version of the TEMPEST code we compute radial particle and heat flux, the Geodesic-Acoustic Mode (GAM), and the development of neoclassical electric field, which we compare with neoclassical theory with a Lorentz collision model. The present work provides a numerical scheme and a new capability for self-consistently studying important aspects of neoclassical transport and rotations in toroidal magnetic fusion devices.

• Research Organization: Lawrence Livermore National Laboratory (LLNL), Livermore, CA
• Sponsoring Organization: USDOE
• DOE Contract Number: W-7405-ENG-48
• OSTI ID: 941379
• Report Number(s): UCRL-JRNL-236541
• Journal Information: Physical Review E, vol. 78, N/A, July 29, 2008, pp. 016406, Journal Name: Physical Review E, vol. 78, N/A, July 29, 2008, pp. 016406 Vol. 78
• Country of Publication: United States
• Language: English

Related Subjects

70 PLASMA PHYSICS AND FUSION TECHNOLOGY
ELECTRIC FIELDS
ELECTRONS
GEOMETRY
HEAT FLUX
IMPLEMENTATION
PHASE SPACE
POISSON EQUATION
SIMULATION
TRANSPORT