Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “capillary 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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Three-Dimensional Free Gas Flow Focuses Basin-Wide Microbial Methane to Concentrated Methane Hydrate Reservoirs in Geological System

Here, we present a systematic model that links the generation, migration, phase partitioning, and accumulation of methane into a closed loop as the sediment is deposited from the seafloor and buried through the base of hydrate stability zone (BHSZ). In our model, methane is generated by biodegradation of organic carbon in muds. Hydrate does not form and methane is not trapped until a coarse-grained layer is deposited, because the small pores prevent hydrate formation in muds. Instead, methane diffuses into sands/silts where methane solidifies into hydrate. As hydrate-bearing sands/silts pass through the BHSZ during sediment burial, methane hydrate dissociates, and releases free gas. The released and the newly generated free gas below the BHSZ concentrates into a vertical/dipping zone with low capillary entry pressure and high permeability and flows upward driven its buoyancy. When free gas reaches the hydrate stability zone (HSZ), capillary forces drive free gas to flow laterally, preferentially enter sands/silts, feed hydrate growth, and elevate hydrate saturation. With three-dimensional focused free gas flow, microbial methane that is generated from a much larger fetch area of the entire basin, both above and below the BHSZ, is concentrated into coarse-grained layers at structural closures for hydrate formation. Our model illustrates how geological evolution, microbial methane generation, and gas flow by buoyancy couple to generate concentrated hydrate deposits in geological system. These insights can be used to explore for high-concentration methane hydrate and are important for understanding the methane budget and carbon cycle under the seafloor.

58 GEOSCIENCES↗
DOE OSTI2024

An integrated approach to derive relative permeability from capillary pressure

Surface tension affects all aspects of fluid flow in porous media. Through measurements of surface tension interaction under multiphase conditions, a relative permeability curve can be determined. Relative permeability is a numerical description of the interaction between two or more fluids and the porous media. It is a critical parameter for various tools that characterize subsurface multiphase flow systems, such as numerical simulation for carbon sequestration, oil and gas development, and groundwater contamination remediation. Therefore, it is critical to get a good statistical distribution of relative permeability in the porous media under study. Empirical formula for determining relative permeability from capillary pressure are already well established but do not provide the needed flexibility that is required to match laboratory-derived relative permeability curves. By expanding the existing methods for calculating relative permeability from capillary pressure data, it is possible to create both two and three-phase relative permeability curves. Mercury intrusion capillary pressure (MICP) data from the Morrow 'B' Sandstone coupled with interfacial tension and contact angle measurements were used to create a suite of relative permeability curves. Furthermore, these curves were then calibrated to a small sample of existing laboratory curves to elucidate common fitting parameters for the formation that were then used to create relative permeability curves from MICP data that does not have an associated laboratory-measured relative permeability curve.

58 GEOSCIENCES↗
DOE OSTI2024

Thermodynamic Analysis of Nano-Electrospray Induced Gas Jets

Nano-electrospray (nES) produces a plume of charged liquid droplets which have drag interactions with the surrounding gas after they are emitted from a source capillary. The resulting induced gas flow has often been considered unimportant, but we have demonstrated that a gas microjet of significant velocity up to 10s of m/s can be produced. In this work we present a thermodynamic framework that enables analysis of gas jet generation from electrosprays and introduce the important metrics for such analysis, effectiveness (a measure of momentum transfer from the electrosprayed aerosol to gas) and efficiency (a measure of energy conversion from electrical energy generating the electrospray to gas kinetic energy). This analytical framework is applicable to any two-phase flows consisting of discrete conservative-force-driven particles which exchange momentum with an inert, otherwise quiescent fluid medium to yield a co-flowing two-phase jet. Here, we apply this framework to sprays of water from nano-electrospray emitters to demonstrate that increasing the applied electrical potential difference, increasing liquid mass flowrate, and decreasing droplet size all can increase electrospray induced gas jet strength, but only the latter two do so while also increasing the momentum transfer effectiveness and energy conversion efficiency.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2024

Optimization of lithium vapor box divertor evaporator location on NSTX-U using SOLPS-ITER

