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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 415 records · Page 23

All sources NASA NTRS DOE OSTI

DOE OSTI2022

High-Precision Characterization of Seismicity from the 2022 Hunga Tonga-Hunga Ha'apai Volcanic Eruption

The earthquake swarm accompanying the January 2022 Hunga Tonga-Hunga Ha'apai (HTHH) volcanic eruption includes a large number of posteruptive moderate-magnitude seismic events and presents a unique opportunity to use remote monitoring methods to characterize and compare seismic activity with other historical caldera-forming eruptions. We compute improved epicentroid locations, magnitudes, and regional moment tensors of seismic events from this earthquake swarm using regional to teleseismic surface-wave cross correlation and waveform modeling. Precise relative locations of 91 seismic events derived from 59,047 intermediate-period Rayleigh- and Love-wave cross-correlation measurements collapse into a small area surrounding the volcano and exhibit a southeastern time-dependent migration. Regional moment tensors and observed waveforms indicate that these events have a similar mechanism and exhibit a strong positive compensated linear vector dipole component. Precise relative magnitudes agree with regional moment tensor moment magnitude ($M_w$) estimates while also showing that event sizes and frequency increase during the days after the eruption followed by a period of several weeks of less frequent seismicity of a similar size. The combined information from visual observation and early geologic models indicate that the observed seismicity may be the result of a complex series of events that occurred after the explosive eruption on 15 January, possibly involving rapid resupply of the magma chamber shortly after the eruption and additional faulting and instability in the following weeks. In addition, we identify and characterize an $M_w$ 4.5 event five days before the paroxysmal explosion on 15 January, indicating that additional seismic events preceding the main eruption could have been identified with improved local monitoring. As a result, our analysis of the HTHH eruption sequence demonstrates the value of potentially utilizing teleseismic surface-wave cross correlation and waveform modeling methods to assist in the detailed analysis of remote volcanic eruption sequences.

58 GEOSCIENCES↗

DOE OSTI2020

Transient Weighted Moving-Average Model of Photovoltaic Module Back-Surface Temperature

Accurate modeling of photovoltaic (PV) performance requires the precise calculation of module temperature. Currently, most temperature models rely on steady-state assumptions that do not account for the transient climatic conditions and thermal mass of the module. On the other hand, complex physics-based transient models are computationally expensive and difficult to parameterize. In order to address this, a new approach to transient thermal modeling was developed, in which the steady-state predictions from previous timesteps are weighted and averaged to accurately predict the module temperature at finer time scales. This model is informed by 3-D finite-element analyses, which are used to calculate the effect of wind speed and module unit mass on module temperature. The model, in application, serves as an added filter over existing steady-state models that smooths out erroneous values that are a result of intermittency in solar resource. Validation of this moving-average model has shown that it can improve the overall PV energy performance model accuracy by as much as 0.58% over steady-state models based on mean absolute error improvements and can significantly reduce the variability between the model predictions and measured temperature times series data.

14 SOLAR ENERGY↗

DOE OSTI2025

dynamics of organic-mineral interactions at the metal oxide-solution interface as studied via binding energetics (Final report)

