Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multiphase flow”

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 361 records · Page 20

Addressing Two-Phase Flow Maldistribution in Microchannel Heat and Mass Exchangers

Multiphase flow phenomena in single micro and minichannels have been widely studied. Characteristics of two-phase flow through a large array of microchannels are investigated here. An air–water mixture is used to represent the two phases flowing through a microchannel array representative of those employed in practical applications. Flow distribution of the air and water flow across 52 parallel microchannels of 0.4 mm hydraulic diameter is visually investigated using high-speed photography. Two microchannel configurations are studied and compared, with mixing features incorporated into the second configuration. Slug and annular flow regimes are observed in the channels. Void fractions and interfacial areas are calculated for each channel from these observations. The flow distribution is tracked at various lengths along the microchannel array sheets. Statistical distributions of void fraction and interfacial area along the microchannel array are measured. The design with mixing features yields improved flow distribution. Void fraction and interfacial area change along the length of the second configuration, indicating a change in fluid distribution among the channels. The void fraction and interfacial area results are used to predict the performance of different microchannel array configurations for heat and mass transfer applications. Results from this study can help inform the design of compact thermal-fluid energy systems.

Engineering↗

Implementing General Framework in MFIX for Radiative Heat Transfer in Gas-Solid Reacting Flows (Final Report)

The current knowledge gap between radiation and multiphase flow communities is the greatest barrier to developing predictive models for radiation, the dominant mode of heat transfer in gas-solid reacting flows. The objectives of this research are (1) to develop and implement a general framework to support the integration of modern radiation models for gas-solid reacting flows; (2) to implement a methodology for developing new multiphase radiation models with accuracy and efficiency commensurate with their importance in energy-related applications, and (3) to demonstrate the accuracy and efficiency of the radiation models under typical gas-solid reacting flow conditions.

20 FOSSIL-FUELED POWER PLANTS↗

Hydrogen, Methane, Brine Flow Behavior, and Saturation in Sandstone Cores During H 2 and CH 4 Injection and Displacement

Large-scale underground hydrogen storage (UHS) is a critical component in the emerging hydrogen economy. Knowledge of multiphase flow behavior involving hydrogen in storage reservoir formations is crucial to characterizing hydrogen transport properties and essential for the deliverability and storage operations of UHS. There are still many gaps in fully understanding hydrogen–methane–brine multiphase phase flow that require further investigation. In this work, H 2 and CH 4 were injected through brine-saturated sandstone cores using a tri-axial core holder system fitted with flow rate meters and pressure transducers, while the effluent gas concentrations were analyzed using an online micro gas chromatograph. Brine displacement, permeability, and gas breakthrough curves were measured. We studied the flow behavior of hydrogen and methane in sandstone cores through testing brine displacement by gas injection and comparing the hydrogen displacement of methane with the methane displacement of hydrogen. We also tested the differences between horizontal and vertical flow in brine displacement. The results showed that brine displacement was more efficient in a core with higher permeability and porosity, resulting in a higher initial gas saturation. A higher gas injection rate brought about faster gas breakthrough measured by pore volume and sharper concentration curves. Hydrogen did not exhibit abnormal flow in the sandstone when the flow was horizontal and downward vertical. Gas overriding was observed in brine displacements when the flow was horizontal, with hydrogen showing this behavior more profoundly compared to methane. Downward vertical gas injection induced higher efficiency brine displacement compared to horizontal displacement and resulted in a higher initial gas saturation in the sandstone cores. These findings address critical knowledge gaps regarding gas flow patterns and displacement behaviors during hydrogen injection and recovery phases in UHS facilities using methane as the cushion gas. The insights from this research offer valuable guidance for optimizing UHS systems, ensuring operational efficiency, and advancing sustainable energy solutions in alignment with decarbonization goals.

08 HYDROGEN↗

Smoothed dissipative particle dynamics with angular momentum conservation

Smoothed dissipative particle dynamics (SDPD) combines two popular mesoscopic techniques, the smoothed particle hydrodynamics and dissipative particle dynamics (DPD) methods, and can be considered as an improved dissipative particle dynamics approach. Despite several advantages of the SDPD method over the conventional DPD model, the original formulation of SDPD by Español and Revenga (2003) [9], lacks angular momentum conservation, leading to unphysical results for problems where the conservation of angular momentum is essential. To overcome this limitation, we extend the SDPD method by introducing a particle spin variable such that local and global angular momentum conservation is restored. The new SDPD formulation (SDPD+a) is directly derived from the Navier–Stokes equation for fluids with spin, while thermal fluctuations are incorporated similarly to the DPD method. We test the new SDPD method and demonstrate that it properly reproduces fluid transport coefficients. Also, SDPD with angular momentum conservation is validated using two problems: (i) the Taylor–Couette flow with two immiscible fluids and (ii) a tank-treading vesicle in shear flow with a viscosity contrast between inner and outer fluids. For both problems, the new SDPD method leads to simulation predictions in agreement with the corresponding analytical theories, while the original SDPD method fails to capture properly physical characteristics of the systems due to violation of angular momentum conservation. In conclusion, the extended SDPD method with angular momentum conservation provides a new approach to tackle fluid problems such as multiphase flows and vesicle/cell suspensions, where the conservation of angular momentum is essential.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Shock-driven dispersal of a corrugated finite-thickness particle layer

