Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “surface complexation model”

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 343 records · Page 19

A computational fluid dynamics model to estimate local quantities in firebrand char oxidation

Firebrand burning is a complex phenomenon that is influenced by several parameters which are difficult to fully explore experimentally. Computational fluid dynamics models capable of predicting local quantities are essential for accurate prediction of char oxidation in firebrands. This article presents a computational fluid dynamics model to estimate firebrand mass loss, diameter change, and surface temperature during char oxidation. The model was validated using previously conducted wind tunnel experiments. These experiments were conducted for firebrands of two different aspect ratios, which were arranged in three different configurations (single, horizontal array, and vertical array), and for four different wind speeds (0.5, 1, 1.5, and 2 m/s). The computational fluid dynamics results were compared with a previous 1 D model. In all the test cases, the computational fluid dynamics model predicted the physical phenomena with significantly improved accuracy compared to a 1 D model. The char oxidation model presented in this article can be coupled with other models to study firebrand generation and trajectory, biomass pyrolysis, fluidized bed reactors, and coal combustion.

Engineering↗

From Micro-environments to Macroscopic Effects: How the Alkaline Hydrogen Evolution Reaction Drives Cu Cathodic Corrosion

Cathodic corrosion of copper (Cu) has posed a significant challenge for over a century, impeding various technological progresses such as electrochemical conversion of CO 2 (eCO 2 RR) into fuels and other value-added carbon products. Here, in this study, employing a combined Density Functional Theory (DFT) and kinetic Monte Carlo (kMC) simulation approach, we delve into the atomistic level mechanism driving this phenomenon in Cu. Our hypothesis posits the pivotal role of alkaline hydrogen evolution reaction (HER) in facilitating cathodic corrosion in Cu. We rigorously develop a pH-dependent hydroxide (OH) adsorption mechanism and calculate the equilibrium OH coverage (𝜃 OH ) at varying pH levels, the thermodynamic stability of subsurface oxygen (O sub ), as well as the Cu-vacancy mediated diffusion of subsurface oxygens (O sub ). Through comprehensive analysis, we establish correlation among various microenvironments, including oxygen diffusion in subsurface layers, pH-dependent OH adsorption, and Cu dissolution into the electrolyte as (Cu-OH) complexes. Furthermore, our investigation explores the correlation between surface coordination environment of active sites and cathodic corrosion of Cu. Finally, by integrating DFT-derived thermodynamic data into a kMC model, we successfully predict the formation of experimentally observed corrosion pits on Cu-surfaces. This combined approach not only advances our fundamental understanding of Cu cathodic corrosion but also offers insights crucial for developing effective corrosion mitigation strategies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Intrinsic 150‐Day Periodicity of the Southern Hemisphere Extratropical Large‐Scale Atmospheric Circulation

Abstract The variability of the Southern Hemisphere (SH) extratropical large‐scale circulation is dominated by the Southern Annular Mode (SAM), whose timescale is extensively used as a key metric in evaluating state‐of‐the‐art climate models. Past observational and theoretical studies suggest that the SAM lacks any internally generated (intrinsic) periodicity. Here, we show, using observations and a climate model hierarchy, that the SAM has an intrinsic 150‐day periodicity. This periodicity is robustly detectable in the power spectra and principal oscillation patterns (aka dynamical mode decomposition) of the zonal‐mean circulation, and in hemispheric‐scale precipitation and ocean surface wind stress. The 150‐day period is consistent with the predictions of a new reduced‐order model for the SAM, which suggests that this periodicity is associated with a complex interaction of turbulent eddies and zonal wind anomalies, as the latter propagate from low to high latitudes. These findings present a rare example of periodic oscillations arising from the internal dynamics of the extratropical turbulent circulations. Based on these findings, we further propose a new metric for evaluating climate models, and show that some of the previously reported shortcomings and improvements in simulating SAM's variability connect to the models' ability in reproducing this periodicity. We argue that this periodicity should be considered in evaluating climate models and understanding the past, current, and projected Southern Hemisphere climate variability.

54 ENVIRONMENTAL SCIENCES↗

New experimental approach to understanding the chemical reactivity of oxide surfaces

