An improved multiphase lattice Boltzmann flux solver with phase interface compression for incompressible multiphase flows

Numerical dissipation is ubiquitous in multiphase flow simulation. This paper introduces a phase interface compression term into the recently developed multiphase lattice Boltzmann flux solver and achieves an excellent interface maintenance. Here, the phase interface compression term only works in the interface region and is solved as the flux in finite volume discretization. At each cell interface, the interfacial compression velocity ur is determined by local reconstruction velocities of the multiphase lattice Boltzmann flux solver, which maintains the consistency of the flux evaluation. Meanwhile, the interfacial order parameter C in the phase interface compression term is obtained by the second order upwind scheme according to the interface normal direction. Numerical validation of the present model has been made by simulating the Zalesak problem, the single vortex problem, Rayleigh–Taylor instability, and bubble rising and coalescence. The obtained results indicate the validity and reliability of the present model.