Search NASA⌕ Search

SEARCH · Search NASA

Results for “interface & surface thermodynamics”

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 91 records · Page 5

Examination of actinide chemistry at solid-water interfaces to support advanced actinide separations (Final Report)

In this work we have examined the thermodynamics of actinide interactions with solid interfaces in the presence and absence of multi-functional complexing ligands. The term “interactions” is a global term to capture many potential reactions including aqueous complexation, sorption/surface complexation, surface-mediated redox reactions, electrostatic attraction and repulsion, surface passivation, surface precipitation, and counter-ion competition. Each of these potential reactions may control actinide interactions with a solid surface. This work is helping to develop a more fundamental understanding of actinide bonding at solid:water interfaces and examine unique surface mediated oxidation/reduction reactions which are not expected based on the reduction potentials of aqueous species. We place a large emphasis on understanding both sorption and desorption processes as a means of understanding the sorption mechanisms. This work is directly aligned with the DOE BES Heavy Element Chemistry programs interests in understanding “the chemical and physical properties of these elements to determine solution, interfacial and solid-state bonding and reactivity”.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

XPS study of interface formation of CVD SiO2 on InSb

The interfacial chemistry of CVD SiO2 films deposited on thin native oxides grown on InSb substrates is examined using X-ray photoemission spectroscopy (XPS) and a relatively benign chemical etching technique for depth profiling. An intensity analysis of XPS spectra is used to derive the compositional structure of the interfaces obtained in the SiO2/native oxide/InSb system. Peak positions in these spectra are used to follow the change in substrate surface potential during the etch sequence, and to establish the chemical nature of the species formed during deposition and subsequent processing. Reaction of the substrate with oxygen resulted in an In-rich native oxide and 1-2 monolayers of excess elemental Sb at the native-oxide/substrate interface, incompletely oxidized silane reduced the native oxide, leaving less than 1 monolayer of elemental In at the SiO2/native oxide interface. Etch removal of this thin In-rich layer leads to a change in the substrate surface potential of 0.06 eV, corresponding to a net increase in positive charge. The results are consistent with simple thermodynamic considerations; they are also compared to previously reported studies of deposited dielectrics on III-V compound semiconductors.

Vasquez, R. P.↗

Dynamics of Liquids in Edges and Corners (DYLCO): IML-2 Experiment for the BDPU

Knowledge of the behavior of fluids possessing free surfaces is important to many fluid systems, particularly in space, where the normally subtle effects of surface wettability play a more dramatic and often surprising role. DYLCO for the IML-2 mission was proposed as a simple experiment to probe the particular behavior of capillary surfaces in containers of irregular cross section. Temperature control was utilized to vary the fluid-solid contact angle, a questionable thermodynamic parameter of the system, small changes in which can dramatically influence the configuration, stability, and flow of a capillary surface. Container shapes, test fluid, and temperature ranges were selected for observing both local changes in interface curvature as well as a global change in fluid orientation due to a critical wetting phenomenon. The experiment hardware performed beyond what was expected and fluid interfaces could be readily digitized post flight to show the dependence of the interface curvature on temperature. For each of the containers tested surfaces were observed which did not satisfy the classic equations for the prediction of interface shape with constant contact angle boundary condition. This is explained by the presence of contact angle hysteresis arising from expansion and contraction of the liquid during the heating and cooling steps of the test procedure. More importantly, surfaces exceeding the critical surface curvature required for critical wetting were measured, yet no wetting was observed. These findings are indeed curious and pose key questions concerning the role of hysteresis for this critical wetting phenomena. The stability of such surfaces was determined numerically and it is shown that stability is enhance (reduced) when a surface is in its 'advancing' ('receding') state, The analysis shows complete instability as the critical wetting condition is reached. The case of ideal dynamic wetting is addressed analytically in detail with results of significant flow characteristics presented in closed form. The solutions indicate a square root of T dependence of the capillary 'rise' rate which is corroborated by drop tower tests. The analysis clearly shows that infinite time is necessary for surfaces to reorient at the critical wetting transition.

Langbein, D.↗

Ensembles of Metastable States Govern Heterogeneous Catalysis on Dynamic Interfaces

