Evaluation of methods using topology and integral geometry to assess wettability

In situ characterization of heterogeneous surface wetting in porous materials

The performance of nano- and micro-porous materials in capturing and releasing fluids, such as during CO2 geo-storage and water/gas removal in fuel cells and electrolyzers, is determined by their wettability in contact with the solid. However, accurately characterizing wettability is challenging due to spatial variations in dynamic forces, chemical heterogeneity, and surface roughness. In situ measurements can potentially measure wettability locally as a contact angle - the angle a denser phase (e.g water) contacts solid in the presence of a second phase (e.g. hydrogen, air, CO2) - but suffer from difficulties in accurately capturing curvatures, contact areas, and contact loops of multiphase fluids. We introduce a novel extended topological method for in situ contact angle measurement and provide a comparative review of current geometric and topological methods, assessing their accuracy on ideal surfaces, porous rocks containing CO2, and water in gas diffusion layers. The new method demonstrates higher accuracy and reliability of in situ measurements for uniformly wetting systems compared to previous topological approaches, while geometric measurements perform best for mixed-wetting domains. This study further provides a comprehensive open-source platform for in situ characterization of wettability in porous materials with implications for gas geo-storage, fuel cells and electrolyzers, filtration, and catalysis.

Universal description of wetting on multiscale surfaces using integral geometry

HYPOTHESIS

Emerging energy-related technologies deal with multiscale hierarchical structures, intricate surface morphology, non-axisymmetric interfaces, and complex contact lines where wetting is difficult to quantify with classical methods. We hypothesise that a universal description of wetting on multiscale surfaces can be developed by using integral geometry coupled to thermodynamic laws. The proposed approach separates the different hierarchy levels of physical description from the thermodynamic description, allowing for a universal description of wetting on multiscale surfaces.

THEORY AND SIMULATIONS

The theoretical framework is presented followed by application to limiting cases of wetting on multiscale surfaces. Limiting cases include those considered in the Wenzel, Cassie-Baxter, and wicking state models. Wetting characterisation of multiscale surfaces is explored by conducting simulations of a fluid droplet on a structurally rough surface and a chemically heterogeneous surface.

FINDINGS

The underlying origin of the classical wetting models is shown to be rooted within the proposed theoretical framework. Integral geometry provides a topological-based wetting metric that is not contingent on any type of wetting state. The wetting metric is demonstrated to account for multiscale features along the common line in a scale consistent way; providing a universal description of wetting for multiscale surfaces.

Determination of contact angles for three-phase flow in porous media using an energy balance

HYPOTHESIS

We define contact angles, θ, during displacement of three fluid phases in a porous medium using energy balance, extending previous work on two-phase flow. We test if this theory can be applied to quantify the three contact angles and wettability order in pore-scale images of three-phase displacement.

THEORY

For three phases labelled 1, 2 and 3, and solid, s, using conservation of energy ignoring viscous dissipation (Δa_1scosθ₁₂-Δa₁₂-ϕκ₁₂ΔS₁)σ₁₂=(Δa_3scosθ₂₃+Δa₂₃-ϕκ₂₃ΔS₃)σ₂₃+Δa₁₃σ₁₃, where ϕ is the porosity, σ is the interfacial tension, a is the specific interfacial area, S is the saturation, and κ is the fluid-fluid interfacial curvature. Δ represents the change during a displacement. The third contact angle, θ₁₃ can be found using the Bartell-Osterhof relationship. The energy balance is also extended to an arbitrary number of phases.

FINDINGS

X-ray imaging of porous media and the fluids within them, at pore-scale resolution, allows the difference terms in the energy balance equation to be measured. This enables wettability, the contact angles, to be determined for complex displacements, to characterize the behaviour, and for input into pore-scale models. Two synchrotron imaging datasets are used to illustrate the approach, comparing the flow of oil, water and gas in a water-wet and an altered-wettability limestone rock sample. We show that in the water-wet case, as expected, water (phase 1) is the most wetting phase, oil (phase 2) is intermediate wet, while gas (phase 3) is most non-wetting with effective contact angles of θ₁₂≈48° and θ₁₃≈44°, while θ₂₃=0 since oil is always present in spreading layers. In contrast, for the altered-wettability case, oil is most wetting, gas is intermediate-wet, while water is most non-wetting with contact angles of θ₁₂=134°±~10°,θ₁₃=119°±~10°, and θ₂₃=66°±~10°.

Contact Angles in Two-Phase Flow Images

In this work, we calculate contact angles in X-ray tomography images of two-phase flow in order to investigate the wettability. Triangulated surfaces, generated using the images, are smoothed to calculate the contact angles. As expected, the angles have a spread rather than being a constant value. We attempt to shed light on sources of the spread by addressing the overlooked mesh corrections prior to smoothing, poorly resolved image features, cluster-based analysis, and local variations of contact angles. We verify the smoothing algorithm by analytical examples with known contact angle and curvature. According to the analytical cases, point-wise and average contact angles, average mean curvature and surface area converge to the analytical values with increased voxel grid resolution. Analytical examples show that these parameters can reliably be calculated for fluid–fluid surfaces composed of roughly 3000 vertices or more equivalent to 1000 pixel2. In an experimental image, by looking into individual interfaces and clusters, we show that contact angles are underestimated for wetting fluid clusters where the fluid–fluid surface is resolved with less than roughly 500 vertices. However, for the fluid–fluid surfaces with at least a few thousand vertices, the mean and standard deviation of angles converge to similar values. Further investigation of local variations of angles along three-phase lines for large clusters revealed that a source of angle variations is anomalies in the solid surface. However, in the places least influenced by such noise, we observed that angles tend to be larger when the line is convex and smaller when the line is concave. We believe this pattern may indicate the significance of line energy in the free energy of the two-phase flow systems.

Characterization of wetting using topological principles

HYPOTHESIS

Understanding wetting behavior is of great importance for natural systems and technological applications. The traditional concept of contact angle, a purely geometrical measure related to curvature, is often used for characterizing the wetting state of a system. It can be determined from Young's equation by applying equilibrium thermodynamics. However, whether contact angle is a representative measure of wetting for systems with significant complexity is unclear. Herein, we hypothesize that topological principles based on the Gauss-Bonnet theorem could yield a robust measure to characterize wetting.

THEORY AND EXPERIMENTS

We introduce a macroscopic contact angle based on the deficit curvature of the fluid interfaces that are imposed by contacts with other immiscible phases. We perform sessile droplet simulations followed by multiphase experiments for porous sintered glass and Bentheimer sandstone to assess the sensitivity and robustness of the topological approach and compare the results to other traditional approaches.

FINDINGS

We show that the presented topological principle is consistent with thermodynamics under the simplest conditions through a variational analysis. Furthermore, we elucidate that at sufficiently high image resolution the proposed topological approach and local contact angle measurements are comparable. While at lower resolutions, the proposed approach provides more accurate results being robust to resolution-based effects. Overall, the presented concepts open new pathways to characterize the wetting state of complex systems and theoretical developments to study multiphase systems.