Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “capillary pressure”

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

Preconditioners for multiphase poromechanics with strong capillarity

This paper aims to enhance the performance of Newton–Krylov solvers for coupled poromechanical problems with two-phase flow. In particular, we investigate the impact of capillary pressure on preconditioning strategies. Capillarity complicates the coupling between the solid deformation and fluid pressure degrees of freedom, as well as increases the nonlinearity of the system. Depending on the capillary pressure relation used in the constitutive formulation, the flow equations may exhibit a spectrum of advection-dominated to diffusion-dominated behavior. We propose preconditioning approaches that account for this behavior and lead to robust numerical performance within a broad range of regimes.

42 ENGINEERING↗
NASA NTRS1975

Capillary flow through heat-pipe wicks

Theoretical expressions are obtained for the capillary-pressure limit and permeability of a fibrous wick in terms of the porosity and fiber diameter. Hysteresis in capillary pressure is included through the introduction of an empirical hysteresis constant. A partial-saturation model based on the statistical distribution of local porosity requires an additional empirical constant, the standard deviation. The theory is compared to results of a beta-ray absorption experiment that measured the liquid content of a partially saturated wick and to results of permeability measurements on partially and fully saturated wicks. A simple wick-weighing experiment is described that yields values for the empirical hysteresis constant and the standard deviation. Theoretical results are used to design an optimum wick.

Eninger, J. E.↗
DOE OSTI2022

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

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

Engineering↗
DOE OSTI2021

Experimental study of three-dimensional CO 2 -water drainage and fracture-matrix interactions in fractured porous media

Estimation of CO 2 storage capacity in fractured porous reservoirs requires a better understanding of the CO 2 -water flow fundamentals in fracture-matrix systems. There are few core-flood experiments on three-dimensional CO 2 -water drainage in a fracture-matrix system, yet none have examined the impacts of the capillary continuity of fractures and fracture-matrix interactions. In this study, twelve drainage experiments were conducted in four fracture-matrix columns, each comprising a vertical stack of a cylindrical rock core, a filter paper serving as an analogue to a horizontal fracture, and a ceramic plate. Here, the core sample was surrounded by open space at the top and circumferential sides to model fractures, allowing for 3-D CO 2 -water drainage in the rock core and displaced water draining across the horizontal fracture. X-ray computed tomography was conducted to visualize the dynamic invasion/drainage processes in four rock core samples showing contrasts in anisotropy, permeability, and heterogeneity. Experimental results show (1) the equilibrium CO 2 saturations in the rock matrix vary from 0.10 to 0.60 at controlled capillary pressures up to 200 kPa, (2) the CO 2 saturations in the matrix increase with water saturation in the horizontal fracture resulting in better capillary continuity for water to drain across the fracture, and (3) the fracture water saturation can be enhanced by the non-uniform fracture capillary pressure and countercurrent flow of CO 2 and water across the fracture-matrix interface. In the case of high fracture water saturation, matrix CO 2 saturation largely depends on matrix anisotropy and heterogeneity. The core-scale experimental results contribute to understand the fracture-matrix interactions and CO 2 storage efficiency in fractured porous media.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2022

On the Stability and Phase Behavior of Titan's Subsurface Liquid Columns

On Titan, liquid hydrocarbon may stay in the subsurface porous permeable crust known as the alkanofer, analogous to water in Earth’s aquifer. In addition to pressure gradient, the subsurface liquid in alkanofers is subject to vertical compositional grading due to the gravity and temperature gradient. The common wisdom is that the liquid would normally stay underground in a stability established by the pressure that increases with depth as observed in aquifers on Earth. However, Titan’s liquids consist of nitrogen and hydrocarbons, mainly methane and ethane, the behavior of which is very sensitive to temperature and pressure. Consequently, the liquid density does not always increase with depth, thus may introduce a reverse density profile that leads to vertical convective instability of the liquid column. If reverse density profiles are present, capillary pressures arising from liquid trapped within small pores in the crust can help with the column stability. The liquid held in the capillaries can seal the space below it thus helping with the stability, unless the overpressure built from underneath becomes larger than the capillary pressure, which causes leakage to allow the liquids to seep upward from the deep. This situation is analogous to hydrocarbon seeps on Earth, where oil and natural gas escape the reservoir and flow slowly through network of cracks to the surface. An algorithm based on an extended Gibbs equation commonly used in petroleum reservoir engineering is employed in this work to produce pressure, density, and compositional profiles for evaluating the stability and phase behavior of Titan’s subsurface fluids.

Sugata Tan↗
NASA NTRS1972

High pressure pulsed capillary viscometry

An analytical and test program was conducted in order to establish the feasibility of a multichamber pulsed-capillary viscometer. The initial design incorporated a piston, ram, and seals which produced measured pulses up to 30,000 psi in the closed chamber system. Pressure pulses from one to ten milliseconds were investigated in a system volume of 1 cuin. Four test fluids: a MIL-L-7808, a 5P4E polyphenyl ether, a MIL-L-23699A, and a synthetic hydrocarbon were examined in the test pressure assembly. The pressure-viscosity coefficient and viscosity delay time were determined for the MIL-L-7808 lubricant tested.

Smith, R. L.↗
DOE OSTI2025

Informing field-scale CO 2 storage simulations with sandbox experiments: The effect of small-scale heterogeneities

Small-scale heterogeneities can significantly affect the fate of the CO 2 plume and trapping during CO 2 migration. We conducted geologic carbon storage field-scale simulations to investigate the impact of small-scale heterogeneities on plume dynamics and trapping performance. Small-scale heterogeneities have been shown to increase the amount of trapped CO 2 during buoyancy-driven flow. The trapped CO 2 saturation is validated by previous sandbox experimental work during buoyancy-driven flow in realistic heterogeneous domains and is implemented through the critical CO 2 saturation parameter (i.e., the first non-zero value in the drainage CO 2 relative permeability curve). Depending on the type and degree of heterogeneity, various critical CO 2 saturation values are exhibited. Furthermore, we investigated the effect of small-scale heterogeneities when multiple capillary pressure models are employed. This study demonstrates that an increase in critical CO 2 saturation reduces the CO 2 plume size and lateral extent, accompanying an increase in residual trapping and a decrease in solubility trapping. Lastly, we show that independent of the capillary pressure model used, an increase in critical saturation leads to similar CO 2 plume dynamics distribution and trapping performance. These results emphasize the importance of quantifying the effect of small-scale heterogeneity as they affect the large-scale behavior of the CO 2 plume.

58 GEOSCIENCES↗
DOE OSTI2020

Link Between CO 2 –Induced Wettability and Pore Architecture Alteration

Changes in pore (throat) size, surface roughness and mineralogy induced by supercritical-CO 2 -water-rock reactions impact petrophysical properties such as porosity, permeability, and especially wettability. Herein, we show that these changes directly impact relative permeability and capillary pressure curves, a fact rarely studied in the literature. In this work, we show that CO 2 contact angle changes emerge after Madison Limestone samples were soaked for 400 hours in CO 2 -enriched brine. Coreflooding results show that the water production rate and cumulative water production increased after the rock was exposed to carbonic acid. Moreover, the mercury capillary pressure decreased in meso- and macro-pores, indicating the increase of size in these pores due to reactions. This compounded wettability and pore network alteration can directly affect CO 2 injectivity, migration and storage capacity. Furthermore, this fundamental insight into CO 2 geological storage processes should aid practitioners to reduce uncertainties in forecasting CO 2 distribution via injection simulation.

58 GEOSCIENCES↗
NASA NTRS1983

Pore size engineering applied to the design of separators for nickel-hydrogen cells and batteries

Pore size engineering in starved alkaline multiplate cells involves adopting techniques to widen the volume tolerance of individual cells. Separators with appropriate pore size distributions and wettability characteristics (capillary pressure considerations) to have wider volume tolerances and an ability to resist dimensional changes in the electrodes were designed. The separators studied for potential use in nickel-hydrogen cells consist of polymeric membranes as well as inorganic microporous mats. In addition to standard measurements, the resistance and distribution of electrolyte as a function of total cell electrolyte content were determined. New composite separators consisting of fibers, particles and/or binders deposited on Zircar cloth were developed in order to engineer the proper capillary pressure characteristics in the separator. These asymmetric separators were prepared from a variety of fibers, particles and binders.

Abbey, K. M.↗
NASA NTRS1983

Pore size engineering applied to the design of separators for nickel-hydrogen cells and batteries

Pore size engineering in starved alkaline multiplate cells involves adopting techniques to widen the volume tolerance of individual cells. Separators with appropriate pore size distributions and wettability characteristics (capillary pressure considerations) to have wider volume tolerances and an ability to resist dimensional changes in the electrodes were designed. The separators studied for potential use in nickel-hydrogen cells consist of polymeric membranes as well as inorganic microporous mats. In addition to standard measurements, the resistance and distribution of electrolyte as a function of total cell electrolyte content were determined. New composite separators consisting of fibers, particles and/or binders deposited on Zircar cloth were developed in order to engineer the proper capillary pressure characteristics in the separator. These asymmetric separators were prepared from a variety of fibers, particles and binders. Previously announced in STAR as N83-24571

Abbey, K. M.↗
DOE OSTI2023

Hydrostratigraphic Region 4 Model - with HSU3-4Pc

Farnsworth Unit (FWU) CO2EOR Eclipse compositional model: Hydrostratigraphic Region 4 model uses the Morrow2 (UNOCAL) relative permeability curve with HSU3-4 capillary pressure assigned homogeneously

Capillary Pressure↗
DOE OSTI2021

Physical and chemical effects of isopropanol exchange in cement-based materials

