Search NASA⌕ Search

SEARCH · Search NASA

Results for “COUPLING NETWORK”

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 1,423 records · Page 79

A multirate mass transfer model to represent the interaction of multicomponent biogeochemical processes between surface water and hyporheic zones (SWAT-MRMT-R 1.0)

Abstract. Surface water quality along river corridors can be modulated by hyporheic zones (HZs) that are ubiquitous and biogeochemically active. Watershed management practices often ignore the potentially important role of HZs as a natural reactor. To investigate the effect of hydrological exchange and biogeochemical processes on the fate of nutrients in surface water and HZs, a novel model, SWAT-MRMT-R, was developed coupling the Soil and Water Assessment Tool (SWAT) watershed model and the reaction module from a flow and reactive transport code (PFLOTRAN). SWAT-MRMT-R simulates concurrent nonlinear multicomponent biogeochemical reactions in both the channel water and its surrounding HZs, connecting the channel water and HZs through hyporheic exchanges using multirate mass transfer (MRMT) representation. Within the model, HZs are conceptualized as transient storage zones with distinguished exchange rates and residence times. The biogeochemical processes within HZs are different from those in the channel water. Hyporheic exchanges are modeled as multiple first-order mass transfers between the channel water and HZs. As a numerical example, SWAT-MRMT-R is applied to the Hanford Reach of the Columbia River, a large river in the United States, focusing on nitrate dynamics in the channel water. Major nitrate contaminants entering the Hanford Reach include those from the legacy waste, irrigation return flows (irrigation water that is not consumed by crops and runs off as point sources to the stream), and groundwater seepage resulting from irrigated agriculture. A two-step reaction sequence for denitrification and an aerobic respiration reaction is assumed to represent the biogeochemical transformations taking place within the HZs. The spatially variable hyporheic exchange rates and residence times in this example are estimated with the basin-scale Networks with EXchange and Subsurface Storage (NEXSS) model. Our simulation results show that (1), given a residence time distribution, how the exchange fluxes to HZs are approximated when using MRMT can significantly change the amount of nitrate consumption in HZs through denitrification and (2) source locations of nitrate have a different impact on surface water quality due to the spatially variable hyporheic exchanges.

58 GEOSCIENCES↗

Adaptive Linear State Estimation for Unbalanced Distribution System

The inclusion of PMU functionality in distribution relays enables the implementation of a linear state estimator (LSE) in Distribution Systems (DS). However, the unbalanced topology and phase coupling in distribution lines necessitate modifications to the LSE formulation. Additionally, the higher fault frequency in distribution systems requires a state estimation approach that is resilient to contingencies. This work proposes an adaptive linear state estimation algorithm tailored for unbalanced distribution systems with single-phase and two-phase laterals. Furthermore, a modified Optimal PMU Placement (OPP) strategy is introduced to ensure full observability in distribution systems with single-phase and two-phase buses. To maintain adaptability to topology changes, the state estimator incorporates circuit breaker status data provided by PMUs, ensuring robust performance during topology changes triggered by faults. The performance of the algorithm is verified on the IEEE 13-bus, 34-bus, and 123-bus systems.

PMUs↗

Identification of Dynamic Force Coefficients for an Additively Manufactured Hermetic Squeeze Film Bearing Support Damper Utilizing a Pass-Through Channel

Abstract The following paper presents breakthrough experimental results for a new hermetic squeeze film damper (HSFD) concept that is integrally designed within an externally pressurized tilting-pad radial gas bearing support. The flexibly damped gas bearing module was designed for a 7.2″ (183 mm) diameter shaft and fabricated using direct metal laser melting (DMLM); also known as additive manufacturing. The bearing and HSFD were sized based on ongoing studies for oil-free supercritical carbon dioxide (sCO2) power turbines in the 8.5 MW–10 MW power range. The development of the new damper concept was motivated by past dynamic testing on HSFD, which generated frequency-dependent stiffness and damping force coefficients. In efforts to eliminate the frequency dependency, a new HSFD architecture was conceived that adds accumulator volumes and a pass-through channel to previously conceived HSFD flow network designs. The other motivation for the work is the need to develop a cost-effective and reliable oil-free bearing technology that is scalable to large power turbomachinery applications. There were several objectives for the following work. The first objective was to successfully design and fabricate a single piece bearing-damper using additive manufacturing, while dimensionally controlling critical design features. The paper discusses the manufacturing steps and shows cut-ups that reveal adequate clearance control capability with internal damper clearances. The second objective was to perform experimental testing with the new HSFD design in efforts to extract stiffness and damping coefficients for excitation frequencies within 20–160 Hz and peak vibration amplitudes between 0.25 mils (6.35 microns) to 1 mil (25.4 microns). The test results for a single HSFD bearing module indicated that the design modifications to the HSFD architecture were successful in eliminating nearly all the frequency dependencies for the stiffness force coefficient. The dynamic tests yielded a stiffness coefficient that varied between 112 klb/in. (19.6 MN/m) and 96 klb/in. (16.8 MN/m). The damping force coefficient however, exhibited relatively more variation with frequency with values residing between 175 lb-s/in. (31 kN-s/m) to 214 lb-s/in. (37 kN-s/m). Finally, the paper advances a three-dimensional fluid-structure interaction (FSI) model using transient finite element analysis (FEA) coupled to a computational fluid dynamics (CFD) model. The FSI analysis performed between 20 Hz and 80 Hz was used to predict the stiffness and damping of the HSFD using a quarter-section model of the damper. The FSI analysis was able to support test results by showing only a 6–7.4% change in the magnitude of force coefficients. Stiffness predictions agree reasonably well with experiments whereas damping is underpredicted.

