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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Quantifying the Effect of Pore‐Size Dependent Wettability on Relative Permeabilities Using Lattice Boltzmann Simulation

Relative permeability is a crucial two-phase property in porous media that can be significantly impacted by wettability conditions. While traditional research has predominantly examined homogeneous wettability, this work explores the less studied pore-size dependent (PSD) wettability, featured by a pore-size dependent wettability distribution. Leveraging high-fidelity Lattice Boltzmann simulations on CT-scanned porous samples, we demonstrate how PSD wettability would impact relative permeability at the pore scale. Our findings reveal that the deviation of relative permeability from the homogeneous wettabilitity induced by PSD wettability can be 5%–20%. The deviation of relative permeability curves increases as the spanning range of the contact angle increases. We also find that this impact is less pronounced as the capillary number increases. By adopting a pore-size-dependent contact angle relationship, our approach provides a more accurate and nuanced understanding of how PSD wettability would impact two-phase flow. These findings discovered at the pore scale may also provide valuable insights on relative permeability at the core to reservoir scales.

25 ENERGY STORAGE↗
DOE OSTI2020

Correlating inter-particle forces and particle shape to shear-induced aggregation/fragmentation and rheology for dilute anisotropic particle suspensions: A complementary study via capillary rheometry and in-situ small and ultra-small angle X-ray scattering

Predicting the stability of suspensions composed of anisotropic particles is difficult due to the complexity and interplay of hydrodynamic and colloidal forces. For example, it is not well understood how non-spherical particle shapes affect suspension-rheological responses. This study presents evidence showing how the stacking of boehmite platelets dispersed in aqueous solution leads to moderate shear thinning, even at low boehmite mass fractions. Wide angle X-ray, small angle X-ray, and ultra-small angle X-ray scattering were used to simultaneously monitor changes to the size and fractal dimensions of boehmite aggregates from 6 – 10,000 Å as the sample was recirculated through an in situ capillary rheometer. The latter also provided simultaneous measurements of the viscosity of the suspensions. Computational fluid dynamics modeling of the apparatus was used to provide a more rigorous analysis of the fluid flow. In addition to shear-induced aggregation/fragmentation correlated to a complicated balance between hydrodynamic and colloidal forces, the results indicate that orientationally-dependent interactions, arising from non-spherical boehmite aggregate sub-units, may explain the relatively large experimental viscosity when the hydrodynamic force is small compared to colloidal forces. Ultimately, the results from this study can provide a foundation for future work investigating particle shape effects on colloidal and hydrodynamics forces.

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

Scaling of two-phase water-steam relative permeability and thermal fluxes in porous media

Two-phase water-steam flow conditions are frequently encountered in many engineering applications, including geothermal reservoirs. Although routine calculations are based on the multiphase Darcy’s law, the role of the topology of the flowing phases at the pore-scale is usually neglected in the estimation of relative permeabilities. Instead, the latter are frequently computed using empirical models like the Corey correlation. In this work, we first apply the model for relative permeabilities based on pore-scale flow regimes developed by Picchi and Battiato (2019), Relative permeability scaling from pore-scale flow regimes, Water Resour. Res. 55, 3215–3233, to scenarios typical of geothermal reservoirs and then extend it by deriving the scaling laws for the transmissibilities and the thermal properties as a function of temperature. First, we discuss the scaling behavior of normalized relative permeabilities in terms of viscosity ratio and capillary number of water-steam systems and, then, we provide a validation of the model against experimental data available in the literature. The model captures well the data trends collected in real 3D porous media. These results suggest that water-steam relative permeabilities follow the same scaling behavior of gas-liquid systems where the non-wetting phase is much less viscous than the wetting phase. Finally, we investigate the impact that relative permeabilities have on heat transfer rates at two-phase flow conditions and the scaling of mass and energy transmissibility and thermal properties of the mixture. An estimation of the exergy carried by the two-phase water-steam mixture is also included.

42 ENGINEERING↗
DOE OSTI2025

Unsteady- and Steady-State Relative Permeability Study with X-ray and Acoustic Monitoring for CO 2 Storage in Deep Saline Aquifers

