Chemical-potential multiphase lattice Boltzmann method with superlarge density ratios

Theoretical model of dynamics and stability of nanobubbles on heterogeneous surfaces

Surface nanobubbles have revealed a new mechanism of gas-liquid-solid interaction at the nanoscale; however, the nanobubble evolution on real substrates is still veiled, because the experimental observation of contact line motions at the nanoscale is too difficult.

HYPOTHESIS

This study proposes a theoretical model to describe the dynamics and stability of nanobubbles on heterogeneous substrates. It simultaneously considers the diffusive equilibrium of the liquid-gas interface and the mechanical equilibrium at the contact line, and introduces a surface energy function to express the substrate's heterogeneity.

VALIDATION

The present model unifies the nanoscale stability and the microscale instability of surface bubbles. The theoretical predictions are highly consistent to the nanobubble morphology on heterogeneous surfaces observed in experiments. As the nanobubbles grow, a lower Laplace pressure leads to weaker gas adsorption, and the mechanical equilibrium can eventually revert to the classical Young-Laplace equation above microscale.

FINDINGS

The analysis results indicate that both the decrease in substrate surface energy and the increase in gas oversaturation are more conducive to the nucleation and growth of surface nanobubbles, leading to larger stable sizes. The larger surface energy barriers result in the stronger pinning, which is beneficial for achieving stability of the pinned bubbles. The present model is able to reproduce the continual behaviors of the three-phase contact line during the nanobubble evolution, e.g., "pinning, depinning, slipping and jumping" induced by the nanoscale defects.

Three-dimensional lattice Boltzmann model with self-tuning equation of state for multiphase flows

In this work, the recent lattice Boltzmann (LB) model with self-tuning equation of state (EOS) [Huang et al., Phys. Rev. E 99, 023303 (2019)2470-004510.1103/PhysRevE.99.023303] is extended to three dimensions for the simulation of multiphase flows, which is based on the standard three-dimensional 27-velocity lattice and multiple-relaxation-time collision operator. To achieve the self-tuning EOS, the equilibrium moment is devised by introducing a built-in variable, and the collision matrix is improved by introducing some velocity-dependent nondiagonal elements. Meanwhile, the additional cubic terms of velocity in recovering the Newtonian viscous stress are eliminated to enhance the numerical accuracy. For modeling multiphase flows, an attractive pairwise interaction force is introduced to mimic the long-range molecular interaction, and a consistent scheme is proposed to compensate for the ɛ^{3}-order discrete lattice effect. Thermodynamic consistency in a strict sense is established for the multiphase LB model with self-tuning EOS, and the wetting condition is also treated in a thermodynamically consistent manner. As a result, the contact angle, surface tension, and interface thickness can be independently adjusted in the present theoretical framework. Numerical tests are first performed to validate the multiphase LB model with self-tuning EOS and the theoretical analyses of bulk and surface thermodynamics. The collision of equal-sized droplets is then simulated to demonstrate the applicability and effectiveness of the present LB model for multiphase flows.

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.

Simplified wetting boundary scheme in phase-field lattice Boltzmann model for wetting phenomena on curved boundaries

In this work, a simplified wetting boundary scheme in the phase-field lattice Boltzmann model is developed for wetting phenomena on curved boundaries. The proposed scheme combines the advantages of the fluid-solid interaction scheme and geometric scheme-easy to implement (no need to interpolate the values of parameters exactly on solid boundaries and find proper characteristic vectors), the value of contact angle can be directly prescribed, and no unphysical spurious mass layer-and avoids mass leakage. Different from previous works, the values of the order parameter gradient on fluid boundary nodes are directly determined according to the geometric formulation rather than indirectly regulated through the order parameters on ghost solid nodes (i.e., ghost contact-line region). For this purpose, two numerical approaches to evaluate the order parameter gradient on fluid boundary nodes are utilized, one with the prevalent isotropic central scheme and the other with a local gradient scheme that utilizes the distribution functions. The simplified wetting boundary schemes with both numerical approaches are validated and compared through several numerical simulations. The results demonstrate that the proposed model has good ability and satisfactory accuracy to simulate wetting phenomena on curved boundaries under large density ratios.

Electrowetting lattice Boltzmann method for micro- and nano-droplet manipulations

