CALCULATION OF 3-DIMENSIONAL MHD EQUILIBRIA WITH ISLANDS AND STOCHASTIC REGIONS

A three-dimensional MHD equilibrium code is described that does not assume the existence of good flux surfaces. Given an initial guess for the magnetic field, the code proceeds by calculating the pressure driven current and then by updating the field using Ampere's law. The numerical algorithm to solve the magnetic differential equation for the pressure driven current is described, and demonstrated for model fields having islands and stochastic regions. The numerical algorithm which solves Ampere's law in three dimensions is also described. Finally, the convergence of the code is illustrated for a particular stellarator equilibrium with no large islands.