This project explores the idea of performing kinetic numerical simulations in the Z inner magnetically insulated transmission line (inner MITL) by reduced physics models such as a guiding center drift kinetic approximation for particles and electrostatic and magnetostatic approximation for the fields. The basic problem explored herein is the generation, formation, and evolution of vortices by electron space charge limited (SCL) emission. The results indicate that for relevant to Z values of peak current and pulse length, these approximations are excellent, while also providing tens to hundreds of times reduction in the computational load. The benefits could be enormous: Implementation of these reduced physics models in present particle-in-cell (PIC) codes could enable them to be routinely used for experimental design while still capturing essential non-thermal (kinetic) physics.
In this study, we examine the effects of the radiation reaction force on electrons in a radial magnetically insulated transmission line (MITL) near a load with peak currents of 60+ MA. More specifically, we study the differences in electron motion and kinetic energy with or without radiation reaction physics using a novel guiding center drift approach that incorporates E $\times$ B and ∇B drifts. A key finding of this study is that an electron's magnetic moment, which would be conserved when radiation reaction physics is not incorporated, can be significantly reduced in magnetic fields on the order of 10,000's T when radiation reaction is included. The reduction of magnetic moment gives rise to a significant reduction in cycloidal kinetic energy as well as a reduction in the electron's ∇B drift.
Bennett, N.; Welch, D.R.; Jennings, C.A.; Yu, E.; Hess, Mark H.; Hutsel, B.T.; Laity, G.; Moore, James M.; Rose, D.V.; Peterson, K.; Cuneo, M.E.
A challenge for the TW-class accelerators driving Z-pinch experiments, such as Sandia National Laboratories' Z machine, is to efficiently couple power from multiple storage banks into a single multi-MA transmission line. The physical processes that lead to current loss are identified in new large-scale, multidimensional simulations of the Z machine. Kinetic models follow the range of physics occurring during a pulse, from vacuum pulse propagation to charged-particle emission and magnetically-insulated current flow to electrode plasma expansion. Simulations demonstrate that current is diverted from the load through a combination of standard transport (uninsulated charged-particle flows) and anomalous transport. Standard transport occurs in regions where the electrode current density is a few 104-105 A/cm2 and current is diverted from the load via transport without magnetic insulation. In regions with electrode current density >106 A/cm2, electrode surface plasmas develop velocity-shear instabilities and a Hall-field-related transport which scales with electron density and may, therefore, lead to increased current loss.
The Magnetized Liner Inertial Fusion concept (MagLIF) [Slutz et al., Phys. Plasmas 17, 056303 (2010)] is being studied on the Z facility at Sandia National Laboratories. Neutron yields greater than 1012 have been achieved with a drive current in the range of 17-18 MA and pure deuterium fuel [Gomez et al., Phys. Rev. Lett. 113, 155003 (2014)]. We show that 2D simulated yields are about twice the best yields obtained on Z and that a likely cause of this difference is the mix of material into the fuel. Mitigation strategies are presented. Previous numerical studies indicate that much larger yields (10-1000 MJ) should be possible with pulsed power machines producing larger drive currents (45-60 MA) than can be produced by the Z machine [Slutz et al., Phys. Plasmas 23, 022702 (2016)]. To test the accuracy of these 2D simulations, we present modifications to MagLIF experiments using the existing Z facility, for which 2D simulations predict a 100-fold enhancement of MagLIF fusion yields and considerable increases in burn temperatures. Experimental verification of these predictions would increase the credibility of predictions at higher drive currents.
A critical component of the magnetically driven implosion experiments at Sandia National Laboratories is the delivery of high-current, 10s of MA, from the Z pulsed power facility to a target. In order to assess the performance of the experiment, it is necessary to measure the current delivered to the target. Recent Magnetized Liner Inertial Fusion (MagLIF) experiments have included velocimetry diagnostics, such as PDV (Photonic Doppler Velocimetry) or Velocity Interferometer System for Any Reflector, in the final power feed section in order to infer the load current as a function of time. However, due to the nonlinear volumetrically distributed magnetic force within a velocimetry flyer, a complete time-dependent load current unfold is typically a time-intensive process and the uncertainties in the unfold can be difficult to assess. In this paper, we discuss how a PDV diagnostic can be simplified to obtain a peak current by sufficiently increasing the thickness of the flyer. This effectively keeps the magnetic force localized to the flyer surface, resulting in fast and highly accurate measurements of the peak load current. In addition, we show the results of experimental peak load current measurements from the PDV diagnostic in recent MagLIF experiments.
We have developed a physics-based transmission-line-circuit model of the Z pulsed-power accelerator. The 33-m-diameter Z machine generates a peak electrical power as high as 85 TW, and delivers as much as 25 MA to a physics load. The circuit model is used to design and analyze experiments conducted on Z. The model consists of 36 networks of transmission-line-circuit elements and resistors that represent each of Zs 36 modules. The model of each module includes a Marx generator, intermediate-energy-storage capacitor, laser-triggered gas switch, pulse-forming line, self-break water switches, and tri-plate transmission lines. The circuit model also includes elements that represent Zs water convolute, vacuum insulator stack, four parallel outer magnetically insulated vacuum transmission lines (MITLs), double-post-hole vacuum convolute, inner vacuum MITL, and physics load. Within the vacuum-transmission-line system the model conducts analytic calculations of current loss. To calculate the loss, the model simulates the following processes: (i) electron emission from MITL cathode surfaces wherever an electric-field threshold has been exceeded; (ii) electron loss in the MITLs before magnetic insulation has been established; (iii) flow of electrons emitted by the outer-MITL cathodes after insulation has been established; (iv) closure of MITL anode-cathode (AK) gaps due to expansion of cathode plasma; (v) energy loss to MITL conductors operated at high lineal current densities; (vi) heating of MITL-anode surfaces due to conduction current and deposition of electron kinetic energy; (vii) negative-space-charge-enhanced ion emission from MITL anode surfaces wherever an anode-surface-temperature threshold has been exceeded; and (viii) closure of MITL AK gaps due to expansion of anode plasma. The circuit model is expected to be most accurate when the fractional current loss is small. We have performed circuit simulations of 52 Z experiments conducted with a variety of accelerator configurations and load-impedance time histories. For these experiments, the apparent fractional current loss varies from 0% to 20%. Results of the circuit simulations agree with data acquired on 52 shots to within 2%.