Search NASA⌕ Search

SEARCH · Search NASA

Results for “interfacial physics”

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

Effects of Heating Configuration and Operating Parameters on Heat Transfer and Interfacial Physics of Microgravity Flow Boiling With Subcooled Inlet Conditions –Experiments Onboard the International Space Station

This study is part of the Flow Boiling and Condensation Experiment (FBCE), a collaborative effort between the Purdue University Boiling and Two-Phase Flow Laboratory (PU-BTPFL) and the NASA Glenn Research Center. The FBCE fitted with the Flow Boiling Module (FBM) was launched to the International Space Station (ISS) in August 2021 and experiments were successfully performed from February to July 2022 to amass a large microgravity-flow-boiling database. This study is focused on heat transfer and flow visualization of microgravity flow boiling of n-Perfluorohexane in a rectangular channel of 5.0 mm height, 2.5 mm width (heated), and 114.6 mm length, with subcooled inlet conditions. High-speed-video photography is utilized to present flow patterns and temporal interfacial behavior. Heat transfer results are presented in the form of flow boiling curves and both parametric curves and streamwise profiles of wall temperature and heat transfer coefficient. Firstly, the parametric effects of mass velocity (199.4 – 3200.0 kg/m 2 s), inlet subcooling (0.2 – 46.0°C), and inlet pressure (124.2 – 176.7 kPa), on the aforementioned aspects are assessed for double-sided heating to establish them for a microgravity environment. Of these three parameters, mass velocity and inlet subcooling mostly determine the microgravity flow boiling behavior, while inlet pressure plays an insignificant role. Flow patterns for double-sided heating are more complex than those for single-sided heating due to interaction between the two vapor layers. Vapor interaction is minimized at high subcoolings and high mass velocities due to strong condensation offered by the subcooled bulk liquid layer separating them. Despite the different flow patterns, both single- and double-sided heating generally result in similar parametric trends and local heat transfer coefficients for similar operating conditions. Flow instabilities manifest as temporal flow anomalies and temperature oscillations, and their severity increases with increasing boiling number. Secondly, the effects of heating configuration are analyzed by comparing and contrasting several aspects of single- and double-sided heating data. The heat fluxes at which onset of nucleate boiling degradation (ONBD) and critical heat flux (CHF) occur are distinctly different for single- and double-sided heating. There exists a threshold inlet subcooling demarcating the dominance of flow acceleration and condensation effects in vapor removal from the near-wall region and replenishment of fresh liquid for boiling. Above the threshold, condensation from the near-wall region is dominant and single-sided heating yields higher heat fluxes, and below it, acceleration is dominant and double-sided yields higher heat fluxes. At mass velocity in the range of 200 – 2400 kg/m 2 s, the threshold inlet subcooling lies in the approximate range of 20 – 30°C (corresponding inlet quality of roughly -0.40 – -0.20).

Microgravity↗

Optical and Electronic Losses Arising from Physically Mixed Interfacial Layers in Perovskite Solar Cells

Perovskite solar cell device performance is affected by optical and electronic losses. To minimize these losses in solar cells, it is important to identify their sources. Here, we report the optical and electronic losses arising from physically mixed interfacial layers between the adjacent component materials in highly efficient two terminal (2T) all-perovskite tandem, single junction wide-bandgap, and single junction narrow-bandgap perovskite-based solar cells. Physically mixed interfacial layers as the sources of optical and electronic losses are identified from spectroscopic ellipsometry measurements and data analysis followed by comparisons of simulated and measured external quantum efficiency spectra. Parasitic absorbance in the physically mixed regions between silver metal electrical contacts and electron transport layers (ETLs) near the back contact and a physical mixture of commercial indium tin oxide and hole transport layers (HTL) near the front electrical contact lead to substantial optical loss. Here, a lower-density void + perovskite nucleation layer formed during perovskite deposition at the interface between the perovskite absorber layer and the HTL causes electronic losses due to incomplete collection of photogenerated carriers likely originating from poor coverage and passivation of the initially nucleating grains.

