Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Design verification”

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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

An approach for spent nuclear fuel containment integrity verification using gas tagging

Verification of containment integrity is required for spent nuclear fuel (SNF) managed by the commercial nuclear industry and U.S. Department of Energy (DOE), especially after extended storage. Certain SNF storage systems, such as the DOE road-ready dry storage system, hold several packaged containments within a welded over-canister. These packaged containments are called Department of Energy Standard Canisters (DOESCs). DOESC leakage identification is challenging because their containment boundary cannot be accessed for testing and their contents (i.e., SNF and fill gas) are often similar. There are concerns that this could result in costly characterization and repackaging operations of DOE road-ready dry storage systems if compromised DOESCs are suspected. Here, to address these concerns, this paper presents an approach for applying a gas tagging process using xenon to uniquely identify compromised inaccessible containments following extended storage. The containments considered for this application are seven DOESCs, each packaged within a single over-canister. Two different SNF loading configurations from the Advanced Test Reactor and Fort Saint Vrain nuclear power plant are considered. These configurations are used to represent research reactor aluminum-clad spent nuclear fuel (ASNF) and TRi-structural ISOtropic (TRISO) SNF types. Results for this application show that for ASNF and TRISO type fuels for which the selected fuels are representative, the volume of taggant required at loading is determined primarily by the lower detection limit and leak rate of taggant from a compromised DOESC, rather than the amount of fission-generated xenon in the loaded fuel. While the application presented is suited for larger leaks, smaller leaks could be detected by modifying certain design parameters. This gas tagging approach can also be applied to other DOE containments and advanced reactor SNF storage systems.

07 - ISOTOPES AND RADIATION SOURCES↗
DOE OSTI2026

The future of subsurface monitoring: AEC’s breakthroughs in CCS technology

Carbon capture and storage (CCS) has emerged as a key solution in the fight against climate change. However, for CCS to succeed, it is crucial to ensure that the sequestered CO2 stays safely trapped underground. The U.S. Department of Energy (DOE) has emphasized the need for advancements in subsurface monitoring, measurement, reporting, and verification. Aside from caprock integrity failure, the other primary failure points usually involve defective cement in the casing annulus of wellbores or plugged and abandoned wells. In addition, many energy producers (e.g., oil and gas, geothermal) and storage and disposal operators (e.g., H2 and water) must deal with the same issue. Poorly placed or degraded cement can create pathways for gas or fluid to escape from casing annuli and in plugged and abandoned or orphan wells, posing environmental risks. Yet, a reliable and cost-effective way to monitor cement and well integrity over multiple decades is still unavailable. Traditional geophysical methods like 4D seismic imaging and surface-based electromagnetic monitoring lack the resolution and accuracy for detecting these types of failures (Vasco et al., 2022; Fawad and Mondol, 2021). Wireline logging is expensive to run continuously and is obtrusive to the operation. While fiber optics can potentially be a solution, its bulkiness can significantly compromise the cement's integrity. To address these challenges, the Advanced Energy Consortium (AEC) at The University of Texas at Austin’s Bureau of Economic Geology (the Bureau) has been pioneering research in subsurface monitoring using its portfolio of distributed autonomous microfabricated sensors for harsh subsurface environments since 2008. A class of these microsensors [System on a Chip (SoC)] can be mixed in cement and permanently placed without compromising the cement column; the sensors would then communicate with each other or a data acquisition (DAQ) master node. Another class of the AEC microsensors can be fully autonomous, with rechargeable micro-batteries capable of exceeding 100°C, flash memory, and, currently, a pressure and temperature sensor. They are designed to circulate in mud, geothermal fluids, U-loops, or pipelines. They can log data into memory and are unobtrusive to operations. Our team has been working on a multi-year DOE-funded project (DE-FE0031856)—supported by $2.95M in federal funding and $0.75M in cost-matching from the AEC—to demonstrate SoC sensor utility for CO2 leakage monitoring in CCS applications. This multi-institutional collaboration developed a novel sensing architecture utilizing radiofrequency (RF) microsensors embedded within the cement sheath. These sensors detect CO2 migration and are interrogated via a Smart Casing Collar (SCC).

58 GEOSCIENCES↗
DOE OSTI2020

Initial Atmospheric Transport of Particles

