Mesoscopic modeling of a two-phase flow in the presence of boundaries: The contact angle

Stick-slip sliding of water drops on chemically heterogeneous surfaces

We present a comprehensive study of water drops sliding down chemically heterogeneous surfaces formed by a periodic pattern of alternating hydrophobic and hydrophilic stripes. Drops are found to undergo a stick-slip motion whose average speed is an order of magnitude smaller than that measured on a homogeneous surface having the same static contact angle. This motion is the result of the periodic deformations of the drop interface when crossing the stripes. Numerical simulations confirm this view and are used to elucidate the principles underlying the experimental observations.

Interaction pressure tensor for a class of multicomponent lattice Boltzmann models

We present a theory to obtain the pressure tensor for a class of nonideal multicomponent lattice Boltzmann models, thus extending the theory presented by X. Shan [Phys. Rev. E 77, 066702 (2008)] for single-component fluids. We obtain the correct form of the pressure tensor directly on the lattice and the resulting equilibrium properties are shown to agree very well with those measured from numerical simulations. Results are compared with those of alternative theories.

Axisymmetric multiphase lattice Boltzmann method

A lattice Boltzmann method for axisymmetric multiphase flows is presented and validated. The method is capable of accurately modeling flows with variable density. We develop the classic Shan-Chen multiphase model [Phys. Rev. E 47, 1815 (1993)] for axisymmetric flows. The model can be used to efficiently simulate single and multiphase flows. The convergence to the axisymmetric Navier-Stokes equations is demonstrated analytically by means of a Chapmann-Enskog expansion and numerically through several test cases. In particular, the model is benchmarked for its accuracy in reproducing the dynamics of the oscillations of an axially symmetric droplet and on the capillary breakup of a viscous liquid thread. Very good quantitative agreement between the numerical solutions and the analytical results is observed.

Effects of the hierarchical structure of rough solid surfaces on the wetting of microdroplets

We used the lattice Boltzmann method to investigate how the hierarchical structure of a rough solid surface, which in this work is modeled as the microstructure (micropillars) covered with nanostructures (nanopillars), affects the contact angle of microdroplets atop of the solid surface and the wetting transition between the Wenzel and Cassie states. Our simulation results show that the Wenzel-to-Cassie state transition can be achieved by decreasing the fluid-solid attraction, increasing the micropillar spacing, or coating the microstructures with nanostructures. For the effect of the hierarchical structure on the contact angle, we find that the micropillars show a negligible effect on the contact angle, but they may affect the sliding angle. In contrast, it is the nanostructure that determines the contact angle. The contact angle increases with the nanopillar length until reaching a maximal value, but its dependence on the nanopillar spacing becomes more complicated. The contact angle may first increase with the nanopillar spacing and then decreases, or decreases monotonously, depending on whether the liquid enters the nanostructure or not. In this work, we also demonstrate in the presence of contact line pinning, that the pinning effect affects the apparent contact angle.

Flow past superhydrophobic surfaces with cosine variation in local slip length

Anisotropic superhydrophobic surfaces have the potential to greatly reduce drag and enhance mixing phenomena in microfluidic devices. Recent work has focused mostly on cases of superhydrophobic stripes. Here, we analyze a relevant situation of cosine variation of the local slip length. We derive approximate formulas for maximal (longitudinal) and minimal (transverse) directional effective slip lengths that are in good agreement with the exact numerical solution and lattice-Boltzmann simulations. Compared to the case of superhydrophobic stripes, the cosine texture can provide a very large effective slip. However, the difference between eigenvalues of the slip-length tensor is smaller, indicating that the flow is less anisotropic.

Scheme for contact angle and its hysteresis in a multiphase lattice Boltzmann method

In this paper, a scheme for specifying contact angle and its hysteresis is incorporated into a multiphase lattice Boltzmann method. The scheme is validated through investigations of the dynamic behaviors of a droplet sliding along two kinds of walls: a smooth (ideal) wall and a rough or chemically inhomogeneous (nonideal) wall. For an ideal wall, the wettability of solid substrates is able to be prescribed. For a nonideal wall, arbitrary contact angle hysteresis can be obtained through adjusting advancing and receding angles. Significantly different phenomena can be recovered for the two kinds of walls. For instance, a droplet on an inclined ideal wall under gravity is impossible to stay stationary. However, the droplet on a nonideal wall may be pinned due to contact angle hysteresis. The steady interface shapes of the droplet on an inclined nonideal wall under gravity or in a shear flow quantitatively agree well with the previous numerical studies. Besides, the complex motion of a droplet creeping like an inchworm could be simulated. The scheme is found suitable for the study of contact line problems with and without contact angle hysteresis.

