Research Article  Open Access
Yong Peng, Bo Wang, Yunfei Mao, "Study on Force Schemes in Pseudopotential Lattice Boltzmann Model for TwoPhase Flows", Mathematical Problems in Engineering, vol. 2018, Article ID 6496379, 9 pages, 2018. https://doi.org/10.1155/2018/6496379
Study on Force Schemes in Pseudopotential Lattice Boltzmann Model for TwoPhase Flows
Abstract
Multiphase flows are very important in industrial application. In present study, the force schemes in the pseudopotential LBM for twophase flows have been compared in detail and the force schemes include ShanChen, EDM, MED, and Guo’s schemes. Numerical simulations confirm that all four schemes are consistent with the Laplace law. For ShanChen scheme, the smaller is, the smaller the surface tension is. However, for other schemes, has no effect on surface tension. When , 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 , the curves of the maximum spurious currents has a minimum value which is corresponding to except for EDM schemes. In the region of , will reduce as decreases. On the other hand, in the area of , will increase as decreases. However, for EDM scheme, will increase as increases.
1. Introduction
Multiphase flows are very important in industrial application [1]. Recently, the lattice Boltzmann method (LBM) has been applied widely for studying the twophase flows and demonstrated its advantages [2–4]. The LBM for twophase flows can be divided into the following four kinds: pseudopotential model [5, 6], free energy model [7], kinetic theory based model [8], and color model [9, 10], in which the pseudopotential LBM is very popular because of its simplicity.
In pseudopotential LBM, the interactions between fluids are simulated by an artificial interparticle potential. So, the force scheme is very important to simulate twophase flows accurately. There are four main kinds of force schemes in the pseudopotential LBM: the first is Shan and Chen’s force scheme [5, 6], the second is the Exact Difference Method (EDM) scheme [11–13], the third is Method of Explicit Derivative (MED) scheme [8], and the last is Guo’s scheme which considered the discrete lattice effects [14]. Shan [15] showed that the spurious current present in LBM is due to the insufficient isotropy of operator and proposed a method to improve it. Shan [16] proposed a general approach of calculating the pressure tensor for nonideal gas LBM. Yu and Fan [17] combined adaptive mesh refinement method and lattice Boltzmann method to improve the twophase flow simulation. It should be noted that the EDM scheme is Galilean invariant and the results obtained by EDM scheme do not depend on relaxation time [11, 18]. Kupershtokh et al. [13] compared three kinds of approximation of the gradient of special potential (“local,” “meanvalue,” and “general”) and showed that the “general” approximation was most precise and stable. Moreover, the maximal density ratio was larger than for the van der Waals Equation of State and up to for the CarnahanStarling Equation of State. Recently, Huang et al. [19] and Sun et al. [20] investigated the performance of different forcing schemes in the pseudopotential LBM. Li et al. [21] carried out a theoretical analysis of the ShanChen and EDM schemes in the pseudopotential LBM. Hu et al. [22] carried out the comparative analysis of different force schemes in pseudopotential LBM. Zheng et al. [23] studied the different force treatments in detail in ShanChen twophase LBM by theoretical analysis.
Based on above analysis, although some studies have been carried out on the force schemes in pseudopotential LBM, but the detailed comparison of different force schemes is scarce in available literature. In present study, a detailed comparison of force schemes including ShanChen, EDM, MED, and Guo’s schemes will be carried out.
2. Pseudopotential Lattice Boltzmann Model for TwoPhase Flows
2.1. ShanChen Scheme
The first pseudopotential lattice Boltzmann model for twophase flows is proposed by Shan and Chen [5, 6] and is expressed as follows:
The equilibrium distribution functions can be calculated by
The density and flow velocity can be obtained by the following:The actual fluid velocity can be defined according to Ginzburg and Adler [24] as
2.2. EDM Scheme
The pseudopotential lattice Boltzmann model with EDM [11, 12] can readwhere , and the equilibrium distribution function is shown as follows:
The real fluid velocity can be obtained from Ginzburg and Adler [24] as follows:
2.3. Guo’s Scheme
The pseudopotential lattice Boltzmann model with Guo’s force scheme [14] can be shown as follows:with .
The density and flow velocity can be obtained by the following [24]:
2.4. MED Scheme
The pseudopotential lattice Boltzmann model with MED force scheme [8] can be shown as follows:with .
The density and flow velocity can be obtained by
3. CarnahanStarling Equation of State (CS EOS)
In order to get relatively large density ratio, the CarnahanStarling Equation of State (CS EOS) is used in the present study [25]:in which , , and and can be obtained bywhere , , are used in the present study.
4. Application Cases
4.1. TwoPhase Separation
In this section, twophase separation will be used to test four schemes of ShanChen, EDM, MED, and Guo. The force term can be calculated by the following:where is weights and, for D2Q9 model, , , .
In the present work, the interaction potential is defined according to the method by Chen et al. [26] as follows:
In computation, , 101 × 101 grids have been used; the periodic boundary conditions are used for all boundaries. The initial density for whole area is critical density, and a random density fluctuation of plus or minus 1% is added at the beginning for whole area. When the computation gets to equilibrium, gas and liquid phases will separate clearly. The simulated results are shown in Figure 1.
It can be known from Figure 1 that all of four schemes perform almost the same and agree well with the analytical solution. Besides, it is clear that the present study agrees with that by Kupershtokh et al. [13].
4.2. Surface Tension
Surface tension is of great importance in the twophase flows and its relationship with bubble radius is consistent with the Laplace law. In this section, Laplace’s law will be used to verify four schemes and the effect of and on the surface tension will be discussed. The simulated results are shown in Table 1 and Figures 2–4.

