Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nuclear data uncertainties”

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.

At least 1,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Deuteron-induced reactions on natural Zr from threshold to 50 MeV: production of 86g Y

Two stacks of thin Zr foils were irradiated with 30 and 50 MeV deuterons, respectively, using the Lawrence Berkeley National Laboratory 88-Inch Cyclotron, and 19 excitation functions for nat Zr(d,x) reactions were measured over a beam energy range of 6.3–47.64 MeV, where the independent cross sections for nat Zr(d,x) 88 Nb and nat Zr(d,x) 86m,g Y were measured for the first time. The well-characterized nat Fe(d,x) 56 Co, nat Ni(d,x) 56 Co, nat Ni(d,x) 58 Co, nat Ni(d,x) 61 Cu, nat Ti(d,x) 46 Sc and nat Ti(d,x) 48 V monitor reactions were used to determine the deuteron beam current throughout the stacks. All cross sections were determined using High Purity Germanium (HPGe) detector γ-ray spectroscopy. A variance minimization technique was employed to simultaneously constrain the deuteron beam currents with multiple monitor reactions, thus reducing systematic uncertainties. An additional 16 channels are reported for reactions on the nickel, titanium, and iron monitor foils, leading to a total of 35 excitation functions, with seven reaction channels reported for the first time in this work. The measured excitation functions are compared to calculations provided by the reaction modeling codes TALYS – 2.0, ALICE – 2020, CoH – 3.5.3 and EMPIRE – 3.2.3, as well as the TENDL – 2023 data library. The degree of agreement between theory and experiments is discussed. The possible production of the important PET radionuclide 86g Y via the nat Zr(d,x) route was critically examined. The physical yields for nat Zr(d,x) 86 Y and other yttrium isotopes produced were calculated and compared to other production pathways. Due to high-level of associated radionuclide impurities, this route cannot deliver 86g Y suitable for medical applications.

86gY↗
DOE OSTI2026

Bayesian discovery of optimal reduced order models from mechanistic and experimental data: A case study of Pd penetration in TRISO fuels using BISON

TRistructural ISOtropic (TRISO) particles rely on a silicon carbide (SiC) layer as the primary structural material and barrier to metallic fission products (FPs) release. Accurate prediction of palladium (Pd) transport and penetration is therefore critical for qualifying TRISO fuels for advanced reactors. The empirical correlation for Pd penetration in BISON is derived from historical particle-fuel data, but cannot explain the large scatter in the experimental data that arises from varying experimental conditions. To aid fuel qualification, we previously developed a mechanistic reduced order model (ROM) using BISON that resolves these dependencies. Here, in this work we build on that mechanistic ROM and perform validation and quantify its uncertainty using Bayesian uncertainty quantification (UQ). calibration against a suite of in-pile and out-of-pile experiments spanning particle compositions, geometries, and operating conditions, and we benchmark it against the empirical correlation. Bayesian UQ identifies influential parameters, calibrates them to data, and yields predictive intervals. Results show that while the empirical correlation can be tuned to fit a single experiment type, it transfers poorly; the mechanistic ROM sustains accuracy with credible uncertainty across disparate conditions. This demonstrates a practical path—via Bayesian UQ applied to mechanistic ROMs—to leverage single-effect experiments for inferring in-reactor behavior and supporting TRISO fuel qualification.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Measurement of the 239 Pu(n,f)/ 235 U(n,f) Cross-Section Ratio with the NIFFTE fission Time Projection Chamber

In this work, the 239 Pu(n,f)/ 235 U(n,f) cross-section ratio has been measured with the fission Time Projection Chamber (fissionTPC) from 100 keV to 100 MeV. The fissionTPC provides three-dimensional reconstruction of fission-fragment ionization profiles, allowing for a precise quantification of measurement uncertainties. The measurement was performed at the Los Alamos Neutron Science Center which provides a pulsed white source of neutrons. The data are recommended to be used as a cross-section ratio shape. A discussion of the status of the absolute normalization and comparisons to ENDF evaluations and previous measurements is included.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Spatially resolved properties of early-type group-dominant galaxies with MUSE: gas content, ionization mechanisms, and metallicity gradients

