Achieving tunable surface tension in the pseudopotential lattice Boltzmann modeling of multiphase flows

Characterizing Dynamic Contact Angle during Gas-Liquid Imbibition in Microchannels by Lattice Boltzmann Method Modeling

Spontaneous imbibition in microchannels is a critical phenomenon in various industrial applications, such as enhanced oil recovery and microfluidic systems. One of the key factors influencing the imbibition process is the dynamic wetting effect, which governs the interaction between the liquid and solid surfaces. This paper improves the original pseudopotential model for interfluid forces by coupling it with the Peng-Robinson equation of state. The model's accuracy is verified through thermodynamic consistency checks, simulations of gas-liquid interfacial tensions, and testing of static equilibrium contact angles. Following model validation, we use it to simulate spontaneous gas-liquid imbibition in microchannels and investigate dynamic contact angle evolution during the process. The results demonstrate that (1) as the microchannel width increases, inertia forces become more significant during the initial imbibition stages, leading to a greater difference between the dynamic and static contact angles. (2) A decrease in fluid-solid interaction strength results in a larger gap between dynamic and static contact angles. (3) Higher interfacial tension strengthens the capillary forces, accelerating the imbibition rate and enlarging the difference between the dynamic and static contact angles. Furthermore, the dynamic contact angle data obtained from our simulations can be used to correct the traditional Lucas-Washburn equation. The corrected equation predicts imbibition distances that closely match the simulation results.

Three-dimensional solidification modeling of various materials using the lattice Boltzmann method with an explicit enthalpy equation

Based on the mesoscopic scale, the lattice Boltzmann method (LBM) with an enthalpy-based model represented in the form of distribution functions is widely used in the liquid-solid phase transition process of energy storage materials due to its direct and relatively accurate characterization of the presence of latent heat of solidification. However, since the enthalpy distribution function itself contains the physical properties of the material, these properties are transferred along with the enthalpy distribution function during the streaming process. This leads to deviations between the enthalpy-based model when simulating the phase transition process of different materials mixed and the actual process. To address this issue, in this paper, we construct an enthalpy-based model for different types of materials. For multiple materials, various forms of enthalpy distribution functions are employed. This method still uses the form of enthalpy distribution functions for collisions and streaming processes among the same type of substance, while for heat transfer between different materials, it avoids the direct transfer of enthalpy distribution functions and instead applies a source term to the enthalpy distribution functions, characterizing the heat transfer between different materials through the energy change before and after mixing based on the temperature. To verify the accuracy of the method proposed in this paper, a detailed solidification model for two different materials is constructed using the example of water droplets solidifying in air, and the results are compared with experimental outcomes. The results of the simulation show that the model constructed in this paper is largely in line with the actual process.

Metastable and unstable hydrodynamics in multiphase lattice Boltzmann

Metastability in liquids is at the foundation of complex phase transformation dynamics such as nucleation and cavitation. Intermolecular interaction details, beyond the equation of state, and thermal hydrodynamic fluctuations play a crucial role. However, most numerical approaches suffer from a slow time and space convergence, thus hindering the convergence to the hydrodynamic limit. This work shows that the Shan-Chen lattice Boltzmann model has the unique capability of simulating the hydrodynamics of the metastable state. The structure factor of density fluctuations is theoretically obtained and numerically verified to a high precision, for all simulated wave vectors, reduced temperatures, and pressures, deep into the metastable region. Such remarkable agreement between the theory and simulations leverages the exact implementation at the lattice level of the mechanical equilibrium condition. The static structure factor is found to consistently diverge as the temperature approaches the critical point or the density approaches the spinodal line at a subcritical temperature. Theoretically predicted critical exponents are observed in both cases. Finally, the phase separation in the unstable branch follows the same pattern, i.e., the generation of interfaces with different topology, as observed in molecular dynamics simulations.

Spontaneous separation and evaporation mechanism of self-rewetting fluid droplets on chemically stripe-patterned surfaces: A lattice Boltzmann study

The evaporation characteristics of self-rewetting fluids have attracted much attention in recent years. However, the evaporation dynamics as well as the underlying evaporation mechanism of self-rewetting fluid droplets has not been well understood. In this paper, we numerically investigate the evaporation performance and the dynamic behavior of self-rewetting fluid droplets on chemically patterned surfaces using a thermal multiphase lattice Boltzmann model with liquid-vapor phase change. First, it is shown that a self-rewetting fluid droplet can spontaneously separate into two droplets during its evaporation on a hydrophilic surface with a hydrophobic stripe, while no separation occurs during the evaporation of a conventional fluid droplet. The positive surface tension gradient of the self-rewetting fluid is found to play an important role in the spontaneous separation of the self-rewetting fluid droplet during the evaporation. Meanwhile, the separation behavior of the self-rewetting fluid droplet can effectively increase the length of the triple-phase contact line, which leads to a significant increase in the evaporation rate as compared with that of a conventional fluid droplet. Moreover, by investigating the evaporation performance of self-rewetting fluid droplets on chemically stripe-patterned surfaces with different values of the widths of the hydrophilic and hydrophobic stripes, it is found that the stripe width and the initial location of the droplet significantly affect the dynamic behavior and the evaporation efficiency of the self-rewetting fluid droplet. For different relative positions between the droplet and the stripes, the droplet may spontaneously separate into two or three droplets and achieve much better evaporation efficiency when the stripe width is within an optimal range.

