Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “state dependent variable”

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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Impact of light output on the timing resolution of organic glass scintillator bars in a dual-ended readout configuration

The effectiveness of detector system modeling and validation depends on the quality of characterization performed on the system. For nuclear nonproliferation applications, which require detectors to address complex and variable field conditions, access to accurate system models is essential. This study presents the timing characterization of organic glass scintillator bars used in an imaging system developed at the University of Michigan. 137 Cs and 60 Co gamma-ray sources were used to determine the coincidence time resolution as a function of measured light output for two Compton scatter events between two organic glass scintillator bars. Timing resolution for two events, each with light output ranging from 100 to 1100 keVee were characterized. The resulting data were fit to a parametric function that was validated to agree with experimental results within ±25 ps. Simulations were performed to establish appropriate calibration points for each organic scintillator, thereby ensuring accurate calibration of light output values during analysis. This work successfully characterized the light output dependence of the coincidence timing resolution of a pair of organic glass detectors for use in an imaging system. Depending on the amount of scintillation light output from the events, the full width at half maximum of the measured coincidence time resolution ranged from 653.2 ± 16.1 ps for low light yields, 100 keVee, to 121.0 ± 8.4 ps for higher yields at 1100 keVee. The insights obtained from the timing resolution behavior will allow for accurate simulation capabilities in future security and verification efforts using scatter-based imaging systems.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Fully implicit crystal plasticity models representing orientations with modified Rodrigues parameters

Here, this work describes a crystal plasticity formulation combining several mathematical, numerical, and implementation choices to produce a highly efficient model. Specifically, the key choices in the implementation are (1) representing orientations with modified Rodrigues parameters, (2) implementing a fully coupled implicit time integration for the elastic stretch, the crystal orientations, and the model internal variables, (3) implementing the model in the NEML2 constitutive modeling framework, based on PyTorch, to vectorize the calculations and port the computation to GPUs and other hardware accelerators, and (4) an exact implementation of the consistent tangent matrix, even for arbitrary coupling to other field variables beyond the displacements, like temperature, neutron fluence, etc. The first two features of the model are, to our knowledge, novel. The paper considers each of these choices individually as well as the final model as a whole. This includes a full description of modified Rodrigues parameters, their advantages over other representations of orientations, the mathematical formulae and tools required to implement a model with modified Rodrigues parameters, and a detailed description of the geometry of the space of modified Rodrigues parameters (in an appendix). It also includes a description of a fully implicit time integration scheme for the orientations and the advantages in representing orientations with modified Rodrigues parameters in implementing such a model. The work then assess, via numerical examples, the advantages of fully coupled implicit time integration versus more common decoupled and explicit time integration schemes. These studies demonstrate the computational advantages of fully coupled integration versus other time integration algorithms, though the performance of the competing models depends on the complexity of the underlying single crystal model. The study concludes by demonstrating that the choice of time integration method affects the sharpness of the predicted texture, with explicit methods for integrating the orientations overestimating texture sharpness and implicit methods underestimating texture sharpness.

Crystal plasticity↗
DOE OSTI2025

Factors That Influence Variability in Stress-Drop Measurements Using Spectral Decomposition and Spectral-Ratio Methods for the 2019 Ridgecrest Earthquake Sequence

Stress drop is a fundamental parameter related to earthquake source physics, but is hard to measure accurately. To better understand how different factors influence stress-drop measurements, we compare two different methods using the Ridgecrest stress-drop validation data set: spectral decomposition (SD) and spectral ratio (SR), each with different processing options. Here, we also examine the influence of spectral complexity on source parameter measurement. Applying the SD method, we find that frequency bandwidth and time-window length could influence spectral magnitude calibration, while depth-dependent attenuation is important to correctly map stress-drop variations. For the SR method, we find that the selected source model has limited influence on the measurements; however, the Boatwright model tends to produce smaller standard deviation and larger magnitude dependence than the Brune model. Variance reduction threshold, frequency bandwidth, and time-window length, if chosen within an appropriate parameter range, have limited influence on source parameter measurement. For both methods, wave type, attenuation correction, and spectral complexity strongly influence the result. The scale factor that quantifies the magnitude dependence of stress drop show large variations with different processing options, and earthquakes with complex source spectra deviating from the Brune-type source models tend to have larger scale factor than earthquakes without complexity. Based on these detailed comparisons, we make a few specific suggestions for data processing workflows that could help future studies of source parameters and interpretations.

58 GEOSCIENCES↗
NASA NTRS1999

A View of the Lunar Interior Through Lunar Laser Range Analysis

