Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interfacial forces”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Thermal buoyancy and Marangoni convection in a two fluid layered system - A numerical study

Thermal buoyancy and surface tension driven convection is numerically investigated in a system with two, immiscible fluids. The geometry investigated, has an open cavity configuration with the lighter fluid situated on top of the heavier fluid, forming a stable layered system. The upper fluid meniscus and the interface are assumed to be flat and undeformable in the calculations. The governing equations and boundary/interface conditions are solved by a control volume based finite difference scheme for two pairs of immiscible fluids; the water-hexadecane system and a so called generic system. The steady state calculations show, that dramatically different flows are predicted when the interfacial tension effects are included or excluded from the system model. These differences are particularly appreciable in surface tension dominated flows, that are typical of microgravity situations. Complex flow patterns, with induced secondary flows are noticed in both the fluids. In general, the overall system heat transfer is found to increase with increases in the thermal buoyancy and surface tension effects, but the behavior of the system flow and thermal fields is not easy to characterize, when different combinations of these forces are considered.

Ramachandran, N.↗
DOE OSTI2025

Binder Effects on Processing, Mechanical Properties, and Performance of Thin Sulfide Solid‐State Electrolytes

All‐solid‐state batteries (ASSBs) offer enhanced safety and energy density compared to conventional lithium‐ion batteries by replacing flammable liquid electrolytes with solid‐state electrolytes (SSEs). Among SSEs, sulfide‐based electrolytes exhibit high ionic conductivity and mechanical deformability, making them promising candidates for next‐generation energy storage. However, their practical implementation is hindered by interfacial instability, mechanical brittleness, and challenges in fabricating ultrathin electrolyte membranes (<30 μm) with robust mechanical integrity. Here, this study systematically examines the influence of polymeric binders—polyisobutylene, hydrogenated nitrile butadiene rubber, and styrene–ethylene–butylene–styrene (SEBS)—on the structural, mechanical, and electrochemical performance of thin sulfide SSE membranes. Key findings reveal that SEBS enables the fabrication of ultrathin, uniform membranes, while binder elasticity significantly affects structural stability during cycling. Operando stack pressure measurements indicate that binder properties directly influence adhesion force for LPSCl particles, influencing their stabilizing cycling period. These results underscore the critical role of polymer binders beyond mechanical reinforcement, positioning them as essential design variables in sulfide SSE engineering. By linking binder chemistry to processability and electrochemical performance, this study provides insights into optimizing sulfide SSEs, advancing their commercial viability in ASSBs.

argyrodite sulfide electrolyte↗
DOE OSTI2025

Directing the adsorption and assembly of laponite nano-discs at oil–water interfaces

Achieving tunable controls on the adsorption and self-assembly of nanoscale building blocks at immiscible fluid interfaces is essential for synthesizing advanced materials, stabilizing emulsions, and the sustainable storage and recovery of fluids. Although extensive efforts have been directed toward resolving the assembly of spherical nanoparticles at water–hydrocarbon interfaces, 2D nanoparticles have received far less attention. In this study, we developed novel controls to direct the adsorption and assembly of 2D laponite nano-discs (diameter = ∼30 nm and thickness = ∼1 nm) at water–heptane interfaces using traces of sodium dodecyl sulfate (SDS) surfactant, demonstrated by advanced in situ small-angle X-ray scattering, atomic force microscopy, spinning drop tensiometer measurements and molecular dynamics simulations. The results show that SDS surfactant displaces the adsorbed laponite nano-discs from the interface toward the aqueous phase. The extent of this displacement increases with SDS concentration such that high SDS concentrations convert laponite-rich interfaces to SDS-rich interfaces free of laponite nano-discs. This transformation is associated with significant alterations in the rheological and nanomechanical properties of the host interface. In this context, reduction in interfacial tension and interfacial stiffness and an increase in the interfacial deformation is observed on increasing the SDS concentration from 0 to 0.1 wt%. The displacement of laponite nano-discs from the interface is driven by strong electrostatic interactions between the hydrophilic SDS group and interfacial water molecules. These studies unlock new insights into the adsorption and assembly of 2D nanoscale particles at water–hydrocarbons interfaces that are relevant for various applications related to energy and environmental science and engineering.

interfacial tension↗
NASA NTRS1971