ABSTRACT With the goal of a thorough investigation of the ionized gas and its origin in early-type group-dominant galaxies, we present archival MUSE data for 18 galaxies from the Complete Local-Volume Groups Sample (CLoGS). These data allowed us to study the spatially resolved warm gas properties, including the morphology of the ionized gas, EW(H α), and kinematics as well as the gas-phase metallicity (12 + log(O/H)) of these systems. In order to distinguish between different ionization mechanisms, we used the emission-line ratios [O iii]/H β and [N ii]/H α in the BPT diagrams and EW(H α). We find that the ionization sources in our sample have variable impacts at different radii; central regions are more influenced by low-luminosity active galactic nuclei, while extended regions of low-ionization nuclear emission-line region-like emission are ionized by other mechanisms with post-asymptotic giant branch stars photoionization likely contributing significantly. We classified our sample into three H α +[N ii] emission morphology types. We calculate the gas-phase metallicity assuming several methods and ionization sources. In general, 12 + log(O/H) decreases with radius from the centre for all galaxies, independently of nebular morphology type, indicating a metallicity gradient in the abundance profiles. Interestingly, the more extended filamentary structures and all extranuclear star-forming regions present shallow metallicity gradients. Within the uncertainties these extended structures can be considered chemically homogeneous. We suggest that group-dominant galaxies in our sample likely acquired their cold gas in the past as a consequence of one or more mechanisms, e.g. gas-clouds or satellite mergers/accretion and/or cooling flows that contribute to the growth of the ionized gas structures.

Lagos, P. (ORCID:0000000223218657)↗
DOE OSTI2020

VADER: A Tool for Criticality Safety Validation

The purpose of criticality safety is to prevent any inadvertent criticality from occurring during the handling or storage of fissile material. Calculations are frequently used to demonstrate that a sufficient subcritical margin exists. Validation is a key aspect of the evaluation process, establishing the suitability, accuracy, and associated uncertainty of the computational method and data to be used for the intended application. The validation process is performed by comparing the results of critical experiments with the calculated results from models of the experiments using the computational method to be validated. Laboratory critical experiments are controlled systems that achieve a k eff of approximately 1 in order to investigate the parameters at which such a critical condition is achieved. The validation parameters that are traditionally applied to safety analysis calculations are the bias and the bias uncertainty . The bias is the deviation of the average k eff of the validation suite from unity. The bias uncertainty accounts for the statistical uncertainty in the bias based on the standard deviation, sample size, and distribution of k eff values of the validation suite. The values of bias and bias uncertainty ensure that the systems predicted to be subcritical by the computational method will indeed be subcritical. The bias and bias uncertainty are often combined to determine an upper subcritical limit (USL) or computational margin that can then be applied to safety analysis calculations. Many methods have been developed by different organizations to calculate the bias and bias uncertainty for various types of criticality analyses. Each of these methods typically requires that the validity of various underpinning statistical assumptions be confirmed to demonstrate that the method is appropriate for the analysis of a given validation suite. An example of the validation decision making flow is shown in Fig.1. As shown in Fig. 1, the analyst performing the validation fits a trend line to the data and performs a test to determine if the trend was a statistically better representation of the data than if it were treated as an uncorrelated sample. If the trend line is a better representation of the data, then the analyst uses any one of a number of trending techniques to determine the bias and bias uncertainty. If a trend is not an appropriate representation of the data, then the analyst proceeds to perform a normality assessment for the data. If the normal assumption can be shown to be acceptable, then the analyst calculates the bias and bias uncertainty with the parametric technique. If the assumption of normality cannot be justified, then the nonparametric technique is used. Once the decision flow has been followed and the appropriate technique has been selected, the bias and bias uncertainty is typically combined with an administrative margin to determine a USL below which calculated values of k eff for safety analysis models can be considered subcritical. The calculations used in each decision are often performed with spreadsheets or with small programs available at various sites performing criticality analyses. Expertise in understanding and interpreting the results must be maintained to perform these calculations. This can often be an error-prone process. Oak Ridge National Laboratory (ORNL) is currently developing the Validation and Data Evaluation Resource (VADER) to simplify and automate the criticality safety validation process and to provide a software quality assurance pedigree to the calculational methods used. This paper discusses the use of the Fulcrum user interface with VADER, the anticipated initial capabilities of VADER to perform validation analyses, and the output from the code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Using Calibrated Water Data for Preliminary Validation of the SRT Code for Advanced Reactors