Heterogeneous catalysis is at the heart of chemical industry. Being able to tune and design efficient catalysts for processes of interest is of an utmost importance, and for this, the molecular-level understanding of heterogeneous catalysts is the first step, and indeed a prime focus of the modern catalysis research. For a long time, a single most thermodynamically stable structure of the catalytic interface attained in reaction conditions had been envisioned as the reactive phase. However, some catalytic interfaces continue to undergo structural dynamics in the steady state, triggered by high temperatures, pressures, and binding and changing reagents. Among particularly dynamic interfaces are such widely-used catalysts as crystalline and amorphous surfaced supporting (sub-)nano metallic clusters. Recently, it became clear that this dynamic fluxionality causes the supported clusters to populate many distinct structural and stoichiometric states in catalytic conditions. Hence, the catalytic interface should be viewed as an evolving statistical ensemble of many (not one) structures. As such, every member in the ensemble contributes to the properties of the catalyst differently, and in proportion to its probability to be populated. This new notion flips the established paradigm and calls for new theory, modeling approaches, operando measurements, and updated design strategies. The statistical ensemble nature of surface-supported sub-nano cluster catalysts can be exemplified by oxide-supported and adsorbate-covered Pt, Pd, Cu, CuPd clusters, catalytic toward oxidative and non-oxidative dehydrogenation. They have access to a variety of 3D and quasi-2D shapes. The compositions of their thermal ensembles are dependent on the cluster size, leading to size-specific catalytic activities and the famous “every atom counts” phenomenon. The support and adsorbates affect catalyst structures, and state of the reacting species causes the ensemble to change in every reaction intermediate. The most stable member of the ensemble dominates the thermodynamic properties of the corresponding intermediate, whereas the kinetics can be determined by more active but less populated metastable catalyst states, and that suggests that many earlier studies might have overlooked the actual active sites. Both effects depend on the relative time-scales of catalyst restructuring and reaction dynamics. The catalyst may routinely operate off-equilibrium. Ensemble phenomena lead to surprising exceptions from established rules of catalysis, such as scaling relations, and Arrhenius behavior. Catalyst deactivation is also an ensemble property, and its extent of mitigation can be predicted through the new paradigm. These findings were enabled by advances in theory, such as global optimization and subsequent utilization of multiple local minima, and pathways sampling, as well as operando catalyst characterization. The fact that the per-site and per-species resolution is needed for the description and predicting of catalyst properties gives theory the central role in catalysis research, as most experiments provide ensemble-average information and cannot detect the crucial minority species that may be responsible for the catalytic activity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Tailoring Interfaces for Enhanced Methanol Production from Photoelectrochemical CO 2 Reduction

Efficient and stable photoelectrochemical reduction of CO 2 into highly reduced liquid fuels remains a formidable challenge, which requires an innovative semiconductor/catalyst interface to tackle. In this study, we introduce a strategy involving the fabrication of a silicon micropillar array structure coated with a superhydrophobic fluorinated carbon layer for the photoelectrochemical conversion of CO 2 into methanol. The pillars increase the electrode surface area, improve catalyst loading and adhesion without compromising light absorption, and help confine gaseous intermediates near the catalyst surface. The superhydrophobic coating passivates parasitic side reactions and further enhances local accumulation of reaction intermediates. Upon one-electron reduction of the molecular catalyst, the semiconductor–catalyst interface changes from adaptive to buried junctions, providing a sufficient thermodynamic driving force for CO 2 reduction. These structures together create a unique microenvironment for effective reduction of CO 2 to methanol, leading to a remarkable Faradaic efficiency reaching 20% together with a partial current density of 3.4 mA cm –2 , surpassing the previous record based on planar silicon photoelectrodes by a notable factor of 17. Furthermore, this work demonstrates a new pathway for enhancing photoelectrocatalytic CO 2 reduction through meticulous interface and microenvironment tailoring and sets a benchmark for both Faradaic efficiency and current density in solar liquid fuel production.

14 SOLAR ENERGY↗

Parallel Grand Canonical Monte Carlo (ParaGrandMC) Simulation Code