Laser ranges between observatories on the Earth and retroreflectors on the Moon started in 1969 and continue to the present. Recent range accuracies are 2 cm while earliest ranges are an C, order of magnitude less certain. Four retroreflectors are ranged: three located at the Apollo 11, 14, and 15 sites and one on the Lunakhod 2 rover. Accurate analysis of the range data determines a number of lunar science parameters. The lunar interior variables include a fluid core parameter. The Lunar Laser Ranging effort is reviewed elsewhere. Many parameters are detected through their influence on rotation. Also detected are solid-body tides and accurate selenocentric reflector locations. Determined through the rotation are moment-of-inertia differences, gravitational harmonics, potential Love number, and dissipation effects due to tides and molten core. The rotation of the Moon is not at its minimum energy state; some recently active process has caused free librations. The moment differences contributed to the recent improvement of the Moon's moment of inertia from the Lunar Prospector gravity field. The Love numbers provide bulk elastic properties. Future possibilities for measurement include oblateness of the core-mantle boundary and core moment. A study of dissipation signatures in the rotation determines tidal Q vs. frequency and concludes that the Moon has a molten core. At 1 month the tidal Q is 37 and at 1 yr it is 60. The core radius is < or = 352 km for Fe and < or = 374 km for the Fe-FeS eutectic. The core detection exceeds 3x its uncertainty. The spin of the core is not aligned with the spin of the mantle and torque arises from the velocity difference at the boundary. Yoder's turbulent boundary layer theory is used to compute the radii. The present heat generation from tides and core interaction is minor compared to radiogenic heating. The heating for ancient times is more interesting. Peale and Cassen investigated lunar tidal heating while the lunar orbit expanded due to tides on Earth. Their calculations predate the measurement of Q and should be multiplied by 3.45 to match the lunar-laser-determined Love number and monthly Q. Tidal-heating computations depend on how fast the lunar orbit evolved and whether the tidal dissipation is localized. Neither is known, but plausible assumptions lead to early central region temperature increases of several hundred degrees. Most of the energy is deposited early in the Moon's history. The turbulent boundary layer theory allows a prediction of energy dissipated at the core-mantle boundary during orbit evolution. Under the assumption that the properties of the early core are the same as at present, the energy dissipated by core-mantle interaction is about the same as for tidal dissipation, but it is deposited in a smaller volume. This source of energy is capable of promoting convection in an early fluid core and driving a dynamo. This is a transient phase with a duration depending on the rate of orbit evolution. Plausible assumptions lead to a duration of a few hundred million years. Thus, the remnant magnetization of many lunar rocks is compatible with a brief global magnetic field powered by dynamical energy dissipation. Analysis of the lunar laser ranges is providing information on lunar geophysics. Future data will improve accuracies of present solution parameters, and several more interior effects should be detectable.

Williams, J. G.↗
DOE OSTI2026

A multi‐variable framework for selecting WRF physics configurations at convection‐permitting scales: An Amazon wet‐season case study

Tropical convection over rainforests modulates atmospheric circulation and the energy and hydrological cycles across multiple scales. However, the scarcity of observations still limits our understanding of these processes. Although numerical models are utilized to investigate atmospheric physical processes, their performance depends on the choice of parameterizations and grid resolution. Here, in this work, we introduce and apply a multi‐variable, multi‐physics framework to evaluate and rank the Advanced Research Weather Research and Forecasting (WRF–ARW) configurations at 1‐km resolution over the central Amazon during the wet season. A 48‐member ensemble combines three land‐surface models (LSMs), four planetary boundary layer (PBL), and four microphysics (MP) schemes. Model performance for seven convection‐related near‐surface and boundary‐layer variables is assessed using Taylor diagrams and the Taylor skill score (TSS), analysis of variance (ANOVA)‐based sensitivity metrics, and non‐parametric rank tests. We then construct a combined, weighted TSS to identify configurations that are comparatively robust across variables. Results showed that most configurations reproduce near‐surface temperature, sensible and latent heat fluxes, and boundary‐layer height reasonably well, whereas humidity and rainfall remain challenging. LSM choice has the strongest impact on the surface fluxes and a secondary influence on near‐surface temperature and humidity, PBL schemes dominate boundary‐layer height, and MP schemes exert the largest control on rainfall. No single configuration is optimal for all variables, but the combination of Noah (LSM), Yonsei University (PBL), and Morrison (MP) emerges as the most robust configuration for this case, with the WRF single‐moment six‐class scheme (WSM6) providing a competitive, computationally cheaper MP alternative. The framework is general and can be applied to other regions, seasons, and convective regimes.

