TERPSICHORE code : linear MHD stability analysis of 3D equilibria
A new 3D linear MHD stability code generalizing ERATO and PEST has been developed . From nested flux surface 3D equilibria it is mapped to Boozer coordinates. Finite Fourier series are used in the two periodic coordinates and finite hybrid elements across the flux surfaces. An eigen value problem in a special block diagonal matrix form is obtained from the MHD energy principle. Two direct matrix solvers are developed, to give the eigen values and eigen vectors using inverse vector iteration. The first is optimized for a 2 processor Cray-2. The second is optimal for machines having 4 or more processors. Constructing the matrix elements is, even when highly optimized, more time consuming than solving the eigensystem. Realistic problems will require up to 100 megawords of memory and 100 minutes of Cray-2 time.
Bibliographic Reference: Paper presented: 30th Annual Meeting, Division of Plasma Physics, organised by the American Physical Society, Hollywood (US), Oct. 31 - Nov. 4, 1988 Article: Bulletin of the American Physical Society, Vol. 33(1988) No.9
Availability: Available from (1) as Paper EN 34287 ORA
Record Number: 198910214 / Last updated on: 1994-12-01
Original language: en
Available languages: en