Ternary free-energy lattice Boltzmann model with tunable surface tensions and contact angles

Dynamics of spreading of an asymmetrically placed droplet near a fluid-fluid interface

Two-dimensional numerical simulations are carried out to study the spreading dynamics of a droplet placed in the vicinity of a fluid-fluid interface. Simulations are performed using the hybrid lattice-Boltzmann technique and the diffuse-interface model by considering three immiscible fluids of the same density and viscosity. In contrast to the well-studied spreading of drops placed symmetrically across fluid-fluid interfaces, this work considers the simultaneous migration, spreading and eventual adsorption of an asymmetrically placed drop. These processes, which are solely driven by interfacial forces, are characterised by monitoring the temporal evolution of geometric parameters, such as the centre of mass, radius and height of the drop, the surface energy of the three interfaces and the associated flow fields inside and outside the droplet. The rate of spreading and rate of adsorption are also calculated to determine the dominant processes that drive the dynamics of the system.

General wetting energy boundary condition in a fully explicit nonideal fluids solver

We present an explicit finite-difference method to simulate the nonideal multiphase fluid flow. The local density and momentum transport are modeled by the Navier-Stokes equations and the pressure is computed by the van der Waals equation of the state. The static droplet and the dynamics of liquid-vapor separation simulations are performed as validations of this numerical scheme. In particular, to maintain the thermodynamic consistency, we propose a general wetting energy boundary condition at the contact line between fluids and the solid boundary. We conduct a series of comparisons between the current boundary condition and the constant contact angle boundary condition as well as the stress-balanced boundary condition. This boundary condition alleviates the instability induced by the constant contact angle boundary condition at θ≈0 and θ≈π. Using this boundary condition, the equilibrium contact angle is correctly recovered and the contact line dynamics are consistent with the simulation by applying a stress-balanced boundary condition. Nevertheless, unlike the stress-balanced boundary condition for which we need to further introduce the interface thickness parameter, the current boundary condition implicitly incorporates the interface thickness information into the wetting energy.

Unrestricted component count in multiphase lattice Boltzmann: A fugacity-based approach

Studies of multiphase fluids utilizing the lattice Boltzmann method (LBM) are typically severely restricted by the number of components or chemical species being modeled. This restriction is particularly pronounced for multiphase systems exhibiting partial miscibility and significant interfacial mass exchange, which is a common occurrence in realistic multiphase systems. Modeling such systems becomes increasingly complex as the number of chemical species increases due to the increased role of molecular interactions and the types of thermodynamic behavior that become possible. The recently introduced fugacity-based LBM [Soomro et al., Phys. Rev. E 107, 015304 (2023)2470-004510.1103/PhysRevE.107.015304] has provided a thermodynamically consistent modeling platform for multicomponent, partially miscible LBM simulations. However, until now, this fugacity-based LB model had lacked a comprehensive demonstration of its ability to accurately reproduce thermodynamic behavior beyond binary mixtures and to remove any restrictions in a number of components for multiphase LBM. In this paper we closely explore these fugacity-based LBM capabilities by showcasing comprehensive, thermodynamically consistent simulations of multiphase mixtures of up to ten chemical components. The paper begins by validating the model against the Young-Laplace equation for a droplet composed of three components. The model is then applied to study mixtures with a range of component numbers from one to six, showing agreement with rigorous thermodynamic predictions and demonstrating linear scaling of computational time with the number of components. We further investigate ternary systems in detail by exploring a wide range of temperature, pressure, and overall composition conditions to produce various characteristic ternary diagrams. In addition, the model is shown to be unrestricted in the number of phases as demonstrated through simulations of a three-component three-phase equilibrium case. The paper concludes by demonstrating simulations of a ten-component, realistic hydrocarbon mixture, achieving excellent agreement with thermodynamics for both flat interface vapor-liquid equilibrium and curved interface spinodal decomposition cases. This study represents a significant expansion of the scope and capabilities of multiphase LBM simulations that encompass multiphase systems of keen interest in engineering.