A research area emerging in the multiphase flow community is the study of shock-driven multiphase instability (SDMI), a gas–particle analog of the traditional fluid-fluid Richtmyer–Meshkov instability (RMI). In this work, we study the interaction of planar air shocks with corrugated glass particle curtains through the use of numerical simulations with an Eulerian–Lagrangian approach. One objective of this study is to compare the simulated particle curtains to a comparable set of shock tube experiments performed to analyze traditional RMI of a gas curtain. The simulations are set to match the experimental shock Mach numbers and perturbation wavelengths (3.6 and 7.2 mm) while also matching the Atwood number of the experiments to the multiphase Atwood number of the simulations. Varying particle diameters are tested in the simulations to explore the impact of particle diameter on the evolution of the particle curtain. This simulation setup allows for a one-to-one comparison between RMI and SDMI under comparable conditions while also allowing for a separate study into the validity of the use of the multiphase Atwood number to compare the single-phase and multiphase instabilities. In particular, we show that the comparison depends on the diameter of the particles (thus, dependent on the Stokes number of the flow). A second objective of this study is to analyze the effect of the initial particle volume fraction on the evolution of the curtain and the behavior of the instability. This is done through analyzing the effect of the multiphase terms of the vorticity evolution equation on the vorticity deposition in SDMI. Also discussed is the effect of the particle diameter on the multiphase generation terms as well as in the baroclinic vorticity generation term in SDMI as the shock passes over the curtain.

42 ENGINEERING↗

Tomographic X-ray particle tracking velocimetry

We investigate the feasibility of in-laboratory tomographic X-ray particle tracking velocimetry (TXPTV) and consider creeping flows with nearly density matched flow tracers. Specifically, in these proof-of-concept experiments we examined a Poiseuille flow, flow through porous media and a multiphase flow with a Taylor bubble. For a full 360° computed tomography (CT) scan we show that the specially selected 60 micron tracer particles could be imaged in less than 3 seconds with a signal-to-noise ratio between the tracers and the fluid of 2.5, sufficient to achieve proper volumetric segmentation at each time step. In the pipe flow, continuous Lagrangian particle trajectories were obtained, after which all the standard techniques used for PTV or PIV (taken at visible wave lengths) could also be employed for TXPTV data. And, with TXPTV we can examine flows inaccessible with visible wave lengths due to opaque media or numerous refractive interfaces. In the case of opaque porous media we were able to observe material accumulation and pore clogging, and for flow with Taylor bubble we can trace the particles and hence obtain velocities in the liquid film between the wall and bubble, with thickness of liquid film itself also simultaneously obtained from the volumetric reconstruction after segmentation. While improvements in scan speed are anticipated due to continuing improvements in CT system components, we show that for the flows examined even the presently available CT systems could yield quantitative flow data with the primary limitation being the quality of available flow tracers.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Theoretical and computational dynamics of a compressible flow

An introduction to the theoretical and computational fluid dynamics of a compressible fluid is presented. The general topics addressed include: thermodynamics and physical properties of compressible fluids; 1D flow of an inviscid compressible fluid; shock waves; fundamental equations of the dynamics of a compressible inviscid non-heat-conducting and radiating fluid, method of small perturbations, linearized theory; 2D subsonic steady potential flow; hodograph and rheograph methods, exact solutions of 2D insentropic steady flow equations, 2D steady transonic and hypersonic flows, method of characteristics, linearized theory of 3D potential flow, nonlinear theory of 3D compressibe flow, anisentropic (rotational) flow of inviscid compressible fluid, electromagnetogasdynamics, multiphase flows, flows of a compressible fluid with transport phenomena.

Pai, Shih-I↗

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↗

Numerical Analysis of Water Harvesting from Saturated Vapor with the Electrospray Water Droplet Injection in The Solar-Driven Desalination Applications

