Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamical balances”

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

Mineralized sclerites in the gorgonian coral Leptogorgia chilensis as a natural jamming system

The soft corals (Cnidaria, Octocorallia), a diverse group of colonial marine invertebrates, can reversibly tune their body stiffness in response to external stimuli. This capability is attributed to their dynamic skeletal systems, which consist of thousands of mineralized skeletal elements, called sclerites, embedded within a gel-like matrix that swells/deswells and unjams/jams the sclerites, thus modulating skeletal stiffness. While sclerite morphology is widely used for species identification, its role in the mechanical performance of a soft coral’s skeletal system is largely unknown. Here, we investigated structure-jamming relationships in sclerite-based skeletal architectures using the red gorgonian octocoral Leptogorgia chilensis as a model system. The sclerites of L. chilensis exhibit a shaft-like geometry with two axial branches and two sets of triradiate side branches, which are aligned with the crystallographic symmetry of the constituent magnesium-containing calcite. By combining multiscale three-dimensional (3D) structural characterization, parametric geometrical modeling, 3D printing, mechanical testing, and discrete element simulations, we demonstrate how sclerite geometry achieves a balanced jamming performance in terms of stiffness, weight, strength, and fracture resistance in comparison to alternative geometries parametrically modified from the native sclerites (e.g., changes in the length and number of side branches). Here, we also found that these performance metrics are achieved through the effective interlocking among side and axial branches, which is further enhanced by the fractal-like microscopic spikes on the branch tips. The findings in this natural jamming system offer insights for designing synthetic mechanotunable material architectures for a wide range of applications, from soft robotics to mechanical dampeners.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Optimisation of the Kaplan hydropower system via PID 2 and digital twin

Here, this paper proposes a proportional–integral-double–derivative (PID 2 ) optimisation method for the Kaplan hydropower system by building a digital twin. The study first uses one multilayer perceptron (MLP) to model the hydroturbine dynamic and then adopts three connected MLPs to model the generator dynamic, both in an open-loop fashion. Inspired by stochastic distribution control (SDC) theory, we regard the training of the turbine's neural network model as a process control problem, and we propose minimising entropy loss to update the network parameters. The next step is to build the digital twin by connecting the neural network models with a PID 2 controller and a lead-lag exciter and run the whole model in a closed-loop fashion. After that, a binary search approach is applied to optimise the PID 2 parameters based on the obtained digital twin model. The simulation results show that the proposed method can reduce the mean square tracking error by more than 90%. Furthermore, the method is extended to jointly optimise the PID 2 controller and excitation system gains through multiobjective optimisation, leveraging Pareto frontier analysis to balance active power and voltage tracking performance. Simulation results confirm the effectiveness of the proposed method, achieving a 83.46% reduction in relative mean square error of active power, a 47.13% reduction in terminal voltage tracking error, and an 82.78% improvement in the overall scalarized objective.

Hydropower system↗
DOE OSTI2025

Scalable Solution-Processed Electrolyte Membranes with Optimized Microstructure for High-Performance Protonic Ceramic Electrochemical Cells

Proton-conducting electrochemical cells (PCECs) are promising for efficient hydrogen production, but achieving dense, uniform, thin electrolyte layers remains a key challenge, particularly for scalable fabrication. Here, we present a solution-processed deposition approach with a mechanistically optimized slurry for uniform electrolyte formation. By tailoring particle size distribution, solid loading, and solvent/additive balance, we regulated wetting behavior and evaporation kinetics of the electrolyte slurry to promote homogeneous electrolyte particle packing. These features facilitate tight grain boundary contact and early stage neck growth during sintering, eliminating residual porosity, and improving mechanical integrity. The resulting ∼15 μm thick electrolyte shows high density, strong electrode adhesion, and stable interfaces outperforming the previously reported spray-based fabricated electrolyte by about 31% at 600 °C in FC mode. Single cells deliver 0.962 W cm –2 at 600 °C in fuel cell mode and 1.31 A cm –2 at 1.3 V in electrolysis mode, maintaining robust performance over 100 h with negligible degradation (≤0.02% h –1 ) in each mode. Scale-up to 2.5 cm diameter substrates confirmed reproducible densification and geometric stability. This work demonstrates a cost-effective, scalable route where control over particle-fluid interactions and drying dynamics enables a superior electrolyte microstructure and high PCEC performance.

dense microstructure↗
DOE OSTI2020

Mineralogy of the deep lower mantle in the presence of H 2 O

