Search NASA⌕ Search

SEARCH · Search NASA

Results for “CAPTURE EFFECT”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

A Modified Network Alteration Model to Predict Quasi-Static Behavior of the Cross-Linked Polymers During Hydrolytic Aging

Abstract The network alteration theory is utilized to capture the effects of hydrolytic aging on the mechanical integrity of the polymer matrix. Hydrolytic degradation generally considered as a result of two simultaneous chemical phenomena; (i) reduction of the cross-links, and (ii) active chain loss. Assuming the alteration of the chain density is identical to the change of the number average molecular weight along the aging trajectory, the strain energy of polymer matrix can be rewritten as a function of deformation, deformation history, storage time, and aging temperature. The proposed model includes only five physically inspired material parameters. Thus, while it is computationally efficient, it shows good agreement with own experimental data, which performed on various range of accelerated aging temperatures and times. With respect to its computational efficiency, simplicity, accuracy, and interpret-ability, the model is a good choice for advanced implementations in FE programs.

Bahrololoumi, Amir↗
DOE OSTI2020

Determination of jet calibration and energy resolution in proton–proton collisions at $\sqrt{s} = 8~\hbox {TeV}$ using the ATLAS detector

The jet energy scale, jet energy resolution, and their systematic uncertainties are measured for jets reconstructed with the ATLAS detector in 2012 using proton–proton data produced at a centre-of-mass energy of 8 TeV with an integrated luminosity of 20fb -1 . Jets are reconstructed from clusters of energy depositions in the ATLAS calorimeters using the anti-k t algorithm. A jet calibration scheme is applied in multiple steps, each addressing specific effects including mitigation of contributions from additional proton–proton collisions, loss of energy in dead material, calorimeter non-compensation, angular biases and other global jet effects. The final calibration step uses several in situ techniques and corrects for residual effects not captured by the initial calibration. Furthermore, these analyses measure both the jet energy scale and resolution by exploiting the transverse momentum balance in γ + jet, Z + jet, dijet, and multijet events. A statistical combination of these measurements is performed. In the central detector region, the derived calibration has a precision better than 1% for jets with transverse momentum 150 GeV < p T < 1500 GeV, and the relative energy resolution is (8.4 ± 0.6)% for p T =100 GeV and (23 ± 2)% for p T =20 GeV. The calibration scheme for jets with radius parameter R = 1.0, for which jets receive a dedicated calibration of the jet mass, is also discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Learning reference governor for cycle-to-cycle combustion control with misfire avoidance in spark-ignition engines at high exhaust gas recirculation–diluted conditions

Cycle-to-cycle feedback control is employed to achieve optimal combustion phasing while maintaining high levels of exhaust gas recirculation by adjusting the spark advance and the exhaust gas recirculation valve position. The control development is based on a control-oriented model that captures the effects of throttle position, exhaust gas recirculation valve position, and spark timing on the combustion phasing. Under the assumption that in-cylinder pressure information is available, an adaptive extended Kalman filter approach is used to estimate the exhaust gas recirculation rate into the intake manifold based on combustion phasing measurements. The estimation algorithm is adaptive since the cycle-to-cycle combustion variability (output covariance) is not known a priori and changes with operating conditions. A linear quadratic regulator controller is designed to maintain optimal combustion phasing while maximizing exhaust gas recirculation levels during load transients coming from throttle tip-in and tip-out commands from the driver. During throttle tip-outs, however, a combination of a high exhaust gas recirculation rate and an overly advanced spark, product of the dynamic response of the system, generates a sequence of misfire events. In this work, an explicit reference governor is used as an add-on scheme to the closed-loop system in order to avoid the violation of the misfire limit. The reference governor is enhanced with model-free learning which enables it to avoid misfires after a learning phase. Experimental results are reported which illustrate the potential of the proposed control strategy for achieving an optimal combustion process during highly diluted conditions for improving fuel efficiency.

42 ENGINEERING↗
DOE OSTI2021

Evaluation of residual gas fraction estimation methods for cycle-to-cycle combustion variability analysis and modeling

