Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “surface conditioning”

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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Toward Hybrid Physics-Machine Learning to improve Land Surface Model predictions

A critical challenge for Land Surface Models (LSMs) is to simulate processes at the surface and the subsurface and their feedbacks to the atmosphere. Even using the same climate forcings, different LSMs predict different surface fluxes and soil moisture conditions due to differences in the formulations of individual processes, parameterizations, and representation of spatial heterogeneity. Ultimately, these differences contribute to the LSM prediction errors and uncertainty. This research seeks to address this challenge by coupling physics-based modeling with state-of-the-art machine learning (ML) techniques to describe complex physical and biogeochemical processes and narrow the gap between model predictions and observations.

58 GEOSCIENCES↗
DOE OSTI2021

SR19021 Tritium Aging of Regenerated LANA.75 (Final Project Report)

The Savannah River Tritium Enterprise (SRTE) has used the metal hydride LaNi 4.25 Al 0.75 (LANA.75) in the Tritium Facilities for over two decades. LANA.75 beds store significant quantities of tritium but have a limited service life due to the radiolytic decay of tritium to He-3 within the metal matrix. It has been shown that heating tritium-aged LANA.75 under vacuum can reverse tritium aging effects, eliminating the heel of trapped hydrogen, and restoring the reversible capacity. Additional investigation is needed to ensure there are no unexpected changes to the hydride before this restoration technique is employed in full scale beds in the Tritium Facilities. This project was to be comprised of three distinct scopes: obtain tritium aging data on the sample regenerated in 2018, regenerate a tritium-aged sample at 600 °C, and perform thermal stability testing on a non-tritiated sample. Isotherms were collected on the previously regenerated sample after approximately two years of tritium aging. Isotherms were collected at 80, 100, and 120 °C. As expected, there was a decrease in the plateau pressure, an increase in plateau slope, and a portion of the “heel” of tritium trapped in the metal had been reestablished. Unexpectedly, it appeared that the plateau had shortened at the higher tritium to metal ratios as well. This is typically seen in older samples. The second scope, to regenerate a second tritium-aged LANA.75 sample, was not completed. A Task Technical and Quality Assurance Plan was written and approved, a high temperature test cell was fabricated, several pre-job briefs were held, and the hydride sample was passivated with air. Despite these successes, the hydride sample was not recovered from the legacy test cell. The third scope was to perform thermal stability testing on a “cold” LANA.75 sample. A non-tritiated sample of LANA.75 was held at 750 °C under vacuum for 200 hours to simulate exposure to multiple regeneration evolutions. Hydride isotherm performance, chemical composition, crystallinity, particle size, and morphology are compared between the pre- and post-regeneration samples. No significant changes were observed in composition, crystallinity, or particle size. Comparison of before and after isotherms showed that performance improved rather than deteriorated during the evolution. Scanning Electron Microscopy (SEM) analysis showed small growths on the particle surface after exposure to regeneration conditions. Additional testing will be required to determine the cause of these growths.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2022

Fundamental Studies on the Recovery of Rare Earth Elements from Coal and Coal Byproducts (Final Product Report)

The objective of this project was to collect fundamental information that can be used to develop disruptive physical and chemical separation methods that may be used to extract the rare earth elements (REEs) from coal-based feedstocks in a highly efficient, cost-effective, and environmentally benign manner. To meet the stated objective, a set of interrelated fundamental studies have been conducted. The project has been subdivided into three broad areas of tasks as follows, (1) measurement of surface forces in the thin liquid films (TLFs) of water confined between two macroscopic surfaces to establish the optimal conditions for efficiently collecting the ultrafine particles of rare earth minerals (REMs) present in the coal-based feedstocks using the flotation and/or the hydrophobic-hydrophilic separation process(es), (2) identification of the fundamental mechanisms by which REE ions adsorb and desorb on clay minerals present in coal byproducts, and the characterization of the ion-adsorption clays (IACs) present in U.S. coals, and (3) developing methods of extracting the REEs from the REMs and IACs present in US coal byproducts in a cost-effective and sustainable manner.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2021