Commercial fusion reactors will be faced with extremely high divertor target heat fluxes that will require mitigation. Simulations of detachment in an NSTX-U scenario projected to have 92 MW/m 2 unmitigated peak target heat flux are presented, which reaches sub-10 MW/m 2 target heat flux using a highly dissipating lithium vapor box divertor design. The lithium vapor box is a detached divertor design which employs lithium vapor evaporation and condensation to contain lithium below the X-point. Previous SOLPS modeling has indicated a lithium vapor box can reduce the heat flux down to 10 MW/m 2 via simultaneous evaporation from the Private Flux Region (PFR) and the Common Flux Region (CFR) sides of the vapor box. It is found here that PFR evaporation has improved access to the separatrix leading to significantly more efficient power dissipation than CFR evaporation. Simulations of target evaporation with an evaporation distribution that is self-consistent with the temperature of a Capillary Porous System with Fast flowing liquid lithium (CPSF) could reach n Li / n e ~ 0.025-0.030 at the Last Closed Flux Surface (LCFS) depending on the liquid metal flow speeds and lithium sputtering yield, while PFR-side evaporation can reach acceptable heat fluxes with n Li / n e ~ 0.038 at the LCFS. However, PFR evaporator performance can be improved if the target is allowed to be hot enough such that it reflects lithium, reaching n Li / n e ~ 0.028 and reducing required lithium evaporation. Ultimately PFR evaporation and target evaporation are found to have similar ability to produce acceptable heat flux solutions with minimal upstream concentration.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Relative permeabilities for two-phase flow through wellbore cement fractures

Multiple fluids are likely to exist in fractures and flow paths associated with leaky wellbores, including liquids (e.g., crude oil) and gases (e.g., gas exsolved from liquid). These fluids occupy and move through different portions of the pore spaces within the fractures depending on many factors, including fluid properties, fracture size, and the amount of the different fluids. Upward leakage of any phase, through the fracture, can contaminate water-bearing formations, create hazardous surface conditions, and compromise the functionality of the wellbore. Early signs of wellbore leaks may be expressed by anomalous pressure behavior at surface monitoring points on cavern storage wells. These pressure anomalies are difficult to interpret, necessitating knowledge of the factors that affect the multiphase flow in fractures and porous media. These parameters are critical to modeling multiphase flow in fractures. This insight can guide further diagnosis and maximize leak remediation. Here, our study focuses on the relationship of the liquid–gas relative permeabilities for representative variable-aperture wellbore cement fracture. To obtain the relative permeability of each phase, two-phase flow tests were conducted where both fluids were flowing simultaneously through a fractured wellbore cement specimen under a range of factors, namely (1) aperture size, (2) capillary numbers, and (3) viscosity ratio. The flow experiments were conducted under a range of confining stresses and flow velocities, using nitrogen gas and silicone oils (of different viscosities) in a specially designed pressure vessel. The sum of gas and oil relative permeabilities were found to be less than one under all conditions, which indicates that the presence of one phase affects the permeability of the other phase, and vice versa. Since the gas phase flow conditions include a significant inertial flow component in addition to viscous flow, the inertial flow coefficients at different saturation states are presented. The factors affecting the relationship between the relative permeabilities are discussed in detail. A new mathematical model for estimating the relative permeability of wellbore cement fracture is presented and experimentally validated.

58 GEOSCIENCES↗
DOE OSTI2022

A 1D–0D–3D coupled model for simulating blood flow and transport processes in breast tissue

In this work, we present mixed dimensional models for simulating blood flow and transport processes in breast tissue and the vascular tree supplying it. These processes are considered, to start from the aortic inlet to the capillaries and tissue of the breast. Large variations in biophysical properties and flow conditions exist in this system necessitating the use of different flow models for different geometries and flow regimes. In total, we consider four different model types. First, a system of 1D nonlinear hyperbolic partial differential equations (PDEs) is considered to simulate blood flow in larger arteries with highly elastic vessel walls. Second, we assign 1D linearized hyperbolic PDEs to model the smaller arteries with stiffer vessel walls. The third model type consists of ODE systems (0D models). It is used to model the arterioles and peripheral circulation. Finally, homogenized 3D porous media models are considered to simulate flow and transport in capillaries and tissue within the breast volume. Sink terms are used to account for the influence of the venous and lymphatic systems. Combining the four model types, we obtain two different 1D–0D–3D coupled models for simulating blood flow and transport processes: The first model results in a fully coupled 1D–0D–3D model covering the complete path from the aorta to the breast combining a generic arterial network with a patient specific breast network and geometry. The second model is a reduced one based on the separation of the generic and patient specific parts. The information from a calibrated fully coupled model is used as inflow condition for the patient specific sub-model allowing a significant computational cost reduction. Several numerical experiments are conducted to calibrate the generic model parameters and to demonstrate realistic flow simulations compared to existing data on blood flow in the human breast and vascular system. Moreover, we use two different breast vasculature and tissue data sets to illustrate the robustness of our reduced sub-model approach.