Transient-forced convection film boiling on an isothermal flat plate

An approach for the solution of transient-forced convection film boiling on an isothermal flat plate using the boundary layer model is developed. The similarity variables are used to convert the governing partial differential equations to ordinary ones. The results of numerical solutions of these ordinary equations indicate that the transient process can be classified as one-dimensional conduction, intermediate, and the steady state regions. The time required for the one-dimensional conduction and the time necessary to attain a steady state condition are obtained. The influence of interfacial shear is seen to be negligible while the Prandtl Number and the ratio (C sub p delta T divided by h sub fg times Pr) have major influence. The use of local similarity approximations for the intermediate regime facilitates prediction of complete boundary layer growth. Using the ratio of time at any instant to the steady state time as abscissa, the curves representing the boundary layer growth can be merged into a single mean curve within 5 percent. Further, the analysis shows that the average rate of heat transfer during transient is 50 to 100 percent higher than those at steady state. The average rate of vapor convected away is 10 to 15 percent lower than at steady state while the average rate of accumulation to form the vapor layer is 1 to 14 times larger. Further, the total heat transferred during transient increases and the evaporation decreases for increasing values of C sub p delta T divided by h sub fg times Pr

Nagendra, H. R.↗
NASA NTRS1984

The role of material properties in adhesion

When two solid surfaces are brought into contact strong adhesive bond forces can develop between the materials. The magnitude of the forces will depend upon the state of the surfaces, cleanliness and the fundamental properties of the two solids, both surface and bulk. Adhesion between solids is addressed from a theoretical consideration of the electronic nature of the surfaces and experimentally relating bond forces to the nature of the interface resulting from solid state contact. Surface properties correlated with adhesion include, atomic or molecular orientation, reconstruction and segregation as well as the chemistry of the surface specie. Where dissimilar solids are in contact the contribution of each is considered as is the role of their interactive chemistry on bond strength. Bulk properties examined include elastic and plastic behavior in the surficial regions, cohesive binding energies, crystal structure, crystallographic orientation and state. Materials examined with respect to interfacial adhesive interactions include metals, alloys, ceramics, polymers and diamond. They are reviewed both in single and polycrystalline form. The surfaces of the contacting solids are studied both in the atomic or molecularly clean state and in the presence of selected surface contaminants.

Buckley, D. H.↗
NASA NTRS1993

Wave effects on ocean-ice interaction in the marginal ice zone

The effects of wave train on ice-ocean interaction in the marginal ice zone are studied through numerical modeling. A coupled two-dimensional ice-ocean model has been developed to include wave effects and wind stress for the predictions of ice edge dynamics. The sea ice model is coupled to the reduced-gravity ocean model through interfacial stresses. The main dynamic balance in the ice momentum is between water-ice stress, wind stress, and wave radiation stresses. By considering the exchange of momentum between waves and ice pack through radiation stress for decaying waves, a parametric study of the effects of wave stress and wind stress on ice edge dynamics has been performed. The numerical results show significant effects from wave action. The ice edge is sharper, and ice edge meanders form in the marginal ice zone owing to forcing by wave action and refraction of swell system after a couple of days. Upwelling at the ice edge and eddy formation can be enhanced by the nonlinear effects of wave action; wave action sharpens the ice edge and can produce ice meandering, which enhances local Ekman pumping and pycnocline anomalies. The resulting ice concentration, pycnocline changes, and flow velocity field are shown to be consistent with previous observations.

Liu, Antony K.↗
NASA NTRS1998

Numerical Simulation of Combustion and Extinction of a Solid Cylinder in Low-Speed Cross Flow

The combustion and extinction behavior of a diffusion flame around a solid fuel cylinder (PMMA) in low-speed forced flow in zero gravity was studied numerically using a quasi-steady gas phase model. This model includes two-dimensional continuity, full Navier Stokes' momentum, energy, and species equations with a one-step overall chemical reaction and second-order finite-rate Arrhenius kinetics. Surface radiation and Arrhenius pyrolysis kinetics are included on the solid fuel surface description and a parameter Phi, representing the percentage of gas-phase conductive heat flux going into the solid, is introduced into the interfacial energy balance boundary condition to complete the description for the quasi-steady gas-phase system. The model was solved numerically using a body-fitted coordinate transformation and the SIMPLE algorithm. The effects of varying freestream velocity and Phi were studied. These parameters have a significant effect on the flame structure and extinction limits. Two flame modes were identified: envelope flame and wake flame. Two kinds of flammability limits were found: quenching at low-flow speeds due to radiative loss and blow-off at high flow speeds due to insufficient gas residence time. A flammability map was constructed showing the existence of maximum Phi above which the solid is not flammable at any freestream velocity.

