Neutron bang times for a series of MagLIF (Magnetic Liner Inertial Fusion) experiments with D2-filled targets have been measured at the Z facility. The emitted neutrons were detected as current-mode pulses in a multichannel, neutron time-of-flight (nTOF) diagnostic with conventional, scintillator-photomultiplier-tube (PMT) detectors. In these experiments, the detectors were fielded at known, fixed distances L (690-2510 cm) from the target, and on three, non-coplanar (but convergent) lines-of-sight (LOS). The primary goal of this diagnostic was to estimate a fiducial time (bang time) relative to an externally generated time-base for synchronizing all the diagnostics in an experiment. Recorded arrival times (A7) of the pulses were characterized experimentally by three numerical methods: a first-moment estimate (centroid) and two nodal measures — Savitzky-Golay (SG) smoothing and a single point peak estimate of the raw data. These times were corrected for internal detector time delays (transit and impulse-response function) — an adjustment that linked the recorded ATs to the corresponding arrival of uncollided neutrons at each detector. The bang time was then estimated by linearly regressing the arrival times against the associated distances to the source; tbang (on the system timescale) was taken as the temporal intercept of the regression equation at distance L = 0. This article reports the analysis for a representative shot #2584 for which (a) the recorded ATs — even without detector corrections — agreed by method in each channel to within 1-2 ns; (b) internal corrections were each ~3 — 5 ns; and (c) a 95% uncertainty (confidence) interval for tbang in this shot was estimated at ±3 ns with 4 degrees of freedom. A secondary goal for this diagnostic was to check that the bang time measurements corresponded to neutrons emitted by the D(d,n)3He reaction in a thermalized DD plasma. According to the theoretical studies by Brysk, such neutrons should be emitted with an isotropic Gaussian distribution of mean kinetic energy $ \overline{E}$ of 2.449 MeV; this energy translates to a mean neutron speed $ \overline{u}$ of 2.160 cm/ns [D. H. Munro, Nuclear Fusion, 56(3) 036001 (2016)]. In the MagLIF series of shots there was no evidence of spatial asymmetry in the time-distance regressions, and it was possible to extract the mean neutron speed from the slope of these fits. In shot 2584 $ \overline{u}$ was estimated at 2.152 cm/ns ± 0.010 cm/ns [95 % confidence, 4 dof] and the mean kinetic energy $ \overline{E}$ (with relativistic corrections) was 2.431 MeV ± 0.022 MeV [95 % confidence, 4 dof] — results supporting the assumption that D-D neutrons were, in fact, measured.
We recently developed a one-dimensional imager of neutrons on the Z facility. The instrument is designed for Magnetized Liner Inertial Fusion (MagLIF) experiments, which produce D-D neutrons yields of ∼3 × 1012. X-ray imaging indicates that the MagLIF stagnation region is a 10-mm long, ∼100-μm diameter column. The small radial extents and present yields precluded useful radial resolution, so a one-dimensional imager was developed. The imaging component is a 100-mm thick tungsten slit; a rolled-edge slit limits variations in the acceptance angle along the source. CR39 was chosen as a detector due to its negligible sensitivity to the bright x-ray environment in Z. A layer of high density poly-ethylene is used to enhance the sensitivity of CR39. We present data from fielding the instrument on Z, demonstrating reliable imaging and track densities consistent with diagnosed yields. For yields ∼3 × 1012, we obtain resolutions of ∼500 μm.
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.
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.
We present results from crystal spectroscopic analysis of silicon aero-gel foams heated by dynamic hohlraums on Z. The dynamic hohlraum on Z creates a radiation source with a 230-eV average temperature over a 2.4-mm diameter. In these experiments silicon aero-gel foams with 10-mg/cm{sup 3} densities and 1.7-mm lengths were placed on both ends of the dynamic hohlraum. Several crystal spectrometers were placed both above and below the z-pinch to diagnose the temperature of the silicon aero-gel foam using the K-shell lines of silicon. The crystal spectrometers were (1) temporally integrated and spatially resolved, (2) temporally resolved and spatially integrated, and (3) both temporally and spatially resolved. The results indicate that the dynamic hohlraum heats the silicon aero-gel to approximately 150-eV at peak power. As the dynamic hohlraum source cools after peak power the silicon aero-gel continues to heat and jets axially at an average velocity of approximately 50-cm/{micro}s. The spectroscopy has also shown that the reason for the up/down asymmetry in radiated power on Z is that tungsten enters the line-of-sight on the bottom of the machine much more than on the top.