Solar-driven seawater desalination systems provide freshwater through environmentally friendly and carbon-neutral processes. Air humidification and dehumidification desalination (HDD) systems extract water directly from moist air. This paper focuses on a new solar-driven seawater desalination system that shares some common traits with HDD systems but in which the air is eliminated from the process. Seawater is vaporized by solar thermal radiation in high-performance solar panels. Water vapor flows to a new electro-condensation chamber and passes through a series of electrosprays, which inject small nuclei of freshwater droplets into the vapor clouds. Here, due to dielectrophoresis and electrohydrodynamic flows, the vapor molecules are captured by the charged droplets, and vapor is condensed at the droplet surface. The water vapor condensation was modeled by using the multiphase flow numerical model. The CFD model, which was also experimentally validated, used SprayFoam solver modules available in the Open FOAM® CFD open-source software. In preliminary results, 802 grams/(hr-m2) (grams per hour per unit area) of water was harvested directly from moist air at 24C and 90% relative humidity. The simulation results indicated that freshwater productivity could increase up to 4,000 grams/(hr-m2) for fully saturated air conditions. This represents about 5 times the water productivity rate of current conventional HDD systems today. In addition, energy consumption could be reduced by 80 percent compared to conventional HDD systems.

Electrospray, electro-condensation, solar, droplet↗

Wettability Effects on Primary Drainage Mechanisms and NAPL Distribution: A Pore-Scale Study

The pore-scale processes governing water drainage behavior in porous media have implications for geoscience multiphase scenarios including carbon capture and storage, contaminant site remediation, oil recovery, and vadose zone processes. However, few studies report directly observed pore-scale water drainage phenomena in 3-D soils. This knowledge gap limits our ability to verify assumptions underlying existing models and develop optimal solutions. As such, this paper utilizes synchrotron X-ray microtomography to present an experimental pore-scale examination of nonaqueous phase liquid (NAPL)/water distribution along a primary drainage front as dense NAPL was injected upward into water wetting (WW) and intermediate wetting (IW) sand-packed columns. Pore-network structures were extracted from imaged data sets and mapped onto segmented NAPL/water data sets which allowed quantitative examinations of wettability impacts on (a) the extent to which NAPL fills individual pore bodies and (b) relationships between pore size and the phase occupying the pore, with both considered as a function of distance (and capillary pressure) relative to the NAPL front. These results revealed that several hypotheses treating IW sand similarly to WW sands are simplistic. IW systems exhibited a sequence of pore filling that deviated from traditional capillary pressure-based model predictions: NAPL invades smaller pores, while larger, adjacent pores are bypassed leaving multipore residual water ganglia. NAPL pore saturations were close to 1 and did not change with capillary pressure in IW systems. Overall, the results illustrate how a relatively small change in operative contact angle alters NAPL distribution during water drainage, with important implications for geoscience multiphase flow scenarios.

58 GEOSCIENCES↗

Effects of Nozzle Radial Position, Separation Angle, and Gas Flow Partitioning on the Mixing, Eye Area, and Wall Shear Stress in Ladles Fitted with Dual Plugs

Mixing time, slag eye area, and wall shear stress in a ladle fitted with dual plugs have been studied as a function of operating variables, namely, gas flow rate, radial nozzle position, separation angle between nozzles and flow rate partitioning. While the mixing time and slag eye area have been experimentally investigated in a scaled water model (with and without a top slag layer), the shear stress on the vessel wall has been studied computationally via a RANS-based turbulent, three-dimensional coupled Eulerian–Lagrangian (VOF-DPM), multiphase flow model. Experimental and computational results have indicated that the mixing, eye area, and wall shear stress depend on the gas flow rate, i.e., while mixing efficiency increases with the increasing gas flow rate, the slag eye area, and wall shear stress, in contrast, become more pronounced with the increasing gas flow, leading to an undesirable operating regime. The radial nozzle position, separation angle between the nozzles and flow rate partitioning at any given flow rate also influence the ladle process performance, influencing the mixing, eye area, and wall shear stress, albeit to a varied degree. Within the range of operating conditions studied, an expected inverse correspondence between mixing time and shear behavior was observed. Experimental and computational results in conjunction have indicated that the best arrangement of porous plugs or gas injection nozzles for superior ladle process performance is dependent on the gas flow rate and is specific to the desired objective (i.e., decreasing the mixing time, or slag eye area and wall shear stress). Hence, a unique gas injection practice cannot and should not be suggested as the optimum for ladle metallurgy in steelmaking. Nevertheless, if the mixing time is the parameter of primary interest, a nozzle configuration with equal flow partitioning (1:1) between the nozzles and identical nozzle radial positions (0.7R/0.7R, 45 deg) should suffice for both low and high gas flow rates. In contrast, a nonidentical nozzle radial position (0.7R/0.5R, 90 deg) and an unequal gas flow rate per nozzle (1:3) appear preferable if both ladle eye and shear stresses, rather than mixing, are the issues of concern.