Engineering↗

The Identity and Chemistry of C 7 H 7 Radicals Observed during Soot Formation

Here we used aerosol mass spectrometry coupled with tunable synchrotron photoionization to measure radical and closed-shell species associated with particle formation in premixed flames and during pyrolysis of butane, ethylene, and methane. We analyzed photoionization (PI) spectra for the C 7 H 7 radical to identify the isomers present during particle formation. For the combustion and pyrolysis of all three fuels, the PI spectra can be fit reasonably well with contributions from four radical isomers: benzyl, tropyl, vinylcyclopentadienyl, and o-tolyl. Although there are significant experimental uncertainties in the isomeric speciation of C7H 7 , the results clearly demonstrate that the isomeric composition of C 7 H 7 strongly depends on the combustion or pyrolysis conditions and the fuel or precursors. Fits to the PI spectra using reference curves for these isomers suggest that all of these isomers may contribute to m/z 91 in butane and methane flames, but only benzyl and vinylcyclopentadienyl contribute to the C 7 H 7 isomer signal in the ethylene flame. Only tropyl and benzyl appear to play a role during pyrolytic particle formation from ethylene, and only tropyl, vinylcyclopentadienyl, and o-tolyl appear to participate during particle formation from butane pyrolysis. There also seems to be a contribution from an isomer with an ionization energy below 7.5 eV for the flames but not for the pyrolysis conditions. Kinetic models with updated and new reactions and rate coefficients for the C 7 H 7 reaction network predict benzyl, tropyl, vinylcyclopentadienyl, and o-tolyl to be the primary C 7 H 7 isomers and predict negligible contributions from other C 7 H 7 isomers. These updated models provide better agreement with the measurements than the original versions of the models but, nonetheless, underpredict the relative concentrations of tropyl, vinylcyclopentadienyl, and o-tolyl in both flames and pyrolysis and overpredict benzyl in pyrolysis. Our results suggest that there are additional important formation pathways for the vinylcyclopentadienyl, tropyl, and o-tolyl radicals and/or loss pathways for the benzyl radical that are currently unaccounted for in the present models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Inversion Reduction Method for Real and Complex Distribution Feeder Models

We report the proliferation of distributed generation on distribution feeders triggers a large number of integration and planning studies. Further, the complexity of distribution feeder models, short simulation time steps, and long simulation horizons rapidly render studies computational burdensome. To mend this issue, we propose a methodology for reducing the number of nodes, loads, generators, line, and transformers of p-phase distribution feeders with unbalanced loads and generation, non-symmetric wire impedance, mutual coupling, shunt capacitance, and changes in voltage and phase. The methodology is derived on a constant power load assumption and employs a Gaussian elimination inversion technique to design the reduced feeder. Compared to previous work by the authors, the inversion reduction takes half the time and voltage errors after reduction are reduced by an order of magnitude. Using a snapshot simulation the reduction is tested on six additional publicly available feeders with a maximum voltage error 0.0075 p.u. regardless of feeder size or complexity, and typical errors on the order of 1 × 10 -4 p.u. Lastly, for a day long quasi-static time series simulation on the UCSD A feeder, errors are shown to increase with changes in loading when a large number of buses removed, but shows less variation for less than 85% of buses removed.

24 POWER TRANSMISSION AND DISTRIBUTION↗

SIW21-95: Hybridizing Synchronous Condensers with Grid-Forming Battery Energy Storage Systems

One of the main challenges associated with deployment of high shares of inverter-based resources (IBRs) in power grid is not only reduced system inertia but also degrading system strength that may cause severe stability impacts. A minimum level of system strength is needed for the power system to remain stable under normal conditions and to return to a steady state condition following a system disturbance. Significant system strength reduction is expected in almost all planned areas for solar and wind generation deployment. Synchronous condensers (SC) have been considered as one main technology to address the system strength issues for the areas with high levels of IBRs. SCs can help improving reliability and resiliency of power system but they do not provide the full range of services needed by power systems for reliable and economic integration of high shares of inverter-coupled variable generation such as PV generation. Services related to active power controls cannot be provided by SCs due to lack of prime mover. Even for provision of reactive power SCs have certain constraints based on their thermal and stability limits. NREL has been conducting research on a hybridized concept that combines SCs with grid forming (GFM) battery energy storage systems (BESS). This super flexible AC transmission system (SuperFACTS) that combines these two technologies in a single plant under the same controller offers a unique scalable set of services to the power system at all levels (transmission, sub-transmission, distribution, islands and isolated microgrids). Depending on use cases, SuperFACTS can be controlled to provide fully dispatchable and flexible operation using energy storage component, provide a full range of existing and future ancillary and reliability services to the grid (similar or better than conventional sources), maintain adequate levels of grid strength and inertia, and provide fault current for proper operation of protection systems. Therefore, this economic and easy to commercialize solution has potential for significant impacts on certain segments of global energy sector. All types of gird services (market based, reliability and resiliency) can be provided by SuperFACTS plants. GFM BESS can act as a self-black start source for each SuperFACTS module, which in turn can act a black start resource for co-located PV and wind power plants, segments of transmission and distribution networks, for conventional power plants, etc. The issue of in-rush currents during black start is addressed by overcurrent capability of SCs. In this paper we describe the results of modeling for SuperFACTS concept.

grid forming↗

Decoding Zeolite Crystallization and Stage III in Nuclear Waste Glasses by Coupled Modeling and Experiments