This project focused on addressing longstanding fundamental and experimental uncertainties on how dissolved organic substances (DOS) interact with metal oxide surface under environmentally relevant conditions. By leveraging a custom-built real-time, in-tandem flow adsorption microcalorimetry-UV-Vis/fluorescence spectroscopy platform, we characterized the binding energetics, kinetics and mechanistic pathways driving DOS-metal oxide interactions at temporal resolution on the order of 1-5 seconds. We studied a diverse suite of model organic compounds/substances – including monocarboxylates (e.g. acetate and benzoate), di-carboxylates (oxalate and succinate), amino acids, amino-based nanparticles and natural organic matter – interacting at the mineral-water interface of structurally- and/or chemically distinct metal oxides (including SiO2, boehmite, ferrihydrite, and γ-Al2O3). Our results indicated that DOS-metal oxide interactions are governed by multi-step reaction pathways, often switching between distinct, resolvable enthalpy- and entropy-driven non-electrostatic or electrostatic configurations. To quantify these interactions, we developed and implemented an analytical workflow that integrates peak deconvolution and Monte-Carlo based error propagation to determine site-specific thermodynamic and kinetic parameters for individual binding/debinding events. In addition to resolving apparent first-order rate constants of each event, we were able to quantify associated apparent equilibrium constants as well as free energy, enthalpy and entropy contribution to the activation and subsequent progression of the binding/debinding process across compounds, compound class and metal oxide surfaces. The kinetic-thermodynamic data produced in this study captured how the interplay between oxide surface reactivity and DOS molecular structure jointly drives binding-debinding dynamics. Notably, that at pH below PZC of the oxide surface, neutral species were heavily involved in monocarboxylate binding, while anionic species drove dicarboxylate binding. Also, that among amino acids 1) positional isomers show distinctive binding characteristics to each other while enantiomers show no significant differences in binding characteristics, 2) molecules that bind via outer-sphere complexation show a larger entropic shift between binding and debinding with no impact on oxide surface while 3) inner-sphere interactions increased anion exchange capacity of the oxide surface. The new insights and data from this work has great potential for improving predictive modeling of carbon dynamics and specifically organic-mineral interactions in environmental and industrial systems.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Urban Land Surface Effects on Summertime Clouds and Moist Convection in Houston Under Different Synoptic Conditions

Urban landscapes modify cloud formation and convection through complex thermodynamic and aerodynamic processes; however, their influence under different synoptic regimes remains poorly understood. This study investigates the impact of the Houston metropolitan area on summertime cloud cover and convective cell characteristics using a combination of satellite observations, radar data, and high-resolution process-based modeling. We isolate urban effects by comparing model simulations with realistic urban land cover against hypothetical scenarios where all urban areas are replaced by rural vegetation. Results reveal that Houston's land cover consistently enhances cloud fraction and convective activity relative to surrounding rural areas, altered by large-scale meteorological forcing. Under weakly forced conditions, enhanced surface heat flux primarily contributes to driving low-level convergence and vertical ascent, leading to over 8% cloud fraction increase between 2 and 6 km over the city. Under strongly forced conditions, urban influences manifest differently depending on cloud type. For non-convective clouds, the city acts as a barrier that decelerates and lifts moist southerly inflow, increasing low cloud cover over the urban core, while decreasing it downwind the city. For convective clouds, both synoptic ascent and urbanization modulate moisture redistribution and cloud structure, producing modest cloud enhancement over the city and slight suppression over the downwind area. Urbanization exerts small changes in the intensity of the convective cells; however, it significantly decreases their duration and traveling distance. This work highlights the importance of accounting for land surface heterogeneity in modeling clouds and precipitation and demonstrates that urban impacts on clouds are highly regime-dependent.

