Linear ideal magnetohydrodynamic stability computations in three dimensional plasmas
The 3-D ideal MHD stability code TERPSICHORE has been developed to investigate the linear ideal magnetohydrodynamic (MHD) stability properties of complicated plasma confinement configurations with the aim of identifying new devices that could become attractive candidates as power producing fusion reactors. The variational energy principle that is minimised is Fourier-decomposed in poloidal and toroidal angular variables that are optimal for the stability calculation. Finite hybrid elements are used in the radial discretisation of the energy integral. In the weak form, the stability problem reduces to a special block pentadiagonal eigenvalue matrix equation which is solved using an inverse vector iteration approach to determine the eigenvalue closest to the null value. All the eigenvalues of the system can be determined via the eigenvalue shift scheme. An application to the advanced toroidal facility, a 12 field period torsatron in operation in Oak Ridge, USA, demonstrates that the device becomes unstable to internal ideal MHD modes with toroidal mode number n = 3 when the plasma pressure reaches about half the value that it was predicted and designed to achieve using standard profiles for the pressure and zero net toroidal current.
Bibliographic Reference: Paper presented: International IMACS Conference on Mathematical Modelling and Applied Mathematics, Moscow (SU), June 18-23, 1990
Availability: Available from (1) as Paper EN 36156 ORA
Record Number: 199110756 / Last updated on: 1994-12-02
Original language: en
Available languages: en