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.

At least 1,261 records · Page 70

Transient heat transfer of superfluid He 4 in nonhomogeneous geometries: Second sound, rarefaction, and thermal layer

Transient heat transfer in superfluid 4 He (He II) is a complex process that involves the interplay of the unique counterflow heat-transfer mode, the emission of second-sound waves, and the creation of quantized vortices. Many past researches focused on homogeneous heat transfer of He II in a uniform channel driven by a planar heater. Here, we report our systematic study of He II transient heat transfer in nonhomogeneous geometries that are pertinent to emergent applications. By solving the He II two-fluid equations of motion coupled with Vinen's equation for vortex-line density, we examine and compare the characteristics of transient heat transfer from planar, cylindrical, and spherical heaters in He II. Our results show that as the heater turns on, an outgoing second-sound pulse emerges, within which the vortex-line density grows rapidly. These vortices attenuate the second sound and result in a heated He II layer in front of the heater, i.e., the thermal layer. In the planar case where the vortices are created throughout the space, the second-sound pulse is continuously attenuated, leading to a thick thermal layer that diffusely spreads following the heat pulse. On the contrary, in the cylindrical and the spherical heater cases, vortices are created mainly in a thin thermal layer near the heater surface. As the heat pulse ends, a rarefaction tail develops following the second-sound pulse, in which the temperature drops. This rarefaction tail can promptly suppress the thermal layer and take away the deposited thermal energy. The effects of the heater size, heat flux, pulse duration, and temperature on the thermal-layer dynamics are discussed. We also show how the peak heat flux for the onset of boiling in He II can be studied in our model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Zero-Valent Iron Permeable Reactive Barrier to Remediate Volatile Organic Compounds in Groundwater - 20209

The US DOE-SRS, the US EPA, and the South Carolina Department of Health and Environmental Control determined it was appropriate to perform a non-time critical removal action at the P-Area Groundwater Operable Unit at the SRS to reduce the mass and downgradient transport of trichloroethylene in the P-Area groundwater plume. Contaminated groundwater discharges to a nearby stream, Steel Creek, within the SRS boundaries, resulting in trichloroethylene concentrations above the maximum contaminant level. Impact to surface water is limited in areal extent and supported by recently collected characterization data. The P-Area Groundwater Operable Unit encompasses the groundwater beneath an industrial area within SRS, P Area, where the P-Reactor once operated. The boundaries of the P-Area Groundwater Operable Unit extend northwest to Steel Creek, northeast toward PAR Pond, and southeast to Meyers Branch. Groundwater in the Upper Three Runs Aquifer of the P-Area Groundwater Operable Unit has been impacted by reactor and facility operations between 1954 and 1991, including tritium and volatile organic compounds. The P-Area surface units contributing to groundwater contamination were remediated as part of the P-Area Operable Unit in 2011. The P-Reactor closure is one of the first of its kind in the DOE Complex and is one of only a few full-sized production reactors in the US to undergo completion of final closure activities. The nature and extent of groundwater contamination was determined using a variety of investigative approaches such as groundwater monitoring wells, direct-push technology, and surface water samples. Groundwater contamination associated with trichloroethylene is primarily exhibited in a narrow plume that extends from the source area at P-Reactor and west to Steel Creek. Maximum contaminant level exceedances in groundwater occur over an area of ∼6.9 hectares for trichloroethylene with concentrations as high as 7.7 milligrams per liter. To the west of the P-Area facility area, the trichloroethylene groundwater plume is controlled by a buried geologic feature, assumed to be an old stream bed, that further narrows the groundwater plume in what has been designated as the 'neck area.' This narrowing of the groundwater plume provides an ideal location for a treatment barrier. The non-time critical removal action alternative chosen is to install a zero-valent iron permeable reactive barrier within the neck area of the trichloroethylene groundwater plume, perpendicular to groundwater flow direction. This technology will provide a treatment barrier that will reduce trichloroethylene groundwater concentrations by 90% and has an anticipated useful life of at least 25 years. A pre-design investigation was performed in the neck area to confirm site lithology, hydrogeology, geochemistry, and extent of trichloroethylene contamination prior to a final design. A treatability study, conducted as part of the pre-design investigation, indicated that the subsurface and groundwater in the PArea Groundwater Operable Unit is compatible with the zero-valent iron and will not lead to excessive buildup from mineralization/precipitation or biofouling. Probabilistic modeling was conducted using field and laboratory data to determine the expected performance of the zero-valent iron permeable reactive barrier. The model simulations indicated that a 3.81-centimeter thick barrier would provide greater than 90% reduction of trichloroethylene groundwater concentration. The final design of the zero-valent iron permeable reactive barrier is a barrier that will extend 80.5 linear meters in a 'zigzag' orientation to best transect the trichloroethylene plume and account for varying groundwater flow. The barrier will be installed from 13.7 meters below ground surface to 41.1 meters below ground surface for 65.8 linear meters and from 13.7 meters below ground surface to 36.6 meters below ground surface for 14.6 linear meters, the base of which is 'keyed' into a low permeability zone. The barrier is designed to a thickness of 10.2 centimeters, which was determined to reduce the trichloroethylene groundwater concentrations by greater than 90% with a safety factor of 2.67. A total of approximately 689 metric tons of zero-valent iron will be injected through 22 injection wells spaced 3.66 meters apart, using guar to suspend the zero-valent iron. Zero-valent iron permeable reactive barrier construction will be monitored through 23 installed resistivity receivers offset 7.32 meters from the zerovalent iron permeable reactive barrier. The zero-valent iron will be energized with a low-voltage 100 Hertz signal during injection and will be monitored using the resistivity receivers to ensure complete coalescence of the zero-valent iron permeable reactive barrier. The zero-valent iron was sized to have a hydraulic conductivity greater than the natural subsurface, thus promoting groundwater flow through the barrier. As contaminated groundwater contacts the zero-valent iron, volatile organic compounds, including trichloroethylene, are immediately degraded to harmless compounds such as ethylene. The performance of the zero-valent iron permeable reactive barrier will be monitored using three upgradient monitoring well clusters, six downgradient monitoring well clusters, and four in-wall monitoring wells. The in-wall monitoring wells will indicate immediate reduction of trichloroethylene mass in the groundwater and allow analyses of the zero-valent iron permeable reactive barrier health. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