A new version of the Initial Atmospheric Transport (IAT) model has been developed. This report fully characterizes the new model. IAT is designed to model heat release into the atmosphere that contains contaminated particles. The heat release forms a buoyant puff or plume. Buoyancy is quickly dissipated in 10's to 100's of seconds through drag and entrainment of ambient air. The final location of contaminant particles in the atmosphere after dissipation of abnormal heat is an important input to long term transportation codes such as the National Oceanic and Atmospheric Administration's (NOAA) Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) code. Ultimately the settling of contaminants after long term transport provides a basis for estimating health consequences. In comparison to the old version, the new model includes explicit simulation of particles, newly formulated equations of motion for a rising puff, and a new plume model. The plume model releases many instances of a single puff along a single trajectory in a way that balances energy released over a longer period than would be appropriate for a single puff. The puff model has additions for vorticity, turbulence, virtual mass, drag, and entrainment. The particle model leads to direct interfacing between IAT and the long-term transport code HYSPLIT. Previously, particle locations were given a pattern based on IAT's puff trajectory with no particle simulation. Now, particles are explicitly simulated, and 3-D spatial distributions of particles are output based on final positions. Ballistic particles are thrown outward at high velocity from explosive events and tend to hit the ground before the end of the simulation. Non-ballistic particles tend to get pulled into the rising puff s ring of vorticity with some escaping to ambient conditions as buoyant energy dissipates. A new model verification and validation effort was made for the new plume model. The Wallops Island Solid Rocket Propellant (SRP) fire plume was extracted to a 3-D transient dataset. After calibration, the model performs better than previous single puff attempts. Even so, error levels in the process used to estimate location and ambient weather conditions approach 30%. This leads to the conclusion that unplanned photogrammetry evaluations of the IAT model are useful but not sufficiently accurate to quantify model parameters. Like previous attempts, the vertical motion entrainment parameters come out much too high in comparison to values for other studies in the literature. It is suspected that a new model of shear entrainment would alleviate this problem. Unfortunately, IAT already has more calibration parameters than the data comparisons thus far can resolve with statistical certainty. This underscores the need to verify and validate the IAT model using detailed CFD studies that allow ideal statistical comparison with the IAT model for hundreds of cases. Even so, validation by comparison to carefully planned experiments (i.e. with acceptable error levels) is also necessary because the dynamics of entrainment are at a sub-grid level for CFD. Though better validation and verification are most important, many new enhancements are still envisioned for future IAT development. These include particle rain out due to condensation and agglomeration of water, direct connection to the North American Mesoscale (NAM) weather dataset, dynamic weather conditions, including ambient particle entrainment into the puff, and detrainment at the puff to ambient mixing interface. In general, including detailed aerosol dynamics is the greatest coming challenge which will require an entire new sub-model in IAT. This is very important if contaminant particles change in such a way that they increase or decrease in their potential harm to the environment.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Quantification of Residual Water in Spent Fuel Dry Storage Canisters Using Mass Spectrometry

The purpose of this report is to document updates on the apparatus to simulate commercial vacuum drying procedures at the Nuclear Energy Work Complex at Sandia National Laboratories. Validation of the extent of water removal in a dry spent nuclear fuel storage system based on drying procedures used at nuclear power plants is needed to close existing technical gaps. Operational conditions leading to incomplete drying may have potential impacts on the fuel, cladding, and other components in the system during subsequent storage and disposal, such as fuel degradation; cladding corrosion, embrittlement, or breaching; and the creation of a flammable environment via radiolysis of water. A general lack of data suitable for model validation of commercial nuclear canister drying processes necessitates well-designed investigations of drying process efficacy and water retention. Scaled tests that incorporate relevant physics and well-controlled boundary conditions are essential to provide insight and guidance to the simulation of prototypic systems undergoing drying processes. This report documents details on the quantification of residual water in the Advanced Drying Cycle Simulator (ADCS), an apparatus built to simulate commercial drying procedures and quantify the amount of residual water remaining in a pressurized water reactor (PWR) fuel assembly after drying. The ADCS was constructed with a prototypic 17×17 PWR fuel skeleton and waterproof heater rods to simulate decay heat. The ADCS is outfitted with thermocouples to measure the thermal response of the ADCS to simulated decay heats and internal helium fill pressures relevant to commercial drying procedures. The ADCS is also instrumented with pressure transducers to measure the pressures and vacuum levels observed during simulated commercial drying. The most unique instrumentation used for quantifying residual water in the ADCS is a Hiden Analytical HPR-30 mass spectrometer (MS), which measures gas compositions of the ADCS internal free volume, based on partial pressures calculated from relative proportions of gas molecules detected by the MS. This report details the methodology used to implement MS measurements in quantifying residual water in the ADCS. This methodology includes the calibration of the HPR-30 MS to a Buck Research Instruments CR-4 chilled mirror hygrometer, which itself is calibrated to a NIST-traceable standard. Data collected by both the MS and the chilled mirror hygrometer from water/helium mixtures ranging from 150 to 500,000 ppmv water in helium were used to generate calibration curves, establishing a source of verification of MS measured water contents. Details regarding water content measurement uncertainties are included in this report, defining the accuracy and verifiability of the HPR-30 MS in measuring residual water content in simulated dry storage canister environments.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Validating Intra-Plate Power Distributions in the Advanced Test Reactor

We report that the neutronics software, HELIOS, was validated in 2015 for performing core reload design and safety analysis of the Advanced Test Reactor. However, when HELIOS was benchmarked against historic fission-wire measurements (i.e., zero-power full-core measurements) a statistically resolved calculation-to-measurement bias was discovered. The azimuthal power along each fuel plate computed by HELIOS has consistently shown to under-predict measurements made by fission-wires in historic zero-power tests near the fuel element side-plates. It was hypothesized during the HELIOS software validation work, that this bias is attributable to local moderation in coolant vents in the side-plates axially just above and below the fission-wires on the fuel-plate edges. This work used detailed MCNP and MC21 models of the side-plate vents to test this hypothesis. By comparing the average azimuthal biases between HELIOS, and 2D and 3D MCNP models, and a 3D MC21 model, it was found that the HELIOS azimuthal bias is not due to the measurement.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

Summary of Graphite Data Stored within NDMAS