14 SOLAR ENERGY↗

Heat Transfer and Interfacial Flow Physics of Microgravity Flow Boiling in Single-Side-Heated Rectangular Channel with Subcooled Inlet Conditions – Experiments Onboard the International Space Station

This study is the culmination of a long-term collaborative effort between researchers from the Purdue University Boiling and Two-Phase Flow Laboratory (PU-BTPFL) and the NASA Glenn Research Center to investigate gravitational effects on flow boiling and flow condensation. The science and design concepts for this large-scale effort were initiated in 2011 and included several studies detailing various aspects of two-phase fluid physics in both Earth gravity and microgravity, culminating in construction of the large-scale experimental facility named “Flow Boiling and Condensation Experiment (FBCE)”. The experiment was launched to the International Space Station (ISS) in August 2021. Following the successful installation of FBCE, equipped with the Flow Boiling Module (FBM), onboard the ISS and completion of several safety checks, flow boiling experiments were performed for five months from February 2022 until July 2022. This resulted in a large flow boiling database covering broad ranges of operating parameters and heating configurations spanning several research objectives. This study investigates microgravity flow boiling of n-perfluorohexane with subcooled inlet in a single-side-heated rectangular channel of dimensions 114.6-mm heated length, 2.5-mm heated width, and 5.0-mm height. Key operating parameters investigated are mass velocity (199.90 – 3200.13 kg/m 2 s), inlet subcooling (0.10 – 45.76°C), and inlet pressure (113.30 – 164.29 kPa). Images and image sequences acquired via high-speed-video are presented to elucidate the interfacial flow physics. To analyze and explain the effects of various parameters in microgravity, heat transfer results are presented as flow boiling curves, streamwise profiles of wall temperature and heat transfer coefficient, and parametric trends of local and averaged heat transfer coefficient. Mass velocity and inlet subcooling significantly influenced most of the aforementioned aspects of flow boiling, whereas effects of inlet pressure were comparatively insignificant. Although the data and observed flow physics might be different, the parametric effects and trends in microgravity are similar to vertical upflow in Earth gravity. Some cases, especially low mass velocities, high heat fluxes, and large degrees of inlet subcooling, experienced temporally anomalous flow behaviors caused by two-phase flow instabilities manifesting as flow reversals and resulted in deviations in overall trends. Severe thermodynamic non-equilibrium is observed throughout the channel. Overall, FBCE’s ISS experiments were successful for subcooled inlet with single-sided heating of rectangular channel, and the collected data well established the various effects on flow boiling physics in highly controlled long-term microgravity conditions.

flow boiling↗

Data and Modeling of Dendrites Subject to A Step Change in Pressure (TDSE)

There is considerable interest in dendritic solidification because of the influence dendrites have in the determination of microstructure, and thereby in the physical properties of cast metals and alloys. Current theories and models of dendritic growth generally couple diffusion effects in the melt with the physics of the interface. Data and subsequent analysis prior of the tip growth speed and radii of thermal succinonitrile dendrites in the near-convection free, on-orbit, free-fall environment demonstrate that these theories yield predictions that are reasonably in agreement with the results of experiment. However, data and analysis for assessing the interfacial physics component of theory are not sufficiently detailed or definitive. To study fundamental aspects of dendritic interface stability, we are measuring and modeling the kinetics and morphology of dendrites as they evolve from one well-defined steady state at a pre-set supercooling, through a transient stage, to a different well-defined steady state. More specifically, we subject succinonitrile dendrites, growing under steady-state conditions, to a rapid change in pressure. This leads to a rapid change in thermal driving force from the corresponding change in both the equilibrium melting temperature due to the Clapeyron effect, and a change in the far-field temperature due to adiabatic temperature changes in the bulk liquid and solid. Subsequently, we observe transformations from a well-characterized initial state into a new steady-state. Initial data reveal that the dendrite tip velocity changes almost as fast as the pressure charges, while the tip radius changes occur more slowly, taking from 10 60 seconds depending on the size of the step change and the final supercooling. Computer modeling of this process shows both agreements and disagreements with the experimental data. In making these observations and measurements, we are gaining new understandings of interfacial dynamics and state-selection physics.

