Over the last two decades, lattice Boltzmann methods have become an increasingly popular tool to compute the flow in complex geometries such as porous media. In addition to single phase simulations allowing, for example, a precise quantification of the permeability of a porous sample, a number of extensions to the lattice Boltzmann method are available which allow to study multiphase and multicomponent flows on a pore scale level. In this article, we give an extensive overview on a number of these diffuse interface models and discuss their advantages and disadvantages. Furthermore, we shortly report on multiphase flows containing solid particles, as well as implementation details and optimization issues.

A new algorithm for imposing the contact angle on solid surfaces is proposed in the Lattice Boltzmann color‐gradient model. The capability and accuracy of this algorithm are validated by simulation of contact angles for a droplet resting on a flat surface and on a cylindrical surface. The color‐gradient model with the proposed contact angle algorithm is then used to study the capillary valve effect in porous media. As a preliminary study, the capillary valve effect is explained by simulating immiscible two‐phase displacement within a single‐pore geometry. It is shown that the capillary valve effect is accurately captured by the present simulations. Further simulations of drainage and imbibition are also conducted to understand the capillary valve effect in an experiment‐matched pore‐network micromodel. The simulated results are found to agree quantitatively with the experimental results reported in literature, except for a few differences which result from the exclusion of contact angle hysteresis in the proposed algorithm.

Droplet dynamics in microfluidic applications is significantly influenced by surfactants. It remains a research challenge to model and simulate droplet behaviour including deformation, breakup and coalescence, especially in the confined microfluidic environment. Here, we propose a hybrid method to simulate interfacial flows with insoluble surfactants. The immiscible two-phase flow is solved by an improved lattice Boltzmann colour-gradient model which incorporates a Marangoni stress resulting from non-uniform interfacial tension, while the convection–diffusion equation which describes the evolution of surfactant concentration in the entire fluid domain is solved by a finite difference method. The lattice Boltzmann and finite difference simulations are coupled through an equation of state, which describes how surfactant concentration influences interfacial tension. Our method is first validated for the surfactant-laden droplet deformation in a three-dimensional (3D) extensional flow and a 2D shear flow, and then applied to investigate the effect of surfactants on droplet dynamics in a 3D shear flow. Numerical results show that, at low capillary numbers, surfactants increase droplet deformation, due to reduced interfacial tension by the average surfactant concentration, and non-uniform effects from non-uniform capillary pressure and Marangoni stresses. The role of surfactants on the critical capillary number ( $Ca_{cr}$ ) of droplet breakup is investigated for various confinements (defined as the ratio of droplet diameter to wall separation) and Reynolds numbers. For clean droplets, $Ca_{cr}$ first decreases and then increases with confinement, and the minimum value of $Ca_{cr}$ is reached at a confinement of 0.5; for surfactant-laden droplets, $Ca_{cr}$ exhibits the same variation in trend for confinements lower than 0.7, but, for higher confinements, $Ca_{cr}$ is almost a constant. The presence of surfactants decreases $Ca_{cr}$ for each confinement, and the decrease is also attributed to the reduction in average interfacial tension and non-uniform effects, which are found to prevent droplet breakup at low confinements but promote breakup at high confinements. In either clean or surfactant-laden cases, $Ca_{cr}$ first remains almost unchanged and then decreases with increasing Reynolds number, and a higher confinement or Reynolds number favours ternary breakup. Finally, we study the collision of two equal-sized droplets in a shear flow in both surfactant-free and surfactant-contaminated systems with the same effective capillary numbers. It is identified that the non-uniform effects in the near-contact interfacial region immobilize the interfaces when two droplets are approaching each other and thus inhibit their coalescence.

