Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “surface complexation modeling”

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.

1,295 records · Page 72
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Terrain-Influenced Winds and Fire-Fire Interactions in Wildland Fire Simulations [Dissertation]

Ensemble-based approaches to prescribed fire planning cannot be supported by computational fluid dynamics based models like FIRETEC and the Wildland-Urban Interface Fire Dynamics Simulator (WFDS) because they are too computationally expensive and cannot leverage large eddy simulation approaches like CAWFE and WRF-SFIRE because they have too coarse of resolution. QUIC-Fire was developed to fill this gap but it cannot currently address complex terrain, that is typical for instance in the Western United States. This dissertation describes a variety of improvements made to QUIC-Fire and its various incorporated algorithms in an effort to make it a viable tool in simulating wildland and prescribed fires on terrain. The modifications made to QUIC-Fire are described in three chapters. The first chapter describes the extension of the diagnostic wind model QUIC-URB, the wind engine of QUIC-Fire, to a terrain-following coordinate system. The terraininfluenced winds it generates are analyzed and compared. In particular, this chapter presents the mathematical derivation of the wind solver leading to a linear system of equations that are solved through the successive over-relaxation method. The model is validated against a standard test used in previous works (the Askervein Hill) and against a new dataset from measurements in the Socorro Mountains, New Mexico. The terrain-following implementation captures the correct phenomenology for the isolated Askervein Hill, with a wind speed up at the top of the hill. The model agrees well with measurements on the upwind side of the peak, but overestimates speed-up on the downwind side of the hill. This is due to the inability of the model to generate flow separation and wake-eddy dynamics. On a common laptop, the divergence-free wind field is obtained in 6 s, making the solver appealing for coupled fire-atmosphere simulations. The Socorro Mountain is highly complex, with many cliff faces, peaks, and valleys. Although the model captures the magnitude and direction of inlet and outlet areas of the domain, it performs rather poorly in the valley region and in the regions near the steep cliffs. Hence, the model shows good agreement with data in areas of open sloped terrain but lacks in areas where flow separation and thermally driven effects may be present (neither effect is addressed in this work). In the second chapter the implementation of the terrain-following version of QUIC-URB into QUIC-Fire, and the necessary changes needed to include terrain are described. No changes to the underlying fire spread algorithm are made other than what is required to correctly account for the inclusion of terrain. Previously published FIRETEC results that use five different topographies that share the same centerline profile are compared to simulation results from the modified QUIC-Fire that use the same topographies and fuels. QUIC-Fire results show overall similar behaviors in terms of how the topographies affect fire shapes and trends in spread rates. Due to the terrain-following version of QUIC-URB being unable to generate flow separations at the crest of hills, fire spread rates in these regions across all topographies are over-predicted when compared to FIRETEC. Lateral fire growth shows similar trends with FIRETEC between topographies but does not capture the increase in spread due to a diagonal interface between grassland and forested fuel region of the domain. These results suggest that there are three algorithms within QUIC-Fire that could use improvement: how flame tilt angle is accounted for, the incorporation of non-local drag effects, and the inclusion of the wake-eddy parameterizations that are used in QUIC-URB. Lastly, the third chapter describes a modification to the initial guess used for the calculation of the QUIC-URB mass-conserved wind solution during fire simulations. The modification is aimed at improving fire-fire interactions in QUIC-Fire simulations. The modification consists of using the solution from the previous timestep as the starting point for the calculation of the solution for the next timestep. Fire-fire interactions is greatly improved by the change but a new source of error is introduced. Due to how plumes are modelled in QUIC-Fire the new solution contains errors where gaps in the plume structure are present. However, these errors are mostly limited to the upper atmosphere, where they do not affect fire behavior at the surface, and their magnitude isn’t significant enough to discount the amount of new fire phenomenology now captured in QUIC-Fire with the change.

58 GEOSCIENCES↗
DOE OSTI2024

Soil Carbon and Nitrogen Elemental Analysis and Middle Infrared Spectroscopy from the Kougarok Fire Complex, Seward Peninsula, Alaska, 2022

In September and October of 2022, soil samples were collected at the Kougarok Fire Complex near mile marker 86 of the Nome-Taylor Highway on the Seward Peninsula of Alaska. This study site was chosen due to its unique fire history, as the Kougarok Fire Complex has experienced multiple wildfires since 1971. As the Arctic warms, the risk of Arctic tundra fires continues to increase due to warmer summer temperatures and higher frequency of lightning. Burned soil carbon or pyrogenic carbon (PyC) is an important component of C cycling after wildfire, and one that is often overlooked in tundra systems where wildfires are historically rare. To better understand PyC signatures and quantify PyC presence in post-regeneration permafrost regions, soils were sampled from soil pits within the 1971, 2002, and 2019 burn sites, as well as two unburned control sites. At each site, three soil pits were dug to the permafrost table. Soil samples were collected with a trowel from the face of each pit at 10 cm increments down to the permafrost table. Thaw depth, maximum vegetation height, vegetation species composition, and O horizon depth were also collected at each soil pit. An HS2 Hydrosense II Handheld Soil Moisture Sensor was used to collect volumetric soil moisture content at each 10 cm sampling increment. This dataset includes one *.csv of middle infrared spectroscopy measurements of soil samples, one *.csv of field observations, one *.csv of carbon and nitrogen analysis of soil samples, and one *.kml of sampling locations.The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic) was a 15-year research effort (2012-2027) to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy’s Office of Biological and Environmental Research.The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska.Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy’s Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Models, data, and scripts associated with “Prediction of Distributed River Sediment Respiration Rates using Community-Generated Data and Machine Learning”