Understanding the mineralogy of the Earth's interior is a prerequisite for unravelling the evolution and dynamics of our planet. Here, we conducted high pressure-temperature experiments mimicking the conditions of the deep lower mantle (DLM, 1800–2890 km in depth) and observed surprising mineralogical transformations in the presence of water. Ferropericlase, (Mg, Fe)O, which is the most abundant oxide mineral in Earth, reacts with H 2 O to form a previously unknown (Mg, Fe)O 2 H x (x ≤ 1) phase. The (Mg, Fe)O 2 H x has a pyrite structure and it coexists with the dominant silicate phases, bridgmanite and post-perovskite. Depending on Mg content and geotherm temperatures, the transformation may occur at 1800 km for (Mg 0.6 Fe 0.4 )O or beyond 2300 km for (Mg 0.7 Fe 0.3 )O. The (Mg, Fe)O 2 H x is an oxygen excess phase that stores an excessive amount of oxygen beyond the charge balance of maximum cation valences (Mg 2+ , Fe 3+ and H + ). This important phase has a number of far-reaching implications including extreme redox inhomogeneity, deep-oxygen reservoirs in the DLM and an internal source for modulating oxygen in the atmosphere.

58 GEOSCIENCES↗
DOE OSTI2022

Resolving Salt-Induced Agglomeration of Laponite Suspensions Using X-ray Photon Correlation Spectroscopy and Molecular Dynamics Simulations

Linking the physics of the relaxation behavior of viscoelastic fluids as they form arrested gel states to the underlying chemical changes is essential for developing predictive controls on the properties of the suspensions. In this study, 3 wt.% laponite suspensions are studied as model systems to probe the influence of salt-induced relaxation behavior arising from the assembly of laponite nanodisks. X-ray Photon Correlation Spectroscopy (XPCS) measurements show that laponite suspensions prepared in the presence of 5 mM concentrations of CaCl2, MgCl2 and CsCl salts accelerate the formation of arrested gel states, with CaCl2 having a significant impact followed by CsCl and MgCl2 salts. The competing effects of ion size and charge on relaxation behavior are noted. For example, the relaxation times of laponite suspensions in the presence of Mg2+ ions are slower compared to Cs+ ions despite the higher charge, suggesting that cation size dominates in this scenario. The faster relaxation behavior of laponite suspensions in the presence of Ca2+ ions compared to Cs+ ions shows that a higher charge dominates the size of the ion. The trends in relaxation behavior are consistent with the cluster formation behavior of laponite suspensions and the electrostatic interactions predicted from MD simulations. Charge balance is achieved by the intercalation of the cations at the negatively charged surfaces of laponite suspensions. These studies show that the arrested gel state of laponite suspensions is accelerated in the presence of salts, with ion sizes and charges having a competing effect on relaxation behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2025

TRAILS Output Files