Tien, J. S.↗
NASA NTRS2000

Effects of G-Jitter on Interfacial Dynamics of Two Miscible Liquids: Application of MIM

We designed an experiment to examine the effects of g-jitter on mixing of two miscible liquids using the Microgravity Vibration Isolation Mount (MIM). The global bifurcation of the interface was observed with the MIM operating alternatively to either transmit the g-jitter, isolate from the g-jitter or to provide controlled vibration levels with well defined amplitude and frequency content. With the MIM in isolation mode, the interface remains stationary indicating buoyancy induced convection is negligibly small such that mixing occurs via intrinsic mass diffusion without the masking effect of vibration driven convection. Analytical and computational results are in agreement with the experimental findings. Operation of the MIM in forced mode with conditions typical of g-jitter shows that vibration induced convective flows can excite instability mechanisms such as Kelvin-Helmholtz to generate large amplitude quasi-stationary waves oriented vertically for various cases with Stokes-Reynolds number in the range of 0.003 to 0.5. The two and four mode quasi-stationary waves are also predicted with a mathematical model. Though unplanned, the effect of a primary thruster filing was captured and shown to cause a catastrophic bifurcation, enhancing local mass transport. In light of the findings, experiments planned for the International Space Station should consider the potential effects of g-jitter.

Duval, Walter M. B.↗
DOE OSTI2022

Local conformations and heterogeneities in structures and dynamics of isotactic polypropylene adsorbed onto carbon fiber

Carbon fiber (CF) reinforced polymers (CFRPs) have experienced widespread use in various industries. One of the most important parameters that controls the macroscopic property of CFRPs is the interface between a polymer matrix and CF. There is growing evidence to suggest the formation of a bound polymer layer (BPL), i.e., polymer chains that physically adsorb on a filler surface. However, this interface is always in contact with the thicker part of a polymer matrix, rendering its understanding a difficult task. Therein, we use CF-reinforced isotactic polypropylene (iPP) as a rational CFRP. To characterize the BPL on the CF surface, we extracted it from the CFRP using solvent-rinsing with p-xylene. The physical and thermal properties of the BPL were characterized by differential scanning calorimetry and thermogravimetric analysis, while its microscopic structures and dynamics were probed by small-angle neutron scattering and quasi-elastic neutron scattering (QENS) techniques. Subsequently, we employed atomistic molecular dynamics (MD) simulations to complement the QENS results above the bulk melting temperature and reveal details that were experimentally inaccessible. In this work, we observed that the degree of crystallinity of the BPL was quite lower than the bulk, while the melting temperature of the BPL remained the same as the bulks. Within the given length and time scales probed by QENS and MD, we also observed that most of the bound chains were mobile, with the formation of a high-density region (less than 1 nm in thickness) near the CF surface. The segmental dynamics of the bound chains probed by both QENS and MD were also much faster than those of the free chains, possibly due to the presence of a free surface region at the topmost surface of the BPL. Furthermore, the MD results demonstrated that the backbone chains and side groups lie nearly flat on the CF surface, which is the driving force for the flattening process of the iPP BPL to overcome the conformational entropy loss in the total free energy.

36 MATERIALS SCIENCE↗
NASA NTRS1994

Separation of crack extension modes in composite delamination problems

