SEARCH · Search NASA
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.
Data sets contain Rossi-alpha and pulsed neutron prompt neutron time decay measurements performed at the Oak Ridge Critical Experiments Facility (ORCEF) between 1960 and 1975. See also ORNL/TM-2019/1155 from April 2019 by John T. Mihalczo, included with these data sets.
The Feynman-alpha method is a neutron noise technique that is used to estimate the prompt neutron period of fissile assemblies. The method and quantity are of widespread interest including in applications such as nuclear criticality safety, safeguards and nonproliferation, and stockpile stewardship; the prompt neutron period may also be used to infer the k eff multiplication factor. The Feynman-alpha method is predicated on time-correlated neutron detections that deviate from a Poisson random variable due to multiplication. Traditionally, such measurements are diagnosed with one-region point kinetics, but two-region models are required when the fissile assembly is reflected. This paper presents a derivation of the two-region point kinetics Feynman-alpha equations based on a double integration of the Rossi-alpha equations, develops novel propagation of measurement uncertainty, and validates the theory. Validation is achieved with organic scintillator measurements of weapons-grade plutonium reflected by various amounts of copper to achieve k eff values of 0.83–0.94 and prompt periods of 5–75 ns. The results demonstrate that Feynman-alpha measurements should use the two-region model instead of the one-region model. The simplified one-region model deviates from the validated two-region models by as much as 10% in the estimate of the prompt neutron period, and the two-region model reduces to the one-region model for small amounts of reflector. The Feynman-alpha estimates of the prompt neutron period are compared to those of the Rossi-alpha approach. The comparative results demonstrate that the Feynman-alpha method is more precise than the Rossi-alpha method and more accurate for k eff < 0.92, whereas the Rossi-alpha method is generally more accurate for higher multiplications. Here, the uncertainty propagation developed in this work should be used for all Feynman-alpha measurements and will therein improve fitting accuracy and appropriate precision estimates.
Nine 7, 11, and 15 in. diam highly enriched uranium (HEU, 93.15 wt % 235 U) metal cylinders were assembled on the vertical assembly machine in the Oak Ridge Critical Experiments Facility (ORCEF) and had 1, 2, or 3 in. thick HLM graphite reflectors on the top and bottom. The experiments, which were performed between April 3, 1970, and February 18, 1971, used 23 operational days at ORCEF. Before those experiments were carried out, unreflected and unmoderated, graphite- and polyethylene-reflected, and polyethylene-moderated HEU metal cylinders had been assembled to obtain delayed criticality at ORCEF in the 1960s and reported by the International Criticality Safety Benchmark Evaluation Project (ICSBEP) at the Nuclear Energy Agency (NEA)*. The data from the nine critical experiments are acceptable for use as criticality safety benchmark experiments for the NEA’s ICSBEP once the uncertainty analysis is completed. Based on previous ICSBEP benchmarks with HEU metal at ORCEF, the uncertainties in k eff are expected to be as low as ±0.0004.
Explore the source record for details and available documents.
For FY20, we had a deliverable to write a report detailing efforts to simultaneously evaluate both the prompt neutron multiplicity, $\overline{ν}$, and the prompt neutron fission neutron spectrum, PFNS, using CGMF. CGMF is the LANL-developed fission fragment decay code that consistently evaporates prompt neutrons and γ rays using the Hauser-Feshbach statistical theory of compound reactions. The decay begins by constructing the initial conditions of the fission fragments, then decaying each one from the excited state by neutrons and γ rays, conserving energy, momentum, spin, and parity in each step of the emission. The initial conditions of the fragments, along with the multiplicity, energy, and direction of each emitted neutron or γ ray, are recorded, allowing for the full reconstruction of the fission event. These event histories allow us to reconstruct average quantities, as well as correlations between observables, that can be compared with experimental or evaluated data. In that initial report, although there was already a favorable comparison between $\overline{ν}$ from CGMF, experiment, and the current ENDF/B-VIII.0 evaluation, we showed that there was still significant work to be done to improve the PFNS from CGMF. Historically, the PFNS is calculated too soft by Hauser-Feshbach fission models, and CGMF is no exception. The incorrect shape presents a significant challenge in fission modeling, including for our understanding of the fission process and for our ability to consistently calculate and predict a variety of prompt fission observables (such as fission fragment initial conditions, neutron and γ-ray multiplicities and energies, and the correlations between all observables). In our companion report, we detail our success in using CGMF to evaluate $\overline{ν}$. Although not included in the optimization explicitly, we also keep the initial conditions of the fission fragments physical, along with reproducing reasonably well the neutron multiplicity distribution. As we would expect from the sensitivities calculations from, the average neutron energies change very little from the $\overline{ν}$ optimization along with the PFNS (as will be shown in Sec. 2.6). The conclusion was that the global and statistical models would have to be investigated instead of just the fission fragment initial conditions (as is sufficient for $\overline{ν}$). This report details those efforts.
We report the Hauser-Feshbach fission-fragment decay (HF 3 D) model is extended to calculate the prompt fission neutron spectrum (PFNS) for the thermal-neutron-induced fission on 235 U, where the evaporated neutrons from all possible fission-fragment pairs are aggregated. By studying model parameter sensitivities on the calculated PFNS as well as nonstatistical behavior of low-lying discrete level spin distribution, we conclude that discrepancies between the aggregation calculation and the experimental PFNS seen at higher neutron emission energies can be attributed to both the primary fission-fragment yield distribution and the possible high spin states that are not predicted by the statistical theory of nuclear structure.
The Hauser-Feshbach fission fragment decay model, HF3D, calculates the statistical decay of fission fragments through both prompt and delayed neutron and γ-ray emissions in a deterministic manner. While previously limited to the calculation of only first-chance fission, the model has recently been extended to include multi-chance fission, up to neutron incident energies of 20 MeV. The deterministic decay takes as input prescission quantities–fission probabilities, pre-fission neutron energies, and the average energy causing fission– and post-scission quantities–yields in mass, charge, total kinetic energy, spin, and parity. From those fission fragment initial conditions, the full decay is followed through both prompt and delayed particle emissions. The evaporation of the prompt neutrons and γ rays is calculated through the Hauser-Feshbach statistical theory, taking into account the competition between neutron and γ-ray emission, conserving energy, spin, and parity. The delayed emission is taken into account using time-independent calculation using decay data. This whole formulation allows for the calculation of prompt neutron and γ-ray properties, such as multiplicities and energy distributions, both independent and cumulative fission yields, and delayed neutron observables, in a consistent framework. Here, we describe the implementation of multi-chance fission into the HF 3 D model, and show an example of prompt and delayed quantities beyond first-chance fission, using the example of neutron-induced fission on 239 Pu. This expansion represents significant progress in consistently modeling the emission of prompt and delayed particles from fissile systems.
This report documents the analysis of dose rate fields from a combination of two sources stemming from the target system of the Oak Ridge National Laboratory (ORNL) Spallation Neutron Source (SNS) Second Target Station (STS). This analysis supports the target drive components’ lifetime estimates and conventional facilities’ preliminary design. The two sources detailed in this report are the radiation generated by the proton beam incident on the target and the activated cooling water in the pipes at the top of the target assembly. The analysis considers the STS target system in operational configuration, with all shielding in place. The dose rate fields are calculated in the target drive room and through the target drive room ceiling. This report calculates the thickness necessary for the target drive room’s high-density concrete ceiling to meet the 0.25 mrem/hr dose limit requirement in the high bay directly above the target drive room. In addition to sizing the target drive room ceiling, this report includes details of the dose field surrounding the target drive assembly to estimate the deposited dose in the target drive assembly components resulting from both streaming and decay of activated water. Prompt neutrons and gammas generated by the proton beam incident on the target primarily stream along the gaps between the target assembly and core vessel shielding. Because the models and analysis detailed in this report are based on the STS preliminary design, future analysis is needed to support the STS final design.
In this article, we perform energy-dependent calculations of independent and cumulative fission product yields for 235 U, 238 U, and 239 Pu in the first chance fission region. Starting with the primary fission fragment distributions taken from available experimental data and analytical functions based on assumptions for the excitation energy and spin-parity distributions, the Hauser-Feshbach statistical decay treatment for fission fragment de-excitation is applied to more than 1,000 fission fragments for the incident neutron energies up to 5 MeV. The calculated independent yields are then used as an input of β-decay calculations to produce the cumulative yield, and summation calculations are performed. Model parameters in these procedures are adjusted by applying the Bayesian technique at the thermal energy for 235 U and 239 Pu and in the fast energy range for 238 U. The calculated fission observable quantities, such as the energy-dependent cumulative yields, and prompt and delayed neutron yields, are compared with available experimental data. We also study the possible impact of the second chance fission opening on the energy dependence of the delayed neutron yield by on the opening of second chance fission on the energy dependence calculation.
Throughout the course of FY21, significant effort was put into investigating models within the LANL developed Hauser-Feshbach fission fragment decay code, CGMF, to understand and potentially solve the long-standing challenge of a too-soft prompt fission neutron spectrum, PFNS. Several inputs and models to CGMF were investigated, including the discrete nuclear levels, the optical model potential, level densities, and the fission fragment initial conditions. Some of the global models within CGMF led to a slight hardening of the neutron spectrum—particularly the likely incomplete discrete levels in through which γ-rays decay—but none of the changes where large enough for the tail of the PFNS to reproduce experimental data. A significant hardening of the spectrum tail was observed when the fission fragment initial conditions were optimized based on their sensitivities to the PFNS data for thermal incident neutrons. In this way, the parameters for the CGMF mass and total kinetic energy distributions, along with the spin cutoff factor were adjusted to better reproduce the experimental PFNS measurements. This optimization hardened the tail of the PFNS slightly but led to unphysical mass distributions for the fission fragments before neutron emission. It was clear from the above that we do not expect to be able to produce an evaluation-quality PFNS with CGMF in the near future. Challenges at thermal will persist–and possibly worsen–with increasing incident energy, where more models are needed to completely describe the fission. Basic-science research funding exceeding the amount available and scope of our NCSP funds would be needed to tackle this decade-long challenge impacting many fission-fragment event generator. And, in fact, Amy Lovell won LDRD ECR funding to do so over the next few years. Therefore, we focused in FY22 on extending evaluation capabilities beyond thermal incident neutrons, to take into account the incident energy dependence of the PFNS and fission fragment initial condition distributions in CGMF. We chose to set up the evaluation methodology to perform PFNS evaluations with CGMF across incident-neutron energies, in order to have it readily available for future NCSP evaluations when the PFNS from CGMF has improved. In this report, we outline the evaluation methodology, along with the results of the optimization, including full model calculations with CGMF using the evaluated parameters.
Advanced instruments and methods need to be developed now to create a technical basis to support the negotiation of future nuclear arms control treaties. One new capability that is anticipated is the ability to confirm either the declared presence or declared absence of high explosive (HE) material in the presence of special nuclear material (SNM). Towards this goal, Passive HE Neutron Inspection (PHENIX) has been developed and demonstrated as a method for confirming the presence or absence of HE in the presence of plutonium. The method exploits the inherent presence of neutrons associated with the decay of plutonium as an internal probe source for performing prompt gamma-ray neutron activation analysis (PGNAA), searching for the presence of HE as revealed by the emission of characteristic gamma rays following neutron absorption in hydrogen and nitrogen which are building blocks of present-day, military-grade HE. Tests using stoichiometrically-correct hemishells of mock HE with plutonium show that a system can be expected to positively confirm the presence or absence of these signatures, supporting determination of HE presence or absence with Pu, in a few hours. To protect other potentially sensitive gamma-ray signatures from a treaty accountable item, an analog information barrier has been conceptualized and tested which physically prevents the collection of gamma-ray spectral data outside of user selected energy windows strategically chosen to view only narrow spectral regions corresponding to the hydrogen (2223.2 keV) and nitrogen (9807.2 keV, 10,318.2 keV, and 10,829.2 keV) PGNAA signatures.
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.
Remote sensing of ionizing radiation has a significant role in waste management, nuclear material management and nonproliferation, and radiation safety. Robotic platforms can surpass the number of tasks that are achieved by humans. With this technique, the operator's radiation exposure can be decreased. Remote sensing allows for the evaluation and monitoring of radiological contamination. Gamma-ray and neutron sensors were integrated onto the robotic platforms. This approach allows for the radiation sensor data to be dynamically tracked and mapped thus enabling further analysis of the radiation flux in temporal and spatial domains. The goal is to complete scheduled tasks while the robot is being irradiated. To achieve this, electronic components must be shielded and radiation hardened. CZT Detector: Cadmium Zinc Telluride (CZT) detector technology has been a promising solution for gamma-ray and x-ray measurements. Detector data is transferred to the Odroid minicomputer that controls and powers the module via the USB. Robot Operating System (ROS) was utilized for data acquisition and data fusion. The Mariscotti method was employed for the spectrum analysis. A function was programmed in ROS for the automatic identification of photopeaks. CLYC Detector: A Cs{sub 2}LiYCl{sub 6}:Ce{sup 3+} (CLYC) detector was used for simultaneous medium-resolution gamma-ray measurements and neutron counting. A 2.54 cm diameter photomultiplier tube (PMT) was equipped with a high voltage supply and a miniature digitizer. Gamma-ray excitation: fast core-to-valence luminescence (CVL) with 1 ns decay constant, and prompt Ce{sup 3+} emission with 50 ns decay constant. Neutron excitation: slow cerium self-trapped excitation (Ce{sup 3+} STE), 1000 ns decay constant. Radiation Source Localization: Maximum Likelihood Estimation (MLE) and gradient-based methods were used to locate the position of a radiation source based on measured radiation intensities. Multi-Particle Transport Code FLUKA: Estimation of radiation damage of the electronic components is important in order to optimize the robot's operational time while it is irradiated. Displacement per atom (DPA) represents the radiation damage in materials exposed to the ionizing radiation. Various shielding layers of different thickness t were analyzed (< 5% statistical error). The model of the controller of the UAS was designed in FLUKA. Conclusion: CZT and CLYC detectors were integrated onto the robotic platforms. Radiation source localization and contour mapping using robotic platforms were studied. Functions for data analysis and fusion were developed in ROS. FLUKA code was utilized to analyze DPA values. Layers of low-density and high-density materials were used to shield the UAS electronics.
Prompt and β-delayed fission observables, such as the average number of prompt and delayed neutrons, the independent and cumulative fission product yields, and the prompt γ-ray energy spectra for the photonuclear reactions on 235,238 U and 239 Pu, are calculated with the Hauser-Feshbach fission fragment decay (HF 3 D) model and compared with available experimental data. Further, in the analysis of neutron-induced fission reactions to the case of photo-induced fission, an excellent reproduction of the delayed neutron yields supports a traditional assumption that the photo fission might be similar to the neutron-induced fission at the same excitation energies regardless of the spin and parity of the fissioning systems.
The neutron-deficient 188 Bi and 188 Po isotopes have been studied by γ-ray spectroscopy using the recoil-decay tagging technique with the Argonne Gas-Filled Analyzer. A new 0.25(5)-μs isomeric state and a prompt cascade formed by 319-, 366-, and 462-keV γ rays have been established on top of the (10¯) α -decaying isomer in 188 Bi. The first excited (2 + ) state in 188 Po was identified, its excitation energy of 242(2) keV continues the nearly constant trend for the first 2 + states in 190, 192,194 Po. Here, the state is most likely a member of a prolate rotational band built on the ground state, albeit mixing with other coexisting configurations cannot be excluded. The new results obtained in the present work provide new information to shape coexistence in bismuth and polonium isotopes near the neutron midshell at N = 104 . In this mass region, a reduction in the prompt γ-ray yield obtained with recoil decay tagging was observed for a few nuclides, and the possible reasons are presented.
The Hauser-Feshbach fission fragment decay model, HF 3 D, which calculates the statistical decay of fission fragments, has been expanded to include multichance fission, up to neutron incident energies of 20 MeV. The deterministic decay takes as input prescission quantities—fission probabilities and the average energy causing fission—and postscission quantities—yields in mass, charge, total kinetic energy, spin, and parity. From these fission fragment initial conditions, the full decay is followed through both prompt and delayed particle emissions, allowing for the calculation of prompt neutron and γ properties, such as multiplicity and energy distributions, both independent and cumulative fission yields, and delayed neutron observables. Here, we describe the implementation of multichance fission in the HF 3 D model, and show an example of prompt and delayed quantities beyond first-chance fission, using the example of neutron-induced fission on 235 U. This expansion represents significant progress in consistently modeling the emission of prompt and delayed particles from fissile systems.
We present the development and application of a collision tracking feature within the OpenMC Monte Carlo particle transport code, designed for diverse applications such as neutron spectroscopy, scatter camera system, neutron noise, and multiplicity counting simulations. This feature enables the tracking of individual particle collisions, with potential applications in nuclear nonproliferation, reactor physics, and nuclear security. Additionally, the feature holds potential for the calibration of neutron detectors, specifically in converting light output into energy deposited within the detectors. The implementation consists of a set of filters—such as reaction type, energy, cell, and material—that constrain the set of collisions that are tracked, extensions to the Python API to enable simple input specification, and support for writing either OpenMC’s native HDF5-based format or the Monte Carlo particle list format. This feature was added to the official OpenMC release in version 0.15.3. In this work, the feature will be applied to showcase scenarios such as time-of-flight simulations, scatter-camera imaging for neutron source localization, neutron-noise analysis to extract integral kinetic parameters such as the prompt decay constant α, and multiplicity counting to estimate the mass of special nuclear materials. Ultimately, this feature aims to expand the application scope of open-source Monte Carlo particle transport codes such as OpenMC.
Abstract Fusion power reactors will generate intense neutron fluxes into plasma-facing and structural materials (SMs). Radiation back-fluxes, generated from neutron–material interactions under these fluxes, can dramatically impact the plasma dynamics, e.g. by seeding runaway electrons during disruptions via Compton scattering of background electrons by wall-emitted gamma radiation. Here, we quantify these back-fluxes, including neutrons, gamma rays, and electrons, using Monte Carlo calculations for a range of SM candidates and first wall (FW) thicknesses. The radiation back-flux magnitudes are remarkably large, with neutron and gamma radiation back-fluxes on the same order of magnitude as the incident fusion neutron flux. Electron back-fluxes are two orders of magnitudes lower, but are emitted at sufficiently high energies to impact the sheath and boundary plasma dynamics. Material configuration plays a key role in determining back-flux magnitudes. The SM chiefly determines the neutron back-flux magnitude, while the FW thickness principally attenuates the gamma ray and electron back-fluxes. In addition to prompt back-fluxes, which are emitted immediately after fusion neutrons impact the surface, significant delayed gamma ray and electron back-fluxes arise from nuclear decay processes in the activated materials. These delayed back-flux magnitudes range from 2% to 7% of the prompt back-fluxes, and remain present during transients when fusion no longer occurs. During disruptions, build-up of delayed gamma radiation back-flux represents potential runaway electron seeding mechanisms, posing additional challenges for disruption mitigation in a power reactor compared with non-nuclear plasma operations. This work highlights the impact of these radiation back-fluxes plasma performance and demonstrates the importance of considering back-flux generation in materials selection for fusion power reactors.