The Graphite Technology Development Project provides data to support the design of graphite core components within specific reactor service conditions of the next generation of high-temperature, gas-cooled nuclear reactors. Physical, mechanical, and thermal properties of nuclear grade graphite were characterized for specimens that were unirradiated, irradiated, and irradiated under various stress conditions. The material properties include diameter, length, mass, density, compressive strength, tensile strength, flexural strength, modulus, resistivity, thermal diffusivity, and thermal expansion coefficient. Baseline graphite specimens are unirradiated from different grades (2114, IG-110, NBG-17, NBG-18, and PCEA) and different types used in different characterization tests (compressive, flexural, tensile, one-inch cylinder, and quarter-inch cylinder). The Advanced Graphite Creep (AGC) irradiation specimens are from a much larger number of grades and cylinder types (creep, piggyback, and pencil). For baseline graphite, characterization data for 7,756 specimens extracted from thirteen graphite billets were captured to the NDMAS database. For AGC experiments, four irradiation campaigns have been completed: AGC-1, AGC-2, AGC-3, and AGC-4. The ongoing HDG-1 (High-Dose Graphite) experiment, which began irradiation with Cycle 168B on August 26, 2020, includes specimens previously irradiated in AGC-2 in addition to the specimens originally destined for AGC-5.. Besides the characterization data, the AGC data includes irradiation monitoring and physics data representing the irradiation conditions of AGC specimens. Currently, all data for AGC-1, AGC-2, and AGC-3 have been captured to the NDMAS database. Only AGC-4 pre-irradiation and irradiation monitoring data have been captured, and HDG-1 pre-irradiation data are in process of being captured for unirradiated specimens. To date, a total of 38,149 material property records have been captured into NDMAS database for baseline and AGC specimens. All characterization data are qualified for use according to their perspective data verification reports. For the AGC irradiation campaigns, a total of 41,502 qualified physics calculation records were added for AGC-1, AGC-2, and AGC-3. Finally, a total of 173,698,505 AGC irradiation monitoring data records (thermocouple temperature, gas flow rate, gas pressure, gas moisture, applied load, and specimen displacement) have been captured to NDMAS; the majority of those records (~94%) are qualified data records.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Simulation Center for Runaway Electron Avoidance and Mitigation (SCREAM SciDAC) (Technical Final Report)

Runaway electrons can severely damage the plasma facing components on ITER during a major disruption and pose a major risk for tokamak fusion. It has been recognized that an adequate disruption mitigation system (DMS) is essential for the safe operation of ITER. The United States is responsible for the design and implementation of the disruption mitigation system on ITER, and in July 2016 the Simulation Center for Runaway Electron Avoidance and Mitigation (SCREAM) was launched by DOE, in a joint Fusion Energy Sciences (FES) and Advanced Scientific Computing Research (ASCR) collaboration. SCREAM was a comprehensive theory and simulation SciDAC center that provided physics guidance in the avoidance and mitigation of runaway electrons, and in tandem with domestic and international experiments, helped establish the qualitative and quantitative bases for safe operational scenarios and viable mitigation techniques. The SCREAM center assembled a national team of experts in runaway electron physics, tokamak disruptions, magnetohydrodynamic (MHD) simulation, and advanced algorithms and computing. The team combined advanced simulation and analysis capability facilitated by direct participation of ASCR SciDAC institutes with theoretical models and code development by FES scientists to focus on the runaway risk for ITER and tokamaks in general. The research scope was focussed on integrated simulations of kinetic runaway electrons, including MHD and fluid models of impurity transport, within a research plan guided by theory. The specific research tasks were (1) establish the fundamental physics of runaway generation, saturation, and dynamical evolution in a tokamak; (2) examine the critical path toward runaway avoidance; and (3) investigate the viability and effectiveness of the leading candidate schemes for runaway mitigation. In all three areas, members of the team carried out scoping studies that established the readiness for rapid and critical advances, especially in the deployment and further development of large-to extreme-scale simulation tools. Our multi-pronged computational approach included (1) relativistic Fokker-Planck solvers with discretization in phase space, (2) self-consistent particle-in-cell techniques, (3) particle-based Monte-Carlo, and (4) MHD-particle hybrid simulations. Cross-check between these different methods provided an additional means for verification and further bolstered the fidelity of our physics prediction. Validation against experimental results brings confidence to the predictive capability for ITER and frequently leads to new ideas for understanding and mitigating the thermal quench driven runaway electron phenomenon.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Development and Validation of Two-Phase Flow Models in MOOSE and Application to Molten Salt Reactors

Two-phase flow in Molten Salt Reactors (MSRs) is important as it impacts reactivity evolution, reactor transient response, and the removal of species dissolved in the molten salt through gas phase transfer. Therefore, accurately predicting the gas distribution and the associated liquid-gas interface area in MSRs is essential for their design and operation. Recently, we integrated two new models into Idaho National Laboratory (INL)’s Multiphysics Object-Oriented Simulation Environment (MOOSE): a multi-D generalization of a mixture drift-flux model and a Euler-Euler model. The Euler-Euler model offers higher fidelity, while the mixture drift-flux model provides greater computational efficiency, which is typically preferred for modeling reactor transients. However, the mixture model's accuracy in capturing void distribution and interfacial area in MSRs still needs to be assessed. This article begins with a description of the mathematical framework for the two-phase models implemented in MOOSE. It then presents validation of these models against relevant experimental data. Finally, both models are applied to the Molten Salt Reactor Experiment case study, analyzing various operational conditions such as different rates of fission product volatilization and diverse cover gas entrainment scenarios at the reactor pump. The article concludes by assessing the suitability of both models for capturing the two-phase flow dynamics critical to MSR operations.

42 - ENGINEERING↗
DOE OSTI2023

Project: Corbomite - Product: ConsoleWorks REACT