In this study, we link the multiphase flow measurements with controls of sedimentary structures (e.g., heterogeneity and anisotropy) on relative permeability to variations in ultrasonic velocities for two deep saline aquifer formations (Entrada and Bluff Sandstones) in the San Juan Basin of the Southwestern USA. The rock specimens were extracted from outcrop sites near Durango, CO, USA. They have distinct differences in grain size, cementation composition, and individual chemical amounts, despite both formations being eolian sandstones. We performed a series of unsteady- and steady-state CO 2 -brine relative permeability experiments under capillary-controlled displacement rates. Unsteady-state experiments were conducted at 71 °C and 9.65 MPa; steady-state experiments were conducted at 85 °C and 22.8 MPa and 89 °C and 24.1 MPa for the Bluff and Entrada Sandstones, respectively. During the unsteady experiments, X-ray computed tomography was used to visualize multiphase flow in porous media and quantify saturations during brine drainage under various flow rates. Scan images and saturation profiles indicate that the CO 2 distribution in the pore volume was strongly impacted by the presence of high-angle cross-laminations, heterogeneous rock structure, and direction of bedding orientation. Those factors contribute to dramatic and quick initial breakthroughs and affect the overall saturation dynamics. Steady-state relative permeability tests were conducted at net flow rates of 1 mL/min for both brine drainage and imbibition. During the experimental steps, the CO 2 fractional flow was increased and decreased for both drainage and imbibition scenarios to mimic the front when CO 2 contacts brine and behind the front when brine enters space previously occupied by CO 2 . It was found that compressional velocity decreased, while shear waves slightly increased as brine saturation decreased. The hysteresis effects for the relative permeability and acoustic velocities were distinct. Furthermore, the CO 2 /brine front stability is quantified by applying a mobility ratio approach to spot saturations at which the boundary line between fluids becomes uneven. The results presented in this work can potentially boost the quality and precision of forecasts for the CO 2 storage projects in which the vertical and horizontal core-scale heterogeneity and anisotropy impact the plume migration within host reservoirs.

carbon dioxide (CO2)↗
DOE OSTI2020

Onset dynamics of air-water menisci on rock fracture surfaces

Information on menisci dynamics and equilibrium contact angles, θ e , is needed for modeling multiphase flow of fluids in geologic formations. The wettability of fracture faces is relevant to applications such as waterflooding in enhanced oil recovery and evaluation of caprock integrity for geologic carbon sequestration. We investigated the vertical capillary rise of air-water menisci on exposed fracture faces for a sandstone, a shale, and three granites. The experimental design approximated that employed in the Wilhelmy plate method. Replicate fracture faces were prepared by fracturing cylindrical rock cores using the Brazilian method and splitting the fractured cores apart along their fracture planes. Meniscus onset was imaged from a cross-sectional perspective using dynamic neutron radiography, and quantified with change point analysis. An existing model for meniscus onset fitted the experimental data very well (median R 2 = 0.96). Capillary rise was typified by t 2/3 scaling (where t is time), followed by a constant equilibrium height, z e . The average time taken for the menisci to reach z e was 0.39 s, with no differences between rock types. The menisci achieved a greater vertical extent on the fracture faces of the sedimentary rocks as compared to those of the granites. Apparent θ e values were calculated from the best fit estimates of z e using the Wilhelmy plate equation for a smooth and planar vertical surface. Despite the fractal roughness of some of the fracture faces, there was a significant linear relationship (R 2 = 0.80) between the apparent θ e values and intrinsic θ e values measured on polished surfaces of the same rocks using the sessile drop method. The regression equation indicated the apparent θ e values were consistently greater than the intrinsic values. This overestimation may be due to averaging of pixels at the air-water-solid interface, a required step in the image analysis. Further exploration of the range of applicability and measurement accuracy of this proposed new technique would be valuable.

58 GEOSCIENCES↗
DOE OSTI2022

Multiscale petrophysical characterization and flow unit classification of the Minnelusa eolian sandstones

Integration of petrophysical and geological information is critical to simulation of subsurface carbon storage (GCS). In this sense, two depositional facies were identified from the core description and well-log interpretation, namely massive (MS) and cross-bedded (CB) facies groups. Additionally, pore-scale characteristics were studied by a combination of techniques, e.g. Nuclear Magnetic Resonance (NMR) and mercury intrusion capillary pressure (MICP). Scanning electron microscope (SEM) and petrographic analyses show that the pore structure is dominantly controlled by the depositional environment and dolomite cementation. NMR-T 2 distributions of MS and CB facies show triple and quadruple modes, respectively. In addition, MICP of high- and low-permeability MS facies samples, and their CB facies group mixtures were collected. The MS sample pore-throat size distribution is uni-modal, while the triple-modal characteristic of the mixtures indicates heterogeneous pore structures at the sub-core scale for CB facies. The reliably estimates of porosity and permeability for both facies groups via NMR techniques and the MLR (Multiple Linear Regression) approach demonstrate the applicability of these techniques to eolian sandstone. Moreover, irreducible water saturation via the T 2 -cutoff method correlates strongly with T 2LM instead of porosity. Finally, the rock quality index and flow zone indicator were calculated based on Combinable Magnetic Resonance (CMR) log interpretations. This provides direct connection to properties measured in the well. Four flow units were classified for both facies groups. Finally, results show that better reservoir quality with significant heterogeneities is observed in the CB facies. This study highlights the importance integrating multiscale petrophysical properties including facies, pore architecture and diagenesis analysis with core- to log-scale property characterization. The results herein validate our reservoir characterization and flow unit classification in eolian reservoirs.