Uranium–Molybdenum Alloy Critical Experiments for the Design of the Health Physics Research Reactor

Clean critical experiments with a uranium-molybdenum alloy (average of 10.1616 wt. % Mo with a density of 17.08 g/cm 3 ) were performed at the Oak Ridge Critical Experiments Facility in 1961 to support the design of the Health Physics Research Reactor (HPRR). The HPRR was similar to the Godiva burst reactor at Los Alamos National Laboratory and was designed to produce 50 microseconds burst of 10 17 fission pulses of radiation for dosimetry measurements, initially in support of the determination of the doses from the nuclear detonations in Japan during World War II. These experiments reported here were used to verify the calculational methods used to design the HPRR. These delayed critical measurements were:1) a solid unreflected and unmoderated 8-in.-dimeter U-Mo cylinder, 2) an unmoderated and unreflected annulus with 8-in.-outside diameter, 2-in.-inside diameter cylinder with a central void, 3) an unmoderated and unreflected annulus with 8-in.-outside diameter, 2-in.-inside diameter cylinder with a central void filled with stainless steel, 4) Same as 3) but with 3-in-thick Plexiglas reflector on top with and without cadmium between the reflector and the U-Mo alloy assembly with steel in the center, and 5) an unmoderated and unreflected annulus which was a modification of the second but with the lower 5 inches of the central hole enlarged to 3.5 in. with various reflector conditions. The reflector conditions were: 1-in.-thick Plexiglas on all outer surfaces-void in the center; 1-in.-thick Plexiglas on all outer surfaces-Plexiglas in the center; 2-in.- thick Plexiglas on radial surface-void in the center; 6-in.-thick Plexiglas on the bottom only-Plexiglas in the center; and 6-in.-thick Plexiglas on bottom, 1-in.-thick on top and on the lower 8.25-cm.-section of the radial surface-void in the center. For some of these reflector conditions 0.025-cm.thick cadmium was located between the reflector and the U-Mo alloy. The uranium contained 93.17 wt. % 235 U. Reflection was a safety concern for this unmoderated and unreflected reactor and reduction of reflection effects was also investigated by insertion of neutron absorber around the U-Mo alloy. The stainless steel 304 contained 18% nickel and 8% chromium and the rest iron. The reflector material was a methacrylate plastic (Plexiglas) containing 5.8 x 10 22 atoms/cm 3 of hydrogen and 3.6 x 10 22 atoms/cm 3 of carbon with a density of 1.20 g/cm 3 . The purpose of this report is to document the experimental information for the measurements performed so that at a later date researchers could perform the required uncertainty and calculational analyses and documentation to use these data for an International Nuclear Criticality Safety Benchmark Program (ICSBEP) or a EURATON Nuclear Energy Agency (NEA) benchmark. The data from the experiments described should be acceptable for use as criticality safety benchmark experiments for the ICSBEP and the NEA nuclear criticality safety benchmark program, once the uncertainty analysis is completed. Based on previous ICSBEP benchmarks with this enriched uranium metal at ORCEF, the uncertainties in k eff could be as low as ±0.0002 for some configurations.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

Crush Gun Impact of Unconfined Explosive with Flat-Faced Projectiles at Nominal Velocity of 60 m/s

A pneumatic gun apparatus was designed and built to deliver crushing impacts at velocities of 0-70 m/s. The apparatus is designed to halt the projectile (if desired) approximately 1 mm before impacting the sapphire imaging surface, preventing a complete “pinch” condition from attaining. The primary diagnostics include high-speed visible videography and infrared imaging to determine spatial temperature rise of the sample during impact.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2022

Capillary Pressure to Relative Permeability: Task 7.3.4 Relative Permeability Analysis Final Report

