skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

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

Conference ·
OSTI ID:6439245

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; TRN: 89-009283
Resource Relation:
Conference: 7. international conference on finite element methods in flow problems, Huntsville, AL, USA, 3 Apr 1989; Other Information: Portions of this document are illegible in microfiche products
Country of Publication:
United States
Language:
English