Recent interest and corresponding progress worldwide regarding advanced nuclear reactors has renewed focus on their performance and related safety assessments. Specifically, the U.S. Nuclear Regulatory Commission has emphasized the importance of mechanistic approaches to source term analysis for advanced reactor licensing applications, which attempt to realistically account radionuclide transport and retention phenomena. Further model development is required due to the numerous and complex physical and chemical phenomena associated with mechanistic source term analyses. Reflecting the need for modeling advancement, Argonne National Laboratory developed a mechanistic source term analysis tool for sodium fast reactors. The Simplified Radionuclide Transport (SRT) code describes fuel pin failure (for simulating the initial condition at the point of fuel pin breach), bubble scrubbing, deposition, leakage and following environmental impact. In the current work, a validation study of the SRT bubble scrubbing model is performed using a water-loop experiment performed at the University of Wisconsin-Madison. Through the analysis, the approach and fundamental bubble scrubbing models in SRT, which examine the removal of aerosols within the bubble as it is transported through a pool, have been widely evaluated. The results of the assessment demonstrate a high level of agreement in the regions of greater aerosol size. In the parameter range of minimal aerosol removal, the simulation slightly underpredicts the experiment results; however, considering the scale of plots and huge uncertainties inherently included, the deviation can be judged to be minor and would produce a conservative result. In addition, uncertainty analysis has been further refined to reflect the experimental distribution of parameters including aerosol sizes, which induces a span of performance for each representative aerosol size. Based upon the initial validation results along with uncertainty effects, SRT is expected to provide meaningful insights for the analysis of bubble scrubbing. Future sodium-loop tests will provide further validation basis.

Kam, Dong Hoon↗
DOE OSTI2020

Observation of a new baryon state in the $Λ^0_bπ^+π^-$ mass spectrum

A new baryon state is observed in the $Λ^0_bπ^+π^-$ mass spectrum with high significance using a data sample of pp collisions, collected with the LHCb detector at centre-of-mass energies $\sqrt{s}$ = 7, 8 and 13 TeV, corresponding to an integrated luminosity of 9 fb -1 . The mass and natural width of the new state are measured to be m = 6072.3 ± 2.9 ± 0.6 ± 0.2 MeV , Γ = 72 ± 11 ± 2 MeV , where the first uncertainty is statistical and the second systematic. The third uncertainty for the mass is due to imprecise knowledge of the $Λ^0_b$ baryon mass. The new state is consistent with the first radial excitation of the $Λ^0_b$ baryon, the Λ b (2S) 0 resonance. Updated measurements of the masses and the upper limits on the natural widths of the previously observed Λ b (5912) 0 and Λ b (5920) 0 states are also reported.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Model 9975 O Ring Fixture LongTerm Leak Performance (FY2021 Status Report)

Leak testing experiments to monitor the aging performance of Viton® GLT and GLT-S O-rings used in the model 9975 shipping package has been ongoing since 2004 at Savannah River National Laboratory. Seventy tests using mock-up 9975 primary containment vessels (PCVs) with GLT O-rings were assembled and heated to temperatures ranging from 200 to 450 °F. Due to material substitution, fourteen tests with GLT-S O-rings were initiated in 2008 and heated to temperatures ranging from 200 to 400 °F. The conditioning temperatures are elevated compared to the calculated maximum O-ring temperature in a 9975 package in storage, 158 °F, to accelerate aging and enable observations of O-ring failures in a reasonable time frame. The fixtures are leak tested periodically and all GLT O-ring fixtures aging at 350 °F or above have failed to maintain a leak-tight seal. Eight GLT O-ring fixtures aging at 300 °F have failed after 2.8 to 5.7 years at temperature while the remaining fixtures at 300 °F were retired from testing following more than five years without failure. There has been one GLT O-ring fixture which failed after 13.4 years aging at 200 °F. However, 19 other GLT O-rings aging at 200 °F have remained leak-tight for over 14 years and remain in test. No failures have yet been observed in GLT O-ring fixtures aging at 270 °F for 9.5 years. All GLT-S O-ring fixtures aging at 300 °F or above have failed their leak test. No failures have yet been observed in GLT-S O-ring fixtures aging at 200 and 250 °F for 11.5 years. Data from the O-ring fixtures are generally consistent with results from compression stress-relaxation testing, and provide confidence in the predictive models based on those results. However, uncertainty exists in extrapolating these elevated temperature results to the lower temperatures of interest for normal storage in K-Area Complex (KAC). Oxygen consumption testing, which includes results from temperatures near KAC normal storage temperatures, is ongoing to provide further confidence and corroborate these extrapolations. The collective data from these test efforts suggest the minimum O-ring service life at KAC normal storage conditions should be at least 34 years for GLT and GLT-S O-rings. Measurement of compression set in O-rings removed from failed fixtures, compared to that from KAC surveillance O-rings, indicate significant margin remains for O-rings still in service in 9975 packages in KAC. Aging and periodic leak testing will continue for the remaining 26 mock-up PCV fixtures.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