Forcing scheme in pseudopotential lattice Boltzmann model for multiphase flows

The pseudopotential lattice Boltzmann (LB) model is a widely used multiphase model in the LB community. In this model, an interaction force, which is usually implemented via a forcing scheme, is employed to mimic the molecular interactions that cause phase segregation. The forcing scheme is therefore expected to play an important role in the pseudoepotential LB model. In this paper, we aim to address some key issues about forcing schemes in the pseudopotential LB model. First, theoretical and numerical analyses will be made for Shan-Chen's forcing scheme [Shan and Chen, Phys. Rev. E 47, 1815 (1993)] and the exact-difference-method forcing scheme [Kupershtokh et al., Comput. Math. Appl. 58, 965 (2009)]. The nature of these two schemes and their recovered macroscopic equations will be shown. Second, through a theoretical analysis, we will reveal the physics behind the phenomenon that different forcing schemes exhibit different performances in the pseudopotential LB model. Moreover, based on the analysis, we will present an improved forcing scheme and numerically demonstrate that the improved scheme can be treated as an alternative approach to achieving thermodynamic consistency in the pseudopotential LB model.

Consistent description of electrohydrodynamics in narrow fluidic confinements in the presence of hydrophobic interactions

Electrohydrodynamics in the presence of hydrophobic interactions in narrow confinements is traditionally represented from a continuum viewpoint by a Navier slip-based conceptual paradigm, in which the slip length carries the sole burden of incorporating the effects of substrate wettability on interfacial electromechanics, precluding any explicit dependence of the interfacial potential distribution on the substrate wettability. Here we show that this traditional way of treating electrokinetics-wettability coupling may lead to serious discrepancies in predicting the resultant transport characteristics as manifested through an effective zeta potential. We suggest that an alternative consistent description of the underlying physics through a free-energy-based formalism, in conjunction with considerations of hydrodynamic and electrical property variations consistent with the pertinent phase-field description, may represent the underlying consequences in a more rational manner, as compared to the traditional slip-based model coupled with a two-layer description. Our studies further reveal that the above discrepancies may not occur solely due to the slip-based route of representing the interfacial wettability, but may be additionally attributed to the act of "discretizing" the interfacial phase fraction distribution through an artificial two-layer route.

Tensorial slip of superhydrophobic channels

We describe a generalization of the tensorial slip boundary condition, originally justified for a thick (compared to texture period) channel, to any channel thickness. The eigenvalues of the effective slip-length tensor, however, in general case become dependent on the gap and cannot be viewed as a local property of the surface, being a global characteristic of the channel. To illustrate the use of the tensor formalism we develop a semianalytical theory of an effective slip in a parallel-plate channel with one superhydrophobic striped and one hydrophilic surface. Our approach is valid for any local slip at the gas sectors and an arbitrary distance between the plates, ranging from a thick to a thin channel. We then present results of lattice Boltzmann simulations to validate the analysis. Our results may be useful for extracting effective slip tensors from global measurements, such as the permeability of a channel, in experiments or simulations.

Forcing term in single-phase and Shan-Chen-type multiphase lattice Boltzmann models

Numerous schemes have been proposed to incorporate a bulk forcing term into the lattice Boltzmann equation. In this paper we present a simple and straightforward comparative analysis of five popular schemes [Shan and Chen, Phys. Rev. E 47, 1815 (1993); Phys Rev Lett. 81, 1618 (1998); He et al., Phys. Rev. E 57, R13 (1998); Guo et al., Phys. Rev. E 65, 046308 (2002); Kupershtokh et al., Comput. Math. Appl. 58, 965 (2009)] in which their differences and similarities are identified. From the analysis we classify the schemes into two groups; the behaviors of the schemes in each group are proven to be identical up to second order. Numerical test simulating the two-dimensional unsteady Taylor-Green vortex flow problem demonstrate that all five schemes are of comparable accuracy for single-phase flow. However, for two-phase flow the situation is different, which is demonstrated by incorporating these schemes into different Shan-Chen-type multiphase models. The forcing scheme in the original Shan-Chen (SC) multiphase model turns out to be inaccurate in terms of the resulting surface tension for different density ratios and relaxation times. In the numerical tests, a typical equation of state and interparticle interactions including next-nearest neighbors were incorporated into the SC model. Our results confirm that the surface-tension values obtained from the original SC lattice Boltzmann method (LBM) simulation depend on the value of the relaxation time τ. For τ<0.7Δt, the surface tension agree well with the analytical solutions. However, when τ>0.7Δt, the surface tension turns out to be systematically larger than the analytical one, exceeding it by more than a factor of 2 for τ=2Δt. In contrast, with the application of the scheme proposed by He et al., the SC LBM produces very accurate surface tensions independent of the value of τ. We also found that the densities of the coexisting liquid and gas can be adjusted to match those at thermodynamic equilibrium if the particle interaction term includes next-nearest-neighbor contributions. The obtained results will be useful for further studies of two-phase flow with high density ratios using the SC LBM approach.