This data package is associated with the publication “Prediction of Distributed River Sediment Respiration Rates using Community-Generated Data and Machine Learning’’ submitted to the Journal of Geophysical Research: Machine Learning and Computation (Scheibe et al. 2024). River sediment respiration observations are expensive and labor intensive to obtain and there is no physical model for predicting this quantity. The Worldwide Hydrobiogeochemisty Observation Network for Dynamic River Systems (WHONDRS) observational data set (Goldman et al.; 2020) is used to train machine learning (ML) models to predict respiration rates at unsampled sites. This repository archives training data, ML models, predictions, and model evaluation results for the purposes of reproducibility of the results in the associated manuscript and community reuse of the ML models trained in this project. One of the key challenges in this work was to find an optimum configuration for machine learning models to work with this feature-rich (i.e. 100+ possible input variables) data set. Here, we used a two-tiered approach to managing the analysis of this complex data set: 1) a stacked ensemble of ML models that can automatically optimize hyperparameters to accelerate the process of model selection and tuning and 2) feature permutation importance to iteratively select the most important features (i.e. inputs) to the ML models. The major elements of this ML workflow are modular, portable, open, and cloud-based, thus making this implementation a potential template for other applications. This data package is associated with the GitHub repository found at Please see the file level metadata (flmd; “sl-archive-whondrs_flmd.csv”) for a list of all files contained in this data package and descriptions for each. Please see the data dictionary (dd; “sl-archive-whondrs_dd.csv”) for a list of all column headers contained within comma separated value (csv) files in this data package and descriptions for each. The GitHub repository is organized into five top-level directories: (1) “input_data” holds the training data for the ML models; (2) “ml_models” holds machine learning models trained on the data in “input_data”; (3) “scripts” contains data preprocessing and postprocessing scripts and intermediate results specific to this data set that bookend the ML workflow; (4) “examples” contains the visualization of the results in this repository including plotting scripts for the manuscript (e.g., model evaluation, FPI results) and scripts for running predictions with the ML models (i.e., reusing the trained ML models); (5) “output_data” holds the overall results of the ML model on that branch. Each trained ML model resides on its own branch in the repository; this means that inputs and outputs can be different branch-to-branch. Furthermore, depending on the number of features used to train the ML models, the preprocessing and postprocessing scripts, and their intermediate results, can also be different branch-to-branch. The “main-*” branches are meant to be starting points (i.e. trunks) for each model branch (i.e. sprouts). Please see the Branch Navigation section in the top-level README.md in the GitHub repository for more details. There is also one hidden directory “.github/workflows”. This hidden directory contains information for how to run the ML workflow as an end-to-end automated GitHub Action but it is not needed for reusing the ML models archived here. Please the top-level README.md in the GitHub repository for more details on the automation.

13C↗
DOE OSTI2022

AI-Enabled Discovery and Physics-Based Optimization of Energy Efficient Processing Strategies for Advanced Turbine Alloys (Final Technical Report)

In this project, the multi-organizational team of academic and industrial researchers from the University of Kentucky an aerospace and energy generation OEM partner has leveraged novel Digital Process Twin (DPT) models of process/structure interactions (i.e., process-induced surface integrity) to advance a paradigm of fully integrated computational materials engineering (ICME). Using efficient process models as the core of a digital process simulator for a reinforcement learning algorithm, the team has integrated industrial data and metrics of structure/performance/energy relationships and manufacturing-related energy metrics to optimize dynamic processing parameters for significantly improved life-cycle energy efficiency of advanced γ-TiAl low-pressure turbine (LPT) alloys, as indicated by a set of design relevant parameters (e.g., residual stresses and scrap rate). The key objective and anticipated outcome of the project was at least a 10% reduction in life-cycle embodied energy for a recently developed, γ-TiAl low-pressure turbine (LPT) alloy and nickel-based superalloy Inconel 718, through the adoption of the proposed AI-enabled process optimization approach. The final project outcomes significantly exceeded this original target, realizing manufacturing-related energy efficiency improvements of more than 130% for TiAl and up to 80% for Inconel 718. Rather than following the prevailing and highly inefficient empirical paradigm, the proposed study demonstrated the feasibility of adopting a digital, physics-based process design and optimization paradigm. The recurring need for manual intervention, rework, reinspection causes significant production bottlenecks and unnecessary expense associated with delivering the requisite component quality. The OEM partner, and turbine industry in general, expect to reap significant cost and resource savings if an AI-optimized set of parameters can be applied to specific machining operations. The technical scope of the proposed project involved the paving of a realistic path towards model-based and AI-enabled Integrated Computational Materials Engineering (ICME), and away from inefficient empirical process optimization and legacy manufacturing practices, which are no longer able to efficiently process novel high-performance turbine alloy materials. The project team will address the fundamental knowledge gap that currently exists within the ICME paradigm with respect to the process/structure/performance/energy impacts of finishing processes. While significant resources have been devoted to the ‘early stages’ of manufacturing, such as alloy design, primary and secondary processing, finishing processes have not been adequately integrated within ICME. To provide an actionable path towards model-based finishing process design (e.g., machining, burnishing, grinding, polishing), we will employ a novel AI-enabled process optimization paradigm, based on a computationally efficient, physics-based process simulator. Through limited experimental work to calibrate and validate our process simulator model via an advanced in-situ characterization technique and process optimization via reinforcement learning, the project will seek to demonstrate a viable alternative to the inefficient ‘legacy’ processing strategies, empirical testing and broad scope machining learning approaches, all of which fail to adequately consider complex process physics. The project team has identified an intermetallic γ-TiAl LPT alloy, which is currently being used as part of the OEM partner’s advanced gas turbine designs. This particular alloy poses significant manufacturing challenges during finishing operations, which limit the degree to which the current turbine design can be manufactured in an energy- and cost-efficient manner. Empirical testing and numerical modeling efforts to optimize processing parameters for γ-TiAl have not been able to resolve these manufacturing challenges, so the proposed physics-based AI-enabled optimization technology would offer a truly novel and transformative capability. The multi-organizational team of academic and industry experts from the UKY and the OEM partner will work together closely to demonstrate the analytical and experimental critical function and characteristic proof of concept of this novel approach.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2023