STOCHASTIC OPTIMIZATION FOR LONG TERM CAPITAL STRUCTURES, SYSTEMS, AND COMPONENTS REFURBISHMENT AND REPLACEMENT

As commercial nuclear power plants (NPPs) pursue extended plant operations in the form of Second License Renewals (SLRs), opportunities exist for these plants to provide capital investments to ensure long-term, safe, and economic performance. Several utilities have already announced their intention to pursue extended operations for one or more of their NPPs via SLR2. The goal of this research is to develop a riskinformed approach to evaluate and prioritize plant capital investments made in preparation for, and during the period of, extended plant operations to support decisions in NPP operations. In order to prioritize project selection via a riskinformed approach we developed a single decision-making tool that integrates safety/reliability, cost, and stochastic optimization models to provide users with data analysis capabilities to more cost effectively manage plant assets. Both stochastic analysis methods—such as Monte Carlo-based sampling strategies—and multi-stage stochastic optimization strategies are employed to provide priority lists to decisionmakers in support of risk-informed decisions. We applied the proposed method to a trial application of projected replacement/refurbishment expenditures for plant capital assets (i.e., Structures, Systems, and Components [SSCs]). The objective is to optimize the SSC replacement/refurbishment schedule in terms of economic constraints, data uncertainties, and SSC reliability data, as well to generate a priority list for maximizing returns on investment.

42 ENGINEERING↗
DOE OSTI2024

The Beta-decay Paul Trap Mk IV: Design and commissioning

Here, the Beta-decay Paul Trap is an open-geometry, linear trap used to measure the decays of 8 Li and 8 B to search for a tensor contribution to the weak interaction. In the latest 8 Li measurement of Burkey et al. (2022), β scattering was the dominant experimental systematic uncertainty. The Beta-decay Paul Trap Mk IV reduces the prevalence of β scattering by a factor of 4 through a redesigned electrode geometry and the use of glassy carbon and graphite as electrode materials. The trap has been constructed and successfully commissioned with 8 Li in a new data campaign that collected 2.6 million triple coincidence events, an increase in statistics by 30% with 4 times less β scattering compared to the previous 8 Li data set.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

DUNE and MINERvA Flux Studies and a Measurement of the Charged-Current Quasielastic Antineutrino Scattering Cross Section with $<E_{\nu}> \sim 6 $ GeV on a CH Target

One of the biggest open questions in physics is the prevalence of matter over anti-matter in the present universe. One of possible answers lies in the violation of charge and parity symmetry in the lepton sector that would favor matter over anti-matter so that the universe becomes dominated by matter over time. Recent results from neutrino experiments like T2K and NOvA have indicated that there is CP violation in the leptonic sector \citep{Abe:2019vii}. However, the precise measurement of CP violation will require an unprecedented level of accuracy in these oscillation experiments. One of the challenges lies with the systematic uncertainties that come with the measurements, which include uncertainties related to the neutrino flux modeling, interaction modeling and detector related systematic uncertainties.\\ The first part of this thesis will go through the challenges related to flux uncertainties. It will discuss the flux uncertainties related to hadron production in the LBN F beamli ne for the DUNE experiment. Hadron production models used in flux simulations vary widely depending on the simulation and choice of physics models used to simulate the hadron production. This thesis will explain the method of using existing hadron production data to constrain these uncertainties, a method which has been used in MINERvA experiment.\\ The second part of the thesis goes through the effect of possible flux mismodeling and MINERvA's approach to address this effect. In doing so, we discovered that the suspected mismodeling could also result from incorrect estimation of the energy scale of the muons. This project not only demonstrated the challenges of flux modeling but also led to a novel approach of using neutrino energy spectra to understand the correlation between detector and neutrino beamline parameters. \\ A high statistics measurement of the anti-neutrino scattering cross section in the charged current quasielastic (CCQE) channel is the final part of the th esis. Th e 2-body interactions in this channel have the advantage of allowing reconstruction of the interaction kinematics from the outgoing muon trajectory which can be measured very well. Double differential cross sections as a function of muon kinematics, one of the deliverables of this analysis, will help future oscillation experiments understand their data. The cross-section as a function of the four-momentum transferred from the leptonic system to the hadronic system ($Q^{2}_{QE}$) can be used to test models used in simulating antineutrino interactions and the nuclear effects that can modify the predicted cross sections. Nuclear effects arise from the complex nuclear environment and both modify the initial scattering and change the fate of the particles produced from in neutrino-nucleus interactions. Study of nuclear effects will help to understand the structure of the atomic nucleus and its impact on neutrino interactions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