A parameterization of sub-grid topographical effects on solar radiation in the E3SM Land Model (version 1.0): implementation and evaluation over the Tibetan Plateau

Abstract. Topography exerts significant influences on the incoming solar radiation at the land surface. A few stand-alone regional and global atmospheric models have included parameterizations for sub-grid topographic effects on solar radiation. However, nearly all Earth system models (ESMs) that participated in the Coupled Model Intercomparison Project (CMIP6) use a plane-parallel (PP) radiative transfer scheme that assumes that the terrain is flat. In this study, we incorporated a well-validated sub-grid topographic (TOP) parameterization in the Energy Exascale Earth System Model (E3SM) Land Model (ELM) version 1.0 to quantify the effects of sub-grid topography on solar radiation flux, including the shadow effects and multi-scattering between adjacent terrain. We studied the role of sub-grid topography by performing ELM simulations with the PP and TOP schemes over the Tibetan Plateau (TP). Additional ELM simulations were performed at multiple spatial resolutions to investigate the role of spatial scale on sub-grid topographic effects on solar radiation. The Moderate Resolution Imaging Spectroradiometer (MODIS) data was used to compare with the ELM simulations. The results show that topography has non-negligible effects on surface energy budget, snow cover, snow depth, and surface temperature over the TP. The absolute differences in surface energy fluxes for net solar radiation, latent heat flux, and sensible heat flux between TOP and PP exceed 20, 10, and 5 W m−2, respectively. The differences in land surface albedo, snow cover fraction, snow depth, and surface temperature between TOP and PP exceed 0.1, 0.1, 10 cm, and 1 K, respectively. The magnitude of the sub-grid topographic effects is dependent on seasons and elevations and is also sensitive to the spatial scales. Although the sub-grid topographic effects on solar radiation are larger with more spatial details at finer spatial scales, they cannot be simply neglected at coarse spatial scales. When compared to MODIS data, incorporating the sub-grid topographic effects overall reduces the biases of ELM in simulating surface energy balance, snow cover, and surface temperature, especially in the high-elevation and snow-covered regions over the TP. The inclusion of sub-grid topographic effects on solar radiation parameterization in ELM will contribute to advancing our understanding of the role of the surface topography on terrestrial processes over complex terrain.

Hao, Dalei (ORCID:0000000334979774)↗

Oxidation Dynamics of Supported Catalytic Cu Clusters: Coupling to Fluxionality

Copper oxide nanoclusters have a wide range of catalytic applications, such as the selective oxidation of hydrocarbons. O 2 binding to the catalyst, activation, and release upon reagent oxidation are key events in these catalytic chemistries. Furthermore, these events are expected to be accompanied by significant structural changes of the Cu clusters, because O atoms integrate into the cluster, rather than bind to its surface. Topping the complexity of the problem, partially oxidized Cu clusters are known to exhibit strong fluxionality and feature diverse and interconverting structures and oxygen contents in conditions of oxidative dehydrogenation (ODH). Hence, a significant dynamic coupling between the “hot” O 2 molecule impacting the cluster at reaction temperatures and the cluster fluxionality can be expected. In this work, we focus on the dynamics of dioxygen integration into a partially oxidized Cu cluster supported on hydroxylated amorphous alumina–a system recently reported to be an exceptionally selective catalyst for cyclohexane ODH with very little CO produced, whose mechanistic underpinnings are of utmost interest. The statistics over a swarm of adsorption and scattering trajectories where O 2 hits various sites on the cluster at reaction temperature shows that the O 2 binding does not only follow the minimal energy paths. O 2 also rarely integrates into the cluster in a single step and instead first binds to a single Cu atom via either an η 1 -O 2 or an η 2 -O 2 mode. Surprisingly, this step often has a higher barrier than the subsequent O 2 integration and dissociation, which in turn take multiple steps and complete the oxidation process. Dynamic trajectories starting from the key transition state of integration of the adsorbed O 2 can also lead to different intermediate structures during or right after the dissociation, due to the energy released from the transition state and the thermal intracluster effects. From these activated O 2 chemisorbed structures, O 2 dissociation occurs with moderate barriers (~0.5 eV), producing multiple final oxidized Cu 4 O 4 states. Hence, a diversity of reaction profiles for the attack of supported Cu cluster by O 2 emerges due to the dynamic effects, with implications for mechanisms, kinetic models, and catalyst design principles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Design and Production for Maximum Structural Efficiency With Respect to Fiber Orientation With Increased Understanding of Hybrid Fiber Flow Behavior

Discontinuous fiber-reinforced thermoplastic composites have gained considerable attention in automotive, aerospace, and other industries, due to their high-rate of production combined with their ability to attain complex and intricate shapes. Among other high-rate thermoplastic manufacturing processes, injection-molding is one of the most common manufacturing methods due to fast production and high surface finishing of complex geometries. Fiber orientation in discontinuous fiber composites plays a pivotal role in determining the mechanical, electrical, and thermomechanical properties, underscoring the necessity to comprehend fiber orientation in injection molded parts. Among different fiber types, glass and carbon fibers are most common in the composite industries. The recent trend of hybrid composites comprising both glass fiber (GF) and carbon fiber (CF) is also gaining importance in the automotive industry. Hybrid fiber options allow designers to optimize the balance between glass and carbon fibers by leveraging the high durability and low cost of GF while the strength and lightweight properties of CF. Consequently, comparing the fiber oriented distribution (FOD) of injection molded composites containing GF, CF, and a hybrid of GF/CF is critical to investigating the local mechanical properties of intricate structures for high-end applications. In Phase I of this project, FOD in injection-molded panels with respect to distance from the gate was analyzed using X-ray computed tomography (X-CT) for GF, CF, and hybrid CF/GF (CGF) reinforced nylon 66. To understand the reason behind the FOD with different fiber types, computational fluid dynamics (CFD) and rheology were performed. Samples were extracted at three locations: near the gate, center, and opposite end. Thickness of the layers of typical skin-shell-core type FOD varies with fiber type and location. GF achieved flow direction alignment (in shell) earlier than viscous CF and CGF near the gate, whereas CF showed the highest flow-direction alignment at the center due to shear induced orientation. At the opposite end, GF experienced more backflow than others indicating faster mold filling owing to its lower viscosity. Hybrid CGF exhibited GF-dominated center and CF-dominated end region. The numerical model used to obtain FOD and rheological predictions for the CF and GF composites served to corroborate the trends observed in the experimental trials. The FOD responses across fiber types and location were reflected in their longitudinal and transverse properties. Only GF showed higher longitudinal modulus over transverse modulus near the gate attributed to rapid alignment, whereas CF and CGF exhibited opposite trend. However, fountain flow enhanced the longitudinal modulus over transverse modulus with the distance for all, particularly for CF. This study offers insights into mold filling behavior of different fibers which are critical in optimizing injection molding conditions for tailored final properties.

36 MATERIALS SCIENCE↗

Response of Subsurface Nitrogen-Cycling Microbial Communities to Environmental Fluctuations (Final Technical Report)

Riparian floodplains are dynamic ecosystems linking terrestrial and riverine systems. These floodplains experience hydrological shifts such as changes in water table height, flooding, and drought and can be ‘hotspots’ of biogeochemical cycling due to shifting sediment moisture (and saturation) and subsurface exchanges of water, nutrients, and other compounds across different sediment layers. Subsurface microbial communities are the primary drivers of biogeochemical processes in floodplains, and thus their structure and function can directly influence both surface and groundwater quality. The microbial nitrogen (N) cycle is particularly important in floodplains as it affects nutrient availability and removal. Two functional guilds of chemoautotrophic (i.e. CO2-fixing) microorganisms are responsible for the first oxidative step of the N cycle, nitrification: ammonia-oxidizing archaea (AOA) and bacteria (AOB) catalyze the oxidation of ammonia to nitrite, while nitrite-oxidizing bacteria (NOB) oxidize nitrite to nitrate. Despite the critical role nitrification plays in N-cycling in both terrestrial and aquatic ecosystems, our understanding of the diversity, ecophysiology, and activity of nitrifying organisms in subsurface floodplain soils/sediments is extremely limited. To help address this critical knowledge gap, the overarching goal of this project was to determine how shifts in key environmental parameters and gradients impact microbial N-cycling communities/processes, with particular emphasis on nitrification, within hydrologically-variable floodplain sediments in the Wind River Basin near Riverton, Wyoming. The three specific objectives of this project were to: (1) to associate in situ environmental drivers of N cycling with distinct functional guilds; (2) determine the guild response to variation in key ecosystem drivers; and (3) develop a dynamic ecosystem model of the microbial N cycle with the Riverton subsurface using community genomic and biogeochemical data collected in the first two objectives. Over the course of this project, we employed both 16S rRNA gene amplicon sequencing and genome-resolved metagenomics to examine the phylogenetic diversity and metabolic potential of subsurface nitrifier communities within 68 samples collected across multiple sites, depths, and time points within the Riverton floodplain, allowing for both spatial and temporal investigations at different scales. This project benefitted tremendously from recent advances in high-throughput sequencing technologies coupled with dramatic improvements in the computational tools and algorithms available for analyzing such large, complex genomic datasets. By pairing these cutting-edge genomic approaches with depth-resolved sampling and detailed geochemical analyses of the Riverton floodplain, we have gained novel insights into the structure and function of subsurface nitrifier communities in relation to both hydrology and biogeochemistry. This project resulted in the most detailed and comprehensive characterization of N-cycling floodplain microbial communities to date and will hopefully inspire and pave the way for future studies using similar approaches in other floodplains. Indeed, such information is critical for understanding subsurface biogeochemical cycling and how elemental stores are altered from perturbations initiated by the water cycle within floodplains. Finally, because of the terrestrial-aquatic nature of the Riverton floodplain, results from this project are also of relevance to disciplines such as soil science, estuarine science, limnology & oceanography, biogeochemistry, geobiology, environmental engineering, as well as genomics and data science.

