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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Roadmap on methods and software for electronic structure based simulations in chemistry and materials

This Roadmap article provides a succinct, comprehensive overview of the state of electronic structure methods and software for molecular and materials simulations. Seventeen distinct sections collect insights by 51 leading scientists in the field. Each contribution addresses the status of a particular area, as well as current challenges and anticipated future advances, with a particular eye towards software related aspects and providing key references for further reading. Foundational sections cover density functional theory and its implementation in real-world simulation frameworks, Green's function based many-body perturbation theory, wave-function based and stochastic electronic structure approaches, relativistic effects and semiempirical electronic structure theory approaches. Subsequent sections cover nuclear quantum effects, real-time propagation of the electronic structure, challenges for computational spectroscopy simulations, and exploration of complex potential energy surfaces. The final sections summarize practical aspects, including computational workflows for complex simulation tasks, the impact of current and future high-performance computing architectures, software engineering practices, education and training to maintain and broaden the community, as well as the status of and needs for electronic structure based modeling from the vantage point of industry environments. Overall, the field of electronic structure software and method development continues to unlock immense opportunities for future scientific discovery, based on the growing ability of computations to reveal complex phenomena, processes and properties that are determined by the make-up of matter at the atomic scale, with high precision.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Thermal Performance of Triply Periodic Minimal Surface Lattice Structures in Single-Phase Dielectric Fluid Cooling of Power Electronics

Additive manufacturing has transformed thermal management by enabling the production of complex, optimized geometries that conventional manufacturing methods cannot achieve. This study investigates the single-phase convective heat transfer performance of gyroid triply periodic minimal surface (TPMS) lattice structures with functional porosity. TPMS structures provide high surface area to volume ratios and are amenable to 3D printing. A gyroid numerical model was created and validated against an existing experimental study with a similar feature size to the investigated geometries. The TPMS structure has a periodic width of 1.6 mm, a length of 10 mm, and a height of 4 mm, with a functional porosity ranging from 0.5 to 0.8, decreasing with distance from the heated surface. Three different flow configurations were examined for an inlet fluid temperature of 70 °C. The inlet velocities range from 0.01 to 1.2 m/s, corresponding to a Reynolds number range of 10–900 with a heat flux of 50 W/cm 2 applied at the base. AmpCool ® AC-110 dielectric fluid (Prandtl number 59.5) was used as the coolant. Thermal performance and friction characteristics were studied for the three flow orientations. The parallel flow configuration was identified as the most efficient for heat removal. A detailed analysis of the numerical results highlights the underlying physics behind the thermal performance differences among the flow configurations.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2025

Planar Collisionless Shock Simulations with the Semi-implicit Particle-in-cell Model FLEKS

This study investigates the applicability of the semi-implicit particle-in-cell code FLexible Exascale Kinetic Simulator (FLEKS) to heliospheric shock simulations. We examine one- and two-dimensional local planar shock simulations, initialized using MHD states with upstream conditions representative of plasmas in the hypersonic, β ∼ 1 regime, for both quasi-perpendicular and quasi-parallel configurations. The refined algorithm in FLEKS proves robust, enabling accurate shock simulations with a grid resolution on the order of the electron inertial length d e . Our simulations successfully capture key shock features, including shock structures (foot, ramp, overshoot, and undershoot), upstream and downstream waves (fast magnetosonic, whistler, Alfvén ion-cyclotron, and mirror modes), and non-Maxwellian particle distributions. Crucially, we find that at least two spatial dimensions are critical for accurately reproducing downstream-wave physics in quasi-perpendicular shocks and capturing the complex dynamics of quasi-parallel shocks, including surface rippling, shocklets, short, large-amplitude magnetic structures, magnetic reconnection, and jets. Furthermore, our parameter studies demonstrate the impact of mass ratio and grid resolution on shock physics. This work provides valuable guidance for selecting appropriate physical and numerical parameters for shock simulations using a semi-implicit PIC method, paving the way for incorporating kinetic shock processes into large-scale collisionless plasma simulations with the MHD-AEPIC model.