A Technical and Economic Assessment of LWR Flexible Operation for Generation and Demand Balancing to Optimize Plant Revenue

With increased penetration of subsidized variable renewable energy (VRE) resources and competition from low natural gas prices, existing light water reactor (LWR) nuclear power plants (NPPs) are struggling to remain economically competitive. This work examines the potential economic competitiveness of various thermal energy storage (TES) technologies when coupled directly or indirectly with a NPP. To highlight their relative economic competitiveness, we contrast several energy storage solutions in stochastic dispatch optimization. We leverage data from recent work analyzing a range of TES technologies with varying capital costs, performance, and technology readiness level (TRL) to establish our case. We explore inserting these technologies into an electricity market with existing nuclear generation and large projected variable renewable energy (VRE) penetration. Although these technologies' projected capital costs may make them unlikely candidates in their current state, this analysis demonstrates a high-fidelity techno-economic analysis of energy storage. Furthermore, as the projected cost of energy storage technologies evolves, this analysis sets a precedent for similar future investigations. One region with projected trends that may be unfavorable for existing nuclear capacity is the New York Independent System Operator (NYISO) market. New York state’s baseload generation has been historically provided by fossil-fired and nuclear assets. However, amid economic pressures from subsidized VREs and low natural gas prices, the state has recently deactivated Indian Point nuclear power plant units 2 and 3. Furthermore, the state plans to meet its zero-emission generation target by 2040 by replacing fossil-fired capacity with significant investments in VRE resources like wind and solar photovoltaic (PV) and battery storage. Increased intermittent resource penetration lowers the baseload power requirement, adding further economic pressure to the state’s three remaining NPPs still in operation. With three NPPs still in operation in New York, this work analyzes potential economic benefits to NPPs on the New York grid when directly or indirectly coupled with various TES technologies. This work requires two modeling steps to analyze the potential economic benefits of various system configurations of the TES directly or indirectly coupled with nuclear. First, this analysis leverages capacity expansion modeling by experts at the Electric Power Research Institute (EPRI). Using their deterministic capacity expansion model, U.S. Regional Economy, Greenhouse Gas, and Energy (US-REGEN), EPRI analysts evaluated the capacity and generation evolution of the New York state energy market under four projection scenarios. These four projection scenarios were developed to represent the potential evolution of the capacity and generation in NYISO from 2015 to 2050 under various economic, technology, and policy constraints. The results from these capacity expansion models are then used as boundary conditions in the second modeling step. The second modeling step uses the Holistic Energy Resource Optimization Network (HERON) for a set of stochastic techno-economic analyses (STEAs) to investigate the potential increase in the economic viability of various configurations of the TES. With no current capacity expansion capabilities, HERON takes the data generated from US-REGEN for 2050 to generate synthetic load, solar, and wind data. Then HERON economically optimizes the capacity and dispatch of the various TES configurations. The potential economic benefit is the differential net present value (NPV) of the TES configurations from the no-TES baseline. As a stochastic techno-economic analysis package, HERON introduces uncertainty into the economic metrics, while US-REGEN trades resolution for reduced computational complexity. Using HERON also allows the modeling of direct thermal coupling, a feature not common in capacity and dispatch models. As expected, with high capital costs, the costs of introducing energy storage for all the technologies considered outweighed the potential economic benefit of this strategy for flexible plant operation. The benefit of this analysis is primarily in demonstrating a workflow that examines innovative solutions to increase NPP revenue via TES coupling. HERON’s stochastic capacity and dispatch optimization process used in this work has proven an effective tool in observing and evaluating the impact of introducing storage technologies in a grid energy system.

25 ENERGY STORAGE↗
DOE OSTI2023

High-resolution absorption spectroscopy of room-temperature and jet-cooled ammonia between 59,000 and 93,000 cm -1