TDi Technologies presents ConsoleWorks REACT, an advanced platform designed to tackle the complexities of cyber and operational risk assessment. This comprehensive solution goes beyond asset-focused approaches by considering the impact of both assets and people have on the security and operation of critical infrastructure, specifically targeting preventing gird mis-operation by evaluating real-time human interaction or commands with critical assets. The hypothesis suggests that by integrating the assessment of user commands into the overall risk assessment process, organizations can make more informed decisions, prioritize resources effectively, and respond promptly to potential threats. This hypothesis forms the basis for the development of a proactive and holistic risk management approach that is more comprehensive and context aware than traditional models which only look at assets, patch levels, configuration and threat intel by collecting that information off the network vs directly from the asset, human, and human interaction all in real-time. ConsoleWorks' unique man-in-the-middle architecture is a key feature that sets it apart in the cybersecurity landscape. This architecture enables the real-time observation, enforcement, and commands or interaction risk transparency of user interaction with critical infrastructure to be risk mitigated and audited as they are aggregated with device and human risk factors for a more comprehensive risk threat score across a device or group of devices. This architecture allows ConsoleWorks to act as a secure intermediary between users and critical assets, monitoring all interactions and ensuring that only authorized commands are sent to the asset to be executed. This not only enhances security but also provides a comprehensive audit trail of all user activities, contributing to compliance efforts and facilitating incident investigation and risk management. By integrating this unique architecture with our comprehensive risk assessment methodology, ConsoleWorks REACT provides a powerful solution for managing cyber and operational risks, enabling organizations to maintain a robust security posture and effectively mitigate potential threats. To that end, the primary objective of this project was to research and develop a robust solution that enables the energy industry to mitigate the risks associated with human actions that can compromise the security or operations of assets critical to energy delivery, generation, transmission and operation. ConsoleWorks REACT plays a pivotal role in achieving this goal by leveraging its unique capabilities to monitor and track all user activity. Through its Zero Trust approach, which emphasizes continuous verification, the platform ensures secure access to assets and serves as the centralized human response and notification platform for addressing cyber and operational issues.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Verification of diesel spray ignition phenomenon in dual-fuel diesel-piloted premixed natural gas engine

Dual-fuel (DF) engines, in which premixed natural gas and air in an open-type combustion chamber is ignited by diesel-fuel pilot sprays, have been more popular for marine use than pre-chamber spark ignition (PCSI) engines because of their superior durability. However, control of ignition and combustion in DF engines is more difficult than in PCSI engines. In this context, here we focus on the ignition stability of n-heptane pilot-fuel jets injected into a compressed premixed charge of natural gas and air at low-load conditions. To aid understanding of the experimental data, chemical-kinetics simulations were carried out in a simplified engine-environment that provided insight into the chemical effects of methane (CH4) on pilot-fuel ignition. The simulations reveal that CH4 has an effect on both stages of n-heptane autoignition: the small, first-stage, cool-flame-type, low-temperature ignition (LTI) and the larger, second-stage, high-temperature ignition (HTI). As the ratio of pilot-fuel to CH4 entrained into the spray decreases, the initial oxidization of CH4 consumes the OH radicals produced by pilot-fuel decomposition during LTI, thereby inhibiting its progression to HTI. Using imaging diagnostics, the spatial and temporal progression of LTI and HTI in DF combustion are measured in a heavy-duty optical engine, and the imaging data are analyzed to understand the cause of severe fluctuations in ignition timing and combustion completeness at low-load conditions. Images of cool-flame and hydroxyl radical (OH*) chemiluminescence serve as indicators of LTI and HTI, respectively. The cycle-to-cycle and spatial variation in ignition extracted from the imaging data are used as key metrics of comparison. The imaging data indicate that the local concentration of the pilot-fuel and the richness of the surrounding natural-gas air mixture are important for LTI and HTI, but in different ways. In particular, higher injection pressures and shorter injection durations increase the mixing rate, leading to lower concentrations of pilot-fuel more quickly, which can inhibit HTI even as LTI remains relatively robust. Decreasing the injection pressure from 80 MPa to 40 MPa and increasing the injection duration from 500 µs to 760 µs maintained constant pilot-fuel mass, while promoting robust transition from LTI to HTI by effectively slowing the mixing rate. This allows enough residence time for the OH radicals, produced by the two-stage ignition chemistry of the pilot-fuel, to accelerate the transition from LTI to HTI before being consumed by CH4 oxidation. Thus from a practical perspective, for a premixed natural gas fuel–air equivalence-ratio, it is possible to improve the “stability” of the combustion process by solely manipulating the pilot-fuel injection parameters while maintaining constant mass of injected pilot-fuel. This allows for tailoring mixing trajectories to offset changes in fuel ignition chemistry, so as to promote a robust transition from LTI to HTI by changing the balance between the local concentration of the pilot-fuel and richness of the premixed natural gas and air. This could prove to be a valuable tool for combustion design to improve fuel efficiency or reduce noise or perhaps even reduce heat-transfer losses by locating early combustion away from in-cylinder walls.

42 ENGINEERING↗
DOE OSTI2022

[Presentation Slides] Code-to-Code Benchmark Study for Thermal Stress Modeling and Preliminary Analysis of the High-temperature Single Heat-Pipe Experiment

