Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Steady state process”

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

New directions and principles for solvent extraction for recovery of lithium from aqueous brines and mineral leachates: A brief review

Increasing demand for lithium for manufacturing of batteries is fueling the unprecedented search for improved recovery and alternative sources. Wider source distribution, lower energy consumption, and greater sustainability make extraction of lithium from brines, both natural and process-derived, an attractive alternative to mineral ores. Solvent extraction, used industrially for production of metals, salts, and pharmaceuticals, has been investigated as a methodology for lithium recovery for several decades. However, industrial application of solvent extraction for lithium recovery has so far been limited. In contrast, direct lithium extraction using adsorbents based on inorganic minerals has rapidly advanced from research to commercialization. A comparison of solvent extraction processes to adsorption highlights these issues and explains the preference for adsorbents. Although the application of solvent extraction has been criticized for use of large amounts of acid, alkali, and organic solvents, steady progress has been made to improve its potential for industrial lithium production, spurred on generally by the advantages of solvent extraction in selectivity and throughput. Previously developed beta-diketone, organophosphate, and crown ether ligands are being adapted and improved. Their novel use with ionic liquids, deep eutectic solvents, and membrane technologies promises to expand capabilities for extraction of lithium from dilute aqueous sources while improving sustainability. Possibilities for further discovery and innovation abound. In this review, we provide a unique perspective from the field of solvent extraction starting with fundamentals such as ion-transfer theory and apply them to understanding lithium selectivity and extraction behavior. In conclusion, the results are cast in the light of the practical realities of developing economical solvent extraction processes.

Brine↗
DOE OSTI2026

Public Reference Data for Megawatt-Scale Hydrogen Electrolysis - NLR Historical Solar PV

The U.S. Department of Energy and National Laboratory of the Rockies (NLR) demonstrate hydrogen electrolysis from variable sources, hydrogen compression and storage, and hydrogen fuel cell power production using megawatt-scale equipment at NLR’s Flatirons Campus as part of the Advanced Research on Integrated Energy Systems (ARIES) research platform. This dataset represents part of that effort and is intended for academic, national laboratory, industrial, and other stakeholders to plan, design, and validate models of megawatt-scale hydrogen technologies and diverse energy infrastructure nationwide. These data provide a baseline for how existing hydrogen electrolysis technologies perform when coupled with various energy technologies. Future datasets will demonstrate how existing hydrogen fuel cell technologies can provide controllable, dispatchable, and variable power output for artificial intelligence data centers and other variable loads. This dataset entry describes the behavior of a 1.25-MW proton exchange membrane MC250 electrolyzer system, manufactured by Nel Hydrogen , [1] when fed historical data generated by the 430-kW, fixed-axis solar photovoltaic (PV) array located at NLR’s Flatirons Campus. (While the electrolyzer balance of plant supports up to 2.5 MW of electrolysis, NLR only has a single 1.25-MW electrolysis stack.) Solar PV power output data for the 2020 calendar year were categorized on a daily basis by total energy generation and standard deviation. Each day was then ranked by these metrics, and the 25th, 50th, and 100th percentiles were selected. The 75th percentile day did not exhibit sufficient variability to make for a valuable experiment. A similar process was used for the related historical wind dataset . [2] The historical days in 2020 that represented these percentiles are Dec. 19, March 29, and May 4, respectively. The entire solar day’s power profile was then fed through the MC250 electrolyzer. Due to its length, the 100th percentile day experiment was split into two parts, and the final 3 hours of the solar day were not captured. These final 3 hours contained no spikes or dips of interest and simply represented a slow decay of input solar power. Also, a single timestamp (13:13:47 on Jan. 14, 2026) was lost in the hydrogen system supervisory control and data acquisition. Finally, during the 25th percentile experiment (solar day Dec. 19, 2020) data recording was lost from 11:00:13 to 11:14:45. The roughly 15 minutes of the solar profile were rerun at the end of the experiment and spliced into this time slot during post-processing. The electrolysis system controls hydrogen production by varying direct current applied to the stack, from a maximum of 3,000 A to a minimum safe operation of 300 A, or 10%. Because the current–voltage characteristic changes as the stack ages and efficiency degrades, the actual minimum safe operating power changes over time. The historical solar profiles were translated from power (kilowatts) to current (amperes) using a curve fit with calibration data and sent to the electrolyzer power supply at 1-Hz frequency. For more details on the statistical analysis process, see the slide deck “Public Reference Data for Megawatt-Scale Hydrogen Electrolysis: NLR Historical Solar PV Analysis and Profile Generation” accessible with this data entry. These datasets report relevant hydrogen balance-of-plant and system data, all captured at 1 Hz, including hydrogen mass production measured with an Emerson Coriolis flow meter. Each .zip file represents a single solar PV electrolysis experiment and is formatted as: {technology}_{percentile}_{scaling factor} For instance, “solarPV-430kW_25_2x.zip” reports the experiment using the 25th percentile solar data from the historical 2020 solar PV dataset, scaled to 200%. Scaling factors were applied to the generated solar PV power output files to more closely match the 1.25-MW capacity of the electrolyzer. Each .zip folder contains the following files: A .csv file containing raw data. An .xlsx file explaining all the fields in the raw data. A .png plot showing the time series of hydrogen production, electrolysis power consumption, and solar power input. A PDF file detailing the historical solar data statistical analysis used to generate the solar profile. An experiment labeled “characterization_200.zip” demonstrates the MC250 electrolyzer steady-state response with 30-minute load steps for a total duration of 5 hours. Finally, a .csv file is provided with all experiments combined into one dataset labeled "combined_solarPV_experiments.csv". [1] nelhydrogen.com/product/mc-series-electrolyser . [2] data.nlr.gov/submissions/316 .

