A multi-species collisional operator for full-F gyrokinetics

A linearized multi-species collision operator has been developed for an efficient implementation in gyrokinetic codes. This operator satisfies the main expected properties: particle, momentum, and energy conservation, and existence of an H-theorem. A gyrokinetic version is then calculated, which involves derivatives with respect to the gyrocenter position, parallel velocity, and magnetic momentum. An isotropic version in the velocity space can be constructed for the specific problem of trace impurities colliding with a main species. A simpler version that involves derivatives with parallel velocity only has been developed. This reduced version has been implemented in the GYSELA gyrokinetic code, and is shown to comply with particle, momentum, and energy conservation laws. Moreover, the interspecies relaxation rates for momentum and energy agree very well with the theoretical values.

[1]  Martin Greenwald,et al.  Quantitative comparison of experimental impurity transport with nonlinear gyrokinetic simulation in an Alcator C-Mod L-mode plasma , 2012 .

[2]  Peter J. Catto,et al.  Linearized gyro‐kinetic equation with collisions , 1977 .

[3]  Laurent Villard,et al.  Gyrokinetic simulations of turbulent transport , 2010 .

[4]  J. Connor,et al.  The neo-classical transport theory of a plasma with multiple ion species , 1973 .

[5]  Marshall N. Rosenbluth,et al.  Plasma Transport in Toroidal Confinement Systems , 1972 .

[6]  M. Barnes,et al.  Linearized model Fokker-Planck collision operators for gyrokinetic simulations. I. Theory , 2008, 0808.1300.

[7]  R. Neu,et al.  Local effects of ECRH on argon transport in L-mode discharges at ASDEX Upgrade , 2011 .

[8]  F. J. Casson,et al.  Impact of the background toroidal rotation on particle and heat turbulent transport in tokamak plasmas , 2009 .

[9]  J. Candy,et al.  Kinetic calculation of neoclassical transport including self-consistent electron and impurity dynamics , 2008 .

[10]  Jet Efda Contributors,et al.  First operation with the JET International Thermonuclear Experimental Reactor-like wall , 2013 .

[11]  Alain J. Brizard.  A guiding-center Fokker–Planck collision operator for nonuniform magnetic fields , 2004 .

[12]  C. Angioni,et al.  Neoclassical transport of heavy impurities with poloidally asymmetric density distribution in tokamaks , 2014 .

[13]  X. Garbet,et al.  Interplay between gyrokinetic turbulence, flows, and collisions: perspectives on transport and poloidal rotation. , 2009, Physical review letters.

[14]  A. Peeters,et al.  Non-adiabatic passing electron response and outward impurity convection in gyrokinetic calculations of impurity transport in ASDEX Upgrade plasmas , 2007 .

[15]  J. Contributors,et al.  Metal impurity transport control in JET H-mode plasmas with central ion cyclotron radiofrequency power injection , 2011 .

[16]  Nonlinear neoclassical transport in a rotating impure plasma with large gradients , 1999 .

[17]  Hideo Sugama,et al.  Linearized Model Collision Operators for Multiple Ion Species Plasmas , 2009 .

[18]  M. Ottaviani,et al.  Effects of density asymmetries on heavy impurity transport in a rotating tokamak plasma , 1998 .

[19]  K. Wenzel,et al.  Neoclassical analysis of impurity transport following transition to improved particle confinement , 1990 .

[20]  F. Hinton,et al.  Theory of plasma transport in toroidal confinement systems , 1976 .

[21]  Alain J. Brizard,et al.  Foundations of Nonlinear Gyrokinetic Theory , 2007 .

[22]  A. Peeters,et al.  Direction of impurity pinch and auxiliary heating in tokamak plasmas. , 2006, Physical review letters.

[23]  J. F. Clarke,et al.  Neoclassical transport theory of a multispecies plasma in the low collision frequency regime , 1976 .

[24]  T. Fülöp,et al.  Effect of poloidal asymmetry on the impurity density profile in tokamak plasmas , 2011 .

[25]  S. P. Hirshman,et al.  Approximate Fokker–Planck collision operator for transport theory applications , 1976 .

[26]  F. Imbeaux,et al.  Anomalous transport of light and heavy impurities in Tore Supra ohmic, weakly sawtoothing plasmas , 2009 .

[27]  Jeff M. Candy,et al.  Full linearized Fokker–Planck collisions in neoclassical transport simulations , 2011 .

[28]  C. Bourdelle,et al.  Experimental electron temperature gradient dependence of heavy impurity transport in fusion devices. , 2010, Physical review letters.

[29]  Laurent Villard,et al.  Synergy between ion temperature gradient turbulence and neoclassical processes in global gyrokinetic particle-in-cell simulations , 2012 .

[30]  X. Garbet,et al.  Neoclassical equilibrium in gyrokinetic simulations , 2009 .

[31]  Shuichi Takamura,et al.  Chapter 4: Power and particle control , 2007 .

[32]  Marshall N. Rosenbluth,et al.  Numerical simulation of ion temperature gradient driven modes in the presence of ion-ion collisions , 1990 .

[33]  J. Contributors,et al.  Tungsten transport in JET H-mode plasmas in hybrid scenario, experimental observations and modelling , 2014 .

[34]  S. P. Hirshman,et al.  Neoclassical transport of a multispecies toroidal plasma in various collisionality regimes , 1977 .

[35]  Y. R. Martin,et al.  Plasma wall interaction and its implication in an all tungsten divertor tokamak , 2007 .

[36]  C. Bourdelle,et al.  Fluid simulations of turbulent impurity transport , 2007 .

[37]  M. Barnes,et al.  Linearized model Fokker-Planck collision operators for gyrokinetic simulations. II. Numerical implementation and tests , 2008, 0809.3945.

[38]  J. Candy,et al.  Pfirsch–Schlüter neoclassical heavy impurity transport in a rotating plasma , 2014 .

[39]  S. Hirshman,et al.  Neoclassical transport of impurities in tokamak plasmas , 1981 .

[40]  A. Lenard,et al.  PLASMA OSCILLATIONS WITH DIFFUSION IN VELOCITY SPACE , 1958 .