SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
This process model is a conceptual process model that simulate the reactive captures of CO2 and the subsequent conversion of captured CO2 to methanol. The model was developed using Aspen Plus.
The future of materials design is rapidly advancing through the combination of quantum mechanics and data-driven modeling. These approaches integrate quantum principles with advanced data analysis, enabling precise insights into material behavior. This talk will highlight recent progress in using these methods for computational design, particularly in high-entropy alloy catalysts, emphasizing the role of hierarchical machine-learning architectures for accurate predictions. Additionally, I will discuss our work on developing machine learning interatomic potentials (MLPs) for single-element metals, metal oxides, and alloys under extreme conditions, focusing on melting behavior and phase properties at high temperatures and pressures. We have also refined our MLP models to capture dynamic surface interactions, such as CO2 and CO adsorption on MgO, using both static and molecular dynamics simulations. These models maintain high accuracy while significantly reducing computational costs compared to first-principles calculations. By enabling efficient and accurate simulations, this work supports broader community adoption, optimizes datasets for materials discovery, and extends the accessible time, size, and environmental conditions beyond the limits of experiments and traditional simulations.
In this work, a decay heat sensitivity analysis capability was developed and implemented in the ORIGEN code of the SCALE nuclear modeling and simulation suite. This capability introduces improved numerical integration schemes, which overcome the challenges associated with accurately modeling the behavior of adjoint nuclide amounts during coarse time steps for both nuclide amount and decay heat sensitivity calculations. This capability significantly improves the accuracy of calculations without compromising computational efficiency compared to the existing method. Extensive verification was conducted for various benchmark problems, including a 238 Pu decay and an irradiation problem involving 135 Xe, evaluated with both coarse and fine time grids. The results show excellent agreement with reference direct perturbation solutions, reaffirming the computational accuracy of the newly proposed numerical integration methods. Furthermore, sensitivity analyses were performed for fission product inventories ( 147 Sm, 150 Sm, 155 Gd) in pressurized water reactor UO 2 and MOX fuel assemblies. These analyses demonstrated that the ORIGEN sensitivity analysis capability can capture detailed sensitivity coefficients and underlying physics in real applications. Additionally, a decay heat sensitivity analysis for high-assay low-enriched uranium fuel, including various initial 235 U enrichment and burnup points, highlights the extended capabilities of SCALE/ORIGEN in comprehensively assessing the factors influencing total decay heat. These advancements in ORIGEN offer valuable insights for reactor analysis, fuel design, and safety assessments, especially in the context of advanced nuclear fuel development and design changes.
Transmission Electron Microscopy (TEM) is a powerful tool for the characterization of materials at the nanoscale; however, its inherent two-dimensional (2D) nature poses significant challenges to accurately measure three-dimensional (3D) properties. We introduce a supervised machine-learning model that predicts 3D structural information, such as sample thickness and curvature, from a series of conventional 2D TEM images. The model, a U-Net convolutional neural network, is trained on a large synthetic dataset generated from dynamical diffraction simulations that model TEM’s complex, nonlinear image formation, accounting for sample thickness and curvature. This physically realistic framework enables exploration of a broad parameter space impractical to sample experimentally. We demonstrate that the trained model has accurate predictions for experimental single-crystal silicon samples, achieving performance comparable to established measurement techniques. This work highlights the critical role of robust, simulation-based training in overcoming the limitations of real-world imaging artifacts and inconsistent sample geometries. By integrating machine learning with numerical simulations, we offer an efficient and scalable framework for quantitative TEM analysis, paving the way for more sophisticated 3D characterization of complex materials.
This project demonstrates, for the first time, the application of simulation-based inference (SBI) techniques to tune neutrino–nucleus interaction models. Using a mock dataset based on the MicroBooNE tuning of the GENIE event generator, our approach employs a Neural Posterior Estimator (NPE) with Masked Autoregressive Flows (MAF) to infer the posterior distributions of key GENIE parameters directly from simulated histograms. The workflow provides a scalable and amortized framework for performing likelihood-free inference in high-dimensional parameter spaces, offering a pathway to more efficient and uncertainty-aware model tuning for next-generation neutrino experiments such as DUNE and SBND.
This report implements a high-fidelity multiphysics modeling framework using the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program tools to track leaching and plating of materials between the fuel salt and structures in Molten Salt Reactors (MSRs). Specifically, this framework is demonstrated in modeling various processes within MSRs including hot and cold leg corrosion, redox potential driven corrosion, colloid precipitation and deposition, and noble metal plating and decay heat modeling. The model integrates neutronics, thermal-hydraulics, depletion, and thermochemistry to simulate the evolving chemical behavior of the fuel salt in interactions with structural materials. Initial findings suggest that hot-cold leg corrosion in MSR systems may become significant if the redox potential of the fuel salt is not sufficiently controlled. Additionally, the impact of noble metal decay heat deposition in large power reactors is explored. The resulting engineering framework for species tracking will be applied in the future toward specific validation work with ongoing experimental efforts focusing on both chloride and fluoride salt loop experiments with in-situ corrosion and redox control instrumentation.
We derive equivalent-circuit models of control- and physical-layer subsystems of grid-forming (GFM) inverter-based resources (IBRs) for electromagnetic-transient (EMT) simulations. Three different primary controllers are considered: droop, virtual synchronous machine (VSM), and dispatchable virtual oscillator control (dVOC). In addition, the models include cascaded voltage- and current-control loops, and LCL output filters. Simulations for a single-inverter setup and for a network of five inverters in a modified IEEE 14-bus topology are presented. The equivalent-circuit models simulated with analog electronic circuit-simulator software (in our case, LTspice) offer the same accuracy with up to 150x lower computational burden compared to block-diagram-based implementations in commercial off-the-shelf EMT software (in our case, MATLAB-Simulink).
Air humidity in indoor spaces plays a critical role in human comfort and health. Dehumidification systems are used for building humidity controls, but they can take significant energy consumption, especially in geographic locations with high outdoor humidity and warm climates. Consequently, there is a growing demand for innovative dehumidification processes that consume minimal energy. Dielectrophoretic air dehumidification represents one such promising approach. However, it has not garnered significant attention due to the absence of engineering models and simulation tools capable of evaluating its performance and limitations at large-scale airflows. A new numerical multiphase CFD model, which is also experimentally validated, is developed in a customized Reacting Foam solver based on OpenFOAM® version 9. The newly developed model seeks to decrease substantial energy consumption and lower costs by leveraging the dielectrophoretic phenomenon to regulate moisture levels in the air. The solver integrates a hybrid Eulerian-Lagrangian framework to track the droplet's trajectory and growth rate while solving the continuum equations for the moist air. An electrospray produces electrically charged droplets, which grow during their in-flight trajectories as water vapor condenses onto their surfaces. The role of electrostatic forces in promoting vapor condensation within a high-gradient electrical field is investigated, and the dielectrophoretic vapor nucleation process on charged water droplets is discussed. The CFD model was validated against results from the literature and from proof-of-concept experiments conducted by the authors, which showed a 2 % air dehumidification with a single electrospray and airflow rate of 5 cubic feet per minute. The simulation results indicated that augmenting the number of electrically charged spray droplets increased the dehumidification of the air to 25 %. The initial mean droplet diameter, the orientation of the injector and relative humidity significantly influence the assessment of dehumidification. As a result, scaling up this approach to larger airflow volumes is identified as a potential future research direction.
Triply Periodic Minimal Surfaces (TPMS) represent a promising geometry for future fuel designs due to their significant surface-to-volume ratio, which facilitates efficient cooling of nuclear fuel, a crucial factor for safety and efficiency. Demonstrating the remarkable capabilities of TPMS fuel requires initial modeling and simulation. This paper presents an implementation of TPMS in the Monte Carlo code OpenMC, enabling reactor physics modeling of TPMS. Here, the primary advantages over traditional methods using CAD files include reduced memory requirements for computations and high-fidelity implementation. This implementation has been tested against CAD files loaded in Serpent2, yielding promising results with low biases in the $k_{\textrm{eff}}$, comparable to biases in the material balance sheet. The implementation presented in this work will be used in future reactor physics computations related to new reactor designs involving TPMS-based fuels.
During high current operation, substantial heterogeneity develops within battery cathodes, particularly when their thickness is large. Heterogeneity relaxation during subsequent rest is important for understanding battery performance under pulsed conditions. Localized charge balancing phenomena within batteries at zero net current are not well understood and merit investigation. In this work, the heterogeneity within cathodes of commercial alkaline Zn–MnO 2 batteries is measured during discharge and monitored during rest using energy dispersive X-ray diffraction (EDXRD). Significant gradients in protonation form during discharge and partially relax under rest. It is demonstrated that the proton gradient relaxation is through local redox activity at zero net current, where local (de)protonation works to redistribute charge across the cathode thickness. To support this redox-based relaxation, a fundamental kinetic study on prismatic MnO 2 cathodes is conducted to determine an appropriate model to describe both discharge and charge kinetics of MnO 2 . These kinetics are incorporated into a computational model to simulate the proton gradient formation and partial relaxation under identical discharge conditions as the operando EDXRD experiments. Model and experimental data are found to be in excellent agreement, correctly predicting localized charge balancing at rest.
Nitrous oxide (N 2 O) is a greenhouse gas and stratospheric ozone-depleting substance with large and growing anthropogenic emissions. Previous studies identified the influx of N 2 O-depleted air from the stratosphere to partly cause the seasonality in tropospheric N 2 O (aN 2 O), but other contributions remain unclear. Here, we combine surface fluxes from eight land and four ocean models from phase 2 of the Nitrogen/N 2 O Model Intercomparison Project with tropospheric transport modeling to simulate aN 2 O at eight remote air sampling sites for modern and pre-industrial periods. Models show general agreement on the seasonal phasing of zonal-average N 2 O fluxes for most sites, but seasonal peak-to-peak amplitudes differ several-fold across models. The modeled seasonal amplitude of surface aN 2 O ranges from 0.25 to 0.80 ppb (interquartile ranges 21%–52% of median) for land, 0.14–0.25 ppb (17%–68%) for ocean, and 0.28–0.77 ppb (23%–52%) for combined flux contributions. The observed seasonal amplitude ranges from 0.34 to 1.08 ppb for these sites. The stratospheric contributions to aN 2 O, inferred by the difference between the surface-troposphere model and observations, show 16%–126% larger amplitudes and minima delayed by ~1 month compared to Northern Hemisphere site observations. Land fluxes and their seasonal amplitude have increased since the pre-industrial era and are projected to grow further under anthropogenic activities. Our results demonstrate the increasing importance of land fluxes for aN 2 O seasonality. Considering the large model spread, in situ aN 2 O observations and atmospheric transport-chemistry models will provide opportunities for constraining terrestrial and oceanic biosphere models, critical for projecting carbon-nitrogen cycles under ongoing global warming.
Abstract Several mechanisms have been proposed to explain why the isotope ratios of precipitation vary in space and time and why they correlate with other climate variables like temperature and precipitation. Here, we argue that this behavior is best understood through the lens of radiative transfer, which treats the depletion of atmospheric vapor transport by precipitation as analogous to the attenuation of light by absorption or scattering. Building on earlier work by Siler et al., we introduce a simple model that uses the equations of radiative transfer to approximate the two-dimensional pattern of the oxygen isotope composition of precipitation ( δ p ) from monthly mean hydrologic variables. The model accurately simulates the spatial and seasonal variability in δ p within a state-of-the-art climate model and permits a simple decomposition of δ p variability into contributions from gradients in evaporation and the length scale of vapor transport. Outside the tropics, δ p is mostly controlled by gradients in evaporation, whose dependence on temperature explains the positive correlation between δ p and temperature (i.e., the temperature effect). At low latitudes, δ p is mostly controlled by gradients in the transport length scale, whose inverse relationship with precipitation explains the negative correlation between δ p and precipitation (i.e., the amount effect). This suggests that the temperature and amount effects are both mostly explained by the variability in upstream rainout, but they reflect distinct mechanisms governing rainout at different latitudes. Significance Statement The isotopic composition of precipitation has long been used to make inferences about past climates based on its observed relationship with precipitation in the tropics and with temperature at higher latitudes. These relationships—known as the “amount effect” and “temperature effect,” respectively—have been attributed to many different mechanisms, most of which are thought to operate at either high or low latitudes but not both. Here, we present a unified framework for interpreting the isotope variability that can explain the latitude dependence of the temperature and amount effects despite making no distinction between high and low latitudes. Although our results are generally consistent with certain interpretations of the amount effect, they suggest that the temperature effect is widely misunderstood.
The MOOSE framework is a foundational capability used by the NEAMS program to create over 15 different simulation tools for advanced nuclear reactors. Due to this ubiquity, improvements to the framework in support of modeling and simulation goals are critical to the program. These improvements can take many forms including optimization, improved user experience, streamlined application programming interfaces (APIs), parallelism, and other new capabilities. The work transcribed in this report was conducted in direct support of the simulation tools and has already been deployed. The capabilities outlined in this report include enabling selective polynomial basis refinement, implementing a custom convergence system, building a scalable preconditioner for saddle-point problems, and much more.
The Delivery Environments Testbeds to Reduce Uncertainty in Simulations and Tests (TRUST) project is a broad project intended to analyze simplified problems experimentally and with modeling and simulation. The purpose of analyzing these simplified problems is to extend solution methods to more complex problems, as well as understand deficiencies and gaps in knowledge of methods currently used in more complex analyses. The TRUST project encompasses several smaller testbeds intended to isolate individual phenomena. The testbed under consideration in this report is the Sensors in Environments: Thermocouples testbed. In previous years, the purpose of this testbed was to quantify uncertainty of thermocouple sensors. To accomplish this, an aluminum plate was placed in a thermal chamber and subject to various types of thermal loading. Thermocouples were placed in various locations on the aluminum plate in various configurations (e.g., embedded in the plate, placed under Kapton tape), and an effort was made to quantify uncertainty in these measurements. Finite element simulations were performed to investigate how sensitive these measurements were to parameters such as the boundary conditions on the plate and material properties. However, a fundamental source of uncertainty in this analysis was the convective heat transfer from the plate. Convective heat transfer is a complex physical phenomenon comprised of a number of interacting sub-processes and is difficult to predict accurately a priori. As such, the main purpose of this testbed in FY25 was to better understand, both experimentally and numerically, the convective heat transfer from the plate. This is a highly applicable problem to several more complex problems, as convective heat transfer occurs in nearly all problems where a body is moving through air. Numerically, this required a two-step approach. First, the air flow in the thermal chamber was in vestigated using computational fluid dynamics. The commercial solver Fluent was used to perform these simulations. From these simulations, a heat transfer coefficient over the surface of the plate was calculated. This heat transfer was then used as boundary conditions for finite element heat transfer simulations within the plate, which were performed using Abaqus. Significant effort was devoted to automating the handoff between these two solvers. Experimentally, previous thermocouple results in the plate were used to validate the time-dependent thermal profiles produced from Abaqus. Further experimental efforts were performed both to help validate the Fluent simulations and to inform its boundary conditions. For example, hot-wire anemometers were used to measure the velocity in the chamber, which would be particularly useful in understanding the chamber inlet velocity. Thermocouple measurements were also taken in the chamber, instead of only on the plate, to serve as validation evidence for the Fluent simulations. Numerical results showed that the Fluent to Abaqus workflow matched previous plate thermocouple measurements well. This type of handoff is useful for more complex experiments, or those that are not able to be examined in as great of detail as this testbed, as it was performed without any experimental input. Experimental results, however, were more mixed. The anemometers proved unreliable, with inconsistent measurements across all anemometers, even at locations that were nearly identical. On the other hand, the thermocouples provided a relatively rich view of the temperature field in the chamber.
The Bering Sea shelf supports a highly productive marine ecosystem that is vulnerable to ocean acidification (OA) due to the cold, carbon-rich waters. Previous observational evidence suggests that bottom waters on the shelf are already seasonally undersaturated with respect to aragonite (i.e. Ω arag <1) and that OA will continue to increase the spatial extent, duration, and intensity of these conditions. Here, we use a regional ocean biogeochemical model to simulate changes in ocean carbon chemistry for the Bering Sea shelf from 1970–2022. Over this timeframe, model results suggest that surface Ω arag decreases by −0.043 per decade and surface pH by −0.014 per decade, comparable to observed global rates of OA. However, bottom water pH decreases at twice the rate of surface pH, while bottom [H + ] decreases at nearly 3 times the rate of surface [H + ]. This amplified bottom water acidification has emerged over the past 25 years and is likely driven by a combination of anthropogenic carbon accumulation and increasing primary productivity and subsurface respiration and remineralization. Due to this enhanced bottom water acidification, the spatial extent of bottom waters with Ω arag <1 has greatly expanded over the past 2 decades, along with pH conditions harmful to red king crab. Interannual variability in surface and bottom Ω arag , pH, and [H + ] has also increased over the past 2 decades, resulting in part from the increased physical climate variability. We also find that the Bering Sea shelf is a net annual carbon sink of 1.1–7.9 Tg C yr −1 , with the range resulting from the difference in the two different atmospheric forcing reanalysis products used. Seasonally, the shelf is a significant carbon sink from April–October but a somewhat weaker carbon source from November–March.
Energy sector planning models determine the aggregate need for new generation, but these models are typically at the state or regional scale and are not equipped to address the wide range of location- and technology-specific issues that are increasingly a factor in power plant siting. These animations demonstrate the aggregate siting suitability of various power plant technology configurations, considering technology-specific factors that can prohibit development. The data presented is from the GRIDCERF (Geospatial Raster Input Data for Capacity Expansion Regional Feasibility) data package. GRIDCERF is a harmonized, open-source geospatial product that can be used to evaluate siting suitability for renewable and non-renewable power plants in the conterminous United States. The animations presented here demonstrate a curated selection of the full suite of technology configurations available. GRIDCERF provides the necessary inputs for models that simulate power plant siting for regional capacity expansion planning such as the Capacity Expansion Regional Feasibility (CERF) model.
This report implements a high-fidelity multiphysics modeling framework using the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program tools to track isotopic species in Molten Salt Reactors (MSRs), with a specific focus on the 91-depletion chain within the Molten Salt Reactor Experiment (MSRE). The model integrates neutronics, thermal-hydraulics, depletion, and thermochemistry to simulate the production, transport, and phase transitions of isotopes under steady-state and transient conditions. The main findings reveal that isotopes such as bromine-91 largely remain in the liquid phase, while others, including krypton-91and yttrium-91, transition to the gas phase, significantly influencing the reactor’s radiological source term. The study also shows that during transients, like a reactivity insertion transient, rapid void formation and the expansion of the liquid-gas interface led to substantial transfers of dissolved isotopes into the gas phase, altering isotope distribution and largely increasing the source term in the off-gas system. Additionally, the research highlights that short-lived isotopes dominate the initial off-gas response during transients, while longer-lived isotopes determine the equilibrium state, underscoring the necessity of dynamic simulations for accurate species tracking and reactor safety analysis. The developed methodology will be applied in the future to the tracking of a larger number of species and introduce other species tracking mechanisms, such as deposition and plating.
The Material Protection, Accounting, and Control Technologies (MPACT) program utilizes modeling and simulation to assess Material Control and Accountability (MC&A) concerns for a variety of nuclear facilities. The Sandia National Laboratories (SNL)-developed Fissile Facility Flow Modeler (F3M) and the Material Accountancy Performance Indicator Toolkit (MAPIT) have historically provided MPACT with the capability to analyze MC&A approaches for nuclear facilities to determine that these facilities meet regulatory requirements. In FY25, improvements on the application of the F3M and MAPIT tools to simulate a generic TRi-structural ISOtropic (TRISO) fuel fabrication facility were successfully completed. The generic TRISO fuel fabrication F3M model captures the entire TRISO fuel fabrication process and is adaptable to any final TRISO fuel form, including spherical pebbles and cylindrical compacts loaded into graphite prismatic blocks. Comprehensive F3M/MAPIT functionality for the generic TRISO fuel fabrication model has been demonstrated. This modeling framework can be applied to support the U.S. Department of Energy and domestic nuclear industry stakeholders in developing MC&A approaches for advanced fuel fabrication facilities via statistical tests that demonstrate compliance to regulatory requirements.