Cycle-to-cycle combustion variability in spark-ignition engines during normal operation is mainly caused by random perturbations of the in-cylinder conditions such as the flow velocity field, homogeneity of the air-fuel distribution, spark energy discharge, and turbulence intensity of the flame front. Such perturbations translate into the variability of the energy released observed at the end of the combustion process. During normal operating conditions, the cycle-to-cycle variability (CCV) of the energy release behaves as random uncorrelated noise. However, during diluted combustion, in either the form of exhaust gas recirculation (EGR) or excess air (lean operation), the CCV tends to increase as dilution increases. Moreover, when the ignition limit is reached at high dilution levels, the combustion CCV is exacerbated by sporadic occurrences of incomplete combustion events, and the uncorrelation assumption no longer holds. The low or null energy released by partial burns and misfires has an impact on the following combustion event due to the residual gas that carries burned and unburned gases, which contributes to the deterministic coupling between engine cycles. Many residual gas fraction estimation methods, however, only address the nominal case where complete combustion occurs and combustion events are uncorrelated. Here we evaluate the efficacy of such methods on capturing the effects of partial burns and misfires on the residual gas estimate for high-EGR operation. The advantages and disadvantages of each method are discussed based on their ability to generate cycle-to-cycle estimates. Finally, a comparison between the different estimation techniques is presented based on their usefulness for control-oriented modeling.

42 ENGINEERING↗
DOE OSTI2024

Conceptual Spacer Design for the ATR GEN I Target for Pu-238 Production in the Advanced Test Reactor at Idaho National Laboratory

The initial target design used for Pu-238 production at Idaho National Laboratory was designed by Oak Ridge National Laboratory to optimize the production of Pu-238 in the High Flux Isotope Reactor (HFIR) and are referred to as HFIR GEN II targets. To take advantage of the Advanced Test Reactor’s (ATR) taller active core region a redesign of the HFIR GEN II targets was needed. It was proposed to stack two HFIR GEN II targets nose to nose about the core center line; however, this resulted in excessive neutron and photon heating in the pellets located in the center. This peak heating was not desirable so three alternative designs were investigated for the ATR GEN I targets. The python-based code, MCNP to ORIGEN2 in Python (MOPY), was used to calculate the heating rates after 40 days of irradiation to capture the effects of each configuration. The purpose of this paper is to document the details of these conceptual design calculations and comparisons for the ATR GEN I targets.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2025

Insights into Slip-Rate Time Functions, Rupture Parameter Correlations, and Ground Motions from Validated Multicycle Earthquake Ruptures

Earthquake strong-motion predictions using kinematic source modeling require knowledge of the slip-rate functions (SRFs) along the rupture and their distinct characteristics in asperities, background (off-asperity) areas and near the surface. Here, in this study we analyzed SRFs from well-validated, self-consistent, and fully dynamic rupture models from earthquake cycles obeying a rate-and-state friction law, from our companion study (Galvez et al., 2021). The shapes of SRFs in asperities are well described by the regularized Yoffe function (RYF), which has only two parameters: rise time T r and smoothing time T s , which control the generation of long- and short-period ground motions, respectively. In background areas, we demonstrate that, in addition to the primary rupture, multiple secondary ruptures may also nucleate from rupture heterogeneities related to asperities, resulting in SRFs with multiple peaks. Because it is impossible to fit a multiple-peak SRF by the single-peak RYF, we describe SRFs in background areas in an effective way by fitting their amplitude spectra with the RYF spectra. Such spectrally effective RYFs capture salient aspects of seismic-wave generation and can be used in rupture generators for strong motion prediction. We found that small T s values correlate with small characteristic weakening distances, large peak slip rates (PSRs), and large rupture velocities. T r values are larger in background areas and smaller in asperities. Within the shallow aseismic zone, T s values approximately quadruple whereas T r values approximately double. Because of this dominant T s increase, PSR values decrease in the near-surface zone. These features indicate that the generation of strong motions by the near-surface portions of the rupture is negligible in the studied scenarios.

Geosciences↗
DOE OSTI2020

Calculation of Friction Factors in Heat Pipes Using CFD in Support of the Sockeye Code

