Contact Angle Distribution Research Articles

Spatial Correlation of Contact Angle and Curvature in Pore‐Space Images

AbstractWe study the in situ distributions of contact angle and oil/brine interface curvature measured within millimeter‐sized rock samples from a producing hydrocarbon carbonate reservoir imaged after waterflooding using X‐ray microtomography. We analyze their spatial correlation combining automated methods for measuring contact angles and interfacial curvature (AlRatrout et al., 2017, https://doi.org/10.1016/j.advwatres.2017.07.018), with a recently developed method for pore‐network extraction (Raeini et al., 2017, https://doi.org/10.1103/PhysRevE.96.013312). The automated methods allow us to study image volumes of diameter approximately 1.9 mm and 1.2 mm long, obtaining hundreds of thousands of values from a data set with 435 million voxels. We calculate the capillary pressure based on the mode oil/brine interface curvature value and associate this value with a nearby throat in the pore space. We demonstrate the capability of our methods to distinguish different wettability states in the samples studied: water‐wet, weakly oil‐wet, and mixed‐wet. The contact angle and oil/brine interface curvature are spatially correlated over approximately the scale of an average pore. There is a wide distribution of contact angles within a single pore. A range of local oil/brine interface curvature is found with both positive and negative values. There is a correlation between interfacial curvature and contact angle in trapped ganglia, with ganglia in water‐wet patches tending to have a positive curvature and oil‐wet regions seeing negative curvature. We observed a weak correlation between average contact angle and pore size, with the larger pores tending to be more oil‐wet.

Wettability in complex porous materials, the mixed-wet state, and its relationship to surface roughness

A quantitative in situ characterization of the impact of surface roughness on wettability in porous media is currently lacking. We use reservoir condition micrometer-resolution X-ray tomography combined with automated methods for the measurement of contact angle, interfacial curvature, and surface roughness to examine fluid/fluid and fluid/solid interfaces inside a porous material. We study oil and water in the pore space of limestone from a giant producing oilfield, acquiring millions of measurements of curvature and contact angle on three millimeter-sized samples. We identify a distinct wetting state with a broad distribution of contact angle at the submillimeter scale with a mix of water-wet and water-repellent regions. Importantly, this state allows both fluid phases to flow simultaneously over a wide range of saturation. We establish that, in media that are largely water wet, the interfacial curvature does not depend on solid surface roughness, quantified as the local deviation from a plane. However, where there has been a significant wettability alteration, rougher surfaces are associated with lower contact angles and higher interfacial curvature. The variation of both contact angle and interfacial curvature increases with the local degree of roughness. We hypothesize that this mixed wettability may also be seen in biological systems to facilitate the simultaneous flow of water and gases; furthermore, wettability-altering agents could be used in both geological systems and material science to design a mixed-wetting state with optimal process performance.

Quasiperiodic granular chains and Hofstadter butterflies.

We study quasiperiodicity-induced localization of waves in strongly precompressed granular chains. We propose three different set-ups, inspired by the Aubry-André (AA) model, of quasiperiodic chains; and we use these models to compare the effects of on-site and off-site quasiperiodicity in nonlinear lattices. When there is purely on-site quasiperiodicity, which we implement in two different ways, we show for a chain of spherical particles that there is a localization transition (as in the original AA model). However, we observe no localization transition in a chain of cylindrical particles in which we incorporate quasiperiodicity in the distribution of contact angles between adjacent cylinders by making the angle periodicity incommensurate with that of the chain. For each of our three models, we compute the Hofstadter spectrum and the associated Minkowski-Bouligand fractal dimension, and we demonstrate that the fractal dimension decreases as one approaches the localization transition (when it exists). We also show, using the chain of cylinders as an example, how to recover the Hofstadter spectrum from the system dynamics. Finally, in a suite of numerical computations, we demonstrate localization and also that there exist regimes of ballistic, superdiffusive, diffusive and subdiffusive transport. Our models provide a flexible set of systems to study quasiperiodicity-induced analogues of Anderson phenomena in granular chains that one can tune controllably from weakly to strongly nonlinear regimes.This article is part of the theme issue 'Nonlinear energy transfer in dynamical and acoustical systems'.

Direct simulation of multiphase flows with modeling of dynamic interface contact angle

We describe a modeling technique for dynamic contact angle between a phase interface and a solid wall using a generalized Navier boundary condition in the context of a front-tracking-based multiphase method. The contact line motion is determined by the generalized Navier slip boundary condition in order to eliminate the infinite shear stress at the contact line. Applying this slip boundary condition only to the interface movement with various slip ratios shows good agreement with experimental results compared to allowing full fluid slip along the solid surface. The interface slip model performs well on grid convergence tests using both the slip ratio and slip length models. A detailed energy analysis was performed to identify changes in kinetic, surface, and potential energies as well as viscous and contact line dissipation with time. A friction coefficient for contact line dissipation was obtained based on the other computed energy terms. Each energy term and the friction coefficient were compared for different grid resolutions. The effect of varying the slip ratio as well as the contact angle distribution versus contact line speed was analyzed. The behavior of drop impact on a solid wall with different advancing and receding angles was investigated. Finally, the proposed dynamic contact model was extended to three dimensions for large-scale parallel calculations. The impact of a droplet on a solid cylinder was simulated to demonstrate the capabilities of the proposing formulation on general solid structures. Widely different contact angles were tested and showed distinctive characteristic behavior clearly.

Image analysis of axisymmetric droplets in wetting experiments: A new tool for the study of 3D droplet geometry and droplet shape reconstruction

A new software tool is developed for droplet image analysis for the study of wetting of solid substrates. The tool extracts information exclusively from images and does not require the use of any properties of the system. Moreover, its applicability covers both axisymmetric and non-axisymmetric droplets. The developed software processes independently droplet images taken from side and top perspectives. Processing of side-view images is made by polynomial fitting to the droplet shape and provides important 2D geometrical features such as contact angles, length, height, contact point coordinates and contour outline. The analysis of top-view images is achieved through active contours (snakes) and yields droplet dimensions, centroid position and contour outline. The combination of synchronized side and top images provides the reconstruction of the droplet 3D shape by means of slices of circular arc shape, which allows estimation of droplet volume and of the distribution of contact angles along its perimeter. The above is an important feature that has not been delivered by other software tools. Experimental results to support the applicability of the new tool are presented for two distinct substrates having different surface properties (glass and Teflon).

Droplet departure modeling and a heat transfer correlation for dropwise flow condensation in hydrophobic mini-channels

Droplet nucleation, growth, coalescence, and departure control dropwise condensation heat transfer. Smaller droplets are associated with higher heat transfer coefficients due to their lower liquid thermal resistances. Unlike quiescent dropwise condensation with gravity-driven droplet departure, droplet departure sizes in flow condensation are governed by flow-droplet shear forces and droplet-solid adhesive forces. This research models droplet departure, droplet size distributions, and heat transfer through single droplets under different flow conditions. Heat transfer through single droplets includes the thermal resistances at the vapor-liquid interface, temperature depression across the curved surface, conduction in the liquid droplet, and conduction through the surface promoter (e.g., Teflon). Droplet size distributions were determined for two ranges using the population balance method and power law function for small and large droplets, respectively. Droplet departure sizes (e.g., 10–500 µm) were derived using force balances between drag forces (obtained using FLUENT) and droplet-solid adhesive forces (determined using a third-order polynomial for contact angle distribution along contact line). The analytical model was compared to experimental flow condensation heat transfer data in a Teflon AF™-coated rectangular mini-gap with hydraulic diameters of 0.95 and 1.8 mm. The correlation was compared against experiments with a steam mass flux range of 35–75 kg/m2 s and quality of 0.2–0.9. There was good agreement between the model and experimental data; without any curving fitting, the mean absolute errors of the heat transfer correlation were 9.6% and 8.8% respectively for the 0.95-mm and 1.8-mm mini-gaps.

Analysis of varying contact angles and load distributions in defective angular contact ball bearing

Abstract Angular contact ball bearing (ACBB) is widely used in rotating machinery for its great external load capacity. The operating condition of ACBB is easily affected by surrounding mechanical components. The variations of contact angles and load distributions can sufficiently reflect the operating condition of ACBB. A localized defect on the raceway will lead to the dramatic change of contact angle and load distribution of ACBB. In this paper, a mechanical model of ACBB with a localized defect on outer raceway is established, in which the effects of centrifugal force and gyroscopic moment on the force equilibrium relation of balls in high speed operating condition are considered. The proposed model is verified by comparing with the existing model in literature. On the basis of the proposed bearing model, the effects of some parameters on the variations of contact angles and load distributions are investigated. The results show that as the operating speed increases, the centrifugal forces of balls increase and the balls are gradually moving towards outer raceway, which results in that the load distributions between balls and outer raceway are much larger than those between balls and inner raceway. In view of the combined external forces and moments acted on the bearing, the variations of contact angles and load distributions become more complicated than those under the single force or moment applied condition. With the increase of the circumferential extent of localized defect, abrupt changes of contact angles and load distributions in defect area increase obviously.

Drop pattern resulting from the breakup of a bidimensional grid of liquid filaments

A rectangular grid formed by liquid filaments on a partially wetting substrate evolves in a series of breakups leading to arrays of drops with different shapes distributed in a rather regular bidimensional pattern. Our study is focused on the configuration produced when two long parallel filaments of silicone oil, which are placed upon a glass substrate previously coated with a fluorinated solution, are crossed perpendicularly by another pair of long parallel filaments. A remarkable feature of this kind of grids is that there are two qualitatively different types of drops. While one set is formed at the crossing points, the rest are consequence of the breakup of shorter filaments formed between the crossings. Here, we analyze the main geometric features of all types of drops, such as shape of the footprint and contact angle distribution along the drop periphery. The formation of a series of short filaments with similar geometric and physical properties allows us to have simultaneously quasi identical experiments to study the subsequent breakups. We develop a simple hydrodynamic model to predict the number of drops that results from a filament of given initial length and width. This model is able to yield the length intervals corresponding to a small number of drops, and its predictions are successfully compared with the experimental data as well as with numerical simulations of the full Navier–Stokes equation that provide a detailed time evolution of the dewetting motion of the filament till the breakup into drops. Finally, the prediction for finite filaments is contrasted with the existing theories for infinite ones.

In situ characterization of mixed-wettability in a\xa0reservoir rock at subsurface conditions

We used X-ray micro-tomography to image the in situ wettability, the distribution of contact angles, at the pore scale in calcite cores from a producing hydrocarbon reservoir at subsurface conditions. The contact angle was measured at hundreds of thousands of points for three samples after twenty pore volumes of brine flooding.We found a wide range of contact angles with values both above and below 90°. The hypothesized cause of wettability alteration by an adsorbed organic layer on surfaces contacted by crude oil after primary drainage was observed with Scanning Electron Microscopy (SEM) and identified using Energy Dispersive X-ray (EDX) analysis. However, not all oil-filled pores were altered towards oil-wet conditions, which suggests that water in surface roughness, or in adjacent micro-porosity, can protect the surface from a strong wettability alteration. The lowest oil recovery was observed for the most oil-wet sample, where the oil remained connected in thin sheet-like layers in the narrower regions of the pore space. The highest recovery was seen for the sample with an average contact angle close to 90°, with an intermediate recovery in a more water-wet system, where the oil was trapped in ganglia in the larger regions of the pore space.

A Pore Network Modelling Study of Fuel Cell Gas Diffusion Layers with Patterned Wettability

Further improvements in power density of polymer electrolyte fuel cells (PEFCs) are required to lower the cost. At high current densities, mass transport related resistances limit the performance due to, mainly, a poor water management within the cell. It is therefore paramount to develop advanced water management strategies to maximize power output. In this context, diffusion layers play a notable role. They must ensure efficient gas transport toward the catalyst layer while removing the electrochemically produced liquid water in the opposite flow direction. A substantial amount of work has been carried out to improve gas diffusion layers (GDLs) and microporous layer over the last years, mainly controlling microstructural features, such as pore sizes, fiber arrangements and porosities, and hydrophobic coating load and distribution [1]. Our group has recently developed an approach consisting of artificially creating hydrophilic patterns throughout the entire GDL thickness using the radiation grafting method. The hydrophilic areas act as low pressure liquid water removal pathways, therefore liberating the remaining areas for the gas to flow throughout less tortuous domains [2]. The performance of PEFCs was significantly increased thanks to the use of the modified GDLs [3]. However, there is still plenty of room for improvement. Capillary pressure is the driving force for liquid water transport and it is influenced by the microstructure and liquid-solid interaction [4]. It therefore stands to reason that, for a given application, further optimization can be achieved when selecting the appropriate diffusion layer microstructure and wettability ratio. While the first experimental data investigating the effect of GDL morphology, contact angle and coating load has been recently published, these experiments require notable infrastructure and time efforts [5]. For that reason, modelling capillary pressure experiments using pore network models can be a valid approach to theoretically assess optimal material parameters. The open-source OpenPNM framework, which has been extensively used in recent fuel cell literature [6], was the software of choice for this investigation. In this talk we will start by presenting the model physics in which a modification of the Washburn equation has been used as required in the intermediate wettability range where this model fails to predict accurately water imbibition and the validation strategy. The definition of coating load and contact angle will be discussed. A good agreement was obtained between model and experiments, specially investigating the influence of coating load, contact angle and pore size distributions. In the second part of the talk, we will provide some material design guidelines based on model output information to maximize the quality of water segregation in GDLs with patterned wettability.

Contact angle entropy and macroscopic friction in noncohesive two-dimensional granular packings.

We study the relationship between the granular contact angle distribution and local particle friction on the macroscopic friction and bulk modulus in noncohesive disk packings. Molecular dynamics in two dimensions are used to simulate uniaxial loading-unloading cycles imposed on the granular packings. While macroscopic Mohr friction depends on the granular pack geometric details, it reaches a stationary limit after a finite number of loading-unloading cycles that render well-defined values for bulk modulus, grain coordination, porosity, and friction. For random packings and for all polydispersities analyzed, we found that as interparticle friction increases, the bulk modulus for the limit cycle decreases linearly, while the mean coordination number is reduced and the porosity increased, also as approximately linear functions. On the other hand, the macroscopic Mohr friction increases in a monotonous trend with interparticle friction. The latter result is compared to a theoretical model that assumes the existence of sliding planes corresponding to definite Mohr-friction values. The simulation results for macroscopic friction are well described by the theoretical model that incorporates the local neighbor angle distribution that can be quantified through the contact angle entropy. As local friction is increased, the limit entropy of the neighbor angle distribution is reduced, thus introducing the geometric component to granular friction. Surprisingly, once the limit cycle is reached, the Mohr friction seems to be insensitive to polydispersity as has been recently reported.

Pore‐scale water dynamics during drying and the impacts of structure and surface wettability

AbstractPlants and microbes secrete mucilage into soil during dry conditions, which can alter soil structure and increase contact angle. Structured soils exhibit a broad pore size distribution with many small and many large pores, and strong capillary forces in narrow pores can retain moisture in soil aggregates. Meanwhile, contact angle determines the water repellency of soils, which can result in suppressed evaporation rates. Although they are often studied independently, both structure and contact angle influence water movement, distribution, and retention in soils. Here drying experiments were conducted using soil micromodels patterned to emulate different aggregation states of a sandy loam soil. Micromodels were treated to exhibit contact angles representative of those in bulk soil (8.4° ± 1.9°) and the rhizosphere (65° ± 9.2°). Drying was simulated using a lattice Boltzmann single‐component, multiphase model. In our experiments, micromodels with higher contact angle surfaces took 4 times longer to completely dry versus micromodels with lower contact angle surfaces. Microstructure influenced drying rate as a function of saturation and controlled the spatial distribution of moisture within micromodels. Lattice Boltzmann simulations accurately predicted pore‐scale moisture retention patterns within micromodels with different structures and contact angles.

X-ray tomography investigations of mono-sized sphere packing structures in cylindrical containers

The structure of mono-sized sphere packings (diameter d) in cylindrical containers (diameter D and height H) both with and without inner cylinders (diameter Di) has been investigated in detail by means of advanced X-ray computed tomography. The geometrical parameters were varied in a wide range; in all experiments 1d vertical vibration was applied. Five experiments were selected with characteristically differing local packing structures. The influence of container geometry, filling and vibration procedures on the formation of regular packings is discussed and a simple correlation is presented to assess whether structured packings occupy a significant fraction of the total packed volume.For a packing with moderate densification, the regular structures are restricted to small wall zones and a random packing exists in the largest part of the packing volume. By selecting appropriate vibration parameters, the zones with regular structures can increase considerably and can persist in the total packed volume. The increasing crystallisation causes an increase of the container packing fraction. For cylinders with H/D ≫1 and moderate D/d, regular structures develop preferentially in radial direction from a hexagonal layer at the concave wall. For H/D<1 and D/d≫1, hexagonal dense structures grow preferentially above the flat bottom plate and can occupy a great portion of the total volume. The role of granular convection on these crystallisation processes has been addressed. Previous statements that the thickness of wall zones is ≈(4–5)d are not generally valid for mono-sized sphere packings; the development of a comprehensive correlation is the task of a future work.Structural details of the packings close to concave, plane and convex walls are analysed via void fraction distributions, sphere centre positions, contact angle distributions, coordination numbers, radial distribution function and Voronoi tessellation. The combination of these methods provides a comprehensive understanding of structural details. Only a few characteristic results are presented; special topics will be the subject of forthcoming publications.

In Situ Local Contact Angle Measurement in a CO2-Brine-Sand System Using Microfocused X-ray CT.

The wettability of porous media is of major interest in a broad range of natural and engineering applications. The wettability of a fluid on a solid surface is usually evaluated by the contact angle between them. While in situ local contact angle measurements are complicated by the topology of porous media, which can make it difficult to use traditional methods, recent advances in microfocused X-ray computed tomography (micro-CT) and image processing techniques have made it possible to measure contact angles on the scale of the pore sizes in such media. However, the effects of ionic strength, CO2 phase, and flow pattern (drainage or imbibition) on pore-scale contact angle distribution are still not clear and have not been reported in detail in previous studies. In this study, we employed a micro-CT scanner for in situ investigation of local contact angles in a CO2-brine-sand system under various conditions. The effects of ionic strength, CO2 phase, and flow pattern on the local contact-angle distribution were examined in detail. The results showed that the local contact angles vary over a wide range as a result of the interaction of surface contaminants, roughness, pore topology, and capillarity. The wettability of a porous surface could thus slowly weaken with increasing ionic strength, and the average contact angle could significantly increase when gaseous CO2 (gCO2) turns into supercritical CO2 (scCO2). Contact angle hysteresis also occurred between drainage and imbibition procedures, and the hysteresis was more significant under gCO2 condition.

Properties of force networks in jammed granular media

Discrete element method simulations are conducted to investigate the effects of applied uniaxial compressive load, and bidisperse particle size distributions on force networks within jammed granular media. The differences between the strong and weak networks are examined through investigating the spatial correlation and distribution of contact angles, and emergence of chainlike structures. The simulation results show that the chainlike structures are more prevalent in the strong network due to the larger cumulative probabilities of contact angles, but not all the contacts belonging to the strong or weak networks are able to constitute the chainlike structures. Although the contacts of coarse-fine particles are dominant for the bidisperse systems, the contacts of coarse–coarse particles dominate the strong network, as well as the linear chainlike structures. Upon increasing the pressure from very low to high, the probability of contact orientations with respect to the compression direction in the strong network increases for contact orientation less than \(60^{\circ }\) and decreases for contact orientation greater than \(60^{\circ }\), while the opposite trends are observed in the weak network. The tails of normalized normal contact forces distributions are quantified by \(\hbox {P}(\hbox {f}) = \hbox {exp}(-\hbox {cf}^{\mathrm{n}})\), and it is found that the value of n depends on the applied pressure and particle size distribution. Statistical analysis shows that the degree of homogeneity of contact force increases with increasing pressure, which is also validated by participation number.

Improving Nanofiber Membrane Characteristics and Membrane Distillation Performance of Heat-Pressed Membranes via Annealing Post-Treatment

Electrospun membranes are gaining interest for use in membrane distillation (MD) due to their high porosity and interconnected pore structure; however, they are still susceptible to wetting during MD operation because of their relatively low liquid entry pressure (LEP). In this study, post-treatment had been applied to improve the LEP, as well as its permeation and salt rejection efficiency. The post-treatment included two continuous procedures: heat-pressing and annealing. In this study, annealing was applied on the membranes that had been heat-pressed. It was found that annealing improved the MD performance as the average flux reached 35 L/m2·h or LMH (&gt;10% improvement of the ones without annealing) while still maintaining 99.99% salt rejection. Further tests on LEP, contact angle, and pore size distribution explain the improvement due to annealing well. Fourier transform infrared spectroscopy and X-ray diffraction analyses of the membranes showed that there was an increase in the crystallinity of the polyvinylidene fluoride-co-hexafluoropropylene (PVDF-HFP) membrane; also, peaks indicating the α phase of polyvinylidene fluoride (PVDF) became noticeable after annealing, indicating some β and amorphous states of polymer were converted into the α phase. The changes were favorable for membrane distillation as the non-polar α phase of PVDF reduces the dipolar attraction force between the membrane and water molecules, and the increase in crystallinity would result in higher thermal stability. The present results indicate the positive effect of the heat-press followed by an annealing post-treatment on the membrane characteristics and MD performance.

A Mixed Wettability Pore Size Distribution Based Mathematical Model for Analyzing Two-Phase Flow in Porous Electrodes

The non-isothermal, two-phase membrane electrodes assembly numerical model previously developed in Part I [J. Electrochem. Soc., 164, 6, F530 (2017)] is validated by comparing to experimentally measured electrochemical performance data under various operating conditions. Water accumulation in catalyst layer (CL) and gas diffusion layer (GDL) are also compared to the neutron and X-ray imaging data and shown to be in agreement. Water fluxes at cathode and anode boundaries and phase change induced flow are analyzed and compared to experimental data. Simulation results indicate that when liquid water is present, reactant transport in the catalyst layer is the key factor limiting fuel cell performance. The model shows that at high relative humidity, 80°C is the optimal operating temperature in order to delay water accumulation without degrading performance due to oxygen dilution by water vapor. The impact of CL and GDL hydrophilic volume fraction, hydrophobic contact angle and pore size distribution on performance are also studied. Results suggest that when liquid water is present, GDL parameters have minimal effect on performance and a CL should have a large hydrophobic contact angle, a low hydrophilic volume fraction, and large pore radius.

In-situ characterization of wettability and pore-scale displacements during two- and three-phase flow in natural porous media

We establish a unique approach to measure in-situ contact angle from micro-CT images acquired during two- and three-phase miniature core-flooding experiments in order to overcome the uncertainties associated with conventional contact angle measurement techniques. The measurements are used to quantify the wettability behavior of the rock and explain pore-level displacement events occurring in three-phase flow. Six two-phase experiments are performed on individual core samples with three pairs of fluids, i.e., oil-brine, gas-oil, and gas-brine, and under two thermodynamic conditions: (a) binary-equilibrated, when only the two respective phases are at equilibrium and (b) ternary-equilibrated, when all three phases are equilibrated and only the two desired fluids are injected into the core. A three-phase experiment set is also performed under ternary-equilibrated conditions, which includes gas injection, a waterflood, and an oilflood process. All experiments are performed on Berea miniature core samples using a nonspreading brine-oil-gas fluid system.We measure receding and advancing contact angles at arc menisci and main terminal menisci for the two-phase binary-equilibrated experiments and characterize contact angle hysteresis for each fluid pair. Contact angle hysteresis values are almost identical for all fluid pairs. The results of the two-phase binary- and ternary-equilibrated experiments show similar contact angle distributions for each fluid pair. Contact angle distributions during the three-phase flow experiment are analyzed to develop new insights into relevant complex displacement mechanisms. The results indicate that, during gas injection, the majority of displacements involving oil and water are oil-to-water events. It is observed that, during the waterflood, both oil-to-gas and gas-to-oil displacement events take place. However, the relative frequency of the former is greater. For the oilflood, gas-water interfaces only slightly hinge in pore elements. Pore-scale fluid occupancy maps and the Bartell–Osterhoff constraint verify the above-mentioned findings.

A note on the distribution of stationary contact angles on organic coatings

Static contact angle to organic coating surfaces is treated as a strictly localized parameter. Contact angle topographies (CAT) are generated for eight organic coatings, typically applied as top coats to maritime and offshore structures. The distribution of the CAT is found to follow a Rosin–Rammler–Sperling–Bennett (RRSB) distribution. Distribution parameters are estimated and compared to conventional static contact angles. It is suggested to consider RRSB distribution parameters to characterize the wetting behavior of organic top coats.

Pore‐scale contact angle measurements of CO 2 –brine–glass beads system using micro‐focused X‐ray computed tomography

Conventional methods for measuring contact angle are usually applied on smooth surfaces. Methods concerning contact angle determinations performed directly on pore surfaces of porous media have rarely been reported. This work approaches the pore-scale measurement of local contact angle in a CO2–brine–glass beads system, using micro-focused X-ray computed tomography (micro-CT). Both drainage and imbibition experiments with 0.1 ml/min injection rate were conducted at 40°C and 8 MPa. The effectiveness of this pore-scale approach is confirmed by comparing the results with the results gathered from traditional sessile drop methodology. Observations indicate that the contact angle hysteresis phenomenon was not so obvious for intermediate-wet glass beads in the employed experimental setting. In real reservoir circumstances, the roughness and capillary variation caused significant deviations in contact angle distribution for both drainage and imbibition, even in rock cores consisting of a single-phase material.