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.
A Novel Relative Permeability Model for Gas and Water Flow in Hydrate-Bearing Sediments With Laboratory and Field-Scale Application
In a producing gas hydrate reservoir the effective porosity available for fluid flow constantly changes with dissociation of gas hydrate. Therefore, accurate prediction of relative permeability using legacy models (e.g. Brooks-Corey (B-C), van Genuchten, etc.) that were developed for conventional oil and gas reservoirs would require empirical parameters to be calibrated at various S h over its range of variation, but such calibrations are precluded because of lack of experimental relative permeability data. This study proposes a new relative permeability model for gas hydrate-bearing media that is a function of maximum capillary pressure, capillary entry pressure, pore size distribution index, residual saturations, hydrate saturation, and four other constants. The three novel features of the proposed model are: (i) requires fitting its six empirical parameters only once using experimental data from any single S h , and the same set of empirical parameters predict relative permeability at all S h , (ii) includes the effect of capillarity, and (iii) includes the effect of pore-size distribution. From practical standpoint, the model can be used to simulate multiphase flow in gas hydrate-bearing sediments where the proposed relative permeability can account for the evolving hydrate saturation. The proposed model is implemented in a numerical simulator and the wall time required to perform simulations using the proposed model is shown to be similar to the time it takes to run same simulations with the B-C model. The proposed model is a step forward towards achieving the goal of physically accurate modeling of multiphase flow for gas hydrate-bearing sediments that accounts for the effect of gas hydrate saturation change on relative permeability.
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.
Direct measurement of capillary blood pressure in the human lip
In this study, we developed and tested a new procedure for measuring microcirculatory blood pressures above heart level in humans. Capillary and postcapillary venule blood pressures were measured directly in 13 human subjects by use of the servonulling micropressure technique adapted for micropuncture of lip capillaries. Pressure waveforms were recorded in 40 separate capillary vessels and 14 separate postcapillary venules over periods ranging from 5 to 64 s. Localization and determination of capillary and postcapillary vessels were ascertained anatomically before pressure measurements. Capillary pressure was 33.2 +/- 1.5 (SE) mm Hg in lips of subjects seated upright. Repeated micropunctures of the same vessel gave an average coefficient of variation of 0.072. Postcapillary venule pressure was 18.9 +/- 1.6 mm Hg. This procedure produces a direct and reproducible means of measuring microvascular blood pressures in a vascular bed above heart level in humans.
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%.
Numerical investigation and experimental validation of water condensation in the gas diffusion layer with different properties
Liquid water in porous media hinders the transport of reactant to catalyst layer, where the electrochemical reactions occur, which in turn, reduces performance and lifetime of a proton exchange membrane fuel cell (PEMFC). In this study, two types of gas diffusion layer (GDL) materials are used to study the two-phase saturation distribution in a PEMFC. The correlations between the effective and anisotropic transport properties and porosity of the two GDLs are solved by the pore scale model. In addition, the anisotropic liquid water permeability and the relationship between saturation pressure and capillary pressure is determined by using the Lattice-Boltzmann method. Fuel cell performance as well as liquid water distribution in the GDL using neutron radiography are obtained for validation. Under wet conditions, the cell with Freudenberg GDL performs better than that using Toray GDL, especially at high current densities. The results of the two-phase model simulation show that the peak water saturation for Toray GDL occurs in the catalyst layer and can reach 40%, while the peak water saturation for Freudenberg GDL occurs inside the GDL directly under the ribs and only reaches 25% saturation. The combined results provide key insights to enable high power density operation of a PEMFC through GDL material optimization.
Two-Phase Fluid Flow Properties of Rough Fractures With Heterogeneous Wettability: Analysis With Lattice Boltzmann Simulations
Fractures are conduits for fluid flow in low-permeability geological formations. Multiphase flow properties of fractures are important in natural processes and in engineering applications such as the evaluation of the sealing capacity of caprocks and productivity of hydrocarbon-bearing tight rocks. Investigations of flow and transport through fractures typically focus on the effects of fracture geometric and mechanical factors such as aperture, roughness, and compressibility. The wettability of the fracture surfaces and its influence on microscale interfacial phenomena and macroscale effective transport properties are seldom studied. In this study, we investigated the effect of heterogeneous wetting properties on the displacement of water by supercritical CO2 through a series of lattice Boltzmann method simulations. The results show the evolution of the CO2 plume within a fracture is controlled by both the roughness of the aperture field and the wetting distribution. We combined these factors into a capillary pressure map that can be related to the macroscopic flow behavior of the fracture. We observed that heterogeneous wetting distributions promote the residual trapping of water where lower capillary pressures allowed for isolated water pockets in higher capillary pressure zones. Analysis of fracture unsteady relative permeability shows the effect of wetting on permeability evolution and provides support for the viscous-coupling relative permeability model. Finally, analysis of the steady-state relative permeability and saturation demonstrates a strong correlation between permeability and the standard deviation of the capillary pressure field. Thus, characterizing the distribution of wetting properties of fractures is crucial to understanding multiphase fracture flow and transport properties.
A multi-dimensional parametric study of variability in multi-phase flow dynamics during geologic CO 2 sequestration accelerated with machine learning
Successful geologic CO 2 storage projects depend on numerical simulations to predict reservoir performance during site selection, injection verification, and post-injection monitoring phases of the project. These numerical simulations solve non-linear sets of coupled partial differential equations, while accounting for multi-phase fluid dynamics on the basis of constitutive equations that are embedded into the solution scheme. As a consequence, individual simulations often require tens to hundreds of hours to complete on high-performance computing clusters. Moreover, laboratory experiments reveal that parametric functions for capillary pressure and relative permeability exhibit substantial variability, even within the same rock type. This combination of computational expense and wide-ranging parametric variability means that there remains substantial uncertainty in the behavior of multi-phase CO 2 -water systems, particularly in the context of feedbacks between relative permeability and capillary pressure. To bridge this knowledge gap, here we develop a novel workflow that utilizes physics-based numerical simulation to train an artificial neural network (ANN) emulator for interrogating the multivariate parameter space that governs both capillary pressure and relative permeability. With this approach, the ANN is trained to emulate both fluid pressure distribution and CO 2 saturation, which are then interrogated quantitatively to generate parametric response surface mappings with high-fidelity resolution. Results from this study initially show that capillary entry pressure is the dominant control on both CO 2 plume geometry and fluid pressure propagation when considering the combined effects of capillary pressure and relative permeability, particularly when phase interference is low and residual CO 2 saturation is high. Moreover, the ANN emulator provides tremendous computational speed-up by computing 2691 individual simulations in several minutes; whereas, the same simulation ensemble would have required ~3 years of simulation time using only physics-based simulation methods (25,000 times speed up).
Relationship of echocardiographic indices to pulmonary capillary wedge pressures in healthy volunteers
OBJECTIVES: We sought to determine the relationship between different echocardiographic indices and pulmonary capillary wedge pressures (PCWP) in normal volunteers. BACKGROUND: Indices based on tissue Doppler (TDE) and color M-mode (CMM) echocardiography have been proposed to reflect left (LV) ventricular filling pressures. These include the ratio of early diastolic transmitral velocity (E) to early myocardial velocity measured by TDE (E') and the ratio of E to the wave propagation velocity (Vp) measured from CMM images. These indices, however, have not been validated in normal individuals. METHODS: We studied seven volunteers during two phases of preload altering maneuvers, baseline, with two stages of lower body negative pressure, and repeat baseline with two stages of volume loading. The PCWP obtained from right heart catheterization was compared with diastolic indices using pulsed Doppler, TDE and CMM echocardiography. RESULTS: The PCWP ranged from 2.2 to 23.5 mm Hg. During preload alterations, significant changes in E and septal E' (both p < 0.05) but not lateral E' or Vp were observed. Furthermore, E, septal E' and E/Vp correlated with PCWP (all r > 0.80) but not combined E and TDE indices (both r < 0.15). Within individuals, a similar linear relationship was observed among E/Vp, E and septal E' (average r > 0.80). CONCLUSIONS: In subjects without heart disease, E, septal E' and E/Vp correlate with PCWP. Because the influence of ventricular relaxation is minimized, the ratio E/Vp may be the best overall index of LV filling pressures.
Heterogeneous multiphase flow properties of volcanic rocks and implications for noble gas transport from underground nuclear explosions
Of interest to the Underground Nuclear Explosion Signatures Experiment are patterns and timing of explosion-generated noble gases that reach the land surface. The impact of potentially simultaneous flow of water and gas on noble gas transport in heterogeneous fractured rock is a current scientific knowledge gap. This article presents field and laboratory data to constrain and justify a triple continua conceptual model with multimodal multiphase fluid flow constitutive equations that represents host rock matrix, natural fractures, and induced fractures from past underground nuclear explosions (UNEs) at Aqueduct and Pahute Mesas, Nevada National Security Site, Nevada, USA. Capillary pressure from mercury intrusion and direct air–water measurements on volcanic tuff core samples exhibit extreme spatial heterogeneity (i.e., variation over multiple orders of magnitude). Petrographic observations indicate that heterogeneity derives from multimodal pore structures in ash-flow tuff components and post-depositional alteration processes. Comparisons of pre- and post-UNE samples reveal different pore size distributions that are due in part to microfractures. Capillary pressure relationships require a multimodal van Genuchten (VG) constitutive model to best fit the data. Relative permeability estimations based on unimodal VG fits to capillary pressure can be different from those based on bimodal VG fits, implying the choice of unimodal vs. bimodal fits may greatly affect flow and transport predictions of noble gas signatures. The range in measured capillary pressure and predicted relative permeability curves for a given lithology and between lithologies highlights the need for future modeling to consider spatially distributed properties.
Correlation of capillary entrance pressure drops with normal stress data.
Pressure drops at entrance of capillary tube correlated with normal stress data for viscoelastic polymer solutions
Comprehensive parametric study of CO 2 sequestration in deep saline aquifers
Carbon dioxide injection in deep saline aquifers is a key method for permanently sequestering anthropogenic CO 2 . Here this study employs a reactive transport model to explore mineral precipitation/dissolution and its impact on reservoir properties in deep saline aquifers. We also assess capillary pressure and relative permeability hysteresis on various CO 2 trapping mechanisms. Results from this study reveal the significant influence of initial brine composition on mineral precipitation/dissolution. The dissolution and precipitation of minerals have different effects around the wellbores compared to the overall reservoir. Additionally, salt concentration (Ca ++ and Mg ++ ) and quartz surface area affect CO 2 mineralization, while Na + impacts halite precipitation, altering flow properties. The effect of capillary pressure is significant, as including the capillary pressure in the simulation case resulted in significantly improved CO 2 trapping, achieving almost total dissolution of the injected CO 2 in around 300 years. This study offers novel insights into the interactions of reservoir minerals, brine properties, and the injected CO 2 .
Constrained Vapor Bubble
The nonisothermal Constrained Vapor Bubble, CVB, is being studied to enhance the understanding of passive systems controlled by interfacial phenomena. The study is multifaceted: 1) it is a basic scientific study in interfacial phenomena, fluid physics and thermodynamics; 2) it is a basic study in thermal transport; and 3) it is a study of a heat exchanger. The research is synergistic in that CVB research requires a microgravity environment and the space program needs thermal control systems like the CVB. Ground based studies are being done as a precursor to flight experiment. The results demonstrate that experimental techniques for the direct measurement of the fundamental operating parameters (temperature, pressure, and interfacial curvature fields) have been developed. Fluid flow and change-of-phase heat transfer are a function of the temperature field and the vapor bubble shape, which can be measured using an Image Analyzing Interferometer. The CVB for a microgravity environment, has various thin film regions that are of both basic and applied interest. Generically, a CVB is formed by underfilling an evacuated enclosure with a liquid. Classification depends on shape and Bond number. The specific CVB discussed herein was formed in a fused silica cell with inside dimensions of 3x3x40 mm and, therefore, can be viewed as a large version of a micro heat pipe. Since the dimensions are relatively large for a passive system, most of the liquid flow occurs under a small capillary pressure difference. Therefore, we can classify the discussed system as a low capillary pressure system. The studies discussed herein were done in a 1-g environment (Bond Number = 3.6) to obtain experience to design a microgravity experiment for a future NASA flight where low capillary pressure systems should prove more useful. The flight experiment is tentatively scheduled for the year 2000. The SCR was passed on September 16, 1997. The RDR is tentatively scheduled for October, 1998.
Transcapillary fluid shifts in head and neck tissues during and after simulated microgravity
To understand the mechanism, magnitude, and time course of facial puffiness that occurs in microgravity, seven male subjects were tilted 6 degrees head down for 8 hr, and all four Starling transcapillary pressures were directly measured before, during, and after tilt. Head-down tilt (HDT) caused facial edema and a significant elevation of microvascular pressures measured in the lower lip: capillary pressures increased from 27.2 +/- 5 mm Hg pre-HDT to 33.9 +/- 1.7 mm Hg by the end of tilt. Subcutaneous and intramuscular interstitial fluid pressures in the neck also increased as a result of HDT, while interstitial fluid colloid osmotic pressures remained unchanged. Plasma colloid osmotic pressures dropped significantly after 4 hr of HDT, suggesting a transition from fluid filtration to absorption in capillary beds between the heart and feet during HDT. After 4 hr of seated recovery from HDT, microvascular pressures remained significantly elevated by 5 to 8 mm Hg above baseline values, despite a significant HDT diuresis and the orthostatic challenge of an upright, seated posture. During the control (baseline) period, urine output was 46.7 ml/hr; during HDT, it was 126.5 ml/hr. These results indicate that facial edema resulting from HDT is primarily caused by elevated capillary pressures and decreased plasma colloid osmotic pressures. Elevation of cephalic capillary pressures sustained for 4 hr after HDT suggests that there is a compensatory vasodilation to maintain microvascular perfusion. The negativity of interstitial fluid pressures above heart level also has implications for the maintenance of tissue fluid balance in upright posture.
Controls on reservoir properties in organic-matter-rich shales: Insights from MICP analysis
Mercury injection capillary pressure (MICP) analysis was used to study organic-matter-rich Devonian shales varying in maturity from marginally mature to postmature to discuss controls on drainage and imbibition in shale reservoirs. Six samples come from the New Albany Shale; their total organic carbon (TOC) ranged from 3.3–8.4 % and vitrinite reflectance (R o ) from 0.54–1.42 %. One Marcellus Shale sample with R o 2.50 % was used to extend the maturity range to the advanced dry gas window. Three coal samples of varying rank were also included for comparison. Porosity of the shales ranges from 1.5–5.0 %, and MICP-derived permeability ranges from 0.17–9.9 mD. Injection capillary pressure curves have a range of entry pressures from 0.055 to0.421 MPa, and no relationship of entry pressure with maturity was revealed. Maturity, however, has an influence on pore size distribution, with the higher-maturity samples having significantly lower volumes of pores with throats less than 10 nm in diameter. Withdrawal capillary pressure curves document a wide range of withdrawal efficiencies from 55 % to almost 100 %, and pores having throats between 3–50 nm in diameter play a major role in the ability of mercury to be released from the pore system. The obtained data provide valuable information about shale properties and their behavior not only in the air-mercury system but also under reservoir conditions.
Simplified Green-Ampt Model, Imbibition-Based Estimates of Permeability, and Implications for Leak-off in Hydraulic Fracturing
Predicting water imbibition into porous materials is important in a wide variety of fields, and hydraulic fracturing of low permeability hydrocarbon reservoirs has emerged as an application that is imposing a large water footprint. Furthermore, reliable predictions of imbibition are needed to better manage water use, yet are challenging because of uncertainties in both the permeability and capillary pressure driving force. Here, this uncertainty is reduced through evaluating correlations between the permeability and the effective capillary pressure associated with the wetting front, P c,f . These correlations allow elimination of P c,f from the Green and Ampt equation and concentrate all uncertainties in fluxes on the effective permeability k . Over a wide range of k and porosities n , imbibition scales approximately with k 1/3 . Although Leverett k 1/4 scaling for predicting P c,f is shown to be inferior when tested with data spanning a wide range of n , it nevertheless predicted imbibition fairly well. From simple imbibition measurements, both the empirical and Leverett scaling approaches allow estimates of k that have root-mean-square deviations of about one order of magnitude relative to measurements that ranged over 10 orders of magnitude in k .
Using Direct Numerical Simulation of Pore-Level Events to Improve Pore-Network Models for Prediction of Residual Trapping of CO2
Direct numerical simulation and pore-network modeling are common approaches to study the physics of two-phase flow through natural rocks. For assessment of the long-term performance of geological sequestration of CO 2 , it is important to model the full drainage-imbibition cycle to provide an accurate estimate of the trapped CO 2 . While direct numerical simulation using pore geometry from micro-CT rock images accurately models two-phase flow physics, it is computationally prohibitive for large rock volumes. On the other hand, pore-network modeling on networks extracted from micro-CT rock images is computationally efficient but utilizes simplified physics in idealized geometric pore elements. This study uses the lattice-Boltzmann method for direct numerical simulation of CO 2 -brine flow in idealized pore elements to develop a new set of pore-level flow models for the pore-body filling and snap-off events in pore-network modeling of imbibition. Lattice-Boltzmann simulations are conducted on typical idealized pore-network configurations, and the interface evolution and local capillary pressure are evaluated to develop modified equations of local threshold capillary pressure of pore elements as a function of shape factor and other geometrical parameters. The modified equations are then incorporated into a quasi-static pore-network flow solver. The modified model is applied on extracted pore-network of sandstone samples, and saturation of residual trapped CO 2 is computed for a drainage-imbibition cycle. The modified model yields different statistics of pore-level events compared with the original model; in particular, the occurrence of snap-off in pore-throats is reduced resulting in a more frontal displacement pattern along the main injection direction. Compared to the original model, the modified model is in closer agreement with the residual trapped CO 2 obtained from core flow experiments and direct numerical simulation.