SEARCH · Engineering Papers

Results for “Interfacial force 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 91 records · Page 5

All sources NASA NTRS DOE OSTI

DOE OSTI2026

From bulk to surface: Structure and dynamics of amorphous alumina from deep potential molecular dynamics

Understanding the atomic-scale structure and dynamics of amorphous oxide surfaces is essential for interpreting their chemical reactivity, mechanical stability, and interfacial behavior, yet direct experimental characterization remains challenging. We employ Deep Potential (DP) molecular dynamics to generate large-scale, ab initio -quality models of amorphous Al 2 O 3 bulk glasses and melt-quenched free surfaces, enabling a quantitative analysis of both structure and relaxation dynamics with statistical confidence inaccessible to direct ab initio simulation. The trained DP model reproduces experimental liquid and glass structure, captures the cooling-rate dependence of the bulk glass transition, and corrects systematic biases in the polyhedral populations predicted by widely used classical force fields. At the free surface, mass density recovers to bulk values over ~10 Å, while local coordination requires a slightly wider subsurface region to fully converge. The outermost layer is oxygen-enriched, exhibits altered polyhedral connectivity with contracted Al–O bonds, and hosts a broad population of under-coordinated motifs (notably AlO 3 and OAl 2 ) whose abundances are governed by glass stability. These under-coordinated surface motifs exhibit distinct vibrational signatures and occur as locally paired Lewis acid and Brønsted base sites consistent with bond-valence compensation, yet remain spatially dispersed rather than aggregating into extended clusters. Despite this pronounced structural heterogeneity, surface relaxation and the glass-transition temperature remain comparable to their bulk counterparts, suggesting that the disordered surface is kinetically stable once formed. Together, these results establish a molecular-level picture of amorphous alumina surfaces and demonstrate the capability of machine-learned potentials to resolve structure–property relationships in disordered oxide interfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Mechanism-Resolved PFM of Ferroionic and Ferroelectric Responses in Thickness-Gradient Hf 0.5 Zr 0.5 O 2 Libraries

Resolving growth mechanisms and thickness evolution of functional properties is one of the key tasks in materials discovery and optimization involving thin-film materials, traditionally requiring significant experimental budgets. Here we introduce the combination of thickness-gradient libraries and automated scanning probe microscopy as a systematic pathway to elucidate growth modes and disentangle ferroelectric and electrochemical contributions in ferroelectric thin films. As a model system, we explore the Hf 0.5 Zr 0.5 O 2 (HZO) gradient thin films grown on La x Sr 1-x MnO 3 (LSMO) bottom electrode thin films. Automated piezoresponse force microscopy, spectroscopy, and lithography reveals that irreversible topographic deformation arises from electrochemical activity at the LSMO surface, whereas reversible phase inversion in HZO reflects ferroelectric switching. Automated topography height-map scans are used to further quantify nucleation density, particle-size evolution, and roughness correlations across the thickness-gradient, demonstrating that improved plume stabilization during growth suppresses interfacial reactions and promotes dense, fine-grained HZO conducive to ferroelectric phase formation. This combined materials-engineering and automated-SPM framework establishes a platform for high-throughput, mechanism-resolved characterization of ferroionic and ferroelectric responses in complex oxide films.

FOS: Physical sciences↗

NASA NTRS1995

Influence of Two-Phase Thermocapillary Flow on Liquid Retention in Microscopic Pores

An important feature of screened propellant acquisition devices is the retention capability or maximum maintainable pressure difference across the porous barrier separating the liquid and gas. Previous experiments with liquid hydrogen showed a marked reduction in retention when the tank containing the device was pressurized with hydrogen vapor. These tests, however did not indicate any appreciable degradation in retention with helium pressurization or direct heating through the screen. The objective of this article is to determine if the thermocapillary convection arising from phase change in the microscopic pores of such screens could cause these disparities in performance. A numerical model of flow in a single pore suggests that the thermocapillary-induced gradient in liquid pressure along the surface can strongly affect surface morphology. In an evaporative environment, this gradient exerts a stabilizing influence on surface curvature, and preserves the momentum balance between the liquid and gas. With condensation, it causes a force imbalance and a destabilizing suction in the middle of the pore that reduces retention. Results also indicate that introducing an inert gas, such as helium, suppresses this retention loss mechanism by lowering thermocapillary circulation and its associated interfacial pressure gradient.