Contact angle hysteresis is an important physical phenomenon omnipresent in nature and various industrial processes, but its effects are not considered in many existing multiphase flow simulations due to modeling complexity. In this work, a multiphase lattice Boltzmann method (LBM) is developed to simulate the contact-line dynamics with consideration of the contact angle hysteresis for a broad range of kinematic viscosity ratios. In this method, the immiscible two-phase flow is described by a color-fluid model, in which the multiple-relaxation-time collision operator is adopted to increase numerical stability and suppress unphysical spurious currents at the contact line. The contact angle hysteresis is introduced using the strategy proposed by Ding and Spelt [Ding and Spelt, J. Fluid Mech. 599, 341 (2008)], and the geometrical wetting boundary condition is enforced to obtain the desired contact angle. This method is first validated by simulations of static contact angle and dynamic capillary intrusion process on ideal (smooth) surfaces. It is then used to simulate the dynamic behavior of a droplet on a nonideal (inhomogeneous) surface subject to a simple shear flow. When the droplet remains pinned on the surface due to hysteresis, the steady interface shapes of the droplet quantitatively agree well with the previous numerical results. Four typical motion modes of contact points, as observed in a recent study, are qualitatively reproduced with varying advancing and receding contact angles. The viscosity ratio is found to have a notable impact on the droplet deformation, breakup, and hysteresis behavior. Finally, this method is extended to simulate the droplet breakup in a microfluidic T junction, with one half of the wall surface ideal and the other half nonideal. Due to the contact angle hysteresis, the droplet asymmetrically breaks up into two daughter droplets with the smaller one in the nonideal branch channel, and the behavior of daughter droplets is significantly different in both branch channels. Also, it is found that the contact angle hysteresis is strengthened with decreasing the viscosity ratio, leading to an earlier droplet breakup and a decrease in the maximum length that the droplet can reach before the breakup. These simulation results manifest that the present multiphase LBM can be a useful substitute to Ba et al. [ Phys. Rev. E 88, 043306 (2013)] for modeling the contact angle hysteresis, and it can be easily implemented with higher computational efficiency.

Double emulsion formation in a hierarchical flow-focusing channel is systematically investigated using a free-energy ternary lattice Boltzmann model. A three-dimensional formation regime diagram is constructed based on the capillary numbers of the inner ( $Ca_{i}$ ), middle ( $Ca_{m}$ ) and outer ( $Ca_{o}$ ) phase fluids. The results show that the formation diagram can be classified into periodic two-step region, periodic one-step region, and non-periodic region. By varying $Ca_{i}$ and $Ca_{m}$ in the two-step formation region, different morphologies are obtained, including the regular double emulsions, decussate regimes with one or two alternate empty droplets, and structures with multiple inner droplets contained in the continuous middle phase thread. Bidisperse behaviours are also frequently encountered in the two-step formation region. In the periodic one-step formation region, scaling laws are proposed for the double emulsion size and for the size ratio between the inner droplet and the overall double emulsion. Furthermore, we show that the interfacial tension ratio can greatly change the morphologies of the obtained emulsion droplets, and the channel geometry plays an important role in changing the formation regimes and the double emulsion sizes. In particular, narrowing the side inlets or the distance between the two side inlets promotes the conversion from the two-step formation regime to the one-step formation regime.

We propose a lattice Boltzmann color-gradient model for immiscible ternary fluid flows, which is applicable to the fluids with a full range of interfacial tensions, especially in near-critical and critical states. An interfacial force for N-phase systems is derived and then introduced into the model using a body force scheme, which helps reduce spurious velocities. A generalized recoloring algorithm is applied to produce phase segregation and ensure immiscibility of three different fluids, where an enhanced form of segregation parameters is derived by considering the existence of Neumann’s triangle and the effect of the equilibrium contact angle in a three-phase junction. The proposed model is first validated by two typical examples, namely, the Young-Laplace test for a compound droplet and the spreading of a droplet between two stratified fluids. It is then used to study the structure and stability of double droplets in a static matrix. Consistent with the theoretical stability diagram, seven possible equilibrium morphologies are successfully reproduced by adjusting the interfacial tension ratio. By simulating near-critical and critical states of double droplets where the outcomes are very sensitive to the model accuracy, we show that the present model is advantageous to three-phase flow simulations and allows for accurate simulation of near-critical and critical states. Finally, we investigate the influence of interfacial tension ratio on the behavior of a compound droplet in a three-dimensional shear flow, and four different deformation and breakup modes are observed.