54 ENVIRONMENTAL SCIENCES↗

Engineered Barrier System R&D and International Collaborations – LANL (FY20): Spent Fuel and Waste Disposition

The United States Department of Energy’s Spent Fuel and Waste Disposition (SFWD) program is investigating the design and safety function of generic nuclear geologic repositories in a variety of geologic settings (salt, argillite, and crystalline rock). Different configurations and loadings of spent nuclear fuel and waste within disposal canisters are also being investigated, some of which have the potential to generate repository temperatures higher than previously considered (i.e., temperatures >100ºC) by foreign and domestic concepts. This report expands on engineered barrier material stability in a high temperature crystalline rock repository through high temperature hydrothermal experiments. Experiments were designed to develop engineered barrier system (EBS) concepts in a hightemperature crystalline environment in 1) bentonite-Grimsel Granodiorite interactions, 2) bentonite-cement reactions, and 3) interaction between waste canister materials and bentonite. Experiment results are applied to understanding long-term repository performance in terms of radionuclide isolation. One hydrothermal experiment was completed in the rocking autoclaves at LANL in FY-20: IEBS-6 (Grimsel Granodiorite + Wyoming bentonite + cured ordinary Portland cement + Grimsel Granodiorite synthetic groundwater, 250ºC/150 bar, 8 weeks). Several other experiments were planned but were delayed due to the COVID-19 pause in laboratory work. In addition, some characterization of the reaction products of experiments conducted in FY-20 was prevented by COVID-19. Mineral phase chemistry was not measured via electron microprobe analyses for IEBS-6. Quantitative X-ray diffraction results from IEBS-6 and HBT-1 were not completed. The missing analyses and discussion of the results will be included in next year’s report. New characterization that was completed includes scanning electron microscopy of reaction products from IEBS-6, quantitative X-ray diffraction results are presented from IEBS-1 through IBES-5, preliminary scanning electron microscope images and chemical analyses for IEBS-6, X-ray diffraction of the clay fraction from all experiments, measurement/imaging of mineral growth on the surface of steel coupons. Major observations pertaining to bentonite stability in a Grimsel Granodiorite environment include the stability of Na-montmorillonite at 250°C, the formation of trace CSH phases, and the formation of bentonite colloids on experiment cooling. The addition of a cured chip of Portland cement to the bentonite-Grimsel system results in slightly higher pH values and the formation of diverse secondary mineral phases that were not observed in the previous experiments (e.g., analcime, garronite, CSH phases). The new characterization efforts related to the interaction of stainless-steel coupons and bentonite clay focused on thickness and mineralogy of phases that formed at the steel surface. In Wyoming bentonite + Grimsel Granodiorite systems, newly formed minerals at the bentonitesteel coupon interface included alteration of the outermost steel edge to Fe,Cr-oxide phases, followed by Fe-rich phyllosilicates (Fe-saponite, chlorite) and interbedded Fe,Cr,Ni-sulfide phases (pentlandite). Hydrothermal experiments were completed to assess uranium-carbonate complexation at conditions relevant to high-temperature disposal. Autoclave solubility experiments were conducted at 150 to 250ºC with a range of carbonate and uranium concentrations. The experiment results were characterized via situ UV-Visible spectroscopy and synchrotron-based in situ XAS techniques. Results show a significant decrease in the stability of uranyl-carbonate complexes at temperatures above 100ºC along with the precipitation of uranium oxides. Further, at T > 200ºC, results show that uranyl-hydroxyl complexes control solubility of uranium instead of the previously predicted uranyl-carbonate species. These results are significant for understanding the mobility of uranium in the EBS, which will likely contain carbonate-rich fluids. International research efforts focused on three main areas: 1) participation in international conferences, 2) building collaborations with foreign repository programs, and 3) the initiation of an experimental program to complement the full-scale HotBENT test at the Grimsel test site. This experiment included Wyoming bentonite + low carbon steel + Grimsel Granodiorite synthetic groundwater and was run at the planned maximum temperature of the HotBENT test (200°C). Complete characterization of reaction products was hindered by disruptions to laboratory work but will be reported in the next FY. The experimental results obtained in FY-20 continue to document the wide-ranging effects of bulk composition and pressure-temperature conditions in the mineralogical and geochemical evolution of a high-temperature repository environment. Concepts developed will be used to inform models of long-term material stability in a generic crystalline rock-hosted repository.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Interfacial Electromechanics Predicts Phase Behavior of 2D Hybrid Halide Perovskites