In the heat-pipe-cooled microreactor, heat pipes remove heat from the reactor core as a passive heat-transfer device, so the fluid circulation is not required for cooling, which can substantially simplify the overall reactor design. However, given the extremely high temperatures in the core region and potentially large temperature gradients across the structure materials, thermal stresses need to be well-analyzed to ensure structural integrity during normal operations and accident scenarios. This presentation slides discuss finite element method-based thermal-stress analysis for the high-temperature single heat-pipe test article in the Single Primary Heat Extraction and Removal Emulator (SPHERE) facility at Idaho National Laboratory (INL), using two commercial software packages, Abaqus and Star-CCM+. A code-to-code benchmark study was performed to crosscheck the model setup and capability of each code and to gain preliminary insights into the potential thermal stress concerns from the current experimental setup. It is observed the significant thermal stresses happen at the inner surface of the heat pipe hole surrounded by electric heaters where the largest temperature gradients appear. The temperature fields have good agreements between Abaqus and Star-CCM+, while the induced thermal stresses show modest deviations probably due to differences of meshing engines used in these two codes. It is found that the local maximum thermal stresses may reach close to the ultimate tensile strength and yield strength of structural material depending on the heater power. Ultimately, the coupled thermal-structural analysis will help guide the current experimental plan and ensure the facility safety for future experimental study.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

HERO CarbonSAFE Phase 2 Project in the Columbia River Basalt Group

The Hermiston, Oregon Basalt CarbonSAFE Phase II project (HERO CarbonSAFE) seeks to accelerate the deployment of commercial carbon dioxide (CO2) storage projects in basaltic rocks. Basalt CO2 storage has several advantages to conventional saline storage reservoirs including 1. The potential for rapid mineralization of CO2, 2. Associated decreases in pressure and CO2 migration risks, 3. Reduced long-term monitoring requirements with respect to plume tracking, 4. Widespread geographic distribution and, 5. Large storage potential due to thickness, porosity, and CO2 interactions with basalt. And for locations such as the Pacific Northwest, Hawaii, Iceland, India and Japan, basalts may offer the only economically feasible option for local CO2 storage. However, there are limited field-scale assessments of CO2 storage in basalt, and current carbon capture utilization and storage (CCUS) permitting and regulatory frameworks were developed for conventional saline reservoirs. HERO CarbonSAFE is designed to address research gaps and uncertainties associated with basalt storage. Specifically, the project will assess the feasibility of CO2 injection in the deep layered basalts, long-term storage (mineralization), practical approaches for large-scale implementation (50+ million metric tons of CO2 over 30 years), lithology-specific risks, and the technoeconomic potential for CO2 storage in basalts. The HERO CarbonSAFE project will assess feasibility of developing a commercial-scale (50+ million metric tons of CO2) geological storage complex within the Columbia River Basalt Group (CRBG), a layered continental flood basalt complex that underlies Calpine’s natural gas-fired Hermiston Power Project (HPP) in Hermiston, OR (Figure 1). Under this 2-year CarbonSAFE Phase II project, the HERO team will conduct a data acquisition campaign that includes drilling a stratigraphic well to a total depth of ~1,500 m into the thick layered basalts proximal to HPP. A comprehensive well logging and hydrologic testing program will be augmented with new core collected from flow zones and sealing units, and comprehensive laboratory testing to help refine the kinetic rates of mineralization. The newly acquired information will be integrated with existing data from regional wells to correlate basalt injection zone properties to develop storage hub/commercial-scale models. Using these models, the project team will evaluate injection scenarios to define the technical and economic potential for storing a minimum of 50 million metric tons of CO2 over a 30-year period, along with a robust sensitivity analysis on key parameters governing reservoir viability for sustainable injection over a commercial project lifetime. Specific technical objectives of HERO are: (1) assessing the reservoir response of a series of stacked layered reservoir flowtop sequences occurring in this area of the CRBG to commercial-scale injection volumes; (2) extending prior efforts by the project team to characterize the deep layered basalts encountered in regional studies, to leverage prior investments by U.S. Department of Energy’s (DOE) Carbon Storage program; (3) leveraging DOE’s mineralization characterization efforts to advance model parametrization for commercial scale injection of CO2 in basalts; (4) conducting risk assessments associated with scaling up to commercial storage hub injection goals, while validating DOE’s National Risk Assessment Partnership (NRAP) tools, to identify potential constraints that would prevent the CRBG from serving as a commercial-scale storage complex; (5) developing mitigation plans to address identified risks; (6) developing a commercial-scale injection and monitoring, verification and accounting (MVA) strategy; (7) utilizing computational models to define and minimize, if possible, the Area of Review (AoR) under Class VI regulations; and (8) developing a robust CO2 management strategy for CRBG that also considers a regional source/sink approach that is responsive to stakeholder needs and industrial demand. Specific institutional objectives are: (1) identifying and developing plans to mitigate the nontechnical challenges associated with the build-out of a commercial-scale storage complex within the CRBG with integrated CO2 sources; (2) implementing the community outreach plan; (3) conducting regulatory research, including a survey of issues related to pore space ownership, MVA and long-term assurance of mineralization-based storage, to support an eventual application for a UIC Class VI permit; (4) advancing the project’s plan for CO2 liability management; and (5) continuing to refine and update the project’s economic model. The final objective is the preparation of a comprehensive Site Characterization Plan that draws upon the technical and institutional feasibility assessments to prepare the project for future commercialization efforts.

58 GEOSCIENCES↗
DOE OSTI2025

Analysis and optimization of seismic monitoring networks with Bayesian optimal experimental design