Metal oxides have been an attractive option for a range of applications, including hydrogen sensors, microelectronics, and catalysis, due to their reactivity and tunability. The properties of metal oxides can vary greatly on their precise surface structure; however, few surface science techniques can achieve atomistic-level determinations of surface structure, and fewer yet can do so for insulator surfaces. Low energy ion beam analysis offers a potential insulator-compatible solution to characterizing the surface structure of metal oxides. As a feasibility study, we apply low energy ion beam analysis to investigate the surface structure of a magnetite single crystal, Fe 3 O 4 (100). We obtain multi-angle maps using both forward-scattering low energy ion scattering (LEIS) and backscattering impact-collision ion scattering spectroscopy (ICISS). Both sets of experimental maps have intensity patterns that reflect the symmetries of the Fe 3 O 4 (100) surface structure. However, analytical interpretation of these intensity patterns to extract details of the surface structure is significantly more complex than previous LEIS and ICISS structural studies of one-component metal crystals, which had far more symmetries to exploit. To gain further insight into the surface structure, we model our experimental measurements with ion-trajectory tracing simulations using molecular dynamics. Our simulations provide a qualitative indication that our experimental measurements agree better with a subsurface cation vacancy model than a distorted bulk model.

36 MATERIALS SCIENCE↗

GEOMWRITER - A Workflow for Automatically Generating MCNP Geometry for Neutron Guide Simulations

This report describes the functionality of the geomwriter.py script, which is a python script for producing complex MCNP CSG geometries from the IGES CAD format. It is specifically written to handle neutron supermirror guides, but can handle more complex objects through utilizing a convex hull algorithm. This work flow reduces the amount of time spent on generating the surface definitions for a complex neutron guide system with hundreds of segments such that changes to the neutron guide system design in McStas can be rapidly incorporated into an MCNP model and integrated with the shielding geometry present in that model.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Urban weather modeling using WRF: linking physical assumptions, code implementation, and observational needs

The Weather Research and Forecasting (WRF) model includes urban schemes that simulate the influence of urban surfaces on the atmosphere using parameterizations for flux, and radiative exchanges. Three core schemes – the Bulk urban parameterization, Single-Layer Urban Canopy Model (SLUCM), and Multi-Layer Urban Canopy Model (MLUCM) – represent increasing levels of complexity. Although the parameterizations within these urban schemes are described in the literature, their specific implementation remains poorly documented, thus slowing down model development efforts. This manuscript provides a roadmap to the three urban schemes in WRF version 4.5.2, presenting equations using the same symbols as in the model code, along with references to code lines, and including graphics and explanations that connect the code to its physical foundations. Our thorough review of the urban parameterizations implemented in WRF version 4.5.2 highlighted a handful of parameters that may introduce discontinuities in simulations: (i) in the SLUCM, a 1 mm h −1 rain rate threshold is employed to switch between two minimum moisture availability parameterizations, thus impacting latent heat flux calculations; (ii) in the SLUCM a threshold is used to partition shortwave radiation into direct and diffuse components; (iii) in all three urban schemes, the bulk Richardson number is employed to select the similarity function, which influences the vertical distribution of heat and momentum. We also identified a highly simplified treatment of the radiative balance on roof surfaces. The implications of these simplifications can be assessed through targeted observations across relevant conditions, including varying precipitation rates, cloud cover, and transitions between stability regimes. Furthermore, the widespread application of the Monin-Obukhov similarity theory in these urban schemes warrants model evaluation under highly stable and unstable conditions and in heterogeneous urban settings with variable land cover and building heights on scales finer than model resolution. To address these challenges, we offer guidance on observational strategies, emphasizing the need for multi-parametric measurements to capture potential compensating biases and multi-height measurements that align with the levels where quantities are diagnostic and prognosed in the model (i.e., the lowest atmospheric level of the WRF model). Finally, our inspection of the code revealed implementation bugs that have now been corrected in WRF versions 4.6.0 and 4.6.1. Sensitivity tests over the Atlanta urban area show that these corrections affect surface temperatures, underscoring the importance of performing rigorous documentation and verification of the implementation of parameterizations in model code.

