Simulation of laminar flow in a three‐dimensional lid‐driven cavity by lattice Boltzmann method

The purpose of this paper is to apply lattice Boltzmann equation method (LBM) with multiple relaxation time (MRT) model, to investigate lid‐driven flow in a three‐dimensional (3D), rectangular cavity, and compare the results with flow in an equivalent two‐dimensional (2D) cavity.

The second‐order MRT model is implemented in a 3D LBM code. The flow structure in cavities of different aspect ratios (0.25‐4) and Reynolds numbers (0.01‐1000) is investigated. The LBM simulation results are compared with those from numerical solution of Navier‐Stokes (NS) equations and with available experimental data.

The 3D simulations demonstrate that 2D models may predict the flow structure reasonably well at low Reynolds numbers, but significant differences with experimental data appear at high Reynolds numbers. Such discrepancy between 2D and 3D results are attributed to the effect of boundary layers near the side‐walls in transverse direction (in 3D), due to which the vorticity in the core‐region is weakened in general. Secondly, owing to the vortex stretching effect present in 3D flow, the vorticity in the transverse plane intensifies whereas that in the lateral plane decays, with increase in Reynolds number. However, on the symmetry‐plane, the flow structure variation with respect to cavity aspect ratio is found to be qualitatively consistent with results of 2D simulations. Secondary flow vortices whose axis is in the direction of the lid‐motion are observed; these are weak at low Reynolds numbers, but become quite strong at high Reynolds numbers.

The findings will be useful in the study of variety of enclosed fluid flows.