This report provides an overview of the Parallel Grand Canonical Monte Carlo (ParaGrandMC) simulation code. This is a highly scalable parallel FORTRAN code for simulating the thermodynamic evolution of metal alloy systems at the atomic level, and predicting the thermodynamic state, phase diagram, chemical composition and mechanical properties. The code is designed to simulate multi-component alloy systems, predict solid-state phase transformations such as austenite-martensite transformations, precipitate formation, recrystallization, capillary effects at interfaces, surface absorption, etc., which can aid the design of novel metallic alloys. While the software is mainly tailored for modeling metal alloys, it can also be used for other types of solid-state systems, and to some degree for liquid or gaseous systems, including multiphase systems forming solid-liquid-gas interfaces.

Vesselin I Yamakov↗

A computational pipeline to generate a synthetic dataset of metal ion sorption to oxides for AI/ML exploration

The charged mineral/electrolyte interfaces are ubiquitous in the surface and subsurface–including the surroundings of the geological disposal sites for radioactive waste. Therefore, understanding how ions interact with charged surfaces is critically important for predicting radionuclide mobility in the case of waste leakage. At present, the Surface Complexation Models (SCMs) are the most successful thermodynamic frameworks to describe ion retention by mineral surfaces. SCMs are interfacial speciation models that account for the effect of the electric field generated by charged surfaces on sorption equilibria. These models have been successfully used to analyze and interpret a broad range of experimental observations including potentiometric and electrokinetic titrations or spectroscopy. Unfortunately, many of the current procedures to solve and fit SCM to experimental data are not optimal, which leads to a non-transferable or non-unique description of interfacial electrostatics and consequently of the strength and extent of ion retention by mineral surfaces. Recent developments in Artificial Intelligence (AI) offer a new avenue to replace SCM solvers and fitting algorithms with trained AI surrogates. Unfortunately, there is a lack of a standardized dataset covering a wide range of SCM parameter values available for AI exploration and training–a gap filled by this study. Here, we described the computational pipeline to generate synthetic SCM data and discussed approaches to transform this dataset into AI-learnable input. First, we used this pipeline to generate a synthetic dataset of electrostatic properties for a broad range of the prototypical oxide/electrolyte interfaces. The next step is to extend this dataset to include complex radionuclide sorption and complexation, and finally, to provide trained AI architectures able to infer SCMs parameter values rapidly from experimental data. Here, we illustrated the AI-surrogate development using the ensemble learning algorithms, such as Random Forest and Gradient Boosting. These surrogate models allow a rapid prediction of the SCM model parameters, do not rely on an initial guess, and guarantee convergence in all cases.

Li, Chunhui↗

Cryogenic Fluid Technologies for Long Duration In-Space Operations

Reliable knowledge of low-gravity cryogenic fluid management behavior is lacking and yet is critical in the areas of storage, distribution, and low-gravity propellant management. The Vision for Space Exploration mission objectives will require the use of high performance cryogenic propellants (hydrogen, oxygen, and methane). Additionally, lunar missions will require success in storing and transferring liquid and gas commodities on the surface. The fundamental challenges associated with the in-space use of cryogens are their susceptibility to environmental heat, their complex thermodynamic and fluid dynamic behavior in low gravity and the uncertainty of the position of the liquid-vapor interface if the propellants are not settled. The Cryogenic Fluid Management (CFM) project is addressing these issues through ground testing and analytical model development, and has crosscutting applications and benefits to virtually all missions requiring in-space operations with cryogens. Such knowledge can significantly reduce or even eliminate tank fluid boil-off losses for long term missions, reduce propellant launch mass and on-orbit margins, and simplify vehicle operations. The Cryogenic Fluid Management (CFM) Project is conducting testing and performing analytical evaluation of several areas to enable NASA s Exploration Vision. This paper discusses the content and progress of the technology focus areas within CFM.

Motil, Susan M.↗

Parallel Grand-Canonical Monte Carlo (ParaGrandMC) User’s Manual Version 2.0