Department of Energy’s Atmospheric System Research (ASR) Program’s Workshop on the Future of Atmospheric Large Eddy Simulation (LES) (Workshop Report)

Large-eddy simulation (LES) is used as a tool to understand physical processes such as turbulence, aerosols, clouds, precipitation, radiation, the interactions among all these, and their interactions with the underlying surface. Over the next 10 years, LES will drive fundamental progress in open scientific questions in these areas as LES is increasingly used to gain understanding of complex interacting physical processes involving atmospheric turbulence. This growth will be driven both by scientific demand and the expansion of computational resources needed to conduct LES, and the form that the growth takes will largely be determined by how computational resources are leveraged for scientific gain. In particular, we suggest that computational resources are likely to be leveraged in two separate but not necessarily distinct ways. On one hand, growth in computational resources will allow LES to be made more routine, that is, performed more frequently, while on the other hand, the computational expense (measured in total floating point operations) afforded to individual LES will expand dramatically, allowing simulations to increase in both domain size and resolution as well as physical detail. Current U.S. Department of Energy (DOE) projects such as LES ARM Symbiotic Simulation and Observation Activity (LASSO) are leading the way in conducting routine LES, building large, public databases that are accessible for data science, sensitivity studies, and training for machine learning. LES will also become more routine as it becomes more accessible for individual researchers to address their scientific questions of interest. Scientific questions addressed by LES over the next 10 years are likely to include cloud organization and aggregation; aerosol cloud interactions and atmospheric chemistry (including geo-engineering); urban-scale LES; atmospheric extreme events, ranging from small-scale severe weather to wildfires; and ocean-wave-atmosphere interactions. Further LES-related research will likely grow significantly in areas related to societal impact studies of air quality and extreme weather events, applications to renewable energy forecasting and resource assessment, and aid in decision-making processes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

First-Principles Simulation of Beam-Induced Processes Underlying Atomic Manipulation in Electron Microscopes

The development of experimental methods and apparatuses capable of promoting atomically precise material manipulations holds great promise for realizing the ultimate limit of feature miniaturization in materials and devices. The ability to modify materials atom by atom is anticipated to usher in new technologies in areas as diverse as separation science, medicine, and quantum information science. Historically, scanning probe-based techniques have been the most prominent approaches in this space. However, these methods are best suited for the manipulation of surface-exposed regions of materials, as the strong perturbations required for bond scission are delivered most effectively to atoms in the near-proximity to the scanning probe. In contrast, convergent electron beams with energies tuned slightly below the threshold for inducing irreversible knock-on damage have recently been employed (within scanning transmission electron microscopy) to promote atomic-scale bond rearrangements in various beam-stable solids. Currently, however, the efficiency and selectivity of beam-induced atomic manipulation processes with focused electron beams are such that long irradiation times are required to induce a desired atomic rearrangement. With a better understanding of the underlying physics dictating the outcome of a given irradiation event, methods can be devised to improve the efficiency of these techniques so that their promise can be fully realized through widespread adoption.To this end, this Account details our recent efforts to develop and apply tractable first-principles simulation approaches for studying the response of materials to electric beam-like external electric potentials applied in real space. We briefly review the concepts and capabilities in the area of atomically precise materials manipulation and review the early demonstrations of accomplishments in this area, focusing on studies using scanned convergent electron beam probes in particular. We expound upon the depth of the challenge and identify critical shortcomings of theoretical methods that have previously been employed in the simulation of beam-induced processes. We then describe the computational methods that we have generalized from the concepts and tools most commonly applied to the study of molecular photochemistry and how our adaptations of these methods can be employed to capture the relevant dynamical phenomena for beam-induced processes ranging from the initial electron scattering to the ensuing multistate reactions. Here, we contextualize these methods within the current state of the art in this area, which has historically focused primarily on the simulation of inelastic image formation in the electron microscope for the purpose of interpreting the results of quantitative electron microscopy experiments. We demonstrate that the spatial distribution of state-specific excitation rates due to the presence of an external (probe) electric charge is inhomogeneous, such that irradiation at particular locations in materials can favor specific electronic transitions (and disallow others). In addition to the potential for excited-state reaction pathways to be accessed through the initial inelastic scattering of the tightly focused electron beam from the targeted atoms, we also identify favorable conditions for the electronically nonadiabatic evolution of the highly vibrationally excited system to open complex multistate reaction pathways. Implications of the early results for understanding the mechanisms and potential routes to improved efficiency and selectivity in beam-induced reactions are discussed. We conclude with a summary of the current state of theory and modeling capabilities in this area and provide our perspective on future directions for theoretical and experimental developments that we view as crucial to advancing the use of convergent electron beams in mode-specific, atomically precise platforms for direct-write materials modifications.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Changing Characteristics of Tropical Extreme Precipitation–Cloud Regimes in Warmer Climates

