Non-Newtonian viscosity in magnetized plasma

Abstract

The particle and momentum balance equations can be solved on concentric circular flux surfaces to determine the effective viscous drag present in a magnetized tokamak plasma in the low aspect ratio limit. An analysis is developed utilizing the first-order Fourier expansion of the poloidal variation of quantities on the flux surface akin to that by Stacey and Sigmar [Phys. Fluids, 28, 9 (1985)]. Expressions to determine the poloidal variations of density, poloidal velocity, toroidal velocity, radial electric field, poloidal electric field, and other radial profiles are presented in a multi-species setting. Using as input experimental data for the flux surface averaged profiles of density, temperature, toroidal current, toroidal momentum injection, and the poloidal and toroidal rotations of at least one species of ion, one may solve the equations numerically for the remaining profiles. The resultant effective viscosities are compared to those predicted by Stacey and Sigmar and Shaing, et al., [Nuclear Fusion, 25, 4 (1985)]. A velocity dependence of the effective viscosity is observed, indicating the presence of non-Newtonian fluid behavior in magnetized plasma.