(a) ShanChen scheme
(b) EDM scheme
(c) Guo’s scheme
(d) MED scheme
(a) EDM scheme
(b) Guo’s scheme
(c) MED scheme
It can be seen from Figures 2 and 3, the pressure difference between the inside and outside the bubble is linearly related to the reciprocal of the radius. This agrees well with the Laplace’s law. Figure 3 shows that has no effect on surface tension for the EDM, MED, and Guo’s schemes.
4.3. Maximum TwoPhase Density Ratio
The study shows that the simulation will become unstable when is low for four schemes, so has a certain range for stable computation. To demonstrate the numerical stability of four schemes and get the minimum with stable computation, twophase separation case is tested by adjusting τ and the results are shown in Figure 5.
It can be seen from Figure 5 that the minimum will become larger as reduces when and it means that achieved maximum ratio of twophase densities will become smaller. For the same , the maximum density ratio of EDM scheme will be larger than those of other schemes for .
4.4. Spurious Currents
The spurious current is one of the important criteria to evaluate pseudopotential models because it will lead to the computation instability. In order to study the effect of four force schemes on spurious current, a series of tests have been carried out. The maximum spurious currents for different and have been shown in Table 2 and Figure 6. It shows that the spurious current will increase as reduces with the same for all of four schemes. That is why computation is instable when is small. The maximum spurious current of MED scheme is slightly larger than that of other schemes.