In this study, we investigated the changing characteristics of climatic scale (monthly) tropical extreme precipitation in warming climates using the Energy Exascale Earth System Model (E3SM). The results are from Atmospheric Model Intercomparison Project (AMIP)-type simulations driven by (a) a control experiment with the present-day sea surface temperature (SST) and CO 2 concentration, (b) P4K, the same as in (a) but with a uniform increase of 4K in the SST globally, and (c) the same as in (a), but with an imposed SST and CO 2 concentration from the outputs of the coupled E3SM forced by a 4xCO 2 concentration. We found that as the surface warmed under P4K and 4xCO 2 , both convective and stratiform rain increased. Importantly, there was an increasing fractional contribution of stratiform rain as a function of the precipitation intensity, with the most extreme but rare events occurring preferentially over land more than the ocean, and more so under 4xCO 2 than P4K. Extreme precipitation was facilitated by increased precipitation efficiency, reflecting accelerated rates of recycling of precipitation cloud water (both liquid and ice phases) in regions with colder anvil cloud tops. Changes in the vertical profiles of clouds, condensation heating, and vertical motions indicate increasing precipitation–cloud–circulation organization from the control and P4K to 4xCO 2 . The results suggest that large-scale ocean warming, that is, P4K, was the primary cause contributing to an organization structure resembling the well-known mesoscale convective system (MCS), with increased extreme precipitation on shorter (hourly to daily) time scales. Additional 4xCO 2 atmospheric radiative heating and dynamically consistent anomalous SST further amplified the MCS organization under P4K. Analyses of the surface moist static energy distribution show that increases in the surface moisture (temperature) under P4K and 4xCO 2 was the key driver leading to enhanced convective instability over tropical ocean (land). However, a fast and large increase in the land surface temperature and lack of available local moisture resulted in a strong reduction in the land surface relative humidity, reflecting severe drying and enhanced convective inhibition (CIN). It is argued that very extreme and rare “record-breaking” precipitation events found over land under P4K, and more so under 4xCO 2 , are likely due to the delayed onset of deep convection, that is, the longer the suppression of deep convection by CIN, the more severe the extreme precipitation when it eventually occurs, due to the release of a large amount of stored surplus convective available potential energy in the lower troposphere during prolonged CIN.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Vegetation transect surveys from 2022 and 2023 within the Kougarok Fire Complex (KFC) on the Seward Peninsula in Alaska, USA

This dataset contains recorded vegetation classifications from 16 transects within the Kougarok Fire Complex (KFC) area in Alaska. Each transect is 50 meters long, and vertical vegetation profiles were sampled every 50cm, resulting in 101 points recorded for each transect. At each point, a vertical rod was inserted in the ground and any living plant or substrate touching the rod was recorded from top to bottom. For each plant, the plant functional type (PFT) and species name was identified and recorded. The file "flmd.xlsx" provides file level meta data for all other included files. The file (transect_metadata.csv) contains descriptive information about each of the 16 transects (field observation year and date; latitude and longitude GPS measurements of the 0, 25, and 50 meter points; and whether that area burned during the 1971, 1997, 2002, 2015, and/or 2019 fires). The files (PFT_names.csv and species_names.csv) map the short hand labels for the plant function type (PFT) and species used in the transect surveys to their full names. The remaining files are named for each transect, and contain the results of the vegetation survey for each transect. These files contain the plant functional type and the species identified at each point. Points along the transect with more than one species at that point are listed vertically from top to bottom in these files.The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research.The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska.Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Water structure and electric fields at the interface of oil droplets

Interfacial water exhibits rich and complex behaviour, playing an important part in chemistry, biology, geology and engineering. However, there is still much debate on the fundamental properties of water at hydrophobic interfaces, such as orientational ordering, the concentration of hydronium and hydroxide, improper hydrogen bonds and the presence of large electric fields. This controversy arises from the challenges in measuring interfacial systems, even with the most advanced experimental techniques and theoretical approaches available. Here we report on an in-solution, interface-selective Raman spectroscopy method using multivariate curve resolution to probe hexadecane-in-water emulsions, aided by a monomer-field theoretical model for Raman spectroscopy. Our results indicate that oil-water emulsion interfaces can exhibit reduced tetrahedral order and weaker hydrogen bonding, along with a substantial population of free hydroxyl groups that experience about 95 cm -1 redshift in their stretching mode compared with planar oil-water interfaces. Given the known electrostatic zeta potential characteristic of oil droplets, we propose the existence of a strong electric field (about 50-90 MV cm -1 ) emanating from the oil phase. This field is inferred indirectly but supported by control experiments and theoretical estimates. These observations are either absent or opposite in the molecular hydrophobic interface formed by small solutes or at planar oil-water interfaces. Instead, water structural disorder and enhanced electric fields emerge as unique features of the mesoscale interface in oil-water emulsions, potentially contributing to the accelerated chemical reactivity observed at hydrophobic-water interfaces.

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

Fostering a Guiding Multiscale Model for the Development of Advanced MgB 2 Hydrogen Storage Materials (Final Technical Report)

