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.
Major evaluation efforts have been completed or are underway for actinide fission observables, including fission product yields, PFNS, neutron multiplicity, and neutron-induced cross sections. Energy-dependence and covariances for fission product yields are being developed due to new measurements and modeling efforts.
We present a new neutron multiplicity counting analysis and measurement method for neutron-shielded fissile material using neutron-capture gamma rays. Neutrons absorbed in shielding produce characteristic gamma rays that preserve the otherwise lost neutron multiplicity signature. Neutron multiplicity counting provides estimates of fission parameters, such as neutron leakage multiplication, spontaneous fissioner (e.g. Pu-240) mass, and (α,n) ratio. Standard neutron multiplicity counting can incorporate the new neutron-capture gamma-ray multiplicity counting technique to characterize previously degenerate or intractable source configurations by maximizing the multiplicity signature. The new method decouples neutron source-detector interferences, such as reflection and thermalization time in the detector, that could improve measurements of the mean neutron lifetime. We also develop a detector prototype for the multiplicity counting of neutron-capture gamma rays and present detector design considerations, such as detection material and shielding, to optimize the detection of the 2.2 MeV hydrogen capture gamma ray. We simulate the prototype neutron- capture gamma-ray multiplicity counter against the BeRP ball in polyethylene shells to inform future measurements.
In this presentation, researchers detail tradeoffs between two simulation suites for thermal neutron scattering data. OXLIMAX + MCViNE provide a more accurate elastic peak shape and NJOY + MCNP have more accurate spectra between phonon mode peaks (most likely due to multiple neutron scattering). Additionally, some outstanding issues with MCViNE are noted, including multiple neutron scattering and non-uniform (Q,E) grid for data.
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 differentiate the event-by-event n-γ multiplicity data from 252 Cf(sf) with respect to the energies of the emitted particles as well as their relative angles of emission. We determine that neutron emission enhances γ-ray emission around 0:7 and 1:2 MeV, but the only directional alignment was observed for E γ ≤ 0:7 MeV and tended to be parallel and antiparallel to neutrons emitted in the same event. The emission of γ rays at other energies was determined to be nearly isotropic. The presence of the emission and alignment enhancements is explained by positive correlations between neutron emission and quadrupole γ-ray emission along rotational bands in the de-exciting fragments. This observation corroborates the hypothesis of positive correlations between the angular momentum of a fragment and its intrinsic excitation energy. The results of this work are especially relevant in view of the recent theoretical and experimental interest in the generation of angular momentum in fission. Specifically, we have determined an alignment of the fragments angular momenta in a direction perpendicular to the direction of motion. We interpret the lack of n-γ angular correlations for fission fragments near closed shells as a weakening of the alignment process for spherical nuclei. Lastly, we have observed that statistical γ rays are emitted isotropically, indicating that the average angular momentum removed by this radiation is small. These results, and the analysis tools presented in this work, represent a stepping stone for future analysis of n-γ emission correlations and their connection to angular momentum properties.
This report documents an evaluation of 235 U prompt fission neutron spectra (PFNS) that is a release candidate for the upcoming U.S. nuclear data library, ENDF/B-VIII.1. This evaluation differs from its predecessor, ENDF/B-VIII.0, mainly by the inclusion of 235 U PFNS measured by the Chi-Nu team of LANL and LLNL. This data set is the first one that covers the 235 U PFNS for continuous incident-neutron energies of 1⁻20 MeV and outgoing-neutron energies from 10 keV⁻10 MeV with high precision. Previous data sets were either measured in a limited energy range or with less precision. Hence, these new Chi-Nu data provide decisive information for the evaluation. The resulting evaluated data correspond well to the new experimental PFNS. The evaluated PFNS also produce average mean energies and 239 Pu/ 235 U PFNS in agreement with associated Chi-Nu data. If one uses the new evaluated data to predict the neutron multiplication factor, k eff , of selected ICSBEP critical assemblies, the differences of simulated values compared to those using ENDF/B-VIII.0 is modest (less than 55 pcm). This difference in k eff can be easily accommodated by changes in the 235 U average prompt fission neutron multiplicity that is currently being re-evaluated. In addition to that, the new PFNS predict on average 235 U LLNL pulsed-sphere neutron-leakage spectra better than ENDF/B-VIIII.0 PFNS.
We reinvestigated the neutron multiplicity yields of Ba-Mo, Ce-Zr, Te-Pd, and Nd-Sr from the spontaneous fission of 252 Cf by (i) using both γ -γ -γ -γ and γ -γ -γ coincidence data, (ii) using up to date level scheme structures, and (iii) crosschecking analogous energy transitions in multiple isotopes, we have achieved higher precision than previous analyses. Particular attention was given to the Ba-Mo pairs where our results clearly confirm that the Ba-Mo yield data have a second hot fission mode where 8–10, and now 11 neutron evaporation channels are observed. These are the first observations of the 11 neutron channel. These 8–11 neutron channels are observed for the first time in the Ce-Zr pairs, but are not observed in other fission pairs. The measured intensities of the second mode in Ba-Mo and Ce-Zr pairs are ~1.5(4)% and ~1.0(3)%, respectively. These high neutron number evaporation modes can be an indication of hyperdeformation and/or octupole deformation in 143-145 Ba and in 146,148 Ce at scission to give rise to such high neutron multiplicities.
The stochastic theory of neutron transport is extended to describe the cumulative distribution of fission numbers and deposited fission energy in a multiplying assembly. Solutions for the probability distributions are obtained using analytical approximations and Monte Carlo simulation in lumped geometry and in symmetric homogeneous and heterogeneous spheres. The results show the development of a power-law tail in the steady state fission number and deposited energy distributions when the medium is critical, independent of the fission neutron multiplicity distribution and domain heterogeneity. In contrast, the asymptotic decay is faster than exponential in subcritical media due to rapid chain extinction and in supercritical media due to the increasing probability of chain divergence. Here, a formal asymptotic analysis of the problem in lumped geometry with an arbitrary fission neutron multiplicity confirms the existence of power-law tails at critical.
These experiments were performed in 1983 between August 1 and September 16 and used 35 operational days of critical facility time at the Babcock & Wilcox critical facility in Lynchburg, Va. These measurements were to assess the capability of the Cf source driven noise analysis method to measure the subcriticality of light water reactor fresh fuel configurations. This report documents the experiments of 1983 that were not reported at that time. These measurement with 2.459 wt. % 235 U enriched uranium oxide fresh fuel pins started at the delayed critical configuration of the fully assembled pressurized water reactor (PWR) core configuration with 1511 parts per million (ppm) natural boron (19.8 wt. % 10 B) in the water. This was a cylindrical configuration of 4961 fuel pins. After the initial measurements at delayed criticality, subcritical measurements by the Cf source driven noise analysis method were performed in which the number of detectors and Cf source locations were varied. Additional boron was added to the water moderator and reflector to reduce the neutron multiplication factor, with resulting boron concentrations of 1561, 1613, 1765, 1880, 2104, 2384, 2975, 3606, and 4303 ppm. After the completion of these measurements the increased, the boron concentration was returned to that at delayed critical (1511 ppm) and the reactivity was decreased by reducing the number of fuel pins in steps to 3713, 2553, 1281, 749, and 333 fuel pins and then to that (289) of a s ingle square 17 x 17 pin array, corresponding to a single PW fuel element. For all measurement that started at delayed criticality and then reduced the neutron multiplication factor to various subcritical states, inverse kinetics rod drop measurements at the higher reactivities were also performed but not reported here. At all subcritical states the break frequency noise analysis data was fitted to obtain the prompt neutron decay constant. Some data presented in this report are from notes not in the logbook. The purpose of this report is to document the experimental information for the measurements performed so that a later date someone could perform the required uncertainty and calculational analyses and documentation to use these data for an International Nuclear Criticality Safety Benchmark Program (INCSBEP) or Nuclear Energy Agency (NEA) benchmark. The data from these measurements is available from the Laboratory Records Department of Oak Ridge National Laboratory and the logbook is available from INCSBEP at Idaho National Laboratory.
Since 1945, Los Alamos National Laboratory (LANL) has performed critical experiments, primarily at the Los Alamos Critical Experiments Facility (LACEF) at Technical Area 18 (TA-18). Since 2011, these experiments have been conducted by the National Criticality Experiments Research Center (NCERC), operated by LANL, at the Device Assembly Facility (DAF) in the Nevada National Security Site (NNSS). These experiments utilize various types of Special Nuclear Material (SNM). Some of these experiments utilize what is referred to as the Rocky Flat Shells, which are called such as they came from the Rocky Flats Plant. These concentric hemi-shells are made of Highly Enriched Uranium (HEU), which is 93 w/o 235U. Fig. 1. Subset of the Rocky Flat Shells LANL is designing a new subcritical hands-on experiment with the goal of achieving a neutron multiplication in the range of 50 to 200, which correlates to Keff values of 0.98 to 0.995. ANSI/ANS-1 is the standard for Conduct of Critical Experiments which governs critical operations at NCERC. Section 3.9 of ANSI/ANS-1 states that when manipulating a critical assembly by hand, the predicted k eff of a known configuration should not exceed 0.95 (a neutron multiplication of 20). This presents a challenge, as this system would have a higher multiplication than 20 and the assembly would have to remain subcritical in normal and credible accident scenarios. To ensure the safety of such an assembly, the different normal and abnormal conditions that could alter the criticality 1 MCNP® and Monte Carlo N-Particle® are registered trademarks owned by Triad National Security, LLC, manager and operator of Los Alamos National Laboratory. Any third party use of such registered marks should be properly of the system must be considered and analyzed. If such a condition is deemed to be credible, it will be further investigated using the Monte Carlo N-Particle (MCNP) transport code. If the results of these simulations show that a k eff larger than one could be possible, modifications to the assembly will be made until such an event is deemed no longer possible.
Nuclear physics advances in the United States and Britain from 1939 to 1945 are described. The Manhattan Project’s work led to an explosion in our knowledge of nuclear science. A conference in April 1943 at Los Alamos provided a simple formula used to compute critical masses and laid out the research program needed to determine the key nuclear constants. In short order, four university accelerators were disassembled and reassembled at Los Alamos, and methods were established to make measurements on extremely small samples owing to the initial lack of availability of enriched 235U and plutonium. I trace the program that measured fission cross sections, fission-emitted neutron multiplicities and their energy spectra, and transport cross sections, comparing the measurements with our best understanding today as embodied in the Evaluated Nuclear Data File ENDF/B-VIII.0. The large nuclear data uncertainties at the beginning of the project, which often exceeded 25% to 50%, were reduced by 1945 often to less than 5% to 10%. Uranium-235 and plutonium-239 fission cross-section assessments in the fast mega-electron-volt range were reduced following more accurate measurements, and the neutron multiplicity $\overline{v}$ increased. By a lucky coincidence of canceling errors, the initial critical mass estimates were close to the final estimated masses. Some images from historical documents from our Los Alamos archives are shown. Many of the original measurements from these early years have not previously been widely available. Through this work, these data have now been archived in the international experimental nuclear reaction data library (EXFOR) in a collaboration with the International Atomic Energy Agency and Brookhaven National Laboratory.
Compensating errors between nuclear data observables in a library can adversely impact application simulations. The primary goal of the EUCLID project (Experiments Underpinned by Computational Learning for Improvements in Nuclear Data) is to reduce compensating errors in nuclear data. A new criticality experiment, described in this work, was designed with the specific target nuclear data of 239 Pu fission, inelastic scattering, elastic scattering, capture, nu-bar, and prompt fission neutron spectrum (PFNS). This work will focus on the design and execution of the EUCLID experiment, performed on the Planet vertical lift critical assembly machine at the National Criticality Experiments Research Center (NCERC). The criticality experiment includes two different configurations with very different geometries: one is cube-like to minimize neutron leakage while the other is slab-like to maximize leakage. Having these two widely varying configurations allows the scattering sensitivities of 239 Pu to the neutron multiplication factor to be greatly changed while minimally impacting the other cross section sensitivities. Both configurations utilize the Pu ZPPR (Zero Power Physics Reactor) plates as fuel. The experiments were designed using a D-Optimality criteria, which is an optimization method minimizing the log-determinant of the adjusted nuclear data covariance for the target reactions. These experiments include not only inference of k eff , as done in all critical benchmark experiments, but several other responses as well, such as neutron multiplication measurements and reaction rate ratios. After analysis of the measured data is complete, adjustment of nuclear data will be performed to assess whether the new experimental data successfully reduced compensating errors.
Spent tri-structural isotropic (TRISO)–based fuels have a strong track record in storage and transportation without documented incidents. This work seeks to reduce uncertainty to aid in more informed spent fuel management of TRISO-based fuels by modeling both fresh and spent pebble bed reactor (PBR) fuel and comparing the results to the regulatory standards from 10 CFR 71. SCALE was used for all modeling due to it having fast and accurate methods for handling PBR fuel modeling, as well as having an efficient method for shielding calculations in monaco with automated variance reduction using importance calculations (MAVRIC), which utilizes the consistent adjoint-driven importance sampling (CADIS) and the forward-weighted consistent adjoint-driven importance sampling (FW-CADIS) methods. KENO-VI was used for all criticality calculations, TSUNAMI was used for uncertainty quantification on k-effective, TRITON and the Oak Ridge isotope generation code (ORIGEN) were both used for depletion of the fuel, and MAVRIC was used for shielding calculations. For criticality assessments, this study focused on the requirement that the value of the neutron multiplication factor, k-effective (k-eff), would not exceed a peak value of 0.95, including uncertainty, with 95% confidence. Criticality was initially examined by modeling fresh fuel from three different designs—HTR-10 fuel, PBMR-400 fuel, and demonstration fuel representative of a TRISO-fueled modern high-temperature gas reactor (HTGR) design, henceforth referred to as Demo HTGR—and placing them into various sized containers with conditions described in 10 CFR 71 to quantify the peak k-eff state. When the peak value of 0.95 k-eff was exceeded, mitigation methods were examined in those scenarios. Burnup credit, pebble displacement in areas of strong neutron multiplication, and random pebble replacement using pebbles of various compositions and replacement fractions were examined. In summary, the criticality of PBR fuels can be well accounted for by restricting container size, taking credit for burnup, or by displacing/replacing pebbles. Uncertainty of the k-eff due to nuclear data uncertainties was recorded at ~0.6644%Δk/k, or roughly 664% mil (pcm). The nuclear data–induced uncertainty was relatively small and should not require significant modification in the design to be accounted for. Revisions to the evaluated nuclear data file values have been shown to have a larger impact than nuclear data–induced uncertainty. For dose rate aspects, U.S. Nuclear Regulatory Commission regulations require a maximum dose rate of 10 millirem per hour (mrem/h) at 2 meters. In examining the dose rate behavior of spent PBR fuel, the representative Demo HTGR fuel was modeled exclusively due to it possessing the highest target burnup of the examined fuels. Equilibrium cycle modeling methods were used to produce a higher-fidelity discharge isotopic composition than simple assumptions, such as reflected pebbles. The discharge composition was used as a source term in the fixed-source transport shielding calculations, and dose rates were calculated at 2 m for the shortest possible cooling time. The low concentration of fuel material led to dose rates that were in line with regulatory limits, despite the high burnup when compared to traditional light water reactor fuels. In conclusion, the methods employed in this study would require more work to further verify and validate and are limited to the criticality and dose rate analyses performed.
After successful fitting of prompt ν for 235 U, 238 U, and 239 Pu(n,f) using CGMF, we have moved on to the minor plutonium isotopes. Minor isotopes pose a greater challenge both because there is less data available and we do not, by default, have parametrizations in CGMF already. Initially, we had planned to mitigate these challenges—and provide consistency within CGMF—by taking a stepped approach to the optimization. Continuing from the 239 Pu work, we would then fit 241Pu(n,f) $\overline{ν}$, where there is also a number of experimental measurements, keeping consistency between the parameters that are included in the calculation for 239 Pu and 241 Pu. With these parametrizations settled, we could consistently optimize 240 Pu(n,f) $\overline{ν}$. Following that step, we would move on to 242 Pu(n,f) (and increasing neutron number) and 238 Pu(n,f) (and decreasing neutron number). By including the fissioning systems in this manner, we should be able to minimize the unknown parameters in CGMF. We were also able to develop systematics for the CGMF input parameters along the Pu isotopic chain. However, after studies in the beginning of FY23, where we found reasonable agreement between CGMF and experimental prompt neutron multiplicities for 238−242 Pu, using a compound mass dependent parametrization from, and from discussions with I. Stetcu and P. Talou, we instead performed the evaluation of all five isotopes simultaneously. For all parameters in CGMF, they are either the same for each compound nucleus, or have a dependence on the compound mass. The details of these updates are given in. This report builds upon those details. In this short report, we first give an example of how we have updated CGMF to include 240 Pu(n,f) and 242 Pu(n,f), keeping consistency with the current 239 Pu and 241 Pu calculations, but without rigorous optimization (Sec. 2). Then, we will show in Section 3 what experimental data exist for the various isotopes to provide some insight into why 241 Pu and 239 Pu are used as anchor points. The evaluation is detailed in Section 4, and the comparison between CGMF calculations using the evaluated parameters and other prompt observables besides average neutron multiplicity are discussed in Section 2.2. Finally, we present conclusions and future work in Section 5.
The covariance committee of CSEWG (Cross Section Evaluation Working Group) established templates of expected measurement uncertainties for neutron-induced total, (n,γ), neutron-induced charged-particle, and (n,xn) reaction cross sections as well as prompt fission neutron spectra, average prompt and total fission neutron multiplicities, and fission yields. Templates provide a list of what uncertainty sources are expected for each measurement type and observable, and suggest typical ranges of these uncertainties and correlations based on a survey of experimental data, associated literature, and feedback from experimenters. Information needed to faithfully include the experimental data in the nuclear-data evaluation process is also provided. These templates could assist (a) experimenters and EXFOR compilers in delivering more complete uncertainties and measurement information relevant for evaluations of new experimental data, and (b) evaluators in achieving a more comprehensive uncertainty quantification for evaluation purposes. This effort might ultimately lead to more realistic evaluated covariances for nuclear-data applications. In this topical issue, we cover the templates coming out of this CSEWG effort–typically, one observable per paper. This paper here prefaces this topical issue by introducing the concept and mathematical framework of templates, discussing potential use cases, and giving an example of how they can be applied (estimating missing experimental uncertainties of 235 U(n,f) average prompt fission neutron multiplicities), and their impact on nuclear-data evaluations.
MicroBooNE, an 85-tonne liquid argon time projection chamber (LArTPC) detector is on-axis to the Booster Neutrino Beam (BNB) beamline facility at Fermi National Accelerator Laboratory. MicroBooNE is elucidating neutrino interactions with argon through cross-section measurements to refine interaction models and reduce uncertainties. In this poster, we present the status of the single and double multi-differential charged current (CC) cross section with zero pions in the final state (CC-0$\pi$) as a function of muon momentum ($0.1<p_\mu<2.0\,\mathrm{GeV/c}$) and the cosine of the muon angle ($-1<\cos\theta_\mu<1$). We present the details of the event selection and cross section extraction along with a set of tests using fake data to establish the robustness of the analysis methodology. We also discuss prospects for a future combined measurement with the Gd-H$_2$O target at the ANNIE experiment, to explore MicroBooNE’s proton multiplicity alongside ANNIE’s neutron multiplicity.
Explore the source record for details and available documents.
This study newly proposes a heterogeneous core calculation method based on local/global iteration using proper orthogonal decomposition (POD). By using the singular value decomposition (SVD) and the low-rank approximation, appropriate POD bases for expanding the neutron flux can be obtained from snapshot data of the neutron flux obtained by fine mesh calculations. By projection using the POD bases, the dimension of the target equation (e.g., discretized neutron diffusion equation) can be dramatically reduced. In the proposed method, POD is effectively applied to each single assembly calculation (local calculation). Furthermore, using the local/global iteration, the effective neutron multiplication factor and the neutron flux distribution in the whole core geometry can be obtained by combining the numerical results of the local calculation for each fuel assembly and the global calculation for the whole core. As a feasibility study, the proposed method is applied to a one-dimensional heterogeneous core analysis, and the accuracy is investigated by changing the total number of POD bases. (authors)