Schmidt, George R.↗

NASA NTRS2013

Effect of Residual Noncondensables on Pressurization and Pressure Control of a Zero-Boil-Off Tank in Microgravity

The Zero-Boil-Off Tank (ZBOT) Experiment is a small-scale experiment that uses a transparent ventless Dewar and a transparent simulant phase-change fluid to study sealed tank pressurization and pressure control with applications to on-surface and in-orbit storage of propellant cryogens. The experiment will be carried out under microgravity conditions aboard the International Space Station in the 2014 timeframe. This paper presents preliminary results from ZBOT's ground-based research that focuses on the effects of residual noncondensable gases in the ullage on both pressurization and pressure reduction trends in the sealed Dewar. Tank pressurization is accomplished through heating of the test cell wall in the wetted and un-wetted regions simultaneously or separately. Pressure control is established through mixing and destratification of the bulk liquid using a temperature controlled forced jet flow with different degrees of liquid jet subcooling. A Two-Dimensional axisymmetric two-phase CFD model for tank pressurization and pressure control is also presented. Numerical prediction of the model are compared to experimental 1g results to both validate the model and also indicate the effect of the noncondensable gas on evolution of pressure and temperature distributions in the ullage during pressurization and pressure control. Microgravity simulations case studies are also performed using the validated model to underscore and delineate the profound effect of the noncondensables on condensation rates and interfacial temperature distributions with serious implications for tank pressure control in reduced gravity.

Marangoni and Natural Convection↗

NASA NTRS2001

The Transient Dendritic Solidification Experiment (TDSE)

Dendritic solidification is a common mode of solidification. It is also an important model problem in nonequilibrium physics and pattern formation physics. Current theories couple the transfer of latent heat with selection mechanisms at the interface. Measurements of succinonitrile (SCN) dendrites in microgravity show reasonable agreement between heat transfer predictions and experiment. However, data and analysis for assessing interfacial physics theories are less definitive. We are studying, and will present data on, transient effects in dendritic growth of SCN. We employ the Clapeyron pressure/melting temperature effect to make a rapid change in a sample's hydrostatic pressure, and thereby rapidly change the specimen's melting temperature, forcing the dendrite to select a new steady-state. These initial measurements show some surprising and non-intuitive effects.

Koss, M. B.↗

NASA NTRS2003

Combustion Of Moving Droplets And Of Droplets Suspended Within A Convective Environment: Transient Numerical Results

The problem considered is that of a single-component liquid fuel (n-heptane) droplet undergoing evaporation and combustion in a hot, convective, low pressure, zero-gravity environment of infinite expanse. For a moving droplet, the relative velocity (U(sub infinity)) between the droplet and freestream is subject to change due to the influence of the drag force on the droplet. For a suspended droplet, the relative velocity is kept constant. The governing equations for the gas-phase and the liquid-phase consist of the unsteady, axisymmetric equations of mass, momentum, species (gas-phase only) and energy conservation. Interfacial conservation equations are employed to couple the two phases. Variable properties are used in the gas- and liquid-phase. Multicomponent diffusion in the gas-phase is accounted for by solving the Stefan-Maxwell equations for the species diffusion velocities. A one-step overall reaction is used to model the combustion. The governing equations are discretized using the finite volume and SIMPLEC methods. A colocated grid is adopted. Hyperbolic tangent stretching functions are used to concentrate grid points near the fore and aft lines of symmetry and at the droplet surface in both the gas- and liquid-phase. The discretization equations are solved using the ADI method with the TDMA used on each line of the two alternating directions. Iterations are performed within each time-step until convergence is achieved. The grid spacing, size of the computational domain and time-step were tested to ensure that all solutions are independent of these parameters. A detailed discussion of the numerical model is given.