Here, we present new high-resolution photoabsorption spectra of ammonia spanning the region between 59,000 cm -1 and 93,000 cm -1 that were recorded by using the Fourier Transform Spectrometer at the Synchrotron SOLEIL. This region extends from just above the Franck-Condon envelope for the $\tilde{A}$ 1 $A$ 2 " ← $\tilde{X}$ 1 $A$ 1 ' transition to well above the NH 3 + $\tilde{X}$ + 2 $A$ 2 " ionization threshold. The spectra were recorded at a measured resolution of 0.23 cm -1 in both a room-temperature cell (293 K) and in a slit-jet supersonic expansion (~70 K). The absolute photoabsorption cross section with an uncertainty of ±5% is also reported for the room-temperature spectrum. The present resolution is a factor of 10 – 100 times higher than in other recently reported broad band spectra of ammonia, and many of the observed bands show partially resolved rotational structure. We have attempted to assign this structure for a number of these bands. The oscillator strengths extracted from the data are in good agreement with previous measurements but, in the case of structured bands, the present higher resolution measurements return higher peak absorption cross sections, that increase further when the sample is cooled. The present higher resolution spectra suggest that a number of previous vibronic band assignments that were based on quantum defect considerations may require some revision. Finally, we discuss the substantial differences between the photoabsorption and photoionization data just above the first ionization threshold.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

A Plutonium Mass Uncertainty Assessment Using a Cherenkov-Based Neutron Multiplicity Water Detector

We present a new analysis of the sensitivity of a water Cherenkov-based neutron multiplicity counter for nondestructive assay of trace quantities of 240 Pu. The multiplicity counter [water-based neutron detector (WaND)] contains 1.0 m 3 of water doped with 0.2% gadolinium in the form of a gadolinium salt (0.4% GdCl 3 ). It can be used for testing the presence of trace quantities of 240 Pu within a sample of unknown material by observing the correlated timing behavior of neutron capture candidates in the detector. The new analysis takes advantage of the fact that the detector is comprised almost entirely of a homogeneous medium, which permits the use of a simple model to reproduce the observed neutron capture timing behavior. Finally, the results of this analysis indicate that the detector has an absolute neutron detection efficiency of 23.4 ± 1.1% and is sensitive to a 240 Pu mass as small as 2.3 mg (at 1σ) after a 1-h data acquisition.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2023

Fuel performance uncertainty quantification and sensitivity analysis in the presence of epistemic and aleatoric sources of uncertainties

Fuel performance modeling and simulation includes many uncertain parameters from models to boundary conditions, manufacturing parameters and material properties. These parameters exhibit large uncertainties and can have an epistemic or aleatoric nature, something that renders fuel performance code-to-code and code-to-measurements comparisons for complex phenomena such as the pellet cladding mechanical interaction (PCMI) very challenging. Additionally, PCMI and other complex phenomena found in fuel performance modeling and simulation induce strong discontinuities and non-linearities that can render difficult to extract meaningful conclusions form uncertainty quantification (UQ) and sensitivity analysis (SA) studies. In this work, we develop and apply a consistent treatment of epistemic and aleatoric uncertainties for both UQ and SA in fuel performance calculations and use historical benchmark-quality measurement data to demonstrate it. More specifically, the developed methodology is applied to the OECD/NEA Multi-physics Pellet Cladding Mechanical Interaction Validation benchmark. A cold ramp test leading to PCMI is modeled. Two measured quantities of interest are considered: the cladding axial elongation during the irradiations and the cladding outer diameter after the cold ramp. The fuel performance code used to perform the simulation is FAST. The developed methodology involves various steps including a Morris screening to decrease the number of uncertain inputs, a nested loop approach for propagating the epistemic and aleatoric sources of uncertainties, and a global SA using Sobol indices. The obtained results indicate that the fuel and cladding thermal conductivities as well as the cladding outer diameter uncertainties are the three inputs having the largest impact on the measured quantities. More importantly, it was found that the epistemic uncertainties can have a significant impact on the measured quantities and can affect the outcome of the global sensitivity analysis.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

DECOVALEX-2023: Task F2 Salt Final Report