Equation-of-state-dependent surface free-energy density for wettability in lattice Boltzmann method

In thermodynamic theory, the liquid-vapor fluids can be described by a single multiphase equation of state and the surface wettability is usually characterized by the surface free-energy density. In this work, we propose an equation-of-state-dependent surface free-energy density for the wettability of the liquid-vapor fluids on a solid surface, which can lead to a simple closed-form analytical expression for the contact angle. Meanwhile, the thermodynamically derived equilibrium condition is equivalent to the geometric formulation of the contact angle. To numerically validate the present surface free-energy density, the mesoscopic multiphase lattice Boltzmann model with self-tuning equation of state, which is strictly consistent with thermodynamic theory, is employed, and the two-dimensional wetting condition treatment is extended to the three-dimensional situation with flat and curved surfaces. Two- and three-dimensional lattice Boltzmann simulations of static droplets on flat and curved surfaces are first performed, and the obtained contact angles agree well with the closed-form analytical expression. Then, the three-dimensional lattice Boltzmann simulation of a moving droplet on an inclined wall, which is vertically and sinusoidally oscillated, is carried out. The dynamic contact angles well satisfy the Cox-Voinov law. The droplet movement regimes are consistent with previous experiments and two-dimensional simulations. The dependence of the droplet overall velocity with respect to the dimensionless oscillation strength is also discussed in detail.

Fugacity-based lattice Boltzmann method for multicomponent multiphase systems

The free-energy model can extend the lattice Boltzmann method to multiphase systems. However, there is a lack of models capable of simulating multicomponent multiphase fluids with partial miscibility. In addition, existing models cannot be generalized to honor thermodynamic information provided by any multicomponent equation of state of choice. In this paper, we introduce a free-energy lattice Boltzmann model where the forcing term is determined by the fugacity of the species, the thermodynamic property that connects species partial pressure to chemical potential calculations. By doing so, we are able to carry out multicomponent multiphase simulations of partially miscible fluids and generalize the methodology for use with any multicomponent equation of state of interest. We test this fugacity-based lattice Boltzmann method for the cases of vapor-liquid equilibrium for two- and three-component mixtures in various temperature and pressure conditions. We demonstrate that the model is able to reliably reproduce phase densities and compositions as predicted by multicomponent thermodynamics and can reproduce different characteristic pressure-composition and temperature-composition envelopes with a high degree of accuracy. We also demonstrate that the model can offer accurate predictions under dynamic conditions.

Spontaneous phase separation of ternary fluid mixtures

We computationally study the spontaneous phase separation of ternary fluid mixtures using the lattice Boltzmann method both when all the surface tensions are equal and when they have different values. To rationalise the phase diagram of possible phase separation mechanisms, previous theoretical works typically rely on analysing the sign of the eigenvalues resulting from a simple linear stability analysis, but we find this does not explain the observed simulation results. Here, we classify the possible separation pathways into four basic mechanisms, and develop a phenomenological model that captures the composition regimes where each mechanism is prevalent. We further highlight that the dominant mechanism in ternary phase separation involves enrichment and instability of the minor component at the fluid-fluid interface, which is absent in the case of binary fluid mixtures.

Phase separation dynamics in deformable droplets