58 GEOSCIENCES↗
DOE OSTI2022

Preparation of illite coated geomaterial microfluidic surfaces: Effect of salinity and heat treatment

Unique characteristics of clay (e.g., high surface area, deficiency in positive charges) make any clay abundant shales a good source of adsorbed gas. However, the presence of clay in porous media can also affect shale rock's producibility. The clay-water interaction can cause significant swelling and fines migration in the formation and impair the overall hydrocarbon recovery from shale. Geomaterial micromodels, developed by functionalizing traditional glass or PDMS surfaces with geomaterials (e.g., calcite, Quartz, clay), can represent the physicochemical properties of the natural porous media and have been used during the last few years to understand and evaluate the solid-fluid physicochemical interactions. Here, this study focuses on developing a clay-coated geomaterial surface and investigates the effect of base fluid's salinity on clay adsorption to the glass surface. Glass capillary tubes and straight channel borosilicate glass micromodels are coated with Illite clay to represent the pore-scale clay chemistry of Caney Shale, a Mississippian unconventional play in Southern Oklahoma, USA. 10 wt% of Illite clay slurries made with brines of four different salinities are used to coat the glass-capillary tubes to understand the effect of salinity in clay adsorption on the glass surface and overall coating quality. To achieve a stable coating and evaluate the impact of heat treatment on coating's stability, straight channel glass flow cells coated with Illite clay are heat-treated at low (25 °C) and high temperature (125 °C), and the results are compared with an untreated coated surface. Flooding tests were carried out with brines of different salinities on the coated surface to evaluate the stability of the coating. The experimental results indicated a strong relationship between the brine's salinity and the adsorption of clay particles on the glass surface. An increase in brine concentration resulted in improved adsorption of clay particles on the glass surface. Experiments involving heat treating the glass surface following the coating demonstrated significant improvement in the stability of the coating.

02 PETROLEUM↗
DOE OSTI2021

PREPARATION OF CLAY-COATED GLASS SURFACES FOR MICROFLUIDICS-BASED STUDY ON CLAY-FLUID INTERACTIONS

Unique characteristics of clay, even though make any clay abundant shales a good source of adsorbed gas, make them equally difficult to produce. The clay-water interaction can cause significant swelling and fines migration in the shale formation, both of which can significantly impair the overall hydrocarbon recovery from shale. Geomaterial micromodels, developed by functionalizing traditional glass or PDMS surface with geomaterials (e.g., calcite, quartz, clay) can represent the physico-chemical properties of the natural porous media and have been used during the last few years to understand and evaluate the solid-fluid physicochemical interactions. This study focuses on developing a clay-coated geomaterial surface and investigates the effect of base fluid's salinity on clay adsorption to the glass surface. Glass capillary tubes and straight channel borosilicate glass micromodels are coated with Illite clay minerals to represent the pore-scale clay chemistry of Caney Shale, a Mississippian unconventional play in Southern Oklahoma, USA. 10 wt.% of Illite clay slurries made with brines of four different salinities are used to coat the glass-capillary tubes to understand the effect of salinity in clay adsorption on the glass surface and overall coating quality. To achieve a stable coating and evaluate the impact of heat treatment on coating's stability, straight channel glass flow cells coated with Illite clay are heat-treated at low (25 °C) and high temperature (125 °C). Dynamic flooding tests were carried out with brines on the coated surface to evaluate the stability of the coating. The experiments' outcome indicated a strong relationship between the brine's salinity and the adsorption of clay particles on the glass surface. An increase in brine concentration resulted in improved adsorption of clay particles on the glass surface. Experiments involving heat treating the glass surface following the coating demonstrated significant improvement in the stability of the coating.

