Full radius linear and nonlinear gyrokinetic simulations for tokamaks and stellarators:: zonal flows, applied E x B flows, trapped electrons and finite beta

被引:64
作者
Villard, L [1 ]
Allfrey, SJ
Bottino, A
Brunetti, M
Falchetto, GL
Grandgirard, V
Hatzky, R
Nührenberg, J
Peeters, AG
Sauter, O
Sorge, S
Vaclavik, J
机构
[1] Ecole Polytech Fed Lausanne, Assoc Euratom Confederat Suisse, Ctr Rech Phys Plasmas, CH-1007 Lausanne, Switzerland
[2] Euratom CEA, Dept Rech Fus Controlee, Cadarache, France
[3] Max Planck Gesell, Rechenzentrum, Garching, Germany
[4] Euratom IPP, Max Planck Inst Plasmaphys, Greifswald, Germany
[5] Euratom IPP, Max Planck Inst Plasmaphys, Garching, Germany
关键词
D O I
10.1088/0029-5515/44/1/019
中图分类号
O35 [流体力学]; O53 [等离子体物理学];
学科分类号
070204 ; 080103 ; 080704 ;
摘要
The aim of this paper is to report on recent advances made in global gyrokinetic simulations of ion temperature gradient (ITG) modes and other microinstabilities. The nonlinear development and saturation of ITG modes and the role of E x B zonal flows are studied with a global nonlinear deltaf formulation that retains parallel nonlinearity and thus allows for a check of the energy conservation property as a means of verifying the quality of the numerical simulation. Due to an optimized loading technique, the conservation property is satisfied with an unprecedented quality well into the nonlinear stage. The zonal component of the perturbation evolves to a quasi-steady state with regions of ITG suppression, strongly reduced radial energy flux and steepened effective temperature profiles alternating with regions of higher ITG mode amplitudes, larger radial energy flux and flattened effective temperature profiles. A semi-Lagrangian approach free of statistical noise is proposed as an alternative to the nonlinear deltaf formulation. An ASDEX-Upgrade experiment with an internal transport barrier is analysed with a global gyrokinetic code that includes trapped electron dynamics. The weakly destabilizing effect of trapped electron dynamics on ITG modes in an axisymmetric bumpy configuration modelling W7-X is shown in global linear simulations that retain the full electron dynamics. Finite beta effects on microinstabilities are investigated with a linear global spectral electromagnetic gyrokinetic formulation. The radial global structure of electromagnetic modes shows a resonant behaviour with rational q values.
引用
收藏
页码:172 / 180
页数:9
相关论文
共 31 条
[1]   The role of radial electric fields in linear and nonlinear gyrokinetic full radius simulations [J].
Allfrey, SJ ;
Bottino, A ;
Sauter, O ;
Villard, L .
NEW JOURNAL OF PHYSICS, 2002, 4 :29.1-29.19
[2]   A revised δf algorithm for nonlinear PIC simulation [J].
Allfrey, SJ ;
Hatzky, R .
COMPUTER PHYSICS COMMUNICATIONS, 2003, 154 (02) :98-104
[3]  
ALLFREY SJ, 2001, P 28 EPS C CONTR F A, V25
[4]  
ALLFREY SJ, 2002, P 29 EPS C PLASM P A, V26
[5]  
ALLFREY SJ, 2002, P JOINT VAR LAUS INT, P171
[6]  
BOTTINO A, 2002, THEORY FUSION PLASMA, P351
[7]  
BOTTINO A, 2002, P 29 EPS C PLASM P A, V26
[8]  
BRUNETTI M, 2002, P 29 EPS C PLASM P A, V26
[9]   Global approach to the spectral problem of microinstabilities in tokamak plasmas using a gyrokinetic model [J].
Brunner, S ;
Fivaz, M ;
Tran, TM ;
Vaclavik, J .
PHYSICS OF PLASMAS, 1998, 5 (11) :3929-3949
[10]   PARTIALLY LINEARIZED ALGORITHMS IN GYROKINETIC PARTICLE SIMULATION [J].
DIMITS, AM ;
LEE, WW .
JOURNAL OF COMPUTATIONAL PHYSICS, 1993, 107 (02) :309-323