multiscale blood flow models↗
DOE OSTI2023

Modeling Geologic Waste Repository Systems Below Residual Saturation

The heat generated by high-level radioactive waste can pose numerical and physical challenges to subsurface flow and transport simulators if the liquid water content in a region near the waste package approaches residual saturation due to evaporation. Here, residual saturation is the fraction of the pore space occupied by liquid water when the hydraulic connectivity through a porous medium is lost, preventing the flow of liquid water. While conventional capillary pressure models represent residual saturation using asymptotically large values of capillary pressure, here, residual saturation is effectively modeled as a tortuosity effect alone. Treating the residual fluid as primarily dead-end pores and adsorbed films, relative permeability is independent of capillary pressure below residual saturation. To test this approach, PFLOTRAN is then used to simulate thermal-hydrological conditions resulting from direct disposal of a dual-purpose canister in unsaturated alluvium using both conventional asymptotic and revised, smooth models. Importantly, while the two models have comparable results over 100 000 years, the number of flow steps required is reduced by approximately 94%.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

A Dynamic Pore Network Model for Imbibition Simulation Considering Corner Film Flow

Wetting films can develop in the corners of angular pores under strong wetting conditions. Modeling the dynamics of corner film remains elusive using direct numerical simulations because of the significant scale difference between main meniscus and corner film flow. In this paper, the modified interacting capillary bundle model (ICB), developed in our previous work to describe accurately corner film dynamics in a single square tube, is incorporated into a single-pressure dynamic pore network model (DPNM) to simulate imbibition in strongly wetting porous media with corner film flow. The traditional pore network is decomposed into several layers of interacting subpore networks where the 0th layer of subpore network simulates the main meniscus flow and higher layers the corner film flow. The fluid flow between different layers is captured by interlayer throats. In addition, the snap-off mechanism caused by the thickening of wetting corner film is considered. The accuracy of the developed model is validated for four cases: spontaneous imbibition in a single square tube, wetting fluid redistribution through corner films under a capillary pressure difference, snap off in a narrow throat connecting two large pores, and imbibition dynamics in a real microfluidic porous geometry. The validated model is then used to simulate both spontaneous and controlled imbibition in a pore network with random pore size distribution. Finally, the interaction between corner film and main meniscus flow in porous media is analyzed from a pore-scale perspective.

58 GEOSCIENCES↗
DOE OSTI2022

A remarkable CHF of 345W/cm 2 is achieved in a wicked-microchannel using HFE-7100

Boiling heat transfer of dielectric fluids is a promising cooling technique for thermal management of microelectronic systems. However, the critical heat flux achievable is generally low, because of the poor thermophysical properties of these fluids. Here, to address this dilemma, we propose a new cooling concept to substantially enhance global liquid supply during phase-change process using enhanced capillary-driven force. Additionally, dedicated vapor pathways are designed among a bank of micro-pillars to facilitate vapor removal. In this work, new wicks comprised of silicon micro-pinfin arrays are explored to significantly enhance the flow boiling heat transfer performance. To examine the functionalities of this wick, experiments on HFE-7100 were carried out with mass velocities varying from 247 to 3,465 kg/m 2 s. To explore the enhancement mechanisms and to analyze the capillary-assisted flow boiling process, visualization studies were conducted. The results indicate that sustainable evaporation induced by wick microstructures and efficient liquid supply are the enhancement mechanisms compared to parallel microchannels with solid walls. It is found that the overall heat transfer coefficient is substantially increased up to 75%. Remarkably, a high critical heat flux (CHF) of approximately 345 W/cm 2 is recorded at G = 3,465 kg/m 2 s at coolant inlet temperature of ~20 °C. Equally importantly, this noticeable enhancement of CHF value is associated with drastically decreased pressure drops compared to microchannels decorated with μ-pinfin fences.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Model Validation and Uncertainty Quantification on the KRUSTY Microreactor Design Using GRIFFIN Neutron Transport Code [Poster]

Argonne National Laboratory (ANL) and INL have developed a GRIFFIN steady state neutronics model for the multiphysics simulations of the Kilopower Reactor Using Sterling TechnologY (KRUSTY) microreactor in the Multiphysics Object Oriented Simulation Environment (MOOSE). The reliability of such deterministic neutronics models can be validated by comparing with computations from Monte Carlo codes (e.g. MCNP, SERPENT, OpenMC, Shift, etc). Furthermore, potential modeling/design improvements can be identified by incorporating uncertainty quantification (UQ), which can be performed by MOOSE’s Stochastic Tools Module (STM). KRUSTY is a prototype for a 5-kW thermal nuclear-powered space reactor. Its primary components consist of nuclear fuel, heat pipes, a control rod, a reflector, and the shielding. The fuel consists of 3 stacked U-7.65Mo cylinders with a hole in the center for the control rod. 8 liquid sodium heat pipes transfer fission energy from the solid fuel block to the Sterling power conversion system where the energy is extracted, and the cooled sodium flows back to the core via capillary action . The movable Boron Carbide control rod regulates the neutron population during startup or when a reactor temperature boost is needed . The beryllium oxide reflector is in 3 places in the reactor; it surrounds the core axially, it lies beneath the core on a platen, and it is present in the shim. The axial and lower reflectors rest on an adjustable stainless-steel platen that moves upward to cover the fuel and help the reactor reach criticality. Lastly, radial stainless steel surrounds the core offering protection from radiation exposure .

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

High Performance Heat Pipe Power Transient Testing at SPHERE Facility

Microreactors are being researched, designed, and built at Idaho National Laboratory (INL). Microreactors are small reactors defined at less than 20MW of power. These reactor concepts are also being looked at throughout industry for various applications. An important aspect of these reactor designs is economic feasibility i.e. lower overnight capital cost. The driving factors for implementing microreactors are quick setup and takedown, minimal operators, and the ability to manufacture them readily and to fit in mid-sized containers for transport. A specific area of research to aid in successful integration of these factors within the designs is passive heat removal of the core’s thermal power. Interest in heat pipes to achieve this passive heat removal has been shown across multiple industry partners. Because of this interest, INL has developed a test facility to facilitate experimental tests for sodium filled heat pipes. INL has developed the Single Primary Heat Extraction and Removal Emulator (SPHERE) facility to run experiments on high performance, sodium filled heat pipes. As mentioned above, heat pipes are passive heat transfer devices. Radially, heat pipes are broken up into an outer wall, a small annular gap, a wick structure, and a centerline gap. They function by utilizing latent heat transfer. Heat pipes are traditionally separated into three regions, an evaporator (heat input), an adiabatic region, and finally a condenser region (heat removal). As heat is being applied to the evaporator, the working fluid undergoes a phase change to a vapor. This phase change causes a differential pressure across the axial length of the pipe driving flow down the center gap of the heat pipe. The vapor flows down past the adiabatic region to the condenser where the heat is removed. This heat removal forces the working fluid to phase change back to a liquid. The wick structure is then utilized to drive the flow back towards the evaporator by capillary forces. This backflow is aided by the annular gap. Because this heat transfer mechanism functions with latent heat transfer, the heat pipe is close to isothermal down the axial length. Heat pipes can operate under a wide range of working fluids. Considerations for these working fluids are primarily driven by operating temperatures amongst other important factors based around overall performance. Sodium filled heat pipes operate from 450°C up to 900°C. This temperature range works well for the current microreactor designs. In conjunction with this experimental capability, INL has developed a modeling software to simulate heat pipe physics within reactor cores. This modeling software is called Sockeye and functions under the established INL Multiphysics Object Oriented Simulation Environment (MOOSE). SPHERE also supports Sockeye development by providing the modeling team with experimental data on an array of setups and operating parameters to support validation efforts. A power transient experiment was performed utilizing the SPHERE facility to continue to aid with Sockeye development. The testing followed a proposed test plan to ramp up and down the temperature of the heat pipe. Sockeye models steady state heat pipe operation with high accuracy, the data provided by the power transient testing aims to assist with the validation efforts and further enhance transient modeling capability of the tool [2].

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

The liquid-conduction, vapor-flow heat pipe model in Sockeye

A single-phase heat pipe flow model implemented in the heat pipe application Sockeye is described. This model solves one-dimensional, compressible flow equations for the vapor phase in the center of a heat pipe, which are coupled to the two-dimensional heat conduction equation for the wick, liquid, and cladding, as well as to an ordinary differential equation tracking the working fluid inventory in the evaporator section of the heat pipe. This model is demonstrated with several test problems, including comparisons to analytic solutions for the vapor flow fields, analytic curves for sonic and capillary limitations of heat pipe operation, and some experimental data. The numerical solution gives excellent agreement for verification problems and good agreement with experimental results. Also, demonstrations show that the model is very robust, allowing for full simulations of heat pipe transients, including frozen startup, sonic-limited (supersonic) flow, and heat pipe shutdown.

42 ENGINEERING↗
DOE OSTI2023

Rapid Droplet Sampling Interface for Low-Volume, High-Throughput Mass Spectrometry Analysis

Here, we present a rapid droplet sampling interface (RDSI) electrospray ionization mass spectrometry (ESI-MS) system as a high-throughput, low-volume, noncontact, and minimal-carryover approach for characterization of liquids. Liquid characterization was achieved by combining droplet ejection with an open-face microflow capillary with a 2.5 μL/min continuous flow of carrier solvent. Through this implementation, single 0.3 nL droplets containing the analyte effectively mix with 4–8 nL of carrier solvent and create a combined electrospray plume. The carrier solvent continuously cleaned the system, eliminating carryover. A sampling rate of 5 Hz was achieved for droplets containing 1 μM propranolol or 5 μM leu-enkephalin with each droplet fully baseline-resolved (138 ± 32 ms baseline-to-baseline). Using a SCIEX API4000 mass spectrometer, a lower limit of quantification (LLOQ) of propranolol was 15 nM, corresponding to 1.16 fg of propranolol in the droplet, and was linear across 3 orders of magnitude. Further, quantitation could be achieved by adding an isotopically labeled internal standard, as done in conventional ESI. Signal transients were faster than the acquisition speed of the mass spectrometer, resulting in artificially high reproducibility of 15–30% RSD droplet-to-droplet. Analyte–solvent mixing ratios could be controlled by adjusting droplet positioning along the open-face capillary with an optimal position about 0.4 mm from the tip end. The range of analyte coverage was exemplified by measures of peptides and drugs in methanol, water, and buffer solutions. In a comparison to the Open Port Sampling Interface (OPSI) implemented on the same system, the RDSI had 78× greater sensitivity, 6× greater throughput and used significantly less carrier solvent.

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

Capillary Pressure to Relative Permeability: Task 7.3.4 Relative Permeability Analysis Final Report

Surface tension affects all aspects of fluid flow in porous media. Through measurements of surface tension interaction under multiphase conditions, a relative permeability relationship can be determined. Relative permeability is a numerical description of the interplay between two or more fluids and the porous media they flow through. It is a critical parameter for various tools used to characterized subsurface multiphase flow systems, such as numerical simulation for oil and gas development, carbon sequestration, and groundwater contamination remediation. Therefore, it is critical to get a good statistic distribution of relative permeability in the porous media under study. Empirical relationships for determining relative permeability from capillary pressure are already well established but do not provide the needed flexibility required to match laboratory derive relative permeability relationships. By expanding the existing methods for calculating relative permeability from capillary pressure data, it is possible to create both two and three-phase relative permeability relationships. Existing laboratory-measured relative permeability data along with mercury intrusion capillary (MICP) data coupled with interfacial tension and contact angle measurements were used to determine the efficacy of this approach to relative permeability curve creation. The relative permeability relationships determined with this method were fit to existing laboratory data to elucidate common fitting parameters that were then used to create relative permeability relationships from MICP data that does not have an associated laboratory-measured relative permeability relationship.

58 GEOSCIENCES↗
DOE OSTI2025

Manufacturing of Continuous Microporous Polymeric Fibers for Separation Applications