Electrowetting has become a widely used tool for manipulating tiny amounts of liquids on surfaces. This paper proposes an electrowetting lattice Boltzmann method for manipulating micro-nano droplets. The hydrodynamics with the nonideal effect is modeled by the chemical-potential multiphase model, in which the phase transition and equilibrium are directly driven by chemical potential. For electrostatics, droplets in the micro-nano scale cannot be considered as equipotential as macroscopic droplets due to the Debye screening effect. Therefore, we linearly discretize the continuous Poisson-Boltzmann equation in a Cartesian coordinate system, and the electric potential distribution is stabilized by iterative computations. The electric potential distribution of droplets at different scales suggests that the electric field can still penetrate micro-nano droplets even with the screening effect. The accuracy of the numerical method is verified by simulating the static equilibrium of the droplet under the applied voltage, and the results show the apparent contact angles agree very well with the Lippmann-Young equation. The microscopic contact angles present some obvious deviations due to the sharp decrease of electric field strength near the three-phase contact point. These are consistent with previously reported experimental and theoretical analyses. Then, the droplet migrations on different electrode structures are simulated, and the results show that droplet speed can be stabilized more quickly due to the more uniform force on the droplet in the closed symmetric electrode structure. Finally, the electrowetting multiphase model is applied to study the lateral rebound of droplets impacting on the electrically heterogeneous surface. The electrostatic force prevents the droplets from contracting on the side which is applied voltage, resulting in the lateral rebound and transport toward the side.

Contact Angle Measurement on Curved Wetting Surfaces in Multiphase Lattice Boltzmann Method

Contact angle is an essential physical quantity that characterizes the wettability of a substrate. Although it is widely used in the studies of surface wetting, capillary phenomena, and moving contact lines, the contact angle measurements in simulations and experiments are still complicated and time-consuming. In this paper, we present an efficient scheme for the measurement of contact angle on curved wetting surfaces in lattice Boltzmann simulations. The measuring results are in excellent agreement with the theoretical predictions without considering the gravity effect. A series of simulations with various drop sizes and surface curvatures confirm that the present scheme is grid-independent. Then, the scheme is verified in gravitational environments by simulating the deformations of sessile and pendent droplets on the curved wetting surface. The numerical results are highly consistent with experimental observations and support the theoretical analysis that the microscopic contact angle is independent of gravity. Furthermore, the method utilizes only the microscopic geometry of the contact angle and does not depend on the droplet profile; therefore, it can be applied to nonaxisymmetric shapes or moving contact lines. The scheme is applied to capture the dynamic contact angle hysteresis on homogeneous or chemically heterogeneous curved surfaces. Importantly, the accurate contact angle measurement enables the dynamic mechanical analysis of moving contact lines. The present measurement is simple and efficient and can be extended to implementations in various multiphase lattice Boltzmann models.

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.

Movable and Focus-Tunable Lens Based on Electrically Controllable Liquid: A Lattice Boltzmann Study

Adjusting the focal length by changing the liquid interface of the liquid lens has become a potential method. In this paper, the lattice-Boltzmann-electrodynamic (LB-ED) method is used to numerically investigate the zooming process of a movable and focus-tunable electrowetting-on-dielectrics (EWOD) liquid lens by combining the LBM chemical potential model and the electrodynamic model. The LB method is used to solve the Navier-Stokes equation, and the Poisson-Boltzmann (PB) equation is introduced to solve the electric field distribution. The experimental results are consistent with the theoretical results of the Lippmann-Young equation. Through the simulation of a liquid lens zoom driven by EWOD, it is found that the lens changes from a convex lens to a concave lens with the voltage increases. The focal length change rate in the convex lens stage gradually increases with voltage. In the concave lens stage, the focal length change rate is opposite to that in the convex lens stage. During the zooming process, the low-viscosity liquid exhibits oscillation, and the high-viscosity liquid appears as overdamping. Additionally, methods were proposed to accelerate lens stabilization at low and high viscosities, achieving speed improvements of about 30% and 50%, respectively. Simulations of lens motion at different viscosities demonstrate that higher-viscosity liquids require higher voltages to achieve the same movement speed.

Free-Energy-Based Discrete Unified Gas Kinetic Scheme for van der Waals Fluid

The multiphase model based on free-energy theory has been experiencing long-term prosperity for its solid foundation and succinct implementation. To identify the main hindrance to developing a free-energy-based discrete unified gas-kinetic scheme (DUGKS), we introduced the classical lattice Boltzmann free-energy model into the DUGKS implemented with different flux reconstruction schemes. It is found that the force imbalance amplified by the reconstruction errors prevents the direct application of the free-energy model to the DUGKS. By coupling the well-balanced free-energy model with the DUGKS, the influences of the amplified force imbalance are entirely removed. Comparative results demonstrated a consistent performance of the well-balanced DUGKS despite the reconstruction schemes utilized. The capability of the DUGKS coupled with the well-balanced free-energy model was quantitatively validated by the coexisting density curves and Laplace's law. In the quiescent droplet test, the magnitude of spurious currents is reduced to a machine accuracy of 10-15. Aside from the excellent performance of the well-balanced DUGKS in predicting steady-state multiphase flows, the spinodal decomposition test and the droplet coalescence test revealed its stability problems in dealing with transient flows. Further improvements are required on this point.