Under specific conditions of pH and temperature, nuclear waste immobilization borosilicate glasses may exhibit a sudden acceleration in their corrosion kinetics (stage III)—a behavior that has been associated with the formation of zeolite crystals. Such accelerated dissolution may compromise the integrity of nuclear wasteforms placed in geological depositories. However, thus far, none of the available models is able to predict the thermodynamic propensity and kinetics of zeolite precipitation as a function of the solution conditions due to (i) a lack of fundamental knowledge regarding the nucleation & growth mechanisms of zeolitic phases, (ii) uncertainty regarding the compositions (types) of zeolites that may form and the rate-limiting step in their precipitation as a function of the solution conditions, and (iii) the complexities that arise due to the vast parametric space (i.e., solution chemistry, temperature, number of secondary phases, etc.) that encompass these systems under conditions of environmental exposure. To resolve these challenges, this project aimed to unambiguously identify the thermodynamic propensity for zeolite precipitation and the kinetics thereof as a function of the solution conditions (composition, pH, and temperature). To achieve this goal: 1) We identified the solution conditions and zeolite phases relevant to nuclear glass dissolution. 2) We performed a series of ab initio molecular dynamics (AIMD) simulations to compute the thermodynamic properties of a group of characteristic zeolites that features a large range of compositions, various hydration levels, a wide range of framework structures, and partial atomic site occupancies. 3) We released a first-of-a-kind self-consistent thermodynamic database that can be used to assess the kinetics and the stability fields of zeolitic phases within a Gibbs energy minimization (GEM) framework. 4) We developed a robust geochemical modeling method allowing us to predict the stability of secondary phases (including zeolites, calcium–silicate–hydrate gels, and clays) upon the dissolution of nuclear waste immobilization glasses. 5) We introduced a model that predicts the dissolution kinetics of a series of borosilicate nuclear waste immobilization glasses in terms of the topology of their atomic network. 6) We investigated the roles of the solution composition on the crystallization kinetics of phillipsite zeolites and tobermorite silicate hydrates. Via PNNL’s collaboration and engagement, this project directly supports DOE’s nuclear waste immobilization activities by offering a technical, science-based foundation that will (i) facilitate predictions of the long-term corrosion rates and extents of existing nuclear waste immobilization glasses to help ensure safe and successful vitrification operations, (ii) inform the development of advanced glass formulations with enhanced durability, and, (iii) enable cost-savings that result from making more decisive and hence less conservative predictions while offering higher levels of nuclear waste embedment in smaller, more compact glass volumes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Nano Graphene-Reinforced Bio-nanocomposites Based on NR/PLA: The Morphological, Thermal and Rheological Perspective

Bio-based graphene-reinforced thermoplastic elastomers (TPE) based on natural rubber (NR), and poly lactic acid (PLA) were successfully prepared via melt blending. The effect of graphene nanosheets (GNS) content and coupling agent were investigated on morphological, thermal and rheological properties of the PLA/NR/GNS bio-nanocomposites. More stable morphology of uniformly dispersed NR phase in the continuous PLA matrix with a relatively narrower diameter-size distribution was achieved in the presence of a compatibilizer. Moreover, transmission electron micrographs showed more improved dispersion of GNS with predominantly exfoliated morphology in the case of ENR-compatibilized PLA/NR/GNS bio-nanocomposites. Crystallization/melting studies revealed that, as the GNS content increased, the cold-crystallization peak, melting peak and glass-to-rubber transition temperature values shifted to higher temperatures and the crystallinities of blends slightly decreased (about 7%) which might be attributed to the reduction of polymeric segments mobility by the restricting effect of GNS. Dynamical mechanical investigation showed that the storage modulus (E′) increased about 25% by the introduction of GNS due to the inherent stiffness of GNS and the occurrence of compatibility in the PLA/NR blend in the presence of ENR. The presence of ENR shifted the T{sub g} of the NR phase and the PLA matrix towards each other which is a characteristic of higher miscibility and compatibility. The thermogravimetry (TGA) and derivative thermogravimetry (DTG) curves revealed higher thermal stability of the compatibilized-PLA/NR blends due to the enhanced interfacial adhesion and the homogenous dispersion of GNSs as direct effects of the presence of ENR. Rheological studies indicated that the formation of the effective GNS-polymer networks by the presence of ENR increased the storage modulus (G′) and complex viscosity (η*) value due to the effectiveness of the GNS in taking loads and restricted molecular motion, respectively. Graphical Abstract: .

36 MATERIALS SCIENCE↗

Ab Initio Computations and Active Thermochemical Tables Hand in Hand: Heats of Formation of Core Combustion Species

Here, the fidelity of combustion simulations is strongly dependent on the accuracy of the underlying thermochemical properties for the core combustion species that arise as intermediates and products in the chemical conversion of most fuels. High level theoretical evaluations are coupled with a wide-ranging implementation of the Active Thermochemical Tables (ATcT) approach to obtain well-validated high fidelity predictions for the 0 K heat of formation for a large set of core combustion species. In particular, high level ab initio electronic structure based predictions are obtained for a set of 348 C, N, O, and H containing species, which corresponds to essentially all core combustion species with 34 or fewer electrons. The theoretical analyses incorporate various high level corrections to base CCSD(T)/cc-pVnZ analyses (n = T or Q) using H 2 , CH 4 , H 2 O, and NH 3 as references. Corrections for the complete-basis-set limit, higher-order excitations, anharmonic zeropoint energy, core–valence, relativistic, and diagonal Born–Oppenheimer effects are ordered in decreasing importance. Independent ATcT values are presented for a subset of 150 species. The accuracy of the theoretical predictions is explored through (i) examination of the magnitude of the various corrections, (ii) comparisons with other high level calculations, and (iii) through comparison with the ATcT values. The estimated 2σ uncertainties of the three methods devised here, ANL0, ANL0-F12, and ANL1, are in the range of ±1.0–1.5 kJ/mol for single-reference and moderately multireference species, for which the calculated higher order excitations are 5 kJ/mol or less. In addition to providing valuable references for combustion simulations, the subsequent inclusion of the current theoretical results into the ATcT thermochemical network is expected to significantly improve the thermochemical knowledge base for less-well studied species.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Searches for New Physics With Muon Conversion at Fermilab and Triboson Production at the LHC

