Engineering Papers⌕ Search

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.

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

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Development of machine learning framework for interface force closures based on bubble tracking data

Interfacial force closures in the two-fluid model play a critical role for the predictive capabilities of void fraction distribution. However, the practices of interfacial force modeling have long been challenged by the inherent physical complexity of the two-phase flows. The rapidly expanding computational capabilities in the recent years have made high-fidelity data from the interface-captured direct numerical simulation become more available, and hence potential for data-driven interfacial force modeling has prevailed. In this work, we established a data-driven modeling framework integrated to the HZDR multiphase Eulerian-Eulerian framework for computational fluid dynamics simulations. The data-driven framework is verified in a benchmark problem, where a feedforward neural network managed to capture the non-linear mapping between bubble Reynolds number and drag coefficient and reproduce the void distribution resulting from the baseline model in the test case. The second focus is on utilizing the bubble tracking data set to form a closure for the bubble drag in the turbulent bubbly flow, in which the drag coefficient is set to be correlated with the bubble Reynolds number and the Eötvös number. Pseudo-steady state filtering in the Frenet Frame was carried out to obtain the drag coefficient from the turbulent bubbly flow data. The performance of the data-driven drag model is also examined through a case study, where improvement of model’s prediction near-wall is regarded necessary. In conclusion, discussion and further plans of investigation are provided.

42 ENGINEERING↗
DOE OSTI2024

Bespoke Liquid/Liquid Interfaces (Final Technical Report)

The goal of DE-SC0001815 was to advance the basic science of liquid:liquid interface formation, to develop a deeper understanding of the mechanisms of phase separation and the essential relationships between solution composition, organization and dynamics that underlie the kinetic regime of solvent extraction. This included learning how interfacial organization and dynamics alters the properties of the primary coordination sphere of ions and the free energy of transport of ions complexes across a phase boundary. We relied primarily upon classical molecular dynamics studies to determine the equilibrium ensembles of these complex systems, but also utilized ab-initio MD and cluster-based density functional theory (DFT) calculations when more detailed investigation of the electronic structure was needed. We continued development of graph-theory based analyses to elucidate hierarchical correlations and expanded into geometric topology methods to quantify the collectively organized structures that can organize at a liquid/liquid interface during solute transport. One of the main conclusions was from the observation of two distinct mechanisms for solute transport - those that derive from amplifications of interfacial heterogeneity and surface roughness, and those wherein surface roughness has been dampened and instead collectively organized macrostructures work to bring solutes into the organic phase. It was our aim to create a concrete chemical model of the underlying driving forces behind interfacial primary and secondary structure formation and to map out the energetic features of solute transport so that tailored liquid/liquid can be developed that have characteristic kinetic features associated with mass transport.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Uncertainty Quantification for Multiphase Computational Fluid Dynamics Closure Relations with a Physics-Informed Bayesian Approach

Multiphase Computational Fluid Dynamics (MCFD) based on the two-fluid model is considered a promising tool to model complex two-phase flow systems. MCFD simulation can predict local flow features without resolving interfacial information. As a result, the MCFD solver relies on closure relations to describe the interaction between the two phases. Those empirical or semi-mechanistic closure relations constitute a major source of uncertainty for MCFD predictions. In this paper, we leverage a physics-informed uncertainty quantification (UQ) approach to inversely quantify the closure relations’ model form uncertainty in a physically consistent manner. This proposed approach considers the model form uncertainty terms as stochastic fields that are additive to the closure relation outputs. Combining dimensionality reduction and Gaussian processes, the posterior distribution of the stochastic fields can be effectively quantified within the Bayesian framework with the support of experimental measurements. As this UQ approach is fully integrated into the MCFD solving process, the physical constraints of the system can be naturally preserved in the UQ results. Here, in a case study of adiabatic bubbly flow, we demonstrate that this UQ approach can quantify the model form uncertainty of the MCFD interfacial force closure relations, thus effectively improving the simulation results with relatively sparse data support.

42 ENGINEERING↗
DOE OSTI2025

Structure and Dynamics of CO 2 at the Air–Water Interface from Classical and Neural Network Potentials