This manual describes the commands and command line options for the Parallel Grand Canonical Monte Carlo version 2.0 (ParaGrandMC.2.0) simulation code. This is a highly scalable parallel FORTRAN 2003 code for simulating the thermodynamic evolution of materials at the atomic level, and predicting their thermodynamic state, phase diagram, chemical composition and mechanical properties. The code is specifically designed to simulate multi-component alloy systems, predict solid-state phase transformations such as austenite-martensite transformations, precipitate formation, recrystallization, capillary effects at interfaces, surface absorption, etc., which can aid the design of novel metallic alloys. While the software is mainly tailored for modeling metal alloys, it can also be used for other types of solid-state systems, and to some degree for liquid or gaseous systems, including multiphase systems forming solid-liquid-gas interfaces. In addition to performing Monte Carlo (MC) simulations, the code can also perform Molecular Dynamics (MD) and Langevin Dynamics (LD) simulations, which can be combined and interchanged with MC for faster and more efficient system evolution. A detailed description of the MC part of the code is provided in the NASA ParaGrandMC report: NASA/CR–2016-219202; http://www.sti.nasa.gov.

High performance computing↗

DFT modelling of explicit solid–solid interfaces in batteries: methods and challenges

Density Functional Theory (DFT) calculations of electrode material properties in high energy density storage devices like lithium batteries have been standard practice for decades. By contrast, DFT modelling of explicit interfaces in batteries arguably lacks universally adopted methodology and needs further conceptual development. In this paper, we focus on solid–solid interfaces, which are ubiquitous not just in all-solid state batteries; liquid-electrolyte-based batteries often rely on thin, solid passivating films on electrode surfaces to function. We use metal anode calculations to illustrate that explicit interface models are critical for elucidating contact potentials, electric fields at interfaces, and kinetic stability with respect to parasitic reactions. The examples emphasize three key challenges: (1) the “dirty” nature of most battery electrode surfaces; (2) voltage calibration and control; and (3) the fact that interfacial structures are governed by kinetics, not thermodynamics. To meet these challenges, developing new computational techniques and importing insights from other electrochemical disciplines will be beneficial.

25 ENERGY STORAGE↗

The Role of Cooperative Interactions Among Surfaces, Solvents, and Reactive Intermediates on Catalysis at Liquid–Solid Interfaces (Final Report DE-SC0020224)

This project established quantitative links between inner‑sphere chemistry (active metal identity, coordination, and zeolite topology) and outer‑sphere organization (solvent identity, hydrogen‑bond networks, and pore condensation) that govern rates, activation barriers, and selectivities for alkene epoxidation and epoxide ring‑opening at solid–liquid and quasi‑liquid–solid interfaces. We deconvoluted contributions from covalent interactions at active sites and noncovalent, solvent‑mediated interactions within pores by pairing well‑defined metal substituted zeolites with controlled solvent environments. We then mapped those contributions onto measurable kinetics (ΔH‡, ΔS‡), adsorption thermodynamics (ITC), and in situ spectroscopy. The transferrable outcomes include a set of design rules that include the following understandings. First, tune silanol ((SiOH)x) density and pore topology to organize solvent networks that selectively stabilize transition states. Second, exploit activity‑coefficient‑normalized rates and adsorption– barrier correlations to diagnose when solvent reorganization rather than surface chemistry limits performance. Third, use partial pore condensation (e.g., acetonitrile, water but also generalizable to other solvents) to elicit liquid‑like stabilization effects even in nominally vapor‑phase reactors. Collectively, these results provide strategies to increase epoxidation rates, improve oxidant utilization (i.e., selectivities), and steer regioselectivity in zeolite‑based catalytic processes relevant to sustainable oxidation chemistry. These outcomes should be transferable to other classes of reactions that proceed in microporous materials and under confinement provided by organized solvents (e.g., electrochemical double layers).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Atomistic thermodynamics and kinetics of dicalcium silicate dissolution

Low-temperature cement manufacturing has garnered academic and industrial attention for its low environmental footprints. However, the sluggish hydration kinetics of the resultant cement affects their early-age strength development. This motivates fundamental studies to unravel the mechanistic picture of the dissolution process and discover science-informed pathways to accelerate hydration. Standard atomistic simulations seldomly exceed a microsecond making them impractical to study slow dissolution processes. We are using rare event sampling techniques, we provide the mechanistic picture of Ca 2+ ion dissolution from a kink site on the dicalcium silicate surface. The Ca 2+ ion dissolution is comprised of two sequential stages: breaking restraints from the kink sites to form a ledge adatom and detaching from the ledge/terrace adatom sites into the solution. The first and second stages feature free energy barriers of ~63 kJ/mol and ~ 29 kJ/mol respectively, making the first stage the rate-limiting step of the entire Ca 2+ dissolution kinetics. Using the reactive flux method, the rate and equilibrium constants for each reaction step are calculated, which yield the Ca 2+ ion activity of ~1.03 × 10 -5 . The diffusion calculations indicate that the surface effects lower the self-diffusion coefficient of Ca 2+ ions at the solid-water interface.

