-
Compared with the lattice Boltzmann equation (LBE) model based on incompressible phase field theory, the LBE based on quasi-incompressible phase field theory has the advantage of local mass conservation. However, previous quasi-incompressible phase-field-based LBE model does not satisfy the well-balance property, resulting in spurious velocities in the vicinity of interface and density profiles inconsistent with thermodynamics.
To address this difficulty, a novel LBE model is developed based on the quasi-incompressible phase-field theory. First, numerical artifacts in the original LBE for the Cahn-Hilliard are analyzed. Based on this analysis, the equilibrium distribution function and source term are reformulated to remove the numerical artifacts, which makes the new LBE realize the well-balance property at discrete level.
The performance of the proposed LBE model is tested by simulating a number of typical two-phase systems. The numerical results of the planar interface and static droplet problems demonstrate that the present model can remove spurious velocities and achieve well-balanced state. Numerical results of the layered Poiseuille flow demonstrate the accuracy of the present model in simulating dynamic two-phase flow problems. The well-balance properties of the LBE model with two different formulations of surface tension ($\boldsymbol{F}_s=-\phi \nabla \mu$ and $\boldsymbol{F}_s=\mu \boldsymbol{\nabla} \phi$) are also investigated. It is found that the formulation of $\boldsymbol{F}_s=\mu \boldsymbol{\nabla} \phi$ cannot eliminate the spurious velocities, while the formulation of $\boldsymbol{F}_s=\mu \boldsymbol{\nabla} \phi$ can achieve the well-balance state. The influences of viscosity formulations of the fluid mixture are also compared. Particularly, four mixing rules are considered. It is shown that the use of step mixing rule gives more accurate results for the layered Poiseuille flow. Finally, we compare the performance of the present quasi-incompressible LBE model and the original fully incompressible LBE model by simulating the phase separation problem, and the results show that the present model can ensure the local mass conservation, while the fully incompressible LBE can yield quite different predictions. -
[1] Guo Z L, Shu C 2013 Lattice Boltzmann Method and Its Application in Engineering (Vol. 3) (Singapore: World Scientific Publishing)
[2] Huang H B, Sukop M, Lu X Y 2015 Multiphase lattice Boltzmann methods: Theory and application (Hoboken, NJ: Wiley-Blackwell)
[3] Yang K, Guo Z L 2016 Phys. Rev. E 93 043303
[4] Zhang C H, Guo Z L, Liang H 2021 Int. J. Numer. Methods Fluids 93 2225
[5] Connington K, Lee T 2012 J. Mech. Sci. Technol. 26 3857
[6] Gong J M, Oshima N, Tabe Y 2019 Comput. Math. Appl. 78 1166
[7] Huang H B, Krafczyk M, Lu X Y 2011 Phys. Rev. E 84 046710
[8] Zhai Q L, Zheng L, Zheng S 2017 Phys. Rev. E 95 023313
[9] Siebert D, Philippi P, Mattila K 2014 Phys. Rev. E 90 053310
[10] Cristea A, Sofonea V 2003 Int. J. Mod. Phys. C 14 1251
[11] Wagner A J 2003 Int. J. Mod. Phys. B 17 193
[12] Shan X 2006 Phys. Rev. E 73 047701
[13] Sbragaglia M, Benzi R, Biferale L, Succi S, Sugiyama K, Toschi F 2007 Phys. Rev. E 75 026702
[14] Guo Z L, Zheng C G, Shi B C 2011 Phys. Rev. E 83 036707
[15] Guo Z L 2021 Phys. Fluids 33 031709
[16] Zhang C H, Guo Z L, Wang L P 2022 Phys. Fluids 34 012110
[17] Zheng L, Zheng S, Zhai Q L 2021 Phys. Fluids 33 092102
[18] Ju L, Liu P Y, Yan B C, Bao J, Sun S Y, Guo Z L 2023 arXiv:2311.10827v1 [physics.flu-dyn]
[19] Lowengrub J, Truskinovsky L 1998 Proc. R. Soc. Lond. A. 454 2617
[20] Shen J, Yang X F, Wang Q 2013 Commun. Comput. Phys. 13 1045
[21] Jacqmin D 1999 Commun. Comput. Phys. 155 96
[22] Deng B, Shi B C, Wang G C 2005 Chin. Phys. Lett. 22 267
Metrics
- Abstract views: 307
- PDF Downloads: 10
- Cited By: 0
下载:
