Equations of state in multiphase lattice Boltzmann method revisited

The single-component multiphase fluids can be described by a single equation of state (EOS), and various EOSs have been employed in the multiphase lattice Boltzmann (LB) method. In this work, we revisit five commonly used EOSs, including the van der Waals EOS, the Redlich-Kwong EOS, the Redlich-Kwong-Soave EOS, the Peng-Robinson EOS, and the Carnahan-Starling EOS. The recent multiphase LB model with self-tuning EOS is employed because of its thermodynamic consistency in a strict sense and clear physical picture at the microscopic level. First, the way to incorporate these multiphase EOSs is proposed. Two scaling factors are introduced to independently adjust the surface tension and interface thickness, and the lattice sound speed is EOS-dependent to ensure the numerical stability. Then, numerical tests are conducted to validate the incorporations of these EOSs and compare their numerical performances. The surface tension and interface thickness are set to the same values for different EOSs in the comparisons. The liquid and gas densities, surface tension, and interface thickness by the LB simulation agree well with the thermodynamic results. The maximum density ratios achieved with different EOSs are at the same level and could be very close to each other when the interface thickness is relatively small. The effects of multiphase EOS, density ratio, and dimensionless relaxation time on the spurious current are discussed in detail. It is interesting to find the van der Waals EOS shows the best numerical performance in reducing the spurious current.