The accurate description of the structure and dynamics of CO 2 at the instantaneous air–water interface, along with the effects of surface fluctuations on the CO 2 -transport processes, is essential for the development of negative emission technologies aimed at minimizing climate change. In this study, we performed molecular dynamics simulations of CO 2 at the air–water interface using neural network potentials (NNPs) trained on ab initio data generated through density-functional-theory-based molecular dynamics simulations. We compared these results with classical force fields to assess their performance in modeling interfacial CO 2 behavior. Our findings revealed that the asymmetric interactions, coupled with thermal surface fluctuations at the air–water interface, significantly influence CO 2 transport into the aqueous phase. The simulations demonstrate that classical force fields underestimate both the free energy of CO 2 transport and the strength of its interactions at the interface compared with the neural network potentials. In conclusion, the free energy and the interfacial dynamics of CO 2 are primarily influenced by the distribution of water within the instantaneous interfacial water layer, responsible for creating an asymmetric intermolecular interaction environment within the interfacial region.

Ab initio molecular dynamics↗
DOE OSTI2026

Mapping Interfacial Solution Forces that Drive Fluorescent Nanoparticle Incorporation into Crystals

Interfacial solution structure governs processes ranging from catalytic and electrochemical reactions to particle aggregation and crystallization. Here, we design bio-inspired nanocomposites by mapping interfacial forces and controlling nanoparticle incorporation during crystal growth. Our analytical model contains measured kinetic barriers for surface approaches and equilibrium binding constants with the crystal surface. We validate this model using fluorescent silica nanoparticles and calcite. Our results show that surface chemistry dictates incorporation pathways: methoxy groups are least favorable, hydroxyls reduce kinetic barriers by interacting favorably with hydration layers, carboxylates bind strongly but must overcome a kinetic barrier, and amines outperform all kinetically and thermodynamically. This framework resolves the interplay between interfacial solution structures, particle forces and dynamics, and growth kinetics, enabling predicative control of nanocomposite formation and multiplexed mixed-particle systems.

