Fast axisymmetric stability calculations using variational techniques

A procedure for treating the axisymmetric (n = 0) stability of diverted plasmas in the presence of arbitrary, but toroidally symmetric, structures and active feedback circuits has been developed and implemented as a module in the TEQ free-boundary equilibrium code. This procedure is based on a variational solution of the ideal MHD normal mode equations. Inertia is ordered small but provides a constraint to allow the calculation of the poloidal and toroidal components of the plasma displacement. Feedback based on flux loop measurements is handled by introducing an adjoint system into the variational principle. Approximately 200 trial functions for the radial component of the plasma displacement and 200 magnetic surfaces are employed to obtain highly accurate estimates of the passive growth rate and the non-rigid eigenfunction. Nevertheless, the method is extremely fast: typically 10-20 sec of Cray 2 CPU time are required to analyze a realistic tokamak configuration. This speed, along with the direct coupling to the MHD equilibrium solver, allows interactive investigations of tokamak axisymmetric stability. Benchmarks with TSC and GATO are presented along with parameter scans for ITER and BPX. The results emphasize the importance of considering non-rigid mode effects which for ITER, yield higher nominal growth rates (non-rigid: 45 Hz, rigid: 25 Hz) and atypical internal inductance dependence (smaller li more unstable)