Galilean Invariant Preconditioned Central Moment Lattice Boltzmann Method without Cubic Velocity Errors for Efficient Steady Flow Simulations

TL;DR

This paper introduces a Galilean invariant preconditioned central moment lattice Boltzmann method that eliminates cubic velocity errors, improving accuracy and efficiency for steady flow simulations.

Contribution

It develops a GI formulation of the preconditioned cascaded central moment LB method that removes cubic velocity errors on standard lattices, with novel correction strategies.

Findings

Elimination of cubic velocity errors in the LB method.

Enhanced accuracy through velocity-dependent corrections.

Accelerated convergence in steady flow simulations.

Abstract

Lattice Boltzmann (LB) models used for the computation of fluid flows represented by the Navier-Stokes (NS) equations on standard lattices can lead to non-Galilean invariant (GI) viscous stress involving cubic velocity errors. This arises from the dependence of their third order diagonal moments on the first order moments for standard lattices, and strategies have recently been introduced to restore GI without such errors using a modified collision operator involving either corrections to the relaxation times or to the moment equilibria. Convergence acceleration in the simulation of steady flows can be achieved by solving the preconditioned NS equations, which contain a preconditioning parameter that alleviates the numerical stiffness. In the present study, we present a GI formulation of the preconditioned cascaded central moment LB method used to solve the preconditioned NS equations, which is free of cubic velocity errors on a standard lattice. A Chapman-Enskog analysis reveals the structure of the spurious non-GI defect terms and it is demonstrated that the anisotropy of the resulting viscous stress is dependent on the preconditioning parameter, in addition to the fluid velocity. It is shown that partial correction to eliminate the cubic velocity defects is achieved by scaling the cubic velocity terms in the off-diagonal third-order moment equilibria with the square of the preconditioning parameter. Furthermore, we develop additional corrections based on the extended moment equilibria involving gradient terms with coefficients dependent locally on the fluid velocity and the preconditioning parameter. Several conclusions are drawn from the analysis of the structure of the non-GI errors and the associated corrections, with particular emphasis on their dependence on the preconditioning parameter. Improvements in accuracy and convergence acceleration are demonstrated.