Towards the simplest simulation of incompressible viscous flows inspired by the lattice Boltzmann method

TL;DR

This paper introduces a simplified, faster, and easier-to-implement simulation method for incompressible viscous flows inspired by the lattice Boltzmann method, reducing complexity while maintaining stability and extending to two-phase flows.

Contribution

The authors propose a minimalistic approach with simpler stabilization terms, eliminating intermediate steps, and extend it to two-phase flow simulations, improving efficiency and ease of implementation.

Findings

The new method is faster than existing approaches.

It effectively simulates two-phase flows with uniform density and viscosity.

Numerical tests confirm stability and accuracy across various geometries.

Abstract

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.