Pore-scale effects during the transition from capillary- to viscosity-dominated flow dynamics within microfluidic porous-like domains.

We perform a numerical and experimental study of immiscible two-phase flows within predominantly 2D transparent PDMS microfluidic domains with disordered pillar-like obstacles, that effectively serve as artificial porous structures. Using a high sensitivity pressure sensor at the flow inlet, we capture experimentally the pressure dynamics under fixed flow rate conditions as the fluid-fluid interface advances within the porous domain, while also monitoring the corresponding phase distribution patterns using optical microscopy. Our experimental study covers 4 orders of magnitude with respect to the injection flow rate and highlights the characteristics of immiscible displacement processes during the transition from the capillarity-controlled interface displacement regime at lower flow rates, where the pores are invaded sequentially in the form of Haines jumps, to the viscosity-dominated regime, where multiple pores are invaded simultaneously. In the capillary regime, we recover a clear correlation between the recorded inlet pressure and the pore-throat diameter invaded by the interface that follows the Young-Laplace equation, while during the transition to the viscous regime such a correlation is no longer evident due to multiple pore-throats being invaded simultaneously (but also due to significant viscous pressure drop along the inlet and outlet channels, that effectively mask capillary effects). The performed experimental study serves for the validation of a robust Level-Set model capable of explicitly tracking interfacial dynamics at sub-pore scale resolutions under identical flow conditions. The numerical model is validated against both well-established theoretical flow models, that account for the effects of viscous and capillary forces on interfacial dynamics, and the experimental results obtained using the developed microfluidic setup over a wide range of capillary numbers. Our results show that the proposed numerical model recovers very well the experimentally observed flow dynamics in terms of phase distribution patterns and inlet pressures, but also the effects of viscous flow on the apparent (i.e. dynamic) contact angles in the vicinity of the pore walls. For the first time in the literature, this work clearly shows that the proposed numerical approach has an undoubtable strong potential to simulate multiphase flow in porous domains over a wide range of Capillary numbers.

Drying in porous media with gravity-stabilized fronts: experimental results.

In a recent paper [Yiotis et al., Phys. Rev. E 85, 046308 (2012)] we developed a model for the drying of porous media in the presence of gravity. It incorporated effects of corner film flow, internal and external mass transfer, and the effect of gravity. Analytical results were derived when gravity opposes drying and hence leads to a stable percolation drying front. In this paper, we test the theory using laboratory experiments. A series of isothermal drying experiments in glass bead packings saturated with volatile hydrocarbons is conducted. The transparent glass cells containing the packing allow for the visual monitoring of the phase distribution patterns below the surface, including the formation of liquid films, as the gaseous phase invades the pore space, and for the control of the thickness of the diffusive mass boundary layer over the packing. The experimental results agree very well with theory, provided that the latter is generalized to account for the effects of corner roundness in the film region (which was neglected in the theoretical part). We demonstrate the existence of an early constant rate period (CRP), which lasts as long as the films saturate the surface of the packing, and of a subsequent falling rate period (FRP), which begins practically after the detachment of the film tips from the external surface. During the CRP, the process is controlled by diffusion within the stagnant gaseous phase in the upper part of the cells, yielding a Stefan tube problem solution. During the FRP, the process is controlled by diffusion within the packing, with a drying rate inversely proportional to the observed position of the film tips in the cell. Theoretical and experimental results compare favorably for a specific value of the roundness of the films, which is found to be constant and equal to 0.2 for various conditions, and verify the theoretical dependence on the capillary Ca(f), Bond Bo, and Sherwood Sh numbers.

Analytical solutions of drying in porous media for gravity-stabilized fronts.

We develop a mathematical model for the drying of porous media in the presence of gravity. The model incorporates effects of corner flow through macroscopic liquid films that form in the cavities of pore walls, mass transfer by diffusion in the dry regions of the medium, external mass transfer over the surface, and the effect of gravity. We consider two different cases: when gravity opposes liquid flow in the corner films and leads to a stable percolation drying front, and when it acts in the opposite direction. In this part, we develop analytical results when the problem can be cast as an equivalent continuum and described as a one-dimensional (1D) problem. This is always the case when gravity acts against drying by opposing corner flow, or when it enhances drying by increasing corner film flow but it is sufficiently small. We obtain results for all relevant variables, including drying rates, extent of the macroscopic film region, and the demarkation of the two different regimes of constant rate period and falling rate period, respectively. The effects of dimensionless variables, such as the bond number, the capillary number, and the Sherwood number for external mass transfer are investigated. When gravity acts to enhance drying, a 1D solution is still possible if an appropriately defined Rayleigh number is above a critical threshold. We derive a linear stability analysis of a model problem under this condition that verifies front stability. Further analysis of this problem, when the Rayleigh number is below critical, requires a pore-network simulator which will be the focus of future work.

Lattice Boltzmann method for Lennard-Jones fluids based on the gradient theory of interfaces.

In the present study we propose a lattice Boltzmann equation (LBE) model derived from density gradient expansions of the discrete BBGKY evolution equations. The model is based on the mechanical approach of the gradient theory of interfaces. The basic input is the radial distribution function, which is related exclusively to the molecular interaction potential, rather than semiempirical equations of state used in previous LBE models. This function can be provided from independent molecular simulations or from approximate theories. Evidently the accuracy of the interaction potential, and thus the radial distribution function, reflects on the accuracy of the thermodynamic properties and consistency of the derived LBE model. We have applied the proposed model to obtain equilibrium bulk and interfacial properties of a Lennard-Jones fluid at different temperatures, T, close to critical, T(c). The results demonstrate that the LBE model is in excellent agreement with gradient theory as well as with independent literature results based on different molecular simulation approaches. Hence the proposed LBE model can recover accurately bulk and interfacial thermodynamics for a Lennard Jones fluid at T/T(c)>0.9.

Thermodynamic consistency of liquid-gas lattice Boltzmann methods: interfacial property issues.

In the present study we examine the thermodynamic consistency of lattice Boltzmann equation (LBE) models that are based on the forcing method by comparing different numerical treatments of the LBE for van der Waals fluids. The different models are applied for the calculation of bulk and interfacial thermodynamic properties at various temperatures. The effect of the interface density gradient parameter, kappa , that controls surface tension, is related explicitly with the fluid characteristics, including temperature, molecular diameter, and lattice spacing, through the employment of a proper intermolecular interaction potential. A comprehensive analysis of the interfacial properties reveals some important shortcomings of the LBE methods when central finite difference schemes are employed in the directional derivative calculations and proposes a proper treatment that ensures thermodynamically consistent interfacial properties in accord with the van der Waals theory. The results are found to be in excellent quantitative agreement with exact results of the van der Waals theory preserving all the major features of the interfacial characteristics of vapor-liquid systems of different shapes and sizes.

Effect of liquid films on the isothermal drying of porous media.

We study the effects of liquid films on the isothermal drying of porous media. They are important for the transport of liquid to an evaporation interface, far from the receding liquid clusters. Through a transformation, the drying problem is mapped to the Laplace equation around the percolation liquid clusters. From its solution, the properties of drying are obtained in terms of the capillary number. Consistent with experimental evidence, film flow is shown to accelerate drying significantly.