Amazon rainforest↗
DOE OSTI2026

Modulation induced by very-large-scale motions on the inclination angle of wall-attached eddies: an atmospheric surface layer experiment

The forward leaning inclination angle, 𝛾, of coherent turbulent structures is a well-known feature of wall-bounded turbulent flows. Although invariant across friction Reynolds numbers within the range 𝑅𝑒 𝜏 =10 3 −10 6 , 𝛾 can vary significantly across turbulent scales within a high-Reynolds-number flow. Very-large-scale motions (VLSMs) are known to induce significant changes in the instantaneous shear profile, which is a conditioning event that could trigger variability in the inclination angle of smaller coherent turbulent structures. Although this aspect has been extensively studied via numerical and laboratory experiments, few studies have explored this feature for a very-high-Reynolds-number atmospheric flow. In this work, the inclination angle of turbulent structures within the atmospheric surface layer at a very high Reynolds number (𝑅𝑒 𝜏 =7.9 ×10 5 ) is investigated by deploying a scanning Doppler light detection and ranging and a super large particle image velocimetry (SLPIV) apparatus. The inclination angle of wall-attached eddies is inferred either from the two-point correlation of streamwise velocity (𝛾 =41.1∘) or with a scale-dependent approach through the spectral linear stochastic estimator (SLSE). The SLSE (and, thus, the scale-dependent inclination angle) is conditionally evaluated based on the high- and low-momentum events induced by VLSMs, both in the streamwise (𝑢′𝑉𝐿𝑆𝑀) and in the vertical (𝑤′ 𝑉𝐿𝑆𝑀) velocity components. As a result, lower inclination angles (𝛾 =30° −50°) are found for 𝑢$^{'}_{VLSM}$ >0 (𝑤$^{'}_{VLSM}$ <0), while higher values (50° −85°) are ascribed to 𝑢$^{'}_{VLSM}$ <0 (𝑤$^{'}_{VLSM}$ >0). This result emphasises the primary role that VLSMs play in shaping the wall-attached eddy geometry, which, in turn, is crucial to determine the Reynolds stress balance within the wall-attached eddy range.

42 ENGINEERING↗
DOE OSTI2025

Explaining drivers of housing prices with nonlinear hedonic regressions

Housing markets play a critical role in shaping the spatial and demographic evolution of urban areas. Simulating housing price dynamics can enhance projections of future urban development outcomes. However, traditional hedonic regressions for housing prices, which neglect nonlinear interactions among explanatory variables, often exhibit limited predictive performance. While machine learning (ML) methods can provide a more flexible representation of the relationships between predictors, they are often regarded as “black boxes” due to their complexity and lack of transparency. Interpretable ML techniques provide a promising route by combining the flexibility of ML methods with approaches to analyze the relationships between inputs and outputs. In this study, we employ interpretable ML to analyze the patterns driving the housing market in Baltimore, Maryland, USA. We train an Artificial Neural Network (ANN) to predict Baltimore housing prices based on structural characteristics (e.g., home size, number of stories) and locational attributes (e.g., distance to the city center). We then conduct sensitivity and Partial Dependence Plot (PDP) analyses to interpret the fitted ANN model. We find that the ML model achieves higher predictive accuracy and explains 16 % more of housing price variance than a traditional linear regression model. The interpretable ML model also reveals more nuanced and realistic nonlinear relationships between housing sales price and predictors as well as interactive effects underlying Baltimore home price dynamics. For instance, while the linear model indicates a steady housing price increase over time, our interpretable ML model detects a post-2008 decline, with smaller properties experiencing the sharpest drop.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Variable Surface Termination and Ligand Passivation of Lead Sulfide Nanocrystals Synthesized with Excess Lead Chloride

The surface termination and ligand passivation of semiconducting nanocrystals (NCs) impact the stability, optical properties, and self-assembly of NCs. In this work, we definitively characterize the surface of lead sulfide (PbS) NCs synthesized from excess PbCl 2 . With a combination of small-angle neutron scattering (SANS), photoluminescence, and 1D and 2D NMR experiments, we show that the surface termination of PbS NCs depends on the presence of PbCl 2 during ligand exchange. When excess PbCl 2 is removed prior to ligand exchange, PbS[RNH 3 + Cl – ] NCs are obtained, which are terminated by a monolayer of PbCl x on the {100} PbS facet and passivated by oleylammonium chloride ligands. On the other hand, when excess PbCl 2 remains in solution during ligand exchange, lead oleate forms and attaches to the {111} PbS facets of PbS[RNH 3 + Cl – ] NCs, creating PbS@PbCl x NCs. PbS@PbCl x NCs are coated in an epitaxial layer of PbCl x on both the {100} PbS and {111} PbS facets, making them 0.3–0.4 nm larger on average than PbS[RNH 3 + Cl – ] NCs with identical absorption wavelengths. Additionally, PbS@PbCl x NCs have a consistently higher photoluminescence quantum yield and longer photoluminescence lifetimes than PbS[RNH 3 + Cl – ] NCs. This study clarifies the surface structure of PbS NCs synthesized from excess PbCl 2 , highlighting ligand exchange strategies and reconciling observations from across the literature.