SUMMARY Monitoring networks increasingly aim to assimilate data from a large number of diverse sensors covering many sensing modalities. Bayesian optimal experimental design (OED) seeks to identify data, sensor configurations or experiments which can optimally reduce uncertainty and hence increase the performance of a monitoring network. Information theory guides OED by formulating the choice of experiment or sensor placement as an optimization problem that maximizes the expected information gain (EIG) about quantities of interest given prior knowledge and models of expected observation data. Therefore, within the context of seismo-acoustic monitoring, we can use Bayesian OED to configure sensor networks by choosing sensor locations, types and fidelity in order to improve our ability to identify and locate seismic sources. In this work, we develop the framework necessary to use Bayesian OED to optimize a sensor network’s ability to locate seismic events from arrival time data of detected seismic phases at the regional-scale. This framework requires five elements: (i) A likelihood function that describes the distribution of detection and traveltime data from the sensor network, (ii) A prior distribution that describes a priori belief about seismic events, (iii) A Bayesian solver that uses a prior and likelihood to identify the posterior distribution of seismic events given the data, (iv) An algorithm to compute EIG about seismic events over a data set of hypothetical prior events, (v) An optimizer that finds a sensor network which maximizes EIG. Once we have developed this framework, we explore many relevant questions to monitoring such as: how to trade off sensor fidelity and earth model uncertainty; how sensor types, number and locations influence uncertainty; and how prior models and constraints influence sensor placement.

58 GEOSCIENCES↗
DOE OSTI2024

Bayesian Analysis of TRISO Fuel: Quantifying Model Inadequacy, Incorporating Lower-Length-Scale Effects, and Developing Parallel Active Learning Capabilities

The U.S. Department of Energy (DOE)’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program aims to develop predictive capabilities by applying computational methods to the analysis and design of advanced reactor and fuel-cycle systems. This program has been providing engineering-scale support for the continued development of BISON, a high-fidelity, high-resolution fuel performance tool. Fuel behavior in nuclear reactors is governed by a complex network of mechanisms that interact with various other physics aspects in the reactor system. Any model developed to represent fuel behavior will likely be idealized, resulting in uncertainties when comparing their predictions against the observed data. In Fiscal Year (FY)-23, we initiated the Uncertainty Quantification (UQ) work by using Bayesian methods to establish a level of model trustworthiness and further improve it, with a particular emphasis on TRI-Structural isOtropic (TRISO) nuclear fuel. This year, we further expanded on that UQ work by investigating an approach to quantifying model inadequacy and accounting for lower-length scale (LLS) effects in TRISO silver (Ag) release modeling. Furthermore, we are implementing parallel active learning capabilities to reduce the computational cost (i.e., required computational resources and elapsed time) of performing UQ. Specifically, we utilized The Kennedy O’Hagan framework for Bayesian uncertainty quantification (KOH) to account for model inadequacy in TRISO Ag release predictions made by BISON. The KOH framework represents an improvement over the standard Bayesian framework used in FY-23. Explicitly accounting for model inadequacy in the Bayesian framework helps establish the level of experimental noise uncertainty in the Advanced Gas Reactor (AGR) data. We compared the inverse UQ results obtained from both the standard Bayesian and KOH frameworks in light of the AGR-2/3/4 data, and also compared the predictive UQ results obtained from these two frameworks in light of the AGR-1 data. Next, we investigated the impact of considering LLS effects in the Ag release simulations. We developed an expanded database of LLS simulated effective diffusivities for Ag, covering a wide range of microstructures and temperatures. Using this database, we developed a framework for incorporating LLS effects into the engineering-scale Ag release UQ. We developed both parametric and non-parametric approaches for bridging the length scales. We then investigated the inverse UQ results in light of the AGR-2/3/4 data and the predictive UQ results in light of the AGR-1 data, and compared the LLS-informed approach and the Arrhenius equation, which does not include microstructure information. Finally, we discussed implementing parallel active learning capabilities in the Multiphysics Object Oriented Simulation Environment (MOOSE)/BISON to reduce the computational cost (i.e., computational resources and elapsed time) of Bayesian UQ. For verification purposes, we first tested these new capabil ities on a species interaction problem. We then demonstrated them on the TRISO Ag release application, showing that parallel active learning capabilities can enhance the accuracy of UQ while also substantially reducing the computational cost in comparison to the reference methods developed in FY-23.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Guidance for Monitoring Passive Groundwater Remedies Over Extended Time Scales