Surface tension affects all aspects of fluid flow in porous media. Through measurements of surface tension interaction under multiphase conditions, a relative permeability relationship can be determined. Relative permeability is a numerical description of the interplay between two or more fluids and the porous media they flow through. It is a critical parameter for various tools used to characterized subsurface multiphase flow systems, such as numerical simulation for oil and gas development, carbon sequestration, and groundwater contamination remediation. Therefore, it is critical to get a good statistic distribution of relative permeability in the porous media under study. Empirical relationships for determining relative permeability from capillary pressure are already well established but do not provide the needed flexibility required to match laboratory derive relative permeability relationships. By expanding the existing methods for calculating relative permeability from capillary pressure data, it is possible to create both two and three-phase relative permeability relationships. Existing laboratory-measured relative permeability data along with mercury intrusion capillary (MICP) data coupled with interfacial tension and contact angle measurements were used to determine the efficacy of this approach to relative permeability curve creation. The relative permeability relationships determined with this method were fit to existing laboratory data to elucidate common fitting parameters that were then used to create relative permeability relationships from MICP data that does not have an associated laboratory-measured relative permeability relationship.

58 GEOSCIENCES↗
DOE OSTI2026

Floating Solar at Lima's Twin Lakes Reservoir

The City of Lima, Ohio is developing a floating photovoltaic (FPV) facility on the surface of Twin Lake Reservoir, adjacent to the City’s water treatment plant. The FPV array is designed with a capacity of 2 megawatts (MW) and will consist of approximately 3,500 solar panels spanning more than four acres of water surface. Based on regional solar conditions, the system will generate 2,000 megawatt-hours (MWh) of electricity annually. While the FPV system will interconnect to the electrical infrastructure serving the water treatment plant, power produced by the system will be net-metered through the utility grid, with energy credits applied to offset the plant’s overall electrical demand rather than directly supplying all plant loads in real time.

14 SOLAR ENERGY↗
DOE OSTI2020

Tritium Transport Phenomena in Molten-Salt Reactors

In this work, we review phenomena relevant to tritium transport in molten-salt reactors, which produce tritium from lithium and beryllium salts at significantly higher levels than other reactor types. Modeling of such phenomena began following MSRE operations, and these early models attempted to predict measured tritium distributions in the MSRE, accounting for turbulent mass-transport processes (using established heat transfer correlations), permeation through a variety of metal structures such as heat exchanger tubes, and transport to and from bubbles introduced into the salt by gas sparging. The models reasonably reproduced the MSRE data, but did so best when the permeability of structures was reduced by about a factor of 1,000. This issue does not appear to have been conclusively resolved, and all of the more recent attempts to model tritium transport in molten salts appear to make use of the same methodology. MSRE remains, however, essentially our only source of integral tritium transport data relevant to MSRs. Here, we generalize the MSRE approach to permeation in order to include potential rate-limiting effects at interfaces, as well the effects of added hydrogen. Appropriately non-dimensionalized, this system of equations identifies two dimensionless numbers whose relative values clearly delineate the conditions under which mass transport, surface effects, permeation, and hydrogen swamping are expected to become rate-limiting. We also describe the preliminary conceptual design of a forced convection FLiBe loop, into which tritium would be introduced for the purpose of providing validation data for such a model. The primary purpose of this is to investigate the coupled transport phenomena described above and identify those that are rate-limiting in MSRs. Additional test section configurations are described that would address other transport phenomena relevant to MSRs, including bubbly flows and graphite interactions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Machine Learning for Improving Surface-Layer-Flux Estimates

Abstract Flows in the atmospheric boundary layer are turbulent, characterized by a large Reynolds number, the existence of a roughness sublayer and the absence of a well-defined viscous layer. Exchanges with the surface are therefore dominated by turbulent fluxes. In numerical models for atmospheric flows, turbulent fluxes must be specified at the surface; however, surface fluxes are not known a priori and therefore must be parametrized. Atmospheric flow models, including global circulation, limited area models, and large-eddy simulation, employ Monin–Obukhov similarity theory (MOST) to parametrize surface fluxes. The MOST approach is a semi-empirical formulation that accounts for atmospheric stability effects through universal stability functions. The stability functions are determined based on limited observations using simple regression as a function of the non-dimensional stability parameter representing a ratio of distance from the surface and the Obukhov length scale (Obukhov in Trudy Inst Theor Geofiz AN SSSR 1:95–115, 1946), $$z/L$$ z / L . However, simple regression cannot capture the relationship between governing parameters and surface-layer structure under the wide range of conditions to which MOST is commonly applied. We therefore develop, train, and test two machine-learning models, an artificial neural network (ANN) and random forest (RF), to estimate surface fluxes of momentum, sensible heat, and moisture based on surface and near-surface observations. To train and test these machine-learning algorithms, we use several years of observations from the Cabauw mast in the Netherlands and from the National Oceanic and Atmospheric Administration’s Field Research Division tower in Idaho. The RF and ANN models outperform MOST. Even when we train the RF and ANN on one set of data and apply them to the second set, they provide more accurate estimates of all of the fluxes compared to MOST. Estimates of sensible heat and moisture fluxes are significantly improved, and model interpretability techniques highlight the logical physical relationships we expect in surface-layer processes.

