Skip to main content
U.S. Department of Energy
Office of Scientific and Technical Information

Numerical calculations using the full MHD equations in toroidal geometry

Journal Article · · J. Comput. Phys.; (United States)

A computer code has been constructed that solves the full magnetohydrodynamic (MHD) equations in toroidal geometry. The code is applicable to toroidal devices, including tokamaks, stellarators, and reversed field pinches. A fully implicit numerical technique is used that allows linear eigenvalues and eigenfunctions to be found in a very few computational steps. Although the present work describes the solution of the linearized equations, generalization of the numerical method to the solution of the nonlinear problem is straightforward. Use of the code is illustrated by calculating the n = 1 instability for a tokamak configuration. The results show the structural changes in the eigenfunctions as the plasma pressure is increased.

Research Organization:
Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831
DOE Contract Number:
AC05-84OR21400
OSTI ID:
5871047
Journal Information:
J. Comput. Phys.; (United States), Journal Name: J. Comput. Phys.; (United States) Vol. 63:1; ISSN JCTPA
Country of Publication:
United States
Language:
English