Consistent pseudopotential interactions in lattice Boltzmann models

We outline a systematic procedure on how to construct and derive interactions for nonideal lattice fluids. Using a mesoscopic approach based on exact lattice theories, we prove the consistency of the resulting diffuse interface theory with continuum thermodynamics. Translated in the framework of the so called "Shan-Chen" model for nonideal lattice fluids, this paper shows how to adjust the associated pseudopotentials to reproduce a free energy model based on a square gradient theory of equilibrium interfaces.

Modelling thermal fluctuations in non-ideal fluids with the lattice Boltzmann method

Recently, we proposed a theoretical framework to include thermal fluctuations into the Lattice Boltzmann (LB) method for non-ideal fluids. Here, we apply a variant thereof to a certain class of force-based non-ideal fluid LB models. We find that ideal-gas-like noise is an exact result of the fluctuation-dissipation theorem in the hydrodynamic regime. It is shown that satisfactory equilibration of the density and fluid momentum can be obtained in a simulation over a wide range of length scales.

Simulations of slip flow on nanobubble-laden surfaces

On microstructured hydrophobic surfaces, geometrical patterns may lead to the appearance of a superhydrophobic state, where gas bubbles at the surface can have a strong impact on the fluid flow along such surfaces. In particular, they can strongly influence a detected slip at the surface. We present two-phase lattice Boltzmann simulations of a flow over structured surfaces with attached gas bubbles and demonstrate how the detected slip depends on the pattern geometry, the bulk pressure, or the shear rate. Since a large slip leads to reduced friction, our results give assistance in the optimization of microchannel flows for large throughput.

Force imbalance in lattice Boltzmann equation for two-phase flows

The capability of modeling and simulating complex interfacial dynamics of multiphase flows has been recognized as one of the main advantages of the lattice Boltzmann equation (LBE). A basic feature of two-phase LBE models, i.e., the force balance condition at the discrete lattice level of LBE, is investigated in this work. An explicit force-balance formulation is derived for a flat interface by analyzing the two-dimensional nine-velocity (D2Q9) two-phase LBE model without invoking the Chapman-Enskog expansion. The result suggests that generally the balance between the interaction force and the pressure does not hold exactly on the discrete lattice due to numerical errors. It is also shown that such force imbalance can lead to some artificial velocities in the vicinity of phase interface.

Shear stress in nonideal fluid lattice Boltzmann simulations

The viscous stress in nonideal fluid lattice Boltzmann methods is investigated theoretically and by simulations. Three representative liquid-gas models are compared in a steady-state situation, where an analytical expression for the viscous stress is available. It is shown that, in the presence of nonideal fluid interactions or strong body forces, the accuracy of the computed viscous stress depends strongly on the underlying model.

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.

Fully reversible transition from Wenzel to Cassie-Baxter states on corrugated superhydrophobic surfaces

Liquid drops on textured surfaces show different dynamical behaviors depending on their wetting states. They are extremely mobile when they are supported by composite solid-liquid-air interfaces (Cassie-Baxter state) and immobile when they fully wet the textured surfaces (Wenzel state). By reversibly switching between these two states, it will be possible to achieve control over the fluid dynamics. Unfortunately, these wetting transitions are usually prevented by surface energy barriers. We demonstrate here a new, simple design paradigm consisting of parallel grooves with an appropriate aspect ratio that allows for the controlled, barrierless, reversible switching of the wetting states upon application of electrowetting. We report a direct observation of the barrierless dynamical pathway for the reversible transitions between the Wenzel (collapsed) and Cassie-Baxter (suspended) states and present a theory that accounts for these transitions, including detailed lattice Boltzmann simulations.