Polyolefin plastic waste, particularly polypropylene, is one of the most prevalent components in the plastic waste stream generated globally. However, only a small fraction of this waste is recycled and reintegrated into second applications. These materials have significant embedded energy that could be used for additional productivity if recycled properly. Therefore, there is a critical need to develop scalable processing techniques with high throughput to effectively recycle it. We report the continuous manufacturing of microporous fibers using waste polypropylene plastic and melt-processable lignin, a plant biomass constituent produced as a byproduct in paper mills or biorefineries, as a template to create pores. Further, these filaments with hierarchical porosity─created by controlled microphase separation during high-speed extensional flow followed by removal of lignin─exhibit exceptional capillary action for hydrophobic liquids. The resulting porous fibers can be used for various separation applications. For example, the oil uptake of the fibers is 9.59 ± 0.82 g/g for applications in oil recovery in bodies of water. In addition, polyethylenimine infiltration within these nanoporous fibers introduces cyclic room temperature sorption and thermal desorption potential of acidic gases such as CO 2 in a bench-scale experiment. The fibers exhibit multiple cycles of CO 2 absorption, with a range of 0.15 to 0.17 mmol/g. Thus, these polyethylenimine-infiltered nanoporous polypropylene fibers can also be used in removing acidic gases from a gas mixture.

36 MATERIALS SCIENCE↗
DOE OSTI2022

The Conventional Gas Diffusion Electrode May Not Be Resistant to Flooding during CO 2 /CO Reduction

The electrochemical CO 2 or CO reduction to chemicals and fuels using renewable energy is a promising way to reduce anthropogenic carbon emissions. The gas diffusion electrode (GDE) design enables low-carbon manufacturing of target products at a current density (e.g., 500 mA cm –2 ) relevant to industrial requirements. However, the long-term stability of the GDE is restricted by poor water management and flooding, resulting in a significant hydrogen evolution reaction (HER) within almost an hour. The optimization of water management in the GDE demands a thorough understanding of the role of the gas diffusion layer (GDL) and the catalyst layer (CL) distinctively. In this work, the hydrophobicity of the GDL and CL is independently adjusted to investigate their influence on gas transport efficiency and water management. The gas transport efficiency is more enhanced with the increase in hydrophobicity of the GDL than the CL. Direct visualization of water distribution by optical microscope and micro-computed tomography demonstrates that the water flow pattern transfers from the stable displacement to capillary fingering as GDL hydrophobicity increases. Unfortunately, only increasing the hydrophobicity is not sufficient to prevent flooding. A revolutionary change in the design of the GDE structure is essential to maintain the long-term stability of CO 2 /CO reduction.

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

Continuous Counter‐Current Microfluidic Liquid–Liquid Extraction Achieved Using a Pair of Wettable Screen Meshes

Continuous counter‐current microfluidic liquid–liquid extraction performs separations by flowing immiscible liquids in opposing directions within a single flow channel. In principle, this flow arrangement enables a large number of theoretical separation units in a small footprint, without using interstage valving, pumping, and phase separation. Despite its potential for excellent separation performance, this microfluidic scheme rarely appears in literature due to the requirement for capillary forces to be greater than hydrodynamic forces for stable flow. We present a novel microfluidic device and flow approaches that overcome this force‐balance challenge, enabling stable, long‐duration continuous counter‐current flow. Additionally, we cover a suite of methodologies for quantifying the performance of the microfluidic device, revealing the number of theoretical equilibrium stages achieved. The enabling technologies include a woven mesh screen‐based microfluidic device architecture that is easily fabricated outside of a clean room, surface functionalization strategies to promote conjugate (organic/aqueous) wettability, flow approaches to eliminate bubbles and carryover, and computer‐aided flow automation with optical measurement of extraction performance. The reported experiments lasted for over 36 h, terminated only at experiment conclusion, where the device still exhibited good performance. Automated Raman spectroscopy was used for solute quantitation of the ternary system tert‐butanol in a toluene/water matrix, a ternary system that was specifically chosen to analyze the device's performance with a small solute partition ratio and to enable in‐line Raman measurements of solute concentrations in both phases. The microfluidic device possessed a 55 mm contact length and a 38.5 µL internal volume. During counter‐current flow, we observed approximately 37 equilibrium stages (37 ± 13) based on a best‐fit of the solute fraction remaining in the aqueous phase using a Kremser Group Method analysis.

36 MATERIALS SCIENCE↗
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↗