08 HYDROGEN↗
DOE OSTI2025

Effect of Zr and Sc additions on coarsening- and creep resistance of AlSi10Mg fabricated by laser powder bed fusion

Microstructure and creep properties are studied in a eutectic AlSi10Mg alloy modified with Zr and Sc additions (Al-9.8Si-0.32Mg-0.70Zr-0.23Sc, wt%) manufactured through laser powder-bed fusion (L-PBF). Three types of Zr/Sc-bearing powders - elemental, master-alloy, and pre-alloyed - are employed in the fabrication process, with the pre-alloyed powders providing the highest incorporation of Zr and Sc in the Al matrix. Here, the as-printed alloy exhibits a fine cellular Al-Si eutectic structure which fragments and coarsens into micron-sized Si particles during aging at 300 ºC, leading to a steady drop in alloy microhardness between 0.1 and 1000 h. Coarsening of the eutectic Si phase during aging is not measurably affected by Zr and Sc in solid solution, which precipitate during aging and increases strength, compensating weakening from Si coarsening between 1 and 200 h at 300 ºC. Atom-probe analysis in the peak-aged condition (96 h/300 °C) confirms the presence of Al 3 (Sc,Zr) secondary nano-precipitates with an average radius of 1.0 nm and some Si solubility. Micron-sized grains are present throughout the alloy, with ultra-fine-grained regions at the melt pool boundaries, neither of which coarsen during long-term aging at 300 ºC. Under creep conditions at 300 °C, the Zr/Sc-bearing alloy with Al 3 (Sc,Zr) nano-precipitates exhibits power-law behavior, with a high apparent stress exponent (n a = 9) and a high threshold stress (σ th = 43 MPa), exhibiting nearly double the strength of a Zr/Sc-free AlSi10Mg control alloy (σ th = 22 MPa).

36 MATERIALS SCIENCE↗
DOE OSTI2020

Low temperature, ambient pressure electrochemical ammonia synthesis in alkaline media — Mechanistic studies and catalyst design. Final Report

Ammonia synthesis via the Haber-Bosch process is a pillar of modern agriculture, which converts the abundant but inert dinitrogen in the atmosphere to nitrogen-based fertilizers. Despite more than a century of optimization, the Haber-Bosch process remains energy intensive and reliant on fossil fuels and produces large amounts of CO 2 . Distributed and modular ammonia synthesis via the electrochemical nitrogen reduction reaction (ENRR) at or close to ambient conditions, powered by renewable electricity is an attractive alternative because it allows as needed production of ammonia, and in turn fertilizers. However, selective ENRR catalysts remain lacking. In this project, we systematically investigated two categories of catalytic materials, i.e., precious metals and transition metal nitrides, in the ENRR. Precious metals exhibit low activity and selectivity in the ENRR, which is in agreement with recent computational studies. Transition metal nitrides, especially VN, has been demonstrated as active and selective ENRR catalysts. We developed the first quantitative method to determine the amount of ammonia produced in the ENRR on N-containing containing catalysts in the membrane electrode assembly configuration. An initial ammonia production rate of 6.6 × 10 -10 mol s -1 mg -1 and a Faradaic efficiency of 6.0% at -0.1 V were reached on VN. After the initial deactivation, VN was shown to be stable in the ENRR for 116 hours with a steady state ammonia production rate of 1.1 × 10 -10 mol s -1 cm -2 . A vanadium oxynitride species was identified as the active phase in the ENRR via a combination of operando and ex situ characterization techniques. In addition, we established that the ENRR proceeded on transition metal nitrides via the Mars van Krevelen mechanism and determined the initial and steady state active site densities by developing a quantitative isotopic exchange method.