This work concerns fracture mechanics modeling of composite delamination problems. In order to predict delamination resistance, an applied stress intensity factor, K, or energy release rate, G, must be compared to a mode-dependent critical value of K or G from experiment. In the interfacial fracture analysis of most applications and some tests, the mode of crack extension is not uniquely defined. It is instead a function of distance from the crack tip due to the oscillating singularity existing at the tip. In this work, a consistent method is presented of extracting crack extension modes in such cases. In particular, use of the virtual crack closure technique (VCCT) to extract modes of crack extension is studied for cases of a crack along the interface between two in-plane orthotropic materials. Modes of crack extension extracted from oscillatory analyses using VCCT are a function of the virtual crack extension length, delta. Most existing efforts to obtain delta-independent modes of crack extension involve changing the analysis in order to eliminate its oscillatory nature. One such method involves changing one or more properties of the layers to make the oscillatory exponent parameter, epsilon, equal zero. Standardized application of this method would require consistent criteria for identifying which properties can be altered without changing the physical aspects of the problem. Another method involves inserting a thin homogeneous layer (typically referred to as a resin interlayer) along the interface and placing the crack within it. The drawbacks of this method are that it requires increased modeling effort and introduces the thickness of the interlayer as an additional length parameter. The approach presented here does not attempt to alter the interfacial fracture analysis to eliminate its oscillatory behavior. Instead, the argument is made that the oscillatory behavior is non-physical and that if its effects were separated from VCCT quantities, then consistent, delta-independent modes of crack extension could be defined. Knowledge of the near-tip fields in a planar orthotropic material interfacial fracture analysis is used to determine the explicit delta dependence of VCCT parameters. Once this delta dependence is determined, energy release rates are defined with this delta dependence factored out. This modified VCCT method is applied to results from two finite element test cases. It is shown that, as predicted, delta-independent modes of crack extension result. The modified VCCT approach shows potential as a consistent method of extracting crack extension modes. It uses the same information from a finite element analysis (i.e., nodal forces and displacements) as the traditional VCCT method does. The A-independent modes extracted using the modified VCCT approach can also be used as guides to test the convergence of finite element solutions.

Beuth, Jack L.↗
DOE OSTI2022

A generalized 3D elastic model for nanoscale, self-assembled oxide-metal thin films with pillar-in-matrix configurations

In recent years, functional oxide-metal based vertically aligned nanocomposite (VAN) thin films have gained interest due to their intriguing physical properties and multifunctionalities stemming from the complex interactions between the two phases in the film and the substrate. In this work, we develop a model for studying the energetics of these thin film systems, including the effects of both lattice mismatch and capillary forces due to interface curvature. Each phase is incorporated into the model using a phase indicator function, and we introduce the capillary forces as body forces using a vector density representation of the interface. The model is implemented using the finite element method to study the deformation of the thin film which is composed of Au nanopillars embedded in a La 0.7 Sr 0.3 MnO 3 (LSMO) matrix on an SrTiO 3 (STO) substrate. Further, the results suggest that the total energy is lowest for random configurations of pillars compared to ordered square and hexagonal lattice configurations, consistent with the random distribution of pillars found in experiments. Furthermore, we find that the interfacial energy dominates the total energy of each configuration, suggesting that interfacial energy in the system is an important design parameter for nanocomposite growth, along with the lattice mismatch.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Geometry and dynamics of lipid membranes: The Scriven-Love number

The equations governing lipid membrane dynamics in planar, spherical, and cylindrical geometries are presented here. Unperturbed and first-order perturbed equations are determined and nondimensionalized. In membrane systems with a nonzero base flow, perturbed in-plane and out-of-plane quantities are found to vary over different length scales. Here, a new dimensionless number, named the Scriven-Love number, and the well-known Föppl–von Kármán number result from a scaling analysis. The Scriven-Love number compares out-of-plane forces arising from the in-plane, intramembrane viscous stresses to the familiar elastic bending forces, while the Föppl–von Kármán number compares tension to bending forces. Both numbers are calculated in past experimental works, and span a wide range of values in various biological processes across different geometries. In situations with large Scriven-Love and Föppl–von Kármán numbers, the dynamical response of a perturbed membrane is dominated by out-of-plane viscous and surface tension forces—with bending forces playing a negligible role. Calculations of non-negligible Scriven-Love numbers in various biological processes and in vitro experiments show in-plane intramembrane viscous flows cannot generally be ignored when analyzing lipid membrane behavior.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Mass transport in a highly immiscible alloy on extended shear deformation

Forced mixing to a single phase or supersaturated solid solution (SSS) and its prerequisite microstructure evolution in immiscible systems has been a focus of research for both fundamental science and a variety of applications. Controlling the formation of SSS by shear deformation assisted processing could enable a material design beyond conventional equilibrium microstructure in immiscible systems. Here, a highly immiscible (mixing enthalpy of ~ 20 kJ·mol -1 ) Cu-50 at. % Cr binary alloy was employed to investigate the microstructure evolution and localized tendencies of SSS during severe shear deformation. The present results demonstrate distinctive microstructure refinement process in each phase and how they lead to localized SSS as a function of shear strain. Preferential dynamic recrystallization occurs in the softer Cu phase due to strain localization, leading to substantial grain refinement. The refinement of Cr phase in the top-most layer, however, is enabled by the progressive evolution of grain lamination, splitting, and fragmentation as a function of shear strain. The eventual SSS is found to be strongly dependent on the local environments that affect the dislocation activity, including the level of microstructure refinement, the interfacial orientation relationship, hardness difference, and supposed stability of oxidation. Ab initio simulations comfirm that it is more favorable to oxidize Cr than Cu at incoherent Cu/Cr interfaces which then limits the mass transport on an incoherent boundary. Our results shed light on the underpinning mechanism for non-equilibrium mass transport in immiscible systems upon severe deformation that can be applicable to a variety of processing techniques aimed at producing immiscible alloys with superior mechanical properties.

36 MATERIALS SCIENCE↗
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↗
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 OSTI2021

Effect of the Segmental Structure of Thermoplastic Polyurethane (Hardness) on the Interfacial Adhesion of Textile-Grade Carbon Fiber Composites

In polymer composites, the fiber–matrix interface is primarily influenced by the surface treatment of the fibers and polymer morphology. Previous studies have investigated the effect of surface treatment of carbon fiber on the mechanical properties of the resulting composites. However, very few studies have explored the chemical structural modification of polymer effect on the fiber–matrix adhesion. In this work, the interfaces of soft thermoplastic polyurethane (S-TPU) and hard-segmented TPU (H-TPU) were investigated through surface, thermal, and mechanical characterization. A textile-grade carbon fiber (TCF) with 1% concentration of epoxy sizing (an emerging material for nonaerospace applications with 450 K filament tows) was used as a reinforcement to investigate the structure–property relationship at the interface. Atomic force microscopy results showed 39% higher surface roughness for S-TPU than for H-TPU. X-ray photoelectron spectroscopy results revealed a 550% increase in C=O content, which can provide multiple hydrogen bonding networks in H-TPU, and these C=O bonds can produce a strong chemical bond at the fiber–matrix interface. Additionally, dynamic mechanical analysis and differential scanning calorimetry results confirmed the presence of hydrogen bonding that enhances the cross-linked density by 232% in H-TPU compared to S-TPU. The improved mechanical properties of H-TPU composites, such as flexural, impact and tensile by 30, 50, and 130% compared to S-TPU composites, prove that the crystallinity and hydrogen bonding significantly increase the load bearing capacity due to the strong interface. The mechanical properties of TCF–TPU composites validate the formation of chemical bonds through a nucleophilic addition reaction at the interface and improve the bond strength of the composites. Thus, tailoring the polyurethane structure broadens the performance of segmented TPU in conjunction with TCF reinforcement, which has implications for cost-effective, high-strength composites.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Computational thermomechanics for crystalline rock. Part II: Chemo-damage-plasticity and healing in strongly anisotropic polycrystals

We present a thermal–mechanical–chemical-phase field model that captures the multi-physical coupling effects of precipitation creeping, crystal plasticity, anisotropic fracture, and crack healing in polycrystalline rock at various temperature and strain-rate regimes. This model is solved via a fast Fourier transfer solver with an operator-split algorithm to update displacement, temperature and phase field, and chemical concentration incrementally. In nuclear waste disposal in salt formation, brine inside the crystal salt may migrate along the grain boundary and cracks due to the gradient of interfacial energy and pressure. This migration has a significant implication on the permeability evolution, creep deformation, and crack healing within rock salt but is difficult to incorporate implicitly via effective medium theories compared with computational homogenization. As such, we introduce a thermodynamic framework and a corresponding computational implementation that explicitly captures the brine diffusion along the grain boundary and crack at the grain scale. Meanwhile, the anisotropic fracture and healing are captured via a high-order phase field that represents the regularized crack region in which a newly derived non-monotonic driving force is used to capture the fracture and healing due to the solution–precipitation. Numerical examples are presented to demonstrate the capacity of the thermodynamic framework to capture the multiphysics material behaviors of rock salt.

42 ENGINEERING↗