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%.
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.
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.
The magnetized liner inertial fusion concept [S. A. Slutz et al., Phys. Plasmas 17, 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.
Using the ideal magnetohydrodynamic model, we calculate the temporal evolution of initial ripples on the boundaries of a planar plasma slab that is subjected to the magneto-Rayleigh-Taylor instability. The plasma slab consists of three regions. We assume that in each region the plasma density is constant with an arbitrary value and the magnetic field is also constant with an arbitrary magnitude and an arbitrary direction parallel to the interfaces. Thus, the instability may be driven by a combination of magnetic pressure and kinetic pressure. The general dispersion relation is derived, together with the feedthrough factor between the two interfaces. The temporal evolution is constructed from the superposition of the eigenmodes. Previously established results are recovered in the various limits. Numerical examples are given on the temporal evolution of ripples on the interfaces of the finite plasma slab.
This Letter presents results from the first fully integrated experiments testing the magnetized liner inertial fusion concept [S.A. Slutz et al., Phys. Plasmas 17, 056303 (2010)], in which a cylinder of deuterium gas with a preimposed axial magnetic field of 10 T is heated by Z beamlet, a 2.5 kJ, 1 TW laser, and magnetically imploded by a 19 MA current with 100 ns rise time on the Z facility. Despite a predicted peak implosion velocity of only 70 km/s, the fuel reaches a stagnation temperature of approximately 3 keV, with Te ≈ Ti, and produces up to 2e12 thermonuclear DD neutrons. In this study, X-ray emission indicates a hot fuel region with full width at half maximum ranging from 60 to 120 μm over a 6 mm height and lasting approximately 2 ns. The number of secondary deuterium-tritium neutrons observed was greater than 1010, indicating significant fuel magnetization given that the estimated radial areal density of the plasma is only 2 mg/cm2.