Numerical solutions of magnetohydrodynamic stability of axisymmetric toroidal plasmas using cubic B-spline finite element method

Conference · Thu Dec 01 04:00:00 EST 1988

OSTI ID:6439245

Cheng, C Z

A nonvariational ideal MHD stability code (NOVA) has been developed. In a general flux coordinate (/psi/, theta, /zeta/) system with an arbitrary Jacobian, the NOVA code employs Fourier expansions in the generalized poloidal angle theta and generalized toroidal angle /zeta/ directions, and cubic-B spline finite elements in the radial /psi/ direction. Extensive comparisons with these variational ideal MHD codes show that the NOVA code converges faster and gives more accurate results. An extended version of NOVA is developed to integrate non-Hermitian eigenmode equations due to energetic particles. The set of non-Hermitian integro-differential eigenmode equations is numerically solved by the NOVA-K code. We have studied the problems of the stabilization of ideal MHD internal kink modes by hot particle pressure and the excitation of ''fishbone'' internal kink modes by resonating with the energetic particle magnetic drift frequency. Comparisons with analytical solutions show that the values of the critical ..beta../sub h/ from the analytical theory can be an order of magnitude different from those computed by the NOVA-K code. 24 refs., 11 figs., 1 tab.

• Research Organization: Princeton Univ., NJ (USA). Plasma Physics Lab.
• DOE Contract Number: AC02-76CH03073
• OSTI ID: 6439245
• Report Number(s): PPPL-2575; CONF-890479-2; ON: DE89004635
• Country of Publication: United States
• Language: English

Related Subjects

70 PLASMA PHYSICS AND FUSION TECHNOLOGY
700107* -- Fusion Energy-- Plasma Research-- Instabilities
AXIAL SYMMETRY
CLOSED PLASMA DEVICES
EQUILIBRIUM
FINITE ELEMENT METHOD
FUNCTIONS
INSTABILITY
MHD EQUILIBRIUM
NUMERICAL SOLUTION
PLASMA INSTABILITY
SPLINE FUNCTIONS
SYMMETRY
THERMONUCLEAR DEVICES
TOKAMAK DEVICES