Phase separation can drive spatial organization of multicomponent mixtures. For instance in developing animal embryos, effective phase separation descriptions have been used to account for the spatial organization of different tissue types. Similarly, separation of different tissue types is also observed in stem cell aggregates, where the emergence of a polar organization can mimic early embryonic axis formation. Here, we describe such aggregates as deformable two-phase fluid droplets, which are suspended in a fluid environment (third phase). Using hybrid finite-volume Lattice-Boltzmann simulations, we numerically explore the out-of-equilibrium routes that can lead to the polar equilibrium state of such a droplet. We focus on the interplay between spinodal decomposition and advection with hydrodynamic flows driven by interface tensions, which we characterize by a Peclet number Pe. Consistent with previous work, for large Pe the coarsening process is generally accelerated. However, for intermediate Pe we observe long-lived, strongly elongated droplets, where both phases form an alternating stripe pattern. We show that these "croissant" states are close to mechanical equilibrium and coarsen only slowly through diffusive fluxes in an Ostwald-ripening-like process. Finally, we show that a surface tension asymmetry between both droplet phases leads to transient, rotationally symmetric states whose resolution leads to flows reminiscent of Marangoni flows. Our work highlights the importance of advection for the phase separation process in finite, deformable systems.

Development of a coupled geophysical-geothermal scheme for quantification of hydrates in gas hydrate-bearing permafrost sediments

Quantification of hydrates in permafrost sediments using conventional seismic techniques has always been a major challenge in the study of the climate-driven evolution of gas hydrate-bearing permafrost sediments due to almost identical acoustic properties of hydrates and ice. In this article, a coupled geophysical-geothermal scheme is developed, for the first time, to predict hydrate saturation in gas hydrate-bearing permafrost sediments by utilising their geophysical and geothermal responses. The scheme includes a geophysical part which interprets the measured elastic wave velocities using a rock-physics model, coupled with a geothermal part, interpreting the measured effective thermal conductivity (ETC) using a new pore-scale model. By conducting a series of sensitivity analyses, it is shown that the ETC model is able to incorporate the effect of the hydrate pore-scale habit and hydrate/ice-forced heave as well as the effect of unfrozen water saturation under frozen conditions. Given that the geophysical and geothermal responses depend on the overburden pressure, the elastic wave velocities and ETC of methane hydrate-bearing permafrost sediment samples were measured at different effective overburden pressures and the results were provided. These experimental data together with the results of our recent study on the geophysical and geothermal responses of gas hydrate-bearing permafrost sediment samples at different hydrate saturations are used to validate the performance of the coupled scheme. By comparing the predicted saturations with those obtained experimentally, it is shown that the coupled scheme is able to quantify the saturation of the co-existing phases with an acceptable accuracy in a wide range of hydrate saturations and at different overburden pressures.

Hysteresis in spreading and retraction of liquid droplets on parallel fiber rails

Wetting and spreading of liquids on fibers occur in many natural and artificial processes. Unlike on a planar substrate, a droplet attached to one or more fibers can assume several different shapes depending on geometrical parameters such as liquid volume and fiber size and distance. This paper presents lattice Boltzmann simulations of the morphology of liquid droplets on two parallel cylindrical fibers. We investigate the final shapes resulting from spreading of an initially spherical droplet deposited on the fibers and from retraction of an initial liquid column deposited between the fibers. We observe three possible equilibrium configurations: barrel-shaped droplet, droplet bridges, and liquid columns. We determine the complete morphology diagram for varying inter-fiber spacings and liquid volumes and find a region of bistability that spans both the column regime and the droplet regime. We further present a simulation protocol that allows one to probe the hysteresis of transitions between different shapes. The results provide insights into energies and forces associated with shape transformations of droplets on fibers that can be used to develop fiber-based materials and microfluidic systems for manipulation of liquids at a small scale.

Modeling ternary fluids in contact with elastic membranes

We present a thermodynamically consistent model of a ternary fluid interacting with elastic membranes. Following a free-energy modeling approach for the fluid phases, we derive the governing equations for the dynamics of the ternary fluid flow and membranes. We also provide the numerical framework for simulating such fluid-structure interaction problems. It is based on the lattice Boltzmann method for the ternary fluid (Eulerian description) and a finite difference representation of the membrane (Lagrangian description). The ternary fluid and membrane solvers are coupled through the immersed boundary method. For validation purposes, we consider the relaxation dynamics of a two-dimensional elastic capsule placed at a fluid-fluid interface. The capsule shapes, resulting from the balance of surface tension and elastic forces, are compared with equilibrium numerical solutions obtained by surface evolver. Furthermore, the Galilean invariance of the proposed model is proven. The proposed approach is versatile, allowing for the simulation of a wide range of geometries. To demonstrate this, we address the problem of a capillary bridge formed between two deformable capsules.

