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.
To support the nuclear fusion program at Sandia National Laboratories (SNL), a consistent and verifiable method to determine fusion ion temperatures needs to be developed. Since the fusion temperature directly affects the width in the spread of neutron energies produced, a measurement of the neutron energy width can yield the fusion temperature. Traditionally, the spread in neutron energies is measured by using time-of-flight to convert a spread in neutron energies at the source to a spread in time at detector. One potential obstacle to using this technique at the Z facility at SNL is the need to shield the neutron detectors from the intense bremsstrahlung produced. The shielding consists of eight inches of lead and the concern is that neutrons will scatter in the lead, artificially broaden the neutron pulse width and lead to an erroneous measurement. To address this issue, experiments were performed at the University of Rochester's Laboratory for Laser Energetics, which demonstrated that a reliable ion temperature measurement can be achieved behind eight inches of lead shielding. To further expand upon this finding, Monte Carlo N-Particle eXtended (MCNPX) was used to simulate the experimental geometric conditions and perform the neutron transport. MCNPX was able to confidently estimate results observed at the University of Rochester.
The conclusions of this report are: (1) 1D and 2D RMHD simulations indicate feasibility of producing high thermonuclear neutron yields in deuterium and DT gas-puff Z-pinches -- (a) Z 1.7 x 10{sup 13} DD neutrons at 70 kV, 13 MA (Z1384); (b) (3 to 6) x 10{sup 14} at 90 kV, 17 MA (Z1422); (c) Predicted for ZR 2 x 10{sup 15} DD and 6 x 10{sup 16} DT neutrons; (2) Theory and modeling issues -- collisionless ions, nonthermal ions; (3) Experimental data on the origin of the neutrons not yet conclusive, need more shots; and (4) Applications -- (a) Fusion 2.5 and 14 MeV neutron source; (b) Pulsed subcritical neutron source with uranium blanket for {approx}10x neutron and {approx}1000x energy multiplication (Smirnov, Feoktistov and Klimov); and (c) Fusion-assisted keV x-ray plasma radiation source.
Progress in understanding the physics of dynamic-hohlraums is reviewed for a system capable of generating 13 TW of axial radiation for high temperature (>200 eV) radiation-flow experiments and ICF capsule implosions.
The lead probe neutron detector was originally designed by Spencer and Jacobs in 1965. The detector is based on lead activation due to the following neutron scattering reactions: {sup 207}Pb(n, n'){sup 207m}Pb and {sup 208}Pb(n, 2n){sup 207m}Pb. Delayed gammas from the metastable state of {sup 207m}Pb are counted using a plastic scintillator. The half-life of {sup 207m}Pb is 0.8 seconds. In the work reported here, MCNP was used to optimize the efficiency of the lead probe by suitably modifying the original geometry. A prototype detector was then built and tested. A 'layer cake' design was investigated in which thin (< 5 mm) layers of lead were sandwiched between thicker ({approx} 1 - 2 cm) layers of scintillator. An optimized 'layer cake' design had Figures of Merit (derived from the code) which were a factor of 3 greater than the original lead probe for DD neutrons, and a factor of 4 greater for DT neutrons, while containing 30% less lead. A smaller scale, 'proof of principle' prototype was built by Bechtel/Nevada to verify the code results. Its response to DD neutrons was measured using the DD dense plasma focus at Texas A&M and it conformed to the predicted performance. A voltage and discriminator sweep was performed to determine optimum sensitivity settings. It was determined that a calibration operating point could be obtained using a {sup 133}Ba 'bolt' as is the case with the original lead probe.
Inertial confinement fusion capsule implosions absorbing up to 35 kJ of x-rays from a {approx}220 eV dynamic hohlraum on the Z accelerator at Sandia National Laboratories have produced thermonuclear D-D neutron yields of (2.6 {+-} 1.3) x 10{sup 10}. Argon spectra confirm a hot fuel with Te {approx} 1 keV and n{sub e} {approx} (1-2) x 10{sup 23} cm{sup -3}. Higher performance implosions will require radiation symmetry control improvements. Capsule implosions in a {approx}70 eV double-Z-pinch-driven secondary hohlraum have been radiographed by 6.7 keV x-rays produced by the Z-beamlet laser (ZBL), demonstrating a drive symmetry of about 3% and control of P{sub 2} radiation asymmetries to {+-}2%. Hemispherical capsule implosions have also been radiographed in Z in preparation for future experiments in fast ignition physics. Z-pinch-driven inertial fusion energy concepts are being developed. The refurbished Z machine (ZR) will begin providing scaling information on capsule and Z-pinch in 2006. The addition of a short pulse capability to ZBL will enable research into fast ignition physics in the combination of ZR and ZBL-petawatt. ZR could provide a test bed to study NIF-relevant double-shell ignition concepts using dynamic hohlraums and advanced symmetry control techniques in the double-pinch hohlraum backlit by ZBL.
Progress in Z-pinch experiments at Sandia's Saturn facility have underscored a need for an absolute yield measurement for DD fusion neutrons. The technique chosen for making this absolute yield measurement was neutron activation of indium metal samples. To calibrate the technique, a 175-keV deuteron beam was allowed to impinge on a 3.0-μm-thick erbium deuteride target, producing neutrons through the 2H(d,n)3He fusion reaction. The neutron flux produced at 0° and incident on nominal 5-g indium samples was determined by the associated particle method. This method employed protons measured from the 2H(d,p)3H reaction to infer the neutron flux produced. After neutron irradiation, the activity of the indium samples was measured with a Ge gamma-ray detector. The total activity of the metastable state 115mIn (336.23 keV) was measured, compared with the total incident flux, and a calibration factor (indium counts/neutron/gram of indium) determined. For completeness, a calibration factor for DT neutrons from the 3H(d,n)4He fusion reaction was also obtained through the measured activity of the metastable state 114mIn(190.29 keV). The experiment and the measured calibration factors for both reactions are described in the paper.
