Numerical study of equilibrium solutions for axisymmetric plasmas with toroidal flow obtained using Solovev approach

Solovev approach of finding equilibrium solutions, which was extended to include the vacuum solutions provided by Zheng, Wooton, and Solano, was found extremely useful for the purpose of shaping studies. Its extension to toroidal equilibria with a general plasma flow was examined theoretically in a companion paper by Chu, Hu and Guo. The only meaningful extension was found for plasmas with a pure toroidal rotation and with a constant Mach number. A set of functions SOLOVEV_ZWSm was obtained which fixed location of the magnetic axis for equilibria with quasi-constant current density profile, with toroidal flow at constant Mach number and with specific heat capacity 1. The set Solovev_ZWSm should have complete shaping capability for plasma shapes with positive curvature at the boundary; but not for plasmas with negative curvature boundary points, i.e. the doublets or bean shaped tokamaks. We report here extensive numerical studies showing the shaping capability of Solovev_ZWSm for plasmas with pure toroidal rotations, including the change in topology of the solution when the rotation mach number changes. Included plasma topology are the sphere (spheromaks); and the tokamaks (including the doublets).