Force approach for the pseudopotential lattice Boltzmann method
The pseudopotential method is one of the most popular extensions of the lattice Boltzmann method (LBM) for phase change and multiphase flow simulation. One attractive feature of the original proposed method consists on its simplicity of adding a force dependent on a nearest-neighbor potential function, which became known as the Shan-Chen interaction force. Some of the well known drawbacks implied by this method involves lack of thermodynamic consistency and impossibility to control the surface tension independently. In order to correct these deficiencies, different approaches were developed in the literature, such as multirange interactions potential, which involves larger stencils than nearest-neighbor approach, and modified forcing schemes. In this work, a strategy is developed to control the liquid-gas density ratio and the surface tension by means of an appropriate interaction force field using only nearest-neighbor interactions. The proposed procedure is devised starting from the desired pressure tensor, which allow for the control of the equilibrium multiphase properties such as liquid-gas coexistence curve and surface tension. Then, it is shown how to derive an external force field able to replicate the effects of this pressure tensor in the macroscopic conservation equations. The final step of our procedure is implementing this external force in the LBE by using the classical Guo forcing scheme. Numerical tests regarding static and dynamic flow conditions were performed. Results obtained from simulations showed good agreement with expected analytical values. Most divergent solution observed was the droplet oscillation period under certain flow conditions, which deviated 9% from expected analytical result. The observed results corroborate that the proposed method is able to replicate the desired macroscopic multiphase behaviour.