36 MATERIALS SCIENCE↗

Simulation insights into wetting properties of hydrogen-brine-clay for hydrogen geo-storage

Hydrogen geo-storage is attracting substantial interdisciplinary interest as a cost-effective and sustainable option for medium- and long-term storage. Hydrogen can be stored underground in diverse formations, including aquifers, salt caverns, and depleted oil and gas reservoirs. The wetting dynamics of the hydrogen-brine-rock system are critical for assessing both structural and residual storage capacities, and ensuring containment safety. Through molecular dynamics simulations, we explore how varying concentrations of cushion gases (CO 2 or CH 4 ) influence the wetting properties of hydrogen-brine-clay systems under geological conditions (15 MPa and 333 K). We employed models of talc and the hydroxylated basal face of kaolinite (kaoOH) as clay substrates. Our findings reveal that the effect of cushion gases on hydrogen-brine-clay wettability is strongly dependent on the clay-brine interactions. Notably, CO 2 and CH 4 reduce the water wettability of talc in hydrogen-brine-talc systems, while exerting no influence on the wettability of hydrogen-brine-kaoOH systems. Detailed analysis of free energy of cavity formation near clay surfaces, clay-brine interfacial tensions, and the Willard-Chandler surface for gas-brine interfaces elucidate the molecular mechanisms underlying wettability changes. Our simulations identify empirical correlations between wetting properties and the average free energy required to perturb a flat interface when clay-brine interactions are less dominant. Here, our thorough thermodynamic analysis of rock-fluid and fluid-fluid interactions, aligning with key experimental observations, underscores the utility of simulated interfacial properties in refining contact angle measurements and predicting experimentally relevant properties. These insights significantly enhance the assessment of gas geo-storage potential. Prospectively, the approaches and findings obtained from this study could form a basis for more advanced multiscale simulations that consider a range of geological and operational variables, potentially guiding the development and improvement of geo-storage systems in general, with a particular focus on hydrogen storage.

25 ENERGY STORAGE↗

A high-order, localized-artificial-diffusivity method for Eulerian simulation of multi-material elastic-plastic deformation with strain hardening

A high-order method for Eulerian simulation of material undergoing large elastic–plastic deformation is developed. Thermodynamically consistent hyperelastic constitutive relations are assumed, facilitating the treatment of solids, liquids, and gases in a unified manner. Here, the method enables the simulation of multi-material interactions using a diffuse interface approach. Numerical capturing of material interfaces, shock waves, contact surfaces, and elastic-plastic strain discontinuities using high-order compact-difference schemes is assisted by Localized Artificial Diffusivity (LAD). In the new setting involving elastic–plastic deformation, the previously established terms for the artificial properties are verified to effectively regularize normal shocks. Additional LAD terms are introduced to the elastic and plastic kinematic equations to regularize shear shocks and other strain discontinuities, improving solution stability. Other important features of the method that improve robustness include the numerical treatment of compatibility terms in the kinematic equations, and the treatment of rotation. Particular emphasis is focused toward new advancements of the methods for plastic-deformation integration and the associated strain hardening of the material, including rate-dependent plasticity. The method is demonstrated on a variety of test problems, including 1-D impacts, a variant of the Shu-Osher problem, a Taylor impact, and a Richtmyer-Meshkov instability between two elastic–plastic solids with strain hardening.

42 ENGINEERING↗

Phase Separation, Capillarity, and Odd-Surface Flows in Chiral Active Matter