We report on several efforts to search for physics beyond the standard model of particle physics at broad energy scales. The Mu2e experiment at Fermilab will search for charged lepton flavor violation via the muon to electron conversion process, which is suppressed in the Standard Model. Mu2e will be operated at a low energy, yet can probe New Physics at very high mass scales (O(1e3 - 1e4 ) TeV). At high energies, the CMS experiment at the CERN LHC continues to deliver an impressive suite of Standard Model measurements and limits on a variety of New Physics signatures. Mu2e is under construction and slated to collect its first physics data in the coming years. This thesis describes work done during the construction phase of Mu2e and focuses on two critical areas: magnetic field modeling and statistical analysis. We describe a novel method for field modeling which we validate using a simulated dataset representing the expected magnetic field in the Detector Solenoid. This method blends a standard least-squares fitting technique that utilizes physically motivated analytical model functions with a novel physics informed network that is constructed to obey Maxwell’s equations. We show the technique can model the field with an accuracy of 10−7 despite the presence of injected noise in the pseudo-measurements at the 10−5 level. We then present preliminary results of the calibration of 3D Hall probes at the sub-10−4 level. These probes will be used to directly measure the Mu2e Detector Solenoid magnetic field on a sparse grid; these measurements serve as the input to the field model fitting. Finally, we describe the first implementation of both an unbinned shape analysis and a Bayesian interpretation applied to Mu2e pseudo-data. Up to 20% tighter limits can be set by the shape analysis compared to a standard cut & count analysis. The AlCap experiment collected data at PSI in 2015 to measure several important quantities related to nuclear muon capture on an aluminum target, which is a significant background process for Mu2e. The neutron emission from muon capture can introduce background hits in the Mu2e detectors and can increase radiation damage in various elements of the apparatus. We present measurements of the neutron group fluence and mean neutron multiplicity for muon capture on aluminum nuclei. Finally, we discuss an analysis of triboson production at CMS using an Effective Field Theory framework. Standard Model triboson production, which was first observed at CMS in 2020, has a relatively small cross section and provides direct access to both anomalous triple gauge couplings and quartic gauge couplings. These couplings, interpreted in the Standard Model Effective Field Theory, are studied in the present work. We target the boosted regime where the background rate is low and yields are enhanced when dimension-6 and dimension-8 Wilson coefficients are non-zero. We do not observe an excess in the data and therefore set bounds on the Wilson coefficients. For dimension-6 coefficients the tightest observed (expected) bounds are set on cW /Λ2 where Λ is the mass scale of new physics; the bounds are [−0.13, 0.12] TeV−2 ([−0.12, 0.12] TeV−2 ) at 95% CL. The tightest bounds in dimension-8 are set on fT,0 / Λ4 ; the observed (expected) bounds at 95% CL are [−0.63, 0.69] TeV−4 ([−0.54, 0.62] TeV−4 ). Additional results are presented which include scenarios where multiple Wilson coefficients are non-zero, the application of signal model clipping to address unitarity violation in Effective Field Theories, and a novel template fit developed for easier reinterpretation of our results.