Unified lattice Boltzmann method with improved schemes for multiphase flow simulation: Application to droplet dynamics under realistic conditions

As a powerful mesoscale approach, the lattice Boltzmann method (LBM) has been widely used for the numerical study of complex multiphase flows. Recently, Luo et al. [Philos. Trans. R. Soc. A: Math. Phys. Eng. Sci. 379, 20200397 (2021)10.1098/rsta.2020.0397] proposed a unified lattice Boltzmann method (ULBM) to integrate the widely used lattice Boltzmann collision operators into a unified framework. In this study, we incorporate additional features into this ULBM in order to simulate multiphase flow under realistic conditions. A nonorthogonal moment set [Fei et al., Phys. Rev. E 97, 053309 (2018)10.1103/PhysRevE.97.053309] and the entropic-multi-relaxation-time (KBC) lattice Boltzmann model are used to construct the collision operator. An extended combined pseudopotential model is proposed to realize multiphase flow simulation at high-density ratio with tunable surface tension over a wide range. The numerical results indicate that the improved ULBM can significantly decrease the spurious velocities and adjust the surface tension without appreciably changing the density ratio. The ULBM is validated through reproducing various droplet dynamics experiments, such as binary droplet collision and droplet impingement on superhydrophobic surfaces. Finally, the extended ULBM is applied to complex droplet dynamics, including droplet pancake bouncing and droplet splashing. The maximum Weber number and Reynolds number in the simulation reach 800 and 7200, respectively, at a density ratio of 1000. The study demonstrates the generality and versatility of ULBM for incorporating schemes to tackle challenging multiphase problems.

Improved three-dimensional thermal multiphase lattice Boltzmann model for liquid-vapor phase change

Modeling liquid-vapor phase change using the lattice Boltzmann (LB) method has attracted significant attention in recent years. In this paper, we propose an improved three-dimensional thermal multiphase LB model for simulating liquid-vapor phase change. The proposed model has the following features. First, it is still within the framework of the thermal LB method using a temperature distribution function and therefore retains the fundamental advantages of the thermal LB method. Second, in the existing thermal LB models for liquid-vapor phase change, the finite-difference computations of the gradient terms ∇·u and ∇T usually require special treatment at boundary nodes, while in the proposed thermal LB model these two terms are calculated locally. Moreover, in some of the existing thermal LB models, the error term ∂_{t_{0}}(Tu) is eliminated by adding local correction terms to the collision process in the moment space, which causes these thermal LB models to be limited to the D2Q9 lattice in two dimensions and the D3Q15 or D3Q19 lattice in three dimensions. Conversely, the proposed model does not suffer from such an error term and therefore the thermal LB equation can be constructed on the D3Q7 lattice, which simplifies the model and improves the computational efficiency. Numerical simulations are carried out to validate the accuracy and efficiency of the proposed thermal multiphase LB model for simulating liquid-vapor phase change.

Phase-field lattice Boltzmann model for two-phase flows with large density ratio

In this work, a lattice Boltzmann (LB) model based on the phase-field method is proposed for simulating large density ratio two-phase flows. An improved multiple-relaxation-time (MRT) LB equation is first developed to solve the conserved Allen-Cahn (AC) equation. By utilizing a nondiagonal relaxation matrix and modifying the equilibrium distribution function and discrete source term, the conserved AC equation can be correctly recovered by the proposed MRT LB equation with no deviation term. Therefore, the calculations of the temporal derivative term in the previous LB models are successfully avoided. Numerical tests demonstrate that satisfactory accuracy can be achieved by the present model to solve the conserved AC equation. What is more, the discrete force term of the MRT LB equation for the incompressible Navier-Stokes equations is also simplified and modified in the present work. An alternative scheme to calculate the gradient terms of the order parameter involved in the discrete force term through the nonequilibrium part of the distribution function is also developed. To validate the ability of the present LB model for simulating large density ratio two-phase flows, series of benchmarks, including two-phase Poiseuille flow, droplet impacting on thin liquid film, and planar Taylor bubble are simulated. It is found that the results predicted by the present LB model agree well with the analytical, numerical, and experimental results.

Shaping the equation of state to improve numerical accuracy and stability of the pseudopotential lattice Boltzmann method