Koss, Matthew B.↗

Insights into interfacial effect and local lithium-ion transport in polycrystalline cathodes of solid-state batteries

Interfacial issues commonly exist in solid-state batteries, and the microstructural complexity combines with the chemical heterogeneity to govern the local interfacial chemistry. The conventional wisdom suggests that “point-to-point” ion diffusion at the interface determines the ion transport kinetics. Here, we show that solid-solid ion transport kinetics are not only impacted by the physical interfacial contact but are also closely associated with the interior local environments within polycrystalline particles. In spite of the initial discrete interfacial contact, solid-state batteries may still display homogeneous lithium-ion transportation owing to the chemical potential force to achieve an ionic-electronic equilibrium. Nevertheless, once the interior local environment within secondary particle is disrupted upon cycling, it triggers charge distribution from homogeneity to heterogeneity and leads to fast capacity fading. Our work highlights the importance of interior local environment within polycrystalline particles for electrochemical reactions in solid-state batteries and provides crucial insights into underlying mechanism in interfacial transport.

25 ENERGY STORAGE↗

The Transient Dendritic Solidification Experiment (TDSE)

Dendritic solidification is a common mode of solidification. It is also an important model problem in nonequilibrium physics and pattern formation physics. Current theories couple the transfer of latent heat with selection mechanisms at the interface. Measurements of succinonitrile (SCN) dendrites in microgravity show reasonable agreement between heat transfer predictions and experiment. However, data and analysis for assessing interfacial physics theories are less definitive. We are studying, and will present data on, transient effects in dendritic growth of SCN. We employ the Clapeyron pressure/melting temperature effect to make a rapid change in a sample's hydrostatic pressure, and thereby rapidly change the specimen's melting temperature, forcing the dendrite to select a new steady-state. These initial measurements show some surprising and non-intuitive effects.

Koss, M. B.↗

Molecular Dynamics Investigation of Nanoscale Hydrophobicity of Polymer Surfaces: What Makes Water Wet?

The wettability of a polymer surface–related to its hydrophobicity or tendency to repel water–can be crucial for determining its utility, such as for a coating or a purification membrane. While wettability is commonly associated with the macroscopic measurement of a contact angle between surface, water, and air, the molecular physics that underlie these macroscopic observations are not fully known, and anticipating relative behavior of different polymers is challenging. To address this gap in molecular-level understanding, we use molecular dynamics simulations to investigate and contrast interactions of water with six chemically distinct polymers: polytetrafluoroethylene, polyethylene, polyvinyl chloride, poly(methyl methacrylate), Nylon 66, and polyvinyl alcohol. We show that several prospective quantitative metrics for hydrophobicity agree well with experimental contact angles. Moreover, the behavior of water in proximity to these polymer surfaces can be distinguished with analysis of interfacial water dynamics, extent of hydrogen bonding, and molecular orientation–even when macroscopic measures of hydrophobicity are similar. The predominant factor dictating wettability is found to be the extent of hydrogen bonding between polymer and water, but the precise manifestation of hydrogen bonding and its impact on surface-water structure varies. In the absence of hydrogen bonding, other molecular interactions and polymer mechanics control hydrophobic ordering. Furthermore, these results provide new insights into how polymer chemistry specifically impacts water-polymer interactions and translates to surface hydrophobicity. Such factors may facilitate the design or processing of polymer surfaces to achieve targeted wetting behavior, and presented analyses can be useful in studying the interfacial physics of other systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Multiphase computational fluid dynamics modeling of reacting flows in absorption columns for carbon capture