Kampa, Cole Erik [Northwestern U. (main)] (ORCID:0↗

Modeling Functional Organic Chemistry in Arctic Rivers: An Idealized Siberian System

Rivers of the Arctic will become ever more important for the global climate, since they carry a majority of continental dissolved organic carbon flux into the rapidly changing polar ocean. Aqueous organics comprise a wide array of functional groups, several of which are likely to impact coastal and open water biophysical properties. Light attenuation, interfacial films, aerosol formation, gas release and momentum exchange can all be cited. We performed Lagrangian kinetic modeling for the evolution of riverine organic chemistry as the molecules in question make their way from the highlands to Arctic outlets. Classes as diverse as the proteins, sugars, lipids, re-condensates, humics, bio-tracers and small volatiles are all included. Our reduced framework constitutes an idealized northward flow driving a major hydrological discharge rate and primarily representing the Russian Lena. Mountainous, high solute and tundra sources are all simulated, and they meet up at several points between soil and delta process reactors. Turnover rates are parameterized beginning with extrapolated coastal values imposed along a limited tributary network, with connections between different terrestrial sub-ecologies. Temporal variation of our total dissolved matter most closely resembles the observations when we focus on the restricted removal and low initial carbon loads, suggesting relatively slow transformation along the water course. Thus, channel combinations and mixing must play a dominant role. Nevertheless, microbial and photochemical losses help determine the final concentrations for most species. Chemical evolution is distinct for the various functionalities, with special contributions from pre- and post-reactivity in soil and delta waters. Several functions are combined linearly to represent the collective chromophoric dissolved matter, characterized here by its absorption. Tributaries carry the signature of lignin phenols to segregate tundra versus taiga sources, and special attention is paid to the early then marine behaviors of low molecular weight volatiles. Heteropolycondensates comprise the largest percentage of reactive carbon in our simulations due to recombination/accumulation, and they tend to be preeminent at the mouth. Outlet concentrations of individual structures such as amino acids and absorbers lie above threshold values for biophysical influence, on the monolayer and light attenuation. The extent of coastal spreading is examined through targeted regional box modeling, relying on salinity and color for calibration. In some cases, plumes reach the scale of peripheral arctic seas, and amplification is expected during upcoming decades. Conclusions are mapped from the Lena to other boreal discharges, and future research questions are outlined regarding the bonding type versus mass release as permafrost degrades. Dynamic aqueous organic coupling is recommended for polar system models, from headwaters to coastal diluent.

54 ENVIRONMENTAL SCIENCES↗

Automated Clean Chemistry for Bulk Analysis of Environmental Swipe Samples – Final Report

Sample preparation methods for mass spectrometry are being automated using commercial-off-the-shelf (COTS) equipment to shorten lengthy and costly manual chemical purification procedures. This development addresses a serious need in the International Atomic Energy Agency’s (IAEA) Network of Analytical Laboratories (NWAL) to increase efficiency in the Bulk Analysis of Environmental Samples for Safeguards program with a method that allows unattended, overnight operation. In collaboration with Elemental Scientific Inc. (ESI), the prepFAST-MC2 was designed based on COTS equipment. It was modified for uranium/plutonium separations using renewable columns packed with Eichrom TEVA® and UTEVA® resins, with a chemical separation method based on the Oak Ridge National Laboratory (ORNL) NWAL chemical procedure. Initial verification of individual columns yielded small elution volumes with consistent elution profiles and good recovery. Combined column calibration demonstrated ample separation without crosscontamination of the eluent. Automated packing and unpacking of the built-in columns initially showed >15% deviation in resin loading by weight, which can lead to inconsistent separations. Optimization of the packing and unpacking methods led to a reduction in the variability of the packed resin to less than 5% daily. Based on the initial experience at ORNL with the prepFAST-MC2, a second-generation model, the prepFAST-SR, was built by ESI. The newly designed prepFAST-SR has several upgrades compared with the original prepFAST-MC2. The upgrades largely addressed issues with the reliability of unpacking the Eichrom columns and with ventilation to ensure adequate removal of acid vapors. Both systems are currently installed in the Ultra-Trace Forensics Science Center at ORNL. Separations of mixed uranium and plutonium samples containing certified reference materials were analyzed by multi-collector inductively coupled plasma mass spectrometry and yielded good results. Current efforts have demonstrated a wider applicability of the prepFAST system. Near-quantitative removal of metal interferences was achieved with the system in both uranium and plutonium separated fractions. Successful system validation was completed with several archived samples. Isotopic results from archived samples and certified reference materials were well within data quality limits for the IAEA NWAL. Additional COTS equipment has been evaluated for its potential to aid the prepFAST-SR system in reducing the time allotments and cleanroom infrastructure requirements for accurate separations. Overall, these efforts help ensure continued efficient and effective operation of the NWAL. Comparison of the amount of personnel time necessary for successful manual vs. automated chemical separations showed a significant decrease in hands-on time from 9.8 hours to 35 minutes for seven samples, respectively. This documented time savings and reduced labor translates to a significant cost savings per sample. Overall, the system will enable faster sample reporting times at reduced cost by limiting personnel hours dedicated to the chemical separation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Confronting Models with Regional CO 2 Observations

This research addressed an important DOE science area: the role of natural disturbances in carbon cycling. Our project included both (a) modeling and field measurements and (b) strengthened ties between DOE national laboratories and universities which continue after the project was completed. Our focus was on understanding the consequences of drought and emissions (fire and urban) to atmospheric trace gas composition (carbon dioxide, carbon monoxide, and methane), including both concentration and stable isotope composition. Measurements of these atmospheric gases were used to infer the spatial and temporal patterns of both sources and sinks in the carbon cycle. As part of this effort, we maintained two long-term monitoring networks along a geographic gradient in Utah and Colorado that spans montane forests, urban regions, and oil/gas fields. In this research, we expanded on a well-accepted atmospheric model (STILT-WRF) as a tool to evaluate fluxes in CLM/CESM-related models. We used the STILT model to provide a strong linkage between point atmospheric measurements and the surface parameterizations/emissions that are part of the CLM/CESM models. Overall, this project benefited and supported the DOE Long Term Mission and Goals in four distinct ways: (a) testing of carbon cycle models, (b) testing and evaluating model mechanisms whereby factors in CLM influence trace gas composition in CESM, (c) acquiring high-quality, long-term data on concentrations of carbon dioxide in western USA ecosystems, and (d) by reducing uncertainties associated with the representation of climate-carbon feedbacks in Earth System models through the development of new methods for evaluating model performance. One postdoctoral associate (Henrique Duarte) and one Ph.D. student (Darek Malia) were trained with this support. Nine peer-reviewed publications and 14 poster presentations were produced with support from this project. This proposal addresses a specific science area in the RFP: the role of natural disturbances in carbon cycling. In particular, we will focus on the consequences of drought and fire at landscape to regional scales for atmospheric trace gas concentration and stable isotope composition (CO 2 , CO, and CH 4 ). We seek to expand on the well accepted atmospheric Stochastic Time-Inverted Lagrangian Transport (STILT) model coupled to the Weather Research and Forecasting (WRF) model as a tool to evaluate stable isotope fluxes in CLM/CESM-related models through existing, long-term data sets and field campaigns where high-precision trace gas concentration and isotope measurements can be made along gradients and transects. Given the confounding impacts of anthropogenic trace gas emissions on landscape-to-regional atmospheric observations, we will also measure and parameterize isotopic characteristics of CO 2 , CO, and CH 4 from anthropogenic or forest sources in order to allow us the opportunity to partition atmospheric trace gases into natural and anthropogenic sources.

54 ENVIRONMENTAL SCIENCES↗

Coupling Visualization, Simulation, and Deep Learning for Ensemble Steering of Complex Energy Models: Preprint

We have developed a framework for the exploration, design, and planning of energy systems that combines interactive visualization with machine-learning based approximations of simulations through a general purpose dataflow API. Our system provides a visual inter- face allowing users to explore an ensemble of energy simulations representing a subset of the complex input parameter space, and spawn new simulations to 'fill in' input regions corresponding to new enegery system scenarios. Unfortunately, many energy simula- tions are far too slow to provide interactive responses. To support interactive feedback, we are developing reduced-form models via machine learning techniques, which provide statistically sound esti- mates of the full simulations at a fraction of the computational cost and which are used as proxies for the full-form models. Fast com- putation and an agile dataflow enhance the engagement with energy simulations, and allow researchers to better allocate computational resources to capture informative relationships within the system and provide a low-cost method for validating and quality-checking large-scale modeling efforts.

97 MATHEMATICS AND COMPUTING↗

HydraGNN_Predictive_GFM_2026 - Ensemble of predictive graph foundation models for atomistic materials modeling

This release contains data and parameters of HydraGNN-based graph foundation models trained as a result of the work published in the pre-print "Exascale Multi-Task Graph Foundation Models for Imbalanced, Multi-Fidelity Atomistic Data" by M. Lupo Pasini et al. (https://arxiv.org/abs/2604.15380). We jointly train on 16 open first-principles datasets (544+ million structures covering 85+ elements) using a multi-task architecture with per-dataset heads and a scalable ADIOS2/DDStore data pipeline. On Frontier, we execute six large-scale DeepHyper hyperparameter optimization campaigns in FP64 and promote the top-performing message-passing models to sustained 2,048-node training, yielding a PaiNN-based lead model. The version of HydraGNN used to generate the outputs provided in this release is HydraGNN v5.0 (https://github.com/ORNL/HydraGNN/releases/tag/v5.0) The list of datasets used for the training of the graph foundation model is the following: 1) Alexandria [1] 2) ANI1x [2] 3) MPTrj [3] 4) Open Catalyst 2020 (OC20) [4] 5) Open Catalyst 2022 (OC22) [5] 6) Open Catalyst 2025 (OC25) [6] 7) Open Direct ir Capture 2023 (ODAC23) [7] 8) Open Materials 2024 (OMat24) [8] 9) Open Molecules 2025 (OMol25) [9] 10) OMol25-neutral (subset of OMol25 that contains only molecules with zero total charge) 11) OMol25-non-neutral (subset of OMol25 that contains only molecules with non-zero total charge) 12) Open Polymers 2026 (OPoly2026) [10] 13) Nabla2DFT [11] 14) QCML [12] 15) QM7X [reference 13] 16) transition1x [14] Dataset references: [1] J. Schmidt et al., “A dataset of 175k stable and metastable materials calculated with the PBEsol and SCAN functionals,” Scientific Data, vol. 9, p. 64, 2022. [2] J. S. Smith et al., “The ANI-1ccx and ANI-1x data sets, coupled-cluster and density functional theory properties for molecules,” Scientific Data, vol. 7, p. 134, 2020. [Online]. Available: https: //www.nature.com/articles/s41597-020-0473-z [3] A. Jain et al., “Commentary: The Materials Project: A materials genome approach to accelerating materials innovation,” APL Materials, vol. 1, no. 1, p. 011002, 07 2013. [Online]. Available: https://doi.org/10.1063/1.4812323 [4] L. Chanussot et al., “Open catalyst 2020 (oc20) dataset and community challenges,” ACS Catalysis, vol. 11, no. 10, pp. 6059–6072, 2021. [Online]. Available: https://doi.org/10.1021/acscatal.0c04525 [5] K. Tran et al., “Open catalyst 2022 (oc22) dataset and challenges for oxidation electrocatalysts,” ACS Catalysis, vol. 13, no. 5, pp. 3066–3084, 2023. [Online]. Available: https://doi.org/10.1021/acscatal.2c05426 [6] S. J. Sahoo et al., “The open catalyst 2025 (oc25) dataset and models for solid-liquid interfaces,” arXiv preprint arXiv:2509.17862, 2025. [Online]. Available: https://arxiv.org/abs/2509.17862 [7] A. Sriram et al., “The open DAC 2023 dataset and challenges for sorbent discovery in direct air capture,” ACS Central Science, vol. 10, no. 5, pp. 923–941, 2024. [8] L. Barroso-Luque et al., “Open materials 2024 (omat24) inorganic materials dataset and models,” 2024. [Online]. Available: https://arxiv.org/abs/2410.12771 [9] D. S. Levine et al., “The open molecules 2025 (OMol25) dataset, evaluations, and models,” 2025. [Online]. Available: https://arxiv.org/abs/2505.08762 [10] D. S. Levine et al., The open polymers 2026 (OPoly26) dataset and evaluations,” arXiv preprint arXiv:2512.23117, 2025. [Online]. Available: https://arxiv.org/abs/2512.23117 [11] K. Khrabrov et al., “Nabla2dft: A universal quantum chemistry dataset of drug-like molecules and a benchmark for neural network potentials,” in NeurIPS 2024 Datasets and Benchmarks Track, 2024. [Online]. Available: https://openreview.net/forum?id=ElUrNM9U8c [12] S. Ganscha et al., “The QCML dataset, quantum chemistry reference data from 33.5M DFT and 14.7B semi-empirical calculations,” Scientific Data, vol. 12, p. 406, 2025. [13] J. Hoja et al., “QM7-X, a comprehensive dataset of quantum-mechanical properties spanning the chemical space of small organic molecules,” Scientific Data, vol. 8, p. 43, 2021. [Online]. Available: https://www.nature.com/articles/s41597-021-00812-2 [14] M. Schreiner et al., “Transition1x - a dataset for building generalizable reactive machine learning potentials,” Scientific Data, vol. 9, p. 779, 2022. The folder "datasets_ADIOS2_format" contains the set of pre-processed datasets in Adaptable I/O System (ADIOS) format (https://www.exascaleproject.org/research-project/adios/) that have been used for the development and training of GFMs in this work. The "datasets_ADIOS2_format" directory contains 2 sub-directories, one for the version "v1" of the datasets and one for the version "v2" of the datasets. The version "v1" of the datasets provides values of the total energy as they are extracted from the original data as it was released by the respective institutions. The version "v2" of the datasets provides values of the energy that have been realigned. The realignment was performed by training a linear regression model that predicts the total energy as a function of the chemical composition of the atomistic structure, and then subtract such prediction from the original value of the total energy. Both folders "v1" and "v2" contain 16 sub-directories, each corresponding to an ADIOS2-formatted dataset The folder "DeepHyper-results" contains the configurational files and model's parameters for all the 186 HPO trials that were successfully completed by the scalable hyperparameter optimization (HPO) runs on Frontier. The content of the folder "DeepHyper-results" I structured as follows: 1) task-list.txt: list of mpnn name, jobid, and deephyper task id 2) gfm_${MPNN}_${JOBID}_0.${TASKID}: run directory with checkpoint files 3) gfm_${MPNN}: deephyper summary directory (*.csv) for each specific MPNN type 4) deephyper-experiment-${JOBID}: output and error logs for each job The file "deephyper-sorted.csv" contains the details of each HydraGNN model built and tested by HPO, obtained by merging the (*.csv) filed from each HPO run executed. Out of all the HPO trials, we selected 10 to continue the training of the respective HydraGNN models. Due to limited computational budget available in the LRN070 allocation we could not complete the training till convergence for all these 10 selected models. The folder "models" contains multiple sub-folders, one per each HydraGNN model trained. Each model sub-folder contains the parameters of each HydraGNN model, with multiple checkpoint-restarts. The list of sub-folders are as follows: 1) multidataset_hpo-BEST1-fp64 2) multidataset_hpo-BEST2-fp64 3) multidataset_hpo-BEST3-fp64 4) multidataset_hpo-BEST4-fp64 5) multidataset_hpo-BEST5-fp64 6) multidataset_hpo-BEST6-fp64 7) multidataset_hpo-BEST7-fp64 8) multidataset_hpo-BEST8-fp64 9) multidataset_hpo-BEST9-fp64 10) multidataset_hpo-BEST10-fp64 Within each one of these folders, additional auxiliary log files are provided with descriptions about how the training proceeded. The lead PaiNN-model is contained inside "multidataset_hpo-BEST6-fp64". The file "mlp_branch_weights" contains the parameters of the multi-layer perceptron (MLP) used to reconcile the predictions of the 16 output decoding heads of the HydragNN architectures. The MLP takes in input the chemical composition of the atomistic structure and predicts averaging weights to linearly mix the predictions of each output decoding head toward consolidating them into a single one. The folder "1.1billion-structure-inference" contains 1.1 billion atomistic structures randomly generated. Each structures is associated with energy and forces predicted with the lead-PaiNN model combined with the MLP model for reconciliation of the multi-branch predictions generated by the 16 output decoding heads. The folder "1.1billion-structure-inference" contains 9,300 (*.tar.gz) subdirectories, one per Frontier compute node used to execute the inference at exascale. Once uncompressed, each (*.tar.gz) subdirectory contains an ADIOS2 (*.bp) file container, where each atomistic structure is stored as a PyTorch-Geometric Data object. The file "export_dataset_environment_variables.sh" contains the environment variables that need to be set before running the HydraGNN code to reproduce the results provided in this dataset release. The code that can be used to load the ADIOS2 files, load HydraGNN models, and run inference is available at: https://github.com/ORNL/HydraGNN/releases/tag/v5.0

36 MATERIALS SCIENCE↗

Numerical analysis of soot emissions from gasoline-ethanol and gasoline-butanol blends under gasoline compression ignition conditions

In the present work, computational fluid dynamics (CFD) simulations of a single-cylinder gasoline compression ignition (GCI) engine were performed to investigate the impact of blending two biofuels, ethanol and n-butanol, with gasoline on the trade-off between combustion phasing and soot emissions under low load conditions. Here, in order to represent market gasoline (RD5-87), a four-component toluene primary reference fuel (TPRF) + ethanol (ETPRF) surrogate (with 20% ethanol by mole; E20) was formulated using a neural network based octane predictor such that the surrogate had the same ethanol content, Research Octane Number (RON) and Octane Sensitivity (S). In addition, a novel skeletal kinetic mechanism for ETPRF and TPRF + n-butanol (BTPRF) blends, incorporating polycyclic aromatic hydrocarbon (PAH) chemistry, was developed. A three-dimensional (3D) engine CFD formulation employing the skeletal mechanism, adaptive mesh refinement (AMR), finite-rate chemistry approach, and hybrid method of moments (HMOM) was adopted to capture the in-cylinder combustion phenomena and soot emissions. The engine CFD model was validated against RD5-87 experimental data for a broad range of start-of-injection (SOI) timings (-21/-27/-36/-45 crank angle degrees (CAD) after top-dead center (aTDC)), with respect to in-cylinder pressure, heat release rate, combustion phasing, and soot emissions. The closed-cycle simulation results were analyzed to elucidate the non-monotonic trend of soot emissions versus SOI timing: SOI-36 > SOI-45 > SOI-21 > SOI-27. Thereafter, the validated CFD model was employed to simulate the combustion of a gasoline-ethanol blend with 45% (by mole) ethanol (E45) and a gasoline-butanol blend with 45% (by mole) n-butanol (B45) under the same operating conditions to study the effects of fuel composition and SOI timing on combustion phasing and soot emissions. The sooting propensity followed the trend: B45 > E20 > E45 at all SOI timings. Overall, it was observed that the autoignition propensity was primarily related to fuel chemistry. On the other hand, sooting propensity showed strong coupling with both fuel chemistry and physical properties, with greater impact of fuel physical properties at advanced SOI timings.

30 DIRECT ENERGY CONVERSION↗

Evaluating the Effectiveness of Soil Profile Rehabilitation for Pluvial Flood Mitigation Through Two-Dimensional Hydrodynamic Modeling

Pluvial flooding, driven by increasingly impervious surfaces and intense storm events, presents a growing challenge for urban areas worldwide. In Baltimore City, MD, USA, climate change, rapid urbanization, and aging stormwater infrastructure are exacerbating flooding impacts, resulting in significant socio-economic consequences. This study evaluated the effectiveness of a soil profile rehabilitation scenario using a 2D hydrodynamic modeling approach for the Tiffany Run watershed, Baltimore City. This study utilized different extreme storm events, a high-resolution (1 m) LiDAR Digital Terrain Model (DTM), building footprints, and hydrological soil data. These datasets were integrated into a fully coupled 2D hydrodynamic model, the City Catchment Analysis Tool (CityCAT), to simulate urban flood dynamics. The pre-soil rehabilitation simulation revealed a maximum water depth of 3.00 m in most areas, with hydrologic soil groups C and D, especially downstream of the study area. The post-soil rehabilitation simulation was targeted at vacant lots and public parcels, accounting for 33.20% of the total area of the watershed. This resulted in a reduced water depth of 2.50 m. Additionally, the baseline runoff coefficient of 0.49 decreased to 0.47 following the rehabilitation, and the model consistently recorded a peak runoff reduction rate of 4.10 across varying rainfall intensities. The validation using a contingency matrix demonstrated true-positive rates of 0.75, 0.50, 0.64, and 0 for the selected events, confirming the model’s capability at capturing real-world flood occurrences.

Baltimore City↗

Understanding Adsorption and Reactions at Aqueous Oxide Interfaces with Neural Network Potential Molecular Dynamics

Chemical processes at metal oxide−water interfaces are of central importance in geochemistry, biology, and energy technologies. A better understanding of these processes would allow us to make a significant step toward optimizing and controlling them, which could in turn lead to broader impacts. Computational modeling is indispensable to accomplishing this task because complexity and disorder often make it difficult to extract atomistic information from experiments. Balancing computational cost and accuracy, simulation schemes based on efficient machine learning representations of the potential energy surface (PES) predicted by ab initio calculations have become increasingly popular over the past decade. In particular, several studies have demonstrated the ability of machine learning models to accurately reproduce the complex ab initio PESs of aqueous oxide interfaces, allowing simulations of systems and processes that are not accessible with ab initio methods. In this Account, we review our recent efforts to understand adsorption processes and reactions at aqueous oxide interfaces using deep potential molecular dynamics (DPMD), a simulation scheme employing deep neural networks (DNNs), which has proven to be quite successful in accurately describing many different systems in the condensed phase. After summarizing the DPMD methodology, we first review our work on the acid−base chemistry of oxide surfaces in contact with water, a fundamental characteristic that controls proton transfer and surface charge at the interface. We focus on the aqueous interface of rutile IrO 2 , an oxide material thus far considered the best catalyst for the oxygen evolution reaction (OER). We show that this interface is characterized by a large fraction of dissociated water and a strong Brønsted acidity of the surface sites, in good agreement with the experimentally measured value of the point of zero proton charge. In our second example, we investigate how the adsorption of organic species from ambient air or water affects the structure and wettability of the aqueous interfaces of TiO 2 , a prototypical photocatalytic material. This is a question that is relevant to understanding the UV-induced hydrophilicity of TiO 2 surfaces, a property at the basis of self-cleaning windows and related applications. Specifically focusing on formic and acetic acids, the two most common atmospheric organic acids, our simulations reveal that these acids control the wettability of TiO 2 largely through acid−base chemistry at the interface rather than chemisorption on the oxide surface, a finding that could help improve the design of self-cleaning surfaces and photocatalytic devices. Finally, we review our recent study of methanol at TiO 2 −water interfaces, a system whose interest is largely motivated by the role of methanol in enhancing photocatalytic hydrogen evolution on TiO 2 . Our simulations provide mechanistic insights into the coupled roles of the organic adsorbate and water at the TiO 2 interface, with implications for how methanol enhances the activity of H 2 evolution.

adsorption↗