plasma astrophysics↗
DOE OSTI2020

Causal Interactions between Southern Ocean Polynyas and High-Latitude Atmosphere–Ocean Variability

Abstract Weddell Sea open-ocean polynyas have been observed to occasionally release heat from the deep ocean to the atmosphere, indicating that their sporadic appearances may be an important feature of high-latitude atmosphere–ocean variability. Yet, observations of the phenomenon are sparse and many standard-resolution models represent these features poorly, if at all. We use a fully coupled, synoptic-scale preindustrial control simulation of the Energy Exascale Earth System Model (E3SMv0-HR) to effectively simulate open-ocean polynyas and investigate their role in the climate system. Our approach employs statistical tests of Granger causality to diagnose local and remote drivers of, and responses to, polynya heat loss on interannual to decadal time scales. First, we find that polynya heat loss Granger causes a persistent increase in surface air temperature over the Weddell Sea, strengthening the local cyclonic wind circulation. Along with responding to polynyas, atmospheric conditions also facilitate their development. When the Southern Ocean experiences a rapid poleward shift in the circumpolar westerlies following a prolonged negative phase of the southern annular mode (SAM), Weddell Sea salinity increases, promoting density destratification and convection in the water column. Finally, we find that the reduction of surface heat fluxes during periods of full ice cover is not fully compensated by ocean heat transport into the high latitudes. This imbalance leads to a buildup of ocean heat content that supplies polynya heat loss. These results disentangle the complex, coupled climate processes that both enable the polynya’s existence and respond to it, providing insights to improve the representation of these highly episodic sea ice features in climate models.

Kaufman, Zachary S.↗
DOE OSTI2023

Insights into Waterflooding in Hydrocarbon-Bearing Nanochannels of Varying Cross Sections from Mesoscopic Multiphase Flow Simulations

Waterflooding is one of the geotechnique used to recover fuel sources from nanoporous geological formations. The scientific understanding of the process that involves the multiphase flow of nanoconfined fluids, however, has lagged, mainly due to the complex nanopore geometries and chemical compositions. To enable benchmarked flow of nanoconfined fluids, the architected geomaterials, such as synthesized mesoporous silica with tunable pore shapes and surface chemical properties, are used for designing and conducting experiments and simulations. This work uses a modified many-body dissipative particle dynamics (mDPD) model with accurately calibrated parameters to perform parametric flow simulations for studying the influences of waterflooding driven power, pore shape, surface roughness, and surface wettability on the multiphase flow in heptane-saturated silica nanochannels. Remarkably, up to 80\% reduction in the effective permeability is found for water-driven heptane flow in a baseline 4.5 nm-wide slit channel, when compared with the Hagen–Poiseuille equation. In the 4.5 nm-wide channels with architected surface roughness, the flow rate is found either higher or lower than the baseline case, depending on the shape and size of cross-sections. High wettability of the solid surface to water is essential for achieving high recovery of heptane, regardless of surface roughness. When the solid surface is less wetting or non-wetting to water, the existence of an optimal waterflooding driven power is found to allow for the highest possible recovery. A detailed analysis on the evolution of the transient water-heptane interface in those nanochannels is presented to elucidate the underlying mechanisms that impact or dictate the multiphase flow behaviors.

58 GEOSCIENCES↗
DOE OSTI2023

Mapping Glacier Basal Sliding Applying Machine Learning

