Forcing scheme in pseudopotential lattice Boltzmann model for multiphase flows

Cavitation bubble collapse near rough solid wall is modeled by the multi-relaxation-time (MRT) pseudopotential lattice Boltzmann (LB) model. The modified forcing scheme, which can achieve LB model’s thermodynamic consistency by tuning a parameter related with the particle interaction range, is adopted to achieve desired stability and density ratio. The bubble collapse near rough solid wall was simulated by the improved MRT pseudopotential LB model. The mechanism of bubble collapse is studied by investigating the bubble profiles, pressure field and velocity field evolution. The eroding effects of collapsing bubble are analyzed in details. It is found that the process and the effect of the interaction between bubble collapse and rough solid wall are affected seriously by the geometry of solid boundary. At the same time, it demonstrates that the MRT pseudopotential LB model is a potential tool for the investigation of the interaction mechanism between the collapsing bubble and complex geometry boundary.

Simulation of two-phase fluid mixture flow in rectangular two-inlet cavity using lattice Boltzmann method

In this work, our recently developed weighted-orthogonal multiple-relaxation-time pseudopotential lattice Boltzmann (PLB) model [J. Tang et al., “Multiphase flow simulation with three-dimensional weighted-orthogonal multiple-relaxation-time pseudopotential lattice Boltzmann model,” Phys. Fluids 33, 123305 (2021)] is further extended to simulate the complex wetting phenomenon on curved surfaces at large density ratios (ρl/ρg∼1000), where a new geometrical formulation scheme is proposed to characterize the wettability of the curved boundary. Compared with the existing geometrical formulation schemes, the significant advantage of the new scheme is that the characteristic vector representing the phase interface is no longer needed, and, thus, the complex calculations induced by the characteristic vector are avoided, which significantly simplifies computations and facilitates the implementation of the geometrical formulation scheme on curved boundaries. Meanwhile, it is applicable to both two-dimensional and three-dimensional (3D) simulations and maintains the feature of setting the contact angle explicitly. Furthermore, the numerical results of four classical wetting phenomenon benchmark cases at large density ratios predicted by the present model agree well with the analytical solutions, numerical results, or experimental results in the literature. It exhibits the capability of the present model coupled with the proposed scheme to simulate the wetting phenomenon involving curved surfaces with good numerical accuracy. Note that, to the author's knowledge, this is the first time that the geometrical formulation scheme has been successfully adopted in the 3D PLB model to simulate the wetting phenomenon on curved surfaces. We believe that this work lays the foundation for further application of the PLB model to the complex wetting phenomenon.

The pseudopotential model in lattice Boltzmann method is widely used to simulate two-phase flows. Thermodynamic consistency is an important issue in the pseudopotential model. To achieve thermodynamic consistency, various force schemes have been proposed for the pseudopotential model. However, their performance under transient and moving conditions remains insufficiently validated. We evaluate existing force schemes and find that most of them suffer from interfacial velocity slip. A Chapman-Enskog analysis reveals that these schemes introduce error terms in the macroscopic equations, leading to the interfacial velocity slip. While other schemes avoid this slip, they fail to achieve thermodynamic consistency. To address this issue, we propose a novel pseudopotential model achieving thermodynamic consistency while preventing interfacial velocity slip. Our model introduces a new pseudopotential force by adding a term 2γ∇^{2}ψ∇ψ+O(∇^{5}) to the original force. Here, ψ represents the pseudopotential and the coefficient γ is used to tune the coexistence densities. The applicability of our model under moving conditions is validated through simulations of two-phase Couette flow and two-phase Poiseuille flow. Numerical results confirm that our model achieves second-order accuracy. The transient results of our model are compared with those of volume of fluid method, validating the applicability of our model under transient conditions. Finally, droplet splashing and liquid film evaporation under shear flow are also simulated to demonstrate our model's capability for simulating two-dimensional flow and interphase mass transfer.