Meteorology & Atmospheric Sciences↗
DOE OSTI2026

Surface and Buried Thermal, and RGB Unexploded Ordnance Data Collection

This document provides a description of a data collection campaign of unexploded ordnance (UXOI) set. The dataset captures a controlled UAV imaging campaign designed to support detection of UXO across varied environmental conditions. Data were collected during three campaigns in Norris and Northeast Knoxville, Tennessee, using RGB, and thermal sensors mounted on Parrot UKR. In total, the dataset contains 9925 images, 26 full-motion video, and approximately 81.99 GB of data, collected across late spring/summer conditions, every hour during sunlight, and multiple surface contexts, including tall grass, short grass, gravel, as well as buried in sand, and other gravel mixtures. The collection was designed to capture thermal and visual variability relevant to UXO detection in agricultural land, bare earth, and subsurface. Review of the imagery showed that ordnance was most detectable during periods of changing solar input, especially approximately 10-60 minutes after sunrise, approximately 20-60 minutes after sunset, and 2-3 min after cloud cover interrupted prolonged solar heating. These conditions increased thermal contrast because many ordnance items retained or released heat differently than the surrounding vegetation and ground surface. This dataset provides a useful resource for developing and evaluating airborne UXO detection methods under realistic field conditions. All ordnance used in the study was inert, and thermal behavior may differ from that of live ordnance. In addition, variation in ordnance type, composition, and placement introduced differences in thermal response that should be considered when interpreting results.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2022

Mechanistic Study of 1,2-Dichloroethane Hydrodechlorination on Cu-Rich Pt–Cu Alloys: Combining Reaction Kinetics Experiments with DFT Calculations and Microkinetic Modeling

Cu-rich Pt–Cu bimetallic catalysts are among the most promising candidates for actively catalyzing the hydrodechlorination of 1,2-dichloroethane (1,2-DCA) toward ethylene production. Combining reaction kinetics experiments with density functional theory (DFT) calculations and mean-field microkinetic modeling, we present a systematic mechanistic study for 1,2-DCA hydrodechlorination on Cu-rich Pt–Cu alloy catalysts. Our DFT (PBE+(TS+SCS)) results suggest that increasing Cu content in the Pt–Cu alloy destabilizes C 2 -species adsorption while stabilizing the binding of atomic chlorine. The reaction energetics of all the elementary steps in the 1,2-DCA reaction network were calculated on a Pt 1 Cu 3 (111) model surface. The DFT results were then used to construct a microkinetic model, and the model-predicted reaction rates were compared with our reaction kinetics experimental results on a Cu-rich SiO 2 -supported Pt–Cu alloy catalyst through a parameter estimation procedure. Both the reaction kinetics experiments and the microkinetic model after parameter adjustments yielded 100% selectivity to ethylene. The microkinetic model pointed to a reaction pathway involving two sequential chlorine-removal steps on the Pt–Cu alloy catalyst, a mechanism distinct from the one previously identified on pure Pt/SiO 2 catalysts, which involved an initial hydrogen-removal step. Adjustments to the DFT-derived parameters indicate the possible formation of chlorine-induced Cu-enriched surface sites during 1,2-DCA hydrodechlorination conditions, sites that are more active than those encountered in the bulk Pt 1 Cu 3 (111) alloy surface. Furthermore, our study offers valuable initial insights on the 1,2-DCA hydrodechlorination reaction mechanism and the nature of the active sites on PtCu bimetallic catalysts.

