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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Exchange bias magnetism in films of NiFe/(Ni,Fe)O nanocrystallite dispersions

Ni{sub 3}Fe/(Ni,Fe)O thin films having a nanocrystallite dispersion morphology were prepared by a reactive ion beam-assisted deposition technique. The crystallite sizes of these dispersion-based films were observed to decrease from 8.4 ± 0.3 nm to 3.4 ± 0.3 nm as the deposition flow-rate increased from 2.78% to 7.89% O{sub 2}/Ar. Thin film composition was determined using selective area electron diffraction images and Multislice simulations. Through a detailed analysis of high resolution transmission electron microscopy images, the nanocrystallites were determined to be Ni{sub 3}Fe (a ferromagnet), NiO, and FeO (both antiferromagnets). It was determined that the interfacial molar Ni{sub 3}Fe ratio in the nanocrystallite dispersions increased slightly at first, then decreased as the oxygen content was increased; at 7.89% O{sub 2}/Ar, the interfacial molar ratio was essentially zero (only NiO and FeO remained). For nanocrystallite dispersion films grown with O{sub 2}/Ar flow-rate greater than 7.89%, no interfacial (intermixed) Ni{sub 3}Fe phase was detected, which resulted in no measurable exchange bias. Comparing the exchange bias field between the nanocrystallite dispersion films at 5 K, we observed a decrease in the magnitude of the exchange bias field as the nanocrystallite size decreased. The exchange bias coupling for all samples measured set in at essentially the same temperature (i.e., the exchange bias blocking temperature). Since the ferromagnetic/anti-ferromagnetic (FM/AFM) contact area in the nanocrystallite dispersion films increased as the nanocrystallite size decreased, the increase in the magnitude of the exchange bias could be attributed to larger regions of defects (vacancies and bond distortions) which occupied a significant portion of the FM/AFM interfaces in the nanocrystallite dispersion films.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Multiple-relaxation-time color-gradient lattice Boltzmann model for simulating two-phase flows with high density ratio

In this paper, we propose a color-gradient lattice Boltzmann (LB) model for simulating two-phase flows with high density ratio and high Reynolds number. The model applies a multi-relaxation-time (MRT) collision operator to enhance the stability of the simulation. A source term, which is derived by the Chapman-Enskog analysis, is added into the MRT LB equation so that the Navier-Stokes equations can be exactly recovered. Also, a new form of the equilibrium density distribution function is used to simplify the source term. To validate the proposed model, steady flows of a static droplet and the layered channel flow are first simulated with density ratios up to 1000. Small values of spurious velocities and interfacial tension errors are found in the static droplet test, and improved profiles of velocity are obtained by the present model in simulating channel flows. Then, two cases of unsteady flows, Rayleigh-Taylor instability and droplet splashing on a thin film, are simulated. In the former case, the density ratio of 3 and Reynolds numbers of 256 and 2048 are considered. The interface shapes and spike/bubble positions are in good agreement with the results of previous studies. In the latter case, the droplet spreading radius is found to obey the power law proposed in previous studies for the density ratio of 100 and Reynolds number up to 500.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Characterization of local parameters of inclined upward air–water bubbly two-phase flows in a round pipe

The present work seeks to develop a detailed experimental database and characterize the local effects of inclination on horizontal and upward bubbly air–water two-phase flows. Experimental conditions are selected to investigate bubbly flows at a variety of orientations while reducing the influence of flow regime transition. Detailed time-averaged data is collected using state-of-the-art four-sensor conductivity probes for measurements across the flow area for five flow conditions at five inclination orientations in a 25.4 mm inner diameter pipe. Using the newly established database, effects of inclination on bubbly flows are analyzed, as are changes in liquid and gas flow rate at different orientations. The observations include (1) bubble migration and agglomeration at the top pipe wall as orientation changes from vertical to horizontal; (2) reducing the angle from vertical reduces the local gas velocity, and correspondingly reduces the global slip from greater than unity at 90° to less than unity for 60° and lower; (3) changes in the bubble size distribution are primarily at the interface between the bubble cluster and the liquid region of the pipe, not coalescence within the bubble cluster. In conclusion, physical mechanisms for these observations are proposed, and the present trends at different orientations support previous experiments in horizontal two-phase flows.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Tracking Local pH Dynamics during Water Electrolysis via In-Line Continuous Flow Raman Spectroscopy