Multiphase lattice Boltzmann method for particle suspensions

A two-dimensional mass conserving lattice Boltzmann method (LBM) has been developed for multiphase (liquid and vapor) flows with solid particles suspended within the liquid and/or vapor phases. The main modification to previous single-phase particle suspension models is the addition of surface (adhesive) forces between the suspended particle and the surrounding fluid. The multiphase dynamics between fluid phases is simulated via the single-component multiphase model of Shan and Chen [Phys. Rev. E 47, 1815 (1993)]. The combined multiphase particle suspension model is first validated and then used to simulate the dynamics of a single-suspended particle on a planar liquid-vapor interface and the interaction between a single particle and a free-standing liquid drop. It is observed that the dynamics of suspended particles near free-standing liquid droplets is affected by spurious velocity currents although the liquid-vapor interface itself is a local energy minimum for particles. Finally, results are presented for capillary interactions between two suspended particles on a liquid-vapor interface subjected to different external forces and for spinodal decomposition of a liquid-vapor mixture in the presence of suspended particles. Qualitative agreements are reached when compared with results of suspended particles in a binary mixture based on multicomponent LBM models.

Contact line dynamics in binary lattice Boltzmann simulations

We show that, when a single relaxation time lattice Boltzmann algorithm is used to solve the hydrodynamic equations of a binary fluid for which the two components have different viscosities, strong spurious velocities in the steady state lead to incorrect results for the equilibrium contact angle. We identify the origins of these spurious currents and demonstrate how the results can be greatly improved by using a lattice Boltzmann method based on a multiple-relaxation-time algorithm. By considering capillary filling we describe the dependence of the advancing contact angle on the interface velocity.

Front pinning in capillary filling of chemically coated channels

The dynamics of capillary filling in the presence of chemically coated heterogeneous boundaries is investigated both theoretically and numerically. In particular, by mapping the equations of front motion onto the dynamics of a dissipative driven oscillator, an analytical criterion for front pinning is derived under the condition of diluteness of the coating spots. The criterion is tested against two-dimensional lattice Boltzmann simulations and found to provide satisfactory agreement as long as the width of the front interface remains much thinner than the typical heterogeneity scale of the chemical coating.

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.

Capillary filling with pseudo-potential binary Lattice-Boltzmann model

We present a systematic study of capillary filling for a binary fluid by using a mesoscopic lattice Boltzmann model for immiscible fluids describing a diffusive interface moving at a given contact angle with respect to the walls. The phenomenological way to impose a given contact angle is analysed. Particular attention is given to the case of complete wetting, that is contact angle equal to zero. Numerical results yield quantitative agreement with the theoretical Washburn's law, provided that the correct ratio of the dynamic viscosities between the two fluids is used. Finally, the presence of precursor films is experienced and it is shown that these films advance in time with a square-root law but with a different prefactor with respect to the bulk interface.

Wall boundary conditions in the lattice Boltzmann equation method for nonideal gases

Wall boundary conditions for the lattice Boltzmann equation (LBE) method for nonideal gases proposed by Lee and Fischer [Phys. Rev. E 74, 046709 (2006)] are examined. The LBE simulations of the contact line are typically contaminated by small but strong counter-rotating parasitic currents near solid surfaces. We find that these currents can be eliminated to round-off if the potential form of the intermolecular force is used with the boundary conditions based on the wall energy approach and the bounce-back rule. Numerical tests confirm the elimination of the parasitic currents in the vicinity of the contact line and the agreement of the calculated density, excess mass, and contact angle at the solid surfaces with theory.

Slip flow over structured surfaces with entrapped microbubbles

On hydrophobic surfaces, roughness may lead to a transition to a superhydrophobic state, where gas bubbles at the surface can have a strong impact on a detected slip. We present two-phase lattice Boltzmann simulations of a Couette flow over structured surfaces with attached gas bubbles. Even though the bubbles add slippery surfaces to the channel, they can cause negative slip to appear due to the increased roughness. The simulation method used allows the bubbles to deform due to viscous stresses. We find a decrease of the detected slip with increasing shear rate which is in contrast to some recent experimental results implicating that bubble deformation cannot account for these experiments. Possible applications of bubble surfaces in microfluidic devices are discussed.

Pressure tensor calculation in a class of nonideal gas lattice Boltzmann models