54 ENVIRONMENTAL SCIENCES↗

Computing water flow through complex landscapes – Part 2: Finding hierarchies in depressions and morphological segmentations

Depressions – inwardly draining regions of digital elevation models – present difficulties for terrain analysis and hydrological modeling. Analogous “depressions” also arise in image processing and morphological segmentation, where they may represent noise, features of interest, or both. Here we provide a new data structure – the depression hierarchy – that captures the full topologic and topographic complexity of depressions in a region. We treat depressions as networks in a way that is analogous to surface-water flow paths, in which individual sub-depressions merge together to form meta-depressions in a process that continues until they begin to drain externally. This hierarchy can be used to selectively fill or breach depressions or to accelerate dynamic models of hydrological flow. Complete, well-commented, open-source code and correctness tests are available on GitHub and Zenodo.

54 ENVIRONMENTAL SCIENCES↗

Dilute Alloys Based on Au, Ag, or Cu for Efficient Catalysis: From Synthesis to Active Sites

The development of new catalyst materials for energy-efficient chemical synthesis is critical as over 80% of industrial processes rely on catalysts, with many of the most energy-intensive processes specifically using heterogeneous catalysis. Catalytic performance is a complex interplay of phenomena involving temperature, pressure, gas composition, surface composition and structure over multiple length and time scales. In response to this complexity, the integrated approach to heterogeneous dilute-alloy catalysis reviewed here brings together materials synthesis, mechanistic surface chemistry, reaction kinetics, in-situ and operando characterization, and theoretical calculations in a coordinated effort to develop design principles to predict and improve catalytic selectivity. Dilute alloy catalysts—in which isolated atoms or small ensembles of the minority metal on the host metal lead to enhanced reactivity while retaining selectivity—are particularly promising as selective catalysts. Several dilute alloy materials using Au, Ag and Cu as the majority host element, including more recently introduced support-free nanoporous metals and oxide-supported nanoparticle "raspberry colloid templated (RCT)" materials, are reviewed for selective oxidation and hydrogenation reactions. Progress in understanding how such dilute alloy catalysts can be used to enhance selectivity of key synthetic reactions is reviewed, including quantitative scaling from model studies to catalytic conditions. The dynamic evolution of catalyst structure and composition studied in surface science and catalytic conditions and their relationship to catalytic function are also discussed, followed by advanced characterization and theoretical modeling that have been developed to determine the distribution of minority metal atoms at or near the surface. Furthermore, the integrated approach demonstrates the success of bridging the divide between fundamental knowledge and design of catalytic processes in complex catalytic systems, which can accelerate the development of new and efficient catalytic processes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

AN INITIAL LOOK AT THE HEAT PIPE RESILIENCY TO REACTOR OPERATION

Several special purpose reactor designs are aiming to utilize heat pipes due to its inherently passive features allowing safe heat removal to the power conversion systems. There is insufficient data on heat pipe survivability in reactor environments, although heat-pipe failures are predicted to have low failure rates. In this paper, we perform a coupled neutronics and thermo-mechanics analyses for a 45 kWth HALEU fueled, hydride moderated homogenous design to evaluate the transient effects during both startup and a heatpipe-failure scenario. The heatpipe-failure study includes both a single failure case and a failure-propagation (cascading) case. While the model only evaluates a homogenous core, the approach could be applied to more complex models for any heat pipe cooled small reactor.

99 GENERAL AND MISCELLANEOUS↗

Short-term impacts of 2017 western North American wildfires on meteorology, the atmosphere’s energy budget, and premature mortality

Abstract Western North American fires have been increasing in magnitude and severity over the last few decades. The complex coupling of fires with the atmospheric energy budget and meteorology creates short-term feedbacks on regional weather altering the amount of pollution to which Americans are exposed. Using a combination of model simulations and observations, this study shows that the severe fires in the summer of 2017 increased atmospheric aerosol concentrations leading to a cooling of the air at the surface, reductions in sensible heat fluxes, and a lowering of the planetary boundary layer height over land. This combination of lower-boundary layer height and increased aerosol pollution from the fires reduces air quality. We estimate that from start of August to end of October 2017, ∼400 premature deaths occurred within the western US as a result of short-term exposure to elevated PM 2.5 from fire smoke. As North America confronts a warming climate with more fires the short-term climate and pollution impacts of increased fire activity should be assessed within policy aimed to minimize impacts of climate change on society.