Overview This data repository contains ZIP files that store compressed versions of the output of running the WaterPaths utility planning and management tool in the DU Re-Evaluation mode (to download the tool, please see this GitHub repository). The tool was used to simulate the six-utility North Carolina Research Triangle problem. Details on the contents of each ZIP file can be seen below. Data details Temporal range: Weekly data for 2,344 weeks from 2015 to 2060 (45 years). Spatial range: Six water utilities in the North Carolina Research Triangle region (0: Chapel Hil/OWASA, 1: Durham, 2: Cary, 3: Raleigh, 4: Pittsboro, and 5: Chatham) File types: CSV and OUT Different solutions available The solution numbers correspond to the different pathway strategies (henceforth referred to as "solutions") discussed in paper's main and supporting text (abstract and link to the paper here). They are as follows: Sol92: The Durham-focused pathway strategy Sol132: The Raleigh-focused pathway strategy Sol140: The regionally-robust pathway strategy Objectives files These files can be accessed by unzipping solXX_objectives_pathways.zip that contains 1,000 Objectives_RDMXX_solsXX_to_XX.csv files. Each CSV file will consist of a row representing all the objective values for that specific solution, while every six columns represents the reliability, restriction frequency, infrastructure net present value ($ mil), peak financial cost, worst-case cost, and unit cost ($ per MG; in that order) for each of the six utilities. There will be 1,000 such files, denoting the performance of the six utilities across the 1,000 deeply uncertain states of the world (DU SOWs). Pathway files These files can be accessed by unzipping solXX_objectives_pathways.zip that contains 1,000 Pathways_sXX_RDMXX.out file. Each OUT corresponds to the set of infrastructure being triggered in a specific DU SOW, and each file will have the name file will consist of four tab-delimited columns that are described as follows: Realization: The realization in which an infrastructure options being triggered utility: The utility currently triggering infrastructure week: The week in which a specific infrastructure option is being triggered infra.: The infrastructure option being triggered If the OUT file contains only the header line, no infrastructure was triggered for that specific DU SOW. Policies files These files can be obtained by unzipping Policies.zip. Each of the 1,000 CSV files within the unzipped folder will contain weekly water use restriction policies for all 1,000 hydroclimatic realizations within a specific DU SOW. The column structure is as follows: 0rest_m: restriction multiplier for utility 0 (values between 0 and 1) 1rest_m: restriction multiplier for utility 1 (values between 0 and 1) 2rest_m: restriction multiplier for utility 2 (values between 0 and 1) 3rest_m: restriction multiplier for utility 3 (values between 0 and 1) 4rest_m: restriction multiplier for utility 4 (values between 0 and 1) 5rest_m: restriction multiplier for utility 5 (values between 0 and 1) 0transf: transfer volume for utility 0 (in MGD) 1transf: transfer volume for utility 1 (in MGD) 2transf: transfer volume for utility 2 (in MGD) 3transf: transfer volume for utility 3 (in MGD) 4transf: transfer volume for utility 4 (in MGD) 5transf: transfer volume for utility 5 (in MGD) Water Sources files These files can be obtained by unzipping WaterSources_subset.zip. Each of the 100 CSV files within the unzipped folder will contain weekly state variables at each water source for all 1,000 hydroclimatic realizations within a specific DU SOW. The column structure is as follows: Xvolume: available water volume from source X (in MGD) Xs_area: surface area of source X (in ACF) Xdemand: demand drawn from a water source from source X (in MGD) Xup_spill: upstream spillage from source X (in MGD) Xww_inflow: wastewater inflow from source X (in MGD) Xcatch_inflow: upstream catchment inflow to source X (in MGD) Xevap: evaporation multiplier for source X (values between 0 and 1) Xds_spill: downstream spillage from source X (in MGD) X_Y_alloc_cap: the allocated capacity from source X to utility Y (values between 0 and 1) X_Y_alloc_dem: the allocated demand from source X to utility Y (values between 0 and 1) Xtrmt_alloc_Y: the allocated treatment capacity from source X to utility Y (values between 0 and 1) Utilities files These files can be obtained by unzipping Utilities_subset.zip. Each of the 100 CSV files within the unzipped folder will contain weekly state variables at each utility for all 1,000 hydroclimatic realizations within a specific DU SOW. The column structure is as follows: Xst_vol: total available storage volume of utility X (in MG) Xcapacity: total storage capacity of utility X (in MG) Xnet_inf: : net inflow for all storage infrastructure for utility X (in MGD) Xst_rof: short term ROF for utility X (values between 0 and 1) Xst_stor_rof: short-term storage ROF for utility X (values between 0 and 1) Xst_trmt_rof: short-term treatment ROF for utility X (values between 0 and 1) Xlt_rof: long-term ROF for utility X (values between 0 and 1) Xlt_stor_rof: long-term storage ROF for utility X (values between 0 and 1) Xlt_trmt_rof: long-term treatment ROF for utility X (values between 0 and 1) Xrest_demand: restricted demand for utility X (in MGD) Xunrest_demand: unrestricted demand for utility X (in MGD) Xunfulf_demand: unfulfilled demand for utility X (in MGD) Xwastewater: wastewater return for utility X (in MGD) Xtreat_capacity: total treatment capacity for utility X (in MG) Xcont_fund: reserve (contingency) fund balance for utility X Xins_pout: insurance payout for utility X (% annual volumetric revenue) Xins_price: insurance price for utility X (% annual volumetric revenue) Xinfra_npv: infrastructure net present value for utility ($mil) Xst_vol: total available storage volume of utility X (in MG) Xdebt_serv: debt service for utility X (usually once per year if the infrastructure is triggered; % annual volumetric revenue) Xstor_vol: total stored volume (in MGD) Xobs_ann_dem: observed annual demand for utility X (in MGD) Xproj_dem: projected annual demand for utility X (in MGD) Xpv_debt_serv: present value of debt service payments for utility X (% annual volumetric revenue) Xgross_rev: gross revenue for utility X ($mil) Acknowledgment IM3 is a multi-institutional effort led by Pacific Northwest National Laboratory and supported by the U.S. Department of Energy's Office of Science as part of research in MultiSector Dynamics, Earth and Environmental Systems Modeling Program.

Artificial Intelligence↗
DOE OSTI2023

Application of a Distributed Element Roughness Model to Additively Manufactured Internal Cooling Channels