Gogos, George↗

NASA NTRS1984

Future Directions of Research in Adhesion and Friction. Status of Understanding

Four main classes of materials are covered: metals, ceramics, elastomers and polymers. In dealing with adhesion a distinction is drawn between adhesion (which often involves thermodynamic concepts) and the observed adhesive strength, here designated as the pull-off force. With metals the level of understanding of the interfacial bond is well understood theoretically and was elegantly studied experimentally. However the role of surface roughness and especially ductility needs to be incorporated into an overall view of metallic adhesion. The effect of contaminant films has not gone far beyond purely descriptive language. The friction of metals is understood in terms of surface topography, adhesion and deformation and there are now promising solutions using slipline-field theory. These do not however, cope very satisfactorily with work hardening. For surfaces sliding at higher speeds there are now good models which deal with asperity instabilities produced by frictional heating.

Tabor, D.↗

NASA NTRS2016

Levitation Technology in International Space Station Research

The International Space Station (ISS) is a unique multidisciplinary orbiting laboratory for science and technology research, enabling discoveries that benefit life on Earth and exploration of the universe. ISS facilities for containerless sample processing in Materials Science experiments include levitation devices with specimen positioning control while reducing containment vessel contamination. For example, ESA's EML (ElectroMagnetic Levitator), is used for melting and solidification of conductive metals, alloys, or semiconductors in ultra-high vacuum, or in high-purity gaseous atmospheres. Sample heating and positioning are accomplished through electromagnetic fields generated by a coil system. EML applications cover investigation of solidification and microstructural formation, evaluation of thermophysical properties of highly reactive metals (whose properties can be very sensitive to contamination), and examination of undercooled liquid metals to understand metastable phase convection and influence convection on structural changes. MSL utilization includes development of novel light-weight, high-performance materials. Another facility, JAXA's ELF (Electrostatic Levitation Furnace), is used to perform high temperature melting while avoiding chemical reactions with crucibles by levitating a sample through Coulomb force. ELF is capable of measuring density, surface tension, and viscosity of samples at high temperatures. One of the initial ELF investigations, Interfacial Energy-1, is aimed at clarification of interfacial phenomena between molten steels and oxide melts with industrial applications in control processes for liquid mixing. In addition to these Materials Science facilities, other ISS investigations that involve levitation employ it for biological research. For example, NASA's "Magnetic 3D Culturing and Bioprinting" investigation uses magnetic levitation for three-dimensional culturing and positioning of magnetized cells to generate spheroid assemblies for biomedical applications. Levitation is also used as a modeled microgravity ground analog in the NASA OsteoOmics ISS investigation, which tests whether magnetic levitation accurately simulates microgravity conditions by studying gravitational regulation of osteoblast and osteoclast genomics and metabolism. Elucidating the cellular mechanisms of bone loss in microgravity contributes to the understanding of bone loss in medical disorders on Earth, which may lead to development of preventive or therapeutic countermeasures. Thus, the ISS state-of-the-art laboratory offers various levitation capability platforms with applications for innovative research in Materials and Life Sciences disciplines, with benefits for humanity.

Guinart-Ramirez, Y.↗

DOE OSTI2026

Heteroepitaxial control of thickness, strain, and domain architecture in few-layer ferroelectric tin monochalcogenides

