SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Nuclear fission is an important process with applications in astrophysics, nuclear reactors, and stockpile stewardship. Despite having been discovered over 80 years ago, fission is not fully described by a predictive model.
Precise knowledge of nuclear fission is a critical component for nuclear energy and stockpile stewardship research. Despite having been discovered over 80 years ago, fission is not fully described by a predictive theory, so thorough and accurate nuclear data libraries are necessary. These libraries provide information about the fission process in studied isotopes that is used in scientific and engineering models.
We explore, using a state-of-the-art simulation code in 3D and late-enough times to witness final observables, the dependence of core-collapse supernova explosions on the nuclear equation of state (EOS). Going beyond questions of explodability, we compare final explosion energies, nucleosynthetic yields, recoil kicks, and gravitational-wave and neutrino signatures using the SFHo and DD2 nuclear EOSs for a 9 M ⊙ /solar-metallicity progenitor star. The DD2 EOS is stiffer and has a lower effective nucleon mass. The result is a more extended protoneutron star (PNS) and lower central densities. As a consequence, the mean neutrino energies, final explosion energy, and recoil kick speed are lower. Moreover, the evolution of PNS convection differs between the two EOS models in significant ways. This translates in part into interestingly altered neutrino “light” curves and noticeably altered gravitational-wave signal strengths and frequency characteristics that may be diagnostic. The faster exploding model (SFHo) yields slightly more neutron-rich ejecta and more species with atomic weights between 60 and 90 and a weak r-process. However, this is merely a preliminary study. The next step is a more comprehensive and multiprogenitor set of 3D supernova simulations for various EOSs to late times when the observables have asymptoted. Such a future investigation will have a direct bearing on the neutron star and black hole birth mass functions and the quest toward a fully quantitative theory of supernova observables.
The goal of this project was to produce modernized nuclear data libraries of light ion fusion reactions and actinide fission product yields. Such libraries are essential input in both stockpile science but also fundamental science (nucleosynthesis, superheavy nuclei). When such libraries are used to model complex processes involving exotic, very short-lived nuclei that have not been measured experimentally, the underlying nuclear data come from theoretical models, the predictive power of which needs to be quantified, together with its uncertainties. By teaming up nuclear theorists and statisticians, we built comprehensive theoretical frameworks based on the most fundamental theories of fusion and fission that we coupled for the first time with a Bayesian framework with machine learning techniques for uncertainty quantification. Our results include: (1) statistical emulators for neutron-alpha scattering, and the fission product mass yields of neutron-induced fission of 239 Pu, which will enable developers to generate random realizations of the corresponding data for uncertainty quantification (UQ) in applications; (2) the first complete and most accurate ab initio prediction of the S-factor for proton capture on 7 Be, which – combined with UQ – may become the recommended value at solar energies in the next review of solar fusion cross sections; (3) statistical emulators for fission mass yields in the neutron-induced fission of 239 Pu; (4) the first quantification of how uncertainties in nuclear forces impact fission properties that are relevant in nucleosynthesis simulations to explain how heavy elements are formed in the Cosmos.
Albergo thermometers with double isotope, isotone and isobar yield ratio pairs with one proton or/and neutron difference are investigated. Without any extra sequential decay correction, a real temperature value of 4.9 ± 0.5 MeV is deduced from the yields of the experimentally reconstructed primary hot intermediate mass fragments (IMFs) from 64 Zn + 112 Sn collisions at 40 MeV/nucleon using the Albergo thermometer for the first time. An experimental sequential decay correction from the apparent temperatures to the real ones for twelve other reaction systems with different neutron-proton (N/Z) asymmetries in the same experiment, 70 Zn, 64 Ni on 112,124 Sn, 58,64 Ni, 197 Au, 232 Th at 40 MeV/nucleon, is performed using an empirical correction factor approach of Tsang et al. with the deduced 4.9 MeV temperature value. The dependence of nuclear temperature on the source N/Z asymmetry is further investigated using these deduced real source temperature values from the present thirteen systems. It is found that the deduced real source temperatures at the present source N/Z range show a rather weak dependence on the source N/Z asymmetry. In conclusion, by comparison between our previous results and those from other independent experiments, a consistent description for the N/Z asymmetry dependence of nuclear temperature is addressed.
Pebble bed reactors (PBRs) can improve the safety and economics of the nuclear energy production. PBRs rely on TRIstructural-ISOtropic (TRISO) fuel pebbles for enhanced fission product retention. Accurate characterization of individual fuel pebbles would enable the validation of computational models, efficient use of TRISO fuel, and improve fuel accountability. Here, we have developed and tested a new neutron multiplicity counter (NMC) based on 192 boron coated straw (BCS) detectors optimized for 235 U assay in TRISO fuel. The new design yielded a singles and doubles neutron detection efficiency of 4.71% and 0.174%, respectively, and a die-away time of 16.7 μs. The NMC has a low intrinsic gamma-ray detection efficiency of 8.71 x 10 –8 at an exposure rate of 80.3 mR/h. In simulation, a high-efficiency version of the NMC encompassing 396 straws was able to estimate the 235 U in a pebble with a relative uncertainty and error both below 2% in 100 s.
The 2-D position sensitive neutron detectors featuring a significantly improved count rate capacity and a sub-mm spatial resolution are in high demand by neutron scattering facilities, especially the next generation high-flux spallation neutron sources (SNSs). They are anticipated to play a crucial role in fully utilizing the unprecedentedly high neutron brightness offered by these sources. In this work, a high rate and high spatial resolution neutron Anger camera has been successfully developed based on a Ce-doped microcolumnar LiI scintillator film and the PETsys TOFPET2 ASIC readout electronics. The performance evaluation was conducted in a standard laboratory environment as well as the High Flux Isotope Reactor (HFIR) at Oak Ridge National Laboratory (ORNL) with regard to light yield, γ-ray sensitivity, neutron detection efficiency, count rate capability, spatial resolution, and position linearity. The LiI:Ce scintillator was found to have a higher light yield and a lower γ-ray sensitivity than GS20 (i.e., 11729 versus 6000 photons/n and 1.67 × 10 −5 versus on the order of 10 −4 ). The neutron detection efficiency was measured to be 63% at 4.2 Å. This camera demonstrated an exceptional count rate capability, up to at least 481 kHz without rate loss. The spatial resolution was as good as 0.47 mm. Here, the position linearity was also enhanced.
Lithium-loaded organic plastic scintillators combine sensitivity to γ rays with the ability to detect both fast and slow neutrons, making them valuable for applications in nuclear security and in basic nuclear and particle physics. The goal of this work is to characterize the neutron response of two custom lithium-loaded organic plastic scintillators developed at Lawrence Livermore National Laboratory. Both are ternary polystyrene-based formulations containing 1.5 wt.% 6 Li salts of isobutyric acid, but they differ in their primary and secondary dye compositions: one uses m-terphenyl as the primary fluor and with Exalite 404 as the wavelength shifter, whereas the other uses 2,5-diphenyloxazole (PPO) and 9,10-diphenyl-anthracene, respectively. The temporal response of the scintillators was measured via time-correlated single photon counting for γ-ray and neutron events. The proton light yield was measured using the double time-of-flight technique from 1.3 to 15 MeV at the 88-Inch Cyclotron at Lawrence Berkeley National Laboratory. For the slow neutron response, an AmBe source moderated with polyethylene was used, and the light output from the 6 Li(n,α)t reaction was characterized. Differences in ionization quenching and temporal response were observed between the two materials with the PPO-containing scintillator exhibiting higher ionization quenching. These results provide performance benchmarks that can guide the design and optimization of future lithium-loaded plastic scintillators for use in basic science and applications.
The most commonly used screens for neutron imaging consist of 6 LiF + ZnS. This type of screen yields the highest light output per detected neutron. For high resolution, gadolinium oxysulfide (GOS, Gadox) screens are employed, which have a much higher detection efficiency, but a light output so much lower than LiF + ZnS that measurements are often limited by photon statistics. Historically, screens using boron as a neutron-sensitive material have not been very successful. However, a new preparation method was introduced recently that produces light output higher than Gadox with detection efficiency greater than LiF + ZnS. Measurements of these new borated screens were performed at the NeXT facility at ILL, Grenoble, in comparison to a high resolution Gadox screen.
Born in the aftermath of core-collapse supernovae, neutron stars contain matter under extraordinary conditions of density and temperature that are difficult to reproduce in the laboratory. In recent years, neutron star observations have begun to yield novel insights into the nature of strongly interacting matter in the high-density regime where current theoretical models are challenged. At the same time, chiral effective field theory has developed into a powerful framework to study nuclear matter properties with quantified uncertainties in the moderate-density regime for modeling neutron stars. In this article, we review recent developments in chiral effective field theory and focus on many-body perturbation theory as a computationally efficient tool for calculating the properties of hot and dense nuclear matter. We also demonstrate how effective field theory enables statistically meaningful comparisons among nuclear theory predictions, nuclear experiments, and observational constraints on the nuclear equation of state.
The deuterium–tritium (DT) fusion reaction, D ( T , γ ) 5 He , has a rare (4.2 × 10 – 5 ) branching ratio that releases a ~16 MeV gamma ray instead of the more common 14.1 MeV fusion neutron. These fusion gamma rays can be used as a complementary source of measuring the net fusion yield of experiments. Fusion gamma rays have the advantage of having negligible Doppler broadening, being isotropic and less down-scattered compared to fusion neutrons. At the National Ignition Facility, we measured the DT fusion gamma rays by the Gamma Reaction History diagnostic that thresholds to only measure >10 MeV gamma rays. After removing a ~12% contribution from neutron capture interactions in deuterium, D ( n , γ ) T , and carbon, 12 C ( n , γ ) 13 C , the deviation between GRH measurements and neutron-based yield measurements has been reduced to 10%. Further improvements are needed for detailed physics comparison.
This presentation discusses the experimental, simulation, and nuclear data methods that were validated for the RPI γ-Multiplicity Detector. When the neutron capture γ-cascade data is well-known, the γ-emission spectra can be accurately calculated using the modified simulation tools. The RPI γ-Multiplicity Detector system is now ready for analysis and recommendations for isotopes with deficiencies in γ-ray data. The presentation also discusses future work which includes developing a method for analyzing and adjusting nuclear data for 59 Co, 55 Mn and other measured isotopes including 181 Ta. Additionally, future work includes comparing experimental γ-emission spectra with MCNP-6.2/DICEBOX simulations for 238 U and 235 U. In summation, new capture and transmission measurements for 54 Fe will help improve resonance parameter evaluation. Neutron capture gamma cascade spectra and yields were measured in the resolved resonance region and compared to evaluations. In addition, the pulsed neutron die-away method was developed as a tool to provide data for validation of TSLs.
To resolve discrepancies in evaluated cross sections among major nuclear data libraries, energy-differential neutron transmission and radiative capture yield of 181 Ta were measured from 0.15 to 100 keV using multiple sample thicknesses. The new measurements provide resolution such that the resolved resonance region (RRR) can be evaluated up to at least 2.5 keV and the unresolved resonance region can be evaluated up to at least 100 keV. The transmission and capture yield measurements were modeled using resonance parameters from three major libraries to assess the predictive capability of each. It was found that JENDL-5.0 performed best in the RRR. Because of the poor performance of the U.S. ENDF/B evaluation, it is recommended that ENDF/B be reevaluated for 181 Ta.
Measurements of physical quantities are often made relative to a reference quantity, simplifying the interpretation of the measurement and its uncertainties. However, if applied without careful considerations, this technique can carry with it a significant and complex systematic error that is usually ignored. For example, measurements of continuous neutron spectra relative to a reference spectrum yield a result that does not properly include the important effects of neutron scattering in the experimental environment. These effects act to distort the measured spectrum shape and create notable errors in spectrum-integrated quantities. While this effect is one of a series of known potential sources of systematic error including, but not limited to target impurities, angular distributions, incident beam flux, and others, the error introduced from differences in environmental scattering effects between the reference and desired spectra are frequently overlooked. In this work, we demonstrate the origin of this error using simulated data corresponding to measurements of the neutron-induced prompt fission neutron spectra of 233 U, 235 U, 238 U, and 239 Pu as well as the spontaneous fission neutron spectrum of 252 Cf. Errors in the spectrum shape can be as high as 10%–15% depending on the combination of incident and outgoing neutron energy. The spectrum integral and average spectrum energy also contain errors that change with incident neutron energy and can be as high as ~2%. Although this study is based on simulated data, we also show a method for proving that this error exists in an experimental data set and we show a correction that can be applied to results of this type to obtain a more accurate result, thereby largely removing systematic errors from environmental scattering effects. The potential for these kinds of errors should be considered in measurements producing continuous neutron spectra.
The CGMF code implements the Hauser-Feshbach statistical nuclear reaction model to follow the de-excitation of fission fragments by successive emissions of prompt neutrons and γ rays. The Monte Carlo technique is used to facilitate the analysis of complex distributions and correlations among the prompt fission observables. Starting from initial configurations for the fission fragments in mass, charge, kinetic energy, excitation energy, spin, and parity, Y(A,Z,KE,U,J,π), CGMF samples neutron and γ-ray probability distributions at each stage of the decay process, conserving energy, spin and parity. Nuclear structure and reaction input data from the RIPL3 library are used to describe fission fragment properties and decay probabilities. Characteristics of prompt fission neutrons, prompt fission γ rays, and independent fission yields can be studied consistently. Correlations in energy, angle and multiplicity among the emitted neutrons and γ rays can be easily analyzed as a function of the emitting fragments. Finally, in its current version, CGMF is limited to certain spontaneous fission reactions on Pu and Cf isotopes, as well as neutron-induced fission reactions up to 20 MeV on U, Np, and Pu isotopes.
The cosmogenic ³⁹Ar(t1∕2 = 268 years) isotope of argon is used for geophysical dating and tracing of underground and ocean water, as well as ice owing to its appropriate half-life and chemical inertness as a noble gas; ³⁹Ar serves also in nuclear weapon test monitoring. We measured for the first time the total cross section of the main ³⁹Ar cosmogenic production reaction in the atmosphere, namely ⁴⁰Ar(𝑛, 2𝑛)39 Ar, using 14.8 ± 0.3 MeV neutrons. The neutrons, produced by a deuterium-tritium generator, impinged on a stainless steel sphere filled with Ar gas highly enriched in the ⁴⁰Ar isotope and were monitored by a stack of fast-neutron activation foils. The reaction yield was measured by atom counting of long-lived ³⁹Ar with noble gas accelerator mass spectrometry and, independently, by decay counting relative to atmospheric argon (³⁹Ar/Ar= 8.12 × 10−16 ). A total ⁴⁰Ar(𝑛, 2𝑛)39 Ar cross section of 610 ± 100 mb was determined at 14.8 ± 0.3 MeV incident neutron energy. This result serves as a benchmark for recent theoretical calculations and evaluations, found to reproduce well the experimental total cross section. We use these energy-dependent theoretical cross sections together with experimental spectra of cosmogenic neutrons at different altitudes to calculate the global average rate of neutron-induced ³⁹Ar atmospheric production, resulting in 770 ± 240 39 Ar atoms/cm²/day. The secular equilibrium between the ³⁹Ar calculated production rate and radioactive decay rate leads to a partial isotopic abundance ³⁹Ar/Ar = (5.9 ± 1.8) × 10−16 , showing that ≈73% of atmospheric ³⁹Ar is produced by cosmogenic neutrons, the remaining part believed to be induced by muons and high-energy 𝛾 rays. The ⁴⁰Ar(𝑛, 2𝑛)³⁹Ar cross section at 14 MeV is also a key parameter for quantifying the anthropogenic contribution to atmospheric 39 Ar produced during the thermonuclear tests of the 1960s. We estimate that anthropogenic ³⁹Ar accounts for roughly 20% of the present atmospheric inventory.
Abstract In this fifth paper of the series, we use the parameterized, spherically symmetric explosion method PUSH to investigate the impact of eight different nuclear equations of state (EOS). We present and discuss the explosion properties and the detailed nucleosynthesis yields, and predict the remnant (neutron star or black hole) for all our simulations. For this, we perform two sets of simulations. First, a complete study of nonrotating stars from 11 to 40 M ⊙ at three different metallicities using the SFHo EOS; and, second, a suite of simulations for four progenitors (16 M ⊙ at three metallicities and 25 M ⊙ at solar metallicity) for eight different nuclear EOS. We compare our predicted explosion energies and yields to observed supernovae and to the metal-poor star HD 84937. We find EOS-dependent differences in the explosion properties and the nucleosynthesis yields. However, when comparing to observations, these differences are not large enough to rule out any EOS considered in this work.
Anomalies observed in the neutron spectral shift of high-yield shots at the National Ignition Facility (NIF) suggest the presence of suprathermal ions [E. P. Hartouni et al., Nat. Phys. 19, 72 (2023)], implying that kinetic effects play a significant role in burning inertial confinement fusion (ICF) plasmas. Furthermore, recent measurements of reaction-in-flight (RIF) neutrons offer a direct probe of the stopping power in the burning fuel region of high-energy alpha particles and up-scattered fuel ions. We have developed the particle-in-cell code PICNIC, an exactly energy-conserving particle-in-cell Monte-Carlo collision (PIC-MCC) code to simulate the burn stage in ICF. We present results from 1D spherical simulations of NIF shot N210808, which was the first to exceed the Lawson criterion for laser fusion. We find that the suprathermal ions generated by large-angle Rutherford and nuclear elastic scattering (NES) with fusion alphas produce an alpha knock-on neutron (AKN) signal consistent with the extent and relative yield of the AKN spectrum identified in ignition experiments at the NIF. Furthermore, we also find that the inclusion of large-angle scattering physics does not explain the anomalously large spectral shift observed in the experiment.