(a)
(b)
(c)
Moreover, it can be seen from Figure 6 that, with a constant , for ShanChen, MED, and Guo’s schemes, the curves of have a minimum value which is corresponding to . During this region of , will reduce as τ decreases. On the other hand, in the area of , will increase as τ decreases. But for EDM scheme, will increase as increases. Moreover, present study is compared with that by Kupershtokh et al. [13] and both of results are similar as shown in Figure 7.
5. Conclusions
In present study, the force schemes in the pseudopotential LBM for twophase flows have been compared in detail and the force schemes include ShanChen, EDM, MED, and Guo’s schemes. The LBM with four schemes have been used to study the twophase separation and surface tension. Besides, the maximum twophase density ratio and spurious currents also are discussed in detail. Based on the above study, the following conclusions can be drawn:
(1) Numerical simulations confirm that all four schemes are consistent with the Laplace law. For the EDM, MED, and Guo’s schemes, has no effect on surface tension.
(2) When , the achieved density ratio will reduce as reduces. For the same , the maximum density ratio of EDM scheme will be greater than that of other schemes .
(3) Except for the EDM scheme, the curves of the maximum spurious currents have a minimum value which is corresponding to for three schemes. In the region of , will reduce as decreases. On the other hand, in the area of , will increase as decreases. But for EDM scheme, will increase as τ increases.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
The first author would like to acknowledge the financial support of the National Natural Science Foundation of China (Grants nos. 51409183, 51579166, and 51611130203). Besides, this work was also supported by the National Key Technologies R&D Program of China (no. 2015BAD24B01).
References
 M. Ishii and T. Hibiki, ThermoFluid Dynamics of TwoPhase Flow, Springer, New York, NY, USA, 2007.
 M. C. Sukop and D. T. Thorne, Lattice Boltzmann Modeling: An Introduction for Geoscientists and Engineers, Springer, Heidelberg, Germany, 2006.
 C. K. Aidun and J. R. Clausen, “Lattice Boltzmann method for complex flow,” Annual Review of Fluid Mechanics, vol. 42, no. 1, pp. 439–472, 2010. View at: Google Scholar
 H. Safari, M. H. Rahimian, and M. Krafczyk, “Consistent simulation of droplet evaporation based on the phasefield multiphase lattice Boltzmann method,” Physical Review E, vol. 90, no. 3, Article ID 033305, 2014. View at: Publisher Site  Google Scholar
 X. Shan and H. Chen, “Lattice Boltzmann model for simulating flows with multiple phases and components,” Physical Review E, vol. 47, no. 3, pp. 1815–1820, 1993. View at: Publisher Site  Google Scholar
 X. Shan and H. Chen, “Simulation of nonideal gases and liquidgas phase transitions by the lattice Boltzmann equation,” Physical Review E, vol. 49, no. 4, pp. 2941–2948, 1994. View at: Publisher Site  Google Scholar
 M. R. Swift, W. R. Osborn, and J. M. Yeomans, “Lattice Boltzmann simulation of nonideal fluids,” Physical Review Letters, vol. 75, no. 5, pp. 830–833, 1995. View at: Publisher Site  Google Scholar
 X. Y. He, X. W. Shan, and G. D. Doolen, “Discrete Boltzmann equation model for nonideal gases,” Physical Review E, vol. 57, no. 1, pp. R13–R16, 1998. View at: Google Scholar
 A. K. Gunstensen, D. H. Rothman, S. Zaleski, and G. Zanetti, “Lattice Boltzmann model of immiscible fluids,” Physical Review A: Atomic, Molecular and Optical Physics, vol. 43, no. 8, pp. 4320–4327, 1991. View at: Publisher Site  Google Scholar
 J. Tölke, S. Freudiger, and M. Krafczyk, “An adaptive scheme using hierarchical grids for lattice Boltzmann multiphase flow simulations,” Computers & Fluids, vol. 35, no. 89, pp. 820–830, 2006. View at: Publisher Site  Google Scholar
 A. L. Kupershtokh, “New method of incorporating a body force term into the lattice Boltzmann equation,” in Proceedings of the 5th International EHD Workshop, pp. 241–246, 2004. View at: Google Scholar
 A. L. Kupershtokh, “Criterion of numerical instability of liquid state in LBE simulations,” Computers and Mathematics with Applications, vol. 59, no. 7, pp. 2236–2245, 2010. View at: Publisher Site  Google Scholar  MathSciNet
 A. L. Kupershtokh, D. A. Medvedev, and D. I. Karpov, “On equations of state in a lattice Boltzmann method,” Computers and Mathematics with Applications, vol. 58, pp. 965–974, 2009. View at: Google Scholar
 Z. Guo, C. Zheng, and B. Shi, “Discrete lattice effects on the forcing term in the lattice Boltzmann method,” Physical Review E, vol. 65, Article ID 046308, 2002. View at: Publisher Site  Google Scholar
 X. Shan, “Analysis and reduction of the spurious current in a class of multiphase lattice Boltzmann models,” Physical Review E, vol. 73, no. 2, Article ID 047701, 2006. View at: Publisher Site  Google Scholar
 X. Shan, “Pressure tensor calculation in a class of nonideal gas lattice Boltzmann models,” Physical Review E, vol. 77, no. 2, Article ID 066702, 2008. View at: Publisher Site  Google Scholar
 Z. Yu and L.S. Fan, “An interaction potential based lattice Boltzmann method with adaptive mesh refinement AMR for twophase flow simulation,” Journal of Computational Physics, vol. 228, no. 17, pp. 6456–6478, 2009. View at: Publisher Site  Google Scholar  MathSciNet
 P. J. Dellar, “Lattice Boltzmann formulation for linear viscoelastic fluids using an abstract second stress,” SIAM Journal on Scientific Computing, vol. 36, no. 6, pp. A2507–A2532, 2014. View at: Publisher Site  Google Scholar  MathSciNet
 H. Huang, M. Krafczyk, and X. Lu, “Forcing term in singlephase and ShanChentype multiphase lattice Boltzmann models,” Physical Review E: Statistical, Nonlinear, and Soft Matter Physics, vol. 84, no. 4, Article ID 046710, 2011. View at: Publisher Site  Google Scholar
 K. Sun, T. Wang, M. Jia, and G. Xiao, “Evaluation of force implementation in pseudopotentialbased multiphase lattice Boltzmann models,” Physica A, vol. 391, no. 15, pp. 3895–3907, 2012. View at: Publisher Site  Google Scholar
 Q. Li, K. H. Luo, and X. J. Li, “Forcing scheme in pseudopotential lattice Boltzmann model for multiphase flows,” Physical Review E, vol. 86, no. 1, Article ID 016709, 2012. View at: Publisher Site  Google Scholar
 A. Hu, L. Li, and R. Uddin, “Force method in a pseudopotential lattice Boltzmann model,” Journal of Computational Physics, vol. 294, pp. 78–89, 2015. View at: Publisher Site  Google Scholar  MathSciNet
 L. Zheng, Q. Zhai, and S. Zheng, “Analysis of force treatment in the pseudopotential lattice Boltzmann equation method,” Physical Review E, vol. 95, no. 4, Article ID 043301, 2017. View at: Publisher Site  Google Scholar
 I. Ginzburg and P. M. Adler, “Boundary flow condition analysis for the threedimensional lattice Boltzmann model,” Journal de Physique II, vol. 4, no. 2, pp. 191–214, 1994. View at: Publisher Site  Google Scholar
 N. F. Carnahan and K. E. Starling, “Intermolecular repulsions and the equation of state for fluids,” AIChE Journal, vol. 18, no. 6, pp. 1184–1189, 1972. View at: Publisher Site  Google Scholar
 S. Chen, X. He, and L.S. Luo, “Computational methods for multiphase flow,” in Lattice Boltzmann Models for Multiphase Flow, A. Prosperetti and G. Tryggvason, Eds., chapter 6, pp. 157–192, Cambridge University Press, Cambridge, UK, 2007. View at: Google Scholar
Copyright
Copyright © 2018 Yong Peng et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.