The performance of electrochemical devices, which play a critical role in decarbonization efforts, is often governed by proton-coupled electron transfer reactions at the electrode–electrolyte interface. These reactions are highly sensitive to the complex and dynamic microenvironment present at the electrode surface. However, characterizing this environment─particularly monitoring interfacial pH and its evolution under reaction conditions─remains challenging, necessitating the development of advanced analytical tools. Here, in this study, we introduce in-line continuous flow Raman spectroscopy (CFRS) as a spectroelectrochemical platform for quantifying interfacial pH swings generated during water-splitting. By monitoring phosphate ion speciation and controlling the hydrodynamics with a flow cell, we measure pH swings as a function of current density, flow rate, and distance from the electrode. Comparison with theoretical models reveals the impact of bulk pH, boundary layer thickness, and bubble dynamics at high current densities. Collectively, these findings establish CFRS as a platform for quantitatively investigating pH dynamics, offering critical insights for advancing electrochemical energy conversion technologies.

Marquez, Raul A. [Univ. of Texas, Austin, TX (Unit↗
DOE OSTI2023

Electrochemical Flow Reactor Design Allows Tunable Mass Transport Conditions for Operando Surface Enhanced Infrared Absorption Spectroscopy

Abstract In situ attenuated total reflection surface enhanced infrared absorption spectroscopy (ATR‐SEIRAS) is often used to investigate the near‐surface electrocatalytic reaction environment. However, there is a gap in directly correlating the near‐surface reaction environment with electrocatalytic reaction rates. To that end, we designed an electrochemical flow reactor for operando electrochemical ATR‐SEIRAS and demonstrate its capability with the CO 2 reduction reaction (CO 2 RR). Roughened gold catalyst thin films are prepared on ATR silicon crystals as a model system to probe local species under CO 2 RR conditions in 0.1 M KHCO 3 . We measured changes in the interfacial CO 2 concentration as a function of applied potential and electrolyte flow rate in operando , allowing us to correlate the changes in reaction rates with the observed CO 2 concentration. Including the choice of the catalyst and electrolyte, coupling hydrodynamic control with ATR‐SEIRAS in this platform enables investigations of how the local microenvironment affects the activity and selectivity of electrochemical reactions.

Avilés Acosta, Jaime E.↗
DOE OSTI2018

The Impact of Pore-Scale Flow Regimes on Upscaling of Immiscible Two-Phase Flow in Porous Media

Empirical or theoretical extensions of Darcy’s law for immiscible two-phase flow have shown significant limitations in properly modeling the flow at the continuum scale. We tackle this problem by proposing a set of upscaled equations based on pore-scale flow regimes, that is, the topology of flowing phases. The incompressible Navier-Stokes equation is upscaled by means of multiple-scale expansions and its closures derived from the mechanical energy balance for different flow regimes at the pore scale. We also derive the applicability conditions of the upscaled equations based on the order of magnitude of relevant dimensionless numbers, that is, Eotvos, Reynolds, capillary, Froude numbers, and the viscosity and density ratio of the system, as well as a set of closures valid for the basic flow regimes of low Eotvos number systems, that is, core-annular and plug and drop traffic flows. We provide analytical expressions for the relative permeability of the wetting and nonwetting phases in different flow regimes and demonstrate that the effect of the flowing-phases topology on the relative permeabilities is significant. Finally, we show that the classical two-phase Darcy law is recovered for a limited range of operative conditions, while specific terms accounting for interfacial and wall interactions should be incorporated to accurately model ganglia or drop traffic flow.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Local Measurement in Annular Flows Using a Two-Sensor Droplet-Capable Conductivity Probe

Two-sensor electrical conductivity probes have been widely used for detailed local measurement in two-phase flow systems. However application of conductivity probes in liquid-dispersed flows has been limited. In the churn-turbulent to annular transition and annular flow regimes, the liquid phase exists as a continuous phase in the film and ligaments, and as a dispersed phase in the droplets. The conventional design of a conductivity probe is not valid because it cannot distinguish the dispersed liquid phase from the gas phase. In this paper, the Two-Sensor Droplet-Capable Conductivity Probe (DCCP-2) is designed in order to extend measurement capability to liquid-dispersed flows. The liquid phase is classified into droplets and ligaments. The interface velocity, volume fraction, interfacial area concentration, as well as the chord length of liquid droplets and ligaments, can be calculated from the raw signal. The radial distribution of these parameters and the probability distribution of chord length and velocity are obtained for downwards flow in a vertical pipe. The results show that although the average volume fraction of droplets is only on the order of 1%, droplets contribute to around 20%-30% to the entrained liquid phase in the gas core, the remainder as ligaments. The interfacial area concentration from droplets and ligaments are in the order of 300 m{sup -1}. (authors)

42 ENGINEERING↗
DOE OSTI2022

Optimizing conformance control for gas injection EOR in unconventional reservoirs

Field observations have shown that interference between wells caused by connected fractures commonly exists in unconventional reservoirs. The well interference effect and conformance issues could lead to minimal enhanced oil recovery response (EOR) in actual operations if they are not properly controlled. A series of data analysis, modeling, and experimental activities were performed in this study to investigate the gas EOR strategies in the Bakken Formation. Employing the embedded discrete fracture modeling technology, a multiple-well, multiple-fracture model was developed to capture the well interference and conformance issues in the Bakken. The model made it possible to investigate the complex flow behavior between wells in the EOR (huff ‘n’ puff) processes. Alternative EOR and conformance control strategies using water and surfactant injections were studied to improve oil production performance based on experimental observations. Results showed that water injection could be used for conformance control in the Bakken because of its relatively low mobility compared to gas; however, water could also block oil flow around the conformance control wells because of the high interfacial tension between oil and water. Adding surfactant to the injection water is more effective for improving gas EOR performance since surfactant reduces the oil–water interfacial tension and contact angle effectively. Up to 9.4% more oil could be produced from the huff ‘n’ puff wells in 2 years of EOR operations when rich gas was injected at only 3 MMscf/d with 6000 psi of surfactant injection in the offset wells for conformance control.

42 ENGINEERING↗
NASA NTRS1985

Modeling of convection phenomena in Bridgman-Stockbarger crystal growth

Thermal convection phenomena in a vertically oriented Bridgman-Stockbarger apparatus were modeled by computer simulations for different gravity conditions, ranging from earth conditions to extremely low gravity, approximate space conditions. The modeling results were obtained by the application of a state-of-the art, transient, multi-dimensional, completely densimetrically coupled, discrete-element computational model which was specifically developed for the simulation of flow, temperature, and species concentration conditions in two-phase (solid-liquid) systems. The computational model was applied to the simulation of the flow and the thermal conditions associated with the convection phenomena in a modified Germanium-Silicon charge enclosed in a stationary fused-silica ampoule. The results clearly indicated that the gravitational field strength influences the characteristics of the coherent vortical flow patterns, interface shape and position, maximum melt velocity, and interfacial normal temperature gradient.

Carlson, F. M.↗
NASA NTRS1991

Fluid seals technology at BHR Group Ltd.

A mathematical model was developed for ESA, which couples the rigid body dynamics in four degrees of freedom of a flexibly mounted stator (with Coulomb damping) of a turbopump face seal with the fluid film behavior. The interfacial film is assumed gaseous, and the model incorporates laminar/turbulent flow, inertia and sonic choking. Seal faces may be wavy and coned, while the rotating ring can be eccentrically mounted, misaligned with the shaft and oscillate axially at arbitrary frequency. The model is currently being coded into software. An experimental program using a high speed cryogenic/hot gas test rig is being conducted to verify the output of the computer program.

Leefe, Simon↗
DOE OSTI2019

The Importance of Viscous and Interfacial Forces in the Hydrodynamics of the Top-Submerged-Lance Furnace

The purpose of this work is to focus on the hydrodynamics of a Top-Submerged-Lance (TSL) smelting furnace, understanding how liquid properties and operational parameters act on key factors of a TSL process, such as splashing, mixing, mass transfer area, and bubble development. A deep knowledge of all those aspects is needed since they all influence the smelting reaction rates; hence the efficiency of the reactor. The characterization and scaling of the TSL gas injection are commonly based on the modified Froude number, the ratio of dynamic and gravitational forces. Detailed literature research reveals a potential weakness of this approach, since it does not consider the effects of viscosity and surface tension. To investigate this question an extensive parametric study was performed applying computational fluid dynamics to cold and non-reactive flows, which provided a broad overview of the physics of the flow. The analysis was performed on fluid dynamic properties (liquid density, liquid viscosity, surface tension) and operational variables (gas volume flow, lance immersion depth). The coupled Level Set—Volume of Fluid model, available in the commercial solver ANSYS Fluent{sup ®}, was used to resolve the gas–liquid interface in the multiphase flow. The results of the work underscore the significance of the viscous and interfacial forces for gas injection in smelting slags, confirming the incompleteness of applying only the Froude number to describe such flows.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Fluid – Fluid Interfacial Area and Its Impact on Relative Permeability - A Pore Network Modeling Study

Abstract Relative permeability (kr) is commonly modeled as an empirical function of phase saturation. Although current empirical models can provide a good match of one or two measured relative permeabilities using saturation alone, they are unable to predict relative permeabilities well when there is hysteresis or when physical properties such as wettability change. Further, current models often result in relative permeability discontinuities that can cause convergence and accuracy problems in simulation. To overcome these problems, recent research has modeled relative permeability as a state function of both saturation (S) and phase connectivity (X). Pore network modeling (PNM) data, however, shows small differences in relative permeability for the same S-X value when approached from a different flow direction. This paper examines the impact of one additional Minkowski parameter (Mecke and Arns, 2005), the fluid-fluid interfacial area, on relative permeability to identify if that satisfactorily explains this discrepancy. We calculate the total fluid-fluid interfacial areas (IA) during two-phase (oil/water) flow in porous media using pore network modeling. The area is calculated from PNM simulations using the areas associated with corners and throats in pore elements of different shapes. The pore network is modeled after a Bentheimer sandstone, using square, triangular, and circular pore shapes. Simulations were conducted for numerous primary drainage and imbibition cycles at a constant contact angle of 0° for the wetting phase. Simultaneous measurements of capillary pressure, relative permeability, saturation, and phase connectivity are made for each displacement. Fluid-fluid interfacial area is calculated from the PNM capillary pressure, the fluid location in the pore elements, and the pore element dimensional data. The results show that differences in the relative permeability at the same (S,X) point is explained well by differences in the fluid-fluid interfacial area (IA). That is, for a larger change in IA at these intersection points, the permeability difference is greater. That difference in relative permeability approaches zero as the difference in IA approaches zero. This confirms that relative permeability can be modeled better as a unique function of S, X, and IA. The results also show that an increase in IA restricts flow decreasing the nonwetting (oil) phase permeability. This decrease is caused by an increase in the throat area fraction compared to the corner area as the total area IA increases. The wetting phase relative permeability, however, shows the inverse trend, in that its relative permeability is greater when IA becomes larger owing to a greater fraction of the total area associated with the corners. The area IA, however, impacts the nonwetting phase relative permeability more than the wetting phase relative permeability. Corner flow improves the wetting phase relative permeability because the wetting phase is continuous there. Finally, a sensitivity analysis shows that relative permeability a is more sensitive to change in S than they are for IA for the case studied implying that if only two parameters are used to model relative permeability it is better to choose S and X.

Engineering↗
DOE OSTI2022

Modeling of vapor-liquid interactions in condensing ejectors

Ejectors are compact mechanical devices that utilize the expansion of a high-pressure primary fluid to entrain and compress a low-pressure secondary fluid by means of momentum transfer between the two streams of fluid. Condensing ejectors feature both momentum and heat transfer through the interaction between the vapor and liquid streams. The goal of this paper is to develop and validate a one-dimensional slug flow model to simulate the vapor–liquid interactions in two types of condensing ejectors – one with primary liquid and secondary vapor flows (Type I), and the other with primary vapor and secondary liquid flows (Type II). Control volume analysis of the mass, momentum, and energy balance in each phase and across the interface was conducted for the ejectors. The slug flow models for both types of ejectors are validated with published results in the literature. The friction coefficient on the inner wall of the ejectors and the interfacial heat transfer coefficient are identified as controlling parameters for the ejector performance in terms of pressure and temperature distributions along the axis. The detailed parametric study shows that the liquid inlet velocity and mixing tube diameter have a significant impact on the performance of the Type I ejectors, and the performance of the Type II ejectors is mainly controlled by the secondary liquid flow inlet temperature.

42 ENGINEERING↗
DOE OSTI2018

A two-phase mixing layer between parallel gas and liquid streams: multiphase turbulence statistics and influence of interfacial instability

The two-phase mixing layer formed between parallel gas and liquid streams is an important fundamental problem in turbulent multiphase flows. The problem is relevant to many industrial applications and natural phenomena, such as air-blast atomizers in fuel injection systems and breaking waves in the ocean. The velocity difference between the gas and liquid streams triggers an interfacial instability which can be convective or absolute depending on the stream properties and injection parameters. In the present study, a direct numerical simulation of a two-phase gas–liquid mixing layer that lie in the absolute instability regime is conducted. A dominant frequency is observed in the simulation and the numerical result agrees well with the prediction from viscous stability theory. As the interfacial wave plays a critical role in turbulence transition and development, the temporal evolution of turbulent fluctuations (such as the enstrophy) also exhibits a similar frequency. To investigate the statistical response of the multiphase turbulence flow, the simulation has been run for a long physical time so that time-averaging can be performed to yield the statistically converged results for Reynolds stresses and the turbulent kinetic energy (TKE) budget. An extensive mesh refinement study using from 8 million to about 4 billions cells has been performed. The turbulent dissipation is shown to be highly demanding on mesh resolution compared with other terms in TKE budget. The results obtained with the finest mesh are shown to be close to converged results of turbulent dissipation which allow us to obtain estimations of the Kolmogorov and Hinze scales. The estimated Kolmogorov scale is found to be similar to the cell size of the finest mesh used here. The computed Hinze scale is significantly larger than the size of droplets observed and does not seem to be a relevant length scale to describe the smallest size of droplets formed in atomization.

Mechanics↗
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 OSTI2015

Liquid jet breakup regimes at supercritical pressures

Previously, a theory has been presented that explains how discrete vapor–liquid interfaces become diminished at certain high-pressure conditions in a manner that leads to well known qualitative trends observed from imaging in a variety of experiments. Rather than surface tension forces, transport processes can dominate over relevant ranges of conditions. In this paper, this framework is now generalized to treat a wide range of fuel-oxidizer combinations in a manner consistent with theories of capillary flows and extended corresponding states theory. Different flow conditions and species-specific molecular properties are shown to produce distinct variations of interfacial structures and local free molecular paths. These variations are shown to occur over the operating ranges in a variety of propulsion and power systems. Despite these variations, the generalized analysis reveals that the envelope of flow conditions at which the transition from classical sprays to diffusion-dominated mixing occurs exhibits a characteristic shape for all liquid–gas combinations. As a result, for alkane-oxidizer mixtures, it explains that these conditions shift to higher pressure flow conditions with increasing carbon number and demonstrates that, instead of widely assumed classical spray atomization, diffusion-dominated mixing may occur under relevant high-pressure conditions in many modern devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Void fraction, bubble size and interfacial area measurements in co-current downflow bubble column reactor with microbubble dispersion

Micro-bubbles dispersed in bubble column reactors have received great interest in recent years, due to their small size, stability, high gas-liquid interfacial area concentrations and longer residence times. The high gas-liquid interfacial area concentrations lead to high mass transfer rates compared to conventional bubble column reactors. In the present work, experiments have been performed in a down-flow bubble column reactor with micro-bubbles generated and dispersed by a novel mechanism to determine the gas-liquid interfacial area concentrations by measuring the void fraction and bubble size distributions. Gamma-ray densitometry has been employed to determine the axial and radial distributions of void fraction and a high speed camera equipped with a borescope is used to measure the axial and radial variations of bubble sizes. Also, the effects of superficial gas and liquid velocities on the two-phase flow characteristics have been investigated. Further, reconstruction techniques of the radial void fraction profiles from the gamma densitometry's chordal measurements are discussed and compared for a bubble column reactor with dispersed micro-bubbles. The results demonstrate that the new bubble generation technique offers high interfacial area concentrations (1,000 to 4,500 m 2 /m 3 ) with sub-millimeter bubbles (500 to 900 µm) and high overall void fractions (10% – 60%) in comparison with previous bubble column reactor designs. The void fraction data was analyzed using slip velocity model and empirical correlation has been proposed to predict the Sauter mean bubble diameter.

09 BIOMASS FUELS↗