Thin-film epitaxy and epitaxial strain have been widely exploited to tune domain configurations, switching behavior, and ferroic properties in conventional three-dimensional ferroelectric thin films; however, its application to controlling the properties of two-dimensional (2D) ferroelectrics has remained largely unexplored. Here, using SnX (X = Se, S) as a model system, we demonstrate heteroepitaxial control of thickness, strain state, and domain architecture in few-layer ferroelectric SnX via growth on monolayer MoS2 van der Waals (vdW) templates. Compared with conventional growth, MoS2-templated heteroepitaxy promotes epitaxial alignment, yielding ultrathin SnSe films with improved crystalline quality, full areal coverage, and enlarged lateral dimensions. Strong interfacial epitaxial coupling induces pronounced in-plane strain and stabilizes a hierarchical ferroelastic domain architecture, in which long-range 90° stripe domains are further subdivided into nanoscale rotational variants, as revealed by scanning transmission electron microscopy and synchrotron X-ray microscopy. Piezoresponse force microscopy and second-harmonic polarimetry confirm robust in-plane polarization, while ferroelectricity in SnSe is established through polarization-electric field hysteresis and nonvolatile ferroelectric resistive switching with on/off ratios approaching 1000. A nonvolatile, switchable ferroelectric diode effect further evidences direct coupling between polarization and charge transport. Notably, ferroelectric switchability exhibits a strong thickness dependence and is preserved only below ∼10 layers. This vdW heteroepitaxial strategy is further extended to ferroelectric SnS. The seamless heteroepitaxial integration of 2D ferroelectrics with CMOS-compatible, wafer-scale MoS2 templates provides a general and scalable route for strain-enabled structural and ferroelectric engineering in emerging memory and low-power optoelectronic applications.

Wang, Yueyin↗

DOE OSTI2021

Learning intermolecular forces at liquid–vapor interfaces

By adopting a perspective informed by contemporary liquid-state theory, we consider how to train an artificial neural network potential to describe inhomogeneous, disordered systems. Here, we find that neural network potentials based on local representations of atomic environments are capable of describing some properties of liquid-vapor interfaces but typically fail for properties that depend on unbalanced long-ranged interactions that build up in the presence of broken translation symmetry. These same interactions cancel in the translationally invariant bulk, allowing local neural network potentials to describe bulk properties correctly. By incorporating explicit models of the slowly varying long-ranged interactions and training neural networks only on the short-ranged components, we can arrive at potentials that robustly recover interfacial properties. We find that local neural network models can sometimes approximate a local molecular field potential to correct for the truncated interactions, but this behavior is variable and hard to learn. Generally, we find that models with explicit electrostatics are easier to train and have higher accuracy. We demonstrate this perspective in a simple model of an asymmetric dipolar fluid, where the exact long-ranged interaction is known, and in an ab initio water model, where it is approximated.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1986

Ice banding as a response of the coupled ice-ocean system to temporally varying winds

This study models formation of ice bands in the marginal ice zones. A one-dimensional coupled ice-ocean model is used in which the ice model is coupled to a reduced gravity ocean model through interfacial stresses. The internal ice stresses are important only at high ice concentrations (90-100 percent); otherwise, the main balance for the ice motion is between the air-ice stress and the ice-water stress, i.e., free drift. The drag coefficients were chosen so that the air-ice momentum flux is 3 times greater than the air-ocean momentum flux. Thus the Ekman transport is larger under the ice than in the open water, so that winds parallel to the ice edge, with the ice on the right, produce upwelling. The upwelling simulation was extended to include temporally varying forcing, which was chosen to vary sinusoidally with a 4-day period. This forcing resembles successive cyclone passings perpendicular to the ice edge. When the oceanic upper layer was thin, which means that the dynamics are strongly nonlinear, the ice bands were formed. The up/downwelling signals do not disappear in wind reversals because of nonlinear advection. This leads to convergences and divergences in oceanic and ice velocities that manifest themselves as ice banding. At least one wind reversal is needed to produce one ice band.

Hakkinen, S.↗

NASA NTRS1983

Motion driven by the interface

Due to the reduction in buoyant forces aboard orbiting spacecraft such as the Space Shuttle, fluid motion driven by gradients in interfacial tension will be important in the processing of materials in space. In this paper, preliminary results from a study of surface tension driven flow in a pendant drop are reported. The drop is heated from above, and the resulting temperature gradients on the drop surface give rise to interfacial tension gradients. These, in turn, drive a circulation in the drop which is made visible by suitable tracers. The velocities are measured using a video technique, and the data on core velocities are found to agree well with results from a predictive theoretical model.

Jayaraj, K.↗

DOE OSTI2023

Quantifying interfacial energetics of 2D semiconductor electrodes using in situ spectroelectrochemistry and many-body theory

