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.
By magnetizing the fusion fuel in inertial confinement fusion (ICF) systems, the required stagnation pressure and density can be relaxed dramatically. This happens because the magnetic field insulates the hot fuel from the cold pusher and traps the charged fusion burn products. This trapping allows the burn products to deposit their energy in the fuel, facilitating plasma self-heating. Here, we report on a comprehensive theory of this trapping in a cylindrical DD plasma magnetized with a purely axial magnetic field. Using this theory, we are able to show that the secondary fusion reactions can be used to infer the magnetic field-radius product, BR, during fusion burn. This parameter, not ρR, is the primary confinement parameter in magnetized ICF. Using this method, we analyze data from recent Magnetized Liner Inertial Fusion experiments conducted on the Z machine at Sandia National Laboratories. We show that in these experiments BR ≈ 0.34(+0.14/-0.06) MG cm, a ∼ 14x increase in BR from the initial value, and confirming that the DD-fusion tritons are magnetized at stagnation. This is the first experimental verification of charged burn product magnetization facilitated by compression of an initial seed magnetic flux.
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.
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.
We present a general methodology to determine the diagnostic sensitivity that is directly applicable to neutron-activation diagnostics fielded on a wide variety of neutron-producing experiments, which include inertial-confinement fusion (ICF), dense plasma focus, and ion beam-driven concepts. This approach includes a combination of several effects: (1) non-isotropic neutron emission; (2) the 1/r2 decrease in neutron fluence in the activation material; (3) the spatially distributed neutron scattering, attenuation, and energy losses due to the fielding environment and activation material itself; and (4) temporally varying neutron emission. As an example, we describe the copper-activation diagnostic used to measure secondary deuterium-tritium fusion-neutron yields on ICF experiments conducted on the pulsed-power Z Accelerator at Sandia National Laboratories. Using this methodology along with results from absolute calibrations and Monte Carlo simulations, we find that for the diagnostic configuration on Z, the diagnostic sensitivity is 0.037% ± 17% counts/neutron per cm2 and is ~ 40% less sensitive than it would be in an ideal geometry due to neutron attenuation, scattering, and energy-loss effects.
There are several machines in this country that produce short bursts of neutrons for various applications. A few examples are the Zmachine, operated by Sandia National Laboratories in Albuquerque, NM; the OMEGA Laser Facility at the University of Rochester in Rochester, NY; and the National Ignition Facility (NIF) operated by the Department of Energy at Lawrence Livermore National Laboratory in Livermore, California. They all incorporate neutron time of flight (nTOF) detectors which measure neutron yield, and the shapes of the waveforms from these detectors contain germane information about the plasma conditions that produce the neutrons. However, the signals can also be %E2%80%9Cclouded%E2%80%9D by a certain fraction of neutrons that scatter off structural components and also arrive at the detectors, thereby making analysis of the plasma conditions more difficult. These detectors operate in current mode - i.e., they have no discrimination, and all the photomultiplier anode charges are integrated rather than counted individually as they are in single event counting. Up to now, there has not been a method for modeling an nTOF detector operating in current mode. MCNPPoliMiwas developed in 2002 to simulate neutron and gammaray detection in a plastic scintillator, which produces a collision data output table about each neutron and photon interaction occurring within the scintillator; however, the postprocessing code which accompanies MCNPPoliMi assumes a detector operating in singleevent counting mode and not current mode. Therefore, the idea for this work had been born: could a new postprocessing code be written to simulate an nTOF detector operating in current mode? And if so, could this process be used to address such issues as the impact of neutron scattering on the primary signal? Also, could it possibly even identify sources of scattering (i.e., structural materials) that could be removed or modified to produce %E2%80%9Ccleaner%E2%80%9D neutron signals? This process was first developed and then applied to the axial neutron time of flight detectors at the ZFacility mentioned above. First, MCNPPoliMi was used to model relevant portions of the facility between the source and the detector locations. To obtain useful statistics, variance reduction was utilized. Then, the resulting collision output table produced by MCNPPoliMi was further analyzed by a MATLAB postprocessing code. This converted the energy deposited by neutron and photon interactions in the plastic scintillator (i.e., nTOF detector) into light output, in units of MeVee%D1%84 (electron equivalent) vs time. The time response of the detector was then folded into the signal via another MATLAB code. The simulated response was then compared with experimental data and shown to be in good agreement. To address the issue of neutron scattering, an %E2%80%9CIdeal Case,%E2%80%9D (i.e., a plastic scintillator was placed at the same distance from the source for each detector location) with no structural components in the problem. This was done to produce as %E2%80%9Cpure%E2%80%9D a neutron signal as possible. The simulated waveform from this %E2%80%9CIdeal Case%E2%80%9D was then compared with the simulated data from the %E2%80%9CFull Scale%E2%80%9D geometry (i.e., the detector at the same location, but with all the structural materials now included). The %E2%80%9CIdeal Case%E2%80%9D was subtracted from the %E2%80%9CFull Scale%E2%80%9D geometry case, and this was determined to be the contribution due to scattering. The time response was deconvolved out of the empirical data, and the contribution due to scattering was then subtracted out of it. A transformation was then made from dN/dt to dN/dE to obtain neutron spectra at two different detector locations.