A multicomponent multiphase (MCMP) pseudopotential lattice Boltzmann (LB) model with large liquid–gas density ratios is proposed for simulating the wetting phenomena. In the proposed model, two layers of neighboring nodes are adopted to calculate the fluid–fluid cohesion force with higher isotropy order. In addition, the different-time-step method is employed to calculate the processes of particle propagation and collision for the two fluid components with a large pseudo-particle mass contrast. It is found that the spurious current is remarkably reduced by employing the higher isotropy order calculation of the fluid–fluid cohesion force. The maximum spurious current appearing at the phase interfaces is evidently influenced by the magnitudes of fluid–fluid and fluid–solid interaction strengths, but weakly affected by the time step ratio. The density ratio analyses show that the liquid–gas density ratio is dependent on both the fluid–fluid interaction strength and the time step ratio. For the liquid–gas flow simulations without solid phase, the maximum liquid–gas density ratio achieved by the present model is higher than 1000:1. However, the obtainable maximum liquid–gas density ratio in the solid–liquid–gas system is lower. Wetting phenomena of droplets contacting smooth/rough solid surfaces and the dynamic process of liquid movement in a capillary tube are simulated to validate the proposed model in different solid–liquid–gas coexisting systems. It is shown that the simulated intrinsic contact angles of droplets on smooth surfaces are in good agreement with those predicted by the constructed LB formula that is related to Young’s equation. The apparent contact angles of droplets on rough surfaces compare reasonably well with the predictions of Cassie’s law. For the simulation of liquid movement in a capillary tube, the linear relation between the liquid–gas interface position and simulation time is observed, which is identical to the analytical prediction. The simulation results regarding the wetting phenomena of droplets on smooth/rough surfaces and the dynamic process of liquid movement in the capillary tube demonstrate the quantitative capability of the proposed model.

Based on the one-component multiphase lattice Boltzmann method, a novel solid–liquid conjugate boiling heat transfer pseudo-potential lattice Boltzmann (LB) model is tentatively proposed in this paper. By respectively introducing the physical property parameters of solids and liquids into the relaxation time τT of the temperature distribution equation, different energy transfer rates in solid, liquid, and vapor regions can be successfully predicted. After verifying the accuracy, stability, and reasonability of this model, the bubble detaching behavior and boiling heat transfer performance on the rectangular cavity structure are analyzed through setting different contact angles of the cavity surface and plane heating surface. The results show that the hydrophobic cavity surface can initialize bubble nucleation earlier and obviously increase the bubble detaching frequency because of its gas-bounding character, while the hydrophilic plane heating surface can restrict the expansion of bubbles and delay the appearance of film boiling. Moreover, for uniform wettability surfaces, the bubble detaching period varies in the quadratic equation with the surface contact angle due to the interaction of surface tension and buoyancy, and there is a minimum detaching period. While for the mixed wettability surfaces, the bubble detaching period also has a minimum value with the decrease in the contact angle the cavity surface, but the average bubble detaching diameter basically does not change with the cavity surface contact angle; moreover, the cavity surface contact angle corresponding to the minimum detaching period also increases with the increase in the plane heating surface contact angle. In addition, for the boiling heat transfer surface with cavity structure, the maximum heat flux and temperature gradient occur on the cavity surface, and the local heat flux of the hydrophobic cavity surface is higher than that of the hydrophilic cavity surface. This work will provide useful help for the further development of the conjugate boiling heat transfer LB model and clarify the mechanism of enhanced boiling heat transfer on a mixed wettability surface.

A pseudopotential multicomponent lattice Boltzmann (LB) model that can account for the real buoyancy effect is proposed to simulate the mass transfer process around a rising bubble. The density profiles at the equilibrium state are determined based on the hydrostatic condition and the zero diffusion flux condition (the balance of chemical potential). Compared with the LB models using effective buoyancy force, the proposed model has three advantages: (1) avoiding the unrealistic distribution of gas components within the bubble due to the upward effective buoyancy force, (2) removing the undesirable diffusion process due to the application of effective buoyancy force, and (3) considering the effect of the pressure gradient on the change of bubble size. In addition, Henry's law, which can be automatically recovered from the multicomponent LB equation, is re-interpreted from the perspective of the balance of chemical potential. Simulation results showed that the diffusion flux non-uniformly distributes over the surface of a rising bubble. The diffusion zone primarily occurs at the top and the lateral side of a rising bubble, whereas the diffusion transport just below the rising bubble is much less significant than its counterpart above the rising bubble. Various bubble shapes and their corresponding diffusion zones have been obtained. Moreover, the correlation between the Sherwood number and the Peclet number derived from the simulation results is consistent with those from previous numerical results. Thus, the proposed LB model is capable of conducting a quantitative analysis of the mass transfer around a rising bubble.