During the RESOLVE project (“High-resolution imaging in subsurface geophysics: development of a multi-instrument platform for interdisciplinary research”), continuous surface displacement and seismic array observations were obtained on Glacier d’Argentière in the French Alps for 35 days in May 2018. The data set is used to perform a detailed study of targeted processes within the highly dynamic cryospheric environment. In particular, the physical processes controlling glacial basal motion are poorly understood and remain challenging to observe directly. Especially in the Alpine region for temperate based glaciers where the ice rapidly responds to changing climatic conditions and thus, processes are strongly intermittent in time and heterogeneous in space. Spatially dense seismic and Global Positioning System (GPS) measurements are analyzed applying machine learning to gain insight into the processes controlling glacial motions of Glacier d’Argentière. Using multiple bandpass-filtered copies of the continuous seismic waveforms, we compute energy-based features, develop a matched field beamforming catalog and include meteorological observations. Features describing the data are analyzed with a gradient boosting decision tree model to directly estimate the GPS displacements from the seismic noise. We posit that features of the seismic noise provide direct access to the dominant parameters that drive displacement on the highly variable and unsteady surface of the glacier. The machine learning model infers daily fluctuations and longer term trends. The results show on-ice displacement rates are strongly modulated by activity at the base of the glacier. The techniques presented provide a new approach to study glacial basal sliding and discover its full complexity.

58 GEOSCIENCES↗
DOE OSTI2025

Surrogate models to optimize plasma-assisted atomic layer deposition in high aspect ratio features

In this work, we explore surrogate models to optimize plasma enhanced atomic layer deposition (PEALD) in high aspect ratio features. In plasma-based processes such as PEALD and atomic layer etching (ALE), surface recombination can dominate the reactivity of plasma species with the surface, which can lead to unfeasibly long exposure times to achieve full conformality inside nanostructures like high aspect ratio vias. Using a synthetic dataset based on simulations of PEALD, we train artificial neural networks to predict saturation times based on cross section thickness data obtained for partially coated conditions. The results obtained show that just two experiments in undersaturated conditions contain enough information to predict saturation times within 10% of the ground truth. A surrogate model trained to determine whether surface recombination dominates the plasma–surface interactions in a PEALD process achieves 99% accuracy. This demonstrates that machine learning can provide a new pathway to accelerate the optimization of PEALD processes in areas such as microelectronics. Our approach can be easily extended to ALE and more complex structures.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Leveraging Additive Manufacturing to Fabricate High Temperature Alloys With Co-Designed Mechanical Properties and Environmental Resistance

A paradigm shift in the traditional sequential design approaches is critically essential to create application-specific hierarchical and multifunctional materials with superior long-term performance for next-generation energy technologies involving extreme environments. In the current work, we aim to leverage the flexibility and geometric/compositional complexity offered by additive manufacturing to demonstrate this new approach by codesigning a compositionally graded Ni-based alloy for molten salts\sCO 2 heat exchangers to enable mitigation of environmental degradation of surfaces exposed to molten halide salts, while simultaneously suppressing the consequent deterioration in mechanical stability. Further, thermokinetic modeling describing the underlying physics of thermally- and environmentally induced spatiotemporal compositional and microstructural evolution will be employed to predict the parameter space of material deposition processes and precisely identify the required composition gradient. Preliminary corrosion and mechanical testing of the dual material demonstrated the potential of the material to replace existing solid solution strengthened materials for this application.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Implications of reduced-complexity aerosol thermodynamics on organic aerosol mass concentration and composition over North America

Atmospheric organic aerosol (OA) mass concentrations can be affected by water uptake through its impact on the gas–particle partitioning of semi-volatile compounds. Current chemical transport models (CTMs) neglect this process. We have implemented the Binary Activity Thermodynamics model coupled to a volatility basis set partitioning scheme in the GEOS-Chem CTM, providing an efficient reduced-complexity OA model that predicts relative-humidity-dependent mixing and partitioning thermodynamics, while limiting the impact on computational efficiency. We provide a quantitative assessment of this water-sensitive OA treatment, focusing on a subdomain over North America. The updated OA scheme predicts a spatiotemporal mean enhancement in surface-level OA mass concentration of 145 % for January 2019 and 76 % for July 2019 compared to GEOS-Chem's most advanced OA scheme. The temporal mean surface-level OA organic mass concentration can increase by up to ∼590 % for January 2019 and ∼280 % for July 2019, with the greatest enhancements occurring over the ocean. The updated OA scheme also quantifies the OA-associated water content. The simulations show how different OA precursors and related OA surrogates contribute and respond to water uptake, including that due to changes in temperature and relative humidity over the diurnal cycle in selected winter and summer months. These results are independent of future CTM improvements involving updates to chemical reaction schemes and emission inventories. Our water-sensitive OA scheme allows for a better representation of the seasonal and regional variations in OA mass concentration in CTMs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