Passive remediation can be appropriate where natural processes and actions such as institutional controls mitigate exposure to contaminated groundwater, achieving remedial action objectives and protectiveness of human health and the environment. Monitored natural attenuation (MNA) is a prevalent passive remediation strategy supported by a regulatory framework and monitoring design guidance. However, long-term passive remedies are usually selected in combination with at least one active remedy, such as source removal, in situ treatment, or pump-and-treat, functioning as a complementary method for achieving remediation objectives and meeting the applicable statutory and regulatory requirements. However, MNA and existing monitoring guidance primarily target situations where the remedial action objectives are met within a few decades. When time scales for passive remediation extend to many decades (50 years or more), a corresponding change in monitoring strategy is needed to adapt to the extended time scale. This document provides guidance for implementing an extended-scale monitoring (ESM) approach appropriate for long-duration passive remediation. Extended-scale is defined in this document with respect to time (i.e., a long duration of remediation) and a large enough physical scale such that downstream receptors will not be impacted within the remediation timeframe. ESM applies to slow-moving groundwater contaminant plumes and emphasizes monitoring primarily for potential exposure pathways. For this approach, the primary monitoring objective is to demonstrate that the plume diminishes before reaching the receptor zone or point of compliance and/or a receptor does not receive concentrations above the compliance limit. While the overall objectives of protecting human health and the environment are the same as for plumes where remediation can occur over a shorter time period, the time scales between decisions are longer and the dynamics of plume evolution are slower. To this end, a scenario-based strategy is described for different plume and source conditions, defining a containment and receptor zones. The containment zone is the area where the risk of exposure to groundwater contamination can be mitigated (e.g., through institutional controls) during the remediation time period. The receptor zone is defined as the area where exposure to groundwater contamination cannot be mitigated and compliance concentration standards must be met. Within the containment zone, slow plume migration may occur, leading to concentrations that exceed compliance standards. However, where distance to the receptor zone is large relative to plume migration and attenuation rate, this approach can be protective of the receptor zone. Selection of a long-duration passive remedy needs to be based on sufficient understanding of contaminant sources, hydrogeology, and contaminant plumes. A strong technical basis, supported by predictive analysis, is recommended to substantiate that contamination is expected to stay within the containment zone during the active remediation and attainment phase of the remedy, and diminish to meet compliance standards within the extended timeframe prior to reaching the receptor zone (e.g., many decades or even centuries). The ESM approach is based on verification of plume behavior and not on detailed plume dynamics. Monitoring is conducted to confirm expected behavior with an emphasis on exposure pathways to verify that plumes remain contained in areas where the protectiveness objectives can be met. ESM should not be adopted if there is significant risk of the plume extending beyond the containment zone. Given the slow movement within the containment zone, less frequent sampling is required relative to approaches used for conventional-scale remediation.

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

Root genetics in the field to understand drought adaptation and carbon sequestration (Final Scientific/Technical Report)

For all crop plants, roots play a critical role in growth. Roots anchor the plants, and are the primary site of nutrient and water uptake. Roots are also the main source of C to soil in the form of root tissues and exudates, and thus greatly influence SOM stocks. To perform these functions, primary roots extend into soil, producing a network of branching roots of characteristic form, known as its root system architecture (RSA). RSA varies among species, and among varieties within a species that are adapted to different environments. Root traits are major targets for the second green revolution because of their potential to improve crop productivity, increase drought tolerance and nutrient acquisition, and increase C capture of soil. Improving the quality of roots in maize will be particularly valuable, since this crop is planted on over 92 million acres annually in the US. The future sustainability of agricultural systems relies on their ability to enhance soil organic matter (SOM) storage and reduce GHG emissions, while maintaining or enhancing productivity. This program had two components, Sensors and Models. For the first component, we designed and built a high-throughput phenotyping platform for root pulling of maize plants. This eliminated the physical labor of manually pulling up plants and reduced the number of personnel required down to one. The standardized pulling mechanism allowed recording force curves during the pulling process, providing additional information. We validated that the maximum force for pulling the root system was well-correlated with the root system mass and provided root crowns for further RSA analysis. These root crowns identified significant correlations with 2D root area and root depth, along with 3D root volume, total root length and number of root tips. We then used this system for field-based studies in maize on the genetics of root system architecture and its relation to nitrogen-use efficiency (NUE), including using lines relevant to the Corteva breeding program. Varieties were also evaluated at Corteva sites in the cornbelt and Danforth farm in Missouri, to establish responses across sites. From these studies we have identified genetic loci associated with root traits and created mutant lines for these loci and correlations of root traits with NUE. For the Models component, we worked to incorporate root and soil characteristics into the MEMS 2.0 soil and ecosystem biogeochemical model. Existing soil C models, such as Century, are unable to represent specific root trait interactions with the soil environment and therefore to accurately forecast the potential C sequestration benefits of root breeding under different climatic and soil type conditions. We have developed the MEMS 2.0 ecosystem biogeochemical model to improve quantification of farm-scale soil carbon and greenhouse gas emissions. The new knowledge and large datasets produced by this project will be used to develop and drive an innovative model capable of forecasting the impacts on soil C stocks and nutrient dynamics. An innovation was to use the empirical data from the field studies (in 1, above) to model genetic variation in nitrogen use efficiencies and soil C input. Our work demonstrated that maize root-derived C rapidly replaces existing soil C and after 3 years of continuous maize, up to 20% of soil organic C in the topsoil (0-15cm) and 3% in the subsoil (15-30cm) was contributed by maize. However, this contribution did not entirely represent a net increase. Root C contribution to soil was affected by maize genetics. We have analyzed soils derived from the CSU field trials for C and N stocks, in the different soil physical fractions represented by the MEMS model, using both physical fractionation with elemental analyses, and Fourier transformed infrared spectroscopy. Data will be used to link crop nitrogen use efficiencies with soil C sequestration and provide data to bridge the field trials with the model development, for verification of model predictions. The project had a number of successful outcomes: we have used the new phenotyping platform to identify new genetic loci that can enhance root phenotypes; we have partnered with multiple maize seed companies phenotype varieties in their breeding programs; we have developed the MEMS model that can help inform industry on the potential for carbon sequestration in the agricultural sector, and which is now available at the CSU Soil Carbon Solutions Center for use.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Requirements and Conceptual Design of Off-gas Systems for the Reprocessing of Metallic Fuels