It has recently been shown that altering the shape of the metastable and unstable branches of an equation of state (EOS) can substantially improve the numerical accuracy of liquid and vapor densities in the pseudopotential lattice Boltzmann method [Peng et al., Phys. Rev. E 101, 063309 (2020)2470-004510.1103/PhysRevE.101.063309]. We found that this approach reduces stability of the method in nonequilibrium conditions and is unstable for bubbles at low reduced temperatures. Here we present an improved method for altering the shape of the metastable and unstable branches of the EOS which remains stable for both equilibrium and nonequilibrium situations and has no issues with bubbles. We also performed a detailed study of the stability of the methods for a droplet impact on a liquid film for reduced temperatures down to 0.35 with Reynolds number of 300. Our approach remained stable for a density ratio of up to 3.38×10^{4}.

Simulations of surface acoustic wave interactions on a sessile droplet using a three-dimensional multiphase lattice Boltzmann model

This study reports the development of a three-dimensional numerical model for acoustic interactions with a microscale sessile droplet under surface acoustic wave (SAW) excitation using the lattice Boltzmann method (LBM). We first validate the model before SAW interactions are added. The results demonstrate good agreement with the analytical results for thermodynamic consistency, Laplace law, static contact angle on a flat surface, and droplet oscillation. We then investigate SAW interactions on the droplet, with resonant frequencies ranging 61.7-250.1 MHz. According to our findings, an increase in wave amplitude elicits an increase in streaming velocity inside the droplet, causing internal mixing, and further increase in wave amplitude leads to pumping and jetting. The boundaries of wave amplitude at various resonant frequencies are predicted for mixing, pumping, and jetting modes. The modeling predictions on the roles of forces (SAW, interfacial tension, inertia, and viscosity) on the dynamics of mixing, pumping, and jetting of a droplet are in good agreement with observations and experimental data. The model is further applied to investigate the effects of SAW substrate surface wettability, viscosity ratio, and interfacial tension on SAW actuation onto the droplet. This work demonstrates the capability of the LBM in the investigation of acoustic wave interactions between SAW and a liquid medium.

A unified lattice Boltzmann model and application to multiphase flows

In this work, we develop a unified lattice Boltzmann model (ULBM) framework that can seamlessly integrate the widely used lattice Boltzmann collision operators, including the Bhatnagar-Gross-Krook or single-relation-time, multiple-relaxation-time, central-moment or cascaded lattice Boltzmann method and multiple entropic operators (KBC). Such a framework clarifies the relations among the existing collision operators and greatly facilitates model comparison and development as well as coding. Importantly, any LB model or treatment constructed for a specific collision operator could be easily adopted by other operators. We demonstrate the flexibility and power of the ULBM framework through three multiphase flow problems: the rheology of an emulsion, splashing of a droplet on a liquid film and dynamics of pool boiling. Further exploration of ULBM for a wide variety of phenomena would be both realistic and beneficial, making the LBM more accessible to non-specialists. This article is part of the theme issue 'Progress in mesoscale methods for fluid dynamics simulation'.

Modified phase-field-based lattice Boltzmann model for incompressible multiphase flows

Based on the phase-field theory, a multiple-relaxation-time (MRT) lattice Boltzmann model is proposed for the immiscible multiphase fluids. In this model, the local Allen-Chan equation is chosen as the target equation to capture the phase interface. Unlike previous MRT schemes, an off-diagonal relaxation matrix is adopted in the present model so that the target phase-field equation can be recovered exactly without any artificial terms. To check the necessity of removing those artificial terms, comparative studies were carried out among different MRT schemes with or without correction. Results show that the artificial terms can be neglected at low March number but will cause unphysical diffusion or interface undulation instability for the relatively large March number cases. The present modified model shows superiority in reducing numerical errors by adjusting the free parameters. As the interface transport coupled to the fluid flow, a pressure-evolution lattice Boltzmann equation is adopted for hydrodynamic properties. Several benchmark cases for multiphase flow were conducted to test the validity of the present model, including the static drop test, Rayleigh-Taylor instability, and single rising bubble test. For the rising bubble simulation at high density ratios, bubble dynamics obtained by the present modified MRT lattice Boltzmann model agree well with those obtained by the FEM-based level set and FEM-based phase-field models.

Structure and isotropy of lattice pressure tensors for multirange potentials

We systematically analyze the tensorial structure of the lattice pressure tensors for a class of multiphase lattice Boltzmann models (LBM) with multirange interactions. Due to lattice discrete effects, we show that the built-in isotropy properties of the lattice interaction forces are not necessarily mirrored in the corresponding lattice pressure tensor. This finding opens a different perspective for constructing forcing schemes, achieving the desired isotropy in the lattice pressure tensors via a suitable choice of multirange potentials. As an immediate application, the obtained LBM forcing schemes are tested via numerical simulations of nonideal equilibrium interfaces and are shown to yield weaker and less spatially extended spurious currents with respect to forcing schemes obtained by forcing isotropy requirements only. From a general perspective, the proposed analysis yields an approach for implementing forcing symmetries, never explored so far in the framework of the Shan-Chen method for LBM. We argue this will be beneficial for future studies of nonideal interfaces.

Lattice Boltzmann Modeling of Drying of Porous Media Considering Contact Angle Hysteresis