Project Goal and Objective. The demand for energy and for an upgraded energy infrastructure has steadily grown, as have the needs for energy independence and alternatives to our reliance on petroleum. Hydrogen is considered the most viable fuels for wide-scale implementation in the near future as it is less-polluting, non-toxic, and has more stored energy than petroleum. It is envisioned that hydrogen can eventually become the prime energy carrier, integrating the transportation, grid, and chemical sectors in a way that improves resiliency, diversifies feedstocks, and affords new economic opportunities. A key remaining challenge is the development materials with enhanced gravimetric and volumetric hydrogen storage capacities that offer a higher performance than compressed gas. These materials would eliminate the need for large-scale compression, thereby dramatically reducing the footprint and cost of gas storage. The high gravimetric and volumetric hydrogen capacities of complex hydrides has prompted an intensive investigation of the potential of this class of materials as hydrogen storage media over the past 25 years. Among the many complex hydrides that have been explored, magnesium borohydride, Mg(BH 4 ) 2 , has been found to possess the best combination of practical thermodynamic properties. These include a gravimetric H 2 density of 14.9 wt% H 2 and thermodynamics for the dehydrogenation of Mg(BH 4 ) 2 to MgB 2 (equation 1) (ΔH° = 39 kJ/mol H 2 , ΔS = 112 J/K mol H 2 ) which lie in the narrow window required Mg(BH 4 ) 2 $\Leftrightarrow$ MgB 2 + 4 H 2 (1) for reversibility under moderate pressure and temperature. However, overcoming the extremely slow kinetics of the reversible release of hydrogen by this material in the solid state is a daunting challenge. At temperatures greater than 400 °C, the borohydride releases up to 14 wt% hydrogen giving MgB 2 . We discovered that the direct re-hydrogenation of MgB 2 to Mg(BH 4 ) 2 can be accomplished under 950 bar H 2 at 400 °C. While this demonstrated that complete reversibility can be achieved, the conditions employed are far too extreme for commercial hydrogen storage applications. More recently, we found through US DOE funded research projects (EERE HyMARC and HySCOR), that hydrogen cycling, can be accomplish at much milder conditions upon modification of the borohydride or boride. Guided by these discoveries these discoveries, the objective of this research project was to obtain key information that will enable the development of a model of reversible hydrogenation of MgB 2 to Mg(BH 4 ) 2 . The ultimate goal of our efforts is to attain a model of this transformation that can be utilized to accelerate development further advanced materials. This project directly follows on discoveries that were made over the course of a US DOE, EERE HyMARC project that was focused on improvement of the hydrogen cycling kinetics of modified MgB 2 . We found that that mechanical milling with graphene results the desired, pronounced kinetic enhancement. The dramatic lowering of the conditions required for the hydrogenation of MgB 2 is a significant step towards overcoming its chemical inertness allowing its development as a practical onboard hydrogen storage material. However, the exact nature of the modification(s) of MgB 2 that is responsible for its activation towards hydrogenation is completely unknown. This situation is not unique, as efforts to develop hydrogen storage materials typically have a narrow focus rather than a comprehensive approach that takes atomic level bonding and structure; molecular dynamics; long range, nano- and mesoscale-structure and their interconnection all into account. The goal of this project was the development of a comprehensive, multi-scale computational model of reversible hydrogenation of MgB 2 to Mg(BH 4 ) 2 that can be utilized for development of higher performance versions of the modified material. Development of the model requires determination of: 1) the bulk, nano-scale, and meso-scale structural changes occurring at elevated pressure following mechano-chemical modification of MgB 2 ; 2) the reaction pathway of the reversible hydrogenation of MgB 2 to Mg(BH 4 ) 2 ; 3) the effect of elevated pressure and mechano-chemical modification on the chemical reaction pathways; 4) the interactions at solid-gas interfaces; and particle surfaces; and 5) the kinetics and thermodynamic parameters associated with each step of the hydrogenation reaction pathway. This investigation required advanced techniques as preliminary, standard XRD, 11 B NMR, and FTIR analysis showed no signs of material modification. In order to gain this level of understanding of modified MgB 2 , required the teaming of a diverse group of experts and state-of-the art experimental capabilities at the University of Hawaii at Manoa (UHM) and collaborating National Laboratories: Craig Jensen , Department of Chemistry (PI and Project Director), solid state, solution, and high pressure NMR spectroscopy; solid-state synthesis; and high pressure hydrogenation (collaboration with SNL); Godwin Severa , Hawaii Natural Energy Institute (co-PI) calorimetry; infrared and Raman spectroscopy (collaboration with NREL); Dera , high pressure X-ray diffraction including in situ experiments (collaboration with ANL); Hope Ishii , Hawaii Institute of Geophysics electron microscopy investigations (collaboration with LBNL); and Joe Brown , Mechanical Engineering , material electronic structure and electric field effects.

08 HYDROGEN↗
DOE OSTI2022

Atomic-Scale Surface Studies of Bulk Metallic Glasses. Final Report

Bulk metallic glasses (BMGs) are of both scientific and technological interest because the absence of periodic atomic arrangements provides them with unique physical, chemical, and mechanical properties. Their high strength, superior elasticity, and an ability to be easily formed into virtually unlimited shapes with feature sizes from centimeters to Angstroms by blow molding and thermoplastic forming makes them an attractive choice for more and more practical applications and products. Due to their complex internal structure, however, experiments that yield insight into their exact atomic arrangements have been scarce. As a result, glass physics is one of the last remaining unexplored fields of materials science despite the scientific and technological importance of glasses in our daily lives, and the question how to characterize and control matter away from equilibrium, as glasses are, was listed as one of five Grand Challenges in a recent DoE report. The main reason for the slow progress in glass physics is the lack of experimental tools that enable access to atomic-scale structural and behavioral information for disordered materials. Such atomic-scale knowledge is mandatory to establish structure-property relationships that ultimately could allow to custom-design alloys featuring specific desired characteristics. With no such relationships available, theory development in glass remains basic and the few that exist are often untested. The aim of this research was to enable progress in our understanding of BMGs by developing a new approach that will allow a meaningful application of local surface science methods to specially prepared BMG samples to obtain a wealth of quantitative information on their atomic arrangements. Key was the availability of specially prepared samples whose surfaces feature large atomically flat terraces despite being entirely amorphous, which we have produced from a Pt 57.5 Cu 14.7 Ni 5.3 P 22.5 alloy (‘Pt-BMG’) both under ambient conditions as well as in ultrahigh vacuum using a unique setup that has been specially developed within this grant. Our approach starts with the in-situ preparation of oxide crystals that are terminated by large terraces, from which exact mirror images out of BMG will be produced using thermoplastic forming (TPF); for the research within this grant, we have successfully used (001)-oriented SrTiO 3 single crystal surfaces as well as (100)-, (110)-, and (111)-oriented single crystals made from LaAlO 3 . Since the resulting BMG replicas display all features of the original crystal with sub-Angstrom fidelity, thereby mimicking the original crystal’s termination by atomically flat terraces without being crystalline themselves, they are ideally suited for further investigation. The following atomic-scale local studies were then carried within this proposal: (i) high-resolution surface imaging and local spectroscopy using scanning probe microscopy, which showed disordered atom-like features and revealed changes in the gradient of the local surface potential on a 1-2 nm length scale; (ii) characterization of atomic-scale plastic flow with affected volumes as low as 1000 atoms, which showed local hardness near or above the theoretically predicted maximum and, once plastic deformation was initiated, homogeneous flow of the atoms involved; (iii) characterization of surface relaxation processes and the onset of crystallization induced by annealing, which showed that upon heating over the material’s glass transition temperature, the surface rearranges and relaxes towards a more stable, denser packed glass, which increases on-terrace surface roughness, while surface tension smoothens step edges; and (iv) studies that investigate the dependence of the material’s mechanical properties and structure on processing parameters, revealing that relaxed glasses get denser, harder, and more elastic. In combination, this information allows to combine structural models with mechanical properties and preparation history, thereby facilitating the development of preparation-structure-property relationships for metallic glasses. With the availability of such information, bulk metallic glasses can be further optimized to be used in more and more applications in industry.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Pattern-enhanced Resonant Soft X-ray Scattering for Operando monitoring of electrochemical solid-liquid interfaces

