Drift-kinetic simulation of neoclassical transport with impurities in tokamaks

Plasmas in modern tokamak experiments contain a significant fraction of impurity ions in addition to the main deuterium background ions. A new multiple ion-species δf particle simulation capability has been developed to self-consistently study the nonlocal effects of impurities on neoclassical transport in toroidal plasmas. A new algorithm for an unlike-particle collision operator, including test-particle and conserving field-particle parts, is described. Effects of the carbon impurity on the main deuterium species heat flux as well as an ambipolar radial electric field in a National Spherical Torus Experiment (NSTX) [M. Ono, S. M. Kaye, Y.-K. M. Peng, Nucl. Fusion 40, 557 (2000)] configuration were studied. A difference between carbon poloidal rotation found from simulation and from conventional theoretical estimates has been investigated and was identified to be a nonlocal finite orbit effect. In the case of large-aspect ratio tokamak configurations with steep toroidal flow profiles, we propose a theoretical model to describe this nonlocal effect. The dominant mechanisms captured by the model are associated with ion parallel velocity modification due to steep toroidal flow and radial electric field profiles. We present simulation results for carbon poloidal velocity in NSTX. Comparisons with neoclassical theory are discussed.