Drying of porous media is governed by a combination of evaporation and movement of the liquid phase within the porous structure. Contact angle hysteresis induced by surface roughness is shown to influence multi-phase flows, such as contact line motion of droplet, phase distribution during drainage and coffee ring formed after droplet drying in constant contact radius mode. However, the influence of contact angle hysteresis on liquid drying in porous media is still an unanswered question. Lattice Boltzmann model (LBM) is an advanced numerical approach increasingly used to study phase change problems including drying. In this paper, based on a geometric formulation scheme to prescribe contact angle, we implement a contact angle hysteresis model within the framework of a two-phase pseudopotential LBM. The capability and accuracy of prescribing and automatically measuring contact angles over a large range are tested and validated by simulating droplets sitting on flat and curved surfaces. Afterward, the proposed contact angle hysteresis model is validated by modeling droplet drying on flat and curved surfaces. Then, drying of two connected capillary tubes is studied, considering the influence of different contact angle hysteresis ranges on drying dynamics. Finally, the model is applied to study drying of a dual-porosity porous medium, where phase distribution and drying rate are compared with and without contact angle hysteresis. The proposed model is shown to be capable of dealing with different contact angle hysteresis ranges accurately and of capturing the physical mechanisms during drying in different porous media including flat and curved geometries.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s11242-021-01644-9.

Improved thermal multiple-relaxation-time lattice Boltzmann model for liquid-vapor phase change

In this paper, an improved thermal multiple-relaxation-time (MRT) lattice Boltzmann (LB) model is proposed for simulating liquid-vapor phase change. A temperature equation is first derived for liquid-vapor phase change, where the latent heat of vaporization is decoupled with the equation of state. Therefore, the latent heat of vaporization can be arbitrarily specified in practice, which significantly improves the flexibility of the present LB model for liquid-vapor phase change. The Laplacian term of temperature is avoided in the proposed temperature equation and the gradient term of temperature is calculated through a local scheme. To solve the temperature equation accurately and efficiently, an improved MRT LB equation with nondiagonal relaxation matrix is developed. The implicit calculation of the temperature, caused by the source term and encountered in previous works, is avoided by approximating the source term with its value at the previous time step. As demonstrated by numerical tests, the results by the present LB model agree well with analytical results, experimental results, or the results by the finite difference method where the fourth-order Runge-Kutta method is employed to implement the discretization of time.

Shape Optimization of a Microhole Surface for Control of Droplet Wettability via the Lattice Boltzmann Method and Response Surface Methodology

The chief aim is to explore the wetting state on a microhole surface and to optimize the shape parameters of a microhole surface. A two-dimensional pseudopotential model was established, and the effects of shapes on the wetting behavior were explored. The shape parameters were optimized via the response surface methodology. The results reveal that the microhole surface can achieve a superhydrophobic state. When the diameter varies from 25 to 200 μm, the droplet is gradually lifted. However, when the diameter of the microhole is too large, the contact angle decreases rapidly. When the microhole diameter increases, relative radii of the x- and y-directions exhibit increasing trends. With the increase of the spacing, the gaps between the microholes are gradually filled with the droplet. When spacing increases, relative radii of x- and y-directions exhibit decreasing trends. The largest contact angle of 171.246° at the diameter of 76 μm and the spacing of 48 μm is observed.

Lattice Boltzmann Solver for Multiphase Flows: Application to High Weber and Reynolds Numbers

The lattice Boltzmann method, now widely used for a variety of applications, has also been extended to model multiphase flows through different formulations. While already applied to many different configurations in low Weber and Reynolds number regimes, applications to higher Weber/Reynolds numbers or larger density/viscosity ratios are still the topic of active research. In this study, through a combination of a decoupled phase-field formulation—the conservative Allen–Cahn equation—and a cumulant-based collision operator for a low-Mach pressure-based flow solver, we present an algorithm that can be used for higher Reynolds/Weber numbers. The algorithm was validated through a variety of test cases, starting with the Rayleigh–Taylor instability in both 2D and 3D, followed by the impact of a droplet on a liquid sheet. In all simulations, the solver correctly captured the flow dynamics andmatched reference results very well. As the final test case, the solver was used to model droplet splashing on a thin liquid sheet in 3D with a density ratio of 1000 and kinematic viscosity ratio of 15, matching the water/air system at We = 8000 and Re = 1000. Results showed that the solver correctly captured the fingering instabilities at the crown rim and their subsequent breakup, in agreement with experimental and numerical observations reported in the literature.

Three-dimensional pseudopotential lattice Boltzmann model for multiphase flows at high density ratio