Unveiling interfaces at sub-nanometer scales is essential for advancing the understanding of complex chemical transformations. However, characterizing solid-liquid interfaces with high dimensional sensitivity and temporal resolution remains challenging, due to their dynamic nature and inaccessibility by conventional probes. Here we present an approach, Pattern-enhanced Resonant Soft X-ray Scattering, to overcome the challenges. Rooted in a “sample-as-optics” philosophy, this technique utilizes precisely engineered line-grating nanopatterns to modulate near-field X-ray illumination, coherently enhancing scattering signals from the line-gratings. We implement the method using Ni line-grating nanopatterns in electrochemical water oxidation. The periodic nanostructures serve as diffractive optical elements to reveal the Ni oxidation gradients and structural dynamics at the electrode-electrolyte interfaces. Finite-element simulations corroborate the observed trends by modeling variations in compositions and structures during electrocatalysis. Through integrating advanced sample design with coherent wave nature of soft X-rays, our approach opens accessible pathways to operando exploring chemical evolution and sub-nanometer dimensional variations simultaneously in electrochemical systems. This non-destructive method is efficient and element-specific, making it valuable for probing chemical and dimensional dynamics with appropriate modeling.

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

MPEX AI Digital Twins

All magnetically confined plasma fusion power plant concepts (Tokamak, Spherical Tokamak, Stellarator, Mirror, ...) must exhaust the heat and plasma from the core confinement region to the material walls. The primary channel for this exhaust is through a plasma divertor which directs plasma along open magnetic field lines to a material target. The Material Plasma Exposure eXperiment (MPEX) illustrated in Figure 1, is a high-power, steady-state linear plasma device designed to produce the plasma material interaction (PMI) conditions of the divertor of future magnetic confinement fusion power plants: energy flux 20MW/m 2 , ion fluence 1031/m 2 , pulse duration 106 sec. These goals of plasma exposure in MPEX are well beyond those achieved in magnetic fusion experimental devices. Successfully achieving these high power steady state conditions for long pulses requires operational control of the heating and particle sources and the plasma flux to the walls and target. The MPEX AI Hot Spot Controller, proposed in this project, will help achieve the operational milestones of MPEX. The MPEX device will begin commissioning at the end of FY26. A smaller proto-MPEX was operated for 14,666 plasma discharges and will resume operation in September of 2025 as proto-MPEX-lite, with reduced capability, to test a new window for the Helicon plasma source. The proto-MPEX data has undergone surrogate modeling with machine learning methods (R. Archibald, 2022 IEEE International Conference on Big Data). This proto-MPEX data will be used to begin development of the AI digital twins described in this white paper. The scientific mission of MPEX is to qualify materials of different composition for use in the high energy and plasma flux conditions of a fusion power plant. The materials exposed in MPEX will in some cases be exposed to high neutron fluxes at other ORNL facilities to measure the changes to their PMI properties. The targets exposed in MPEX will be transported under vacuum to a Surface Analysis Station (SAS). The SAS will be equipped with the following diagnostics: Focused Ion Beam (FIB) for trench milling, 100-400 angstrom resolution scanning electron microscope (SEM), surface mapping x-ray spectrometer, high resolution camera, and a future upgrade to a laser induced breakdown spectroscopy quadruple mass spectrometer (LIBS-QMS). The MPEX experiments will generate diverse pre- and post-exposure measurement data of detailed material properties down to the crystal grain level in 3D for post-exposure assessment of PMI damage (e.g. cracking, melting, erosion and redeposition of the material). Physics models for the PMI, and how the material composition and manufacturing impact its performance under high energy plasma exposure, need to be validated with MPEX data to guide the selection of new candidate materials. Our vision for the MPEX AI Digital Twins project is to supply experimental and physics model simulation data to train Artificial Intelligence (AI) models for data processing, analysis, operational control, PMI and materials simulation to maximize the scientific output of the MPEX device. Ultimately, an AI digital twin of MPEX material assessment metrics for tested and synthetic material types with simulated PMI will be trained by the AI Modeling Teams on the experimental and physics simulation data submitted to the American Science Cloud by this project. A purely empirical search for the best material is inefficient given the finite number of samples that can be tested on MPEX. In order to expand the material properties database for training the MPEX Material Assessment AI Digital Twin, and to gain physics understanding of the PMI processes, physics models of the material properties and PMI processes are required. The physics simulations provide detailed simulation data, like impact angles for plasma ions, sputtering yields, transport of the ionized sputtered target material in the plasma, and redeposition locations. This simulation data expands the measurement data for deeper physics understanding. The experimental data is essential to validate the PMI and material structure simulation models. The validated models can then be used to generate new simulation data of MPEX material assessments for synthetic material compositions that have not been exposed in MPEX. These predictive simulations, plus the whole experimental dataset, will be used to train the MPEX Material Assessment AI Digital Twin allowing a rapid generative AI search for new materials with reduced PMI damage by interpolating the domain of the training set. These new optimum materials can be simulated with the physics codes and/or tested in MPEX. The ability of AI neural networks to interpolate multi-dimensional parameter spaces and generate virtual data is exploited for a more efficient search for optimum materials. The advent of the Transformational AI Models Consortium (TAIMC) is an opportunity to engage with state of the art private and public AI developers to achieve the goals of the AI digital twins and AI accelerated physics models proposed in this project. Our partners at ORNL from the Advance Scientific Computing Research (ASCR) organization will collaborate in accelerating the integrated plasma material interaction simulation framework. This simulation framework will provide a platform for generating simulation data across a range of physical fidelities, including hybrid methods that produce multi-fidelity results. This data will be leveraged for AI model development, both for generation of surrogates and the automation of simulation campaigns. A part of the research below will include collaborative efforts with the TAIMC to (i) adapt data storage approaches to ensure AI-readiness, (ii) provide a protypical exemplar to inform and exercise constructed workflows, and (iii) generate and share data, using the TAIMC unified AI data standard, for foundational models that will be trained from multiple sources across the DOE complex. We will also collaborate with the TAIMC, as well as the planned AI modeling teams, to develop approaches for reducing the cost of data generation. These include tailored multi-fidelity approaches as well as fine-tuning strategies to augment general, large-scale foundational models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