1,2-Dichloroethane↗
DOE OSTI2021

CFD/FEA Co-Simulation Framework for Analysis of the Thermal Barrier Coating Design and Its Impact on the HD Diesel Engine Performance

Thermal barrier coatings (TBCs) have been investigated both experimentally and through simulation for mixing controlled combustion (MCC) concepts as a method for reducing heat transfer losses and increasing cycle efficiency, but it is still a very active research area. Early studies were inconclusive, with different groups discovering obstacles to realizing the theoretical potential. Nuanced papers have shown that coating material properties, thickness, microstructure, and surface morphology/roughness all can impact the efficacy of the thermal barrier coating and must be accounted for. Adding to the complexities, a strong spatial and temporal heat flux inhomogeneity exists for mixing controlled combustion (diesel) imposed onto the surfaces from the impinging flame jets. In support of the United States Department of Energy SuperTruck II program goal to achieve 55% brake thermal efficiency on a heavy-duty diesel engines, this study sought to develop a deeper insight into the inhomogeneous heat flux from mixing controlled combustion on thermal barrier coatings and to infer concrete guidance for designing coatings. To that end, a co-simulation approach was developed that couples high-fidelity computational fluid dynamics (CFD) modeling of in-cylinder processes and combustion, and finite element analysis (FEA) modeling of the thermal barrier-coated and metal engine components to resolve spatial and temporal thermal boundary conditions. The models interface at the surface of the combustion chamber; FEA modeling predicts the spatially resolved surface temperature profile, while CFD develops insights into the effect of the thermal barrier coating on the combustion process and the boundary conditions on the gas side. The paper demonstrates the capability of the framework to estimate cycle impacts of the temperature swing at the surface, as well as identify critical locations on the piston/thermal barrier coating that exhibit the highest charge temperature and highest heat fluxes. In addition, the FEA results include predictions of thermal stresses, thus enabling insight into factors affecting coating durability. An example of the capability of the framework is provided to illustrate its use for investigating novel coatings and provide deeper insights to guide future coating design.

42 ENGINEERING↗
DOE OSTI2024

Contrasting Effects of Catecholate and Hydroxamate Siderophores on Molybdenite Dissolution

Molybdenum (Mo) is essential for many enzymes but is often sequestered within minerals, rendering it not readily bioavailable. Metallophores, metabolites secreted by microorganisms and plants, promote mineral dissolution to increase metal bioavailability. However, interactions between metallophores and Mo-bearing minerals remain unclear. Here, in this study, catecholate protochelin and hydroxamate desferrioxamine B (DFOB) were utilized to examine their effects on dissolution of the common Mo-bearing mineral, molybdenite (MoS 2 ), under both oxic and anoxic conditions. Protochelin promoted molybdenite dissolution under oxic conditions with the formation of MoO 3 on the surface and Mo-siderophore complexes in solution. This was attributed to air-oxidation of both molybdenite and protochelin, as evidenced by lack of dissolution under anoxic conditions but enhanced dissolution by either pre-oxidized protochelin or pre-oxidized molybdenite. Liquid chromatography-mass spectroscopy (LC-MS), X-ray photoelectron spectroscopy (XPS), and Time-of-flight secondary ion mass spectrometry (ToF-SIMS) analyses revealed degradation of protochelin and adsorptions of its byproducts on molybdenite surface to promote dissolution. Conversely, DFOB inhibited molybdenite dissolution under both oxic and anoxic conditions, likely attributed to surface adsorption of DFOB and its weak complexation with Mo(VI) at the circumneutral pH. This work highlights the need to consider the balance between promoting and inhibitory effects of different metallophores on Mo-mineral dissolution.

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

Off-Stoichiometric Restructuring and Sliding Dynamics of Hexagonal Boron Nitride Edges in Conditions of Oxidative Dehydrogenation of Propane

