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

TL;DR

This paper introduces a new multiple-relaxation-time lattice Boltzmann method that accurately recovers volume-averaged Navier-Stokes equations, reducing spurious velocities and improving consistency for multiphase fluid simulations.

Contribution

It develops a consistent LB scheme incorporating a provisional equation of state and penalty source term, enhancing accuracy and stability over traditional density-based methods.

Findings

Achieves second-order accuracy in recovering VANSE.

Effectively handles large gradients in void fraction fields.

Reduces spurious velocities and improves numerical stability.

Abstract

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.