Implementation of contact angles in pseudopotential lattice Boltzmann simulations with curved boundaries

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.

Implementation of contact line motion based on the phase-field lattice Boltzmann method

This paper proposes a strategy to implement the free-energy-based wetting boundary condition within the phase-field lattice Boltzmann method. The greatest advantage of the proposed method is that the implementation of contact line motion can be significantly simplified while still maintaining good accuracy. For this purpose, the liquid-solid free energy is treated as a part of the chemical potential instead of the boundary condition, thus avoiding complicated interpolations with irregular geometries. Several numerical testing cases, including droplet spreading processes on the idea flat, inclined, and curved boundaries, are conducted, and the results demonstrate that the proposed method has good ability and satisfactory accuracy to simulate contact line motions.

Mesoscopic Model for Disjoining Pressure Effects in Nanoscale Thin Liquid Films and Evaporating Extended Meniscuses

Disjoining pressure effect is the key to describe contact line dynamics, micro/nanoscale liquid-vapor phase change heat transfer, and liquid transport in nanopores. In this paper, by combining a mesoscopic approach for nanoscale liquid-vapor interfacial transport and a mean-field approximation of the long-range solid-fluid molecular interaction, a mesoscopic model for the disjoining pressure effect in nanoscale thin liquid films is proposed. The capability of this model to delineate the disjoining pressure effect is validated. We demonstrate that the Hamaker constant determined from our model agrees very well with molecular dynamics (MD) simulation and that the transient evaporation/condensation mass flux predicted by this mesoscopic model is also consistent with the kinetic theory. Using this model, we investigate the characteristics of the evaporating extended meniscus in a nanochannel. The nonevaporating film region, the evaporating thin-film region, and the intrinsic meniscus region are successfully captured by our model. Our results suggest that the apparent contact angle and thickness of the nonevaporating liquid film are self-tuned according to the evaporation rate, and a higher evaporation rate results a in larger apparent contact angle and thinner nonevaporating liquid film. We also show that disjoining pressure plays a dominant role in the nonevaporating film region and suppresses the evaporation in this region, while capillary pressure dominates the intrinsic meniscus region. Strong evaporation takes place in the thin-film region, and both the disjoining pressure and capillary pressure contribute to the total pressure difference that delivers the liquid from the intrinsic meniscus region to the evaporating thin-film region, compensating for the liquid mass loss due to strong evaporation. Our work provides a new avenue for investigating thin liquid film spreading, liquid transport in nanopores, and microscopic liquid-vapor phase change heat/mass transfer mechanisms near the three-phase contact line region.

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.

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.

Coalescence-Induced Bubble Departure: Effects of Dynamic Contact Angles

Coalescence-induced bubble departure is a common phenomenon in boiling and gas evolution reactions, which has significant impacts on the heat/mass transport. In this work, we systematically investigate the effects of dynamic contact angles on the coalescence and departure processes of two equal-sized bubbles. A critical contact angle (θ_cr) of 76° is determined for an ideal surface on the basis of a surface energy analysis, beyond which the coalesced bubble does not depart from the wall. Using 3D multi-relaxation-time (MRT) lattice Boltzmann simulations, we demonstrate that the advancing contact angle mainly governs the movement of the outer side of the contact lines, and the increase of the advancing contact angle may delay or even prevent the departure of the coalesced bubble. On the other hand, the receding contact angle dominates the motion of the inner side of the contact lines, and the decrease of the receding contact angle facilitates the departure of the coalesced bubble. We identify a regime map for the coalescence-induced bubble departure with respect to the contact angles, which includes four regions: the all-departure region, the advancing contact angle dominated region, the receding contact angle dominated region, and the nondeparture region. Numerically simulated critical contact angles that separate the above-mentioned regions agree well with theoretical analyses. The results of this study will contribute to the manipulation of bubble behaviors and the optimal design of working surfaces in a variety of energy systems involving boiling and gas-evolving reaction processes.

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.

Pore-Scale Study on Convective Drying of Porous Media

In this work, a numerical model for isothermal liquid-vapor phase change (evaporation) of the two-component air-water system is proposed based on the pseudopotential lattice Boltzmann method. Through the Chapman-Enskog multiscale analysis, we show that the model can correctly recover the macroscopic governing equations of the multicomponent multiphase system with a built-in binary diffusion mechanism. The model is verified based on the two-component Stefan problem where the measured binary diffusivity is consistent with theoretical analysis. The model is then applied to convective drying of a dual-porosity porous medium at the pore scale. The simulation captures a classical transition in the drying process of porous media, from the constant rate period (CRP, first phase) showing significant capillary pumping from large to small pores, to the falling rate period (FRP, second phase) with the liquid front receding in small pores. It is found that, in the CRP, the evaporation rate increases with the inflow Reynolds number (Re), while in the FRP, the evaporation curves almost collapse at different Res. The underlying mechanism is elucidated by introducing an effective Péclet number (Pe). It is shown that convection is dominant in the CRP and diffusion in the FRP, as evidenced by Pe > 1 and Pe < 1, respectively. We also find a log-law dependence of the average evaporation rate on the inflow Re in the CRP regime. The present work provides new insights into the drying physics of porous media and its direct modeling at the pore scale.

Comparative investigation of a lattice Boltzmann boundary treatment of multiphase mass transport with heterogeneous chemical reactions

Multiphase reactive transport in porous media is an important component of many natural and engineering processes. In the present study, boundary schemes for the continuum species transport-lattice Boltzmann (CST-LB) mass transport model and the multicomponent pseudopotential model are proposed to simulate heterogeneous chemical reactions in a multiphase system. For the CST-LB model, a lattice-interface-tracking scheme for the heterogeneous chemical reaction boundary is provided. Meanwhile, a local-average virtual density boundary scheme for the multicomponent pseudopotential model is formulated based on the work of Li et al. [Li, Yu, and Luo, Phys. Rev. E 100, 053313 (2019)10.1103/PhysRevE.100.053313]. With these boundary treatments, a numerical implementation is put forward that couples the multiphase fluid flow, interfacial species transport, heterogeneous chemical reactions, and porous matrix structural evolution. A series of comparison benchmark cases are investigated to evaluate the numerical performance for different pseudopotential wetting boundary treatments, and an application case of multiphase dissolution in porous media is conducted to validate the present models' ability to solve complex problems. By applying the present LB models with reasonable boundary treatments, multiphase reactive transport in various natural or engineering scenarios can be simulated accurately.

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.

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'.

Fluid-wall interactions in pseudopotential lattice Boltzmann models

Designing proper fluid-wall interaction forces to achieve proper wetting conditions is an important area of interest in pseudopotential lattice Boltzmann models. In this paper, we propose a modified fluid-wall interaction force that applies for pseudopotential models of both single-component fluids and partially miscible multicomponent fluids, such as hydrocarbon mixtures. A reliable correlation that predicts the resulting liquid contact angle on a flat solid surface is also proposed. This correlation works well over a wide variety of pseudopotential lattice Boltzmann models and thermodynamic conditions.

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.

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.