Hot carrier extraction occurs in 2D semiconductor photoelectrochemical cells. Boosting the energy efficiency of hot carrier-based photoelectrochemical cells requires maximizing the hot carrier extraction rate relative to the cooling rate. One could expect to tune the hot carrier extraction rate constant (k ET ) via a Marcus–Gerischer relationship, where k ET depends exponentially on ΔG°' (the standard driving force for interfacial electron transfer). ΔG°' is defined as the energy level difference between a semiconductor's conduction/valence band (CB/VB) minima/maxima and the redox potential of reactant molecules in solution. A major challenge in the electrochemistry community is that conventional approaches to quantify ΔG°' for bulk semiconductors (e.g., Mott–Schottky measurements) cannot be directly applied to ultrathin 2D electrodes. The specific problem is that enormous electronic bandgap changes (>0.5 eV) and CB/VB edge movement take place upon illuminating or applying a potential to a 2D semiconductor electrode. Here, we develop an in situ absorbance spectroscopy approach to quantify interfacial energetics of 2D semiconductor/electrolyte interfaces using a minimal many-body model. Our results show that band edge movement in monolayer MoS 2 is significant (0.2–0.5 eV) over a narrow range of applied potentials (0.2–0.3 V). Such large band edge shifts could change k ET by a factor of 10–100, which has important consequences for practical solar energy conversion applications. We discuss the current experimental and theoretical knowledge gaps that must be addressed to minimize the error in the proposed optical approach.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Modeling Equilibrium Solid–Liquid Interfaces under Effective Constant Chemical Potential Using Machine Learning Interatomic Potentials

The chemical potential (μ) of species in solution is essential for understanding various chemical processes at interfaces. Molecular dynamics (MD) simulations, constrained by fixed compositions, cannot maintain constant chemical potential with reference to a targeted concentration or chemical potential under nonequilibrium or dynamic conditions, as solute species can migrate to the interface and deplete (or enrich) the bulk due to solute-interface interactions. In this study, we introduce a simple and computationally efficient approach named iterative quasi-constant chemical potential molecular dynamics (iqCμMD) simulation, which helps simulate targeted molar concentrations of species in solution. iqCμMD overcomes the limitations of conventional MD by adjusting the number of species in the solution to reach a target bulk concentration (chemical potential), which allows simulation of the interface under the bulk conditions comparable to experiment. We demonstrate our approach using machine learning interatomic potential (MLIP)-based MD simulations of the Na 2 SO 4,aq –graphene interface, and to show the transferability of our approach, we also perform classical force field-based MD simulations of NaCl aq –air and NaCl aq –graphite interfaces, which produce comparable results to previous CμMD simulations. Our results also show that the iqCμMD approach efficiently achieves the desired bulk ion concentration within two iterations, and by utilizing MLIPs, we can achieve converged results using relatively small-scale simulations compared to previous CμMD simulations. By combining iqCμMD with MLIP-driven simulations, solid–liquid interfaces can be modeled under an effective constant chemical potential with DFT-level accuracy. Here, we show that iqCμMD offers a robust and simple computational framework for constant chemical potential simulations, as its only requirement is to be able to converge interfacial simulations with a measurable bulk region.

Chemical structure↗

NASA NTRS1983

A tire contact solution technique

An efficient method for calculating the contact boundary and interfacial pressure distribution was developed. This solution technique utilizes the discrete Fourier transform to establish an influence coefficient matrix for the portion of the pressurized tire surface that may be in the contact region. This matrix is used in a linear algebra algorithm to determine the contact boundary and the array of forces within the boundary that are necessary to hold the tire in equilibrium against a specified contact surface. The algorithm also determines the normal and tangential displacements of those points on the tire surface that are included in the influence coefficient matrix. Displacements within and outside the contact region are calculated. The solution technique is implemented with a finite-element tire model that is based on orthotropic, nonlinear shell of revolution elements which can respond to nonaxisymmetric loads. A sample contact solution is presented.

Tielking, J. T.↗

NASA NTRS1994

Interfacial force field characterization of a constrained vapor bubble thermosyphon using IAI