Abstract Design for cooling effectiveness in turbine blades relies on accurate models for dynamic losses and heat transfer of internal cooling passages. Metal additive manufacturing (AM) has expanded the design space for these configurations, but can give rise to large-scale roughness features. The range of roughness length scales in these systems makes morphology resolved computational fluid dynamics (CFD) impractical. However, volumetric roughness models can be leveraged, as they have computational costs orders of magnitude lower. In this work, a discrete element roughness model (DERM), based on the double-averaged Navier–Stokes equations, is presented and applied to additively manufactured rough channels, representative of gas turbine blade cooling passages. Unique to this formulation of DERM is a generalized sheltering-based treatment of drag, a two-layer model for spatially averaged Reynolds stresses, and explicit treatment of dispersion. Six different AM rough surface channel configurations are studied, with roughness trough to peak sizes ranging from 15% to 60% nominal channel passage half-width, and the roughness Reynolds number ranges from Rek = 60 to 300. DERM predictions for spatially and temporally averaged mean flow quantities are compared to previously reported direct numerical simulation results. Good agreement in the mean velocity profiles, stress balances, and drag partitions are observed. While DERM models are typically calibrated to specific deterministic roughness morphologies at comparatively small roughness Reynolds numbers, the present more generalized DERM formulation has wider applicability. Here, it is demonstrated that the model can accommodate random roughness of large scale, typical of AM.

Engineering↗
DOE OSTI2024

River–aquifer interactions enhancing evapotranspiration in a semiarid riparian zone: A modelling study

The hydrologic flows across the river–aquifer interface play an important role in groundwater dynamics and biogeochemical reactions within the subsurface; however, little is known about the effects of river–aquifer interactions on land surface processes. In this study, we developed a fully coupled three-dimensional (3D) land surface and subsurface model at a high resolution (~1 km) that accounts for high-frequency hydrologic exchange flow conditions to investigate how river–aquifer interactions modulate surface water budgets in the Upper Columbia-Priest Rapids watershed, a typical semiarid watershed located in the northwestern United States where river stage fluctuates in response to reservoir releases changing. Our results show that the spatiotemporal dynamics of river–aquifer interactions are highly heterogeneous, driven mainly by river-stage fluctuations. Adding 6.64 × 10 6 m 3 year –1 of water over the watershed from the river to groundwater owing to the lateral flow, river–aquifer interactions led to an increase in soil evaporation and transpiration supplied by higher soil moisture content, particularly in deeper subsurface. In a hypothetic future scenarios where a 5-m rise in river stage was assumed, the hydrologic flow exchange rates were intensified, resulting in higher surface water over the entire watershed. Overall, lateral flow induced by river–aquifer exchanges leads to an increase in evapotranspiration of ~75% in the historical period and of ~83% in the hypothetical future scenario. Finally, our study demonstrates the potential of coupled model as an effective tool for understanding river–aquifer–land surface interactions, and indicates that river–aquifer interactions fundamentally alter the water balance of the riparian zone for the semiarid watershed and will likely become more frequent and intense in the future under the effects of climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Sustainable Functional Epoxies through Boric Acid Templating

Thermoset polymers (e.g. epoxies, vulcanizable rubbers, polyurethanes, etc.) are crosslinked materials with excellent thermal, chemical, and mechanical stability; these properties make thermoset materials attractive for use in harsh applications and environments. Unfortunately, material robustness means that these materials persist in the environment with very slow degradation over long periods of time. Balancing the benefits of material performance with sustainability is a challenge in need of novel solutions. Here, we aimed to address this challenge by incorporating boronic acid-amine complexes into epoxy thermoset chemistries, facilitating degradation of the material under pH neutral to alkaline conditions; in this scenario, water acts as an initiator to remove boron species, creating a porous structure with an enhanced surface area that makes the material more amenable to environmental degradation. Furthermore, the expulsion of the boron leaves the residual pores rich in amines which can be exploited for CO 2 absorption or other functionalization. We demonstrated the formation of novel boron species from neat mixing of amine compounds with boric acid, including one complex that appears highly stable under nitrogen atmosphere up to 600 °C. While degradation of the materials under static, alkaline conditions (our “trigger”) was inconclusive at the time of this writing, dynamic conditions appeared more promising. Additionally, we showed that increasing boronic acid content created materials more resistant to thermal degradation, thus improving performance under typical high temperature use conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Performance Testing of a Moving-Bed Gasifier Using Coal, Biomass, and Waste Plastic Blends with Washed and Unwashed Legacy Coals and Other Waste Fuels to Generate White Hydrogen

