This project has progressed several physics models in the EMPIRE plasma simulation code to achieve higher fidelity simulations of high-current diode operation in pulsed-power accelerators. In this report, we present details for the following major work products: (1) a set of verification problems covering all key processes involved in gap closure physics has been designed; this suite has facilitated feature vetting and overall model maturation, (2) a new EMPIRE exemplar has been developed: the Radiographic Integrated Test Stand 6 (RITS-6) diode, and (3) An exemplar for the Saturn accelerator exemplar was enabled by the models matured under this work to self-consistently simulate further into the diode pulse than previously possible (bipolar flow regime). These developments have lead to the highest confidence EMPIRE power flow predictions of the Saturn accelerator to date. Additionally, three modeling approaches for simulating electrode plasmas have been investigated. We report on these results and provide recommendations.
This milestone work baselines electromagnetic particle-in-cell capability of the EMPIRE plasma simulation code to model key processes germane to the physics of electrode plasmas arising in magnetically-insulated transmission lines operating at or near 20 MA. This evaluation is done so through the provision of benchmark verification problems designed to exercise the individual and combined physics models on a small-scale surrogate geometry for the final-feed-to-load region of the Z accelerator under representative operating conditions. In this report, we overview our test designs, and present a portfolio of simulation results along with performance assessments which altogether establish state-of-the-art. In particular, two main verification categories are covered this report: (1) Z-relevant desorption physics (Temkin isotherm), and (2) two approaches to simulate electrode plasma creation and dynamics (automatic creation versus self-consistent creation through direct simulation Monte Carlo collisions).
The Z accelerator at Sandia National Laboratories conducts z-pinch experiments at 26 MA in support of DOE missions in stockpile stewardship, dynamic materials, fusion, and other basic sciences. Increasing the current delivered to the z-pinch would extend our reach in each of these disciplines. To achieve increases in current and accelerator efficiency, a fraction of Z’s shots are set aside for research into transmission-line power flow. These shots, with supporting simulations and theory, are incorporated into this Advanced Diagnostics milestone report. The efficiency of Z is reduced as some portion of the total current is shunted across the transmission-line gaps prior to the load. This is referred to as “current loss”. Electrode plasmas have long been implicated in this process, so the bulk of dedicated power-flow experiments are designed to measure the plasma environment. The experimental analyses are enhanced by simulations conducted using realistic hardware and Z voltage pulses. In the same way that diagnostics are continually being improved for sensitivity and resolution, the modeling capability is continually being improved to provide faster and more realistic simulations. The specifics of the experimental hardware, diagnostics, simulations, and algorithm developments are provided in this report. The combined analysis of simulation and data confirms that electrode plasmas have the most detrimental impact on current delivery. Experiments over the last three years have tested the theoretical current-loss mechanisms of enhanced ion current, plasma gap closure, and Hall-related current. These mechanisms are not mutually exclusive and may be coincident in the final feed as well as in upstream transmission lines. The final-feed geometries tested here, however, observe lower-density plasmas without dominant ion currents which is consistent with a Hall-related current. The picture of plasma formation and transport formed from experiment and simulation is informing hardware designs being fielded on Z now and being proposed for the Next-Generation Pulsed Power (NGPP) facility. In this picture, the strong magnetic fields that heat the electrodes above particle emission thresholds also confine the charged particles near the surface. Some portion of the plasmas thus formed is transported into the transmission-line gap under the force of the electric field, with aid from plasma instabilities. The gap plasmas are then transported towards the load by a cross-field drift, where they accumulate and contribute to a likely Hall-related cross-gap current. The achievements in experimental execution, model validation, and physical analysis presented in this report set the stage for continued progress in power flow and load diagnostics on Z. The planned shot schedule for Z and Mykonos will provide data for extrapolation to higher current to ensure the predicted performance and efficiency of a NGPP facility.
We explore the use of reduced physics models for efficient kinetic particle simulations of space charge limited (SCL) emission in inner magnetically insulated transmission lines (inner MITLs), with application to Sandia National Laboratories' Z machine. We propose a drift kinetic (guiding center) model of electron motion in place of a fully kinetic model and electrostatic-magnetostatic fields in place of electromagnetic fields. The validity of these approximations is suggested by the operational parameters of the Z machine, namely, current pulse lengths of order 100 ns compared with Larmor periods typically smaller than 10-11 s, typical Larmor radii of a few (tens) of microns (magnetic fields of tens to hundreds of Tesla) compared with MITL dimensions of a few centimeters, and transient time of light waves along the inner MITL of order a fraction of a nanosecond. Guiding center orbits eliminate the fast electron gyromotion, which enables the use of tens to hundreds of times larger time steps in the numerical particle advance. Electrostatic-magnetostatic fields eliminate the Courant-Friedrichs-Lewy (CFL) numerical stability limit on the time step and allow the use of higher grid resolutions or, alternatively, larger time steps in the fields advance. Overall, potential computational cost savings of tens to hundreds of times exists. The applicability of the reduced physics models is examined on two problems. First, in the simulation of space charge limited emission of electrons from the cathode surface due to high electric fields in a radial inner MITL geometry with a short load. In particular, it is shown that a drift kinetic-based particle-in-cell (PIC) model with electrostatic-magnetostatic fields is able to accurately reproduce well-known physics of electron vortex formation, spatially and temporally. Second, deeper understanding is gained of the mechanism behind vortex formation in this MITL geometry by considering an exemplar problem of an electron block of charge. This simpler setup reveals that the main mechanism of vortex formation can be attributed to pure drift motion of the electrons, that is, the (fully kinetic) gyromotion of the electrons is inessential to the process. This exemplar problem also suggests a correlation of the spatial dimensions of vortices to the thickness of the electron layer, as observed in SCL simulations. It also confirms that the electromagnetic nature of the fields does not play an essential role. Finally, an improved hybrid fully kinetic and drift kinetic model for electron motion is proposed, as means of capturing finite Larmor radius (FLR) effects; the particular FLR physics that is missed by the drift kinetic model is the particle-wall interaction. By initializing SCL emitted electrons as fully kinetic and later transitioning them to drift kinetic, according to simple criteria, the accuracy of SCL simulations can be improved, while preserving the potential for computational efficiency.
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 article, we derive the vacuum electric fields within specific cylindrically symmetric magnetically insulated transmission lines (MITLs) in the limit of an infinite speed of light for an arbitrary time-dependent current. We focus our attention on two types of MITLs: the radial MITL and a spherically curved MITL. We then simulate the motion of charged particles, such as electrons, present in these MITLs due to the vacuum fields. In general, the motion of charged particles due to the vacuum fields is highly nonlinear since the fields are nonlinear functions of spatial coordinates and depend on an arbitrary time-dependent current drive. Using guiding center theory, however, one can describe the gross particle kinetics using a combination of $\textbf {E} \times \textbf {B}$ and $\nabla B$ drifts. In addition, we compare our approximate inner MITL field models and particle kinetics with those from a fully electromagnetic simulation code. We find that the agreement between the approximate model and the electromagnetic simulations is excellent.
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, Nicole; Cuneo, Michael E.; Yu, Edmund; Jennings, Christopher A.; Laity, George; Hutsel, Brian T.; Peterson, Kyle; Welch, Dale R.; Rose, David V.; Hess, Mark H.; Moore, James M.
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.
We report that 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. Additionally, we show the results of experimental peak load current measurements from the PDV diagnostic in recent MagLIF experiments.