42 ENGINEERING↗
DOE OSTI2023

A comparative investigation of the effect of gas type on foam strength and flow behavior in tight carbonates

Foam generation technique has been practically applied to overcome gas mobility issue during enhanced oil recovery (EOR) processes. Incorporated with the utilized gas, the capability of foaming agent to generate persistent foam determines the flow performance under reservoir conditions. Here, in this study, we first identified the surfactants with excellent phase behavior and foaming properties and, second, evaluated their bulk foam performance with different gases, including N 2 , CO 2 , and hydrocarbon gas (CH 4 ) under high-pressure and high-temperature conditions. In bulk, the foam generated with the high-density gas, i.e., CO 2 , exhibited the greatest foamability but the lowest foam stability. With amidoamine oxide-based surfactant, CH 4 -foam was found to be more stable than N 2 -foam, while the opposite was observed for sulfobetaine-based surfactant. The most effective surfactant was later used to perform flow experiments on two tight carbonate rocks, including Minnesota Norther Cream Buff and Edwards limestone, with permeabilities of ~0.8 and 16 mD, respectively. The core-flooding experiments showed that the optimum quality of CO 2 -foam was lower than that of N 2 - and CH 4 -foams. For a given flow rate and fixed foam quality, the steady-state strength followed the descending order of CH 4 -, N 2 -, and CO 2 -foam. However, in general, the difference in foam strength between N 2 - and CH 4 -foams was insignificant when the total flow rate was varied. We also found that the variation in rock permeability had a greater effect on the strength of the CO 2 -foam compared to those generated with N 2 and CH 4 . In the presence of oil, the shear-thinning behavior of all foams was found more prominent. The increase in oil saturation had a more adverse effect on the strength of N 2 - and CH 4 -foams, particularly in the higher permeability rock. The results indicated that the foam strength, which is key in the effectiveness of foam EOR applications, is greatly influenced by several factors such as gas density, gas solubility in water, the presence of oil, and rock permeability.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Climate Warming Alters Nutrient Storage in Seasonally Dry Forests: Insights From a 2,300 m Elevation Gradient

Understanding potential response of forest carbon (C) and nutrient storage to warming is important for climate mitigation policies. Unfortunately, those responses are difficult to predict in seasonally dry forests, in part, because ecosystem processes are highly sensitive to both changes in temperature and precipitation. We investigated how warming might alter stocks of C, nitrogen (N), and phosphorus (P) in vegetation and the entire regolith (soil + weathered bedrock or “saprock”) using a space-for-time substitution along a bioclimatic gradient in the Sierra Nevada, California. The pine-oak and mixed-conifer forests between 1,160–2,015 m elevation have more optimal climates (not too dry or hot) for ecosystem productivity, soil weathering, and cycling of essential elements than the oak savannah (405 m) and subalpine forest (2,700 m). We found decreases in overstory vegetation nutrient stocks with decreasing elevation because of enhanced water limitation and greater occurrence of disturbances. Stocks of C, N, and P in the entire regolith peaked at the pine-oak and mixed-conifer forests across the bioclimatic gradient, driven by thicker regolith profiles and greater nutrient input rates. These observations suggest long-term warming will decrease ecosystem nutrient storage at the warmer, transitional pine-oak zone, but will increase nutrient storage at the colder, subalpine zone. Assuming steady-state conditions, we found the mean residence time of ecosystem C decreased with projected rising air temperatures and increased following a major drought event across the bioclimatic gradient. Our study emphasizes potentially elevation-dependent changes in nutrient storage and C persistence with warming in seasonally dry forests.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

SAS4A/SASSYS-1 Commercial Grade Dedication Example Report for a Generic Sodium Pool-Type Fast Reactor Application