Factors controlling the pinning force of liquid droplets on liquid infused surfaces

Liquid infused surfaces with partially wetting lubricants have recently been exploited for numerous intriguing applications, such as for droplet manipulation, droplet collection and spontaneous motion. When partially wetting lubricants are used, the pinning force is a key factor that can strongly affect droplet mobility. Here, we derive an analytical prediction for contact angle hysteresis in the limit where the meniscus size is much smaller than the droplet, and numerically study how it is controlled by the solid fraction, the lubricant wetting angles, and the various fluid surface tensions. We further relate the contact angle hysteresis and the pinning force experienced by a droplet on a liquid infused surface, and our predictions for the critical sliding angles are consistent with existing experimental observations. Finally, we discuss why a droplet on a liquid infused surface with partially wetting lubricants typically experiences stronger pinning compared to a droplet on a classical superhydrophobic surface.

Lattice Boltzmann model for ternary fluids with solid particles

On the basis of phase-field theory, we develop a lattice Boltzmann model for ternary fluids containing solid. We develop a modified bounce-back method to describe the interactions between the solid and N-phase (N≥2) fluids. We derive a wetting boundary condition for three-phase flows from the point of mass conservation and propose a scheme for implementing the wetting condition on curved boundaries. We develop a diffuse interface method to compute the capillary force acting on the moving solid objects at the ternary fluids-sold contact lines. In addition, this model can deal with problems involving high density and viscosity contrasts. The proposed method is examined through several test cases. We test the modified bounce-back scheme, wetting boundary condition, and capillary force model in three different cases, and the numerical results agree well with the analytical solutions. Finally, we apply the model to two three-dimensional problems to assess its numerical accuracy and stability.

Reduction-consistent phase-field lattice Boltzmann equation for N immiscible incompressible fluids

In this paper, we develop a reduction-consistent conservative phase-field method for interface-capturing among N (N≥ 2) immiscible fluids, which is governed by conservative Allen-Cahn equation (CACE); here the reduction-consistent property is that if only M (1≤M≤N-1) immiscible fluids are present in a N-phase system, the governing equations for N immiscible fluids must reduce to the corresponding M immiscible fluids system. Then we propose a reduction-consistent lattice Boltzmann equation (LBE) method for solving N immiscible incompressible fluids with high density and viscosity contrasts. Some numerical simulations are carried out to validate the present LBE such as stationary droplets, spreading of a liquid lens, and spinodal decomposition together with the reduction-consistent property, and the numerical results predicted by present LBE are in good agreement with the analytical solutions/other numerical results, which also demonstrate the reduction-consistent property by present LBE.

Multiphase flows of N immiscible incompressible fluids: Conservative Allen-Cahn equation and lattice Boltzmann equation method

In this paper, we develop a conservative phase-field method for interface-capturing among N (N≥2) immiscible fluids, the evolution of the fluid-fluid interface is captured by conservative Allen-Cahn equation (CACE), and the interface force of N immiscible fluids is incorporated to Navier-Stokes equation (NSE) by chemical potential form. Accordingly, we propose a lattice Boltzmann equation (LBE) method for solving N (N≥2) immiscible incompressible NSE and CACE at high density and viscosity contrasts. Numerical simulations including stationary droplets, Rayleigh-Taylor instability, spreading of liquid lenses, and spinodal decompositions are carried out to show the accuracy and capability of present LBE, and the results show that the predictions by use of the present LBE agree well with the analytical solutions and/or other numerical results.

Multicomponent flow on curved surfaces: A vielbein lattice Boltzmann approach