First-principles derived computational fluid dynamics (CFD) simulations have been proposed as a fundamental tool for investigating solvent-based CO 2 absorption in packed columns due to their ability to accurately represent the underlying nonlinear, multiscale dynamics. Numerous studies have previously utilized such CFD simulations to investigate hydrodynamics of columns with structured and random packings by assessing the key hydrodynamic metrics such as the interfacial and wetted areas. While mapping such metrics for different conditions is essential to the optimization of absorption columns, it is not sufficient, as the CO 2 capture rate depends also on the coupled, nonlinear dynamics from the underlying chemical reaction kinetics, thermodynamics, and heat-transfer rates. In this work, we present detailed CFD simulation results augmented by incorporating the effects of interfacial physical mass transfer of CO 2 , heat release from chemical reaction kinetics, and thermophysical property variations from resulting temperature gradients. We demonstrate the applicability of the proposed approach in numerically assessing the performance of packed columns by evaluating key hydrodynamic quantities, CO 2 absorption rates, and temperature rise in a reference column with packings that are structurally similar to the Sulzer Mellapak™ 250.Y packing, for different solvent inflow velocities and temperatures. Predictions from simulation results are found to be consistent with the trends in experimental observations from the literature, suggesting that the predictive capabilities of the simulation framework can be leveraged to guide the future development of absorber-column designs and optimized process flowsheets.

Absorption columns↗

Pressurization of a Flightweight, Liquid Hydrogen Tank: Evaporation and Condensation at a Liquid Vapor Interface

Evaporation and condensation at a liquidvapor interface is important for long-term, in-space cryogenic propellant storage. Yet the current understanding of interfacial physics does not predict behavior or evaporation condensation rates. The proposed paper will present a physical model, based on the 1-D Heat equation and Schrages equation which demonstrates thin thermal layers at the fluidvapor interface.

Cryogenic hydrogen↗

Experimental Heat Transfer Results and Flow Visualization of Vertical Upflow Boiling in Earth Gravity with Subcooled Inlet Conditions – In Preparation for Experiments Onboard the International Space Station

Since 2012, researchers at the Purdue University Boiling and Two-Phase Flow Laboratory (PU-BTPFL) and NASA Glenn Research Center have been collaborating on a long-term effort to study flow boiling and condensation in microgravity. The ultimate goal has been to develop the Flow Boiling and Condensation Experiment (FBCE) for the International Space Station (ISS). Based on the findings from prior flow boiling experiments both at different orientations in Earth gravity and onboard parabolic flights simulating short durations of microgravity, a final refined experiment design, construction, and operating procedure have been arrived at for long-duration microgravity flow boiling experiments onboard the ISS. This study investigates flow boiling of n-Perfluorohexane with subcooled inlet in a rectangular channel of dimensions 114.6 mm heated length, 2.5 mm width, and 5 mm height. These pre-launch experiments (Mission Sequence Testing) were conducted in vertical upflow orientation in Earth gravity using the same experimental rig that was launched to the ISS in August 2021. The various operating parameters varied are heating configuration (single- and double-sided), mass velocity (180 – 3200 kg/m2s), inlet subcooling (+0 – 32°C, encompassing both highly subcooled and near-saturated inlet conditions), and inlet pressure (119 – 191 kPa). High-speed video flow visualization images are presented to explain the two-phase interfacial physics within the channel’s heated section. Heat transfer results in terms of flow boiling curves, streamwise profiles of wall temperature and heat transfer coefficient, and averaged heat transfer coefficients are analyzed and parametric effects elucidated. Severe temporal thermodynamic equilibrium is observed for near-saturated inlet at very low velocities. Nucleate boiling degradation starts at larger heat fluxes for single-sided heating than double sided at low mass velocities with highly subcooled inlet, and conversely at high mass velocities with near-saturated inlet. Nucleate boiling degradation can be delayed to higher heat fluxes by highly subcooling the inlet and increasing mass velocity. The entire local heat transfer coefficient profiles are degraded at higher heat fluxes for near-saturated inlet, but only the downstream part for highly subcooled inlet. This study also confirmed reliability of the upcoming ISS experimental data for subcooled inlet conditions and the collected Earth-gravity data will be used for comparison against the ISS data.