In this study, we extend the pseudopotential lattice Boltzmann model proposed by Huang and Wu [J. Comput. Phys. 327, 121 (2016)10.1016/j.jcp.2016.09.030] to a three-dimensional model for practical simulations of multiphase flows with high density ratio. In this model, an additional source term is introduced into the evolution function, and the performed high-order Chapman-Enskog analysis demonstrates that the Navier-Stokes equations with accurate pressure tensor are recovered. Also, an alternative geometric formulation is developed to obtain various contact angles and an iteration scheme is involved in the initialization to improve the stability of the model. Theoretical and numerical investigations both validate that the thermodynamic consistency and tuning surface tension independently of density ratio is achieved through varying the two free parameters in the source term. Numerical simulations of droplet wetting indicate that a large degree range of contact angles can be precisely realized with the implementation of the wetting boundary scheme. Further dynamic examinations of droplet impingement on a thin film and a dry surface also verify the stability and capability of the proposed pseudopotential lattice Boltzmann model.

Force approach for the pseudopotential lattice Boltzmann method

One attractive feature of the original pseudopotential method consists on its simplicity of adding a force dependent on a nearest-neighbor potential function. In order to improve the method, regarding thermodynamic consistency and control of surface tension, different approaches were developed in the literature, such as multirange interactions potential and modified forcing schemes. In this work, a strategy to combine these enhancements with an appropriate interaction force field using only nearest-neighbor interactions is devised, starting from the desired pressure tensor. The final step of our procedure is implementing this external force by using the classical Guo forcing scheme. Numerical tests regarding static and dynamic flow conditions were performed. Static tests showed that current procedure is suitable to control the surface tension and phase densities. Based on thermodynamic principles, it is devised a solution for phase densities in a droplet, which states explicitly dependence on the surface tension and interface curvature. A comparison with numerical results suggest a physical inconsistency in the pseudopotential method. This fact is not commonly discussed in the literature, since most of studies are limited to the Maxwell equal area rule. However, this inconsistency is shown to be dependent on the equation of state (EOS), and its effects can be mitigated by an appropriate choice of Carnahan-Starling EOS parameters. Also, a droplet oscillation test was performed, and the most divergent solution under certain flow conditions deviated 7.5% from the expected analytical result. At the end, a droplet impact test against a solid wall was performed to verify the method stability, and it was possible to reach stable simulation results with density ratio of almost 2400 and Reynolds number of Re=373. The observed results corroborate that the proposed method is able to replicate the desired macroscopic multiphase behavior.

Multipseudopotential interaction models for thermal lattice Boltzmann method simulations

In this work, in the first instance, the multipseudopotential interaction (MPI) model's capabilities are extended for hydrodynamic simulations. This is achieved by combining MPI with the multiple-relaxation-time collision operator and with surface tension modification methods. A method of approaching thermodynamic consistency is also proposed, which consists of splitting the ɛ_{j} term into separate terms. One of these terms is used in the calculation of the interparticle force, and the second one is used in the forcing scheme. Secondly, MPI is combined with thermal models in order to simulate droplet evaporation and bubble nucleation in pool boiling. Thermal coupling is implemented using a double distribution function thermal model and a hybrid thermal model. It is found that MPI thermal models obey the D^{2}-law closely for droplet evaporation. MPI is also found to correctly simulate bubble nucleation and departure from the heating element during nucleate pool boiling. It can be suggested that MPI thermal models are comparatively better suited to thermal simulations at low reduced temperatures than single pseudopotential interaction models, although such cases remain very challenging. Droplet evaporation simulations are carried out at a reduced temperature (T_{r}) of 0.6 by setting the parameters in the Peng-Robinson equation of state to a=1/6272 and b=1/168.

Simulation of Boiling Heat Transfer at Different Reduced Temperatures with an Improved Pseudopotential Lattice Boltzmann Method

The pseudopotential Lattice Boltzmann Method has attracted much attention in the recent years for the simulation of boiling heat transfer. Many studies have been published recently for the simulation of the bubble cycle (nucleation, growth and departure from a heated surface). This paper puts forward two-dimensional simulations of bubble nucleation, growth and departure using an improved pseudopotential Lattice Boltzmann Model from the literature at different reduced temperatures, Tr=0.76 and Tr=0.86. Two different models using the Bhatnagar–Gross–Krook (BGK) and the Multiple-Relaxation-Time (MRT) collision operators with appropriate forcing schemes are used. The results for pool boiling show that the bubbles exhibit axial symmetry during growth and departure. Numerical results of departure diameter and release period for pool boiling are compared against empirical correlations from the literature by varying the gravitational acceleration. Reasonable agreement is observed. Nucleate boiling trends with heat flux are also captured by the simulations. Numerical results of flow boiling simulations are compared by varying the Reynolds number for both reduced temperatures with the MRT model. It was found that the departure diamenter and release period decreases with the increase of the Reynolds number. These results are a direct effect of the drag force. Proper conclusions are commented at the end of the paper.

Attainment of rigorous thermodynamic consistency and surface tension in single-component pseudopotential lattice Boltzmann models via a customized equation of state