Drying cement-based materials is a necessary step to characterize the microstructure by microscopy, adsorption, or mercury intrusion. The isopropanol (IPA) exchange method is commonly used to replace the pore solution prior to drying, as it reduces the capillary pressure and thereby helps to preserve the microstructure of the hydrated cementitious materials. However, some physical and chemical effects of IPA on hydration products have been reported. These effects cannot be completely avoided, but can be reduced by shortening the exchange duration, in particular for effects on ettringite. This study carried out experiments with different exchange durations and IPA removal techniques to address research gaps in the literature. For 1-mm cubes of cement paste with IPA diffusion coefficient of 1 × 10{sup −11} m{sup 2}/s, the compared exchange durations are 5 min, 15 min (value from the literature), 40 min (when the mixture at the sample center achieves 95% of the drop in surface tension from water to IPA), 80 min (when the IPA/water mixture at the sample center reaches an azeotrope), 5 h (when over 99% IPA is replaced at the sample center), and 24 h (prolonged exchange). After IPA exchange, samples were either directly dried at 40 °C with flushing N{sub 2} or washed by diethyl ether prior to N{sub 2} drying. Even though the exchange duration of 40 and 80 min, in theory, can reduce the damage from capillary pressure, the nitrogen adsorption results do not show any advantage in preventing microstructural alteration. The 5-h exchange provided results similar to the 24-h exchange, but the latter induced slightly more serious chemical effects resulting from the longer contact of IPA with hydration products. In addition, the compared IPA removal techniques do not show significant differences. Therefore, this study suggests using small samples (~1 mm) for microstructural studies following IPA exchange for a period corresponding to the time when 99% of the water (mole fraction) at the sample center is replaced.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Adhesion-force micro-scale study of desiccating granular material

Experiments on five-, four-, three- and two-wet-hydrophilic-grain clusters were performed to investigate evolution of adhesion of granular media during drying on the micro-scale. The experiments show that the adhesion-force of a cluster initially grows at most to three times the original value before decreasing to zero by the end of evaporation. The adhesion-force is composed of capillary pressure force acting over the liquid/solid contact surface area, and surface tension forces acting over the three-phase contact perimeter length. This is in contrast with most macro-scale phenomenological models, in which the only desaturation process variables affecting strength are suction and saturation. Both the contact surface area and contact perimeter length are reduced to zero upon complete liquid evaporation. The morphology of an evaporating water body evolves through slow flow controlled by evaporation rate, interrupted by various modes of fast air entry, which are non-equilibrium jumps of liquid/gas interfaces (Haines jumps). The instabilities involve large adhesion force discontinuities and substantial water mass reconfiguration with water flow in an extremely short time, which makes the process transient. The reconfigurations can reduce the original multi-grain water clusters to four-, three- and two-grain clusters by way of three different instability modes: of thin-sheet instability, or meniscus snap-through instability, depending on the sign of the Gauss curvature of the liquid surface, or finally, for two-grain bridges only, a liquid wire pinch-off. For larger meso-scale assemblies, however, the global adhesion-force evolution is little affected by the jumps. The air entries are potential sites for drying cracks. The (approximately) calculated capillary pressure for two- and three-grain clusters, in no cases is seen to reach high values, predicted from water retention curves.

Engineering↗
DOE OSTI2023

Hydrostratigraphic Region 4 Model - with HSUPc

Farnsworth Unit (FWU) CO2EOR Eclipse compositional model: Hydrostratigraphic Region 4 model uses the Morrow2 (UNOCAL) relative permeability curve with capillary pressure assigned heterogeneously by hydrostratigraphic unit

Capillary Pressure↗
DOE OSTI2021

Spontaneous imbibition in a square tube with corner films: theoretical model and numerical simulation

Spontaneous imbibition in an angular tube with corner films is a fundamental problem in many scientific and engineering processes. In this paper, a modified interacting capillary bundle model is developed to describe the liquid imbibition dynamics in a square tube with corner films. The square tube is decomposed into several interacting sub-capillaries and the local capillary pressure in each sub-capillary is derived based on the specific shape of its meniscus. The conductance of each sub-capillary is calculated using single-phase lattice Boltzmann simulation. The modified interacting capillary bundle model and color-gradient lattice Boltzmann method are used to simulate the liquid imbibition dynamics in the square tube with different fluid properties. The predictions by the modified interacting capillary bundle model match well with the lattice Boltzmann simulation results for different conditions, demonstrating the accuracy and robustness of the interacting capillary bundle model to describe the imbibition dynamics with corner films. In addition, the interacting capillary bundle model is helpful to investigate the mechanisms during spontaneous imbibition and the influences of fluid viscosity, surface tension, wetting phase contact angle and gravity on imbibition dynamics. Finally, a universal scaling law of imbibition dynamics for the main meniscus is developed and the scaling law for arc meniscus is also analyzed.

42 ENGINEERING↗
DOE OSTI2023

Hydrostratigraphic Region 1 Model

Farnsworth Unit (FWU) CO2EOR Eclipse compositional model: Hydrostratigraphic Region 1 model uses the HS1 relative permeability curves assigned heterogeneously by hydrostratigraphic unit. The model with capillary pressure applies the HS1 curves heterogeneously by hydrostratigraphic unit.

Capillary Pressure↗