flow boiling↗

Interfacial Strength and Physical Properties of Functionalized Graphene - Epoxy Nanocomposites

The toughness and coefficient of thermal expansion of a series of functionalized graphene sheet - epoxy nanocomposites are investigated. Functionalized graphene sheets are produced by splitting graphite oxide into single graphene sheets through a rapid thermal expansion process. These graphene sheets contain approx. 10% oxygen due to the presence of hydroxide, epoxide, and carboxyl functional groups which assist in chemical bond formation with the epoxy matrix. Intrinsic surface functionality is used to graft alkyl amine chains on the graphene sheets, and the addition of excess hardener insures covalent bonding between the epoxide matrix and graphene sheets. Considerable improvement in the epoxy dimensional stability is obtained. An increase in nanocomposite toughness is observed in some cases.

Miller, Sandi G.↗

Engineering the Microstructure and Morphology of Explosive Films via Control of Interfacial Energy

Physical vapor deposition of organic explosives enables growth of polycrystalline films with a unique microstructure and morphology compared to the bulk material. This study demonstrates the ability to control crystal orientation and porosity in pentaerythritol tetranitrate films by varying the interfacial energy between the substrate and the vapor-deposited explosive. Additionally, variation in density, porosity, surface roughness, and optical properties is achieved in the explosive film, with significant implications for initiation sensitivity and detonation performance of the explosive material. Various surface science techniques, including angle-resolved X-ray photoelectron spectroscopy and multiliquid contact angle analysis, are utilized to characterize interfacial characteristics between the substrate and explosive film. Optical microscopy and scanning electron microscopy of pentaerythritol tetranitrate surfaces and fracture cross sections illustrate the difference in morphology evolution and the microstructure achieved through surface energy modification. X-ray diffraction studies with the Tilt-A-Whirl three-dimensional pole figure rendering and texture analysis software suite reveal that high surface energy substrates result in a preferred (110) out-of-plane orientation of pentaerythritol tetranitrate crystallites and denser films. Low surface energy substrates create more randomly textured pentaerythritol tetranitrate and lead to nanoscale porosity and lower density films. This work furthers the scientific basis for interfacial engineering of polycrystalline organic explosive films through control of surface energy, enabling future study of dynamic and reactive detonative phenomena at the microscale. Results of this study also have potential applications to active pharmaceutical ingredients, stimuli-responsive polymer films, organic thin film transistors, and other areas.

36 MATERIALS SCIENCE↗

Anion Effects on the Interfacial Structure and Bulk Physical Properties in Choline-Based Hydrogen-Bonded Electrolytes

Electrolytes consisting of choline (Ch) salts with acetate (Ac) and oxalate (Ox) anions were studied in mixtures with ethylene glycol (EG), in comparison to chloride (Cl)- and bis(trifluorosulfonyl)imide (TFSI)-containing systems to understand the impact of anion size and charge density on the bulk electrolyte and interfacial properties. These represent hydrogen-bonded concentrated electrolytes where the Ch salts are the H-bond acceptors (HBAs) and EG is the H-bond donor (HBD) in 1:2 and 1:4 HBA:HBD ratios. Increased H-bonding in Ch 2 Ox with EG compacts the liquid compared to ChAc, consistent with the increased dissociation energies of the HBA:HBD clusters with Ox –2 calculated by density functional theory (DFT), thus leading to denser and more viscous liquids. Electrochemical impedance spectroscopy (EIS) and surface-enhanced Raman spectroscopy (SERS) were used to investigate interfacial charge density and the electrode surface species. Differential capacitance curves demonstrate weak dependence on potential except for ChAc:EG, where Ac – is found to adsorb on the Au surface, in particular at +0.055 V vs Ag quasi-reference while no specific adsorption was observed on glassy carbon (GC). In situ SERS measurement confirms the formation of a surface layer in the ChAc:EG systems resulting from a Au–acetate interaction. Despite the structural similarities, Ox –2 did not show specific binding or any dependence of capacitance on potential as a result of strongly solvating EG molecules via H-bonds. Furthermore, this study shows that the charge density of the HBAs impacts both bulk properties and the electrical double-layer formation in concentrated H-bonded electrolytes that is relevant to electrochemical kinetics and morphology in electrodeposition processes and energy storage applications.