In the U.S., a key component of the commercialization of advanced reactors is completion of a license application, which must ultimately be approved by the Nuclear Regulatory Commission (NRC). The approval of the license application by the NRC is contingent on satisfactory demonstration of the design basis and the response of the advanced reactor design to transient and accident scenarios using accepted codes and methods. This report describes the qualification and dedication requirements that the advanced reactor safety analysis system software SAS4A/SASSYS-1 are expected to need to fulfill to be used for sodium-cooled pool-type fast reactor licensing. The qualification and dedication requirements are identified through performance critical characteristics and evaluation model acceptance criteria representative of the advanced reactor design considered for licensing. This document captures, additionally, the verification process developed to demonstrate that the software fulfills the qualification and dedication requirements for a generic sodium-cooled pool-type fast reactor as part of the commercial grade dedication process. Like most software that has primarily existed in the research and development space, the most significant challenge facing SAS4A/SASSYS-1 for use in a licensing framework is the availability of a documentation basis describing the code pedigree. SAS4A/SASSYS-1 has been used for licensing of the fast flux test facility (FFTF) and the JOYO sodium-cooled fast reactor in Japan, as well as the design of the CRBR Plant. However, the historical verification and validation (V&V) activities supporting SAS4A/SASSYS-1 development do not align with modern software quality assurance (SQA) and V&V requirements. Two approaches to use of SAS4A/SASSYS-1 in a commercial licensing framework have been identified: commercial-grade dedication (CGD) and software qualification. The methods and requirements prescribed in the ASME NQA-1-2008/2009 Standard and Regulatory Guide 1.203 on the evaluation model development and assessment process (EMDAP) have been used as guidance to define the CGD and qualification processes, respectively. A qualification and dedication requirements matrix has been developed which utilizes fundamental software verification. In this process, software verification is defined as a software quality process aimed at defining software requirement specifications, developing software design documentation, and performing and documenting acceptance testing of the code against requirements. A key element of software qualification and dedication includes determination of software acceptance with respect to critical characteristics relevant to the functional requirements of the software. To assist with identification of cross-cutting transient phenomena and functional requirements, domestic SFR vendor designs have been reviewed to identify a reference SFR design. For this report, the reference design is defined as a pool-type reactor with metal alloy fuel, a liquid-metal intermediate heat transport system, and passive decay heat rejection systems. Given this reference, a series of high-level cross-cutting phenomena was identified for a general class of single-fault undercooling or reactivity insertion transients that scopes the design basis space, with the goal of assisting with prioritization of documentation development efforts for key transient models in SAS: 1) Reactivity feedback response prior to scram; 2) System-wide thermal inertia; 3) Transition in natural circulation flow regime in heat removal systems; 4) Decay heat generation; 5) Steady-state fuel characterization; 5) Clad/fuel behavior at elevated temperatures; 6) Point kinetics and decay heat; 7) Pump coastdown behavior; 8) Core flow redistribution in loss of forced convection; 9) Pool stratification. As a demonstration of CGD of SAS4A/SASSYS-1 for a sodium pool reactor, a software qualification and dedication gap analysis as it relates to code documentation has been performed. This effort leverages the framework established as part of the SAS4A/SASSYS-1 SQA Program. This CGD demonstration provides a framework that vendors can build upon to demonstrate the applicability of the SAS4A/SASSYS-1 software for licensing a sodium-cooled pool-type fast reactor.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

Rules for the crystallite size and dislocation density evolution in phases during α-ω transformation in Zr under high-pressure and severe plastic flow

The first in-situ X-ray diffraction (XRD) study of the evolution of the crystallite size and dislocation density in phases during plastic strain-induced phase transformation (PT) is performed utilizing α-ω PT in strongly pre-deformed commercially pure Zr as an important example. Rough diamond anvils (rough-DA) are introduced to intensify all occurring processes during heterogeneous compression of Zr in a diamond anvil cell (DAC). The main rule is found that during α-ω PT the crystallite size and dislocation density in ω-Zr depend solely on the volume fraction of ω-Zr and are independent of pressure, plastic strain tensor, its path, and initial nanostructure. Crystallite size in ω-Zr increases from 10 to 60 nm during the PT, while dislocation density reduces from 1.83×10 15 /m 2 to 0.6×10 15 /m 2 . Rough-DA produce a steady nanostructure in α-Zr before PT with smaller crystallite size and larger dislocation density than smooth-DA, leading to a reduction of the minimum pressure for α-ω PT to a record value 0.67 GPa, 9 times smaller than under hydrostatic loading and 5.1 times lower than the phase equilibrium pressure. In addition to strain, the kinetics of strain-induced PT unexpectedly depends on time. Also, strain-controlled part of kinetics is zero order, in contrast to the first-order kinetics with smooth-DA. The obtained results open a new window for understanding the mutual effects of nanostructure evolution and PT during severe plastic flow in various technological and natural processes. In conclusion, they may bring up economic strategies of producing nanocomposites and single-phase nanostructured materials with optimal properties.

Crystalline size↗
DOE OSTI2025

Alpha Decay Chains as Thermal Power Sources: Analysis and Applications for RTGs

Radioactive sources can provide power in remote and environmentally harsh locations such as the arctic or space. The generators powered by such sources are rugged and can withstand extreme temperatures, lack of sunlight, and require no human intervention for multiple years. Radioisotopes are used in thermoelectric generators to provide power at remote sites and deep in space. Isotopes like Pu-238, Cm-244, and Am-241 are used in these generators by NASA for power in space probes and spacecrafts. These power sources deliver a steady supply of energy over extended periods of time. Alpha particles created during decay do not travel far in a material. Their kinetic energy is transferred to heat that we can then convert into energy. Unlike beta and gamma decay, the slower-moving alpha particles stop in the material, making their energy available for use. Energy from these natural decay processes provides a reliable source of power. Spontaneous fission is rare and unreliable, and unlike induced fission processes, alpha decay occurs naturally and does not require external management or ignition. The ideal properties of an isotope for use as a power source depend upon the intended use. For use in an Arctic research base over a period of several years, but less than a decade, an isotope that provides high power output over a shorter lifespan may be the most suitable option. Whereas, for deep space missions where a consistent power source for decades or perhaps more than 100 years is needed that would require a very different isotope. One with a much longer half-life that would provide consistent power throughout that time and survive in that state in for these extended periods of time. These examples represent two extreme sides in terms of time frames. By analyzing the power produced by different radioactive decay processes over time, we can evaluate the suitability of various isotope decay chains for specific uses. Some unstable isotopes undergo a series of radioactive decays, transforming into different isotopes at each step and resulting in a stable isotope. The lists of isotopes in these decay processes are known as decay chains. Some of these chains, illustrated in the figures below, are currently being investigated for use in radioisotope thermoelectric generators (RTGs) designed for a range of operational durations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Scaling of atomic layer etching of SiO 2 in fluorocarbon plasmas: Transient etching and surface roughness

Fabricating sub-10 nm microelectronics places plasma processing precision at atomic dimensions. Atomic layer etching (ALE) is a cyclic plasma process used in semiconductor fabrication that has the potential to remove a single layer of atoms during each cycle. In self-limiting ideal ALE, a single monolayer of a material is consistently removed in each cycle, typically expressed as EPC (etch per cycle). In plasma ALE of dielectrics, such as SiO 2 and Si 3 N 4 , using fluorocarbon gas mixtures, etching proceeds through deposition of a thin polymer layer and the process is not strictly self-terminating. As a result, EPC is highly process dependent and particularly sensitive to the thickness of the polymer layer. In this paper, results are discussed from a computational investigation of the ALE of SiO 2 on flat surfaces and in short trenches using capacitively coupled plasmas consisting of a deposition step (fluorocarbon plasma) and an etch step (argon plasma). We found that ALE performance is a delicate balance between deposition of polymer during the first half cycle and etching (with polymer removal) during the second half cycle. In the absence of complete removal of the overlying polymer in each cycle, ALE may be transient as the polymer thickness grows with each cycle with a reduction in EPC until the thickness is too large to enable further etching. Small and statistical amounts of polymer left from a previous cycle can produce statistical variation in polymer thickness on the next cycle, which in turn can lead to a spatially dependent EPC and ALE roughness. Based on synergy between T i (sputtering time) and T p (passivation time), dielectric ALE can be described as having three modes: deposition, roughening surface (transitioning to etch-stop), and smooth surface with steady-state EPC.

Materials Science↗
DOE OSTI2021

Loss and Isotopic Fractionation of Alkali Elements during Diffusion-Limited Evaporation from Molten Silicate: Theory and Experiments

Moderately volatile elements (MVEs) are variably depleted in planetary bodies, reflecting the imprints of nebular and planetary processes. Among MVEs, Na, K, and Rb are excellent tracers for unraveling the history of MVE depletion in planetary bodies because they have similar geochemical behaviors but can be chemically fractionated by evaporation and condensation processes. Furthermore, K and Rb are amenable to high-precision isotopic analyses, which can help constrain the conditions of evaporation and condensation. To quantitatively understand why Na, K, and Rb are depleted in planetary bodies, we have carried out vacuum evaporation experiments from basaltic melt at 1200 and 1400 °C to study their evaporation kinetics and isotopic fractionations. We chose this composition because it is relevant to evaporation from small differentiated planetesimals. The Rb isotopic compositions of the evaporation residues were measured by multicollector inductively coupled plasma mass spectrometry (MC-ICPMS), and the K isotopic compositions were measured along profiles across the residues by secondary ion mass spectrometry (SIMS). In the 1400 °C run products, we found that the concentrations of both K and Rb in the run products decreased from core to rim, which was accompanied by a heavy K isotope enrichment near the surface. This indicates that, in this run, evaporation was limited by diffusion. To use those data quantitatively, we derive analytical equations that describe the evaporation rate and isotopic fractionation associated with diffusion-limited evaporation from a sphere, slab, and cylinder in transient and quasi-steady state regimes. This model is used to tease out the roles that diffusive transport in the melt and evaporation at the melt/gas interface play in setting the elemental depletion and isotopic composition of the residue. Under our experimental conditions, volatility decreases in the order of Na, Rb, and K. Using our experimental results in a thermodynamic model, we have estimated the product γΓ of activity coefficients × evaporation coefficients of Na, Rb, and K. The measured isotopic compositions of the residues are well explained using Rayleigh distillations, whereby the relative volatilities of K and Rb isotopes are given by the square root of their masses. We use our results and previously published data to predict how K and Rb could have been lost as a function of temperature, melt composition, oxygen fugacity, and saturation degree relevant to Vesta’s building blocks. We find that the K and Rb depletions, K/Rb elemental fractionation, and δ 41 K and δ 87 Rb isotopic fractionations of Vesta (as sampled by howardite-eucrite-diogenite (HED) meteorites) are best explained by evaporation of submillimeter size objects for 0.1-10 years at moderate temperatures (~1050 °C) in a medium ~98.8% saturated.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Network Uncertainty Quantification for Analysis of Multi-Component Systems

To impact physical mechanical system design decisions and realize the full promise of high-fidelity computational tools, simulation results must be integrated at the earliest stages of the design process. This is particularly challenging when dealing with uncertainty and optimizing for system-level performance metrics, as full-system models (often notoriously expensive and time-consuming to develop) are generally required to propagate uncertainties to system-level quantities of interest. Methods for propagating parameter and boundary condition uncertainty in networks of interconnected components hold promise for enabling design under uncertainty in real-world applications. These methods avoid the need for time consuming mesh generation of full-system geometries when changes are made to components or subassemblies. Additionally, they explicitly tie full-system model predictions to component/subassembly validation data which is valuable for qualification. These methods work by leveraging the fact that many engineered systems are inherently modular, being comprised of a hierarchy of components and subassemblies that are individually modified or replaced to define new system designs. By doing so, these methods enable rapid model development and the incorporation of uncertainty quantification earlier in the design process. The resulting formulation of the uncertainty propagation problem is iterative. We express the system model as a network of interconnected component models, which exchange solution information at component boundaries. We present a pair of approaches for propagating uncertainty in this type of decomposed system and provide implementations in the form of an open-source software library. We demonstrate these tools on a variety of applications and demonstrate the impact of problem-specific details on the performance and accuracy of the resulting UQ analysis. This work represents the most comprehensive investigation of these network uncertainty propagation methods to date.

42 ENGINEERING↗
DOE OSTI2025

Lower Snake River Dams Contribution to Grid Services

Hydroelectric generation and water storage have long been components of the clean energy mix, providing both reliable steady output and operational flexibility. As variable renewable energy sources such as wind and solar increasingly replace traditional generation, the role of all flexible resources—including hydropower—in balancing supply and demand continues to evolve. This study examined the contribution of the Lower Snake River (LSR) Dam plants to Bonneville Power Administration grid services in maintaining power system reliability within the Western Interconnection. By analyzing publicly available data, the study evaluated various reliability services through performance metrics including energy capacity, balancing and ramping, voltage and reactive power support, frequency response, and transmission impact. Results indicated that the LSR plants deliver services as expected based on their size, contributing to the balancing process, ramping capabilities, and operational reserves, particularly during peak load conditions and weather events, while also providing measurable frequency and voltage support to the grid.

13 HYDRO ENERGY↗
DOE OSTI2020

Design and analysis of a wake steering controller with wind direction variability

Wind farm control strategies are being developed to mitigate wake losses in wind farms, increasing energy production. Wake steering is a type of wind farm control in which a wind turbine's yaw position is misaligned from the wind direction, causing its wake to deflect away from downstream turbines. Current modeling tools used to optimize and estimate energy gains from wake steering are designed to represent wakes for fixed wind directions. However, wake steering controllers must operate in dynamic wind conditions and a turbine's yaw position cannot perfectly track changing wind directions. Research has been conducted on robust wake steering control optimized for variable wind directions. In this paper, the design and analysis of a wake steering controller with wind direction variability is presented for a two-turbine array using the FLOw Redirection and Induction in Steady State (FLORIS) control-oriented wake model. First, the authors propose a method for modeling the turbulent and low-frequency components of the wind direction, where the slowly varying wind direction serves as the relevant input to the wake model. Next, we explain a procedure for finding optimal yaw offsets for dynamic wind conditions considering both wind direction and yaw position uncertainty. We then performed simulations with the optimal yaw offsets applied using a realistic yaw offset controller in conjunction with a baseline yaw controller, showing good agreement with the predicted energy gain using the probabilistic model. Using the Gaussian wake model in FLORIS as an example, we compared the performance of yaw offset controllers optimized for static and dynamic wind conditions for different turbine spacings and turbulence intensity values, assuming uniformly distributed wind directions. For a spacing of five rotor diameters and a turbulence intensity of 10 %, robust yaw offsets optimized for variable wind directions yielded an energy gain equivalent to 3.24 % of wake losses recovered, compared to 1.42 % of wake losses recovered with yaw offsets optimized for static wind directions. In general, accounting for wind direction variability in the yaw offset optimization process was found to improve energy production more as the separation distance increased, whereas the relative improvement remained roughly the same for the range of turbulence intensity values considered.

17 WIND ENERGY↗
DOE OSTI2020

Two-Phase Fluid Flow Properties of Rough Fractures With Heterogeneous Wettability: Analysis With Lattice Boltzmann Simulations

Fractures are conduits for fluid flow in low-permeability geological formations. Multiphase flow properties of fractures are important in natural processes and in engineering applications such as the evaluation of the sealing capacity of caprocks and productivity of hydrocarbon-bearing tight rocks. Investigations of flow and transport through fractures typically focus on the effects of fracture geometric and mechanical factors such as aperture, roughness, and compressibility. The wettability of the fracture surfaces and its influence on microscale interfacial phenomena and macroscale effective transport properties are seldom studied. In this study, we investigated the effect of heterogeneous wetting properties on the displacement of water by supercritical CO2 through a series of lattice Boltzmann method simulations. The results show the evolution of the CO2 plume within a fracture is controlled by both the roughness of the aperture field and the wetting distribution. We combined these factors into a capillary pressure map that can be related to the macroscopic flow behavior of the fracture. We observed that heterogeneous wetting distributions promote the residual trapping of water where lower capillary pressures allowed for isolated water pockets in higher capillary pressure zones. Analysis of fracture unsteady relative permeability shows the effect of wetting on permeability evolution and provides support for the viscous-coupling relative permeability model. Finally, analysis of the steady-state relative permeability and saturation demonstrates a strong correlation between permeability and the standard deviation of the capillary pressure field. Thus, characterizing the distribution of wetting properties of fractures is crucial to understanding multiphase fracture flow and transport properties.

58 GEOSCIENCES↗
DOE OSTI2025

Data for Wilson and Megonigal (2025), "Nitrate reduction across soils transitioning from coastal forest to wetland are hotspots for denitrification"

Sea level rise drives spatial migration of coastal ecosystems and can lead to the accelerated replacement of coastal forests with tidal wetlands. Soil biogeochemical cycles in steady-state upland and wetland ecosystems are well studied, but pathways and rates in rapidly changing ecosystems are largely unconstrained. Wilson and Megonigal (2025) performed a one-time sampling and a subsequent incubation experiment, and characterized the reduction of reactive nitrogen (N) via denitrification and dissimilatory nitrate reduction to ammonia. Sampling was done at four sites where coastal deciduous forest is undergoing ecosystem state change and becoming wetland throughout the Chesapeake Bay, USA. The COMPASS-FME project (http://compass.pnnl.gov) established the sites sampled in this study in 2022–2023.This dataset consists of:* Isotope-labeled incubation results comparing nitrate reduction rates across transects spanning upland, transition, and wetland; and* Ancillary porewater chemistry data.All files in this dataset are plain text, comma-separated value (CSV), and no special software is required to read them.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Isotopic Substitution Reveals the Importance of Aluminate Diffusion Dynamics in Gibbsite (Al(OH) 3 ) Crystallization from Alkaline Aqueous Solution

Understanding molecular-scale factors governing the precipitation of aluminum hydroxides, such as gibbsite, under alkaline conditions is important for the formation of laterite deposits, as well as aluminum processing. However, mechanisms enabling tetrahedral aluminate ions to assemble into octahedral sites of the gibbsite lattice remain unclear. Formation of oligomeric complexes has been hypothesized as a critical intermediate step. Here, we report a study of gibbsite solubility in highly alkaline solutions using deuterium substitution to probe equilibrium and kinetic factors that could affect oligomeric intermediate formation, including the reactivity and the diffusivity of aluminate ions. When substituting sodium hydroxide with sodium deuteroxide, solution analysis shows a nearly 40% and 50% increase in gibbsite solubility in 2.4 and 3.3 mol·kg –1 total sodium solutions, respectively. Raman spectroscopy indicated that monomeric aluminate ions are the predominant species in solution irrespective of deuteration. However, both 27 Al and 23 Na nuclear magnetic resonance (NMR) spectroscopy revealed significant differences in chemical shifts in deuterated solutions, and analysis of 1 H, 23 Na, and 27 Al diffusion coefficients with pulsed-field gradient, stimulated echo NMR spectroscopy shows a decrease in ion diffusivity with increasing total deuterium in solution, consistent with the notion of increasing strength in the hydrogen/deuterium bonding network. In conclusion, because the relative change in 1 H diffusion coefficients are commensurate with the difference in apparent solubility constants, there is evidence for an oligomeric intermediate whose steady-state concentration is maintained by the collision frequency of aluminate ions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Nuclear Materials Packaging, Transportation, and Systems Analysis Group Software Quality Assurance Plan: ANSYS Mechanical Finite Element Analysis Software Version 2023R1

ANSYS Inc. develops and markets engineering simulation software and services used in the aerospace, automotive, manufacturing, electronics, biomedical, energy, defense, and many other industries. ANSYS is dedicated to engineering simulation and is the world’s leading software provider. ANSYS was founded in 1970 and is headquartered in Canonsburg, Pennsylvania. ANSYS provides an engineering analysis tool combining structural, thermal, computational fluid dynamics, acoustic, and electromagnetic simulation capabilities. ANSYS has two main programs, which use the same solvers: (1) Mechanical APDL (ANSYS Design Parametric Language), a Fortran-based coding platform, and (2) ANSYS Workbench, which uses a graphical user interface to aid in finite element analysis implementation. This plan covers both APDL and Workbench. The ANSYS computer program is a large-scale, multipurpose finite element program that can be used to solve several classes of engineering analyses. The analysis capabilities of ANSYS include the ability to solve static and dynamic structural analyses, steady-state and transient heat transfer problems, mode-frequency and buckling eigenvalue problems, static or time-varying magnetic analyses, and various types of field and coupled-field applications. The program contains many special features that allow nonlinearities or secondary effects such as plasticity, large strain, hyperelasticity, creep, swelling, large deflections, contact, stress stiffening, temperature dependency, material anisotropy, and radiation to be included in the solution. As ANSYS has been developed, other special capabilities such as substructuring, submodeling, random vibration, kinetostatics, kinetodynamics, free convection fluid analysis, acoustics, magnetics, piezoelectrics, coupled-field analysis, and design optimization have been added to the program. These capabilities contribute further to making ANSYS a multipurpose analysis tool for varied engineering disciplines. The ANSYS program has been in commercial use for over 50 years and has been used extensively in the aerospace, automotive, construction, electronic, energy services, manufacturing, nuclear, plastics, oil, and steel industries. Additionally, many consulting firms and hundreds of universities have used ANSYS for analysis, research, and educational purposes. ANSYS is recognized worldwide as one of the most widely used and capable programs of its type. Ansys design analysis software is the first created within a quality system with ISO 9001 certification, the internationally accepted quality standard. Product development, testing, maintenance and support processes also meet the United States Nuclear Regulatory Commission's quality requirements, as they have for nearly four decades. The Quality Assurance Service Agreement is suitable for the customers working in the nuclear industry who need to meet specific federal regulations including 10CRF50 Appendix B and provisions of 10CFR21. ANSYS has retained its original International Organization for Standardization (ISO) 9001 accreditation certificate since1995-05-04, It’s current certificate is valid until 2027-05-29.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