We develop and implement a finite difference lattice Boltzmann scheme to study multicomponent flows on curved surfaces, coupling the continuity and Navier-Stokes equations with the Cahn-Hilliard equation to track the evolution of the binary fluid interfaces. The standard lattice Boltzmann method relies on regular Cartesian grids, which makes it generally unsuitable to study flow problems on curved surfaces. To alleviate this limitation, we use a vielbein formalism to write the Boltzmann equation on an arbitrary geometry, and solve the evolution of the fluid distribution functions using a finite difference method. Focusing on the torus geometry as an example of a curved surface, we demonstrate drift motions of fluid droplets and stripes embedded on the surface of the torus. Interestingly, they migrate in opposite directions: fluid droplets to the outer side while fluid stripes to the inner side of the torus. For the latter we demonstrate that the global minimum configuration is unique for small stripe widths, but it becomes bistable for large stripe widths. Our simulations are also in agreement with analytical predictions for the Laplace pressure of the fluid stripes, and their damped oscillatory motion as they approach equilibrium configurations, capturing the corresponding decay timescale and oscillation frequency. Finally, we simulate the coarsening dynamics of phase separating binary fluids in the hydrodynamics and diffusive regimes for tori of various shapes, and compare the results against those for a flat two-dimensional surface. Our finite difference lattice Boltzmann scheme can be extended to other surfaces and coupled to other dynamical equations, opening up a vast range of applications involving complex flows on curved geometries.

Wetting boundaries for a ternary high-density-ratio lattice Boltzmann method

We extend a recently proposed ternary free-energy lattice Boltzmann model with high density contrast [Phys. Rev. Lett. 120, 234501 (2018)PRLTAO0031-900710.1103/PhysRevLett.120.234501] by incorporating wetting boundaries at solid walls. The approaches are based on forcing and geometric schemes, with implementations optimized for ternary (and, more generally, higher-order multicomponent) models. Advantages and disadvantages of each method are addressed by performing both static and dynamic tests, including the capillary filling dynamics of a liquid displacing the gas phase and the self-propelled motion of a train of drops. Furthermore, we measure dynamic angles and show that the slip length critically depends on the equilibrium value of the contact angles and whether it belongs to liquid-liquid or liquid-gas interfaces. These results validate the model capabilities of simulating complex ternary fluid dynamic problems near solid boundaries, for example, drop impact solid substrates covered by a lubricant layer.

Phase-field-theory-based lattice Boltzmann equation method for N immiscible incompressible fluids

From the phase field theory, we develop a lattice Boltzmann equation (LBE) method for N (N≥2) immiscible incompressible fluids, and the Cahn-Hilliard equation, which could capture the interfaces between different phases, is also solved by LBE for an N-phase system. In this model, the interface force of N immiscible incompressible fluids is incorporated by chemical potential form, and the fluid-fluid surface tensions could be directly calculated and independently tuned. Numerical simulations including two stationary droplets, spreading of a liquid lens with and without gravity and two immiscible liquid lenses, and phase separation are conducted to validate the present LBE, and numerical results show that the predictions by LBE agree well with the analytical solutions and other numerical results.

Lattice Boltzmann simulation of immiscible three-phase flows with contact-line dynamics

A multiphase lattice Boltzmann method is developed to simulate immiscible three-phase flows with contact-line dynamics. In this method, the immiscible three-phase flow is modeled by a multiple-relaxation-time color-gradient model, which not only allows for a full range of interfacial tensions but also can produce viscosity-independent results especially when the fluid-surface interactions are considered. To achieve the desired contact angles, a weighted contact angle model is utilized to obtain a relatively smooth transition of contact angle for each fluid, which is enforced through a geometrical wetting condition. This method is first validated by simulations of a Janus droplet resting on a surface for various contact angles and fluid properties and dynamic capillary filling of ternary fluids with different viscosity ratios. It is then used to simulate a Janus droplet on a substrate subject to Poiseuille flow. Results show that the droplet may undergo three typical modes, namely, two stable deformation modes and breakup mode, which depend not only on the inlet velocity but also on the fluid viscosity. The terminal velocity of moving droplet increases linearly with the inlet velocity in both stable modes only when three fluids do not differ much in their viscosities.