25 ENERGY STORAGE↗

Characterizing the Geometry and Quantifying the Impact of Nanoscopic Electrocatalyst/Semiconductor Interfaces under Solar Water Splitting Conditions

The materials that are receiving the most attention in photoelectrochemical water splitting are metallic nanoparticle electrocatalysts (np-EC) attached to the surface of a semiconductor (SC) light absorber. In these multicomponent systems, the interface between the semiconductor and electrocatalysts critically affects performance. However, the np-EC/SC interface remains poorly understood as it is complex on atomic scales, dynamic under reaction conditions, and inaccessible to direct experimental probes. This contribution sheds light on how the electrocatalyst/semiconductor interface evolves under reaction conditions by investigating the behavior of nickel electrocatalysts (as nanoparticles and films) deposited on silicon semiconductors. Rigorous electrochemical experiments, interfacial atomistic characterization, and computational modeling are combined to demonstrate critical links between the atomistic features of the interface and the overall performance. It is shown that electrolyte-induced atomistic changes to the interface lead to (1) modulation of the charge carrier fluxes and a dramatic decrease in the electron/hole recombination rates and (2) a change in the barrier height of the interface. Additionally, the critical roles of nonidealities and electrocatalyst coverage due to interfacial geometry are explored. Each of these factors must be considered to optimize the design of metal/semiconductor interfaces which are broadly applicable to photoelectrocatalysis and photovoltaic research.

36 MATERIALS SCIENCE↗

Wickability-optimized textured liquid-desiccant air dehumidifiers for independent moisture management in energy-efficient buildings

Liquid-desiccant-based air conditioning systems are envisioned to enable independent humidity management, thereby improving the energy efficiency of future buildings. Existing liquid-desiccant-based air conditioning concepts, however, suffer from a poor liquid flow distribution deteriorating moisture removal rate. They are consequently flooded with the liquid-desiccant solution, which significantly degrades the energy efficiency of the dehumidification process. Here, in this work, the capillary forces and wickability effect of textured air dehumidifier surfaces are altered to minimize the liquid-desiccant flow rate of the fully wetted state, thereby transforming the physics of interfacial desiccant flow distribution. Consequently, the wickability-optimized air dehumidifier surface maximizes both moisture removal rate and dehumidification energy efficiency. It was interestingly found that the length scale of a textured air dehumidifier surface concept is optimized at an intermediate pattern density. Dry solid-air menisci appear at length scales exceeding the optimum pattern distance while the effective liquid–air interfacial area is reduced at smaller length scales, both of which degrade the moisture removal rate. At the optimum pattern density, the effective liquid–air interfacial area increases with the solution flow rate, thereby increasing the dehumidification rate. At a water vapor pressure potential of 3 kPa and a solution flow rate of 2.8 g/s, experimental results indicated a moisture removal rate of 0.16 g/m2-s for a textured surface concept with a capillary length scale of 3 mm, a 28% improvement compared with that of smooth-plate dehumidifier surfaces. A high moisture removal rate of the textured surface at a low desiccant flow rate led to a high thermal efficiency of 0.75 at a water vapor pressure potential of 5.6 kPa and a LiBr flow rate of 2.8 g/s. The insights gained from the present study accelerate the development of advanced textured surface concepts for next-generation liquid-desiccant-based air dehumidification systems offering independent humidity management for future energy-efficient buildings.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Investigation of γ - ( U , Z r ) structural properties and its interfacial properties with liquid sodium using ab initio molecular dynamics