In nonideal gas lattice Boltzmann (LB) models, obtaining the correct form of the pressure tensor is essential in determining many of the statistical mechanical properties such as the surface tension and the density profile. Here we outline a general approach for calculating the pressure tensor in LB models with interactions beyond nearest neighbors. The statistical mechanical properties calculated from such a pressure tensor are shown to agree very well with those measured from numerical experiments. Comparisons with alternative theories are also made.

Proposed approximation for contact angles in Shan-and-Chen-type multicomponent multiphase lattice Boltzmann models

We propose a method for approximating the adhesion parameters in the Shan and Chen multicomponent, multiphase lattice Boltzmann model that leads to the desired fluid-solid contact angle. The method is a straightforward application of Young's equation with substitution of the Shan and Chen cohesion parameter and a density factor for the fluid-fluid interfacial tension, and the adhesion parameters for the corresponding fluid-solid interfacial tensions.

Spontaneous breakdown of superhydrophobicity

In some cases water droplets can completely wet microstructured superhydrophobic surfaces. The dynamics of this rapid process is analyzed by ultrahigh-speed imaging. Depending on the scales of the microstructure, the wetting fronts propagate smoothly and circularly or-more interestingly-in a stepwise manner, leading to a growing square-shaped wetted area: entering a new row perpendicular to the direction of front propagation takes milliseconds, whereas once this has happened, the row itself fills in microseconds ("zipping"). Numerical simulations confirm this view and are in quantitative agreement with the experiments.

Investigations of drop impact on dry walls with a lattice-Boltzmann model

In this work, axisymmetric computations of drop impingement on dry walls are presented. The two-phase model employed is an axisymmetric lattice-Boltzmann model. Computations are performed in the parametric range of Weber number We of 7 to 8770, Ohnesorge number Oh of 0.02 to 0.137, and drop-wall equilibrium contact angle theta(eq) of 35 degrees to 150 degrees . Deposition and rebound outcomes are reported. In deposition, the different stages of drop evolution including spread, recoil and equilibration are reproduced and studied. Comparisons made with experimentally reported data of temporal evolution of the spread factor and the dynamic evolution of the contact angle show good agreement. Rebound is observed on non-wetting surfaces. The transition between deposition and rebound is shown to be influenced by the impact We, Oh, and advancing and receding static contact angles. Apart from impingement outcomes, the influence of We and Oh on the dynamic contact angle is investigated.

Droplets on inclined rough surfaces

The behaviour of liquid droplets on inclined heterogeneous surfaces was simulated by the lattice-Boltzmann method using the Shan-Chen multiphase model. The effect of topography of the surface on the contact angle hysteresis was investigated. It is shown in particular, by using anisotropic rough surfaces, how surface topography and thereby the continuity of the three-phase contact line, affect this hysteresis. Our results clearly indicate that the superhydrophobicity of a surface cannot be judged by the contact angle alone.

Generalized lattice Boltzmann method with multirange pseudopotential

The physical behavior of a class of mesoscopic models for multiphase flows is analyzed in details near interfaces. In particular, an extended pseudopotential method is developed, which permits to tune the equation of state and surface tension independently of each other. The spurious velocity contributions of this extended model are shown to vanish in the limit of high grid refinement and/or high order isotropy. Higher order schemes to implement self-consistent forcings are rigorously computed for 2d and 3d models. The extended scenario developed in this work clarifies the theoretical foundations of the Shan-Chen methodology for the lattice Boltzmann method and enhances its applicability and flexibility to the simulation of multiphase flows to density ratios up to O(100).

Surface roughness-hydrophobicity coupling in microchannel and nanochannel flows

An approach based on a lattice version of the Boltzmann kinetic equation for describing multiphase flows in nano- and microcorrugated devices is proposed. We specialize it to describe the wetting-dewetting transition of fluids in the presence of nanoscopic grooves etched on the boundaries. This approach permits us to retain the essential supramolecular details of fluid-solid interactions without surrendering--actually boosting--the computational efficiency of continuum methods. The method is used to analyze the importance of conspiring effects between hydrophobicity and roughness on the global mass flow rate of the microchannel. In particular we show that smart surfaces can be tailored to yield very different mass throughput by changing the bulk pressure. The mesoscopic method is also validated quantitatively against the molecular dynamics results of [Cottin-Bizonne, Nat. Mater. 2, 237 (2003)].