Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Interfacial flows”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

An SECM-Based Spot Analysis for Redoxmer-Electrode Kinetics: Identifying Redox Asymmetries on Model Graphitic Carbon Interfaces

The fundamental process in non-aqueous redox flow battery (NRFB) operation revolves around electron transfer (ET) between a current collector electrode and redox-active organic molecules (redoxmers) in solution. Here, we present an approach utilizing scanning electrochemical microscopy (SECM) to evaluate interfacial ET kinetics between redoxmers and various electrode materials of interest at desired locations. This spot-analysis method relies on the measurement of heterogeneous electron transfer rate constants (k f or k b ) as a function of applied potential (E-E 0 '). As demonstrated by COMSOL simulations, this method enables the quantification of Butler-Volmer kinetic parameters, the standard heterogeneous rate constant, k 0 , and the transfer coefficient, α. Our method enabled the identification of inherent asymmetries in the ET kinetics arising during the reduction of ferrocene-based redoxmers, compared to their oxidation which displayed faster rate constants. Similar behavior was observed on a wide variety of carbon electrodes such as multi-layer graphene, highly ordered pyrolytic graphite, glassy carbon, and chemical vapor deposition-grown graphite films. However, aqueous systems and Pt do not exhibit such kinetic effects. Our analysis suggests that differential adsorption of the redoxmers is insufficient to account for our observations. Displaying a greater versatility than conventional electroanalytical methods, we demonstrate the operation of our spot analysis at concentrations up to 100 mM of redoxmer over graphite films. Looking forward, our method can be used to assess non-idealities in a variety of redoxmer/electrode/solvent systems with quantitative evaluation of kinetics for applications in redox-flow battery research.

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

Nanocrystal-to-Ligand Interfacial Thermal Transport: Plasmonic Indium Tin Oxide to Perylene Diimide Adsorbates

The generation and dissipation of heat in nanocrystals upon optical excitation is a process which may either limit or enhance their performance in certain applications. Using transient absorption spectroscopy, we tracked the flow of heat from plasmonic tin-doped indium oxide (ITO) nanocrystals to surface adsorbed perylenediimide (PDI) molecules upon excitation of the ITO near-infrared plasmon resonance. Here, we rationalize the derivative line shapes observed in the transient absorption bleach features are the result of thermal transfer from the ITO core to surface ligands based on temperature-dependent static absorption studies of the PDI molecules adsorbed on plasmonic ITO. Through a series of pump power dependent measurements, we demonstrate that the PDI heating time is largely fluence independent at the powers measured, while the overall recovery time increases. Elucidating thermal transfer rates at the nanocrystal and organic ligand interface is key for applications such as plasmon mediated photocatalysis, where heating may obfuscate hot carrier transfer processes.

heat dissipation↗
DOE OSTI2024

Dual material fused filament fabrication of composite core‐shell structures with improved impact resistance and interfacial adhesion

The mechanical performance of parts produced by fused filament fabrication (FFF) has been limited due to the presence of voids and poor interlayer welding. Recent advancements in FFF have enabled the fabrication of void-free objects with strong interlayer welding through the use of semicrystalline polymer shells such as high-density polyethylene (HDPE) along with a high viscosity core polymer like acrylonitrile-butadiene-styrene (ABS). Furthermore, the zero-shear viscosity (η0) of ABS is three orders of magnitude higher than HDPE making the ABS-HDPE core-shell configuration preferable. ABS holds the shape by preventing bulk flow and part bending while HDPE promotes full surface contact across the layers. Most polymers, however, are immiscible which causes a weak weld line along the core-shell interface. Herein, maleic anhydride (MAH) is grafted to the butadiene segment of the ABS core thereby compatibilizing the interface with HDPE, improving the interfacial adhesion. Attenuated total reflectance-Fourier transform infrared spectroscopy was employed to confirm successful grafting. Using a custom-made die affording the core-shell structure, the ABS-g-MAH is shown to improve the impact resistance by 253% and 16% compared to neat HDPE and ABS specimens, respectively. Additionally, a 10% increase compared to the unmodified ABS-HDPE core-shell configuration is observed.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Seamlessly joining length scales: From atomistic thermal graphs to anisotropic continuum conductivity

Thermal transport in complex solids is governed by local structure, defects, and anisotropy, yet most continuum models still rely on oversimplified and homogenized conductivities. Here, we bridge atomistic and continuum descriptions by building finite element (FE) models directly from the site-projected thermal conductivity (SPTC), an atomic-level decomposition of the Green–Kubo thermal conductivity. We introduce a toolkit, the “Simulator Collection for Atomic-to-Continuum Scales (SCACS)”, which uses a graph neural network to predict SPTC on large atomic structures, coarse-grains these fields into anisotropic conductivity tensors, and embeds them into the heat-flow FE equation with a customized, anisotropy-aware adaptive mesh refinement scheme. Applied to silicon nanostructures, the resulting FE models act as representative volume elements, reproduce bulk conductivities, and capture interfacial and defect-driven anisotropy while maintaining thermodynamic consistency. Additionally, SCACS predicts experimental conductance trends and fields. This work demonstrates a general route for transferring atomistic transport information into device-scale thermal simulations with physics-based approximations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

Simulating Bubble Deformations Using the Intersection Marker (ISM) Front-Tracking method

Surface and volume tracking methods have been under continuous development over the last 50+ years since the first MAC (Marker and Cell) and SLIC (Simple Line Interface Calculation) methods appeared in the 1960's. The sustained interest in simulating multiphase flow stem from their many practical applications in the physical world; from biological systems, to chemical systems, to boiling dynamics and atmospheric heat and mass transfer. The ability to capture the details of interface evolution as well as the physics of interfacial heat and mass transfer has thus been a tempting goal for many numerical researchers. Numerical analysts, apart from the well-established and popular VOF PLIC (Volume of Fluid, Piecewise Linear Interface Construction) method, can now call upon even more computationally intensive options such as the LS (Level Set) method, the FT (Front Tracking) method, the MoF (Moment of Fluid) method and hybrid methods.

42 ENGINEERING↗
DOE OSTI2014

Self-sustained turbulence without dynamical forcing: A two-dimensional study of a bistable interstellar medium

In this paper, the nonlinear evolution of a bistable interstellar medium is investigated using two-dimensional simulations with a realistic cooling rate, thermal conduction, and physical viscosity. The calculations are performed using periodic boundary conditions without any external dynamical forcing. As the initial condition, a spatially uniform unstable gas under thermal equilibrium is considered. At the initial stage, the unstable gas quickly segregates into two phases: cold neutral medium (CNM) and warm neutral medium (WNM). Then, self-sustained turbulence with velocity dispersion of 0.1-0.2 km s{sup –1} is observed in which the CNM moves around in the WNM. We find that the interfacial medium (IFM) between the CNM and WNM plays an important role in sustaining the turbulence. The self-sustaining mechanism can be divided into two steps. First, thermal conduction drives fast flows streaming into concave CNM surfaces toward the WNM. The kinetic energy of the fast flows in the IFM is incorporated into that of the CNM through the phase transition. Second, turbulence inside the CNM deforms interfaces and forms other concave CNM surfaces, leading to fast flows in the IFM. This drives the first step again and a cycle is established by which turbulent motions are self-sustained.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

On the consistency of scale among experiments, theory, and simulation

As a tool for addressing problems of scale, we consider an evolving approach known as the thermodynamically constrained averaging theory (TCAT), which has broad applicability to hydrology. We consider the case of modeling of two-fluid-phase flow in porous media, and we focus on issues of scale as they relate to various measures of pressure, capillary pressure, and state equations needed to produce solvable models. We apply TCAT to perform physics-based data assimilation to understand how the internal behavior influences the macroscale state of two-fluid porous medium systems. A microfluidic experimental method and a lattice Boltzmann simulation method are used to examine a key deficiency associated with standard approaches. In a hydrologic process such as evaporation, the water content will ultimately be reduced below the irreducible wetting-phase saturation determined from experiments. This is problematic since the derived closure relationships cannot predict the associated capillary pressures for these states. We demonstrate that the irreducible wetting-phase saturation is an artifact of the experimental design, caused by the fact that the boundary pressure difference does not approximate the true capillary pressure. Using averaging methods, we compute the true capillary pressure for fluid configurations at and below the irreducible wetting-phase saturation. Results of our analysis include a state function for the capillary pressure expressed as a function of fluid saturation and interfacial area.

13 HYDRO ENERGY↗
DOE OSTI2014

Mechanical resilience and cementitious processes in Imperial Roman architectural mortar

The pyroclastic aggregate concrete of Trajan’s Markets (110 CE), now Museo Fori Imperiali in Rome, has absorbed energy from seismic ground shaking and long-term foundation settlement for nearly two millenia while remaining largely intact at the structural scale. The scientific basis of this exceptional service record is explored through computed tomography of fracture surfaces and synchroton X-ray microdiffraction analyses of a reproduction of the standardized hydrated lime–volcanic ash mortar that binds decimeter-sized tuff and brick aggregate in the conglomeratic concrete. The mortar reproduction gains fracture toughness over 180 d through progressive coalescence of calcium–aluminum-silicate–hydrate (C-A-S-H) cementing binder with Ca/(Si+Al) ≈ 0.8–0.9 and crystallization of strätlingite and siliceous hydrogarnet (katoite) at ≥90 d, after pozzolanic consumption of hydrated lime was complete. Platey strätlingite crystals toughen interfacial zones along scoria perimeters and impede macroscale propagation of crack segments. In the 1,900 year old mortar, C-A-S-H has low Ca/(Si+Al) ≈ 0.45–0.75. Dense clusters of 2- to 30-µm strätlingite plates further reinforce interfacial zones, the weakest link of modern cement-based concrete, and the cementitious matrix. These crystals formed during long-term autogeneous reaction of dissolved calcite from lime and the alkali-rich scoriae groundmass, clay mineral (halloysite), and zeolite (phillipsite and chabazite) surface textures from the Pozzolane Rosse pyroclastic flow, erupted from the nearby Alban Hills volcano. The clast-supported conglomeratic fabric of the concrete presents further resistance to fracture propagation at the structural scale.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Pore-scale modelling of multiphase reactive flow: application to mineral dissolution with production of CO 2

A micro-continuum approach is proposed to simulate the dissolution of solid minerals at the pore scale in the presence of multiple fluid phases. The approach employs an extended Darcy–Brinkman–Stokes formulation that accounts for the interfacial tension between the two immiscible fluid phases and the moving contact line at the mineral surface. The simulation framework is validated using an experimental microfluidic device that provides time-lapse images of the dissolution dynamics. The set-up involves a single-calcite crystal and the subsequent generation of CO 2 bubbles in the domain. The dissolution of the calcite crystal and the production of gas during the acidizing process are analysed. We then show that the production of CO 2 bubbles during the injection of acid in a carbonate formation may limit the overall dissolution rate and prevent the emergence of wormholes.

bio-inspired, mechanical behavior, carbon sequestr↗
DOE OSTI2018

Geometric state function for two-fluid flow in porous media

Models that describe two-fluid flow in porous media suffer from a widely recognized problem that the constitutive relationships used to predict capillary pressure as a function of the fluid saturation are nonunique, thus requiring a hysteretic description. As an alternative to the traditional perspective, we consider a geometric description of the capillary pressure, which relates the average mean curvature, the fluid saturation, the interfacial area between fluids, and the Euler characteristic. The state equation is formulated using notions from algebraic topology and cast in terms of measures of the macroscale state. Synchrotron-based x-ray microcomputed tomography and high-resolution pore-scale simulation is applied to examine the uniqueness of the proposed relationship for six different porous media. We show that the geometric state function is able to characterize the microscopic fluid configurations that result from a wide range of simulated flow conditions in an averaged sense. The geometric state function can serve as a closure relationship within macroscale models to effectively remove hysteretic behavior attributed to the arrangement of fluids within a porous medium. As a result, this provides a critical missing component needed to enable a new generation of higher fidelity models to describe two-fluid flow in porous media.

42 ENGINEERING↗
DOE OSTI2016

Tracking interface and common curve dynamics for two-fluid flow in porous media

Pore-scale studies of multiphase flow in porous medium systems can be used to understand transport mechanisms and quantitatively determine closure relations that better incorporate microscale physics into macroscale models. Multiphase flow simulators constructed using the lattice Boltzmann method provide a means to conduct such studies, including both the equilibrium and dynamic aspects. Moving, storing, and analyzing the large state space presents a computational challenge when highly-resolved models are applied. We present an approach to simulate multiphase flow processes in which in-situ analysis is applied to track multiphase flow dynamics at high temporal resolution. We compute a comprehensive set of measures of the phase distributions and the system dynamics, which can be used to aid fundamental understanding and inform closure relations for macroscale models. The measures computed include microscale point representations and macroscale averages of fluid saturations, the pressure and velocity of the fluid phases, interfacial areas, interfacial curvatures, interface and common curve velocities, interfacial orientation tensors, phase velocities and the contact angle between the fluid-fluid interface and the solid surface. Test cases are studied to validate the approach and illustrate how measures of system state can be obtained and used to inform macroscopic theory.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1988

Formula Gives Better Contact-Resistance Values

Lateral currents in contact strips taken into account. Four-terminal test structures added to intergrated circuits to enable measurement of interfacial resistivities of contacts between thin conducting layers. Thin-film model simplified quasi-two-dimensional potential model that accounts adequately for complicated three-dimensional, nonuniform current densitites. Effects of nonuniformity caused by lateral current flow in strips summarized in equivalent resistance Rs and voltage Vs.

Lieneweg, Udo↗
NASA NTRS1994

Convective flow effects on protein crystal growth

A high-resolution microscopic interferometric setup for the monitoring of protein morphologies has been developed. Growth or dissolution of a crystal can be resolved with a long-term depth resolution of 200 A and a lateral resolution of 2 microns. This capability of simultaneously monitoring the interfacial displacement with high local depth resolution has yielded several novel results. We have found with lysozyme that (1) the normal growth rate is oscillatory, and (2) depending on the impurity content of the solution, the growth step density is either greater or lower at the periphery of a facet than in its center. The repartitioning of Na plus and Cl minus ions between lysozyme solutions and crystals was studied for a wide range of crystallization conditions. A nucleation-growth-repartitioning model was developed, to interpret the large body of data in unified way. The results strongly suggest that (1) the ion to lysozyne ratio in the crystal depends mostly on kinetic rather than crystallographic parameters, and (2) lysozyme crystals possess a salt-rich core with a diameter electron microscopy results appear to confirm this finding, which could have far-reaching consequences for x-ray diffraction studies. A computational model for diffusive-convective transport in protein crystallization has been applied to a realistic growth cell geometry, taking into account the findings of the above repartitioning studies and our kinetics data for the growth of lysozyme. The results show that even in the small cell employed, protein concentration nonuniformities and gravity-driven solutal convection can be significant. The calculated convection velocities are of the same order to magnitude as those found in earlier experiments. As expected, convective transport, i.e., at Og, lysozyme crystal growth remains kinetically limited. The salt distribution in the crystal is predicted to be non-uniform at both 1g and 0g, as a consequence of protein depletion in the solution. Static and dynamic light scattering studies in undersaturated and supersaturated solutions have been performed. Diffusivities in undersaturated solutions, were found to vary with lysozyme concentrations. Depending on the salt concentration, the diffusivities either increase or decrease. Interestingly, the corresponding static scattering intensities behave oppositely, Our current analysis indicates that these changes are inconsistent with aggregation in undersaturated solutions. However, the data are compatible with concentration-dependent changes of the interactions between protein and salt.

Rosenberger, Franz↗
DOE OSTI2018

Simultaneous “Clean-and-Repair” of Surfaces Using Smart Droplets

An autonomous self-healing system, inspired by transportation processes inherent to biology, is described for materials transportation and repair. The selected model system combines inorganic nanoparticles (NPs) on damaged substrates with functional emulsion droplets that pick up the particles from pristine portions of the substrate and deposit them into damaged regions. The droplets are stabilized by polymer surfactants containing phosphorylcholine groups, a polymer composition selected to impart surfactant properties for droplet stabilization as well as fouling resistance to prevent irreversible droplet adsorption on the substrates. Both the NP pickup (cleaning) and drop off (repair) steps are conducted in a system driven by an imposed flow and characterized by fluorescence microscopy. To evaluate and optimize the efficiency of this NP transportation process, the effect of both the chemical composition of the polymer surfactant and the NP surface chemistry is investigated in this work. Interfacial interactions proved enabling for these NP transportation processes, specifically those involving NP/droplet, NP/substrate, and droplet/substrate interactions. Ultimately, droplets capable of both picking up and dropping off NPs are realized by adjusting fluid/fluid and fluid/substrate interactions, with electrostatic interactions between NPs and droplets proving most effective.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Stability-based optimization of ACRT for the growth of CZT by the traveling heater method

Numerical modeling is employed to study the effects of several accelerated crucible rotation technique (ACRT) rotation schedules on the growth of cadmium zinc telluride (CZT) by the traveling heater method (THM). In conjunction with these analyses, a measure is developed to estimate the amplitude of disturbances associated with cellular interface growth morphologies that develop over time and across the melt-solid interface. Furthermore, this disturbance amplitude measure is used to assess the ability of ACRT to stabilize the growth interface and thereby minimize inclusion formation processes. Notably, our analyses reveal that interfacial instability is not well correlated with melt undercooling, explaining why classical ACRT approaches to reduce undercooling by mixing may not stabilize growth. Computations show that interface stability is best achieved via rotation schedules that accentuate outward Ekman flows during spin-up and, during spin-down, eliminate Taylor-Görtler flows and minimize inward Ekman flows.

36 MATERIALS SCIENCE↗
NASA NTRS1984

Diffusive-convective vapor transport across horizontal and inclined rectangular enclosures

The enhancement of vapor transport across horizontal and inclined rectangular (two-dimensional) enclosures by expansive convection is studied numerically for a range in Grashof number that extends up to boundary layer flow. Comparison is made with predictions of the simplifying model and resulting analytical treatment of Klosse and Ullersma (1973) for various aspect ratios of the enclosure. The influence of varying Schmidt and Prandtl numbers is investigated. Solutions for interfacial (growth) flux and temperature distributions are given, showing practically important nonuniformities.

Markham, B. L.↗
DOE OSTI2024

Impact of Pressure-Dependent Interfacial Tension and Contact Angle on Capillary Heterogeneity Trapping of CO2 in Storage Aquifers

Summary Carbon dioxide (CO2) capillary trapping increases the total amount of CO2 that can be effectively immobilized in storage aquifers. This trapping, manifesting itself as accumulated CO2 columns at a continuum scale, is because of capillary threshold effects that occur below low-permeability barriers. Considering that capillary pressure is dictated by heterogeneous pore throat size, the trapped CO2 column height and associated CO2 saturation will vary spatially within a storage aquifer. This variation will be influenced by two pressure-dependent interfacial parameters—CO2/brine interfacial tension (IFT) and CO2/brine/rock contact angle. Our objective is to understand how the pressure dependence of these two parameters affects the heterogeneity of capillary trapped CO2 at a continuum scale. Our conceptual model is a 1D two-zone system with the upper zone being a flow barrier (low permeability) and the lower zone being a flow path (high permeability). The inputs to this model include microfacies-dependent capillary pressure vs. saturation curves and permeability values. The input capillary pressure curves were collected in the literature that represents carbonate microfacies (e.g., dolograinstone) in a prevalent formation in the Permian Basin. We then used the Leverett j-function to scale the capillary pressure curve for the two zones that are assigned with the same or different microfacies. During scaling, we considered the influence of pressure on both the IFT and contact angle of CO2/brine/dolomite systems. We varied the zone permeability contrast ratio from 2 to 50. We then assumed capillary gravity equilibriums and calculated the CO2 saturation buildup corresponding to various trapped CO2 column heights. The CO2 saturation buildup is defined as the CO2 saturation in the lower layer minus that in the upper one. We found that the saturation buildup can be doubled when varying pressure in a storage aquifer, after considering pressure-dependent IFT and contact angles. Thus, assuming these two parameters to be constant across such aquifers would cause large errors in the quantification of capillary trapping of CO2. The whole study demonstrates the importance of considering pressure-dependent interfacial properties in predicting the vertical distribution of capillary trapped CO2. It has important implications in developing a better understanding of leakage risks and consequent storage safety.

Engineering↗