Study on C–S and P–R EOS in pseudo-potential lattice Boltzmann model for two-phase flows

2017 ◽  
Vol 28 (09) ◽  
pp. 1750120 ◽  
Author(s):  
Yong Peng ◽  
Yun Fei Mao ◽  
Bo Wang ◽  
Bo Xie

Equations of State (EOS) is crucial in simulating multiphase flows by the pseudo-potential lattice Boltzmann method (LBM). In the present study, the Peng and Robinson (P–R) and Carnahan and Starling (C–S) EOS in the pseudo-potential LBM with Exact Difference Method (EDM) scheme for two-phase flows have been compared. Both of P–R and C–S EOS have been used to study the two-phase separation, surface tension, the maximum two-phase density ratio and spurious currents. The study shows that both of P–R and C–S EOS agree with the analytical solutions although P–R EOS may perform better. The prediction of liquid phase by P–R EOS is more accurate than that of air phase and the contrary is true for C–S EOS. Predictions by both of EOS conform with the Laplace’s law. Besides, adjustment of surface tension is achieved by adjusting [Formula: see text]. The P–R EOS can achieve larger maximum density ratio than C–S EOS under the same [Formula: see text]. Besides, no matter the C–S EOS or the P–R EOS, if [Formula: see text] tends to 0.5, the computation is prone to numerical instability. The maximum spurious current for P–R is larger than that of C–S. The multiple-relaxation-time LBM still can improve obviously the numerical stability and can achieve larger maximum density ratio.

2018 ◽  
Vol 2018 ◽  
pp. 1-9 ◽  
Author(s):  
Yong Peng ◽  
Bo Wang ◽  
Yunfei Mao

Multiphase flows are very important in industrial application. In present study, the force schemes in the pseudopotential LBM for two-phase flows have been compared in detail and the force schemes include Shan-Chen, EDM, MED, and Guo’s schemes. Numerical simulations confirm that all four schemes are consistent with the Laplace law. For Shan-Chen scheme, the smaller τ is, the smaller the surface tension is. However, for other schemes, τ has no effect on surface tension. When 0.6<τ≤1, the achieved density ratio will reduce as τ reduces. During this range of τ, the maximum density ratio of EDM scheme will be greater than that of other schemes. For a constant T, the curves of the maximum spurious currents (u′) has a minimum value which is corresponding to τ′ except for EDM schemes. In the region of τ′<τ≤1, u′ will reduce as τ decreases. On the other hand, in the area of τ≤τ′, u′ will increase as τ decreases. However, for EDM scheme, u′ will increase as τ increases.


Author(s):  
Takeshi Seta ◽  
Kenichi Okui ◽  
Eisyun Takegoshi

We propose a lattice Boltzmann model capable of simulating nucleation. This LBM modifies a pseudo-potential so that it recovers a full set of hydrodynamic equations for two-phase flows based on the van der Waals-Cahn-Hilliard free energy theory through the Chapman-Enskog expansion procedure. Numerical measurements of thermal conductivity and of surface tension agree well with theoretical predictions. Simulations of phase transition, nucleation, pool boiling are carried out. They demonstrate that the model is applicable to two-phase flows with thermal effects. Using finite difference Lattice Boltzmann method ensures numerical stability of the scheme.


Fluids ◽  
2021 ◽  
Vol 6 (12) ◽  
pp. 434
Author(s):  
Assetbek Ashirbekov ◽  
Bagdagul Kabdenova ◽  
Ernesto Monaco ◽  
Luis R. Rojas-Solórzano

The original Shan-Chen’s pseudopotential Lattice Boltzmann Model (LBM) has continuously evolved during the past two decades. However, despite its capability to simulate multiphase flows, the model still faces challenges when applied to multicomponent-multiphase flows in complex geometries with a moderately high-density ratio. Furthermore, classical cubic equations of state usually incorporated into the model cannot accurately predict fluid thermodynamics in the near-critical region. This paper addresses these issues by incorporating a crossover Peng–Robinson equation of state into LBM and further improving the model to consider the density and the critical temperature differences between the CO2 and water during the injection of the CO2 in a water-saturated 2D homogeneous porous medium. The numerical model is first validated by analyzing the supercritical CO2 penetration into a single narrow channel initially filled with H2O, depicting the fundamental role of the driving pressure gradient to overcome the capillary resistance in near one and higher density ratios. Significant differences are observed by extending the model to the injection of CO2 into a 2D homogeneous porous medium when using a flat versus a curved inlet velocity profile.


Sign in / Sign up

Export Citation Format

Share Document