Anions↗
NASA NTRSDate not supplied

Exploring Radiation, Water and Energy Closure in Key Deep Tropical Regions Over the Last 3 Decades

The closure of global energy and water budgets at various time scales depends heavily upon the on the processes of water and energy exchange in the deep tropics. A key process observed is the change in convection regimes from disorganized convection to organized deep convection. During these periods an imbalance in key tropical regions is observed of water budget terms between observations and atmospheric reanalysis systems, implying that processes are not simulated well. For instance, periods of imbalance in the water budget may be at least partly attributed to precipitation estimates due to the transition between cloud and precipitation systems in which ice cloud affects microwave radar estimates (Lietmann-Niimi et al. 2024, submitted). A recent paper by Hsaio et al (2023, pre-published version) found that the transition process is characterized by wind shear and longwave cloud radiative feedback signatures as well. This is just one example how these processes also affect the radiative fluxes and thus the energy budget at the TOA. In this presentation, we assess the net radiative budget and components in the deep tropics by comparing multiple data products (GEWEX SRB Rel4IP, CERES EBAF/SYN1deg Ed4.2, ISSCP FH, etc.) and describing the temporal variability of energy balance terms across the deep tropics for the period from 1988 to near present in six key oceanic regions. These six key tropical oceanic regions are West Indian Ocean, the south Indian Ocean east, tropical western, central and eastern Pacific Ocean, and Tropical Atlantic Ocean. In these regions we assess the sensitivity of the fluxes as a function of fluctuations in cloud types in response to various larger scale atmospheric processes (e.g. El Niño, Indian Ocean Dipole). We assess this variability of the radiative flux anomalies (including the net TOA, surface and the atmospheric divergence fluxes) against water vapor divergence, cloud properties (include ISCCP “Weather States, Tselioudis et al. 2021) and larger scale wind shear to assess the uncertainties of energy closure in these regions. These are contrasted against the modeled radiative fluxes from ERA-5 and MERRA-2. Ocean buoy measurements of radiative fluxes are utilized to help assess data radiative flux uncertainty over the nearly 40 year period. In general, the variability for overlap periods of these various data products agrees well, but there are significant differences in the net fluxes that vary according to the rendering of surface, atmospheric, aerosol and cloud properties. These differences also affect closure terms and inferences on systematic uncertainties within the observational estimates. We conclude with recommendations for continuing work

Radiation budget↗
DOE OSTI2025

Ensemble cure kinetics network (ECK-Net): A method to derive cure kinetics of thermosetting resin

This paper introduces an Ensemble Cure Kinetics Network (ECK-Net), a neural network (NN)–based framework for modeling the cure kinetics of thermosetting resins within a phenomenological context. ECK-Net replaces traditional analytic models, which require extensive chemical insight and multiple isothermal/non-isothermal experiments, with a data-driven surrogate that maps nonlinear relationships between temperature, degree of cure, and reaction rate from differential scanning calorimetry data. The proposed approach predicts input-dependent kinetic coefficients of a generalized nth-order reaction equation rather than reaction rates directly, enabling a single unified model to represent various epoxy systems without relying on iso-conversional analysis or predefined functional forms. To ensure robustness, multiple independently trained networks under different random initializations are blended through an ensemble strategy, effectively mitigating the stochastic variability inherent to neural networks. The framework is validated using experimental datasets from multiple resin systems, including aerospace-grade materials (Toray 3900-2, Cycom 5320-1, and Hexcel 8552) and a windmill-grade resin (RIMR 035c). The model accurately reproduces the temporal evolution of the degree of cure under manufacturers’ recommended cure cycles across all tested resins systems, yielding Pearson’s correlation coefficients of 0.992, 0.994, 0.993, 0.997, respectively. To demonstrate process-level applicability, the trained network was implemented within the Abaqus environment to simulate out-of-autoclave (OOA) curing process of the CFRP panel composed of Toray T830H-6K/3900-2D prepreg. The simulation results showed excellent agreement with experimental temperature response (maximum peak temperature, simulation: 189.6 °C, experiment: 188.5 °C) and the final degree of cure (simulation: 0.948, experiment: 0.960 ± 0.013), confirming ECK-Net’s capability as a reliable alternative to conventional cure kinetics modeling methods.