Quasi-two-dimensional (2D) mixed-cation hybrid halide perovskites (A' 2 A N–1 M N X 3N+1 ; A' = large organic molecule with cationic group, A = [Cs + , CH 3 NH 3 + , HC(NH 2 ) 2 + ], M = [Pb, Sn, Ge], X = [I – , Br – , Cl – ]) have rapidly emerged as candidates to improve the structural stability and device lifetime of 3D perovskite semiconductor devices under operating conditions. The addition of the large A' cation to the traditional AMX 3 structure introduces several synthetic degrees of freedom and breaks M–X bonds, giving rise to peculiar critical phase behavior in the phase space of these complex materials. In this work, we propose a thermodynamic model parametrized by first-principles calculations to generate the phase diagram of 2D and quasi-2D perovskites (q-2DPKs) based on the mechanics and electrostatics of the interface between the A' cations and the metal halide octahedral network. Focusing on the most commonly studied methylammonium lead iodide system where A' is n-butylammonium (BA; CH 3 (CH 2 ) 3 NH 3 + ), we find that the apparent difficulty in synthesizing phase-pure samples with a stoichiometric index N > 5 can be attributed to the energetic competition between repulsion of opposing interfacial dipole layers and mechanical relaxation induced by interfacial stress. Our model shows quantitative agreement with experimental observations of the maximum phase-pure stoichiometric index (N crit ) and explains the nonmonotonic evolution of the lattice parameters with increasing stoichiometric index (N). This model is generalizable to the entire family of q-2DPKs and can guide the design of photovoltaic and optical materials that combine the structural stability of the q-2DPKs while retaining the charge carrier properties of their 3D counterparts.

14 SOLAR ENERGY↗

Evaluation of native Earth system model output with ESMValTool v2.6.0

Earth system models (ESMs) are state-of-the-art climate models that allow numerical simulations of the past, present-day, and future climate. To extend our understanding of the Earth system and improve climate change projections, the complexity of ESMs heavily increased over the last decades. As a consequence, the amount and volume of data provided by ESMs has increased considerably. Innovative tools for a comprehensive model evaluation and analysis are required to assess the performance of these increasingly complex ESMs against observations or reanalyses. One of these tools is the Earth System Model Evaluation Tool (ESMValTool), a community diagnostic and performance metrics tool for the evaluation of ESMs. Input data for ESMValTool needs to be formatted according to the CMOR (Climate Model Output Rewriter) standard, a process that is usually referred to as “CMORization”. While this is a quasi-standard for large model intercomparison projects like the Coupled Model Intercomparison Project (CMIP), this complicates the application of ESMValTool to non-CMOR-compliant climate model output. In this paper, we describe an extension of ESMValTool introduced in v2.6.0 that allows seamless reading and processing of “native” climate model output, i.e., operational output produced by running the climate model through the standard workflow of the corresponding modeling institute. This is achieved by an extension of ESMValTool's preprocessing pipeline that performs a CMOR-like reformatting of the native model output during runtime. Thus, the rich collection of diagnostics provided by ESMValTool is now fully available for these models. For models that use unstructured grids, a further preprocessing step required to apply many common diagnostics is regridding to a regular latitude–longitude grid. Extensions to ESMValTool's regridding functions described here allow for more flexible interpolation schemes that can be used on unstructured grids. Currently, ESMValTool supports nearest-neighbor, bilinear, and first-order conservative regridding from unstructured grids to regular grids. Example applications of this new native model support are the evaluation of new model setups against predecessor versions, assessing of the performance of different simulations against observations, CMORization of native model data for contributions to model intercomparison projects, and monitoring of running climate model simulations. For the latter, new general-purpose diagnostics have been added to ESMValTool that are able to plot a wide range of variable types. Currently, five climate models are supported: CESM2 (experimental; at the moment, only surface variables are available), EC-Earth3, EMAC, ICON, and IPSL-CM6. As the framework for the CMOR-like reformatting of native model output described here is implemented in a general way, support for other climate models can be easily added.

58 GEOSCIENCES↗

Feature engineering descriptors, transforms, and machine learning for grain boundaries and variable-sized atom clusters