Boron-containing materials, such as hexagonal boron nitride (h-BN), recently shown to be active and selective catalysts for the oxidative dehydrogenation of propane (ODHP), have been shown to undergo significant surface oxyfunctionalization and restructuring. Although experimental ex situ studies have probed the change in chemical environment on the surface, the structural evolution of it under varying reaction conditions has not been established. Herein, we perform global optimization structure search with a grand canonical genetic algorithm to explore the chemical space of off-stoichiometric restructuring of the h-BN surface under ambient as well as ODHP-relevant conditions. A grand canonical ensemble representation of the surface is established, and the predicted 11B solid-state NMR spectra are consistent with previous experimental reports. In addition, we investigated the relative sliding of h-BN sheets and how it influences the surface chemistry with ab initio molecular dynamics simulations. Furthermore, the B–O linkages on the edges are found to be significantly strained during the sliding, causing the metastable sliding configurations to have higher reactivity toward the activation of propane and water.

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

SURFACE FINISHING AND ELECTROLESS NICKEL PLATING OF ADDITIVELY MANUFACTURED (AM) METAL COMPONENTS

This study investigates the application of electroless nickel deposition on additively manufactured stainless steel samples. Current additive manufacturing (AM) technologies produce metal components with a rough surface. Rough surfaces generally exhibit fatigue characteristics, increasing the probability of initiating a crack or fracture to the printed part. For this reason, the direct use of as-produced parts in a finished product cannot be actualized, which presents a challenge. Post-processing of the AM parts is therefore required to smoothen the surface. This study analyzes chempolish (CP) and electropolish (EP) surface finishing techniques for post-processing AM stainless steel components CP has a great advantage in creating uniform, smooth surfaces regardless of size or part geometry EP creates an extremely smooth surface, which reduces the surface roughness to the sub-micrometer level. In this study, we also investigate nickel deposition on EP, CP, and as-built AM components using electroless nickel solutions. Electroless nickel plating is a method of alloy treatment designed to increase manufactured component's hardness and surface resistance to the unrelenting environment. The electroless nickel plating process is more straightforward than its counterpart electroplating.. We use low-phosphorus (2-5% P), medium-phosphorus (6-9% P), and high-phosphorus (10-13% P). These Ni deposition experiments were optimized using the L9 Taguchi design of experiments (TDOE), which compromises the prosperous content in the solution, surface finish, plane of the geometry, and bath temperature. The pre-and post-processed surface of the AM parts was characterized by KEYENCE Digital MicroscopeVHX-7000 and Phenom XL Desktop SEM. The experimental results show that electroless nickel deposition produces uniform Ni coating on the additively manufactured components up to 20 μm per hour. Mechanical properties of as-built and Ni coated AM samples were analyzed by applying a standard 10 N scratch test. Nickel coated AM samples were up to two times scratch resistant compared to the as-built samples. This study suggests electroless nickel plating is a robust viable option for surface hardening and finishing AM components for various applications and operating conditions.

Keywords: additive manufacturing, fatigue, chempol↗
DOE OSTI2020

Coverage-Dependent Adsorption of Hydrogen on Fe(100): Determining Catalytically Relevant Surface Structures via Lattice Gas Models

Hydrogen adatoms are a critical surface species for several reactions catalyzed by Fe surfaces such as Fischer–Tropsch synthesis, ammonia synthesis, and the hydrodeoxygenation of biomass-derived oxygenates. Parameterizing the energetics for H/Fe in terms of both coverage and configuration space can significantly aid in the development of multiscale models as well as provide atomic level insight into the dominant surface structures present under realistic reaction conditions. Here, we construct a lattice gas model for H/Fe(100), where the lateral interactions are determined from first-principles using density functional theory. Using 950 symmetrically unique H/Fe(100) configurations, we generate a cluster expansion with a predictive accuracy in terms of surface energy of 3.8 meV/site over a coverage range from 0 to 3 monolayers. Ten electronic ground state structures are identified from this thorough scan (including the structures at 0 and 3 monolayers), which were subsequently used to generate ab initio phase diagrams under a range of temperatures and pressures. Under reaction conditions typical of Fischer–Tropsch synthesis, ammonia synthesis, and biomass oxygenate hydrodeoxygenation, we find that the 1.0 monolayer structure is dominant. Furthermore, examination of the total H–H lateral interactions for the H/Fe(100) electronic ground state structures shows that H/Fe(100) can be accurately modeled via a mean-field ideal lattice gas model for coverages less than 1.0 monolayers. Altogether, this work enables the incorporation of H–H lateral interactions on Fe(100) into multiscale models, via either mean-field or site-dependent techniques, and provides atomic insight into the catalytically relevant H/Fe(100) structures for a range of heterogeneous reactions.