VENTSAR V3.0: A Python GUI for Estimating Contaminant Concentrations on or Near Buildings Due to Building Effects and Plume Rise and For Calculating Inhalation and Plume Shine Doses

VENTSAR: Originated as VENTAX (Smith and Weber 1983) on the IBM Mainframe at the Savannah River Site (SRS) as a Fortran program. Updated to VENTSAR XL V1.0 (Simpkins 1997) as a spreadsheet version using macros created in Microsoft Excel. Later updated to VENTSAR XL V2.0 (Dixon 2018) as a spreadsheet executing a modified version of the original VENTAX Fortran code. Estimates contaminant concentrations on or near a building from a release at a nearby location. Calculates concentrations for a given meteorological exceedance probability or for a given stability and wind speed combination. Can model a single building with or without a penthouse on top or a ground location from either a stack or ground release. Plume rise can be considered. Contaminant releases can be chemical or radioactive with downwind concentrations determined at user-specified distances. Wind passing over and around buildings creates a complicated dispersion pattern. Air-intake vents may be located on building roofs or near the ground downwind of a release source. An estimation of pollutant concentrations on or near a structure is important in determining expected pollutant levels. Meteorological data are selected based on a specific area of the site. Fortran-based VENTAX able to make fast, complex mathematical calculations. Not user-friendly VENTSAR XL V1.0 no longer supported due to obsolete Excel macros. VENTSAR XL V2.0 an interim solution using an Excel spreadsheet to interface with the modified VENTAX Fortran code. Requires user to have Microsoft Excel. Not an intuitive interface for someone unfamiliar with VENTSAR. Does not perform dose calculations. VENTSAR V3.0 goal to combine the strengths of previous versions of VENTSAR into a single, powerful yet user-friendly program with a Graphical User Interface (GUI). Not dependent upon a specific program or plug-ins. Self-contained package to be used on any computer running Microsoft Windows. User-friendly, intuitive GUI created in Python for parameter input. Executes same modified VENTAX Fortran code as VENTSAR XL V2.0. Outputs formatted text file of completed calculations for easy review and dissemination. Performs inhalation and plume shine dose calculations for up to 11 user-selected nuclides from a nuclide dose factor library of almost 500 nuclides. 12 test cases were created for verification of V1.0 and V2.0. Same test cases run in V3.0 to verify correct performance. V3.0 produced similar results to V1.0. Confirmed GUI did not alter calculations in any way. Comparisons of the test case results confirm the V3.0 GUI does not influence the VENTSAR calculations. Simply passes same input parameters to VENTAX Fortran code for execution. Provides end user an easy-to-use, intuitive tool to quickly make building effect and plume rise calculations, as well as, inhalation and plume shine dose calculations. No experience with or working knowledge of Fortran, Excel, command line, or macros required. Self-contained package allows VENTSAR be deployed to any Windows computer without the need for specific programs or plug-ins.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

Reducing Long‐Standing Surface Ozone Overestimation in Earth System Modeling by High‐Resolution Simulation and Dry Deposition Improvement

The overestimation of surface ozone concentration in low‐resolution global atmospheric chemistry and climate models has been a long‐standing issue. We first update the ozone dry deposition scheme in both high‐ (0.25°) and low‐resolution (1°) Community Earth System Model (CESM) version 1.3 runs, by adding the effects of leaf area index and correcting the sunlit and shaded fractions of stomatal resistances. With this update, 5‐year‐long summer simulations (2015–2019) using the low‐resolution CESM still exhibit substantial ozone overestimation (by 6.0–16.2 ppbv) over the U.S., Europe, eastern China, and ozone pollution hotspots. The ozone dry deposition scheme is further improved by adjusting the leaf cuticle conductance, reducing the mean ozone bias by 19%, and increasing the model resolution further reduces the ozone overestimation by 43%. We elucidate the mechanism by which model grid spacing influences simulated ozone, revealing distinctive pathways in urban versus rural areas. In rural areas, grid spacing mainly affects daytime ozone levels, where additional NO x emissions from nearby urban areas result in an ozone boost and overestimation in low‐resolution simulations. In contrast, over urban areas, daytime ozone overestimation follows a similar mechanism due to the influence of volatile organic compounds from surrounding rural areas. However, nighttime ozone overestimation is closely linked to weakened NO titration owing to the redistribution of urban NO x to rural areas. Additionally, stratosphere‐troposphere exchange may also contribute to reducing ozone bias in high‐resolution simulations, warranting further investigation. This optimized high‐resolution CESM may enhance understanding of ozone formation mechanisms, sources, and changes in a warming climate.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Soil Temperature and Moisture within the Kougarok Fire Complex, Kougarok Road Mile Marker 86, Seward Peninsula, Alaska, 2019-2023

