Czelusniak LE, Mapelli VP, Wagner AJ, Cabezas-Gómez L. Shaping the equation of state to improve numerical accuracy and stability of the pseudopotential lattice Boltzmann method.
Phys Rev E 2022;
105:015303. [PMID:
35193246 DOI:
10.1103/physreve.105.015303]
[Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/15/2021] [Accepted: 12/09/2021] [Indexed: 11/07/2022]
Abstract
It has recently been shown that altering the shape of the metastable and unstable branches of an equation of state (EOS) can substantially improve the numerical accuracy of liquid and vapor densities in the pseudopotential lattice Boltzmann method [Peng et al., Phys. Rev. E 101, 063309 (2020)2470-004510.1103/PhysRevE.101.063309]. We found that this approach reduces stability of the method in nonequilibrium conditions and is unstable for bubbles at low reduced temperatures. Here we present an improved method for altering the shape of the metastable and unstable branches of the EOS which remains stable for both equilibrium and nonequilibrium situations and has no issues with bubbles. We also performed a detailed study of the stability of the methods for a droplet impact on a liquid film for reduced temperatures down to 0.35 with Reynolds number of 300. Our approach remained stable for a density ratio of up to 3.38×10^{4}.
Collapse