Imploding to equilibrium of helically symmetric theta pinches

The time-dependent, single-fluid, dissipative magnetohydrodynamic equations are solved in helical coordinates (r,phi), where phi = THETA-kz, k = 2π/L and L is the periodicity length in the z-direction. The two-dimensional numerical calculations simulate theta pinches which have an l = 1 helical field added to them. Given the applied magnetic fields and the initial state of the plasma, we study the time evolution of the system. The plasma is found to experience two kinds of oscillations, occurring on different time scales. These are the radial compression oscillations, and the slower helical oscillations of the plasma column. The plasma motion is followed until these oscillations disappear and an equilibrium is nearly reached. Hence given the amplitude and the rise time of the applied magnetic fields, the calculations allow one to relate the initial state of a cold, homogeneous plasma to its final equilibrium state where it is heated and compressed