A three-dimensional high-order flux reconstruction lattice boltzmann flux solver for incompressible laminar and turbulent flows

Abstract

A three-dimensional (3D) high-order flux reconstruction lattice Boltzmann flux solver (FR-LBFS) is developed in this paper for accurate and efficient simulations of incompressible laminar and turbulent flows. The original lattice Boltzmann flux solver (LBFS) is a second-order finite volume method (FVM) which combines the advantages of conventional Navier-Stokes (N-S) solvers and lattice Boltzmann equation (LBE) solvers. But it is not convenient to construct high-order FVM, and the compactness is always a bottleneck. The present work separates LBFS from FVM and combines LBFS with high-order FR scheme. Thus, the FR-LBFS inherits the superiority of LBFS and FR and shares the following attractive features: (i) a fully-explicit arbitrary order compact method, (ii) weakly compressible (WC) model for unsteady incompressible flows and (iii) evaluating inviscid and viscous fluxes simultaneously and no particular technique requirement, such as local discontinuous Galerkin method (LDG), for viscous term discretization. To validate the current method, various 3D incompressible laminar isothermal and thermal numerical experiments are presented first. Good agreements are achieved between the current results and the previous experimental and numerical data whilst the present high-order solver adopt coarse meshes. Furthermore, the ability of the proposed method to perform accurate and stable computations of incompressible decaying homogeneous isotropic turbulent flows and turbulent channel flows with different mesh resolutions and accuracy orders are investigated. The results show that the present 3D FR-LBFS is a promising tool for under-resolved or implicit large eddy simulation (ILES) of incompressible turbulent flows.