Predictive Model of Supercooled Water Droplet Pinning/Repulsion Impacting a Superhydrophobic Surface: The Role of the Gas-Liquid Interface Temperature.

Dynamical analysis of an impacting liquid drop on superhydrophobic surfaces is mostly carried out by evaluating the droplet contact time and maximum spreading diameter. In this study, we present a general transient model of the droplet spreading diameter developed from the previously defined mass-spring model for bouncing drops. The effect of viscosity was also considered in the model by definition of a dash-pot term extracted from experiments on various viscous liquid droplets on a superhydrophobic surface. Furthermore, the resultant shear force of the stagnation air flow was also considered with the help of the classical Homann flow approach. It was clearly shown that the proposed model predicts the maximum spreading diameter and droplet contact time very well. On the other hand, where stagnation air flow is present in contradiction to the theoretical model, the droplet contact time was reduced as a function of both droplet Weber numbers and incoming air velocities. Indeed, the reduction in the droplet contact time (e.g., 35% at a droplet Weber number of up to 140) was justified by the presence of a formed thin air layer underneath the impacting drop on the superhydrophobic surface (i.e., full slip condition). Finally, the droplet wetting model was also further developed to account for low temperature through the incorporation of classical nucleation theory. Homogeneous ice nucleation was integrated into the model through the concept of the reduction of the supercooled water drop surface tension as a function of the gas-liquid interface temperature, which was directly correlated with the Nusselt number of incoming air flow. It was shown that the experimental results was qualitatively predicted by the proposed model under all supercooling conditions (i.e., from -10 to -30 °C).

Pore-Scale Modeling of the Effect of Wettability on Two-Phase Flow Properties for Newtonian and Non-Newtonian Fluids.

The Darcy-scale properties of reservoir rocks, such as capillary pressure and relative permeability, are controlled by multiphase flow properties at the pore scale. In the present paper, we implement a volume of fluid (VOF) method coupled with a physically based dynamic contact angle to perform pore-scale simulation of two-phase flow within a porous medium. The numerical model is based on the resolution of the Navier-Stokes equations as well as a phase fraction equation incorporating a dynamic contact angle model with wetting hysteresis effect. After the model is validated for a single phase, a two-phase flow simulation is performed on both a Newtonian and a non-Newtonian fluid; the latter consists of a polymer solution displaying a shear-thinning power law viscosity. To investigate the effects of contact angle hysteresis and the non-Newtonian nature of the fluid, simulations of both drainage and imbibition are carried out in order to analyze water and oil saturation-particularly critical parameters such as initial water saturation (Swi) and residual oil saturation (Sor) are assessed in terms of wettability. Additionally, the model sensitivities to the consistency factor (χ), the flow behavior index (n), and the advancing and receding contact angles are tested. Interestingly, the model correctly retrieves the variation in Sor and wettability and predicts behavior over a wide range of contact angles that are difficult to probe experimentally.

The Effect of Heat Transfer and Polymer Concentration on Non-Newtonian Fluid from Pore-Scale Simulation of Rock X-ray Micro-CT.

Most of the pore-scale imaging and simulations of non-Newtonian fluid are based on the simplifying geometry of network modeling and overlook the fluid rheology and heat transfer. In the present paper, we developed a non-isothermal and non-Newtonian numerical model of the flow properties at pore-scale by simulation of the 3D micro-CT images using a Finite Volume Method (FVM). The numerical model is based on the resolution of the momentum and energy conservation equations. Owing to an adaptive mesh generation technique and appropriate boundary conditions, rock permeability and mobility are accurately computed. A temperature and concentration-dependent power-law viscosity model in line with the experimental measurement of the fluid rheology is adopted. The model is first applied at isothermal condition to 2 benchmark samples, namely Fontainebleau sandstone and Grosmont carbonate, and is found to be in good agreement with the Lattice Boltzmann method (LBM). Finally, at non-isothermal conditions, an effective mobility is introduced that enables to perform a numerical sensitivity study to fluid rheology, heat transfer, and operating conditions. While the mobility seems to evolve linearly with polymer concentration in agreement with a derived theoretical model, the effect of the temperature seems negligible by comparison. However, a sharp contrast is found between carbonate and sandstone under the effect of a constant temperature gradient. Besides concerning the flow index and consistency factor, a master curve is derived when normalizing the mobility for both the carbonate and the sandstone.