Drop Dynamics on Liquid-Infused Surfaces: The Role of the Lubricant Ridge

We employ a free-energy lattice-Boltzmann method to study the dynamics of a ternary fluid system consisting of a liquid drop driven by a body force across a regularly textured substrate, infused by a lubricating liquid. We focus on the case of partial wetting lubricants and observe a rich interplay between contact line pinning and viscous dissipation at the lubricant ridge, which become dominant at large and small apparent angles, respectively. Our numerical investigations further demonstrate that the relative importance of viscous dissipation at the lubricant ridge depends on the drop to lubricant viscosity ratio, as well as on the shape of the wetting ridge.

Ternary Free-Energy Entropic Lattice Boltzmann Model with a High Density Ratio

A thermodynamically consistent free energy model for fluid flows comprised of one gas and two liquid components is presented and implemented using the entropic lattice Boltzmann scheme. The model allows a high density ratio, up to the order of O(10^{3}), between the liquid and gas phases, and a broad range of surface tension ratios, covering partial wetting states where Neumann triangles are formed, and full wetting states where complete encapsulation of one of the fluid components is observed. We further demonstrate that we can capture the bouncing, adhesive, and insertive regimes for the binary collisions between immiscible droplets suspended in air. Our approach opens up a vast range of multiphase flow applications involving one gas and several liquid components.

Lattice Boltzmann model for three-phase viscoelastic fluid flow

A lattice Boltzmann (LB) framework is developed for simulation of three-phase viscoelastic fluid flows in complex geometries. This model is based on a Rothman-Keller type model for immiscible multiphase flows which ensures mass conservation of each component in porous media even for a high density ratio. To account for the viscoelastic effects, the Maxwell constitutive relation is correctly introduced into the momentum equation, which leads to a modified lattice Boltzmann evolution equation for Maxwell fluids by removing the normal but excess viscous term. Our simulation tests indicate that this excess viscous term may induce significant errors. After three benchmark cases, the displacement processes of oil by dispersed polymer are studied as a typical example of three-phase viscoelastic fluid flow. The results show that increasing either the polymer intrinsic viscosity or the elastic modulus will enhance the oil recovery.

Numerical simulation of three-component multiphase flows at high density and viscosity ratios using lattice Boltzmann methods

In this paper, we propose a multiphase lattice Boltzmann model for numerical simulation of ternary flows at high density and viscosity ratios free from spurious velocities. The proposed scheme, which is based on the phase-field modeling, employs the Cahn-Hilliard theory to track the interfaces among three different fluid components. Several benchmarks, such as the spreading of a liquid lens, binary droplets, and head-on collision of two droplets in binary- and ternary-fluid systems, are conducted to assess the reliability and accuracy of the model. The proposed model can successfully simulate both partial and total spreadings while reducing the parasitic currents to the machine precision.

Fluid-structure interaction with the entropic lattice Boltzmann method

We propose a fluid-structure interaction (FSI) scheme using the entropic multi-relaxation time lattice Boltzmann (KBC) model for the fluid domain in combination with a nonlinear finite element solver for the structural part. We show the validity of the proposed scheme for various challenging setups by comparison to literature data. Beyond validation, we extend the KBC model to multiphase flows and couple it with a finite element method (FEM) solver. Robustness and viability of the entropic multi-relaxation time model for complex FSI applications is shown by simulations of droplet impact on elastic superhydrophobic surfaces.

Lattice Boltzmann simulations of heat transfer in fully developed periodic incompressible flows