The isothermal profiles of the extended meniscus in a quartz cuvette were measured in a gravitational field using IAI (image analyzing interferometer) which is based on computer enhanced video microscopy of the naturally occurring interference fringes. The experimental results for heptane and pentane menisci were analyzed using the extended Young-Laplace Equation. These isothermal results characterized the interfacial force field in-situ at the start of the heat transfer experiments by quantifying the dispersion constant for the specific liquid-solid system. The experimentally obtained values of the disjoining pressures and the dispersion constants are compared to the subsequent non-isothermal experiments because one of the major variables in the heat sink capability of the CVBT is the dispersion constant. In all previous studies of micro heat pipes the value of the dispersion constant has been 'guesstimated'. The major advantages of the current glass cell is the ability to view the extended meniscus at all times. Experimentally, we find that the extended Young-Laplace Equation is an excellent model for for the force field at the solid-liquid vapor interfaces.

Dasgupta, Sunando↗

NASA NTRS1990

Free convection in the Matian atmosphere

The 'free convective' regime for the Martian atmospheric boundary layer (ABL) was investigated. This state occurs when the mean windspeed at the top of the ABL drops below some critical value U(sub c) and positive buoyant forces are present. Such forces can arise either from vertical temperature or water vapor gradients across the atmospheric surface layer. During free convection, buoyant forces drive narrow plumes that ascend to the inversion height with a return circulation consisting of broad slower-moving downdraughts. Horizontal pressure, temperature, windspeed, and water vapor fluctuations resulting form this circulation pattern can be quite large adjacent to the ground (within the surface layer). The local turbulent fluctuations cause non-zero mean surface stresses, sensible heat fluxes, and latent heat fluxes, even when the mean regional windspeed is zero. Although motions above the surface layer are insensitive to the nature of the surface, the sensible and latent heat fluxes are primarily controlled by processes within the interfacial sublayer immediately adjacent to the ground during free convection. Thus the distinction between aerodynamically smooth and rough airflow within the interfacial sublayer is more important than for the more typical situation where the mean regional windspeed is greater than U(sub c). Buoyant forces associated with water vapor gradients are particularly large on Mars at low pressures and high temperatures when the surface relative humidity is 100 percent, enhancing the likelihood of free convection under these conditions. On this basis, Ingersol postulated the evaporative heat losses from an icy surface on Mars at 237 K and current pressures would exceed the available net radiative flux at the surface, thus prohibiting ice from melting at low atmospheric pressures. Schumann has developed equations describing the horizontal fluctuations and mean vertical gradients occurring during free convection. Schumann's model was generalized to include convection driven by water vapor gradients and to include the effects of circulation above both aerodynamically smooth and rough surfaces.

Clow, G. D.↗

DOE OSTI2022

Electrochemically-Induced Phase Transformations in Battery Storage Compounds (Final Technical Report)

Compounds of interest for ion storage in advanced batteries frequently exhibit phase transformations as the working ion concentration varies. Under large electrochemical driving forces inherent to practical use, systems are often driven far from equilibrium. This program combines experiments and theory to understand the phase transition behavior of ion insertion compounds when electrochemically driven far from equilibrium. As model systems, we focus on alkaline metal phosphates AMPO4 (A = alkali; M = first row transition metal) of olivine structure, which are both technologically interesting and ideally suited for fundamental study due to the ability to systematically tune transformation strain, and along with it, the phase transformation pathway. Behavior in compositions having large transformation strains (~15 vol%) requiring plasticity for strain accommodation is emphasized. Experimental techniques include operando characterization of structure while simultaneously varying electrokinetic parameters, and high resolution microscopy of nanoscale and interfacial phenomena. Phase-field modeling is used to model the thermodynamics and kinetics of competing transformation pathways, extended to include the effects of plasticity, and integrated with porous electrode kinetic theory to treat multi-particle effects. Success in this project will lead to an ability to design ion storage compounds with predictable transformation pathways, electrochemical kinetics, capacity utilization, and durability. New technologically important compounds may also be discovered.

36 MATERIALS SCIENCE↗