“Lagrangian disks” in M-theory

While the study of bordered (pseudo-)holomorphic curves with boundary on Lagrangian submanifolds has a long history, a similar problem that involves (special) Lagrangian submanifolds with boundary on complex surfaces appears to be largely overlooked in both physics and math literature. We relate this problem to geometry of coassociative submanifolds in G 2 holonomy spaces and to Spin(7) metrics on 8-manifolds with T 2 fibrations. As an application to physics, we propose a large class of brane models in type IIA string theory that generalize brane brick models on the one hand and 2d theories T[M 4 ] on the other.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Influence of ethylenediaminetetraacetic acid on the long-term oxidation state distribution of plutonium

Here, spectrophotometry was used to study the effect of EDTA on plutonium oxidation state distribution as a function of time, pH, and ligand-to-metal ratio (L/M) under anoxic conditions. Novel Pu(V)-EDTA absorption bands were identified at 571, 993, 1105, and 1150 nm with molar absorption coefficients of 15 ± 1, 6 ± 1, 10 ± 1, and 10 ± 1 cm –1 M –1 , respectively. Pu(V)-EDTA spectral changes occurred at L/M < 1, indicating only Pu V O 2 (EDTA) 3- formed with logK = 3.6 ± 0.3. Time-resolved experiments showed EDTA drastically increased the Pu(V/VI) reduction rate, which we propose is driven by amine lone-pair electron donation and the oxidative decarboxylation of EDTA. Oxidation of Pu(III)-EDTA to Pu(IV)-EDTA occurred on a slower time scale (110–237 days) than previously reported (<15 min) and is hypothesized to be radiolysis driven. Pu(V/VI)-EDTA and Pu(III)-EDTA both approached Pu(IV)-EDTA stabilization over time, yet Pu(V/VI)-EDTA solubility data was ≥ 1.0 log 10 units higher than predicted by Pu(IV)-EDTA solubility models, indicating that current thermodynamic models are incomplete. Ultimately, the data show EDTA preferentially stabilizes Pu(IV) over time regardless of initial oxidation state, but Pu(V)-EDTA can persist under environmentally-relevant conditions, emphasizing the need to continue investigating redox reactions, speciation, and behavior of these complexes to support the transuranic waste disposal and surface remediation/containment efforts.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2022

Utilization of coupled eigenmodes in Akiyama atomic force microscopy probes for bimodal multifrequency sensing

Akiyama atomic force microscopy probes represent a unique means of combining several of the desirable properties of tuning fork and cantilever probe designs. As a hybridized mechanical resonator, the vibrational characteristics of Akiyama probes result from a complex coupling between the intrinsic vibrational eigenmodes of its constituent tuning fork and bridging cantilever components. Here, through a combination of finite element analysis modeling and experimental measurements of the thermal vibrations of Akiyama probes, we identify a complex series of vibrational eigenmodes and measure their frequencies, quality factors, and spring constants. We then demonstrate the viability of Akiyama probes to perform bimodal multi-frequency force sensing by performing a multimodal measurement of a surface's nanoscale photothermal response using photo-induced force microscopy imaging techniques. Further performing a parametric search over alternative Akiyama probe geometries, we propose two modified probe designs to enhance the capability of Akiyama probes to perform sensitive bimodal multifrequency force sensing measurements.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2025

Comparing multi-source urban flood indicators: satellite, simulation, and citizen-reported data