Flow and heat transfer in periodic structures are of great interest for many applications. In this paper, we carefully examine the periodic features of fully developed periodic incompressible thermal flows, and incorporate them in the lattice Boltzmann method (LBM) for flow and heat transfer simulations. Two numerical approaches, the distribution modification (DM) approach and the source term (ST) approach, are proposed; and they can both be used for periodic thermal flows with constant wall temperature (CWT) and surface heat flux boundary conditions. However, the DM approach might be more efficient, especially for CWT systems since the ST approach requires calculations of the streamwise temperature gradient at all lattice nodes. Several example simulations are conducted, including flows through flat and wavy channels and flows through a square array with circular cylinders. Results are compared to analytical solutions, previous studies, and our own LBM calculations using different simulation techniques (i.e., the one-module simulation vs. the two-module simulation, and the DM approach vs. the ST approach) with good agreement. These simple, however, representative simulations demonstrate the accuracy and usefulness of our proposed LBM methods for future thermal periodic flow simulations.

Modeling mass transfer and reaction of dilute solutes in a ternary phase system by the lattice Boltzmann method

In this work, we propose a general approach for modeling mass transfer and reaction of dilute solute(s) in incompressible three-phase flows by introducing a collision operator in lattice Boltzmann (LB) method. An LB equation was used to simulate the solute dynamics among three different fluids, in which the newly expanded collision operator was used to depict the interface behavior of dilute solute(s). The multiscale analysis showed that the presented model can recover the macroscopic transport equations derived from the Maxwell-Stefan equation for dilute solutes in three-phase systems. Compared with the analytical equation of state of solute and dynamic behavior, these results are proven to constitute a generalized framework to simulate solute distributions in three-phase flows, including compound soluble in one phase, compound adsorbed on single-interface, compound in two phases, and solute soluble in three phases. Moreover, numerical simulations of benchmark cases, such as phase decomposition, multilayered planar interfaces, and liquid lens, were performed to test the stability and efficiency of the model. Finally, the multiphase mass transfer and reaction in Janus droplet transport in a straight microchannel were well reproduced.

Apparent contact angle and contact angle hysteresis on liquid infused surfaces

We theoretically investigate the apparent contact angle and contact angle hysteresis of a droplet placed on a liquid infused surface. We show that the apparent contact angle is not uniquely defined by material parameters, but also has a dependence on the relative size between the droplet and its surrounding wetting ridge formed by the infusing liquid. We derive a closed form expression for the contact angle in the limit of vanishing wetting ridge, and compute the correction for small but finite ridge, which corresponds to an effective line tension term. We also predict contact angle hysteresis on liquid infused surfaces generated by the pinning of the contact lines by the surface corrugations. Our analytical expressions for both the apparent contact angle and contact angle hysteresis can be interpreted as 'weighted sums' between the contact angles of the infusing liquid relative to the droplet and surrounding gas phases, where the weighting coefficients are given by ratios of the fluid surface tensions.

Comparative study of the lattice Boltzmann models for Allen-Cahn and Cahn-Hilliard equations

In this paper, a comparative study of the lattice Boltzmann (LB) models for the Allen-Cahn (A-C) and Cahn-Hilliard (C-H) equations is conducted. To this end, a new LB model for the A-C equation is first proposed, where the equilibrium distribution function and the source term distribution function are delicately designed to recover the A-C equation correctly. The gradient term in this model can be computed by the nonequilibrium part of the distribution function such that the collision process can be implemented locally. Then a detailed numerical study on several classical problems is performed to give a comparison between the present model for the A-C equation and the recently developed LB model [H. Liang et al., Phys. Rev. E 89, 053320 (2014)PLEEE81539-375510.1103/PhysRevE.89.053320] for the C-H equation in terms of tracking the interface of two-phase flow. The results show that the present LB model for the A-C equation is more accurate and more stable, and also has a second-order convergence rate in space, while the convergence rate of the previous LB model for the C-H equation is only about 1.5.