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.
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%.
Recent Magnetized Liner Inertial Fusion experiments at the Sandia National Laboratories Z pulsed power facility have featured a PDV (Photonic Doppler Velocimetry) diagnostic in the final power feed section for measuring load current. In this paper, we report on an anomalous pressure that is detected on this PDV diagnostic very early in time during the current ramp. Early time load currents that are greater than both B-dot upstream current measurements and existing Z machine circuit models by at least 1 MA would be necessary to describe the measured early time velocity of the PDV flyer. This leads us to infer that the pressure producing the early time PDV flyer motion cannot be attributed to the magnetic pressure of the load current but rather to an anomalous pressure. Using the MHD code ALEGRA, we are able to compute a time-dependent anomalous pressure function, which when added to the magnetic pressure of the load current, yields simulated flyer velocities that are in excellent agreement with the PDV measurement. As a result, we also provide plausible explanations for what could be the origin of the anomalous pressure.
Recent Magnetized Liner Inertial Fusion experiments at the Sandia National Laboratories Z pulsed power facility have featured a PDV (Photonic Doppler Velocimetry) diagnostic in the final power feed section for measuring load current. In this paper, we report on an anomalous pressure that is detected on this PDV diagnostic very early in time during the current ramp. Early time load currents that are greater than both B-dot upstream current measurements and existing Z machine circuit models by at least 1 MA would be necessary to describe the measured early time velocity of the PDV flyer. This leads us to infer that the pressure producing the early time PDV flyer motion cannot be attributed to the magnetic pressure of the load current but rather to an anomalous pressure. Using the MHD code ALEGRA, we are able to compute a time-dependent anomalous pressure function, which when added to the magnetic pressure of the load current, yields simulated flyer velocities that are in excellent agreement with the PDV measurement. We also provide plausible explanations for what could be the origin of the anomalous pressure.
Several magnetized liner inertial fusion (MagLIF) experiments have been conducted on the Z accelerator at Sandia National Laboratories since late 2013. Measurements of the primary DD (2.45 MeV) neutrons for these experiments suggest that the neutron production is thermonuclear. Primary DD yields up to 3e12 with ion temperatures ∼2-3 keV have been achieved. Measurements of the secondary DT (14 MeV) neutrons indicate that the fuel is significantly magnetized. Measurements of down-scattered neutrons from the beryllium liner suggest ρRliner∼1g/cm2. Neutron bang times, estimated from neutron time-of-flight (nTOF) measurements, coincide with peak x-ray production. Plans to improve and expand the Z neutron diagnostic suite include neutron burn-history diagnostics, increased sensitivity and higher precision nTOF detectors, and neutron recoil-based yield and spectral measurements.
We report on the progress made to date for a Laboratory Directed Research and Development (LDRD) project aimed at diagnosing magnetic flux compression on the Z pulsed-power accelerator (0-20 MA in 100 ns). Each experiment consisted of an initially solid Be or Al liner (cylindrical tube), which was imploded using the Z accelerator's drive current (0-20 MA in 100 ns). The imploding liner compresses a 10-T axial seed field, B z ( 0 ) , supplied by an independently driven Helmholtz coil pair. Assuming perfect flux conservation, the axial field amplification should be well described by B z ( t ) = B z ( 0 ) x [ R ( 0 ) / R ( t )] 2 , where R is the liner's inner surface radius. With perfect flux conservation, B z ( t ) and dB z / dt values exceeding 10 4 T and 10 12 T/s, respectively, are expected. These large values, the diminishing liner volume, and the harsh environment on Z, make it particularly challenging to measure these fields. We report on our latest efforts to do so using three primary techniques: (1) micro B-dot probes to measure the fringe fields associated with flux compression, (2) streaked visible Zeeman absorption spectroscopy, and (3) fiber-based Faraday rotation. We also mention two new techniques that make use of the neutron diagnostics suite on Z. These techniques were not developed under this LDRD, but they could influence how we prioritize our efforts to diagnose magnetic flux compression on Z in the future. The first technique is based on the yield ratio of secondary DT to primary DD reactions. The second technique makes use of the secondary DT neutron time-of-flight energy spectra. Both of these techniques have been used successfully to infer the degree of magnetization at stagnation in fully integrated Magnetized Liner Inertial Fusion (MagLIF) experiments on Z [P. F. Schmit et al. , Phys. Rev. Lett. 113 , 155004 (2014); P. F. Knapp et al. , Phys. Plasmas, 22 , 056312 (2015)]. Finally, we present some recent developments for designing and fabricating novel micro B-dot probes to measure B z ( t ) inside of an imploding liner. In one approach, the micro B-dot loops were fabricated on a printed circuit board (PCB). The PCB was then soldered to off-the-shelf 0.020- inch-diameter semi-rigid coaxial cables, which were terminated with standard SMA connectors. These probes were recently tested using the COBRA pulsed power generator (0-1 MA in 100 ns) at Cornell University. In another approach, we are planning to use new multi-material 3D printing capabilities to fabricate novel micro B-dot packages. In the near future, we plan to 3D print these probes and then test them on the COBRA generator. With successful operation demonstrated at 1-MA, we will then make plans to use these probes on a 20-MA Z experiment.
Velocity Interferometer System for Any Reflector (VISAR) [Barker and Hollenbach, J. Appl. Phys. 43, 4669 (1972)] is a well-known diagnostic that is employed on many shock physics and pulsed-power experiments. With the VISAR diagnostic, the velocity on the surface of any metal flyer can be found. For most experiments employing VISAR, either a kinetic pressure [Grady, Mech. Mater. 29, 181 (1998)] or a magnetic pressure [Lemke et al., Intl J. Impact Eng. 38, 480 (2011)] drives the motion of the flyer. Moreover, reliable prediction of the time-dependent pressure is often a critical component to understanding the physics of these experiments. Although VISAR can provide a precise measurement of a flyer’s surface velocity, the real challenge of this diagnostic implementation is using this velocity to unfold the time-dependent pressure. The purpose of this study is to elucidate a new method for quickly and reliably unfolding VISAR data.
Magnetized Liner Inertial Fusion experiments performed at Sandia's Z facility have demonstrated significant thermonuclear fusion neutron yields (∼1012 DD neutrons) from multi-keV deuterium plasmas inertially confined by slow (∼10 cm/μs), stable, cylindrical implosions. Effective magnetic confinement of charged fusion reactants and products is signaled by high secondary DT neutron yields above 1010. Analysis of extensive power, imaging, and spectroscopic x-ray measurements provides a detailed picture of ∼3 keV temperatures, 0.3 g/cm3 densities, gradients, and mix in the fuel and liner over the 1-2 ns stagnation duration.
In this study, the magnetized liner inertial fusion concept [S. A. Slutz et al., Phys. Plasmas17, 056303 (2010)] utilizes a magnetic field and laser heating to relax the pressure requirements of inertial confinement fusion. The first experiments to test the concept [M. R. Gomez et al., Phys. Rev. Lett. 113, 155003 (2014)] were conducted utilizing the 19 MA, 100 ns Z machine, the 2.5 kJ, 1 TW Z Beamlet laser, and the 10 T Applied B-field on Z system. Despite an estimated implosion velocity of only 70 km/s in these experiments, electron and ion temperatures at stagnation were as high as 3 keV, and thermonuclear deuterium-deuterium neutron yields up to 2 × 1012 have been produced. X-ray emission from the fuel at stagnation had widths ranging from 50 to 110 μm over a roughly 80% of the axial extent of the target (6–8 mm) and lasted approximately 2 ns. X-ray yields from these experiments are consistent with a stagnation density of the hot fuel equal to 0.2–0.4 g/cm3. In these experiments, up to 5 ×1010 secondary deuterium-tritium neutrons were produced. Given that the areal density of the plasma was approximately 1–2 mg/cm2, this indicates the stagnation plasma was significantly magnetized, which is consistent with the anisotropy observed in the deuterium-tritium neutron spectra. Control experiments where the laser and/or magnetic field were not utilized failed to produce stagnation temperatures greater than 1 keV and primary deuterium-deuterium yields greater than 1010. An additional control experiment where the fuel contained a sufficient dopant fraction to substantially increase radiative losses also failed to produce a relevant stagnation temperature. The results of these experiments are consistent with a thermonuclear neutron source.
The standard approaches to inertial confinement fusion (ICF) rely on implosion velocities greater than 300 km/s and spherical convergence to achieve the high fuel temperatures (T > 4 keV) and areal densities (ρr > 0.3 g/cm2) required for ignition1. Such high velocities are achieved by heating the outside surface of a spherical capsuleeither directly with a large number of laser beams (Direct Drive) or with x-rays generated within a hohlraum (Indirect Drive). A much more energetically efficient approach is to use the magnetic pressure generated by a pulsed power machine to directly drive an implosion. In this approach 5-10% of the stored energy can be converted to the implosion of a metal tube generally referred to as a “liner”. However, the implosion velocity is not very high 70-100 km/s and the convergence is cylindrical (rather than spherical) making it more difficult to achieve the high temperatures and areal densities needed for ignition.