A heat pipe is two-phase heat transfer device which relies on surface tension and capillary pressure to provide a very efficient heat transfer mechanism. Currently, there is interest in designing micro nuclear power reactors using high temperature heat pipes to provide totally passive cooling. To support this, Sockeye is being developed as a heat pipe simulation capability. To support the development of Sockeye, individual phenomena must be tested. This work focuses on applying a CFD model to evaluate friction factors for use in Sockeye. A CFD model has been developed in Nek5000 which simulates each component fluid, i.e. sodium liquid and vapor, independently and captures the effects of the wick on friction factor. The model tests the effect of the wick permeability and the thickness of the wick occupied by each fluid component. The thickness of the wick occupied by each component can be directly related to the local vapor volume fraction from the Sockeye heat pipe model. Results from the CFD model indicate that friction factor in the liquid flow decreases with increasing wick permeability, while the friction factor in the vapor flow increases with increasing wick permeability. In the liquid, this is likely a direct result of the increased cross sectional area available for the flow. For the vapor, this is explained as the wick acting like a rough wall, while the increased area has only a negligible impact. Results for the tests of the thickness of the wick occupied by each fluid showed a thresholding behavior for both liquid and vapor. When the thickness of the wick occupied by each fluid was increased, it increased the effect of the wick, but only up to a certain point. Once a certain thickness was reached, the friction factors remained constant. When the fluid occupied only a very small fraction of the wick, friction factors returned to nominal values.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Simulating Boundary Region Plasma Instabilities

The boundary regions of low-temperature plasmas are known to be susceptible to kinetic instabilities, which can affect the energies and fluxes of particles directed at the material boundary. For example, both the ion acoustic instability as well as an instability near the electron plasma frequency have been observed. Particle-in-cell (PIC) simulation is a tool that, alongside experiments, can capture the effects these instabilities have on the particle distribution functions. Ultimately, simulations can determine under what conditions these effects are significant by comparing to theoretical predictions and explore conditions unamenable to experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A New 1-step Arrhenius Decomposition Model for PBX 9501: Calibration with Isochoric-Adiabatic High-Activation-Energy Asymptotics and Uncertainty Quantification

Here we present a new 1-step Arrhenius thermal decomposition model of PBX 9501. As opposed to previous studies, we use high-energy asymptotic theory, and previously published equations of state to parameterize the model to experimental critical-time data. The model parameters were fit to experimental data using a bootstrapping procedure to capture the effects of experimental uncertainty. Subsequent implementation of the model in a chemicalkinetics ODE solver shows good agreement of times to ignition vs initial temperature in isochoric thermal explosion calculations with the results of the high-energy asymptotic theory. We also compare the resultant ZND structures to those of one of our previously calibrated burn models.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2023

Development of HPC based phase field simulations tool for modification of alloy morphology to enhance material properties during additive manufacturing (AM) process

Currently, additive manufacturing (AM) technology is extensively used to manufacture aerospace components because of the unique advantages of the technique in producing custom designs, complex geometries in a near-net fashion that significantly reduces the manufacturing cost, and improved energy efficiency. However, the thermodynamic properties of existing aerospace alloys and the thermal conditions that prevail in most of the powder-bed based additive manufacturing processes lead to undesirable columnar / columnar-dendritic microstructures that promote solidification cracking during the process, and the development of anisotropy in the mechanical properties of the as-built components. Time consuming, energy-intensive post AM treatments are required to recover the mechanical properties of the wrought components for which the existing alloy compositions were designed for. However, novel alloy compositions, that exploit the AM thermal conditions to produce fine, equiaxed grain structures in the as-built condition, with mechanical properties equal to or exceeding those of the wrought work-horse alloy, Ti-6Al-4V has recently been demonstrated in Ti-Cu alloys. Significant energy savings up to 66% can be obtained for a typical aerospace component through digitally designing the AM process, process parameters, and alloy composition. The current project focused on developing a fundamental understanding of the evolution of the columnar to equiaxed transition (CET) occurring during solidification under AM thermal conditions in Ti-Cu and Ti-Cu-X alloys using high-fidelity phase field (PF) simulations using the leadership class computing facilities that exist in the national laboratories. The simulations were able to capture the effect of alloy and process conditions on CET, and clearly showed the beneficial effect of a ternary solute addition to Ti-Cu binary alloys on CET. The simulations indicated that bulk nucleation in the liquid ahead of the solidifying epitaxial front was promoted by a large, transient, thermal undercooling promoted by the mismatch between the growth rate of the dendrite tips in the epitaxial front and the imposed solidification rate. AM experiments performed at Raytheon Technologies Research Center using the same thermal conditions used in the PF simulations indicated the formation of equiaxed grains, although the grain size was bigger than the ones reported in the literature. On the other hand, the equiaxed grain size predicted by PF simulations were smaller than the grain size of Ti-Cu alloys reported in the literature. The industry will follow up the existing work in a future program that would involve further optimization of the process for the Ti-Cu-X alloys to demonstrate the application of the process and the alloy to a specific aerospace component.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Additive Manufacturing Flaw Assessment Methodology