Composite curing↗
DOE OSTI2026

Constraining the Aerosol Effects on Deep Convective Clouds by Considering the Coupling Between Clouds and the Planetary Boundary Layer

Several mechanisms have been proposed for the aerosol invigoration effect. Although their principles are well established, their actual magnitudes and roles in cloud development remain uncertain and debatable. This uncertainty partly stems from observational-based studies, in which it has been challenging to separate the co-variability between aerosols and meteorology. Addressing this problem requires large data samples. To this end, this study employs the Atmospheric Radiation Measurement data set expanding to 16 years (some up to 17 years, compared to 10 years in previous work) in the U.S. Southern Great Plains. It also conducts a more careful and rigorous analysis to isolate the influences of convective available potential energy (CAPE) and synoptic patterns to address a previously raised concern. We incorporated a new key process affecting aerosol-cloud interaction: cloud-surface coupling. The state/degree of the coupling relationship turns out to play an important role in the invigoration effect. Our analysis reinforces earlier findings of a robust positive relationship between cloud thickness and aerosol loading across CAPE percentiles—but only under cloud-surface coupled conditions. The increase in cloud thickness with aerosol loading is most pronounced in coupled clouds with high CAPE and bases below 1 km. Coupled clouds with bases below 1 km thicken between 1 and 4 km, depending on the CAPE percentile. Decoupled clouds show no such systematic changes. Synoptic patterns also lead to different strengths of the invigoration effect. Clean and polluted air masses are predominantly associated with northerly and southerly winds, respectively, with a stronger invigoration effect in cleaner air masses.

Geosciences↗
NASA NTRS1962

BOILING HEAT TRANSFER TO LIQUID HYDROGEN AND NITROGEN IN FORCED FLOW

Boiling heat transfer to liquid hydrogen and nitrogen was investigated experimentally. Results are presented from a study of bulk boiling inside a cylindrical tube under vertically upward forced-flow conditions. A 0.555-inch-inside-diameter and 16 (1/8)-inch-long electrically heated stainless-steel tube was used. The range of variables studied for hydrogen were mass velocity of 2850 to 17,000 pounds per hour per square foot, local heat flux of 3600 to 40,000 Btu per hour per square foot, inlet pressure of 30 to 74 pounds per square inch absolute, and inlet subcooling of 0° to 9° R. Nitrogen test conditions were mass velocity of 15,000 to 56,000 pounds per hour per square foot, local heat flux of 2300 to 40,000 Btu per hour per square foot, inlet pressure of 47 to 56 pounds per square inch absolute, and inlet sub cooling of 1° to 6° R. The axial distribution of the tube-wall temperatures is presented. A transition in the type of boiling heat transfer was obtained. The critical heat flux corresponding to this transition was determined over a range of flow and heating rates and local qualities. At specific combinations of flow and transition location, a range of critical-heat-flux values was obtained and maximum values were determined. The maximum critical heat flux increased with increasing fluid-flow rate and decreased with increasing length of tube before transition. Similar variations of the maximum critical heat flux have been reported for water. The tube-inner-wall temperatures upstream of transition were essentially uniform and were only slightly greater than the fluid saturation temperature. The Wall-temperature profiles downstream of transition generally resembled those obtained in film-boiling studies and appeared to be strongly dependent upon local quality at the point of transition. Maximum wall temperatures of 9000 and 18000 R were obtained with hydrogen and nitrogen, respectively. Fluctuations of pressure, flow rate, and temperature occurred during some of the boiling tests. Under some conditions, maximum critical-heat-flux values were attained during steady-state operation with fluctuations. In other cases the fluctuations became uncontrolled, and critical-flux values less than the maximum values were obtained upon restabilization of the test conditions. No measurable pressure drop across the test section was obtained at any condition.

LIQUID HYDROGEN↗
NASA NTRS2004

A Meteorological Model's Dependence on Radiation Update Frequency