The lack of thermodynamic consistency is a well-recognized problem in the single-component pseudopotential lattice Boltzmann models which prevents them from replicating accurate liquid and vapor phase densities; i.e., current models remain unable to exactly match coexisting density values predicted by the associated thermodynamic model. Most of the previous efforts had attempted to solve this problem by introducing tuning parameters, whose determination required empirical trial and error until acceptable thermodynamic consistency was achieved. In this study, we show that the problem can be alternatively solved by properly designing customized equations of state (EOSs) that replace any cubic EOS of choice during the computation of effective mass used in Shan-Chen forces. A two-parameter cubic-shaped customized EOS is introduced. Contrary to previous efforts, customization parameters in the new EOS are nonempirical and are rather derived from solving the integral mechanical stability equation, which neglects the need for any type of tuning for the attainment of rigorous thermodynamic consistency. The proposed approach reduces the errors of the coexisting densities and saturated pressure in the simulation to a maximum of 0.01% within the liquid-vapor density ratio range from O(1) to O(10^{4}), which had not been achieved in any of the previous tuning-based efforts. A straightforward way for achieving the desired surface tension via the customized EOS is also provided.

Effects of hysteresis window on contact angle hysteresis behaviour at large Bond number

Contact angle hysteresis, defined as the difference between advancing and receding contact angles, is an important phenomenon in multiphase flow on a wetting surface. In this study, a modified pseudo-potential lattice Boltzmann (LB) multiphase model with tunable surface tension is proposed, which is further coupled with the geometrical formulation contact angle scheme to investigate the motion of droplets invoking the contact angle hysteresis. We focus on the dynamic behaviour of droplets driven by a body force at the Bond number ranging from 1 to 6, which is defined as the ratio of the body force to the capillary force. The droplet morphology change is examined by varying (i) the Bond number and (ii) the hysteresis window. Results show the droplet morphology evolution can be classified into different stages, including stretch, relaxation, and equilibrium. The droplet oscillation phenomenon at large Bond numbers at the equilibrium stage is observed for the first time. In addition, it is found that such oscillation can lead to the breakup and/or coalescence of droplets when the surface waves spread on the top of the droplet. Furthermore, there is slight oscillation of the normalized length, width and height at the equilibrium stage for the neutral hysteresis window while more dramatic oscillation will appear for the hydrophobic hysteresis window.

Hybridized method of pseudopotential lattice Boltzmann and cubic-plus-association equation of state assesses thermodynamic characteristics of associating fluids

It is crucial to properly describe the associating fluids in terms of phase equilibrium behaviors, which are needed for design, operation, and optimization of various chemical and energy processes. Pseudopotential lattice Boltzmann method (LBM) appears to be a reliable and efficient approach to study thermodynamic behaviors and phase transition of complex fluid systems. However, when cubic equations of state (EOSs) are incorporated into single-component multiphase LBM, simulation results are not well matched with experimental data. This study presents the utilization of cubic-plus-association (CPA) EOS in the LBM structure to obtain more accurate modeling results for associating fluids. An approach based on the global search optimization algorithm is introduced to find the optimal association parameters of CPA EOS for water and primary alcohols in the lattice units. The thermodynamic consistency is verified by the Maxwell construction and is also improved by the forcing scheme of [Q. Li, K. H. Luo, and X. J. Li, Phys. Rev. E 86, 016709 (2012)10.1103/PhysRevE.86.016709]. The spurious velocity is reduced with increasing isotropy in the gradient operator. Furthermore, an extended version of CPA EOS is introduced, which increases the system stability at low reduced temperatures. There is a very good match between the LBM results and experimental data, confirming the reliability of the model developed in the present study. The introduced approach has potential to be employed for simulating transport phenomena and interfacial characteristics of associating fluids in porous systems.

Hydrodynamic behavior of the pseudopotential lattice Boltzmann method for interfacial flows

The lattice Boltzmann method (LBM) is routinely employed in the simulation of complex multiphase flows comprising bulk phases separated by nonideal interfaces. The LBM is intrinsically mesoscale with a hydrodynamic equivalence popularly set by the Chapman-Enskog analysis, requiring that fields slowly vary in space and time. The latter assumptions become questionable close to interfaces where the method is also known to be affected by spurious nonhydrodynamical contributions. This calls for quantitative hydrodynamical checks. In this paper, we analyze the hydrodynamic behavior of the LBM pseudopotential models for the problem of the breakup of a liquid ligament triggered by the Plateau-Rayleigh instability. Simulations are performed at fixed interface thickness, while increasing the ligament radius, i.e., in the "sharp interface" limit. The influence of different LBM collision operators is also assessed. We find that different distributions of spurious currents along the interface may change the outcome of the pseudopotential model simulations quite sensibly, which suggests that a proper fine-tuning of pseudopotential models in time-dependent problems is needed before the utilization in concrete applications. Taken all together, we argue that the results of the proposed paper provide a valuable insight for engineering pseudopotential model applications involving the hydrodynamics of liquid jets.

Eliminating cubic terms in the pseudopotential lattice Boltzmann model for multiphase flow

