Capillary filling with pseudo-potential binary Lattice-Boltzmann model

Color-gradient lattice Boltzmann model for immiscible fluids with density contrast

We present a color-gradient-based lattice Boltzmann model for immiscible fluids with a large density contrast. The model employs the velocity-based equilibrium distribution function, initially proposed for the phase-field-based model by Zu and He [Phys. Rev. E 87, 043301 (2013)1539-375510.1103/PhysRevE.87.043301], with a modification necessary to satisfy the kinematic condition at the interface. Different from the existing color-gradient models, the present model allows to specify interface mobility that is independent of the fluid density ratio. Further, we provide a unified framework, which uses the recursive representation of the lattice Boltzmann equation, to derive the governing equations of the system. The emergent color dynamics thus obtained, through an analysis of the segregation operator, is shown to obey the locally conservative Allen-Cahn equation. We use a series of benchmarks, which include a stationary drop, a layered Poiseuille flow, translation of a drop under a forced velocity field, the Rayleigh-Taylor instability, and the capillary intrusion test to demonstrate the model's ability in dealing with complex flow problems.

Spontaneous imbibition in parallel layers of packed beads

The imbibition of a wetting fluid in a homogeneous porous medium follows the diffusion-like behavior described by Washburn. The impregnation of a two-layered porous medium by a wetting fluid due to capillary action has been previously described to have two fronts, one saturating the medium and the other, leading front, which propagates in finer pores. Here, we report that the leading front is governed by the porous structure and is not always in the finer pores. Based on the experiments in a layered porous medium of permeability varying perpendicular to the direction of flow, we show that the permeability of the adjacent layers plays a significant role in determining the leading front amongst the layers. We have also developed an analytical model which describes the flow dynamics in the layered porous medium. The model predicts the condition for which the leading front in the larger pores is followed by the front in the finer pores. This condition is also verified experimentally.

Capillary Displacement of Viscous Liquids in Geometries with Axial Variations

Axial variations in geometry and presence of viscous displaced fluid are known to alter the diffusive-dynamics of capillary imbibition of a wetting liquid. We here show that the coupled effect of axially varying capillary geometry and finite viscosity of the displaced fluid can lead to significant variations in both short and long time dynamics of imbibition. Based on a theoretical model and lattice Boltzmann simulations, we analyze capillary displacement of a viscous liquid in straight and diverging capillaries. At short times, the imbibition length scales proportionally with time as opposed to the diffusive-dynamics of imbibition of a single wetting liquid. Whereas, at long times, geometry-dependent power-law behavior occurs which qualitatively resembles single liquid imbibition. The distance at which the crossover between these two regimes occurs depends strongly on the viscosities of the imbibing and the displaced liquid. Additionally, our simulations show that the early time imbibition dynamics are also affected by the dynamic contact angle of the meniscus.

Soft matter in hard confinement: phase transition thermodynamics, structure, texture, diffusion and flow in nanoporous media

Spatial confinement in nanoporous media affects the structure, thermodynamics and mobility of molecular soft matter often markedly. This article reviews thermodynamic equilibrium phenomena, such as physisorption, capillary condensation, crystallisation, self-diffusion, and structural phase transitions as well as selected aspects of the emerging field of spatially confined, non-equilibrium physics, i.e. the rheology of liquids, capillarity-driven flow phenomena, and imbibition front broadening in nanoporous materials. The observations in the nanoscale systems are related to the corresponding bulk phenomenologies. The complexity of the confined molecular species is varied from simple building blocks, like noble gas atoms, normal alkanes and alcohols to liquid crystals, polymers, ionic liquids, proteins and water. Mostly, experiments with mesoporous solids of alumina, gold, carbon, silica, and silicon with pore diameters ranging from a few up to 50 nm are presented. The observed peculiarities of nanopore-confined condensed matter are also discussed with regard to applications. A particular emphasis is put on texture formation upon crystallisation in nanoporous media, a topic both of high fundamental interest and of increasing nanotechnological importance, e.g. for the synthesis of organic/inorganic hybrid materials by melt infiltration, the usage of nanoporous solids in crystal nucleation or in template-assisted electrochemical deposition of nano structures.

