Driven, steady-state RFP computations

Abstract

The pseudospectral three-dimensional MHD code of Dahlburg et al. (1986 and 1987) is used to compute the dynamical behavior of a channel of magnetofluid carrying an axial current and magnetic flux. This situation contains the essential MHD behavior of the reversed-field pinch (RFP). An externally imposed electric field is applied to an initially current-free magnetofluid and drives currents that rise and eventually fluctuate about values corresponding to pinch ratios Theta of about 1.3, 2.2, and 4.5. A period of violent turbulence leads to an approximately force-free core, surrounded by an active MHD boundary layer that is not force-free. A steady state is reached that can apparently be sustained indefinitely (for several hundred Alfven transit times or longer). The turbulence level and time variability in the steady state increase with increasing Theta. The average toroidal magnetic field at the wall reverses for Theta = 2.2 and 4.5, but not for Theta = 1.3. Negative toroidal current filaments are observed. The Lundquist numbers are of the order of a few hundred.