Bernstein, Diana N. (ORCID:0000000308279336)↗

Ring Model Development and Validation for Prismatic HTGR Core Thermal-hydraulics and Safety Analysis

Because of the complex core geometry, prismatic high temperature gas-cooled reactors (prismatic HTGRs) often exhibit complex thermal fluid behaviors during both normal operating and transient conditions. Most HTGR designs rely on passive safety system for decay heat removal, such as the reactor cavity cooling system (RCCS). During postulated accidents like Pressurized Conduction Cooldown (PCC) event, the decay heat is first radially transferred from the core region to the reactor vessel outer surface, then to the RCCS cooling panels. The peak fuel temperature is controlled by heat transfer mechanisms with two distinctive characteristic length scales, i.e., the core-wise effective heat conduction and local heat conduction in the fuel pellet scale. As both length scales are essential to determine the fuel temperature, from the modeling perspective, computer codes must be able to capture heat transfer in both scales. This is challenging for both computational fluid dynamics (CFD) tools because of extremely large amount of computation resources required, and for system analysis codes because of the challenge to model the complex core geometry. Under the support of DOE-NE’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, efforts have been pursued to support HTGR technology development and its modeling and simulation needs. There is a particular need for advanced modeling and simulation tools to predict thermal-fluid behavior during safety-related transients. In our previous studies, the ring model was adopted in SAM and further developed to simulate the normal operating condition and a PCC event using the MHTGR-350 design of General Atomics as the reference design. This current work represents a continuation of these previous efforts, and the focus is to critically review and examine simplifications and assumptions made to develop the ring model, and to perform code validation using experimental data from an integral-effect test facility, the High Temperature Test Facility (HTTF) at the Oregon State University. In this study, the test PG-27 from the HTTF test suite was selected for code validation purpose. The test PG-27 is a transient test designed to simulate the PCC event of the MHTGR design. Very good agreements between SAM prediction and experimental measurements were found in both coolant and solid structure temperatures during the transient.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Mechanistic understanding of speciated oxide growth in high entropy alloys

Abstract Complex multi-element alloys are gaining prominence for structural applications, supplementing steels, and superalloys. Understanding the impact of each element on alloy surfaces due to oxidation is vital in maintaining material integrity. This study investigates oxidation mechanisms in these alloys using a model five-element equiatomic CoCrFeNiMn alloy, in a controlled oxygen environment. The oxidation-induced surface changes correlate with each element’s interactive tendencies with the environment, guided by thermodynamics. Initial oxidation stages follow atomic size and redox potential, with the latter becoming dominant over time, causing composition inversion. The study employs in-situ atom probe tomography, transmission electron microscopy, and X-ray absorption near-edge structure techniques to elucidate the oxidation process and surface oxide structure evolution. Our findings deconvolute the mechanism for compositional and structural changes in the oxide film and will pave the way for a predictive design of complex alloys with improved resistance to oxidation under extreme conditions.

36 MATERIALS SCIENCE↗

Connection between the physicochemical characteristics of amorphous carbon thin films and their electrochemical properties

Connecting a material's surface chemistry with its electrocatalytic performance is one of the major questions in analytical electrochemistry. This is especially important in many sensor applications where analytes from complex media need to be measured. Unfortunately, today this connection is still largely missing except perhaps for the most simple ideal model systems. Here we present an approach that can be used to obtain insights about this missing connection and apply it to the case of carbon nanomaterials. In this paper we show that by combining advanced computational techniques augmented by machine learning methods with x-ray absorption spectroscopy (XAS) and electrochemical measurements, it is possible to obtain a deeper understanding of the correlation between local surface chemistry and electrochemical performance. As a test case we show how by computationally assessing the growth of amorphous carbon (a-C) thin films at the atomic level, we can create computational structural motifs that may in turn be used to deconvolute the XAS data from the real samples resulting in local chemical information. Then, by carrying out electrochemical measurements on the same samples from which x-ray spectra were measured and that were further characterized computationally, it is possible to gain insight into the interplay between the local surface chemistry and electrochemical performance. To demonstrate this methodology, we proceed as follows: after assessing the basic electrochemical properties of a-C films, we investigate the effect of short HNO 3 treatment on the sensitivity of these electrodes towards an inner sphere redox probe dopamine to gain knowledge about the influence of altered surface chemistry to observed electrochemical performance. Overall, these results pave the way towards a more general assessment of electrocatalysis in different systems and provide the first steps towards data driven tailoring of electrode surfaces to gain optimal performance in a given application.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Review—Thermal Safety Management in Li-Ion Batteries: Current Issues and Perspectives

