Intrinsic toroidal rotation driven by turbulent and neoclassical processes in tokamak plasmas from global total-f gyrokinetic simulations

Gyrokinetic tokamak plasmas can exhibit intrinsic toroidal rotation driven by the residual stress $\Pi^{\rm res}$. While most studies have attributed $\Pi^{\rm res}$ to the parallel-momentum flux from the turbulent $\boldsymbol{E}\times\boldsymbol{B}$ drift, the neoclassical drift-orbit contribution (denoted $\Pi_\parallel^D$) is often neglected. Here, we use the global total-$f$ gyrokinetic code XGC to study $\Pi^{\rm res}$ in the core and the edge of a DIII-D H-mode plasma. Numerical results showed that $\Pi_\parallel^D$ makes up a significant portion of $\Pi^{\rm res}$; specifically, in the core it is higher than the neoclassical level in the presence of turbulence, while in the edge it forms a counter-current peak even without turbulence. Using a recently developed ``orbit-flux'' formulation, we show that its higher-than-neoclassical level in the core is from the divergence of the turbulent flux, while its counter-current peak in the edge is mainly from the neoclassical collisional ion orbit loss. These results suggest that $\Pi_\parallel^D$ may be important for the study of intrinsic toroidal rotation.