Abstract Obtaining microscopic structure-property relationships for grain boundaries is challenging due to their complex atomic structures. Recent efforts use machine learning to derive these relationships, but the way the atomic grain boundary structure is represented can have a significant impact on the predictions. Key steps for property prediction common to grain boundaries and other variable-sized atom clustered structures include: (1) describing the atomic structure as a feature matrix, (2) transforming the variable-sized feature matrix to a fixed length common to all structures, and (3) applying a machine learning algorithm to predict properties from the transformed matrices. We examine how these steps and different combinations of engineered features impact the accuracy of grain boundary energy predictions using a database of over 7000 grain boundaries. Additionally, we assess how different engineered features support interpretability, offering insights into the physics of the structure-property relationships.

36 MATERIALS SCIENCE↗

Stochastic 3D reconstruction of cracked polycrystalline NMC particles using 2D SEM data

Li-ion battery performance is strongly influenced by the 3D microstructure of its cathode particles. Cracks within these particles develop during calendaring and cycling, reducing connectivity but increasing reactive surface, making their impact on battery performance complex. Understanding these contradictory effects requires a quantitative link between particle morphology and battery performance. However, informative 3D imaging techniques are time-consuming, costly and rarely available, such that analyses often have to rely on 2D image data. This paper presents a novel stereological approach for generating virtual 3D cathode particles exhibiting crack networks that are statistically equivalent to those observed in 2D sections of experimentally measured particles. Consequently, 2D image data suffices for deriving a full 3D characterization of cracked cathodes particles. Such virtually generated 3D particles could serve as geometry input for spatially resolved electro-chemo-mechanical simulations to enhance our understanding of structure-property relationships of cathodes in Li-ion batteries.

36 MATERIALS SCIENCE↗

The Challenge of Achieving Site Closure of a Complex Groundwater Plume - 20531

Over the last three decades, DOE has made substantial progress in nearly every area of nuclear waste cleanup and has closed and transitioned 91 of the 107 sites in the EM complex. Several of the remaining groundwater plumes at Hanford, Los Alamos, Paducah, Oak Ridge, and Savannah River are complex and tend to have unique site or contaminant characteristics which provide an almost impossible challenge to remediate to regulatory levels. As a result, these complex sites will require a large investment to remediate. In order for DOE to achieve regulatory closure, a technically sound case must be developed and presented to regulators and stakeholders that clearly demonstrates that remedial goals have been or will be irreversibly achieved or that there is an alternate acceptable end-state. Five lines of evidence are required: - A conceptual site model that is as complete as possible, with the knowledge that it may continue to evolve with time; - Demonstration that trends in contaminant concentrations, or flux, versus time indicate that remedial goals will be achieved on schedule or within a reasonable time thereafter and that trends that deviate from this path are explainable by the conceptual model; - Demonstration that remobilization of attenuated contaminants is unlikely to occur in the future; - Explicit consideration of remedial contingencies, should deviations in the path toward remedial goals occur; - Development of an effective long-term monitoring plan to detect any conditions that might alter the path toward remedial goals. SRS has made significant progress in the remediation of contamination associated with the F-area and H-Area Seepage Basins, and it is currently developing a strategic closure path that addresses each of these lines of evidence and identifies both basic and applied scientific issues that will need to be addressed in order to obtain regulatory closure. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

The Role of Methylaryl Radicals in the Growth of Polycyclic Aromatic Hydrocarbons: The Formation of Five-Membered Rings

The regions of the C 13 H 11 potential energy surface (PES) related to the unimolecular isomerization and decomposition of the 1-methylbiphenylyl radical and accessed by the 1-/2-methylnaphthyl + C 2 H 2 reactions have been explored by ab initio G3(MP2,CC)//B3LYP/6-311G(d,p) calculations. The kinetics of these reactions relevant to the growth of polycyclic aromatic hydrocarbons (PAH) under high-temperature conditions in circumstellar envelopes and in combustion flames has been studied employing the RRKM-Master Equation approach. The unimolecular reaction of 1-methylbiphenylyl proceeding via a five-membered ring closure followed by H elimination is predicted to be very fast, on a submicrosecond scale above 1000 K and to result in the formation of an embedded five-membered ring in the 9H-fluorene product. The 1-/2-methylnaphthyl + C 2 H 2 reaction mechanism involves acetylene addition to the radical on the methylene group followed by a six- or five-membered ring closure and aromatization via an H atom loss. Despite of the complexity of the C 13 H 11 PES, these straightforward pathways are dominant in the high-temperature regime (above ~1000 K), with the prevailing products being phenalene, with a significant contribution of 1H-cyclopenta(a)naphthalene, for 1-methylnaphthyl + C 2 H 2 , and 1H-cyclopenta(b)naphthalene and 3H-cyclopenta(a)naphthalene, for 2-methylnaphthyl + C 2 H 2 . The methylnaphthyl reactions with acetylene represent a clean source of the three-ring PAHs, but they are relatively slow owing to the high entrance barriers of ~10 kcal/mol, with the rate constants of about an order of magnitude lower as compared to those for naphthyl + allene and σ-aryl + C 2 H 2 . Here, the 1-methylnaphthyl + C 2 H 2 and 2-methylnaphthyl + C 2 H 2 reactions represent prototypes for PAH growth by an extra six- and five-membered ring on a zigzag edge or a corner of PAH and the generated modified Arrhenius expressions are recommended for kinetic modeling of PAH expansion by the mechanism of acetylene addition to methylaryl radicals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Quaternary slip rates and most recent surface rupture of the Bitterroot fault, western Montana