Capillary filling dynamics of viscoelastic fluids

We consider the filling of a capillary by a viscoelastic fluid described by the Phan-Thien-Tanner (PTT) constitutive behavior. By considering both vertical capillary filling and horizontal capillary filling, we demarcate the role played by gravity and fluid rheology towards long-time oscillations in the capillary penetration depth. We also consider the isothermal filling of the capillary for a closed channel and thus bring out the fundamental differences in the nature of capillary filling for PTT and Newtonian fluids for closed channels in comparison to open channels. Through a scaling analysis, we highlight a distinct viscoelastic regime in the horizontal capillary filling which is in contrast to the Washburn scaling seen in the case of Newtonian fluids. Such an analysis with a very general constitutive behavior is therefore expected to shed light on many areas of microfluidics which focus on biofluids that are often well described by the PTT constitutive behavior.

Gravity-driven thin liquid films over topographical substrates

We investigate the time-dependent evolution of thin liquid films over inclined substrates using a multi-component lattice Boltzmann algorithm. Substrates with and without grooves are considered and the effects of the inclination angle on the dynamics and the coating of the substrates are studied. Our results indicate that the dynamics is enhanced and the ridge height and its displacement are increased by increasing the inclination angle. However, by increasing the inclination angle the maximum depth that can be successfully coated is reduced. Also, although for any given groove depth the width should be larger than a critical value for successful coating, the critical width decreases for smaller inclination angles. For different inclination angles we derive and report the critical sizes of the grooves for successful coating of the substrates.

Effective fluid front of the moving meniscus in capillary

Traditionally, the meniscus bottom is taken as the fluid front when tracking the fluid motion in capillary, but in simulation studies, the thus-calculated motion curve deviates notably from the modified Lucas-Washburn equation. Here, we report that, by considering a volume equivalent of the meniscus part, the motion of the equivalent front agrees very well with the theoretical prediction; furthermore, such an effective fluid front can be directly represented by a specific position of the meniscus, which is independent of the capillary radius. These findings provide an accurate and practical method for describing the motion of the fluid front in capillary.

Thin liquid film flow over substrates with two topographical features

A multicomponent lattice Boltzmann scheme is used to investigate the surface coating of substrates with two topographical features by a gravity-driven thin liquid film. The considered topographies are U- and V-shaped grooves and mounds. For the case of substrates with two grooves, our results indicate that for each of the grooves there is a critical width such that if the groove width is larger than the critical width, the groove can be coated successfully. The critical width of each groove depends on the capillary number, the contact angle, the geometry, and the depth of that groove. The second groove critical width depends on, in addition, the geometry and the depth of the first groove; for two grooves with the same geometries and depths, it is at least equal to that of the first groove. If the second groove width lies between the critical widths, the second groove still can be coated successfully on the condition that the distance between the grooves is considered larger than a critical distance. For considered contact angles and capillary numbers our results indicate that the critical distance is a convex function of the capillary number and the contact angle. Our study also reveals similar results for the case of substrates with a mound and a groove.

Multicomponent interparticle-potential lattice Boltzmann model for fluids with large viscosity ratios

