Nonlinear magnetohydrodynamic stability

The computer code developed by Bauer et al. (1978) for the study of the magnetohydrodynamic equilibrium and stability of a plasma in toroidal geometry is extended so that the growth rates of instabilities may be estimated more accurately. The original code, which is based on the variational principle of ideal magnetohydrodynamics, is upgraded by the introduction of a nonlinear formula for the growth rate of an unstable mode which acts as a quantitative measure of instability that is important in estimating numerical errors. The revised code has been applied to the determination of the nonlinear saturation, ballooning modes and beta limits for tokamaks, stellarators and torsatrons.