Three-dimensional cascaded lattice Boltzmann method: improved implementation and consistent forcing scheme

TL;DR

This paper introduces an improved 3D cascaded lattice Boltzmann method with simplified implementation and a consistent forcing scheme, enhancing accuracy, efficiency, and applicability for force-driven flow simulations.

Contribution

It presents a simplified central moment set and a multi-relaxation-time framework for 3D CLBM, extending the forcing scheme for better force incorporation and reduced computational cost.

Findings

Reduced computational cost compared to non-orthogonal CLBM

High accuracy and convergence in force-driven flow simulations

Enhanced consistency with the no-slip boundary condition

Abstract

Cascaded or central-moment-based lattice Boltzmann method (CLBM) proposed in [Geier \textit{et al.}, Phys. Rev. E \textbf{63}, 066705 (2006)] possesses very good numerical stability. However, two constraints exist in three-dimensional (3D) CLBM simulations. Firstly, the conventional implementation for 3D CLBM involves cumbersome operations and requires much higher computational cost compared to the single-relaxation-time (SRT) LBM. Secondly, it is a challenge to accurately incorporate a general force field into the 3D CLBM. In this paper, we present an improved method to implement CLBM in 3D. The main strategy is to adopt a simplified central moment set, and carry out the central-moment-based collision operator based on a general multi-relaxation-time (GMRT) framework. Next, the recently proposed consistent forcing scheme in CLBM [L. Fei and K. H. Luo, Phys. Rev. E \textbf{96}, 053307 (2017)] is extended to incorporate a general force field into 3D CLBM. Compared with the recently developed non-orthogonal CLBM [A. D. Rosis, Phys. Rev. E \textbf{95}, 013310 (2017)], our implementation is proved to reduce the computational cost significantly. The inconsistency of adopting the discrete equilibrium distribution functions (EDFs) in the non-orthogonal CLBM is revealed and discussed. The 3D CLBM developed here in conjunction with the consistent forcing scheme is verified through numerical simulations of several canonical force-driven flows, highlighting very good properties in terms of accuracy, convergence and consistency with the nonslip rule. Finally, the techniques developed here for 3D CLBM can be applied to make the implementation and execution of 3D MRT-LBM much more efficient.