This work focuses on an improved multicomponent interparticle-potential lattice Boltzmann model. The model results in viscosity-independent equilibrium densities and is capable of simulating kinematic viscosity ratios greater than 1000. External forces are incorporated into the discrete Boltzmann equation, rather than through an equilibrium velocity shift as in the original Shan and Chen (hereafter, SC) model. The model also requires the derivation of a momentum conserving effective velocity, which is substituted into the equilibrium distribution function and applies to both the single- and multiple-relaxation-time formulations. Additionally, higher-order isotropy is used in the calculation of the fluid-fluid interaction forces to reduce the magnitude of spurious currents (i.e., numerical errors) in the vicinity of interfaces. First, we compare the model to the SC model for static bubble simulations. We demonstrate that the model results in viscosity-independent equilibrium bubble densities for a wide range of kinematic viscosities, which is not the case for the SC model. Furthermore, we show that the model is capable of simulating stable bubbles for kinematic viscosity ratios greater than 1000 (when higher-order isotropy is used), whereas the SC model is known to be limited to kinematic viscosity ratios on the order of 10. Next we verify the model for surface tension via Laplace's law and show that the model results in the same surface tension values for a range of kinematic viscosities and kinematic viscosity ratios of 10, 100, and 1000. The model is also verified for layered cocurrent flow though parallel plates. We show that the simulated velocity profiles preserve continuity at the interface for kinematic viscosity ratios ranging from 0.001 to 1000 and that the model accurately predicts nonwetting and wetting phase relative permeability for kinematic viscosity ratios of 0.01 to 100.

Modeling the transport of nanoparticle-filled binary fluids through micropores

Understanding the transport of multicomponent fluids through porous medium is of great importance for a number of technological applications, ranging from ink jet printing and the production of textiles to enhanced oil recovery. The process of capillary filling is relatively well understood for a single-component fluid; much less attention, however, has been devoted to investigating capillary filling processes that involve multiphase fluids, and especially nanoparticle-filled fluids. Here, we examine the behavior of binary fluids containing nanoparticles that are driven by capillary forces to fill well-defined pores or microchannels. To carry out these studies, we use a hybrid computational approach that combines the lattice Boltzmann model for binary fluids with a Brownian dynamics model for the nanoparticles. This hybrid approach allows us to capture the interactions among the fluids, nanoparticles, and pore walls. We show that the nanoparticles can dynamically alter the interfacial tension between the two fluids and the contact angle at the pore walls; this, in turn, strongly affects the dynamics of the capillary filling. We demonstrate that by tailoring the wetting properties of the nanoparticles, one can effectively control the filling velocities. Our findings provide fundamental insights into the dynamics of this complex multicomponent system, as well as potential guidelines for a number of technological processes that involve capillary filling with nanoparticles in porous media.

Spontaneous imbibition in disordered porous solids: a theoretical study of helium in silica aerogels

We present a theoretical study of spontaneous imbibition of liquid (4)He in silica aerogels focusing on the effect of porosity on the fluid dynamical behavior. We adopt a coarse-grained three-dimensional lattice-gas description like in previous studies of gas adsorption and capillary condensation and use a dynamical mean-field theory, assuming that capillary disorder predominates over permeability disorder as in recent phase-field models of spontaneous imbibition. Our results reveal a remarkable connection between imbibition and adsorption as also suggested by recent experiments. The imbibition front is always preceded by a precursor film, and the classical Lucas-Washburn √t scaling law is generally recovered, although some deviations may exist at large porosity. Moreover, the interface roughening is modified by wetting and confinement effects. Our results suggest that the interpretation of the recent experiments should be revised.

Capillary rise of water in hydrophilic nanopores

We report on the capillary rise of water in three-dimensional networks of hydrophilic silica pores with 3.5 nm and 5 nm mean radii, respectively (porous Vycor monoliths). We find classical square root of time Lucas-Washburn laws for the imbibition dynamics over the entire capillary rise times of up to 16 h investigated. Provided we assume two preadsorbed strongly bound layers of water molecules resting at the silica walls, which corresponds to a negative velocity slip length of -0.5 nm for water flow in silica nanopores, we can describe the filling process by a retained fluidity and capillarity of water in the pore center. This anticipated partitioning in two dynamic components reflects the structural-thermodynamic partitioning in strongly silica bound water layers and capillary condensed water in the pore center which is documented by sorption isotherm measurements.