An evolution fatigue data and flaw tolerance of components produced using the Powder Bed Fusion (PBF) Additive Manufacturing (AM) process is documented in this report. Initial differences in fatigue data for AM components compared to smooth bar fatigue data indicated a need for a more detailed analysis of AM data available in technical literature. The investigation was initiated to support the development of a fatigue analysis methodology for AM components to support the of codification of AM technology for pressure equipment. The project was initiated to collect and analyze stainless steel 304L and 316L AM fatigue test data and corresponding process and quality information to develop S-N and E-N based fatigue data representation. Additional AM fatigue test data including Inconel Ti-6-4 and aluminum alloys were also considered for comparison purposes Metallic AM parts tend to contain various forms of defects distributed throughout the part. If an AM part is subjected to fatigue loading in service, a fatigue analysis needs to be performed during the design process to ensure an acceptable service life for the part. Post-process machining and polishing do not to improve fatigue resistance in any significant degree. The low cycle fatigue regime is of particular interest to this project in support of flaw acceptance criteria currently under development by ASME’s BPTCS/BNCS. Internal defects become exposed as external surface defects during machining for the machining of the AM part to final dimensions. This implies that as long as inherent AM defects are within a controlled limit in terms of both size and distribution characteristics, the corresponding fatigue test data in terms of either S-N (stress life) or E-N (strain life) can be investigated and characterized to establish fatigue properties of AM parts for design and fatigue evaluation purposes. The resulting S-N or E-N curves and their scatter bands can be used to derive fatigue design allowable stress values by capturing the effects of distributed discontinuities within an acceptable limit.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Failure Analysis for Molten Salt Thermal Energy Storage Tanks for In-Service CSP Plants

Thermal Energy Storage (TES) is a fundamental component in concentrating solar power (CSP) plants to increase the plant's dispatchability, capacity factor, while reducing the levelized cost of electricity. In central receivers CSP plants, nitrate molten salts have been used for several years for operation temperatures of up to 565 degrees C. Despite many efforts to advance nitrate salt to higher operation temperatures (even considering a replacement with molten chloride salts) to achieve higher energy conversion efficiencies, the 565 degrees C temperature is currently considered the state-of-the art. Although molten salt tanks have been broadly deployed in commercial CSP plants worldwide, several failures have been reported in these tanks after a few months or years of operation, causing significant economic loss and mistrust in CSP technologies. Most of these failures are associated with the infancy of the technology and multiple issues related to tank design, fabrication, commissioning, and aggressive operation. A technical standard dedicated to the design and fabrication of molten nitrate TES tanks does not exist today. Current in-service molten salt tanks have been generally designed based on the American Petroleum Institute's (API) 650 and ASME Section II standards. The API 650 code provides guidelines for dimensions and fabrication for oil storage tanks up to 260 degrees C. The ASME standard provides allowable stress values for various materials at a range of temperatures and conditions. Both standards seem to be limited for molten salt TES tanks where high temperatures, thermal cycling, and transient conditions are expected. In 2020, NREL released the Concentrating Solar Power Best Practices Study (NREL/TP-5500-75763) that summarized multiple issues in CSP plants, along with potential alternatives and recommendations to address those issues based on information collected from participants representing about 80% of operating CSP plants in the world. One of the recommendations from this study was the development of accurate and validated models to evaluate the plant's transient operation, capable of capturing the effect of short-term clouds and operator response, while being flexible in being adapted to various spatial and temporal resource data. The "Failure Analysis for Molten Salt Thermal Energy Tanks for In-Service CSP Plants" project was inspired on this recommendation and was focused on (1) the development and validation of a physics-based model for a representative, commercial-scale molten salt tank, (2) performing simulations to evaluate the behavior of the tank as a function of typical plant operation conditions, (3) understanding tank failures mechanisms, (4) determining the residual stress and distortion in the tank floor after welding fabrication and evaluating their impact in the stresses developed in the tank during operation, (5) assessing the impact of key operation parameters on the temperature and stress distribution, (6) conduct a preliminary evaluation of design features to reduce stress and improve tank's reliability, and (7) estimate the tank's service life based on the stress developed under diverse operation scenarios. From the analysis conducted in the project and presented in this report, it was found that maximum stresses surpassing the yield strength point of the stainless steel (SS) 347H are developed on the tank floor near the perimeter. These large stresses are strongly influenced by the initial residual stresses and distortion of the tank floor after welding fabrication. During operation, large stresses are developed in the tank floor at high operation temperatures with large salt inventory levels during transient operation. High stresses are also related to elevated temperature gradients in the tank floor that could be attributed to insufficient mixing within the salt inflow and the salt inventory. Based on the analysis, creep is the predominant failure mechanism. However, the large stress levels could favor the plastic deformation into buckles, and crack formation due to stress relaxation cracking during cycle operation. A lifetime below 3 years was estimated for the typical plant operation conditions and a specific initial residual stress and deformation distribution of the tank floor. The estimated life agrees with the service time to failure reported in several commercial molten salt tanks. Desing and operation guidelines can be extracted from the analysis presented in this report, which could be adopted by tank manufacturers and CSP operators to advance toward an ultimate solution for tank failures by reducing residual and operational stresses to achieve a tank service life of more than 30 years. Addressing failures in molten salt TES tanks is fundamental for the CSP industry's survivability, but it is also important for other industrial and power generation applications using this technology, including nuclear and concentrating solar thermal.

