Intrinsic flow in plasma physics is a long-standing puzzle, since it is difficult to understand its origin without contradiction to momentum conservation in conventional wisdom. It is proved that the electromagnetic turbulent acceleration as a candidate for intrinsic parallel flow generation driven by pressure gradient along the total magnetic field line does not contradict momentum conservation. The conserved quantity corresponding to axial symmetry is the total gyrocenter parallel canonical momentum carried by both species or the total gyrocenter parallel momentum including the ion gyrocenter kinematic momentum and electromagnetic fields momentum, but not the ion kinematic momentum, or even the ion parallel flow. A conservation equation of total parallel momentum including the ion particles' kinematic momentum and electromagnetic fields momentum is also presented.

The quasilinear intrinsic parallel flow drive including parallel residual stress, kinetic stress, cross Maxwell stress and parallel turbulent acceleration by electromagnetic ion temperature gradient (ITG) turbulence is calculated analytically using electromagnetic gyrokinetic theory. Both the kinetic stress and cross Maxwell stress also enter the mean parallel flow velocity equation via their divergence, as for the usual residual stress. The turbulent acceleration driven by ion pressure gradient along the total magnetic field (including equilibrium magnetic field and fluctuating radial magnetic field) cannot be written as a divergence of stress, and so should be treated as a local source/sink. All these terms can provide intrinsic parallel rotation drive.Electromagnetic effects reduce the non-resonant electrostatic stress force and even reverse it, but enhance the resonant stress force. Both the non-resonant and resonant turbulent acceleration terms are also enhanced by electromagnetic effects. The possible implications of our results for experimental observations are discussed.

Both the parallel residual stress and parallel turbulent acceleration driven by electrostatic collisionsless trapped electron mode (CTEM) turbulence are calculated analytically using

The mean parallel current density evolution equation is presented using electromagnetic (EM) gyrokinetic equation. There exist two types of intrinsic current driving mechanisms resulted from EM electron temperature gradient (ETG) turbulence. The first type is the divergence of residual turbulent flux including a residual stress-like term and a kinetic stress-like term. The second type is named as residual turbulent source, which is driven by the correlation between density and parallel electric field fluctuations. The intrinsic current density driven by the residual turbulent source is negligible as compared to that driven by the residual turbulent flux.The ratio of intrinsic current density driven by EM ETG turbulence to the background bootstrap current density is estimated. The local intrinsic current density driven by the residual turbulent flux for mesoscale variation of turbulent flux can reach about 80% of the bootstrap current density in the core region of ITER standard scenario, but there is no net intrinsic current on a global scale. Based on this, the local intrinsic current driven by EM micro-turbulence and its effects on local modification of the profile of safety factor may be needed to be carefully taken into account in the future device with high which is the ratio between electron pressure to the magnetic pressure.

Understanding the generation of intrinsic rotation in tokamak plasmas is crucial for future fusion reactors such as ITER. We proposed a new mechanism named turbulent acceleration for the origin of the intrinsic parallel rotation based on gyrokinetic theory. The turbulent acceleration acts as a local source or sink of parallel rotation, i.e., volume force, which is different from the divergence of residual stress, i.e., surface force. However, the order of magnitude of turbulent acceleration can be comparable to that of the divergence of residual stress for electrostatic ion temperature gradient (ITG) turbulence. A possible theoretical explanation for the experimental observation of electron cyclotron heating induced decrease of co-current rotation was also proposed via comparison between the turbulent acceleration driven by ITG turbulence and that driven by collisionless trapped electron mode turbulence. We also extended this theory to electromagnetic ITG turbulence and investigated the electromagnetic effects on intrinsic parallel rotation drive. Finally, we demonstrated that the presence of turbulent acceleration does not conflict with momentum conservation.