Numerical optimisation of the pseudopotential-based lattice Boltzmann method

Pseudopotential-based lattice Boltzmann models provide an efficient means to simulate multiphase flows. Although thoroughly analysed and extended, they still suffer from various numerical deficiencies, such as stability issues at large density and viscosity ratios as well as spurious velocities in interfacial regions.In this work, we show how to optimise the model numerically with regard to the latter. This optimised scheme is analysed in static as well as oscillating droplet simulations and is compared to different higher-order isotropic schemes.We found that the optimised scheme substantially lowers spurious velocities and competes well with comparable schemes in both test cases.