The sheared-flow stabilized (SFS) Z-pinch has demonstrated long-lived plasmas with fusion-relevant parameters. This Letter presents the first experimental results demonstrating sustained, quasi-steady-state neutron production from the Fusion Z-pinch Experiment (FuZE), operated with a mixture of 20% deuterium/80% hydrogen by volume. Neutron emissions lasting approximately 5 µs are reproducibly observed with pinch currents of approximately 200 kA during an approximately 16 µs period of plasma quiescence. The average neutron yield is estimated to be (1.25 ± 0.45) × 10 5 neutrons/pulse and scales with the square of the deuterium concentration. Coincident with the neutron signal, plasma temperatures of 1 − 2 keV, and densities of approximately 10 17 cm −3 with 0.3 cm pinch radii are measured with fully-integrated diagnostics.

The first fully kinetic particle-in-cell (PIC) simulations of sheared flow stabilized Z-pinch plasmas show the suppression of the sausage instability by shear, ∂rvz ≠ 0, with flow Mach numbers ≲1, consistent with experimental observations. Experimental investigations of sheared-flow stabilized Z-pinches demonstrated stability for 10 s of microseconds, over 1000 Alfvén radial transit times, in quasi steady-state plasmas that are an intermediate between conventional inertial and magnetic confinement systems. The observed stability coincides with the presence of radial shear in axial flow profiles with peak speeds less than Mach 1, and experiments are underway to validate scaling this design to fusion conditions. The experimentally observed stability agrees with models of m = 1 kink mode suppression by sheared flows, but existing models of the m = 0 sausage mode underestimate the efficacy of sheared flow stabilization. These models rely on fluid approximations and find that stabilization requires flows ranging from Mach 1.7 to 4.3, and in some cases, stabilization is not reproduced in the models. This is faster than the measured flows in long-lived plasmas and would necessitate substantial energy convection out of the Z-pinch and the need to drive and sustain supersonic flows in future devices. The MHD models typically used in the literature are invalid in the high-temperature, high-current environments desirable for many Z-pinch applications, and they ignore large Larmor radius effects and viscous dissipation which are known to impact Z-pinch stability. PIC simulations can capture all these effects as well as kinetic instabilities that could influence the performance of high-temperature sheared flow stabilized Z-pinch plasmas. The PIC simulations presented here show the suppression and damping of m = 0 modes by sheared flows ∂rvz = 0.75vA/r0 with flow Mach numbers ≲1. Equivalent stability occurs under plasma conditions ranging from the limits of present-day experimental capabilities to the projected conditions of a sheared flow stabilized Z-pinch reactor.

A kinetic electrostatic eigenvalue equation for the lower-hybrid drift instability (LHDI) in a thin Harris current sheet with a guide field is derived based on the gyrokinetic electron and fully kinetic ion(GeFi) description. Three-dimensional nonlocal analyses are carried out to investigate the influence of a guide field on the stabilization of the LHDI by finite parallel wavenumber, k∥. Detailed stability properties are first analyzed locally, and then as a nonlocal eigenvalue problem. Our results indicate that at large equilibrium drift velocities, the LHDI is further destabilized by finite k∥ in the short-wavelength domain. This is demonstrated in a local stability analysis and confirmed by the peak in the eigenfunction amplitude. We find the most unstable modes localized at the current sheet edges, and our results agree well with simulations employing the GeFi code developed by Lin et al. [Plasma Phys. Controlled Fusion 47, 657 (2005); Plasma Phys. Controlled Fusion 53, 054013 (2011)].

The eigenmode stability properties of three-dimensional lower-hybrid-drift-instabilities (LHDI) in a Harris current sheet with a small but finite guide magnetic field have been systematically studied by employing the gyrokinetic electron and fully kinetic ion (GeFi) particle-in-cell (PIC) simulation model with a realistic ion-to-electron mass ratio mi/me. In contrast to the fully kinetic PIC simulation scheme, the fast electron cyclotron motion and plasma oscillations are systematically removed in the GeFi model, and hence one can employ the realistic mi/me. The GeFi simulations are benchmarked against and show excellent agreement with both the fully kinetic PIC simulation and the analytical eigenmode theory. Our studies indicate that, for small wavenumbers, ky, along the current direction, the most unstable eigenmodes are peaked at the location where k→·B→=0, consistent with previous analytical and simulation studies. Here, B→ is the equilibrium magnetic field and k→ is the wavevector perpendicular to the nonuniformity direction. As ky increases, however, the most unstable eigenmodes are found to be peaked at k→·B→≠0. In addition, the simulation results indicate that varying mi/me, the current sheet width, and the guide magnetic field can affect the stability of LHDI. Simulations with the varying mass ratio confirm the lower hybrid frequency and wave number scalings.