Liu, Ye [Pacific Northwest National Laboratory (PN↗

DOE OSTI2020

Coexistence of Surface Ferromagnetism and a Gapless Topological State in MnBi 2 Te 4

Surface magnetism and its correlation with the electronic structure are critical to understanding the topological surface state in the intrinsic magnetic topological insulator MnBi 2 Te 4 . Here, using static and time resolved angle-resolved photoemission spectroscopy (ARPES), we find a significant ARPES intensity change together with a gap opening on a Rashba-like conduction band. Comparison with a model simulation strongly indicates that the surface magnetism on cleaved MnBi 2 Te 4 is the same as its bulk state. Finally, the inability of surface ferromagnetism to open a gap in the topological surface state uncovers the novel complexity of MnBi 2 Te 4 that may be responsible for the low quantum anomalous Hall temperature of exfoliated MnBi 2 Te 4 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

Partitioning of Complex Fluids at Mineral Surfaces

This report summarizes the results obtained during the LDRD project entitled "Partitioning of Complex Fluids at Mineral Interfaces." This research addressed fundamental aspects of such interfaces, which are relevant to energy-water applications in the subsurface, including fossil energy extraction and carbon sequestration. This project directly addresses the problem of selectivity of complex fluid components at mineral-fluid interfaces, where complex fluids are defined as a mixture of hydrophobic and hydrophilic components: e.g., water, aqueous ions, polar/nonpolar organic compounds. Specifically, this project investigates how adsorption selectivity varies with surface properties and fluid composition. Both experimental and molecular modeling techniques were used to better understand trends in surface wettability on mineral surfaces. The experimental techniques spanned the macroscale (contact angle measurements) to the nanoscale (cryogenic electronic microscopy and vibrational spectroscopy). We focused on an anionic surfactant and a well-characterized mineral phase representative of clay phases present in oil- and gas-producing shale deposits. Collectively, the results consistently demonstrate that the presence of surfactant in the aqueous fluid significantly affects the mineral-fluid interfacial structure. Experimental and molecular modeling results reveal details of the surfactant structure at the interface, and how this structure varies with surfactant coverage and fluid composition.

04 OIL SHALES AND TAR SANDS↗

DOE OSTI2021

Assessment of a tomography-informed polyhedral discrete element modelling approach for complex-shaped granular woody biomass in stress consolidation

The design of handling equipment for granular biomass primarily requires experiments. Discrete element models (DEM) can provide designers with detailed insight into the behaviour of granular materials. However, granular biomass comprising complex-shaped particles is difficult to model with DEM. A tomography-informed DEM approach and an exhaustive assessment of the ability of the approach to predict the bulk behaviour of milled pines using experimental data is presented. Nano-CT scan was conducted to obtain 3D particle surface geometries as the basis for particle shape approximation by a polyhedral model and a sphero-polyhedral model. These models were applied in the simulation of a compressibility test. Our parametric study showed that particle Young's modulus and restitution coefficient are the two main properties influencing the simulated bulk behaviour of the DEM particles and that the level of approximating DEM particle shape for real particles is critical for accurately replicating the bulk behaviour of milled pines. The polyhedral model demonstrated better suitability than its sphero-polyhedral counterpart for modelling pine particles. The polyhedral model calibrated in a compressibility test was then applied in the simulations of a friction test without additional parameter tuning. Interlocking was dominant in the shear of bulk pine particles. Remarkably, the polyhedral model predicted the frictional behaviour of the pine particles when compared to the experimental results. The limitations of the model, as well as possible ways for enhancement, are discussed. This work provided novel insights into the suitability of complex-shaped DEM models for granular woody biomass.

Biomass flow↗

DOE OSTI2021

An FEM-Based Peridynamic Model for Failure Analysis of Unidirectional Fiber-Reinforced Laminates

To predict the mixed damage modes of unidirectional fiber-reinforced polymer (FRP) laminates under dynamic loading, an FEM-based peridynamic model is introduced in this paper. Based on its geometric structure and material composition, a long fiber lamina is considered a transversely isotropic medium as a result of homogenization at the meso-scale. The laminated structure is modeled by stacking surface mesh layers with arbitrary fiber angles along the thickness direction. The peridynamic bonds between Gauss points connect the separated elements. These bonds are classified as inner-layer bonds and inter-layer bonds. To represent the anisotropy of a laminate, the micro-elastic modulus of the inner-layer and inter-layer bonds is calculated from the anisotropic engineering material constants separately. To capture complex failure behaviors of laminate structures, an empirical damage model is proposed for the tension/ compression breakage of peridynamic bonds. This damage model can control the in-plane and delamination failure process. Finally, benchmark tests are conducted to validate the elastic response of laminates under dynamic loading. In terms of damage analysis, the proposed model can capture the complex damage modes and resistive force of laminate structures.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Mechanical Actuation via Homeomorphic Transformations of Topological Solitons within Polymer Coatings

Abstract Topological solitons are currently under investigation for their exotic properties, especially in nonlinear physics, optics, and material sciences. However, challenges of robust generation and limited stability over time have hindered their practical uses. To address this issue, an approach is developed to form structured arrays of solitons in films of polymerizable liquid crystals. Their complex molecular architecture is preserved by in situ photopolymerization forming a stable liquid crystal network. Most excitingly, their properties are advanced to include responsiveness functions. When thermally actuated, these topological solitons mediate the reconfiguration of surface topographies. Complex shape changes occur depending on the intrinsic complex spatial distribution of the director, which may even lead to full shape inversion and topographical changes as high as ≈40% of the initial thickness. Conversely, the shape changes provide information on the initial director profile, which is consistent with the mathematical model. The soliton‐containing polymer coatings are applicable in multiple domains, ranging from tunable optics to haptics, and from shape‐coupled sensing systems to temperature‐coupled heat management.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Mechanisms and kinetics of citrate-promoted dissolution of a uranyl phosphate mineral

The formation of uranyl phosphate precipitates is a remediation strategy to limit the mobility of uranium in contaminated soils. However, exposure to organic ligands, like the plant exudate citrate, can remobilize the uranium. The purpose of this study is to provide a more thorough comprehension of citrate-promoted dissolution of a uranyl phosphate mineral, chernikovite [(H 3 O)(UO 2 )(PO 4 )·3H 2 O], by determining the extent of uranium release from chernikovite at a wide range of citrate concentrations. Here we have quantified the kinetics of dissolution and proposed potential mechanisms of chernikovite dissolution to gain a better understanding of the fate of uranyl phosphate precipitates in the environment. Batch dissolution and continuously stirred tank reactor (CSTR) experiments indicate that increasing citrate concentrations from 0.1 mM to 50 mM increases the concentration of dissolved uranium. However, at citrate concentrations of 10 mM and greater, the effectiveness of the ligand to enhance the dissolution decreases. Density functional theory modeling and Raman spectroscopy indicate a strong interaction between citrate and the uranyl ions at the surface of the uranyl phosphate. The interaction between citrate and uranyl suggests the formation of an alteration layer of uranyl-citrate surface complexes, which could impede uranium release at citrate concentrations of 10 mM or higher. Together, these results indicate that the citrate-promoted dissolution of chernikovite is a fast reaction that is hindered by a combination of surface saturation and a secondary-phase precipitation reaction at higher concentrations of citrate and reaction times greater than or equal to 72 min based on the hydraulic residence time of CSTR experiments. These findings improve our understanding of the stability of uranyl phosphate for use in sequestration of uranium groundwater plumes at contaminated sites.

58 GEOSCIENCES↗

DOE OSTI2023

Unraveling the 2021 Central Tennessee flood event using a hierarchical multi-model inundation modeling framework

Flood prediction systems need hierarchical atmospheric, hydrologic, and hydraulic models to predict rainfall, runoff, streamflow, and floodplain inundation. The accuracy of such systems depends on the error propagation through the modeling chain, sensitivity to input data, and choice of models. In this study, we used multiple precipitation forcings (hindcast and forecast) to drive hydrologic and hydrodynamic models to analyze the impacts of various drivers on the estimates of flood inundation depth and extent. We implement this framework to unravel the August 2021 extreme flooding event that occurred in Central Tennessee, USA. We used two radar-based quantitative precipitation estimates (STAGE4 and MRMS) as well as quantitative precipitation forecasts (QPF) from the National Weather Service Weather Prediction Center (WPC) to drive a series of models in the hierarchical framework, including the Variable Infiltration Capacity (VIC) land surface model, the Routing Application for Parallel Computation of Discharge (RAPID) river routing model, and the AutoRoute and TRITON inundation models. An evaluation with observed high-water marks demonstrates that the framework can reasonably simulate flood inundation. Despite the complex error propagation mechanism of the modeling chain, we show that inundation estimates are most sensitive to rainfall estimates. Most notably, QPF significantly underestimates flood magnitudes and inundations leading to unanticipated severe flooding for all stakeholders involved in the event. Finally, we discuss the implications of the hydrodynamic modeling framework for real-time flood forecasting.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Observing the three-dimensional dynamics of supported metal complexes

Dynamics are intricately linked with activity and selectivity when it comes to catalysis, as noted for instance in the enzymatic principles of induced fit and allostery, and yet the range of motions heterogeneous catalytic sites are able to undergo is poorly understood. Solid-state nuclear magnetic resonance (NMR) spectroscopy is perhaps the only tool capable of probing the rapid conformational dynamics found in heterogeneous catalysts but has historically been restricted by its low sensitivity, limiting the detail with which structures can be resolved. In this work, we apply solid-state NMR and dynamic nuclear polarization, in combination with density functional theory modeling, to reveal the high-resolution structure and motional freedom of a scandium supported complex in three dimensions. The results are contrasted with the study of the analogous homoleptic complex in the crystalline state, highlighting the impacts that surface structure may have on the dynamics of supported complexes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Ion emission from warm dense matter produced by irradiation with a soft x-ray free-electron laser

Here we report on an experiment performed at the FLASH2 free-electron laser (FEL) aimed at producing warm dense matter via soft x-ray isochoric heating. In the experiment, we focus on study of the ions emitted during the soft x-ray ablation process using time-of-flight electron multipliers and a shifted Maxwell–Boltzmann velocity distribution model. We find that most emitted ions are thermal, but that some impurities chemisorbed on the target surface, such as protons, are accelerated by the electrostatic field created in the plasma by escaped electrons. The morphology of the complex crater structure indicates the presence of several ion groups with varying temperatures. We find that the ion sound velocity is controlled by the ion temperature and show how the ion yield depends on the FEL radiation attenuation length in different materials.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

First-principles design of a single-atom–alloy propane dehydrogenation catalyst

The complexity of heterogeneous catalysts means that a priori design of new catalytic materials is difficult, but the well-defined nature of single-atom–alloy catalysts has made it feasible to perform unambiguous theoretical modeling and precise surface science experiments. Herein we report the theory-led discovery of a rhodium-copper (RhCu) single-atom–alloy catalyst for propane dehydrogenation to propene. Although Rh is not generally considered for alkane dehydrogenation, first-principles calculations revealed that Rh atoms disperse in Cu and exhibit low carbon-hydrogen bond activation barriers. Surface science experiments confirmed these predictions, and together these results informed the design of a highly active, selective, and coke-resistant RhCu nanoparticle catalyst that enables low-temperature nonoxidative propane dehydrogenation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Scalable Hot-Water-Repellent Superhydrophobicity via Thermal Insulation

Superhydrophobic surfaces, which rely on a combination of surface texture and chemistry, often lose their repellent behavior when contacted by hot water (≳40 °C) because the impinging hot water replaces the requisite air layer within the surface texture via evaporation and recondensation. In contrast to previous approaches targeting this condensation-induced failure mode that rely on intricately tailored surface structures or complex chemical treatments, we present a scalable approach based on thermal design: the multilayered insulated superhydrophobic (MISH) coating mitigates condensation-induced failure by preventing heat transfer. Superhydrophobicity is retained at impinging water temperatures up to 90 °C, with durability demonstrated via long-term (>1 million impacts) droplet impingement experiments. We explain the mechanism for this approach with a detailed thermal model; the model reveals that the underlying physical behavior is self-similar across coating parameters and impinging fluid temperatures. Finally, the MISH coating accommodates curved geometries and large surfaces, and it is over 4 orders of magnitude less expensive than cleanroomnanofabricated alternatives, indicating promise for practical use in the energy sector, chemical processing, and the food and medical industries.

coating materials↗

DOE OSTI2021

Interface creation on a mixed-terminated perovskite surface

In the field of complex oxide heterostructures, understanding of the initial substrate surface can be critical to fundamental studies regarding the development of emergent properties at the film–substrate interface. For this reason, a considerable amount of effort has gone into the development of techniques to achieve surfaces with single termination for a variety of perovskite single crystals. However, a decisive understanding of how an interface is created when epitaxial growth occurs on a mixed terminated surface remains lacking. Employing in situ synchrotron X-ray scattering during thin film growth by molecular beam epitaxy, we investigate the initial stages of growth on (LaAlO 3 )(Sr 2 AlTaO 6 ) (001) substrates with mixed termination. Using LaNiO 3 as a model system, we find that the surface layers of the substrate restructure during deposition such that while a NiO 2 monolayer is weakly bound, a LaO monolayer bonds strongly, effectively incorporating with (Al, Ta)O 2 from the surface and forming a La(Al, Ta)O 3 ultrathin film.  With regard to the synthesis of high-quality complex oxide heterostructures, it is important that the substrate meet an array of criteria related to the crystal structure, lattice parameter, thermodynamic properties, electrical and optical properties, and so on. In particular, the surface of the substrate should be atomically smooth and exhibit the step-terrace structure ideal for epitaxial thin film growth. Perovskite oxides with the formula ABO 3 have two possible (001) surfaces: AO (A-site termination) and BO 2 (B-site termination). The ability to consistently achieve one of these surface terminations on SrTiO 3 , i.e., TiO 2 , has helped SrTiO 3 become the predominant oxide in the field.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Plasmon coupling in topological insulator multilayers

Topological insulators (TIs) house two-dimensional Dirac plasmons on their surfaces. In TI thin films, plasmons on the top and bottom surfaces couple electrostatically, giving rise to an acoustic and an optical plasmon mode. By extension, a superlattice comprising TI layers and trivially-insulating layers could house multiple complex plasmon modes which could be used to create a new type of Dirac metamaterial. In this paper, we synthesize TI superlattices, fabricate them into stripe arrays, characterize their optical and plasmonic properties, and model the samples using transfer matrices. Furthermore, we excite plasmon modes that couple to the phonons in the superlattice, resulting in hybrid plasmon phonon polaritons. These modes are modeled using an analytical Fano resonance model and the extracted resonant positions are reproduced by transfer matrix modeling. We also excite an epsilon-near-zero mode in the top dielectric material (Bi 0.5 In 0.5 ) 2 Se 3 . Understanding the behavior of the polariton modes in this complex system will lend insight into the many-body interaction in low-dimensional systems.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Calculation of the Heat Transfer Coefficient in the Outer Body for a Rotational Detonation

Unsteady heat transfer characterization on the combustion surfaces of Rotational Detonation Engines (RDE) is not well understood. It is generally thought that the complex nature of the unsteady, reacting, compressible fluid flow inside the combustion anulus of the RDE causes the convective heat transfer coefficient to be significantly higher than it is in other applications. Empirical models that have been used to analyze this strictly apply to steady flows where dimensionless groups can be employed. The reacting flow fields that are characteristic of RDEs are inherently compressible, three dimensional, unsteady, and turbulent, having properties that change by orders of magnitude throughout the flow field. They will therefore contain multiple length scales and time scales operating everywhere in the flow during all times. Because of this it is not likely that the RDE flow fields will lend themselves to explanation using simple dimensionless parameters. The dimensionless groups have meaning only in situations where length scales and time scales are singular and well defined. In spite of this it may be possible to get a relatively good idea of what the convection heat transfer coefficient is. In this work a numerical study is performed where the inside wall surface temperature distribution in the RDE outer body is systematically changed over a given range that would be characteristic of the start-up flows inside an RDE. For each case, temperature distributions inside the outer containment wall of the RDE was calculated and compared with experimental data. The closet match can then be used to directly calculate the convection heat transfer coefficient on the inside surface of the RDE.

VanOsdol, John↗