Zhang, Mingyi [Pacific Northwest National Laborato↗
DOE OSTI2026

Validated Reactive Force Field Quantifies MXene Interfacial Properties, Mechanics, and Thermal Transport

MXenes combine rich surface chemistry, mechanical strength, and high conductivity for a multitude of emerging applications. Predictive modeling supports accelerated materials designs and has been limited by the absence of validated and transferable force fields. Here, we introduce an interpretable, reactive INTERFACE force field (IFF and IFF-R) for Ti 3 C 2 T x MXenes that is trained based on chemical knowledge and achieves quantitative agreement with experiments across lattice parameters (<0.5%), density (<0.2%), liquid contact angles, Raman spectra, and the in-plane elastic modulus (∼320 GPa). The models cover surface terminations from hydroxyl (−OH) to fluorine (−F) groups and are extensible to other chemistries. We introduce pH-resolved surface chemistry and identify dopamine adsorption mechanisms at MXene–aqueous interfaces supported by QCM-D and UV–Vis experiments. The data reveal coplanar and perpendicular binding modes and concentration-dependent multilayer assembly. We predict previously inaccessible properties, including termination-dependent cleavage energies, interlayer shear moduli and dynamic shear failure, nanoindentation and brittle fracture, anisotropic in-plane and out-of-plane thermal conductivities, including the role of defects. Agreement with available experimental data is consistently close and exceeds DFT accuracy across the benchmark properties examined. The IFF/IFF-R model is compatible with CHARMM, AMBER, OPLS, and CVFF force fields for simulations of MXenes with diverse surface terminations, electrolyte interfaces, biointerfaces, and polymer composites without additional parameters. Parameter sets, 3D models, and analysis scripts are provided for community use. The validated, reactive, and transferable IFF framework facilitates predictive design of MXene-based films, membranes, sensing interfaces, and composites.

MXene↗
DOE OSTI2025

Model form and sensitivity analysis of CALPHAD-based nucleation models in b-stabilized Ti alloys

Accurate prediction of α-phase nucleation and growth in β-stabilized titanium alloys is crucial for designing heat treatments to optimize mechanical properties in additively manufactured lightweight components. Ideally, predictions of nucleation and growth would incorporate both top-down observations of past experimental heat treatments and bottom-up modeling of phase transformations; however, the appropriate method of combining these information sources is not self-evident. Combining top-down and bottom-up information requires a unified form of model that can connect between spatiotemporal scales, as well as sets of fitting parameters that can be identified by each data source. The selection of which parameters to fit to which data source can be made based on expert opinion, or by performing a sensitivity analysis. In solid-solid nucleation, direct observation of the nucleation and growth process is challenging. Most data on the heat treatment-controlled phase transformations are not in-situ. To predict the process and outcome of the nucleation, growth and coarsening of precipitates, theoretical models of the nucleation pathway are used to bridge the gap. Many sources of uncertainty affect the modeling of this nucleation process. It can be influenced by small variations in the thermomechanical processing history, chemical composition, and initial microstructure. If molecular dynamics (MD) simulations are used to determine thermodynamic quantities and inform CALPHAD modeling, additional uncertainty can be introduced and accounted for using Bayesian methods. Top-down uncertainties require additional steps to quantify. The influence of nucleation model form on the sensitivity of predictions to input parameters and physical conditions is the focus of this study. Classical nucleation theory (CNT) allows modeling to formulate the nucleation as homogeneous or, more commonly, heterogeneous. Non-classical nucleation models are also increasingly explored as a means of reconciling top-down and bottom-up data. In this study, the sensitivity of the intragranular nucleation of α in a β-annealed, slow-cooled aging (BASCA) heat treatment of β-stabilized Ti5553 alloy is explored using CNT and both heterogeneous and homogeneous assumptions. The Kampmann-Wagner Numerical model of precipitate nucleation and growth is employed. Using open-source tools (pyCalphad and thermodynamic modeling of TiMo as a surrogate system, a sensitivity analysis is performed to measure variations in key parameters, including chemical driving force, interfacial energy, and diffusivity, as they relate to predictions of precipitate number density. The inclusion of top-down and bottom-up data in selection of nucleation model form is discussed.

Rodriguez Negron, A. M.↗
DOE OSTI2026

Deriving effective electrode–ion interactions from free-energy profiles at electrochemical interfaces

Understanding ion adsorption at electrified metal–electrolyte interfaces is essential for accurate modeling of electrochemical systems. Here, in this study, we systematically investigate the free energy profiles of Na + , Cl − , and F − ions at the Au(111)–water interface using enhanced sampling molecular dynamics with both classical force fields and machine-learned interatomic potentials (MLIPs). Our classical metadynamics results reveal a strong dependence of predicted ion adsorption on the Lennard-Jones parameters, highlighting that—without due care—standard mixing rules can lead to qualitatively incorrect descriptions of ion–metal interactions. We present a systematic methodology for tuning the cross term LJ parameters to control adsorption energetics in agreement with more accurate models. As a surrogate for an ab initio model, we employed the recently released Universal Models for Atoms MLIP, which validates classical trends and displays strong specific adsorption for chloride, weak adsorption for fluoride, and no specific adsorption for sodium, in agreement with experimental and theoretical expectations. By integrating molecular-level adsorption free energies into continuum models of the electric double layer, we show that specific ion adsorption substantially alters the interfacial ion population, the potential of zero charge, and the differential capacitance of the system. Our results underscore the critical importance of force field parameterization and advanced interatomic potentials for the predictive modeling of ion-specific effects at electrified interfaces and provide a robust framework for bridging molecular simulations and continuum electrochemical models.

Roncoroni, Fabrice [Lawrence Berkeley National Lab↗
DOE OSTI2022

Relation between interfacial shear and friction force in 2D materials

Understanding the interfacial properties between an atomic layer and its substrate is of key interest at both the fundamental and technological levels. From Fermi level pinning to strain engineering and superlubricity, the interaction between a single atomic layer and its substrate governs electronic, mechanical and chemical properties. Here, we measure the hardly accessible interfacial transverse shear modulus of an atomic layer on a substrate. By performing measurements on bulk graphite, and on epitaxial graphene films on SiC with different stacking orders and twisting, as well as in the presence of intercalated hydrogen, we find that the interfacial transverse shear modulus is critically controlled by the stacking order and the atomic layer–substrate interaction. Importantly, we demonstrate that this modulus is a pivotal measurable property to control and predict sliding friction in supported two-dimensional materials. The experiments demonstrate a reciprocal relationship between friction force per unit contact area and interfacial shear modulus. As a result, the same relationship emerges from simulations with simple friction models, where the atomic layer– substrate interaction controls the shear stiffness and therefore the resulting friction dissipation.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Charge Separation at Organic Interfaces with Near-Zero Energy Offset: A Step toward Designing Inorganic-like Organic Semiconductors

For organic semiconductors, it is often presumed that an energy level offset at the donor–acceptor (D-A) interface is required to provide the driving force for charge separation (CS). This energy level offset unavoidably leads to a voltage loss in organic photovoltaics. In this work, by using zinc phthalocyanine (ZnPc) and fluorinated zinc phthalocyanine (F 4 ZnPc) as a model D/A interface, we found that spontaneous CS, with an enthalpy increase of ∼0.3–0.4 eV, can occur even with an interfacial energy offset as small as ∼0.1 eV. This enthalpy-increase CS process is driven by entropy. The entropic driving force can be enhanced by two factors: (1) a point-like spatial contact between the delocalized electron and hole wave function in the charge transfer exciton; (2) a small band bending near the interface originating from long-range electrostatic interaction. Our work demonstrates that effective CS can occur at interfaces with a near-zero energy level offset, which means that the energy loss at the D/A interface can be avoided.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Tunable Topological Phonon for Next-generation Quantum Transduction

This project aims to understand the critical factors that determine transduction performance of topological phonons across an oxide perovskite/tungsten diselenide heterojunction. We will investigate mechanisms of their propagation and interfacial coupling using modeling and machine learning approaches. Methods include density functional theory, molecular dynamics, numerical transport simulations, and active learning for building up a training dataset for force field development.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Numerical simulation of evaporating wavy falling liquid films in laminar gas streams

Numerical simulations are performed to investigate the interfacial heat and mass transfer from evaporating wavy falling liquid films in interaction with laminar gas streams. The OpenFOAM solver has been used to conduct the simulations where the liquid-gas interface is resolved using the Volume of Fluid method of solving for three phases, i.e., liquid, vapor, and air. The configuration considered is a falling liquid (water) film on a heated vertical plate with a confined laminar moist air (gas) flow that is either (a) co-current or (b) counter-current to the downward liquid flow. The evaporation at the liquid-gas interface is driven by the interfacial gradient of the vapor mass fraction. Interfacial waves are triggered using a monochromatic forcing disturbance that leads to sinusoidal or solitary waves forming at the liquid-gas interface under respective forcing frequencies. Further, the numerical model is validated with the available experimental data. The results show nearly a 15% enhancement in time-averaged Sherwood number (Sh) due to film waviness (sinusoidal or solitary) at the lower volumetric gas flow rate, Q g = +50 (co-current) and Q g = -50 (counter-current). This enhancement in the Sh for both the waves further increases by 11% with Q g = +800 and 196% with Q g = -800. A closer examination of the mass transfer process over a wave demonstrates that with Q g = +50, the concentration of the gas side streamlines at the trough locations of the wave leads to higher values of Sh at these locations. However, with Q g = +800, although the overall Sh increases, vortices appear at the wave trough locations, leading to a corresponding decrease in the local Sh values. Correlations are proposed for predicting Sh under co-current and counter-current gas flow effects.

36 MATERIALS SCIENCE↗
DOE OSTI2025

LiNi0.8Mn0.1Co0.1O2 Thin Films Prepared by Polymer-Assisted Deposition for the Study of Cathode-Electrolyte Interphases in Lithium-Ion Batteries

High-nickel layered oxide cathodes such as LiNi0.8Mn0.1Co0.1O2 (NMC811) are critical for next-generation lithium-ion batteries (LIBs) due to their superior energy density and reduced reliance on cobalt. However, many Ni-rich cathodes suffer from rapid capacity fade and structural instability originating from complex interfacial reactions at the cathode-electrolyte interface. Traditional composite electrodes exhibit degradation mechanisms that are challenging to quantitatively understand due to additives, including binders and carbon black. In this study, we demonstrate a new synthesis approach for binder- and additive-free NMC811 thin films using polymer-assisted deposition (PAD). PAD-NMC811 are model thin-film cathodes for investigating interfacial phenomena that can be obscured in composite cathodes. Structural and chemical characterization by X-ray diffraction, soft X-ray absorption spectroscopy, and atomic force microscopy show that PAD-NMC811 films possess high phase purity, crystallinity, chemical homogeneity, and morphological uniformity. Electrochemical analyses using cyclic voltammetry and galvanostatic cycling revealed electrochemical behavior consistent with that of composite electrodes, along with a moderate capacity fade indicative of cathode-electrolyte interphase (CEI) formation. Our findings illustrate the effectiveness of PAD synthesis of thin films tailored for detailed mechanistic studies, which offer critical insights into CEI evolution and cathode degradation pathways.

25 ENERGY STORAGE↗
DOE OSTI2024

Role of Dynamic Polarization Interactions in the Electrical Double Layer at Calcite (104) Interfaces with Aqueous Solutions

Reactions at mineral interfaces with aqueous solutions control many geochemical and biogeochemical processes in the Earth’s critical zone. At the molecular level, insights into important properties such as the structure of the electrical double layer (EDL) at specific mineral interfaces continue to improve due to the increasing fidelity of laboratory instrumentation and computational approaches. However, molecular simulation approaches suffer from limited reach into relevant scales of time, length, and system complexity. To span this gap, a novel hybrid approach that couples first principle plane-wave density functional theory (DFT) with classical DFT (cDFT) is demonstrated and applied to calcite (104) interfaces with various electrolytes. In this approach, a region of interest described using DFT interacts with the surrounding medium described using cDFT to arrive at a self-consistent ground state. Benchmarking against experimental observations and entirely first principle DFT simulations demonstrates that this hybrid model efficiently encompasses the key short-range and collective interactions in the EDL. Simulations of calcite (104)/solution interfaces reveal the key static and dynamic polarization interactions that give rise to structuring of ions and water. Ion hydration interactions have the strongest effect on the depth of the first minimum in the density distribution of counterions at the surface, and the position and width of the first density peak is largely determined by the strength of ion-correlation forces. Finer details of ion distributions are controlled by mutual polarization of the calcite surface and interfacial electrolyte. Finally, this new ability to efficiently and rigorously predict EDL structure at mineral surfaces in contact with complex solutions paves the way to accurately modeling sorption, nucleation, dissolution, and growth in realistic systems.

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

Differences in the Interfacial Mechanical Properties of Thiophosphate and Argyrodite Solid Electrolytes and Their Composites

Interfacial mechanics are a significant contributor to the performance and degradation of solid-state batteries. Spatially resolved measurements of interfacial properties are extremely important to effectively model and understand the electrochemical behavior. Herein, we report the interfacial properties of thiophosphate (Li 3 PS 4 )- and argyrodite (Li 6 PS 5 Cl)-type solid electrolytes. Using atomic force microscopy, we showcase the differences in the surface morphology as well as adhesion of these materials. Additionally we investigate solvent-less processing of hybrid electrolytes using UV-assisted curing. Physical, chemical, and structural characterizations of the materials highlight the differences in the surface morphology, chemical makeup, and distribution of the inorganic phases between the argyrodite and thiophosphate solid electrolytes.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Homogeneous ice nucleation in an ab initio machine-learning model of water

Molecular simulations have provided valuable insight into the microscopic mechanisms underlying homogeneous ice nucleation. While empirical models have been used extensively to study this phenomenon, simulations based on first-principles calculations have so far proven prohibitively expensive. Here, we circumvent this difficulty by using an efficient machine-learning model trained on density-functional theory energies and forces. We compute nucleation rates at atmospheric pressure, over a broad range of supercoolings, using the seeding technique and systems of up to hundreds of thousands of atoms simulated with ab initio accuracy. The key quantity provided by the seeding technique is the size of the critical cluster (i.e., a size such that the cluster has equal probabilities of growing or melting at the given supersaturation), which is used together with the equations of classical nucleation theory to compute nucleation rates. We find that nucleation rates for our model at moderate supercoolings are in good agreement with experimental measurements within the error of our calculation. We also study the impact of properties such as the thermodynamic driving force, interfacial free energy, and stacking disorder on the calculated rates.

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

Driving forces for particle-based crystallization: From experiments to theory and simulations

The multi-step crystallization processes involving the formation of stable building blocks that subsequently assemble into a crystal are ubiquitous in mineral formation and biomineralization and are particularly attractive in materials synthesis. Utilizing these pathways offers the approach to overcoming the restrictions on the expression of various crystal faces imposed by the interfacial energy during monomer-by-monomer growth to unlock the breadth of architectures with unique properties. Further, controlling particle-based crystallization proved challenging despite its promise due to the complex interdependence of interfacial forces and their non-linear dependence on synthesis parameters. Here, the status of the development of state-of-the-art approaches to measuring interparticle forces and predictive theoretical models of particle-based crystallization are reviewed.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Crystallization pathways and interfacial drivers for the formation of hierarchical architectures

The development of structural hierarchy on various length scales during crystallization process is ubiquitous in biological systems and is also observed in synthetic nanomaterials. The driving forces for the formations of complex architectures range from local interfacial interactions, that modify interfacial speciation, local supersaturation, and nucleation barriers, to macroscopic interparticle forces. Although it is enticing to interpret the formation of hierarchical architectures as the assembly of independently nucleated building blocks, often crystallization pathways follow monomer-by-monomer addition with structural complexity arising from interfacial chemical coupling and strongly correlated fluctuation dynamics in the electric double layers. Here, the mechanism of the development of structural hierarchy through heterogeneous nucleation, coupled interfacial nucleation and assembly, and oriented attachment of independently nucleated particles is discussed. In this article, the emphasis is made on the discussion of the underlying interfacial forces and chemical coupling that drives crystallization pathways towards the formation of structural hierarchy.

36 MATERIALS SCIENCE↗