The Bitterroot fault is a ~100-km-long Quaternary active normal fault that bounds the eastern margin of the north–south-trending Bitterroot Mountains and accommodates extension near the Intermountain Seismic Belt. Here, new detailed mapping using high-resolution topographic data derived from light detection and ranging (lidar) along the southern Bitterroot Range documents multiple generations of fault scarps in Holocene–Pleistocene deposits with vertical off sets that increase in magnitude with age. Fault mapping indicates a complex fault geometry characterized by an en echelon pattern of discontinuous segments of 41–78° east-dipping normal faults that appear to cut the older Eocene detachment fault, and locally 70–88° west-dipping antithetic normal faults. 10 Be cosmogenic exposure dating provides in situ age control for 32 surface boulders (>1 m) sampled in glacial deposits. Near Como Dam, two Pinedale-age glacial moraine sequences yield peak age distributions of 15.0 ± 0.4 ka and 16.4 ± 0.6 ka as apparent exposure ages (e = 0), and 15.4 ± 0.4 ka and 16.8 ± 0.6 ka based on the maximum allowed boulder surface erosion rate (e = 2 mm/ka). Vertical separation of 3.5 ± 0.2 m across the ~16–17 ka glacial moraine off set by the Bitterroot fault scarp yields a fault slip rate of 0.2–0.3 mm/yr. Glacial Lake Missoula highstand shorelines, inset into the ~15 ka glacial moraine and vertically off set 4.6 ± 1.6 m by an antithetic strand of the Bitterroot fault, yield fault slip rates of 0.2–0.5 mm/yr that overlap with fault slip rates on the main strand near Lake Como (0.2–0.3 mm/yr). At the Ward Creek fan located ~15 km to the north of Lake Como, two glacial debris fan sequences yield peak age distributions of 16.6 ± 0.4 ka and 62.8 ± 1.7 ka (e = 0), and 17.0 ± 0.4 ka and 69.9 ± 2.2 ka (e = 2 mm/ka). Vertical separations of 2.4 ± 0.3 m and 4.5 ± 0.2 m on the ~17 ka and ~63–70 ka fan surfaces off set by the Bitterroot fault yield fault slip rates of 0.2–0.3 mm/yr and 0.1 mm/yr, respectively. Our results indicate broadly consistent fault slip rates for the main fault segments at Lake Como (0.2–0.3 mm/yr) and the Ward Creek fan (0.1–0.3 mm/yr) with an along-strike range of 0.1–0.3 mm/yr for the southern Bitterroot fault. Fault scaling relations and evidence of multiple late Quaternary fault surface ruptures suggest the Bitterroot fault could produce a Mw ~7.2 earthquake. Structural model constraints and our slip rate results indicate both high-angle or low-angle fault geometries are possible at depth. A seismogenic low-angle fault model could generate a larger earthquake of Mw >7.2. Earthquake history is unknown for the Bitterroot fault, but fault scarps in young glacial deposits demonstrate its seismogenic potential. Data from this study suggest seismic hazards from the Bitterroot fault may pose a high level of risk to the Missoula metropolitan area, the State’s second most populous region, and major infrastructures across the Missoula and Bitterroot Valleys.

58 GEOSCIENCES↗

Durability of plasma sprayed Thermal Barrier Coatings with controlled properties part Ⅱ: Effects of geometrical curvature

The present investigation elucidates the effects of substrate curvature on the durability of Air Plasma Sprayed (APS) Thermal Barrier Coatings (TBCs). Traditionally, planar disk specimens are utilized in Furnace Cycle Testing (FCT). In most cases, delamination is initiated at the disk's free edge and then propagates along the TBC-bond coat interface. However, in turbine components (e.g., blades/vanes), they lack significant free edges and the coatings are deposited on non-flat surfaces of varying radii of curvature. These geometrical discrepancies imply significant differences in the stress states as compared to planar disk specimens. Therefore, the part geometry inevitably affects the TBC failure mechanisms in the turbine system, and consequently influences TBC durability. Nevertheless, these effects have not been fully explored in the past, despite the anecdotal knowledge which suggests the leading edge of turbine blades to be one of the most common failure/spallation locations. Here in this study, representative TBC systems were deposited onto superalloy disks with flat substrate and rods with curved substrate. The experimental results from FCT suggest TBC durability on rods is significantly lower than when sprayed on disks. Furthermore, the durability of TBCs on rods appears to increase with higher porosity, which is contradictory to reported trends for TBCs on disks. In addition, the effects of various bond coats are not consistent with those observed in disks. To clarify the underlying effects of curved geometries, detailed stress analyses were performed. They revealed a state of thermal stresses which is unique to curved geometries. It was found that TBC on rods experience tensile radial stresses as well as tensile hoop stresses during the cooling. These stresses result in a complex interplay in failure processes in rods, which explains the different observed trends in the durability from those obtained with disks.