02 PETROLEUM↗
DOE OSTI2024

Hierarchical Composites Patterned via 3D Printed Cellular Fluidics

Additive manufacturing of freeform structures containing multiple materials with deterministic spatial arrangement and interactions remains a challenge for most 3D printing processes, due to complex fabrication tool requirements and limitations in printability of some material classes. Here, in this paper, a versatile method is reported to produce architected composites using the concept of cellular fluidics, in which lattices of unit cells are used as templating scaffolds to guide flowable infill materials in a programmed spatial pattern, upon which they are cured in place to produce a deterministically ordered multimaterial solid. The lattice design relies on the unit cell size, type, strut diameter, surface wetting, and distribution of cellular structures to control liquid flow and retention. Individual unit cells are tuned to achieve reliable infilling and combined into higher-order architectures to achieve multiscale composite materials with disparate mechanical properties, including those considered non-printable. Lattice design considerations for leveraging capillary phenomena and demonstrate several methods of patterning polymers in 3D-printed cellular fluidic structures are presented. The concept of tuning the compressive response of an architected composite using a flexible-elastomer as the lattice and a stiff-epoxy as the infill material is illustrated.

36 MATERIALS SCIENCE↗
DOE OSTI2025

CO2 hydrate crystal thickening, morphology, and Raman spectroscopy in a microfluidic device

Gas hydrates are a solid, crystalline form of water that often form at low temperatures and high pressures. Carbon dioxide (CO2) hydrates may form during carbon dioxide capture and storage (CCS) processes. These solid compounds may form in CO2 pipelines, potentially leading to a full blockage and process shutdown for plug removal. On the other hand, formation of CO2 hydrates may be desired for CO2 capture and separation. In either case, understanding the growth behavior and nature of the hydrates is vital to managing these CCS processes. Using a high-pressure, transparent microfluidic reactor, the crystalline film thickening of CO2 hydrates was observed and measured through visual microscopy and Raman spectroscopy. The impact of subcooling, pressure, and CO2 flow rate was investigated, and only CO2 flow rate was found to have a significant impact on the overall thickness of the film. Visual observations and Raman spectroscopy measurements confirmed that two distinct hydrate layers formed during thickening, one which was more porous than the other. The capillary-like channels in the porous layer indicated a mechanism for mass transfer of water through the hydrate layer. A model was developed based on this observation, and it was fit to the thickening data in order to obtain mass transfer coefficients. Results of this study can be applied to CO2 hydrate formation in pipelines and near porous media used for CO2 capture.

