学术文献库

The lattice Boltzmann method (LBM) has gained increasing popularity in incompressible viscous flow simulations, but it uses many more variables than necessary. This defect was overcome by a recent approach that solves the more actual macroscopic equations obtained through Taylor series expansion analysis of the lattice Boltzmann equations [Lu et al., J. Comp. Phys., 415, 109546 (2020)]. The key is to keep some small additional terms (SATs) to stabilize the numerical solution of the weakly compressible Navier-Stokes equations. However, there are many SATs that complicate the implementation of their method. Based on some analyses and numerous tests, we ultimately pinpoint two essential ingredients for stable simulations: (1) suitable density (pressure) diffusion added to the continuity equation; (2) proper numerical dissipation related to the velocity divergence added to the momentum equations. Then, we propose a simplified method that is not only easier to implement but noticeably faster than the original method and the LBM. It contains much simpler SATs that only involve the density (pressure) derivatives and it requires no intermediate steps or variables. Besides, it is extended for two-phase flows with uniform density and viscosity. Several test cases, including some two-phase problems under two dimensional, axisymmetric and three dimensional geometries, are presented to demonstrate its capability. This work may help pave the way for the simplest simulation of incompressible viscous flows on collocated grids based on the artificial compressibility methodology.