Numerical weather models are used to simulate circulations in the atmosphere including clouds and precipitation by applying a set of mathematical equations over a three-dimensional grid. The grid is composed of discrete points at which the meteorological variables are defined. As computing power continues to rise these models are being used at finer grid spacing, but they must still cover a wide range of scales. Some of the physics that must be accounted for in the model cannot be explicitly resolved, and their effects, therefore, must be estimated or "parameterized". Some of these parameterizations are computationally expensive. To alleviate the problem, they are not always updated at the time resolution of the model with the assumption being that the impact will be small. In this study, a coupled land-atmosphere model is used to assess the impact of less frequent updates of the computationally expensive radiation physics for a case on June 6, 2002, that occurred during a field experiment over the central plains known as International H20 Project (IHOP). The model was tested using both the original conditions, which were dry, and with modified conditions wherein moisture was added to the lower part of the atmosphere to produce clouds and precipitation (i.e., a wet case). For each of the conditions (i.e., dry and wet), four set of experiments were conducted wherein the model was run for a period of 24 hours and the radiation fields (including both incoming solar and outgoing longwave) were updated every 1, 3, 10, and 100 time steps. Statistical tests indicated that average quantities of surface variables for both the dry and wet cases were the same for the various update frequencies. However, spatially the results could be quite different especially in the wet case after it began to rain. The near-surface wind field was found to be different most of the time even for the dry case. In the wet case, rain intensities and average vertical profiles of heating associated with cloudy areas were found to differ for the various radiation update frequencies. The latter implies that the mean state of the model could be different as a result of not updating the radiation fields every time step and has important implications for longer term climate studies

Eastman, Joseph L.↗
DOE OSTI2021

Digital Twin Model for Advanced Manufacture of a Rotating Detonation Engine Injector

A digital twin material model (DTMM) of an additive manufacturing (AM) process was created to advance the state of the art in rotating detonation engine (RDE) injector design. Current RDE injectors are designed with large pressure drops, enabling a stable and repeatable combustion process. However, this comes at the cost of system efficiency. For the technology to transition to commercial fossil-based power generation, it is important to develop injectors with reduced flow losses. Low-loss injectors are difficult to design and manufacture with conventional manufacturing techniques. AM enables new design options, but the AM manufacturing process must be thoroughly understood to result in a robust design. A DTMM provides the necessary insight by defining the cause-effect relationships between process parameters, microstructure features, and properties. Therefore, a DTMM to support the design and manufacturing process was developed and applied to the design of a new additively manufactured low-loss injector. The injector combustion behavior was characterized through hot-fire tests, and mechanical performance was compared to the DTMM predictions. The two project goals were the successful development of the DTMM and the demonstration of an improved RDE injector design. The RDE injector design and DTMM developments occurred on parallel but dependent paths. The injector was designed to reduce pressure drop by increasing the cross-sectional flow area ratio between the injector air passages and the combustor annulus. This resulted in less structural material, raising the concern that thin members would be susceptible to high-cycle fatigue (HCF) under the periodic loading inherent to an RDE. It was most important for the DTMM to predict behavior in these features; therefore, the injector design concept guided the material thicknesses used in fatigue tests. The DTMM development started by manufacturing a series of coupons over the range of possible AM process variations. A design-of-experiment approach was used to select which process variable combinations gave the most efficient coverage relevant to the injector design space. The microstructure in each of these coupons was characterized, and then computational methods were used to create a numerical model of the correlation between process variables and microstructure. Next, a set of HCF samples were tested to calibrate existing models that map microstructure to HCF performance. Together, these two links formed the DTMM that calculates HCF behavior from AM process variables. Two injector prototypes were additively manufactured. The first injector design strategy aggressively pursued low-loss performance by substantially increasing the oxidizer flow area. The combination of manufacturing lead times and the fatigue testing schedule meant that the DTMM was not available when building this first prototype. Therefore, its process parameters were chosen based on a manual review of the available coupon data. This prototype was built successfully and evaluated in 58 combustion tests. Sustained detonation was achieved with remarkably reduced pressure loss, and some tests even displayed pressure loss characteristics similar to conventional gas turbine combustors. This achieved the project goal of improving RDE injector design. The second injector was manufactured according to the optimized parameters predicted by the DTMM. The flow area modifications of this injector were less aggressive than the first injector since demonstrating low pressure loss was not an objective of the second hot-fire test series. Rather, the test objective was to cause high cycle fatigue failure in the part due to periodic loading from the rotating detonation wave. The observed number of cycles to failure was to be compared to the number predicted by the DTMM and thereby assess the utility of the DTMM in component design. However, the required level of vibration was not obtained during combustion. Therefore, high cycle fatigue was not experienced in the hot-fire tests of the second injector. Fatigue data was obtained by further testing the second injector in a conventional HCF test apparatus. The injector demonstrated HCF strength above the DTMM prediction. In fact, it did not fail and testing was only discontinued due to reaching the end of the period of performance. This points to some success in the project’s primary goal of successfully developing and applying the DTMM to a component design. Implementing the DTMM recommendations for optimal processing parameters led to a part with acceptable properties. The DTMM was also shown to be an efficient correlator of data and to provide insight into the relationship between process settings, microstructure, and property performance. However, the failure of the DTMM prediction to match the experimental result of the injector fatigue test also points to the need to include significantly more data in the model development. In this project, coupons made with identical processing parameters exhibited drastically different properties from each other and from the injector part, which clearly influences the accuracy of a model that predicts performance based on parameters. Uncertainties in the build process must be quantified to develop more robust models. A denser and broader matrix of coupon process and geometry variations, several repeated builds of every point, more in-situ build process measurements, and direct observation of tensile and HCF sample microstructure (as opposed to separate microstructure specimens) are recommendations to improve future AM modeling efforts.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2026

Temperature‐Dependent Crystallization in Two‐Step Perovskite Deposition Revealed by In Situ GIWAXS and Machine Learning‐Guided Analysis

The performance and stability of perovskite solar cells are strongly governed by the crystallization behavior of their active layer. In two-step sequential deposition, early-stage film formation plays a decisive role in determining final phase purity and device quality. Guided by a data-driven analysis of nearly 39 000 devices in the FAIR perovskite database, we identified solvent-mediated quenching and thermal processing as key variables affecting power conversion efficiency (PCE), particularly in two-step fabrication. Here, to investigate these effects in real time, we designed and implemented a custom-built, temperature-controlled spin-coating system, enabling precise thermal modulation during precursor deposition. Using this platform, we performed in situ GIWAXS measurements to study the crystallization dynamics of FA 0.5 MA 0.5 PbI 3 films over a temperature range of 30°C–90°C. Our results reveal a non-monotonic relationship between spin-coating temperature and α-phase formation, governed by the interplay between precursor interdiffusion, PbI 2 crystallinity, and δ-phase suppression. The custom thermal control enabled us to isolate and quantify these competing effects during the earliest stages of film formation, providing mechanistic insight into how spin-coating temperature governs both phase purity and kinetic pathways in two-step perovskite systems. Temperature-dependent SEM and photovoltaic device measurements further demonstrate that early-stage crystallization pathways directly translate into differences in morphology, charge-transport continuity, and device performance. These findings inform targeted strategies for optimizing deposition protocols to balance rapid nucleation, phase stability, and device performance.

Saadawy, Ahmed [King Fahd University of Petroleum ↗
DOE OSTI2026

Flavor, transverse momentum, and azimuthal dependence of charged pion multiplicities in semi-inclusive deep-inelastic scattering with 10.6 GeV electrons

Measurements of semi-inclusive deep-inelastic scattering multiplicities for 𝜋 + and 𝜋 − from proton and deuteron targets are reported on a grid of hadron kinematic variables 𝑧, 𝑃 𝑇 , and 𝜙* for leptonic kinematic variables in the range 0.3<𝑥<0.6 and 3<𝑄 2 <5GeV 2 . Data were acquired in 2018 and 2019 at Jefferson Lab Hall C with a 10.6 GeV electron beam impinging on 10-cm-long liquid hydrogen and deuterium targets. Scattered electrons and charged pions were detected in the High Momentum Spectrometer and Super High Momentum Spectrometer, respectively. The multiplicities were fitted for each bin in (𝑥,𝑄 2 ,𝑧,𝑃𝑡) to extract the 𝜙*—independent 𝑀 0 and the azimuthal modulations ⟨cos⁡(𝜙*)⟩ and ⟨cos⁡(2⁢𝜙*)⟩. The 𝑃𝑡 dependence of the 𝑀 0 results was found to be remarkably consistent for the four cases studied: 𝑒⁢𝑝→𝑒⁢𝜋 + ⁢𝑋, 𝑒⁢𝑝→𝑒⁢𝜋 − ⁢𝑋, 𝑒⁢𝑑→𝑒⁢𝜋 + ⁢𝑋, 𝑒⁢𝑑→𝑒⁢𝜋 − ⁢𝑋 over the range 0GeV<𝑃 𝑡 <0.4GeV, as were the multiplicities evaluated near 𝜙*=180∘ over the extended range 0GeV<𝑃 𝑡 <0.7GeV. The Gaussian widths of the 𝑃 𝑡 dependence exhibit a quadratic increase with 𝑧. The cos⁡(𝜙*) modulations were found to be consistent with zero for 𝜋 + , in agreement with previous world data, while the 𝜋 − moments were, in many cases, significantly greater than zero. The cos⁡(2⁢𝜙*) modulations were found to be consistent with zero. The higher statistical precision of this dataset of about 20 000 individual multiplicity values, compared with previously published data, should allow improved determinations of quark transverse momentum distributions and higher twist contributions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Cost of Fish Exclusion and Passage Technologies for Hydropower

Hydropower represents a reliable source of renewable energy and accounts for approximately 7% of the total electrical generation in the United States. Future expansion of hydropower is likely to be in the form of either smaller new stream development projects or powering existing non-powered dams. For these new projects to be successful, careful analysis of risks, costs, and uncertainty to offset reduced power production as well as ensuring the protection and safe passage of migratory fish to gain public support, will be required. Exclusion and passage are two common approaches to protect fish from entrainment and impingement at hydropower facilities. The thresholds for entrainment risk and requirements for exclusion and passage often differ depending on the species involved, the characteristics of the facility, and the goals of stakeholders. While the costs associated with environmental mitigations represent a large proportion of the total costs required for the licensing of hydropower facilities, little quantitative information is present within the literature regarding the specific costs of fish exclusion and passage. Working with FOA awardee Natel Energy, scientists at Oak Ridge National Laboratory were tasked with assessing the capital construction costs for downstream fish exclusion and passage infrastructure. This report used keyword searches of an existing environmental mitigation cost data set and manual extraction of additional cost data associated with protection, mitigation, and enhancement (PM&E) measures related to positive barrier screening and passage from regulatory licensing documents available in the Federal Energy Regulatory Commission (FERC) eLibrary. This approach yielded a total of 50 PM&E mitigation measures with estimated capital construction costs pertaining to positive barrier screens, 142 pertaining to passage studies, and 26 pertaining to passage-related studies. PM&E measures associated with positive barrier screens represented <10% of the 171 total FERC project dockets available in the data set. These data were highly skewed toward conventional relicensing projects, as <7% were associated with new stream development (NSD) projects. Results from these data indicate highly variable costs associated with fish screening, with flow-normalized costs one to two orders of magnitude higher for screening with the highest exclusion capability (≤0.09 in. spacing) compared with coarser screening (1 to 2 in.). Furthermore, estimated capital costs of passage infrastructure were positively related to the scale of the project based on installed capacity for some, but not all, types of passage. These data provide an initial baseline for estimating exclusion and passage costs for hydropower development and may help developers consider options for more fish-friendly generation technologies, though gaps remain relating to a lack of data, particularly for NSD projects. More data may still be available within the FERC eLibrary, but significant effort will be required to manually identify and extract the data for future analyses.

13 HYDRO ENERGY↗
DOE OSTI2025

Simulation insights into wetting properties of hydrogen-brine-clay for hydrogen geo-storage

Hydrogen geo-storage is attracting substantial interdisciplinary interest as a cost-effective and sustainable option for medium- and long-term storage. Hydrogen can be stored underground in diverse formations, including aquifers, salt caverns, and depleted oil and gas reservoirs. The wetting dynamics of the hydrogen-brine-rock system are critical for assessing both structural and residual storage capacities, and ensuring containment safety. Through molecular dynamics simulations, we explore how varying concentrations of cushion gases (CO 2 or CH 4 ) influence the wetting properties of hydrogen-brine-clay systems under geological conditions (15 MPa and 333 K). We employed models of talc and the hydroxylated basal face of kaolinite (kaoOH) as clay substrates. Our findings reveal that the effect of cushion gases on hydrogen-brine-clay wettability is strongly dependent on the clay-brine interactions. Notably, CO 2 and CH 4 reduce the water wettability of talc in hydrogen-brine-talc systems, while exerting no influence on the wettability of hydrogen-brine-kaoOH systems. Detailed analysis of free energy of cavity formation near clay surfaces, clay-brine interfacial tensions, and the Willard-Chandler surface for gas-brine interfaces elucidate the molecular mechanisms underlying wettability changes. Our simulations identify empirical correlations between wetting properties and the average free energy required to perturb a flat interface when clay-brine interactions are less dominant. Here, our thorough thermodynamic analysis of rock-fluid and fluid-fluid interactions, aligning with key experimental observations, underscores the utility of simulated interfacial properties in refining contact angle measurements and predicting experimentally relevant properties. These insights significantly enhance the assessment of gas geo-storage potential. Prospectively, the approaches and findings obtained from this study could form a basis for more advanced multiscale simulations that consider a range of geological and operational variables, potentially guiding the development and improvement of geo-storage systems in general, with a particular focus on hydrogen storage.

25 ENERGY STORAGE↗