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.
This work validates a method of estimating sensitivities that utilized the MCNP PERT card and tallies. The agreement of the PERT card and Brute Force results proves the validity of calculating nuclear data sensitivities by use of the MCNP PERT card, especially at higher energies. The overall agreement will be easier to understand once an uncertainty analysis is completed.
The calculation of the spectrum of antineutrinos ($\bar{v}_e$) from a reactor is a complicated problem requiring several nuclear data and physics inputs. Many of these have been investigated thoroughly to improve calculations and properly account for uncertainties. One input which has heretofore escaped consideration is the fission-yield distribution between ground and isomeric states. Here, in this work, we explore the effect of incorporating newly evaluated isomeric yield ratios (IYR) for 43 fission products into summation calculations and identify the disproportionate importance of certain isotopes, particularly at higher energies. Our analysis indicates that updated IYRs contribute to a significant increase in the $\bar{v}_e$ spectrum around and above 7 MeV, with increases of more than 50% at higher energies. Through a detailed sensitivity study, we highlight a number of isotopes for which the IYR has a substantial effect on the $\bar{v}_e$ spectrum. This work stresses the critical role of isomeric yields in calculations of reactor $\bar{v}_e$ spectra and points to the necessity for their accurate experimental determination, especially for key fission products, in order to refine our understanding and address the observed anomalies between measured and calculated $\bar{v}_e$ spectra.
The interplay of quantum anomalies with strong magnetic fields and vorticity in chiral systems could lead to novel transport phenomena, such as the chiral magnetic effect (CME), the chiral magnetic wave (CMW), and the chiral vortical effect (CVE). In high-energy nuclear collisions, these chiral effects may survive the expansion of a quark–gluon plasma fireball and be detected in experiments. The experimental searches for the CME, the CMW, and the CVE have aroused extensive interest over the past couple of decades. The main goal of this article is to review the latest experimental progress in the search for these novel chiral transport phenomena at the Relativistic Heavy Ion Collider at Brookhaven National Laboratory and the Large Hadron Collider at CERN. Future programs to help reduce uncertainties and facilitate the interpretation of the data are also discussed.
Forthcoming in the next decade of neutrino physics advancements is the construction of the Deep Underground Neutrino Experiment (DUNE), a long baseline neutrino oscillation experiment which has the potential to answer fundamental questions such as why there is more matter than antimatter in the universe, and ultimately why we exist at all. A prototype of this experiment, ProtoDUNE, sits in a charged particle test beam at the European Organization for Nuclear Research (CERN) and uses Liquid Argon time projection chambers to image interactions of these particles. ProtoDUNE Single-Phase (ProtoDUNESP) had positively charged hadrons as well as positrons in its initial 2018 beam data run. Because ProtoDUNE is a surface-based time projection chamber, it experiences a constant influx of cosmic rays. These cosmic rays ionize the liquid argon within the detector, resulting in an accumulation of slow moving positive ions towards the cathode. This local charge distribution results in electric field non uniformities, which can have a significant effect on the accuracy of the event reconstruction. This paper outlines a monte carlo simulation study which adapts methodology from the MicroBooNE experiment, another neutrino oscillation experiment, to derive electric field correction maps which can correct for distortions of the electric field within the detector. These correction maps accurately estimate the electric field throughout the detector volume and correct for uncertainties in reconstructed event data enabling for physics discoveries.
Using a simple geometric framework with a realistic nuclear density distribution, we fit published HERMES data to determine fundamental properties of hadronization using the nuclear medium as a spatial analyzer. Our approach uses a fit to the transverse momentum broadening observable and the hadronic multiplicity ratio; the simultaneous fit to two different observables strongly constrains the outcome. Using the known sizes of the target nuclei, we extract the color lifetime, finding a z h -dependent range of values from 2 to 8 fm/c for these data. We also extract estimates for the qˆ transport coefficient characterizing the strength of the interaction between the quark and the cold nuclear medium, finding an average value of 0.035±0.011 GeV 2 /fm. With a three-parameter model we obtain satisfactory fits to the data with a goodness-of-fit parameter χ 2 /dof of 1.1 or less. In a secondary fit of the results from that model we independently find a value for the Lund String Model string tension of 1.00±0.05 GeV/fm. We evaluated the sensitivity for extracting quark energy loss and effective in-medium hadronic cross sections using four-parameter variants of the model, finding large uncertainties in both cases. Our results suggest that hadronic interaction of forming hadrons in the nuclear medium is the primary dynamical cause of meson attenuation in the HERMES data, with quark energy loss playing a more minor role.
When performing a validation exercise for determining the upper subcritical limit of a nuclear criticality safety application, an analyst should select and perform a statistical analysis on a population of benchmark experiments that are neutronically similar to the application. The size of this population should be sufficiently large such that the statistical analysis has a high degree of confidence that the bias plus bias uncertainty (calculational margin) has been accurately quantified. A complication arises because many benchmark experiments share common components, leading to correlations in their measured effective multiplication factors. Correlations between benchmark experiments within the population reduces its predictive power. This motivates the need for methods that consider benchmark experiment correlations and ensure adequate statistical significance of results. The Whisper code is a statistical analysis pack- age that incorporates nuclear data sensitivity coefficients from MCNP to assess benchmark experiment similarity and then performs an extreme-value analysis to estimate the bias plus bias uncertainty. The original methodology in Whisper does not consider the effect of benchmark experiment correlations when making this estimation, and this summary proposes the uniformly ordered binary decision algorithm to address this shortcoming. The original methodology in Whisper computes similarity coefficients ck for an application compared to all benchmark experiments in its library and develops weighting factors for a selected population proportional to the ck values. The methodology can be interpreted as statistically emulating a validation exercise for a particular application where the weighting factors may be viewed as the likelihood that an analyst would include a particular benchmark experiment within the population. The effective sample size of the population is the expected or mean number of benchmark experiments in the population. The uniformly ordered binary decision algorithm identifies clusters of correlated benchmark experiments within the population and then computes adjusted weighting factors based on the magnitude of the correlation coefficients within the cluster to compute a reduced effective sample size accounting for the lower information content because of correlations. Benchmark experiments within the cluster are ordered randomly with equal probability and probabilistic decisions are made as to whether a benchmark. experiment within the cluster should treated as redundant with a previous one; if two redundant benchmark experiments are included, then the conservative worst case bias plus bias uncertainty is used and the pair is counted as a single benchmark experiment in the population. Results are provided for HEU solutions in a research version of the Whisper software using benchmark experiment correlations provided by DICE, the Database for the International Criticality Safety Benchmark Evaluation Project (ICSBEP). These show that there can be a significant increase in the bias plus bias uncertainty because the effective sample size is reduced, and therefore the algorithm, needing to meet sample size requirements, expands the benchmark experiment population by accepting less similar benchmark experiments that would have otherwise not been included.
Gekko is an optimization suite in Python that solves optimization problems involving mixed-integer, nonlinear, and differential equations. The purpose of this study is to integrate common Machine Learning (ML) algorithms such as Gaussian Process Regression (GPR), support vector regression (SVR), and artificial neural network (ANN) models into Gekko to solve data based optimization problems. Uncertainty quantification (UQ) is used alongside ML for better decision making. These methods include ensemble methods, model-specific methods, conformal predictions, and the delta method. An optimization problem involving nuclear waste vitrification is presented to demonstrate the benefit of ML in this field. ML models are compared against the current partial quadratic mixture (PQM) model in an optimization problem in Gekko. GPR with conformal uncertainty was chosen as the best substitute model as it had a lower mean squared error of 0.0025 compared to 0.018 and more confidently predicted a higher waste loading of 37.5 wt% compared to 34 wt%. The example problem shows that these tools can be used in similar industry settings where easier use and better performance is needed over classical approaches. Future works with these tools include expanding them with other regression models and UQ methods, and exploration into other optimization problems or dynamic control.
239 Pu(n,f) nu-bar evaluated with CGMF from E inc = 0.1-20 MeV. Changes compared to VIII.0: Prior: CGMF model included via Kalman and sensitivities of CGMF model parameters to nu-bar. Evaluation technique: Kalman including correction for PPP. Experimental data: Nearly all data that Phil took into account (I rejected: Huanqiao, Johnstone, Leroy, Nesterov, Smirenkin); New UQ for all experimental data; Marini; No correlations between unc. of different exp., except for Cf-252(sf) nu-bar uncertainty cross-correlating all uncertainties.
239 Pu(n,f) nu-bar evaluated with CGMF from E inc = 0.1-20 MeV. Changes compared to VIII.0: Prior: CGMF model included via Kalman and sensitivities of CGMF model parameters to nu-bar. Evaluation technique: Kalman including correction for PPP. Experimental data: Nearly all data that Phil took into account (I rejected: Huanqiao, Johnstone, Leroy, Nesterov, Smirenkin); New UQ for all experimental data; Marini; No correlations between unc. of different exp., except for Cf-252(sf) nu-bar uncertainty cross-correlating all uncertainties.
Here, the presence of 240 Pu in nuclear fuels for reactors has resulted in high uncertainties in the results of reactor and nuclear transmutation calculations because of deficiencies in 240 Pu-related nuclear data. Specifically for the prompt fission neutron spectrum (PFNS) of 240 Pu, there is only one neutron-induced, (n,f), measurement at 0.85 MeV incident neutron energy and only one complete spontaneous fission, (sf), measurement. This limited availability of data does not sufficiently guide nuclear data evaluations of these quantities. Here we report on a measurement of both the 240 Pu(sf) and the 240 Pu(n,f) PFNS, both over the emitted neutron energy range of 0.79–10.0 MeV, and from incident neutron energies of 1.0–20.0 MeV for the (n,f) reaction. Measurements were made with a hemispherical array of liquid scintillators at the high-energy Los Alamos Neutron Science Center white neutron source at the Weapons Neutron Research facility as part of the joint LANL-LLNL Chi-Nu experimental campaign to measure actinide fission neutron spectra. These measurements are the first of their kind, and provide clear experimental evidence for second-chance fission, third chance fission, and pre-equilibrium neutron emission processes in neutron-induced fission of 240 Pu, and are the first ever measurements above 1 MeV incident neutron energy.
The presence of 240 Pu in nuclear fuels for reactors has resulted in high uncertainties in the results of reactor and nuclear transmutation calculations because of deficiencies in 240 Pu-related nuclear data. Specifically for the prompt fission neutron spectrum (PFNS) of 240 Pu, there is only one neutron-induced, ($n,f$), measurement at 0.85 MeV incident neutron energy and only one complete spontaneous fission, ($sf$), measurement. This limited availability of data does not sufficiently guide nuclear data evaluations of these quantities. Here we report on a measurement of both the 240 Pu($sf$) and the 240 Pu($n,f$) PFNS, both over the emitted neutron energy range of 0.79–10.0 MeV, and from incident neutron energies of 1.0–20.0 MeV for the ($n,f$) reaction. Measurements were made with a hemispherical array of liquid scintillators at the high-energy Los Alamos Neutron Science Center white neutron source at the Weapons Neutron Research facility as part of the joint LANL-LLNL Chi-Nu experimental campaign to measure actinide fission neutron spectra. These measurements are the first of their kind, and provide clear experimental evidence for second-chance fission, third-chance fission, and pre-equilibrium neutron emission processes in neutron-induced fission of 240 Pu, and are the first ever measurements above 1 MeV incident neutron energy.
The calculation of the spectrum of antineutrinos ( $\overline{V}$ e ) from a reactor is a complicated problem requiring several nuclear data and physics inputs. Many of these have been investigated thoroughly to improve calculations and properly account for uncertainties. One input which has heretofore escaped consideration is the fission yield distribution between ground and isomeric states. In this work, we explore the effect of incorporating newly evaluated isomeric yield ratios (IYR) for 43 fission products into summation calculations and identify the disproportionate importance of certain isotopes, particularly at higher energies. Our analysis indicates that updated IYRs contribute to a significant increase in the $\overline{V}$ e spectrum around and above 7 MeV, with increases of more than 50% at higher energies. Through a detailed sensitivity study, we highlight a number of isotopes for which the IYR has a substantial effect on the $\overline{V}$ e spectrum. This work stresses the critical role of isomeric yields in calculations of reactor $\overline{V}$ e spectra and points to the necessity for their accurate experimental determination, especially for key fission products, in order to refine our understanding and address the observed anomalies between measured and calculated $\overline{V}$ e spectra.
Understanding the harmful effects of galactic cosmic rays (GCRs) on space exploration requires a substantial amount of nuclear data. Specifically, the interaction of energetic GCR charged particles with spacecraft materials generates secondary radiations that, through energy deposition, can harm astronauts and electronic systems. By identifying the gaps in our knowledge of the relevant nuclear data—such as interaction cross sections—and identifying ways to fill those gaps—with measurements, compilations, evaluations, dissemination, reaction modeling, sensitivity studies, and uncertainty quantification—the safety and viability of space exploration can be improved. This work surveys the state of the art in this interdisciplinary field and identifies promising collaborative research topics that have significant potential to advance our understanding of the effects of the space radiation environment on space exploration.
Recent observations of neutron stars with gravitational waves and X-ray timing provide unprecedented access to the equation of state (EoS) of cold dense matter at densities difficult to realize in terrestrial experiments. At the same time, predictions for the EoS equipped with reliable uncertainty estimates from chiral effective field theory ( x EFT) allow us to bound our theoretical ignorance. In this work, we analyze astrophysical data using a nonparametric representation of the neutron-star EoS conditioned on x EFT to directly constrain the underlying physical properties of the compact objects without introducing modeling systematics. We discuss how the data alone constrain the EoS at high densities when we condition on x EFT at low densities. Here, we also demonstrate how to exploit astrophysical data to directly test the predictions of x EFT for the EoS up to twice nuclear saturation density, in order to estimate the density at which these predictions might break down. We nd that the existence of massive pulsars, gravitational waves from GW170817, and NICER observations of PSR J0030+0451 favor x EFT predictions for the EoS up to nuclear saturation density over a more agnostic analysis by as much as a factor of 7 for the quantum Monte Carlo (QMC) calculations used in this work. While x EFT predictions using QMC are fully consistent with gravitational-wave data up to twice nuclear saturation density, NICER observations suggest that the EoS stiffens relative to these predictions at or slightly above nuclear saturation density. Additionally, for these QMC calculations, we marginalize over the uncertainty in the density at which x EFT begins to break down, constraining the radius of a 1.4M $\bigodot$ neutron star to R 1.4 = $11.40^{+1.38}_{–1.04}$ ($12.54^{+0.71}_{–0.63}$) km and the pressure at twice nuclear saturation density to p(2n sat ) = $14.2^{+18.1}_{–8.4}$ ($28.7^{+15.3}_{–15.0}$) MeV=fm 3 with massive pulsar and gravitational-wave (and NICER) data.
Critical heat flux (CHF) is a key quantity in nuclear system modeling due to its impact on heat transfer, safety margins, and reactor performance. This study develops and validates an uncertainty-aware hybrid modeling approach that combines machine learning with physics-based models to predict CHF in cases of dryout. The Biasi and Bowring empirical correlations were paired with three ML uncertainty quantification (UQ) techniques: deep neural network (DNN) ensembles, Bayesian neural networks (BNNs), and deep Gaussian processes (DGPs). A pure ML model without a base model was evaluated for comparison. Model performance was assessed under plentiful (7,350 points) and limited (9 points) training data scenarios using parity, uncertainty distributions, and calibration curves. Results show that the Biasi hybrid DNN ensemble achieved the best overall performance, with a mean absolute relative error of 1.846%, and well-calibrated uncertainty estimates. The BNN-based hybrids showed slightly higher error (2.14%) but superior uncertainty calibration. DGP models underperformed, with over 6% error and poor uncertainty calibration. All hybrid models outperformed pure machine learning configurations, demonstrating resistance against data scarcity. These findings indicate that hybrid modeling significantly improves predictive accuracy, interpretability, and resilience to data scarcity. The integration of uncertainty awareness provides actionable confidence in CHF predictions, which is vital for safety-critical decisions in nuclear applications. This hybrid approach offers a viable pathway for deploying ML models in reactor analysis tools while preserving domain knowledge and physical consistency.
Understanding the harmful effects of galactic cosmic rays (GCRs) on space exploration requires a substantial amount of nuclear data. Specifically, the interaction of energetic GCR charged particles with spacecraft materials generates secondary radiations that, through energy deposition, can harm astronauts and electronic systems. By identifying the gaps in our knowledge of the relevant nuclear data – interaction cross sections – and identifying ways to fill those gaps – with measurements, compilations, evaluations, disseminations, reaction modeling, sensitivity studies, and uncertainty quantification – the safety and viability of space exploration can be improved. This work surveys the state of the art in this interdisciplinary field and identifies promising collaborative research topics that have significant potential to advance our understanding of the effects of the space radiation environment on space exploration.
Here we present a method for measurement analyses based on probabilistic deep neural networks that provide several advantages over conventional analyses with phenomenological models. These include predicting physical quantities directly from data, the rapid generation of statistically robust uncertainties, and the ability to bypass some parameters that may induce ambiguities and complications in data analysis. As deep learning methods make predictions through “black boxes,” the uncertainty quantification is typically challenging. We use a probabilistic framework that provides thorough uncertainty quantification and is straightforward to follow in practice. With the network architecture based on the Transformer, we demonstrate the current method for predicting nuclear resonance parameters from scattering data using the phenomenological R-matrix model.
The SCALE code package offers several nuclear data libraries to support Monte Carlo (MC) transport, as well as MC-based derivation of $\kappa$ eff sensitivity and uncertainty (S/U) data. The CSAS sequence using the KENO MC code can utilize continuous-energy (CE) cross sections, or pre-generated multigroup (MG) cross section libraries. The use of MG libraries introduces bias into calculations in exchange for faster transport solutions. The TSUNAMI-3D sequence also utilizes KENO MC calculations. TSUNAMI-3D has two CE calculational methods: the Iterated Fission Probability (IFP) method, and the Contribution-Linked eigenvalue sensitivity/Uncertainty estimation via Tracklength importance CHaracterization (CLUTCH) method. Previous work has shown poor agreement between CLUTCH and confirmatory direct perturbation calculations in specific applications, e.g., fissionable and polyethylene reflectors