Approaches for thermal management of lithium-ion (Li-ion) batteries do not always keep pace with advances in energy storage and power delivering capabilities. Root-cause analysis and empirical evidence indicate that thermal runaway (TR) in cells and cell-to-cell thermal propagation are due to adverse changes in physical and chemical characteristics internal to the cell. However, industry widely uses battery management systems (BMS) originally designed for aqueous-based batteries to manage Li-ion batteries. Even the "best" BMS that monitor both voltage and outside-surface temperature of each cell are not capable of preventing TR or TR propagation, because voltage and surface-mounted temperature sensors do not track fast-emerging adverse events inside a cell. Most BMS typically include a few thermistors mounted on select cells to monitor their surface temperature. Technology to track intra-cell changes that are TR precursors is becoming available. Simultaneously, the complex pathways resulting in cell-to-cell TR propagation are being successfully modelled and mapped. Innovative solutions to prevent TR and thermal propagation are being advanced. These include modern BMS for rapid monitoring the internal health of each individual cell and physical as well as chemical methods to reduce the deleterious effects of rapid cell-to-cell heat and material transport in case of TR.

25 ENERGY STORAGE↗

Benchmarking Near-Surface Winds in the HRRR Analyses Using Multisource Observations over Complex Terrain in the Southeastern United States

Wind energy plays a crucial role in sustainable power generation, yet its full potential in the southeast United States (SEUS) remains underexplored. This study advances wind resource assessment in the SEUS using existing regional modeling and multisource observations. We find that terrain complexity, influenced by orography and forest canopy, significantly impacts the accuracy of modeled wind speed. While the High-Resolution Rapid Refresh (HRRR) model effectively simulates wind profiles over flat, nonforested terrain, larger errors are produced above forest canopies, and particularly for those stands in hilly and mountainous terrain where there are no observations aloft that could be used by HRRR’s data assimilation software. Seasonal variations, e.g., changes in the leaf area index, further complicate the relationship between terrain complexity and simulation errors. Here, our findings emphasize the critical need for comprehensive wind profile measurements, extending from below the canopy to height above the canopy and through the lower planetary boundary layer, to fully quantify model performance over and near forests, and for improving wind resource assessment, planning, and management. This study provides valuable insights not only for further model development but also for future field campaign deployment with the goal of further improving our understanding of wind resources in the region.

Complex terrain↗

Seasonal Cycle of Idealized Polar Clouds: Large Eddy Simulations Driven by a GCM

The uncertainty in polar cloud feedbacks calls for process understanding of the cloud response to climate warming. As an initial step toward improved process understanding, we investigate the seasonal cycle of polar clouds in the current climate by adopting a novel modeling framework using large eddy simulations (LES), which explicitly resolve cloud dynamics. Resolved horizontal and vertical advection of heat and moisture from an idealized GCM are prescribed as forcing in the LES. The LES are also forced with prescribed sea ice thickness, but surface temperature, atmospheric temperature, and moisture evolve freely without nudging. A semigray radiative transfer scheme without water vapor and cloud feedbacks allows the GCM and LES to achieve closed energy budgets more easily than would be possible with more complex schemes. This enables the mean states in the two models to be consistently compared, without the added complications from interaction with more comprehensive radiation. Furthermore, we show that the LES closely follow the GCM seasonal cycle, and the seasonal cycle of low-level clouds in the LES resembles observations: maximum cloud liquid occurs in late summer and early autumn, and winter clouds are dominated by ice in the upper troposphere. Large-scale advection of moisture provides the main source of water vapor for the liquid-containing clouds in summer, while a temperature advection peak in winter makes the atmosphere relatively dry and reduces cloud condensate. The framework we develop and employ can be used broadly for studying cloud processes and the response of polar clouds to climate warming.