14 SOLAR ENERGY↗
DOE OSTI2024

Calibration and validation of the foundation for a multiphase strength model for tin

In this work, the Common Model of Multi-phase Strength and Equation of State (CMMP) model was applied to tin. Specifically, calibrations of the strength-specific elements of the CMMP foundation were developed with a combination of experiments and theory, and then the model was validated experimentally. The first element of the foundation is a 10 multi-phase analytic treatment of the melt temperature and the shear modulus for the solid phases. These models were parameterized for each phase based on ab initio calculations using the software VASP (Vienna Ab initio Simulations Package) based on density functional theory (DFT). The shear modulus model for the ambient β phase was validated with ultrasonic sound speed measurements as a function of pressure and temperature. The second element of the foundation is a viscoplastic strength model for the β phase, upon which strength for inaccessible higher-pressure phases can be scaled as necessary. The stress-strain response of tin was measured at strain rates of 10 -3 to 3 x 10 3 s -1 and temperatures ranging from 87 to 373 K. The Preston-Tonks-Wallace (PTW) strength model was fit to that data using Bayesian model calibration. For validation, six forward and two reverse Taylor impact experiments were performed at different velocities to measure large plastic deformation of tin at strain rates up to 10 5 s -1 . The PTW model accurately predicted the deformed shapes of the cylinders, with modest discrepancies attributed to the inability 20 of PTW to capture the effects of twinning and dynamic recrystallization. Some material in the simulations of higher velocity Taylor cylinders reached the melting temperature, thus testing the multiphase model because of the presence of a second phase, the liquid. In simulations using a traditional modeling approach, the abrupt reduction of strength upon melt resulted in poor predictions of the deformed shape and non-physical temperatures. With CMMP, the most deformed material points evolved gradually to a mixed solid-liquid but never fully liquid state, never fully lost strength, 25 stayed at the melt temperature as the latent heat of fusion was absorbed, and predicted the deformed shape well.

36 MATERIALS SCIENCE↗
DOE OSTI2024

First-Principles Studies of Tritium Species Diffusivity Across the Interface of Nickel-Plated Zircaloy-4

