Consistent multiple-relaxation-time lattice Boltzmann method for the volume averaged Navier-Stokes equations

Recently, we notice that a pressure-based lattice Boltzmann (LB) method was established to recover the volume-averaged Navier-Stokes equations (VANSE), which serve as the cornerstone of various fluid-solid multiphase models. It decouples the pressure from density and exhibits excellent numerical performance, however, the widely adopted density-based LB scheme still suffers from significant spurious velocities and inconsistency with VANSE. To remedy this issue, a multiple-relaxation-time LB method is devised in this work, which incorporates a provisional equation of state in an adjusted density equilibrium distribution to decouple the void fraction from density. The Galilean invariance of the recovered VANSE is guaranteed by introducing a penalty source term in moment space, effectively eliminating unwanted numerical errors. Through the Chapman-Enskog analysis and detailed numerical validations, this novel method is proved to be capable of recovering VANSE with second-order accuracy consistently, and well-suited for handling void fraction fields with large gradients and spatiotemporal distributions.
Comment: 23 pages