Multiphase curved boundary condition in lattice Boltzmann method

The boundary treatment is fundamental for modeling fluid flows especially in the lattice Boltzmann method; the curved boundary conditions effectively improve the accuracy of single-phase simulations with complex-geometry boundaries. However, the conventional curved boundary conditions usually cause dramatic mass leakage or increase when they are directly used for multiphase flow simulations. We find that the principal reason for this is the absence of a nonideal effect in the curved boundary conditions, followed by a calculation error. In this paper, incorporating the nonideal effect into the linear interpolation scheme and compensating for the interpolating error, we propose a multiphase curved boundary condition to treat the wetting boundaries with complex geometries. A series of static and dynamic multiphase simulations with large density ratio verify that the present scheme is accurate and ensures mass conservation.

Alternative wetting boundary condition for the chemical-potential-based free-energy lattice Boltzmann model

The free-energy lattice Boltzmann (LB) method is a multiphase LB approach based on the thermodynamic theory. Compared with traditional free-energy LB models, which employ a nonideal thermodynamic pressure tensor, the chemical-potential-based free-energy LB model has attracted much attention in recent years as it avoids computing the thermodynamic pressure tensor and its divergence. In this paper, we propose an improved wetting boundary condition for the chemical-potential-based free-energy LB model. Different from the original wetting boundary condition in the literature, the improved wetting boundary condition utilizes a surface chemical potential that is compatible with the chemical potential of the fluid domain. Accordingly, the thermodynamic consistency of the chemical-potential-based free-energy LB model can be retained by the improved wetting boundary condition. Numerical simulations are performed for droplets resting on flat and cylindrical surfaces with different contact angles. The numerical results show that the improved wetting boundary condition yields more reasonable results and the maximum spurious velocities are found to be smaller by 2 ∼ 3 orders of magnitude than those produced by the original wetting boundary condition.

Lateral Drop Rebound on a Hydrophobic and Chemically Heterogeneous Surface

In this paper, the process of a drop rebounding from a hydrophobic and chemically heterogeneous surface is investigated using the multiphase lattice Boltzmann method. The bounce behavior of drops is dependent on the degree of hydrophobicity and heterogeneity of the surface. When the surface is homogeneous, the drop rebounds vertically and the height increases with the enhancement of the surface hydrophobicity. For the heterogeneous surface with two different hydrophobic parts, the drop rebounds laterally toward the lower hydrophobic side. Because the high hydrophobic side exerts the stronger unbalanced Young's force on the contact line compared with the low hydrophobic side, the difference of the forces results in the asymmetrical rebound. A phase diagram displays the rebound numbers of a drop impacting on the various chemically homogeneous and heterogeneous surfaces. A simply quantitative estimation is made to predict the rebound heights and flying times through the contact angles of the surface. This work promotes the understanding of the rebound mechanism of a drop impacting on a chemically heterogeneous surface and provides a guiding strategy for the precise control of the lateral behavior of rebounding drops by hydrophobic and heterogeneous surfaces.

Achieving thermodynamic consistency in a class of free-energy multiphase lattice Boltzmann models

The free-energy lattice Boltzmann (LB) model is one of the major multiphase models in the LB community. The present study is focused on a class of free-energy LB models in which the divergence of thermodynamic pressure tensor or its equivalent form expressed by the chemical potential is incorporated into the LB equation via a forcing term. Although this class of free-energy LB models may be thermodynamically consistent at the continuum level, it suffers from thermodynamic inconsistency at the discrete lattice level owing to numerical errors [Guo et al., Phys. Rev. E 83, 036707 (2010)10.1103/PhysRevE.83.036707]. The numerical error term mainly includes two parts: one comes from the discrete gradient operator and the other can be identified in a high-order Chapman-Enskog analysis. In this paper, we propose an improved scheme to eliminate the thermodynamic inconsistency of the aforementioned class of free-energy LB models. The improved scheme is constructed by modifying the equation of state of the standard LB equation, through which the discretization of ∇(ρc_{s}^{2}) is no longer involved in the force calculation and then the numerical errors can be significantly reduced. Numerical simulations are subsequently performed to validate the proposed scheme. The numerical results show that the improved scheme is capable of eliminating the thermodynamic inconsistency and can significantly reduce the spurious currents in comparison with the standard forcing-based free-energy LB model.