Daily averages of soil temperature and moisture measured once every hour at different heights located at Intensive Monitoring Stations within the Kougarok Fire Complex, Kougarok Road Mile Marker 86 site. Data were retrieved annually from 2019-2023. Package contains 21 *.CSV data files plus a file level metadata *.CSV, data dictionary *.CSV, data file inventory *.CSV, and sensor location site map *.JPG. Data files have header rows, NaN fields indicate invalid or missing data, and negative vertical offsets are above ground.The Kougarok tundra fire complex (KFC) is located north of Nome and the Kigluaik Mountains, near Quartz Creek and the Kougarok River. The site is accessed by foot from the end of the Nome-Taylor Highway (mile marker 86; also called the Kougarok or Beam Road). The KFC burned in six major fires in the decades since 1950 (Alaska Interagency Coordination Center, unpublished data). Lightning ignited five of these fires (1971, 1997, 2015, and 2019) and one was human caused (2002). The mosaic of overlapping fire scars allows for the study of repeat fires in the tundra which, until recently, was not a common phenomenon outside the boreal forest in Alaska. Our reference unburned tundra fire site is south of the KFC located at mile marker 80 of the Nome-Taylor Highway.The two most recent fires are the Mingvk Lake (2015; 21,698 acres burned from 7/27/2015 to 9/28/2015) and Garfield Creek (2019; 422 acres burned from 7/31/2019 to 8/20/19). The Mingvk Lake fire scar includes areas that burned 1-4x (1971, 1997, 2002), while the entirety of the Garfield Creek fire scar has burned 2x previously (1971, 2002).Previous research at the KFC focused on permafrost (Liljedahl et al. 2007; Narita et al. 2015; Iwahana et al. 2016; Tsuyuzaki, Iwahana, and Saito 2017) and vegetation (Narita et al. 2015; Hollingsworth et al. 2021) response to fire. The central Seward Peninsula is characterized by continuous permafrost with a thickness of 15 to 30 m and a mean active layer thickness of 56 cm (Hinzman et al. 2003). Sloping hills with mixed shrub–tussock tundra and tussock tundra vegetation in the uplands are characteristic of the region. Three micrometeorological towers near the Kougarok field site recorded a mean annual temperature of −2.4°C, mean January temperature of −23.1°C, mean July temperature of +11°C, and mean summer rainfall (June–August) of 94 mm from 2000 to 2006 (Liljedahl et al. 2007).The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research.The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska.Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Monthly averages of ED2 model simulations initialized with airborne lidar structure, Jan 1981-Dec 2018, Brazilian Amazon

Deforestation and forest degradation (selective logging, fires, fragmentation) have impacted nearly 40% of the original extent of the Brazilian Amazon, and have markedly impacted forest structure across the region. To date, few studies analysed how shifts in forest structure from degradation influence the forest sensitivity to climate extremes, because of the complex interactions between forest structure and micro-environmental conditions. To address this knowledge gap, we carried out a series of simulations across the Brazilian Amazon using the Ecosystem Demography Model (ED2), using observed forest structure derived from 541 airborne lidar transects (375 ha each) and two scenarios representing forest recovery and expansion of degradation to investigate how shifts in forest structure impact ecosystem function under near-average and extreme climate conditions, as part of the manuscript Longo et al 2025 "Degradation and Deforestation Increase the Sensitivity of the Amazon Forest to Climate Extremes". This dataset provides the output results from the ED2 model simulations for the three simulations at monthly time scales, in NetCDF format. For all simulations, we used bias-corrected hourly reanalyses (WFDE5) for most meteorological drivers, except for precipitation, which was obtained from CHIRPS. The meteorological drivers used in the study span 38 years (Jan 1981–Dec 2018). The output results correspond to the last 38 years of simulation (one full cycle of meteorological drivers), in which ED2 simulations used static stand structure (i.e., the forest structure was held constant). The following files are provided:ED2_emean_Global_R004_BrAmaz_s1c0t0l0f0.nc. This corresponds to the Control simulation. The forest structure was obtained from the airborne lidar.ED2_emean_Global_R005_BrAmaz_s1c0t1l1f0.nc. This corresponds to the Degraded simulation. The forest structure was obtained from a spin-up simulation initialized with airborne lidar and a scenario that expanded deforestation and selective logging across the Amazon.ED2_emean_Global_R006_BrAmaz_s1c0t1l0f0.nc. This corresponds to the Recovery simulation. The forest structure was obtained from a spin-up simulation initialized with airborne lidar and a scenario that completely halted deforestation and degradation, allowing degraded forests to recover for 38 years.We also provide file ED2_zones_R004_BrAmaz_s1c0t0l0f0.nc, which classifies each grid cell into zones used in the reference manuscript: 1: Southeast. 2: South. 3: West. 4: Central. 5: Northeast. 6: North. 7: Northwest". Index 0 corresponds to grid cells excluded from sub-region analyses because they were dominated by flooded forests, deforestation, and naturally non-forest vegetation.

54 ENVIRONMENTAL SCIENCES↗