Study of up–down poloidal density asymmetry of high- $Z$ impurities with the new impurity version of XGCa

Addition of multispecies impurity ions to the total-f gyrokinetic particle-in-cell code XGCa is reported, including a cross-verification of neoclassical physics against the NEO code. This new version of the neoclassical gyrokinetic code XGCa is used to benchmark and confirm the previous reduced-equation-based prediction that high- $Z$ impurity particles in the Pfirsch–Schlüter regime can exhibit a significant level of up–down poloidal asymmetry, through the large parallel friction force, and thus influence the radial plasma transport significantly. The study is performed in a plasma with weak toroidal rotation. In comparison, when the impurity particles are in the plateau regime, the up–down poloidal asymmetry becomes weak, with the parallel friction force becoming weaker than the parallel viscous force. It is also found that the linearization of the perturbed distribution function, based on the small poloidal asymmetry assumption, can become invalid. Using the numerical data from XGCa, each term in the parallel fluid force-balance equation have been analysed to find that both the main ions and the electrons respond to the poloidal potential variation adiabatically when the high- $Z$ tungsten possesses large poloidal variation.

[1]  Choong-Seock Chang,et al.  Full-f XGC1 gyrokinetic study of improved ion energy confinement from impurity stabilization of ITG turbulence , 2017 .

[2]  P. Helander.  Bifurcated neoclassical particle transport , 1998 .

[3]  Theoretical explanation of the poloidal asymmetry of impurity spectral line emission in collisional tokamak Plasmas , 1979 .

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

[5]  J. Rice,et al.  Parallel transport studies of high-Z impurities in the core of Alcator C-Mod plasmas , 2013 .

[6]  F. Hinton,et al.  Effect of finite aspect ratio on the neoclassical ion thermal conductivity in the banana regime , 1982 .

[7]  Alcator C-Mod Team,et al.  Tungsten impurity transport experiments in Alcator C-Mod to address high priority research and development for ITER , 2015 .

[8]  O. Sauter,et al.  Neoclassical conductivity and bootstrap current formulas for general axisymmetric equilibria and arbitrary collisionality regime , 1999 .

[9]  Virginie Grandgirard,et al.  Self-consistent gyrokinetic modeling of neoclassical and turbulent impurity transport , 2018 .

[10]  R. Churchill,et al.  Neutral recycling effects on ITG turbulence , 2017 .

[11]  A. Polevoi,et al.  W transport and accumulation control in the termination phase of JET H-mode discharges and implications for ITER , 2018, Plasma Physics and Controlled Fusion.

[12]  Robert Hager,et al.  Gyrokinetic neoclassical study of the bootstrap current in the tokamak edge pedestal with fully non-linear Coulomb collisions , 2016 .

[13]  F. Hinton,et al.  Neoclassical ion transport in rotating axisymmetric plasmas , 1985 .

[14]  Patrick H. Worley,et al.  A fast low-to-high confinement mode bifurcation dynamics in the boundary-plasma gyrokinetic code XGC1 , 2018 .

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

[16]  Bruce I. Cohen,et al.  A suitable boundary condition for bounded plasma simulation without sheath resolution , 1993 .

[17]  R. Hazeltine,et al.  Impurity transport in the collisional regime for large poloidal variations , 1980 .

[18]  Robert Hager,et al.  A new hybrid-Lagrangian numerical scheme for gyrokinetic simulation of tokamak edge plasma , 2016, J. Comput. Phys..

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

[20]  Patrick H. Worley,et al.  A fully non-linear multi-species Fokker-Planck-Landau collision operator for simulation of fusion plasma , 2016, J. Comput. Phys..

[21]  Matthew Logan Reinke.  Experimental tests of parallel impurity transport theory in Tokamak plasmas , 2011 .

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

[23]  J. Contributors,et al.  High Z neoclassical transport: Application and limitation of analytical formulae for modelling JET experimental parameters , 2018 .

[24]  Choong-Seock Chang,et al.  Gyroaveraging operations using adaptive matrix operators , 2018, 1802.09125.

[25]  C. Chang,et al.  Enhancement of neoclassical transport coefficients by a poloidal electric field in tokamaks , 1983 .

[26]  J. Contributors,et al.  The impact of poloidal asymmetries on tungsten transport in the core of JET H-mode plasmas , 2015 .

[27]  S. Lanthaler,et al.  Overview of the JET results in support to ITER , 2017, Nuclear Fusion.

[28]  J. Manickam,et al.  Numerical study of tokamak equilibria with arbitrary flow , 2004 .

[29]  Szymon Suckewer,et al.  Vertical poloidal asymmetries of low-Z element radiation in the PDX tokamak , 1983 .

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

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

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

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