In this study, we will calculate the 3H2 and 3H2O diffusion mechanism in the Ni coated layer and the dominant 3H species dissolution into the Zircaloy-4 getter. We will optimize the Ni/Zr interface and introduce Sn and O impurities. To achieve the objectives, we will perform the investigation by using spin-polarized density functional theory (DFT) methods to capture the effects of Ni magnetic properties on 3H diffusion. We will introduce NiO/Ni(OH) on the Ni side closer to the interface. Then, the 3H diffusion pathways and barriers from Ni side to the zircalory-4 will be explored. We will investigate how NiO/Ni(OH) affects the 3H diffusion across the interface.

Paudel, Hari P.↗
DOE OSTI2024

DECOVALEX-2023: Task F1 Final Report

DECOVALEX-2023 Task F is a comparison of models and methods for post-closure performance assessment (PA) of a deep geologic repository for radioactive waste. The general aims of Task F are to build confidence in the models, methods, and software used for PA and to stimulate additional research and development in PA methodologies. The task objectives are to motivate development of PA modelling skills and capabilities, to examine the influence of model choices on calculated repository performance, and to compare the uncertainties introduced by model choices to other sources of uncertainty. Task F involves no actual experiment or site. It is a PA modelling exercise that requires the conceptual development of hypothetical repository designs and geologic settings. Because three of the teams were interested in salt and the rest of the teams were interested in crystalline rock, Task F was split into two branches: Task F1 for crystalline rock and Task F2 for salt. This report is for Task F1, crystalline rock. Teams from seven countries (Canada, Czech Republic, Germany, Korea, Sweden, Taiwan, and United States) participated in Task F1. The teams worked together to define the features, events, and processes of the reference case repository and established a set of performance measures. In addition, they defined a set of benchmark problems designed to test and compare modelling capabilities for fracture flow and transport at different scales. The repository design and benchmark problems are documented in a Task Specification that evolved over time as the group honed the specifications. The benchmark problems verified that each team can aptly model flow and transport in fractured media in 1-, 2-, and 3-dimensions. Two general approaches were used for the 3-dimensional benchmarks: discrete fracture network (DFN) and equivalent continuous porous medium (ECPM). DFN modelling involves explicit meshing of each fracture while ECPM modelling aims to capture the effective porosity and directional permeability of each cell in a space-filling mesh as affected by intersecting fractures. In some models, a combination of the two is used, i.e., DFN for large known fractures and ECPM for the rest of the domain. Transport is solved by using either the advection-dispersion equation or particle tracking. Although some variation is observed among model breakthrough curves in the benchmark problems, there is strong agreement in breakthrough behaviour up to at least the 75 th percentile for all benchmarks. At the 90 th percentile, breakthrough results show larger differences, suggesting several models retain substantially higher fractions of tracer in regions of slower moving water. In addition to the flow and transport benchmarks, several teams completed the source term benchmark, verifying capabilities for modelling radionuclide decay and ingrowth, waste package breach, instant release fractions, fuel matrix degradation rates, and radionuclide solubility limitations. The reference case is conceptualized as a generic spent fuel repository at a depth of 450 m in fractured crystalline rock. The repository has 50 parallel backfilled drifts, each with 50 deposition holes 6 m apart. Each deposition hole contains a 4-PWR waste package and bentonite buffer. The rock domain is 5 km in length, 2 km in width, and 1 km in depth. It has 6 deterministic fractured deformation zones and a multitude of stochastic fractures. Teams generally used the ECPM approach for the entire rock or a hybrid approach in which the deterministic fracture zones are modelled with a DFN and the rest of the rock is modelled by ECPM. Of the reference case problems specified, only the results of the initial reference case problem are compared in this report. The initial problem focuses on transport from the deposition holes to the surface, i.e., it neglects waste package performance. Tracers are released at all waste package locations at time zero and tracked for their releases to the near field and ground surface. The water fluxes calculated at the ground surface entry and exit regions of the domain are similar for all models except for two that have considerably lower fluxes. For tracer transport, large differences are observed among models in the magnitude of tracer transported. Much of the difference appears to be due to how the repository is implemented and hence the different degrees of repository simplification. Models that exclude the drifts, buffer, and backfill from the domain tend to show greater release of tracers and radionuclides from the repository. The initial study presented here indicates that major differences in modelling important processes within the repository (e.g., diffusion through buffer and backfill) can produce broadly different release and transport results, especially when those processes are excluded. Even for the models that included all specified features, events, and processes, the results show significant differences and demonstrate the importance of examining multiple modelling approaches in performance assessment. The differences in results observed in this study are expected to motivate teams to either increase complexity in future versions of the reference case models or to improve methods to account for the effects of simplified features and processes. Either way, future improvements in these models are expected to produce results that more closely agree.

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