36 MATERIALS SCIENCE↗

Command of active and responsive elastomers by topological defects and patterns

The project resulted in the development of stimuli-responsive elastomer coatings formed by photopolymerized liquid crystal molecules. The molecular orientation of the liquid crystal elastomers is coupled to rubber-like elasticity. The project developed an approach to produce complex patterns of molecular orientation in elastomers by employing photopatterning technique based on plasmonic metamasks that convert unpolarized incoming light into a transmitted light beam with spatially-varying linear polarization. The polarization-modulated light beam aligns the film of azodye molecules, which serves as the substrate for liquid crystal elastomer coating. The alignment of the azodye molecules follows the polarization pattern with a phase shift of the alignment direction by 90 degrees. The azodye substrate imposes orientation of molecules in the adjacent liquid crystal layer which is preserved after the monomers are polymerized into a liquid crystal elastomer. By using differently aligned substrates one could produce various patterns of molecular orientation in the elastomer, which dictate the mechanical properties of the elastomers and their response to external stimuli such as temperature, humidity, and ultraviolet irradiation. In particular, the predesigned molecular orientation produces deterministically defined topography response of the liquid crystal coatings, in which their free surface changes from flat to either locally elevated or locally compressed, depending on whether the pre-inscribed molecular orientation is predesigned with a splay or bend deformation. The research established that the mechanism of this deterministic relation between the dynamic topography and pre-inscribed molecular orientation is rooted in the change of the tensorial order parameter that characterizes the degree of orientational order, and in the contraction/expansion of polymer molecules in response to the order parameter changes. In particular, lowering the order causes contraction of polymer globules. The work demonstrated that the dynamic topography of elastomer coatings with three-dimensional pattern of orientational order includes up and down motion of the material (along the normal to the coating) and also shifts in lateral directions. We demonstrate that the deterministic relationship between the orientational patterns and variations of coatings’ profile is caused by forces that are mathematically equivalent to active forces in the system of “pullers” or “pushers” in active matter, demonstrating universality of the out-of-equilibrium behavior. The project resulted in theoretical models that predictively describe the dynamic response of coatings. The topography of liquid crystal elastomers is sensitive to stimuli such as temperature, ultraviolet irradiation, and humidity. The applicability of the patterned approach to produce dynamic coatings was expanded from nematic elastomers to their smectic counterparts. The project combined a battery of imaging techniques to unveil the dynamic properties of the coatings, ranging from the state-of-the-art dynamic holographic microscopy to fluorescent confocal polarizing microscopy and polscope microscopy. The research uncovered potential applications of dynamic elastomer coatings, such as control placement of colloidal particles, engineering of living tissues with alignment patterns of cells that follow the orientational order of elastomer coatings. The approach developed in the project can be used to design programmable dynamic coatings with functionalities that mimic biological tissues such as skin.

36 MATERIALS SCIENCE↗

Simultaneous x-ray and neutron 4D tomographic study of drying-driven hydro-mechanical behavior of cement-based materials at moderate temperatures

The drying of cement-based materials is intimately related to their durability, which has significant economic, social and environmental repercussions. The evolution of the saturation of the pore network and the associated drying shrinkage are in fact leading causes of cracking and of the ensuing penetration of aggressive chemicals. This process is highly heterogeneous, due to the thermo-hydric spatial gradients developing in the material from the exposed surfaces to its core and because of local effects, driven by the intrinsically heterogeneous micro-structure (e.g., by the distribution of pores and aggregates). It follows that macroscopic, sample-scale measurements cannot fully disclose the complexity of the underlying processes. In the last few decades, significant advances in full-field techniques have allowed an unprecedented insight into these local processes. For cement-based materials, x-ray and neutron tomography lend themselves as ideal, and highly complementary, tools for the study of their thermo-hydro-mechanical behavior. Notably, the high sensitivity to density variations of x-ray imaging gives access to the developments of fractures, in 4D (3D + time). On the other hand, neutron tomography allows the study of the evolution of the moisture field in 4D, thanks to its high hydrogen sensitivity. The combination of these two techniques provides a unique insight in thermo-hydro-mechanical couplings, e.g., the effect of cracks on the water content field. This contribution presents novel 5D datasets (3D tomographies along time, plus truly simultaneous x-ray and Neutron rapid acquisitions) in-operando of a cement paste and of a concrete sample heated at moderate temperatures (up to 140{sup ∘}C). The analysis of this 5D data-set (once aligned in time and across modalities) allows for example a quantification of the 4D moisture profiles which were found to predict an overall water loss at hydric equilibrium coherent with the corresponding analytical analysis. In the cement paste sample, the x-ray dataset captures the evolution of an extensive cracking network, opening and propagation toward the core of the sample. A novel analysis procedure is here proposed which allows the extraction of these fractures and the analysis of their interplay with local drying as captured through neutron imaging. This for example reveals the depth of penetration of drying in the vicinity of the fractures along time, which is essential for the assessment and calibration of hydro-mechanical coupled models.

36 MATERIALS SCIENCE↗