It is well recognized that there exist additional cubic terms of velocity in the lattice Boltzmann (LB) model based on the standard lattice. In this work, elimination of these cubic terms in the pseudopotential LB model for multiphase flow is investigated, where the force term and density gradient are considered. By retaining high-order (≥3) Hermite terms in the equilibrium distribution function and the discrete force term, as well as introducing correction terms in the LB equation, the additional cubic terms of velocity are entirely eliminated. With this technique, the computational simplicity of the pseudopotential LB model is well maintained. Numerical tests, including stationary and moving flat and circular interface problems, are carried out to show the effects of such cubic terms on the simulation of multiphase flow. It is found that the elimination of additional cubic terms is beneficial to reduce the numerical error, especially when the velocity is relatively large. Numerical results also suggest that these cubic terms mainly take effect in the interfacial region and that the density-gradient-related cubic terms are more important than the other cubic terms for multiphase flow.

Analysis of force treatment in the pseudopotential lattice Boltzmann equation method

In this paper, different force treatments are analyzed in detail for a pseudopotential lattice Boltzmann equation (LBE), and the contribution of third-order error terms to pressure tensor with a force scheme is analyzed by a higher-order Chapman-Enskog expansion technique. From the theoretical analysis, the performance of the original force treatment of Shan-Chen (SC), Ladd, Guo et al., and the exact difference method (EDM) are ɛ_{Ladd}<ɛ_{Guo}<ɛ_{EDM}≤ɛ_{SC} with the relaxation time τ≥1, while ɛ_{Ladd}<ɛ_{Guo}<ɛ_{SC}<ɛ_{EDM} with τ<1; here ɛ is a parameter related to the mechanical stability and the subscripts are the corresponding force scheme. To be consistent with the thermodynamic theory, a force term is introduced to modify the coefficients in the pressure tensor. Some numerical simulations are conducted to show that the predictions of modified force treatment of the pseudopotential LBE are all in good agreement with the analytical solution and other predictions.

Cascaded lattice Boltzmann method with improved forcing scheme for large-density-ratio multiphase flow at high Reynolds and Weber numbers

A recently developed forcing scheme has allowed the pseudopotential multiphase lattice Boltzmann method to correctly reproduce coexistence curves, while expanding its range to lower surface tensions and arbitrarily high density ratios [Lycett-Brown and Luo, Phys. Rev. E 91, 023305 (2015)PLEEE81539-375510.1103/PhysRevE.91.023305]. Here, a third-order Chapman-Enskog analysis is used to extend this result from the single-relaxation-time collision operator, to a multiple-relaxation-time cascaded collision operator, whose additional relaxation rates allow a significant increase in stability. Numerical results confirm that the proposed scheme enables almost independent control of density ratio, surface tension, interface width, viscosity, and the additional relaxation rates of the cascaded collision operator. This allows simulation of large density ratio flows at simultaneously high Reynolds and Weber numbers, which is demonstrated through binary collisions of water droplets in air (with density ratio up to 1000, Reynolds number 6200 and Weber number 440). This model represents a significant improvement in multiphase flow simulation by the pseudopotential lattice Boltzmann method in which real-world parameters are finally achievable.

Pinning-Depinning Mechanism of the Contact Line during Evaporation on Chemically Patterned Surfaces: A Lattice Boltzmann Study

In this paper, the pinning and depinning mechanism of the contact line during droplet evaporation on chemically stripe-patterned surfaces is numerically investigated using a thermal multiphase lattice Boltzmann (LB) model with liquid-vapor phase change. A local force balance in the context of diffuse interfaces is introduced to explain the equilibrium states of droplets on chemically patterned surfaces. It is shown that when the contact line is pinned on a hydrophobic-hydrophilic boundary, different contact angles can be interpreted as the variation of the length of the contact line occupied by each component. The stick-slip-jump behavior of evaporating droplets on chemically patterned surfaces is well captured by the LB simulations. Particularly, a slow movement of the contact line is clearly observed during the stick (pinning) mode, which shows that the pinning of the contact line during droplet evaporation on chemically stripe-patterned surfaces is actually a dynamic pinning process and the dynamic equilibrium is achieved by the self-adjustment of the contact lines occupied by each component. Moreover, it is shown that when the surface tension varies with the temperature, the Marangoni effect has an important influence on the depinning of the contact line, which occurs when the horizontal component (toward the center of the droplet) of the force caused by the Marangoni stress overcomes the unbalanced Young's force toward the outside.

Contact angle adjustment in equation-of-state-based pseudopotential model

The single component pseudopotential lattice Boltzmann model has been widely applied in multiphase simulation due to its simplicity and stability. In many studies, it has been claimed that this model can be stable for density ratios larger than 1000. However, the application of the model is still limited to small density ratios when the contact angle is considered. The reason is that the original contact angle adjustment method influences the stability of the model. Moreover, simulation results in the present work show that, by applying the original contact angle adjustment method, the density distribution near the wall is artificially changed, and the contact angle is dependent on the surface tension. Hence, it is very inconvenient to apply this method with a fixed contact angle, and the accuracy of the model cannot be guaranteed. To solve these problems, a contact angle adjustment method based on the geometry analysis is proposed and numerically compared with the original method. Simulation results show that, with our contact angle adjustment method, the stability of the model is highly improved when the density ratio is relatively large, and it is independent of the surface tension.