In this study, the elastic properties, structural parameters, sound velocity, and Debye temperature of γ–(U,Zr) were computed using ab initio molecular dynamics (AIMD) at temperatures between 1000 K and 1400 K and for Zr content between 0 at. % and 100 at. %. UZr is used as a metallic fuel for Sodium Fast Reactors (SFRs). The study of the mechanical and thermal behavior of these alloys leads to a better data-informed fuel design. The bulk modulus, shear modulus, Young's modulus, and Poisson's ratio were calculated from the elastic constants and their dependence on Zr content and temperature was investigated, comparing the results with previous computational work and the available experimental data in the literature. Interfacial properties between UZr (up to 32 at. % which typically exists in nuclear fuel) and liquid sodium are also of interest due to the presence of a sodium bond between the fuel and the cladding in metallic nuclear fuel. The interfacial energy between γ–(U,Zr) and liquid sodium, the surface tension of liquid sodium, and the work of adhesion were computed at different temperatures and Zr concentrations. It was demonstrated that γ–(U,Zr) is completely wetted by liquid sodium at all the investigated temperatures and Zr concentrations. Finally, this work provides the basis for the determination of interfacial resistances in SFRs and their implementation into heat transfer fuel performance simulations, which will be the subject of future work.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Pore-Scale Dynamics of Liquid CO2–Water Displacement in 2D Axisymmetric Porous Micromodels Under Strong Drainage and Weak Imbibition Conditions: High-Speed μPIV Measurements

Resolving pore-scale transient flow dynamics is crucial to understanding the physics underlying multiphase flow in porous media and informing large-scale predictive models. Surface properties of the porous matrix play an important role in controlling such physics, yet interfacial mechanisms remain poorly understood, in part due to a lack of direct observations. This study reports on an experimental investigation of the pore-scale flow dynamics of liquid CO 2 and water in two-dimensional (2D) circular porous micromodels with different surface characteristics employing high-speed microscopic particle image velocimetry (μPIV). The design of the micromodel minimized side boundary effects due to the limited size of the domain. The high-speed μPIV technique resolved the spatial and temporal dynamics of multiphase flow of CO 2 and water under reservoir-relevant conditions, for both drainage and imbibition scenarios. When CO 2 displaced water in a hydrophilic micromodel (i.e., drainage), unstable capillary fingering occurred and the pore flow was dominated by successive pore-scale burst events (i.e., Haines jumps). When the same experiment was repeated in a nearly neutral wetting micromodel (i.e., weak imbibition), flow instability and fluctuations were virtually eliminated, leading to a more compact displacement pattern. Energy balance analysis indicates that the conversion efficiency between surface energy and external work is less than 30%, and that kinetic energy is a disproportionately smaller contributor to the energy budget. This is true even during a Haines jump event, which induces velocities typically two orders of magnitude higher than the bulk velocity. These novel measurements further enabled direct observations of the meniscus displacement, revealing a significant alteration of the pore filling mechanisms during drainage and imbibition. While the former typically featured burst events, which often occur only at one of the several throats connecting a pore, the latter is typically dominated by a cooperative filling mechanism involving simultaneous invasion of a pore from multiple throats. This cooperative filling mechanism leads to merging of two interfaces and releases surface energy, causing instantaneous high-speed events that are similar, yet fundamentally different from, burst events. Finally, pore-scale velocity fields were statistically analyzed to provide a quantitative measure of the role of capillary effects in these pore flows.

Li, Yaofa↗