An improved pseudopotential lattice–Boltzmann model was proposed for simulating multiphase flow dynamics to describe fuel droplets, and its thermodynamic consistency was tested against the Peng–Robinson equation of state. The studied liquid fuels included paraffinic hydrocarbons with a different number of carbon atoms (C[Formula: see text]–C[Formula: see text]), methanol (CH[Formula: see text]OH), hydrogen (H[Formula: see text]), ammonia (NH[Formula: see text]), and water (H[Formula: see text]O). To improve accuracy and reduce the magnitude of the spurious currents, the multi-relaxation times collision operator was implemented and the forcing term was computed using the hybrid pseudopotential interaction force with an eighth-order isotropic degree. The pseudopotential lattice–Boltzmann model accurately predicted the equilibrium densities and captured satisfactorily the thermodynamic vapor-liquid coexistence curve given by the analytical solution of the Peng–Robinson equation of state for acentric factors ranging from [Formula: see text]0.22 to 0.56, keeping the maximum average error for the liquid and vapor branches below 0.8% and 3.7%, respectively. Nevertheless, Peng–Robinson was found to be insufficiently accurate to replicate the actual thermodynamic state, especially for H[Formula: see text]O and CH[Formula: see text]OH, for which the results strongly deviated from the experimental vapor-liquid equilibrium densities and reached average errors for the vapor phase of nearly 28%. Furthermore, the surface tension ([Formula: see text]) was retrieved using the multiphase pseudopotential lattice–Boltzmann results and served to verify the thermodynamic consistency of the pseudopotential lattice–Boltzmann with respect to the parachor model. Lastly, the pseudopotential lattice–Boltzmann model was also shown to predict accurately the transient behavior of oscillating droplets. Overall, the enhanced model satisfactorily predicted the properties and behavior of the substances for a wide range of conditions.

A lattice Boltzmann model for the conjugate heat transfer

The dynamics of inkjet deposition in square microcavities are investigated utilizing a three-dimensional multi-relaxation-time pseudopotential lattice Boltzmann (LB) model with large density ratios. A geometric scheme is considered within the pseudopotential LBM framework to obtain the desired contact angles. The effects of wettability, density ratios, droplet viscosity and impact velocity are explored to reveal the droplet–microcavity interactions. With the contact angles of microcavity increasing, the physical outcomes including the crown-like shape with a small round dot, circular hollow core, uniform film and convex film are identified and analyzed. At a lower density ratio ρr = 11.6, the surrounding denser gas resists the droplet recoiling flow resulting in an increasing hollow core. The appropriate higher droplet viscosity and decreasing impact velocity are preferred which could eliminate the hollow core in the recoiling phase and accelerate the inkjet deposition process straightforward. The revelation of droplet-microcavity dynamics is beneficial for optimizing inkjet deposition process and fabricating uniform OLEDs panels.

The multiphase lattice Boltzmann (LB) models based on pairwise interactions show great potential for simulating multiphase flows due to the conceptual and computational simplicity. Although the dynamics of multiphase flows are reproduced by the pairwise interaction force, the gradient of density (or effective density, i.e., pseudopotential) is implicitly involved in these models via the specialized forcing scheme or the consistent scheme for ɛ^{3}-order term. This work focuses on the calculation of density gradient in this class of multiphase LB models. Theoretical analyses are first carried out to reveal the involvement and calculation of density gradient. On the basis of a low Mach number approximation, an improved scheme is then proposed to calculate the density gradient for the recent LB model with self-tuning equation of state. Analytical and numerical calculations show that the improved scheme is more accurate and can help to reduce the numerical error when the reduced temperature is relatively low.

The present paper analyzes the typical unified lattice Boltzmann (LB) models for different convection-diffusion (CD) problems in multiphase systems. The CD problems in multiphase systems can be roughly classified into three groups: CD problems with a continuous scalar value and a continuous flux, a discontinuous scalar value and a continuous flux, a continuous scalar value and a discontinuous flux. The characteristics of the corresponding unified LB models for the three kinds of CD problems are analyzed and the equivalence between the LB models based on different perspectives or numerical schemes is revealed. Finally, a comprehensive multiphase LB model (CMLBM) capable of solving different isotropic and anisotropic CD problems in multiphase systems is proposed. Four typical CD problems in multiphase systems are calculated to validate the CMLBM; the results show that it performs well against the typical isotropic and anisotropic CD problems in multiphase systems.

Numerical analyses and optimization of tubular thermochemical heat storage reactors using axisymmetric thermal lattice Boltzmann model

In this paper, we develop an efficient and alternative lattice Boltzmann (LB) model for simulating immiscible incompressible N-phase flows (N≥2) based on the Cahn-Hilliard phase field theory. In order to facilitate the design of LB model and reduce the calculation of the gradient term, the governing equations of the N-phase system are reformulated, and they satisfy the conservation of mass, momentum and the second law of thermodynamics. In the present model, (N-1) LB equations are employed to capture the interface, and another LB equation is used to solve the Navier-Stokes (N-S) equations, where a new distribution function for the total force is delicately designed to reduce the calculation of the gradient term. The developed model is first validated by two classical benchmark problems, including the tests of static droplets and the spreading of a liquid lens, the simulation results show that the current LB model is accurate and efficient for simulating incompressible N-phase fluid flows. To further demonstrate the capability of the LB model, two numerical simulations, including dynamics of droplet collision for four fluid phases and dynamics of droplets and interfaces for five fluid phases, are performed to test the developed model. The results show that the present model can successfully handle complex interactions among N (N≥2) immiscible incompressible flows.

Lattice Boltzmann model for incompressible flows through porous media with time-fractional effects