Urban flooding arises from complex mechanisms, making it challenging to capture accurately with a single detection method. This study evaluates three complementary approaches to detect flooding across three Chicago neighborhoods: (i) Sentinel-1 synthetic aperture radar (SAR), offering weather-independent, high-resolution (10 m) imagery of surface inundation; (ii) the storm water management model (SWMM), simulating combined sewer overflow and drainage performance; and (iii) citizen-generated 311 service requests, capturing observed flooding impacts. By analyzing six storms ranging from severe to mild, we examine how each source uniquely contributes to identifying urban flood events. SAR imagery effectively identifies standing water but can miss brief flooding due to satellite revisit constraints. SWMM provides detailed insights into system-wide drainage behavior yet may underestimate localized street-level flooding. Meanwhile, 311 calls reflect real-world flooding impacts but are vulnerable to underreporting. Statistical overlap analysis highlights chronic flood hotspots repeatedly identified across multiple detection methods, indicating persistent infrastructure and topographic vulnerabilities. Temporal analysis further reveals that while SWMM flooding aligns closely with rainfall peaks, 311 calls typically precede or persist beyond these peaks. Our findings emphasize the value of using satellite observations, hydrological modeling, and resident-reported data in a complementary manner to better interpret patterns in flood timing, severity, and spatial distribution—providing insights that can inform targeted infrastructure improvements and contribute to urban flood resilience planning.

311↗
DOE OSTI2021

Three-Dimensional Site Effects at Mortandad Canyon, NM: I. Validation against ambient noise recordings of topography effects

Observations from large earthquakes have shown that the presence of a strong topographic relief can significantly aggravate the catastrophic consequences of strong seismic motions. Structures on the tops of hills, ridges, and canyons, for example, have suffered systematically greater damage than similar structures at the hill bases or on level ground. Various earthquake recordings have confirmed the macroseismic observations, indicating systematic amplification of seismic motion over convex topographies, deamplification over concave topographic features, and complex amplification and deamplification patterns on hill slopes that result in significant differential motion. Despite the documented evidence of their role in elevating seismic risk, however, topography effects are not accounted for in the majority of seismic code provisions and microzonation studies, which are based instead on models of flat earth surface.

58 GEOSCIENCES↗
DOE OSTI2023

Arctic Aerosol Sources and Mixing States Field Campaign Report

The Arctic is warming at a faster rate than anywhere else on Earth, with rapidly shrinking sea ice extent transforming the region. Depending on chemical and physical properties, atmospheric aerosols directly scatter and/or absorb radiation, serve as cloud droplet and/or ice crystal nuclei, and/or reduce the reflectiveness of the snow surface, thereby altering the atmospheric energy budget. There is a wide spread in the magnitude of simulated arctic aerosol radiative forcing, and significant differences in aerosol concentration levels and seasonal cycles often exist between models and observations. Increasing local natural and anthropogenic emissions are significant, with uncertain climate impacts due to complex feedbacks. Model evaluations, however, are limited by the dearth of arctic aerosol observations available and an inadequate understanding of arctic aerosol processes. The majority of previous arctic aerosol observations have been made through intensive spring/summer field campaigns, with few intensive measurement studies focused on the fall-winter transition, a period when freeze-up is occurring later and thinning sea ice is resulting in wintertime ice fracturing. Aerosol monitoring at arctic coastal stations has provided knowledge of long-term seasonal trends in aerosol concentrations. The completed ARM field campaign addresses observational and knowledge gaps through detailed aerosol size and chemical composition measurements during the fall-winter transition in the coastal Arctic and through the entire annual cycle in the central Arctic. The observations are improving our understanding of the sources and processes controlling the aerosol population in the rapidly changing Arctic.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Fundamental Insights into Cathode Stability: Linking Compositional Tuning and Local Coordination in Complex Metal Oxides under Aqueous Transformations

