3-D nonlinear evolution of MHD instabilities

The nonlinear evolution of ideal MHD internal instabilities is investigated in straight cylindrical geometry by means of a 3-D initial-value computer code. These instabilities are characterized by pairs of velocity vortex cells rolling off each other and helically twisted down the plasma column. The cells persist until the poloidal velocity saturates at a few tenths of the Alfven velocity. The nonlinear phase is characterized by convection around these essentially fixed vortex cells. For example, the initially centrally peaked temperature profile is convected out and around to form an annulus of high temperature surrounding a small region of lower temperature. Weak, centrally localized instabilities do not alter the edge of the plasma. Strong, large-scale instabilities, resulting from a stronger longitudinal equilibrium current, drive the plasma against the wall. After three examples of instability are analyzed in detail, the numerical methods and their verification are discussed.