A New Flux Coordinates-Based Solver for Fixed-Boundary Tokamak Equilibrium with Toroidal Flow

The plasma equilibrium plays a crucial role in nuclear fusion studies, serving as the foundation for various aspects of fusion research, including plasma stability, transport, and current drive. In this paper, a new Grad-Shafranov equation solver is developed for the fixed-boundary plasma equilibria with toroidal flow. This solver utilizes the pressure profile, safety factor profile (not current profile), and any two profiles of the toroidal angular velocity, plasma temperature, and square of the Mach number as inputs. The numerical results obtained by this solver exhibit good agreement with known analytic solution under identical parameters, and the potential applications of the solver are demonstrated through several numerical equilibria with toroidal flow. It is very convenient to apply this code to simulate the tokamak equilibrium with a smooth plasma shape. In addition, the effect of toroidal flow on the plasma equilibria is investigated as a simple application. The results reveal a notable outward shift in the contour profiles of magnetic flux surface, density, pressure, and temperature induced by toroidal flow.