Numerical calculations using the full MHD equations in toroidal geometry (Q1072850)

A computer code has been constructed that solves the full magnetohydrodynamic 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. Althought 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.