Challenges in quantifying unparameterized spatial uncertainties in deep geologic repositories for nuclear waste

Spatially heterogeneous uncertainties are prevalent in geophysical modeling applications, such as probabilistic post-closure performance assessment (PA) of deep geologic repositories for nuclear waste. Such uncertainties are often highly influential to model outputs, so it is desirable to identify the most important mechanisms by which they influence model predictions. However, these uncertainties are often unparameterized in the sense that there is no set of parameters that can be specified to yield a particular realization of the uncertainty. Additionally, the uncertainty is not intrinsically endowed with a parameterization that captures a realization’s mechanistic influence on model outputs. Therefore, in this work we present a novel methodology to develop and assess a set of proxy variables that aim to represent this influence. We show how they can be computed, downselected, and incorporated into the construction of statistical surrogate models mapping model inputs to outputs. We present our methodology in the context of a motivating application problem in deep geologic repository PA and discuss the challenges in capturing the effects of these spatial heterogeneities in uncertainty analyses.

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

Optimize the Autonomous Sensor Fish Device for Understanding Interactions of Aquatic Animals with Marine and Hydrokinetic and Hydro Systems (CRADA 476) Final Report

This project at the Pacific Northwest National Laboratory (PNNL) enhances our understanding of how aquatic life interacts with marine renewable energy (MRE) and hydropower systems through the development and testing of advanced sensor technology. The Marine Sensor Fish (MSF) and its smaller variants, the Sensor Fish Mini (SF Mini) and Flexible Sensor Fish (FSF), measure critical environmental stressors, such as shear forces, pressure changes, and strike impacts, that are essential for assessing the risks posed by energy devices, including tidal turbines, hydropower installations, and other MRE infrastructure. Rigorous field testing has shown that these tools are effective in capturing detailed data on turbulence, pressure variations, strike impacts, and other stressors near these devices, providing valuable insights into the conditions faced by aquatic species. Economically, these tools offer renewable energy developers a practical solution for streamlining environmental assessments and ensuring regulatory compliance, ultimately reducing costs associated with evaluating and mitigating ecological impacts. This research also benefits the public by promoting the growth of renewable energy sources that protect aquatic ecosystems, advancing both sustainable energy production and environmental stewardship.

13 HYDRO ENERGY↗
DOE OSTI2026

Model Parameter Development for Complex Materials: Species-Specific Diffusion Barriers in 316 Stainless Steel from Systematic DFT Calculations

Vacancy-mediated diffusion barriers in 316 stainless steel have been systematically calculated using density functional theory to provide essential parameters for mesoscale microstructure evolution models. A statistical sampling approach employing 210 nudged elastic band calculations across multiple special quasi-random structures captures the effects of local chemical environments in this concentrated alloy. The computational methodology addresses challenges specific to chemically disordered systems, including proper magnetic treatment throughout multi-step calculations and validation against experimental structural properties. The calculated activation barriers reveal clear species-dependent diffusion behavior with the hierarchy Ni >> Fe ˜ Cr >> Mo. Nickel exhibits the highest barriers (0.74–1.31 eV, mean 1.045 eV), confirming its role as the slowest-diffusing major component. Iron and chromium show similar moderate barriers averaging 0.587 eV and 0.522 eV, respectively. Remarkably, molybdenum demonstrates exceptionally low barriers (0.12–0.28 eV, mean 0.194 eV), suggesting much higher mobility than previously recognized and potentially significant implications for precipitation kinetics and microstructure evolution. The barrier ranges remain consistent across different 316 SS compositions, supporting parameter transferability for modeling applications. The overall mean barrier of 0.64 eV provides a practical approximation for phase field simulations, while species-specific values enable detailed treatments of diffusion-controlled processes. This systematic approach establishes a validated framework for generating diffusion parameters in other concentrated alloys where experimental data are limited, while providing the first systematic set of species-specific barriers for predictive modeling of 316 stainless steel microstructure evolution.

36 MATERIALS SCIENCE↗