54 ENVIRONMENTAL SCIENCES↗

Quantifying Groundwater Response and Uncertainty in Beaver‐Influenced Mountainous Floodplains Using Machine Learning‐Based Model Calibration

Abstract Beavers ( Castor canadensis ) alter river corridor hydrology by creating ponds and inundating floodplains, and thereby improving surface water storage. However, the impact of inundation on groundwater, particularly in mountainous alluvial floodplains with permeable gravel/cobble layers overlain by a soil layer, remains uncertain. Numerical modeling across various floodplain structures considers topographic and sediment complexity and multidirectional flow, linking inundation to groundwater response. This study develops a model‐data integration workflow to address uncertainty in groundwater response to beaver‐induced inundations in a mountainous alluvial floodplain in the Upper Colorado River Basin. Uncertain factors include seasonal hydrologic dynamics, hydraulic conductivities, floodplain structures, and meteorological forcings. We employed an ensemble of groundwater models, based on geophysical and hydrologic data, with machine learning‐based calibration using a neural density estimator. This allowed us to quantify the vertical flux from the soil layer to the permeable gravel bed, the down‐valley underflow within the gravel bed, and their ratios. Results show a significant increase in the vertical flux relative to down‐valley underflow, from 2 during dry pond periods to 20 during wet periods, serving as an analogy for conditions without and with beaver ponds. The study highlights the influence of floodplain structure on groundwater storage, water balance, and water quality impacted by beaver ponds. A thick gravel bed layer, with a large down‐valley underflow, minimizes the effect of beaver‐induced inundation on water quality. We emphasize the need for field‐scale measurements of floodplain structure and improved characterization of evapotranspiration changes to reduce uncertainty in groundwater response. Plain Language Summary Beavers change the flow of water in river corridors by creating ponds, expanding wetlands, and flooding floodplains. This increases surface water area, promotes plant growth, and enhances biodiversity. However, the impact of this flooding on groundwater flow is not well understood, especially in mountainous areas with gravel layers where water moves easily beneath soil. In this study, we used numerical modeling to investigate how beaver ponds influence groundwater in a mountainous floodplain of the Upper Colorado River Basin. We adapted a machine learning method to validate our numerical models using multiple field data sets. Our findings show that beaver ponds significantly increase vertical water flow from the soil to the gravel during wet periods, compared to when the ponds are fully drained. The study also highlights the importance of floodplain structure in controlling both water flow in gravel layers along the river direction and vertical flow from the soil to the gravel with the presence of beavers. To reduce uncertainty in groundwater response, we emphasize the need for more field‐scale measurements of floodplain structure, hydraulic properties, and evapotranspiration changes. Key Points Floodplain structures and hydraulic conductivities are important for groundwater response with beaver ponds in mountainous floodplains Large down‐valley underflow in permeability‐stratified floodplains reduces beaver‐induced impacts on groundwater storage and water quality Machine learning‐based model calibration methods are effective for estimating posterior distributions of groundwater model parameters

Wang, Lijing↗

Competing electronic states emerging on polar surfaces

Excess charge on polar surfaces of ionic compounds is commonly described by the two-dimensional electron gas (2DEG) model, a homogeneous distribution of charge, spatially-confined in a few atomic layers. Here, by combining scanning probe microscopy with density functional theory calculations, we show that excess charge on the polar TaO 2 termination of KTaO 3 (001) forms more complex electronic states with different degrees of spatial and electronic localization: charge density waves (CDW) coexist with strongly-localized electron polarons and bipolarons. These surface electronic reconstructions, originating from the combined action of electron-lattice interaction and electronic correlation, are energetically more favorable than the 2DEG solution. They exhibit distinct spectroscopy signals and impact on the surface properties, as manifested by a local suppression of ferroelectric distortions.

36 MATERIALS SCIENCE↗