The objective of this effort, primarily funded by the United States Department of Energy (DOE), and led by the Electric Power Research Institute, Inc. (EPRI), with support by Hamilton Maurer International (HMI) and Sotacarbo S.p.A. (Sotacarbo), has been to qualify coal, biomass, and plastic waste blends based on performance testing of selected fuel pellet compositions in a pilot-scale updraft moving-bed (UDMB) gasifier. The testing provided relevant data to advance the commercial-scale design of the moving-bed gasifier to be able to successfully use these feedstocks to produce hydrogen. In particular, the effects of waste plastics on feedstock development (i.e., blending and pelletizing) and the resulting products (i.e., syngas compositions, organic condensate production, and ash characteristics) are the focus. The gasifier used for testing is HMI’s moving-bed gasifier, which has been proven capable of gasifying nearly all coal ranks. It has also shown the ability in prior testing work to gasify wood chips (biomass). However, mixtures of these fuels with plastic wastes have not been prepared and gasified together. The three feedstocks were densified and pelletized by California Pellet Mill (CPM) to meet the feedstock size required by Sotacarbo’s 30mm ID UDMB gasifier, under contract to HMI. The technical tasks and results from this two-year research project included: (1) Feed Procurement and Preparation: Nine different tri-fuel pellets were prepared from varying compositions of fresh mined PRB coal, corn stover biomass, and car fluff waste plastics. Tri-fuel pellets were produced by CPM and shipped to Sotacarbo’s test facility in Carbonia, Sardinia, Italy. (2) Test Plan Development: A test plan was created to define the test runs to be performed. The test plan detailed the different UDMB gasification tests to be performed in Sotacarbo’s 12-inch ID pilot scale gasifier, the process monitoring instrumentation used, and the extractive samples recovered for analysis of the total gasification process mass and energy balance. (3) Gasifier Testing: Nine different gasification runs were performed in the pilot-scale gasifier at Sotacarbo using nine different fuel feedstock compositions generated from varying mixtures of PRB coal, biomass, and plastic wastes. The testing generated performance data on gasification reaction efficiency and performance, yielding relevant data for models used to scale up the gasifier design. This task also included work to refurbish and reassemble the pilot gasifier at Sotacarbo and perform a baseline 100% PRB coal run. (4) Data Analysis and Reporting: Review of the data, determination of figures of merit, and interpretation of the results are reported in the project’s final report, published in March 2024. The results show that all tri-fuel pellets gasified well and maintained structural integrity throughout the gasification process. The syngas generated can be shifted to hydrogen by using commercial syngas shifting technologies. (5) High Fidelity computational fluid dynamics (CFD) Simulation: The National Energy Technology Laboratory (NETL) team performed CFD simulations of the UDMB gasifier for two of the tri-fuel pellets gasified in Sotacarbo’s pilot scale gasifier. The kinetic mechanisms for the pyrolysis of each constituent, PRB coal, corn stover biomass, and waste plastics are based on thermogravimetric analysis performed by Sotacarbo. The gasification model was validated by comparing the predicted syngas composition at the exit of the gasifier with the measured syngas composition. In addition, the reactor’s measured internal temperature profile agreed well with the predicted internal reactor temperature profile. These results validate that the model can be used to predict the performance of the updraft moving bed gasifier for different feedstocks and operating conditions. This paper summarizes the results of the completed work in which the pelletizing procedure was validated to ensure the viability of the tri-fuel pellets for the gasification runs performed at Sotacarbo’s 30 mm UDMB gasifier. The gasification performance data from this series of nine runs will enable modeling of a full-scale HMI industrial scale gasifier supporting both combined heat and power, and Hydrogen production from coal (both fresh mined and legacy) combined with various biomass and waste plastics. Additionally, plans and progress on a follow-up project, being executed by the same project team, will be presented. In this project, a total of twenty (20) different feedstocks are being prepared from varying compositions of biomass (both woody biomass and corn stover) with a mixture of legacy coal waste, plastic waste, and refuse-derived fuel (RDF). The testing will provide information on gasification reaction efficiency/performance, yielding relevant data for models used to scale up the gasifier design to 50 megawatt electric (MWe) (equivalent hydrogen production). Tests will also be performed on a bench-scale fluidized-bed gasifier for comparison purposes. The results of this testing will be used to specify the range of feedstock blends that can be successfully gasified as well as quantify gasifier outputs based on specific blends.

08 HYDROGEN↗
DOE OSTI2025

CO 2 Capture Characteristics of Hyperbranched Poly(alkylene imine): A Molecular Dynamics Simulation Approach

This study explores the CO 2 capture characteristics of hyperbranched poly(ethylenimine) (HB-PEI) and poly- (propyleneimine) (HB-PPI) through molecular dynamics simulations using density functional theory-calibrated force fields. Key features such as density, free volume, glass transition temperature, CO 2 /H 2 O distribution, and molecular diffusion are systematically investigated to elucidate structure−function relationships under dry and hydrated conditions. HB-PEI demonstrates a slightly higher density and lower free volume compared to HB-PPI yet shows superior CO 2 capture due to the high amine concentration. Glass transition analysis indicates a higher thermal mobility in HBPEI, enhancing the CO 2 diffusivity. Pair correlation and coordination analyses confirm a stronger affinity of CO 2 with primary and secondary amines, particularly in hydrated environments where water competes with CO 2 for binding sites. Despite its more compact structure, HB-PEI outperformed HB-PPI in CO 2 and H 2 O transport, as confirmed by higher diffusion coefficients across all hydration levels. These findings highlight a critical balance among polymer architecture, amine accessibility, and hydration in designing next-generation solid amine sorbents for efficient direct air capture applications.

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