Multipseudopotential interaction: A consistent study of cubic equations of state in lattice Boltzmann models

A method is developed to analytically and consistently implement cubic equations of state into the recently proposed multipseudopotential interaction (MPI) scheme in the class of two-phase lattice Boltzmann (LB) models [S. Khajepor, J. Wen, and B. Chen, Phys. Rev. E 91, 023301 (2015)]10.1103/PhysRevE.91.023301. An MPI forcing term is applied to reduce the constraints on the mathematical shape of the thermodynamically consistent pseudopotentials; this allows the parameters of the MPI forces to be determined analytically without the need of curve fitting or trial and error methods. Attraction and repulsion parts of equations of state (EOSs), representing underlying molecular interactions, are modeled by individual pseudopotentials. Four EOSs, van der Waals, Carnahan-Starling, Peng-Robinson, and Soave-Redlich-Kwong, are investigated and the results show that the developed MPI-LB system can satisfactorily recover the thermodynamic states of interest. The phase interface is predicted analytically and controlled via EOS parameters independently and its effect on the vapor-liquid equilibrium system is studied. The scheme is highly stable to very high density ratios and the accuracy of the results can be enhanced by increasing the interface resolution. The MPI drop is evaluated with regard to surface tension, spurious velocities, isotropy, dynamic behavior, and the stability dependence on the relaxation time.

Lattice Boltzmann modeling of self-propelled Leidenfrost droplets on ratchet surfaces

In this paper, the self-propelled motion of Leidenfrost droplets on ratchet surfaces is numerically investigated using a thermal multiphase lattice Boltzmann model with liquid-vapor phase change. The capability of the model for simulating evaporation is validated via the D(2) law. Using the model, we first study the performances of Leidenfrost droplets on horizontal ratchet surfaces. It is numerically shown that the motion of self-propelled Leidenfrost droplets on ratchet surfaces is owing to the asymmetry of the ratchets and the vapor flows beneath the droplets. It is found that the Leidenfrost droplets move in the direction toward the slowly inclined side from the ratchet peaks, which agrees with the direction of droplet motion in experiments [Linke et al., Phys. Rev. Lett., 2006, 96, 154502]. Moreover, the influences of the ratchet aspect ratio are investigated. For the considered ratchet surfaces, a critical value of the ratchet aspect ratio is approximately found, which corresponds to the maximum droplet moving velocity. Furthermore, the processes that the Leidenfrost droplets climb uphill on inclined ratchet surfaces are also studied. Numerical results show that the maximum inclination angle at which a Leidenfrost droplet can still climb uphill successfully is affected by the initial radius of the droplet.

Improved forcing scheme in pseudopotential lattice Boltzmann methods for multiphase flow at arbitrarily high density ratios

The pseudopotential lattice Boltzmann method has been widely used to simulate many multiphase flow applications. However, there still exist problems with reproducing realistic values of density ratio and surface tension. In this study, a higher-order analysis of a general forcing term is derived. A forcing scheme is then constructed for the pseudopotential method that is able to accurately reproduce the full range of coexistence curves. As a result, multiphase flow of arbitrarily high density ratios independent of the surface tension can be simulated. Furthermore, the interface width can be tuned to allow for grid refinement and systematic error reduction. Numerical results confirm that the proposed scheme enables independent control of density ratio, surface tension, and interface width simultaneously.

Effect of the forcing term in the pseudopotential lattice Boltzmann modeling of thermal flows

The pseudopotential lattice Boltzmann (LB) model is a popular model in the LB community for simulating multiphase flows. Recently, several thermal LB models, which are based on the pseudopotential LB model and constructed within the framework of the double-distribution-function LB method, were proposed to simulate thermal multiphase flows [G. Házi and A. Márkus, Phys. Rev. E 77, 026305 (2008); L. Biferale, P. Perlekar, M. Sbragaglia, and F. Toschi, Phys. Rev. Lett. 108, 104502 (2012); S. Gong and P. Cheng, Int. J. Heat Mass Transfer 55, 4923 (2012); M. R. Kamali et al., Phys. Rev. E 88, 033302 (2013)]. The objective of the present paper is to show that the effect of the forcing term on the temperature equation must be eliminated in the pseudopotential LB modeling of thermal flows. First, the effect of the forcing term on the temperature equation is shown via the Chapman-Enskog analysis. For comparison, alternative treatments that are free from the forcing-term effect are provided. Subsequently, numerical investigations are performed for two benchmark tests. The numerical results clearly show that the existence of the forcing-term effect will lead to significant numerical errors in the pseudopotential LB modeling of thermal flows.