An assessment has been conducted to determine how key regulations regarding volatile radionuclide emissions to the atmosphere may apply to the off-gas streams associated with electrochemical reprocessing. The scope of this assessment was based upon a generic electrochemical reprocessing scheme with a throughput rate of 200 MTIHM/y applied to metallic fuel discharged from a sodium fast reactor (SFR), but the findings are able to be translated to other advanced nuclear scenarios as merited. Air dispersion modeling was performed using the EPA CAP-88 model and evaluated the uncontrolled decontamination factors (DFs) that would be required to achieve regulatory compliance with the dose-based limits set forth by EPA regulation 40 CFR 190.10(a). These DFs were compared to those required by fuel cycle–based limits set forth by EPA regulation 40 CFR 190.10(b). Two theoretical sites with disparate climatological conditions were selected for air dispersion modeling (Idaho and Tennessee). The radionuclides modeled included 3 H, 85 Kr, 129 I, and selected alpha-emitting transuranic isotopes (referred to here as 239 Pu-TRU <1y ). It was found that the fuel cycle-based limits in 40 CFR 190.10(b) are most restrictive for 85 Kr and 239 Pu-TRU <1y , with DFs of 3 and 6.1E+09, respectively. The dose-based limit as derived from 40 CFR 190.10(a) could require mitigation of tritium in some scenarios, with an estimated DF of about 3 for the reference scenarios. The fuel cycle-based limit for 129 I resulted in a DF of about 240 for the reference scenario. The need for iodine mitigation based on dose to the public depended upon the physical form of iodine as either particulate or vapor-phase species. Emission of iodine from the facility as a vapor necessitated DFs of about 2 but emission as a particulate would require DFs >6,000 to meet thyroid dose-based limits. Effects of physical form on needed iodine mitigation are significant, but the understanding of speciation of iodine both during electrochemical reprocessing and after release to the atmosphere is limited. The electrochemical processing unit operations were evaluated to identify potential release points for the volatile radionuclides and to assess the potential for retention of the radionuclides within the process (thus decreasing the need for mitigation). Mitigation strategies for 3 H, 85 Kr, 129 I, and 239 Pu-TRU <1y were identified. In all cases, there are reasonably achievable pathways to regulatory compliance, although in some cases additional R&D is merited to verify the chemical speciation of these isotopes and to develop and demonstrate potential treatment technologies for this application. Whether or not additional off-gas controls (beyond common operations such as HEPA filtration and oxygen and moisture control) are needed for any of these regulated or volatile radionuclides depends on the (a) type of facility (NRC-regulated or DOE), (b) used fuel process rate, (c) used fuel burnup and composition, (d) speciation and retention of volatile radionuclides in the process and in the cell gas cleanup system, (e) site-specific parameters such as location, meteorology, stack height, and site boundaries, and (f) levels of conservatism and safety factors used in assessing compliance to air emissions regulations. Performance of this assessment revealed several areas where information is lacking or additional research is required in order to better determine if or what kinds of off-gas control might be needed. First, and most significantly, the understanding of the chemical speciation and physical form and partitioning of iodine during electrochemical processing operations is lacking and prevents the ability to accurately assess the potential iodine mitigation requirements. Future research in this area should be multifaceted and include thermodynamic modeling of iodine speciation in different process steps, experiments to quantify the kinetics of vapor-phase and melt-phase transitions, bench-scale experiments to determine the potential chemical and physical form of iodine emissions from the electrorefining process, and verification of iodine behavior with experiments utilizing operational facilities. Similarly, an improved understanding of iodine behavior in the environment after release from the facility stack will be required to refine dose estimations, as particulate and vapor-phase emissions result in significantly different doses to the MEI.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Toward engineering lattice structures with the material point method (MPM)

This study examines the potential of two variants of the material point method—the generalized interpolation material point (GIMP) and dual domain material point (DDMP) methods—in developing a robust computational framework for engineering lattice structures under different loading conditions. The study begins with assessing the ability of the two methods in predicting elastic buckling phenomena using column geometries with and without initial geometric imperfections. The results indicate that both methods effectively capture buckling phenomena when initial geometric imperfections are introduced. After this verification step, we create several models of tetrahedral lattice structures with varying strut diameter and orientation and subject them to quasi-static loading. We then validate the numerical results using laboratory test results. The results show that, while both methods accurately predict load–displacement curves in the pre-buckling regime, their predictive capabilities diminish in the post-buckling regime. Through visual comparison between the numerical and experimental deformed shapes, it appears that the discrepancies between model and experimental results are attributed to initial geometric imperfections in the lattices that occurred during 3D printing. We then establish a second set of lattice models where different types of initial geometric imperfections are considered. The results from these models show that imperfections have a negligible influence in the pre-buckling regime but affect the behavior considerably in the post-buckling regime. As a final step in this work, we subject the lattice models to impact loading and employ hypothetical soft and stiff materials. These results show that the lattice stiffness, which depends on material stiffness, strut diameter, and orientation, significantly influences the ability of a lattice structure to resist impact. In particular, we find that a stiffer lattice (i.e., one made with a stiff material and thicker struts) is capable of absorbing more energy than a softer one during impact. Although material nonlinearities, inelasticity, and detailed contact formulations are not considered in this study, the findings obtained herein lay the groundwork for engineering lattice structures under extreme loading conditions through a simulation-driven framework based on particle-based methods.

97 MATHEMATICS AND COMPUTING↗