36 MATERIALS SCIENCE↗

Three-Dimensional Bubble Fluidics in Architected Porous Media

Gas bubble flows in porous media often exhibit complex and seemingly unpredictable behaviors that are difficult to control. This lack of control limits the ability to design effective devices which manage multiphase flows. Here, we show how the design of 3D printed pores can deterministically control the flow path of an injected gas stream. Open cell structures can be designed to shape the gas/liquid interface with fidelity to control how the two phases are distributed throughout a porous material. The distributed gas volume is free to interact physically and chemically with the surrounding liquid phase, an effect we exploit to create a logical control gate to redirect flows within a lattice. This also allows us to design architectures for reactive capture and aerating bioreactors, resulting in patterned boundaries which can make more effective use of the liquid and gas reagents.

3D microfluidics↗

Construction of an Exact Pressure-Equilibrium Scheme for the Five-Equation Two-Phase Flow Model With Thermal Relaxation

Numerical simulation of compressible multiphase flows based on the four-equation (homogeneous relaxation) model is known to suffer from two fundamental difficulties with (a) wave propagation and (b) pressure equilibrium preservation. First, the mixture sound speed exhibits non-monotonic dependency with respect to the volume fraction, which leads to robustness issues in the resolution of shocks and acoustic wave propagation across two-phase regions. This difficulty can be mitigated by solving Allaire’s five-equation model augmented with infinitely fast phasic temperature equilibrium, from which solutions of the four-equation model can be recovered. However, when temperature is non-uniform, this augmented five-equation formulation still fails to preserve pressure equilibrium across material interfaces. In this work, we propose a fully conservative numerical scheme that exactly preserves pressure equilibrium at the discrete level for the augmented five-equation model, for arbitrary initial distributions of temperature and volume fraction. Combined with the monotonic sound speed property of the five-equation formulation, the proposed pressure-equilibrium preserving scheme significantly improves robustness in the presence of strong multiphase interactions, including shock–interface interactions and advection of material interfaces.

ESG↗

New insights on the complex dynamics of two-phase flow in porous media under intermediate-wet conditions

Multiphase flow in porous media is important in a number of environmental and industrial applications such as soil remediation, CO 2 sequestration, and enhanced oil recovery. Wetting properties control flow of immiscible fluids in porous media and fluids distribution in the pore space. In contrast to the strong and weak wet conditions, pore-scale physics of immiscible displacement under intermediatewet conditions is less understood. This study reports the results of a series of two-dimensional highresolution direct numerical simulations with the aim of understanding the pore-scale dynamics of twophase immiscible fluid flow under intermediate-wet conditions. Our results show that for intermediatewet porous media, pore geometry has a strong influence on interface dynamics, leading to co-existence of concave and convex interfaces. Intermediate wettability leads to various interfacial movements which are not identified under imbibition or drainage conditions. These pore-scale events significantly influence macro-scale flow behaviour causing the counter-intuitive decline in recovery of the defending fluid from weak imbibition to intermediate-wet conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

On the dynamics and kinematics of two-fluid-phase flow in porous media

A model formulated in terms of both conservation and kinematic equations for phases and interfaces in two-fluid-phase flow in a porous medium system is summarized. Macroscale kinematic equations are derived as extensions of averaging theorems and do not rely on conservation principles. Models based on both conservation and kinematic equations can describe multiphase flow with varying fidelity.When only phase-based equations are considered, a model similar in form to the traditional model for two-fluid-phase flow results. When interface conservation and kinematic equations are also included, a novel formulation results that naturally includes evolution equations that express dynamic changes in fluid saturations, pressures, the capillary pressure, and the fluid-fluid interfacial area density in a two-fluid-system. This dynamic equation set is unique to this work, and the importance of the modeled physics is shown through both microfluidic experiments and high-resolution lattice Boltzmann simulations. The validation work shows that the relaxation of interface distribution and shape toward an equilibrium state is a slow process relativeto the time scale typically allowed for a system to approach an apparent equilibrium state based upon observations of fluid saturations and external pressure measurements. Consequently, most pressure-saturation data intended to denote an equilibrium state are likely a sampling from a dynamic system under-going changes of interfacial curvatures that are not typically monitored. The results confirm the importance of kinematic analysis in combination with conservation equations for faithful modeling of system physics.

54 ENVIRONMENTAL SCIENCES↗