Active phase separations evade canonical thermodynamic descriptions and have thus challenged our understanding of coexistence and interfacial phenomena. Considerable progress has been made towards a nonequilibrium theoretical description of these traditionally thermodynamic concepts. Spatial parity symmetry is conspicuously assumed in much of this progress, despite the ubiquity of chirality in experimentally realized systems. Here, in this Letter, we derive a theory for the phase coexistence and interfacial fluctuations of a system that microscopically violates spatial parity. We find suppression of the phase separation as chirality is increased as well as the development of steady-state currents tangential to the interface dividing the phases. These odd flows are irrelevant to stationary interfacial properties, with stability, capillary fluctuations, and surface area minimization determined entirely by the capillary surface tension. Using large-scale Brownian dynamics simulations, we find excellent agreement with our theoretical scaling predictions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Kinetically versus thermodynamically controlled factors governing elementary pathways of GaP(111) surface oxidation

GaP and related III-V semiconductors have attracted interest as high-efficiency photoelectrochemical electrodes for solar water splitting. However, their efficacy is linked to the presence, identity, and integrity of native surface oxides, which are structurally and chemically complex and evolve during operation. Using ambient pressure X-ray photoelectron spectroscopy (APXPS) coupled with ab initio simulations, we track key chemical motifs expressed during evolution of GaP(111) surface oxides, their associated reaction kinetics, and their correlations with electronic properties. We identify two distinct thermal regimes corresponding to kinetically and thermodynamically controlled oxidation. Below 600 K, exposure to O2 generates kinetically facile Ga-O-Ga configurations, whereas higher temperatures cause activated oxygen to insert into Ga–P bonds as part of a thermodynamically driven transformation into a complex, heterogeneous 3D network of surface POx (1≤x ≤ 4) groups and Ga2O3 species, the latter of which eventually dominates upon depletion of surface phosphorus. Here our study highlights the critical competition between kinetic and thermodynamic factors during GaP oxidation, yielding insights for fabricating stable III-P-based photoelectrodes with precisely engineered surface properties.

08 HYDROGEN↗

Near-field models and simulations of pellet ablation in tokamaks

Here, detailed numerical studies of the ablation of a single neon pellet in the plasma disruption mitigation parameter space have been performed. Simulations were carried out using FronTier, a hydrodynamic/low magnetic Reynolds number MHD code with explicit tracking of material interfaces. FronTier's physics models resolve the pellet surface ablation and the formation of a dense, cold cloud of ablated material, the deposition of energy from hot plasma electrons, the expansion of the ablation cloud along magnetic field lines, and the radiation losses. A local thermodynamic equilibrium model based on Saha equations has been used to resolve atomic processes in the cloud and Redlich–Kwong corrections to the ideal gas equation of state for cold and dense gases, which have been used near the pellet surface. The FronTier pellet code is a next generation of the code described. It has been validated against the semi-analytic improved neutral gas shielding model in the 1D spherically symmetric approximation. Main results include quantification of the influence of atomic processes and Redlich–Kwong corrections on the pellet ablation in spherically symmetric approximation and verification of analytic scaling laws in a broad range of pellet and plasma parameters. Using axially symmetric MHD simulations, the properties of ablation channels and the reduction of pellet ablation rates in magnetic fields of increasing strength have been studied. While the main emphasis has been given to neon pellets for the plasma disruption mitigation, selected results on deuterium fueling pellets have also been presented.

97 MATHEMATICS AND COMPUTING↗

Observations of the liquid/solid interface in low-gravity melting

Time-lapsed photography of the liquid/solid interface of a melting ice cylinder was taken on Skylab 3 over a period of three hours. The same experiment was simulated on earth such that morphological and thermodynamic differences could be noted. A study of the returned color film clearly shows the dominance of surface tension effects in low-gravity melting. In the Skylab experiment, the ends of the ice cylinder melted first with the water being driven by surface tension onto the cylindrical surfaces. At any time, the principle of minimum surface area governs the overall appearance of the water-ice globule which changed from a cylindrical to a spherical shape. The latent heat of melting in low-gravity is supplied only by radiation (81%) and conduction (19%); whereas in one-g, the convective (55%) and radiative (38%) mode of heat transfer dominates over the conductive portion (7%). Information is also provided on containerless melting and heat transfer in space in the absence of convective air currents.

Otto, G. H.↗