The subject of Task F of DECOVALEX-2023 concerns performance assessment modelling of radioactive waste disposal in deep mined repositories. The primary objectives of Task F are to build confidence in the models, methods, and software used for performance assessment (PA) of deep geologic nuclear waste repositories, and/or to bring to the fore additional research and development needed to improve PA methodologies. In Task F2- (salt), these objectives have been accomplished through staged development and comparison of the models and methods used by participating teams in their PA frameworks. Coupled-process submodels and deterministic simulations of the entire PA model for a reference scenario for waste disposal in domal salt have been conducted. The task specification has been updated continuously since the initiation of the project to reflect the staged development of the conceptual repository model and performance metrics. Thermal, hydrological, mechanical, and chemical properties of individual components of the engineered and natural system were chosen for relevance by participating teams. The salt reference case system was characterized using data and measurements collected at relevant underground research laboratories (URLs), field sites, and simulation results from teams with specialized modelling capability. Participating teams made a wide range of model assumptions from compartmentalized networks to full 3D models of the salt formation. No single contributed model includes full-fidelity representation of all the features, events, and processes (FEPs) detailed in the task specification, but almost all features and processes are represented in at least one model. Despite differences in the modelling strategies developed by participating teams, all models indicate that salt compaction and radionuclide diffusion are key processes in the repository, and for the FEPs and model scenario considered, little of the disposed radionuclides will migrate beyond the repository seal over the 100,000 year simulations. In general, the model output quantities have the largest differences over the short term and near the waste. The models tend to be more similar further from waste and at later time. Disparities between the models are believed to be due to differing simplifications from the task specification, some of which are chosen simplifications to reduce complexity, and some are restrictions imposed by the modelling tools. A second round of this task has been accepted for DECOVALEX-2027 in conjunction with Task F1 on crystalline PA modelling. The future round includes waste package heating, improved modelling of salt creep closure, additional comparisons of coupled-process sub-models, and the impact of repository engineering design on radionuclide migration in the repository. Participants will also propose and finalize a set of uncertain inputs for the reference case simulations, propagate these uncertainties in a set of realizations, and conduct sensitivity analyses on the simulation results.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Measurement of material isotopics and atom number ratio with α-particle spectroscopy for a NIFFTE fission Time Projection Chamber actinide target

In this work, we present the results of a measurement of isotopic concentrations and atomic number ratio of a double-sided actinide target using α-spectroscopy and mass spectrometry. The double-sided actinide target, with predominantly 239 Pu on one side and 235 U on the other, was used in the fission Time Projection Chamber (fissionTPC) for a measurement of the neutron-induced fission cross-section ratio between the two isotopes. The measured atomic number ratio is needed to extract an absolute measurement fission cross-section ratio. The 239 Pu/ 235 U atom number ratio was measured with a combination of mass spectrometry and α-spectroscopy with a planar silicon detector achieving uncertainties of less than 1%. Different strategies for estimating isotopic concentration from the α-spectrum are presented to demonstrate the potential of these methods for non-destructive target assay. We found that a combination of fitting spectra with constraints from mass spectrometry, and summing counts in a region of the spectrum provided the most consistent results with the lowest uncertainty.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Integration and Demonstration of Monitoring, Modeling, and Prediction of DV-1 Amendment Performance at the Bench Scale: DV-1 Amendment Demonstration

During fiscal years 2024 and 2025, the U.S. Department of Energy’s Hanford Field Office commissioned Pacific Northwest National Laboratory to conduct applied research aimed at reducing the cost, time, and uncertainty associated with in situ treatment of vadose zone contaminants at the Hanford Site. This report outlines the integration of three key research efforts into a meso-scale demonstration designed to advance field-scale solutions that aim to (1) optimize the delivery of chemical amendments to contaminated soils, (2) reduce uncertainty in amendment delivery performance assessment using advanced monitoring techniques, and (3) provide real-time insights into when and where amendment-induced precipitation reactions occur in the subsurface. To achieve these objectives, the tank-scale (~ 1 cubic meter) Geophysical Imaging of Flow and Transport (GIFT) system was developed. GIFT enables experimental testing of amendment delivery while incorporating automated multi-modal monitoring approaches, including pressure measurements, direct fluid sampling, and remote time-lapse geophysical imaging. The data generated from these monitoring techniques will serve as inputs for a generative artificial-intelligence-driven digital twin – a numerical simulation model designed to honor observed data while quantifying uncertainty in simulation accuracy. Using this simulator, researchers will refine an amendment injection strategy to maximize delivery efficiency within a low-permeability soil zone. Monitoring data will be interpreted through simulated outputs to enhance understanding of the injection process. The efficacy of this integrated approach will be evaluated through direct sampling at the conclusion of the experiment.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