HydraGNN_Predictive_GFM_2024 - Ensemble of predictive graph foundation models for ground state atomistic materials modeling

We provide the ensemble of fifteen pre-trained graph foundation models (GFMs) for atomistic materials modeling applications. Each one of the fifteen GFMs has been trained on five open-source datasets that (once aggregated) amount to over 154 million atomistic structures, which cover over two-thirds of the natural elements of the periodic table and that comprises a broad set of organic and inorganic compounds. This vast set of atomistic structures comprises ground state configurations that are dynamically stable (i.e., equilibrated structures with atomic forces approximately close to zero values) as well as dynamically unstable structures (i.e., non-equilibrium structures with non-negligible non-zero values of atomic forces). The ensemble of datasets aggregated does NOT include excited states. The datasets have been curated to remove atomistic structures with spectral norm of the force tensor above 100 eV/angstrom. Moreover, a linear term of the energy was computed for each dataset using a linear regression model that uses the chemical concentration of each natural element as regressor. The linear term predicted by the linear regression model has been subtracted from each original energy value to perform a re-alignment of the energy values across different electronic structures approximation theories performed to generate the diverse multi-source, multi-fidelity datasets. The folder "ADIOS_files" 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 "ADIOS_files" directory contains 6 sub-directories named as follows: - ANI1x-v3.bp - MPTrj-v3.bp - OC2020-20M-v3.bp - OC2020-v3.bp - OC2022-v3.bp - qm7x-v3.bp Each sub-directory contains the pre-processed datasets converted in Adaptable I/O System (ADIOS) format (https://www.exascaleproject.org/research-project/adios/) that have been used to the development, training, and performance testing of the ensemble go predictive graph foundation models. Each GFM was developed using HydraGNN (https://github.com/ORNL/HydraGNN) as underlying graph neural network (GNN) architecture. The multi-task learning (MTL) capability of HydraGNN was used to simultaneously train the GFMs on labeled values for direct predictions of energy (a total system property of an atomistic structure that measures the chemical stability) and atomic forces (an atomic level property of an atomistic structure that measures the dynamical stability). The hyper parameters of the GFM have been tuned using scalable hyperparameter optimization (HPO) algorithms implemented in the software DeepHyper (https://github.com/deephyper/deephyper). The pre-training of each HPO trial was performed using distributed data parallelism (DDP) to scale the training across 128 compute nodes of the exascale OLCF supercomputer Frontier. Each HPO trial was trained only for 10 epochs and an early stopping was performed to avoid wasting significant computational resources on GNN architectures that were clearly underperforming. For each HPO trial, the 'omnistat' tool developed by (AMD Research - Advanced Micro Device) was used to measure the total energy consumption in kWh. The ensemble of GFMs was obtained by selecting the fifteen best performing HPO trials. Four models have been selected for their clear advantage in accuracy, and these are the GFMs with IDs 229, 156, 147, 260. Additional eleven models have been selected based on judicious balance between accuracy and energy consumption needed for training, and these are the GFMs with IDs 165, 78, 137, 1, 175, 171, 181, 67, 179, 167, 351. Each selected GFM of the ensemble was continued to cumulate a total of at most 30 epochs. In some cases, the total number of epochs actually performed was les than 30 due to two combined factors: (1) the size of the GFM (i.e., the number of model parameters to train) and (2) the total wall-clock time for which the computational resources could be allocated on OLCF-Frontier. The "Ensemble_of_models" directory contains 15 sub-directories named as follows: - gfm_0.229 - gfm_0.156 - gfm_0.147 - gfm_0.260 - gfm_0.165 - gfm_0.78 - gfm_0.137 - gfm_0.1 - gfm_0.175 - gfm_0.171 - gfm_0.181 - gfm_0.67 - gfm_0.179 - gfm_0.167 - gfm_0.351 Each one of these sub-directories refers to one of the fifteen HPO trials that have been selected to continue the pre-training with at most 30 epochs. With each sub-directory associated with a specific HPO trial, the following files can be found: - config.json: file for argument parsing to develop and train an HydraGNN architecture - gfm_0.ID_epoch_N.pk: file with model parameters for HPO ID trial after N epochs of training The ensemble of fifteen GFM architectures was used for (1) ensemble averaging to stabilize the predictions of energy and atomic forces after pre-training for post-processing analysis and (2) ensemble uncertainty quantification (UQ). The code used to develop, pre-train, and load the pre-trained models for post-processing analysis is available on the ORNL-GitHub at the following link: https://github.com/ORNL/HydraGNN/tree/Predictive_GFM_2024

36 MATERIALS SCIENCE↗
DOE OSTI2025

The role of thermal stratification on the co‐spectral properties of momentum transport above an Amazonian forest

The influence of thermal stratification on the turbulent kinetic energy balance has been widely studied; however, its influence on the turbulent stress remains less explored in the presence of tall vegetated canopies and less ideal micrometeorological conditions. Here, the impact of thermal stratification on turbulent momentum flux is considered in the roughness sublayer (RSL) and the atmospheric surface layer (ASL) using the Amazon Tall Tower Observatory (ATTO) in Brazil. A scalewise co‐spectral budget (CSB) model is developed using standard closure schemes for the pressure–velocity decorrelation. The CSB revealed that the co‐spectrum ${F}_{wu}\left({k}_x\right)$ between longitudinal (u') and vertical (w') velocity fluctuations is impacted by the energy spectrum of the vertical velocity ${E}_{ww}\left({k}_x\right)$ and the much less studied longitudinal heat‐flux co‐spectrum ${F}_{u{\theta}_{\mathrm{v}}}\left({k}_x\right)$, where ${\theta}_{\mathrm{v}}^{\prime }$ are temperature fluctuations and ${k}_x$ is the longitudinal wavenumber. Under stable, very stable, and dynamic–convective conditions, the scaling exponent ${F}_{wu}\left({k}_x\right)$ in for the inertial subrange (ISR) scales is dominated by ${F}_{u{\theta}_{\mathrm{v}}}\left({k}_x\right)$ instead of ${E}_{ww}\left({k}_x\right)$. A near ${k}_x^{-7/3}$scaling in ${F}_{u{\theta}_{\mathrm{v}}}\left({k}_x\right)$ robust to large variations in thermal stratification is found, whereas the Kolmogorov ISR scaling for ${E}_{ww}\left({k}_x\right)\sim {k}_x^{-5/3}$ is not found. The scale‐dependent decorrelation time between u' and w' is dominated by ${\epsilon}^{-1/3}{k}_x^{-2/3}$ in the ISR, but is nearly constant for eddies larger than the vertical velocity integral scale, regardless of stability. Implications of these findings for generalized stability correction functions that are based on the turbulent stress budget instead of the turbulent kinetic energy budget are discussed.

canopy turbulence↗
DOE OSTI2022

Decadal forest soil respiration following stem girdling

Soil respiration (Rs) represents the largest flux of carbon (C) from forests to the atmosphere, but the long-term influence of phloem-disrupting disturbance on Rs is poorly understood, limiting robust forecasts of ecosystem C balance. Using a decade of observations from the Forest accelerated succession experiment (FASET), we examined relationships among Rs, soil temperature, soil moisture, and leaf area index (LAI) following the stem girdling-induced mortality of 40% of all canopy trees within a 39-ha area. Mean annual Rs declined by about 20% relative to the control two years after disturbance, but recovered to near pre-disturbance values within five years; this reduction correlated with LAI losses and lower Rs temperature sensitivity (i.e., Q 10 ), with the latter counteracting soil warming caused by partial canopy defoliation. These observations are consistent with progressive reductions in belowground labile C causing reductions in Rs. In conclusion, we conclude that the effects of stem girdling on Rs (1) were not immediate, occurring two years after the treatment, (2) were primarily influenced by biotic rather than soil microclimate changes, and (3) persisted for nearly a decade but were temporally dynamic, underscoring the value of long-term experiments.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

City-Scale Building Anthropogenic Heating during Heat Waves

More frequent and longer duration heat waves have been observed worldwide and are recognized as a serious threat to human health and the stability of electrical grids. Past studies have identified a positive feedback between heat waves and urban heat island effects. Anthropogenic heat emissions from buildings have a crucial impact on the urban environment, and hence it is critical to understand the interactive effects of urban microclimate and building heat emissions in terms of the urban energy balance. Here we developed a coupled-simulation approach to quantify these effects, mapping urban environmental data generated by the mesoscale Weather Research and Forecasting (WRF) coupled to Urban Canopy Model (UCM) to urban building energy models (UBEM). We conducted a case study in the city of Los Angeles, California, during a five-day heat wave event in September 2009. We analyzed the surge in city-scale building heat emission and energy use during the extreme heat event. We first simulated the urban microclimate at a high resolution (500 m by 500 m) using WRF-UCM. We then generated grid-level building heat emission profiles and aggregated them using prototype building energy models informed by spatially disaggregated urban land use and urban building density data. The spatial patterns of anthropogenic heat discharge from the building sector were analyzed, and the quantitative relationship with weather conditions and urban land-use dynamics were assessed at the grid level. The simulation results indicate that the dispersion of anthropogenic heat from urban buildings to the urban environment increases by up to 20% on average and varies significantly, both in time and space, during the heat wave event. The heat dispersion from the air-conditioning heat rejection contributes most (86.5%) of the total waste heat from the buildings to the urban environment. We also found that the waste heat discharge in inland, dense urban districts is more sensitive to extreme events than it is in coastal or suburban areas. The generated anthropogenic heat profiles can be used in urban microclimate models to provide a more accurate estimation of urban air temperature rises during heat waves.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

Coupling plant litter quantity to a novel metric for litter quality explains C storage changes in a thawing permafrost peatland

Abstract Permafrost thaw is a major potential feedback source to climate change as it can drive the increased release of greenhouse gases carbon dioxide (CO 2 ) and methane (CH 4 ). This carbon release from the decomposition of thawing soil organic material can be mitigated by increased net primary productivity (NPP) caused by warming, increasing atmospheric CO 2 , and plant community transition. However, the net effect on C storage also depends on how these plant community changes alter plant litter quantity, quality, and decomposition rates. Predicting decomposition rates based on litter quality remains challenging, but a promising new way forward is to incorporate measures of the energetic favorability to soil microbes of plant biomass decomposition. We asked how the variation in one such measure, the nominal oxidation state of carbon (NOSC), interacts with changing quantities of plant material inputs to influence the net C balance of a thawing permafrost peatland. We found: (1) Plant productivity (NPP) increased post‐thaw, but instead of contributing to increased standing biomass, it increased plant biomass turnover via increased litter inputs to soil; (2) Plant litter thermodynamic favorability (NOSC) and decomposition rate both increased post‐thaw, despite limited changes in bulk C:N ratios; (3) these increases caused the higher NPP to cycle more rapidly through both plants and soil, contributing to higher CO 2 and CH 4 fluxes from decomposition. Thus, the increased C‐storage expected from higher productivity was limited and the high global warming potential of CH 4 contributed a net positive warming effect. Although post‐thaw peatlands are currently C sinks due to high NPP offsetting high CO 2 release, this status is very sensitive to the plant community's litter input rate and quality. Integration of novel bioavailability metrics based on litter chemistry, including NOSC, into studies of ecosystem dynamics, is needed to improve the understanding of controls on arctic C stocks under continued ecosystem transition.

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

An open-source framework for balancing computational speed and fidelity in production cost models

Studies of bulk power system operations need to incorporate uncertainty and sensitivity analyses, especially around exposure to weather and climate variability and extremes, but this remains a computational modeling challenge. Commercial production cost models (PCMs) have shorter runtimes, but also important limitations (opacity, license restrictions) that do not fully support stochastic simulation. Open-source PCMs represent a potential solution. They allow for multiple, simultaneous runs in high-performance computing environments and offer flexibility in model parameterization. Yet, developers must balance computational speed (i.e. runtime) with model fidelity (i.e. accuracy). In this paper, we present Grid Operations (GO), a framework for instantiating open-source, scale-adaptive PCMs. GO allows users to search across parameter spaces to identify model versions that appropriately balance computational speed and fidelity based on experimental needs and resource limits. Results provide generalizable insights on how to navigate the fidelity and computational speed tradeoff through parameter selection. We show that models with coarser network topologies can accurately mimic market operations, sometimes better than higher-resolution models. It is thus possible to conduct large simulation experiments that characterize operational risks related to climate and weather extremes while maintaining sufficient model accuracy.

42 ENGINEERING↗
DOE OSTI2024

Dataset for: Price Controls for Scarcity Events in Real-Time and Transactive Energy Systems

Real time pricing (RTP) is often promoted as a mechanism to improve the economic efficiency of the electricity system. However, many regulators have been hesitant to adopt RTP due to concerns about exposing customers to extreme price swings. To balance these concerns, this paper proposes a methodology for establishing price controls, based on the supply of demand-side flexibility in the system. As an illustrative example, we measure price responsiveness using an agent-based simulation model that is representative of the ERCOT market. The model is composed of a distribution feeder that has 250 customers with active agents controlling their HVAC systems in response to the historical ERCOT RTP with an artificially added high-price event. These agents are subjected to increasing electricity prices during the event, which we then use to create a supply curve for demand-side resources in our modeled scarcity event. We set potential price caps at points on the supply curve where customers’ have exhausted their flexible capacity. Using historical prices, we examine the systemic costs of these price caps, and present regulatory options for recouping them. Utilities and regulators interested in limiting consumer risk from dynamic pricing can utilize these methods to develop rate structures and encourage conservation. The attached data upload allows for the duplication or modification of the analysis performed in this study.

Kerby, Jessica R↗