Wadsworth, Lindsey [Colorado School of Mines, Gold↗
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 OSTI2020

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.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Millisecond timescale reactions observed via X-ray spectroscopy in a 3D microfabricated fused silica mixer

Determination of electronic structures during chemical reactions remains challenging in studies which involve reactions in the millisecond timescale, toxic chemicals, and/or anaerobic conditions. Here, a three-dimensionally (3D) microfabricated microfluidic mixer platform that is compatible with time-resolved X-ray absorption and emission spectroscopy (XAS and XES, respectively) is presented. This platform, to initiate reactions and study their progression, mixes a high flow rate (0.50–1.5 ml min –1 ) sheath stream with a low-flow-rate (5–90 µl min –1 ) sample stream within a monolithic fused silica chip. The chip geometry enables hydro­dynamic focusing of the sample stream in 3D and sample widths as small as 5 µm. The chip is also connected to a polyimide capillary downstream to enable sample stream deceleration, expansion, and X-ray detection. In this capillary, sample widths of 50 µm are demonstrated. Further, convection–diffusion-reaction models of the mixer are presented. The models are experimentally validated using confocal epifluorescence microscopy and XAS/XES measurements of a ferricyanide and ascorbic acid reaction. The models additionally enable prediction of the residence time and residence time uncertainty of reactive species as well as mixing times. Residence times (from initiation of mixing to the point of X-ray detection) during sample stream expansion as small as 2.1 ± 0.3 ms are also demonstrated. Importantly, an exploration of the mixer operational space reveals a theoretical minimum mixing time of 0.91 ms. The proposed platform is applicable to the determination of the electronic structure of conventionally inaccessible reaction intermediates.

3D microfabrication↗
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

INVESTIGATION OF ROCK-FLUID INTERACTIONS USING GEOMATERIAL MICROFLUIDICS

Illite and Illite-Smectite clays are reportedly the Caney Shale's most prevalent clay minerals. The swelling and fines migration of clays may occur when water-based fluids contact the clay minerals. This may result in the reduction of pore space and permeability, and wettability alteration. Therefore, we quantified the swelling potentials, fines migration and wettability alteration of Illite and Illite-Smectite clays upon exposure to model brines and a produced water, and a model oil and a crude oil. We measured density, viscosity, pH, TDS, and interfacial tension as they are important in multiphase flow and displacement characteristics. Geomaterial microfluidics is recently being used to study physicochemical interactions of solid-fluid systems. It facilitates visualizing the surface behavior upon exposure to various fluids. This study includes preparation and characterization of Illite-Smectite clay-coated glass capillary tubes. Illite and Illite-Smectite clay-coated geomaterial microfluidic surfaces were prepared to investigate the effects of first contact fluid, brine salinity, and aging on wettability. In this work, advancing and receding contact angles of model brine-air, model oil-air, produced water-air, crude oil-air, model brine-model oil systems, and produced water-crude oil systems have been measured in the untreated, and clay-coated microfluidic channels of 1000 µm width and 50 µm depth. The aging effect was studied by retaining the non-aqueous and aqueous solutions in the corresponding microfluidic channels for approximately 24 h in the case of channels first contacted by aqueous and non-aqueous solutions, respectively. The advancing and receding contact angles of Caney shale-crude oil-produced water systems were also measured. The experimental findings indicate the swelling potential, fines migration, and wettability alteration of Illite and Illite-Smectite clays. The results signify the influence of first contact fluid, brine salinity, and aging on wettability of the untreated and clay-coated microfluidic surfaces. The experiments also demonstrate the wetting nature of the Caney shale samples from reservoir and non-reservoir zones.

02 PETROLEUM↗
DOE OSTI2020

Calculation of Hydraulic Conductivity for a Porous Medium with Staggered Array of Solid Grains - 20087

For a porous medium with periodic array of solid grains. the flow of an incompressible viscous fluid in the pore space is considered with an aim of calculating the permeability of the medium. The solid grains are of ellipsoidal shape. For this purpose, a unit-cell boundary-value problem (Stokes problem) is solved numerically. With the assumption of periodicity of the medium structure and the fluid properties (in fact constant), the flow field and pressure distributions are determined numerically with inertial effects ignored. By applying the theory of homogenization to the viscous flow through pore space, a microscale boundary-value problem is defined. After solving the problem, the permeability is determined by taking the micro-cell average of the velocity components, Various porous structures are chosen to allow a wide range of porosity values. The permeability increases mildly for porosity values in the lower range. But it increases sharply for porosity values in the upper range of the porosity. The permeability of the present study agrees with the Kozeny-Carman relation excellently for smaller porosity. Certain discrepancy between the calculated permeability in this study and the Kozeny-Carman relation appears and increases with porosity for larger porosity. It is believed that the discrepancy is due to the replacement of the pores by a bundle of capillaries, which is unrealistic, in the Kozeny-Carman relation. (authors)

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

Impact of organic solutes on capillary phenomena in water-CO2-quartz systems

Hypothesis: The migration of supercritical CO2 (scCO2) injected into underground reservoirs as part of carbon capture and storage is influenced by organic contamination affecting mineral wettability. Molecular dynamics (MD) simulations of relevant systems that incorporate representative organic solutes allow detailed investigation of changes in fundamental interfacial and capillary properties. Experiments: We use MD simulations to explore the effects of four organic solutes (quinoline, decanoic acid, coronene, sorgoleone) on the wettability of quartz by water in the presence of scCO2. We examine the impacts of polar, alkyl, and aromatic moieties as well as fluid flow velocity at elevated temperatures and pressures. Findings: Organic molecules accumulate at the water-CO2 interface, where they distribute according to their size and functional groups. Certain organics penetrate the adsorbed water film at the quartz-CO2 interface, revealing two modes of hydrogen bonding between polar organic functional group, water, and quartz surface –OH groups. Interfacial energies and contact angles are minimally impacted by organic adsorption at the water-CO2 interface, possibly due to simultaneous CO2 desorption. Non- equilibrium MD simulations reveal that flow-induced redistribution of organic compounds modulates the radii of curvature of the advancing and receding water-CO2 interfaces. Our results indicate that organic adsorption on water surfaces is likely ubiquitous during multi-phase flow in soils and sedimen- tary rocks, with implications for the mobilization and transport of organic compounds.

36 MATERIALS SCIENCE↗