Publications

Abstract not provided.

Abstract not provided.

Abstract not provided.

Kinetic simulations of Sandia National Laboratories' Z machine are conducted to understand particle transport in the highly magnetized environment of a multi-MA accelerator. Joule heating leads to the rapid formation of electrode surface plasmas. These plasmas are implicated in reducing accelerator efficiency by diverting current away from the load [M.R. Gomez et al., Phys. Rev. Accel. Beams 20, 010401 (2017)PRABCJ2469-988810.1103/PhysRevAccelBeams.20.010401, N. Bennett et al., Phys. Rev. Accel. Beams 22, 120401 (2019)PRABCJ2469-988810.1103/PhysRevAccelBeams.22.120401]. The fully-relativistic, electromagnetic simulations presented in this paper show that particles emitted in a space-charge-limited manner, in the absence of plasma, are magnetically insulated. However, in the presence of plasma, particles are transported across the magnetic field in spite of being only weakly collisional. The simulated cross-gap currents are well-approximated by the Hall current in the generalized Ohm's law. The Hall conductivities are calculated using the simulated particle densities and energies, and the parameters that increase the Hall current are related to transmission line inductance. Analogous to the generalized Ohm's law, we extend the derivation of the magnetized diffusion coefficients to include the coupling of perpendicular components. These yield a Hall diffusion rate, which is equivalent to the empirical Bohm diffusion.

The rudiments of a particle-based single-fluid two-temperature magnetohydrodynamic (MHD) algorithm have been outlined in Thoma et al. (2013). The extension of this algorithm to include the effect of Hall physics is described in this paper. An implicit leapfrog version of the algorithm, which allows timesteps large compared to the resistive decay time and other relevant timescales, has recently been added to a hybrid particle-in-cell code. In standard MHD the Hall term in the generalized Ohm’s law can often be neglected when the Hall parameter is small. This term must, however, be retained in regimes where it is non-negligible. The retention of displacement current in Maxwell’s equations avoids the numerical difficulties associated with the whistler mode, which are encountered in standard explicit Hall-MHD codes, and allows the algorithm to be incorporated into hybrid particle-in-cell codes, for which particles may migrate from a kinetic to fluid to MHD description based upon local ambient plasma conditions. A highly-coupled implicit Hall-MHD formalism is presented, in which displacement current can either be retained or neglected. Even when displacement current is neglected, the highly-coupled implicit formalism avoids the restrictive timesteps for the whistler mode in explicit Hall-MHD codes. A comparison of numerical and analytic dispersion analysis demonstrates the feasibility of this approach and establishes relevant constraints to assure numerical stability. The implementation of the algorithm is described, and test simulation results in 1D and 2D in both linear and nonlinear regimes are presented.

This report describes the high-level accomplishments from the Plasma Science and Engineering Grand Challenge LDRD at Sandia National Laboratories. The Laboratory has a need to demonstrate predictive capabilities to model plasma phenomena in order to rapidly accelerate engineering development in several mission areas. The purpose of this Grand Challenge LDRD was to advance the fundamental models, methods, and algorithms along with supporting electrode science foundation to enable a revolutionary shift towards predictive plasma engineering design principles. This project integrated the SNL knowledge base in computer science, plasma physics, materials science, applied mathematics, and relevant application engineering to establish new cross-laboratory collaborations on these topics. As an initial exemplar, this project focused efforts on improving multi-scale modeling capabilities that are utilized to predict the electrical power delivery on large-scale pulsed power accelerators. Specifically, this LDRD was structured into three primary research thrusts that, when integrated, enable complex simulations of these devices: (1) the exploration of multi-scale models describing the desorption of contaminants from pulsed power electrodes, (2) the development of improved algorithms and code technologies to treat the multi-physics phenomena required to predict device performance, and (3) the creation of a rigorous verification and validation infrastructure to evaluate the codes and models across a range of challenge problems. These components were integrated into initial demonstrations of the largest simulations of multi-level vacuum power flow completed to-date, executed on the leading HPC computing machines available in the NNSA complex today. These preliminary studies indicate relevant pulsed power engineering design simulations can now be completed in (of order) several days, a significant improvement over pre-LDRD levels of performance.

Abstract not provided.

Abstract not provided.

Two-dimensional electromagnetic (EM) particle-in-cell (PIC) simulations of a radial magnetically-insulated-transmission-line are presented and compared to the model of E. M. Waisman, M. P. Desjarlais, and M. E. Cuneo [Phys. Rev. Accel. Beams 22, 030402 (2019) in the “high-enhancement” (WDC-HE) limit. The simulations use quasi-equilibrium current and voltage values based on the Sandia National Laboratories Z accelerator, with prescribed injection of an electron sheath that gives electron density profiles qualitatively similar to those used in the WDC-HE model. We find that the WDC-HE model accurately predicts the quasiequilibrium ion current losses in the EM PIC simulations for a wide range of current and voltage values. For the case of two ion species where one is magnetically insulated by the ambient magnetic field and the other is not, the charge of the lighter insulated species in the anode-cathode gap can modify the electric field profile, reducing the ion current density enhancement for the heavier ion species. On the other hand, for multiple ion species, when the lighter ions are not magnetically insulated and are a significant fraction of the anode plasma, they dominate the current loss, producing loss currents which are a significant fraction of the lighter ion WDC values. The observation of this effect in the present work is new to the field and may significantly impact the analysis of ion current losses in the Z machine inner MITL and convolute.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Current loss in magnetically insulated transmission lines (MITLs) was investigated using data from experiments conducted on Z and Mykonos. Data from experiments conducted on Z were used to optimize an ion diode current loss model that has been implemented into the transmission line circuit model of Z. Details on the current loss model and comparisons to data from Z experiments have been previously published in a peer-reviewed journal [Hutsel, et al., Phys. Rev. Accel. Beams 21, 030401]. Dedicated power flow experiments conducted on Mykonos investigated current loss in a millimeter-scale anode-cathode gap MITL operated at lineal current densities greater than 410 kA/cm and with electric field stresses in excess of 240 kV/cm where it is expected that both anode and cathode plasmas are formed. The experiment MITLs were exposed to varying vacuum conditions; including vacuum pressure at shot time, time under vacuum, and vacuum storage protocols. The results indicate that the vacuum conditions have an effect on current loss in high lineal current density MITLs.

Abstract not provided.

Abstract not provided.

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%.

Abstract not provided.

Abstract not provided.

A series of simulations and experiments to resolve questions about the operation of arrays of closely spaced small aspect ratio rod pinches has been performed. Design and postshot analysis of the experimental results are supported by 3-D particle-in-cell simulations. Both simulations and experiments support these conclusions. Penetration of current to the interior of the array appears to be efficient, as the current on the center rods is essentially equal to the current on the outer rods. Current loss in the feed due to the formation of magnetic nulls was avoided in these experiments by design of the feed surface of the cathode and control of the gap to keep the electric fields on the cathode below the emission threshold. Some asymmetry in the electron flow to the rod was observed, but the flow appeared to symmetrize as it reached the end of the rod. Interaction between the rod pinches can be controlled to allow the stable and consistent operation of arrays of rod pinches.

Abstract not provided.

Abstract not provided.

Abstract not provided.