09 BIOMASS FUELS↗
DOE OSTI2020

Calculation of the detonation state of HN 3 with quantum accuracy

HN 3 is a unique liquid energetic material that exhibits ultrafast detonation chemistry and a transition to metallic states during detonation. We combine the Chebyshev interaction model for efficient simulation (ChIMES) many-body reactive force field and the extended-Lagrangian multiscale shock technique molecular dynamics method to calculate the detonation properties of HN 3 with the accuracy of Kohn–Sham density-functional theory. ChIMES is based on a Chebyshev polynomial expansion and can accurately reproduce density-functional theory molecular dynamics (DFT-MD) simulations for a wide range of unreactive and decomposition conditions of liquid HN 3 . We show that addition of random displacement configurations and the energies of gas-phase equilibrium products in the training set allows ChIMES to efficiently explore the complex potential energy surface. Schemes for selecting force field parameters and the inclusion of stress tensor and energy data in the training set are examined. Structural and dynamical properties and chemistry predictions for the resulting models are benchmarked against DFT-MD. We demonstrate that the inclusion of explicit four-body energy terms is necessary to capture the potential energy surface across a wide range of conditions. Our results generally retain the accuracy of DFT-MD while yielding a high degree of computational efficiency, allowing simulations to approach orders of magnitude larger time and spatial scales. The techniques and recipes for MD model creation we present allow for direct simulation of nanosecond shock compression experiments and calculation of the detonation properties of materials with the accuracy of Kohn–Sham density-functional theory.

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

SNF Interim Storage Canister Corrosion and Surface Environment Investigations (FY2020 Status Report)

This progress report describes work performed during FY20 at Sandia National Laboratories (SNL) to assess the localized corrosion performance of container/cask materials used in the interim storage of spent nuclear fuel (SNF). Of particular concern is stress corrosion cracking (SCC), by which a through-wall crack could potentially form in a canister outer wall over time intervals that are shorter than possible dry storage times. Work in FY20 further defined our understanding of the potential chemical and physical environment present on canister surfaces, evaluated the relationship between the environment and the resultant corrosion that occurs, and initiated crack growth rate testing under relevant environmental conditions. In FY20, work to define dry storage canister surface environments included several tasks. First, collection of dust deposition specimens from independent spent fuel storage installation (ISFSI) site locations helped to establish a more complete understanding of the potential chemical environment formed on the canister. Second, the predicted evolution of canister surface relative humidity RH) values was estimated using ISFSI site weather data and the horizontal canister thermal model used by the SNL probabilistic SCC model. These calculations determined that for typical ISFSI weather conditions, seasalt deliquescence to produce MgCl2-rich brines could occur in less than 20 years at the coolest locations on the canister surface, and, even after nearly 300 years, conditions for NaCl deliquescence (75% RH) are not reached. This work illustrates the importance of understanding the stability of MgCl2-rich brines on the heated canister surface, and the potential impact of brine composition on corrosion processes, including pitting and stress corrosion cracking. In an additional study, the description of the canister surface environment was refined in order to define more realistic corrosion testing environments including diurnal cycles, soluble salt chemistries, and inert mineral particles. The potential impacts of these phenomena on canister corrosion are being evaluated experimentally. Finally, work over the past few years to evaluate the stability of magnesium chloride brines continued in FY20. MgCl2 degassing experiments were carried out, confirming that MgCl2 brines slowly degas HCl on heated surfaces, converting to less deliquescent magnesium hydroxychloride phases and potentially leading to brine dryout.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