Compositional tuning of complex metal oxides in Li-ion battery materials influences their performance as well as their end-of-life behavior, in particular, the tendency to release toxic metal cations in aqueous solution. We modeled ternary variants of a parent LiCoO 2 delafossite structure by varying the metal identity and relative amounts. This yielded ten model formulations of Li(A 4/6 B 1/6 C 1/6 )O 2 , where the material is enriched with the A metal and doped with B and C, with Ni, Mn, Co, Fe, Al, V, and Ti as constituent metals. To assess their stability in aqueous conditions, metal release energetics were calculated using a combination of Density Functional Theory calculations and thermodynamics. Metal release in ternary oxides is dictated by subtle variations in the coordination environment of the leaving group. To identify governing chemical features across diverse compositions with varying local coordination environments, we leverage random forest regression and descriptor importance analysis. A key result is that metal–oxygen orbital hybridization, quantified using a projected density-of-states-derived descriptor, H d/p , provides a physically grounded measure of interaction strength that governs metal release energetics. This refined perspective goes beyond conventional oxidation state considerations and offers more robust insights for materials science. Finally, we model defect surface-bound O 2 dimer formation as a proxy for reactive oxygen species (ROS) generation. The results show that Ni-rich compositions more readily stabilize spin-polarized O 2 dimers, corroborating experimental reports of an increased ROS-driven biological response. In conclusion, our results establish a compositional and electronic basis for metal release and surface oxygen reactivity that form a rationale for complex metal oxide design principles.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Evaluation of global teleconnections in CMIP6 climate projections using complex networks

In climatological research, the evaluation of climate models is one of the central research subjects. As an expression of large-scale dynamical processes, global teleconnections play a major role in interannual to decadal climate variability. Their realistic representation is an indispensable requirement for the simulation of climate change, both natural and anthropogenic. Therefore, the evaluation of global teleconnections is of utmost importance when assessing the physical plausibility of climate projections. We present an application of the graph-theoretical analysis tool δ-MAPS, which constructs complex networks on the basis of spatio-temporal gridded data sets, here sea surface temperature and geopotential height at 500 hPa. Complex networks complement more traditional methods in the analysis of climate variability, like the classification of circulation regimes or empirical orthogonal functions, assuming a new non-linear perspective. While doing so, a number of technical tools and metrics, borrowed from different fields of data science, are implemented into the δ-MAPS framework in order to overcome specific challenges posed by our target problem. Those are trend empirical orthogonal functions (EOFs), distance correlation and distance multicorrelation, and the structural similarity index. δ-MAPS is a two-stage algorithm. In the first place, it assembles grid cells with highly coherent temporal evolution into so-called domains. In a second step, the teleconnections between the domains are inferred by means of the non-linear distance correlation. We construct 2 unipartite and 1 bipartite network for 22 historical CMIP6 climate projections and 2 century-long coupled reanalyses (CERA-20C and 20CRv3). Potential non-stationarity is taken into account by the use of moving time windows. The networks derived from projection data are compared to those from reanalyses. Our results indicate that no single climate projection outperforms all others in every aspect of the evaluation. But there are indeed models which tend to perform better/worse in many aspects. Differences in model performance are generally low within the geopotential height unipartite networks but higher in sea surface temperature and most pronounced in the bipartite network representing the interaction between ocean and atmosphere.

58 GEOSCIENCES↗
DOE OSTI2020

Pb Sorption at the Barite (001)–Water Interface

Sorption of ions at the mineral-water interface is an important factor that determines the fate of toxic metals in the environment. Here, we use barite as a model substrate to understand the interaction of toxic-metal lead (Pb) with ionic crystals. The coverage and location of Pb sorbed to the (001) surface was measured as a function of aqueous Pb concentration ([Pb] aq ) using in situ specular resonant anomalous X-ray reflectivity (RAXR) to determine the sorption capacity and process. In this work, the results show that Pb sorption occurs via incorporation (primarily within the top barite layer similar to 3 angstrom in depth) and adsorption (mostly as an inner-sphere complex at similar to 2 angstrom in height) simultaneously. Both the incorporated and adsorbed Pb coverages increase with increasing [Pb] aq up to [Pb] aq similar to 200 mu M, above which the adsorbed fraction increases more rapidly than the incorporated fraction. This enhanced adsorption has a height distribution that is further extended (>= 15 angstrom from the surface) than that observed in lower [Pb] aq . This change in distribution is interpreted as arising from additional sorption of outer-sphere species or Pb-bearing phases precipitated on the surface. Desorption experiments in Pb-free solutions show that the incorporated fraction is more resistant to removal than the adsorbed fraction, consistent with the speciation-dependent stabilities premised in the classical sorption models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