A Semi-Implicit Scheme of Lattice Boltzmann Method for Two Dimensional Cavity Flow Simulation

The calculation sequence of collision, propagation and macroscopic variables is not very clear in lattice Boltzmann method (LBM) code implementation. According to the definition, three steps should be computed on all nodes respectively, which mean three loops are needed. While the “pull” scheme makes the only one loop possible for coding, this is called semi-implicit scheme in this study. The accuracy and efficiency of semi-implicit scheme are discussed in detail through the simulation of lid-driven cavity flow. Non-equilibrium extrapolation scheme is adopted on the boundary of simulation area. The results are compared with two classic articles, which show that semi-implicit scheme has good agreement with the classic scheme. When Re is less than 3000, the iterations steps of semi-scheme can be decreased by about 30% though comparing the semi-implicit scheme with standard scheme containing three loops. As the Re increases into more than 3400, the standard scheme is not converged. On the contrary, the iterations of semi-implicit scheme are approximately linear to Re.

[1]  Kejian Wu,et al.  SHIFT: an implementation for lattice Boltzmann simulation in low-porosity porous media. , 2010, Physical review. E, Statistical, nonlinear, and soft matter physics.

[2]  John G. Hagedorn,et al.  Large Scale Simulations of Single and Multi-Component Flow in Porous Media | NIST , 1999 .

[3]  P. Bhatnagar,et al.  A Model for Collision Processes in Gases. I. Small Amplitude Processes in Charged and Neutral One-Component Systems , 1954 .

[4]  B. Shi,et al.  Non-equilibrium extrapolation method for velocity and pressure boundary conditions in the lattice Boltzmann method , 2002 .

[5]  U. Ghia,et al.  High-Re solutions for incompressible flow using the Navier-Stokes equations and a multigrid method , 1982 .

[6]  Y. Qian,et al.  Lattice BGK Models for Navier-Stokes Equation , 1992 .

[7]  Gerhard Wellein,et al.  On the single processor performance of simple lattice Boltzmann kernels , 2006 .

[8]  John G. Hagedorn,et al.  Large-scale simulations of single- and multicomponent flow in porous media , 1999, Optics & Photonics.

[9]  Jonas Tölke,et al.  Implementation of a Lattice Boltzmann kernel using the Compute Unified Device Architecture developed by nVIDIA , 2009, Comput. Vis. Sci..

[10]  Sauro Succi,et al.  Massively Parallel Lattice-Boltzmann Simulation of Turbulent Channel Flow , 1997 .

[11]  Christian Obrecht,et al.  LBM based flow simulation using GPU computing processor , 2010, Comput. Math. Appl..

[12]  Panayotis Skordos,et al.  Parallel simulation of subsonic fluid dynamics on a cluster of workstations , 1994, Proceedings of the Fourth IEEE International Symposium on High Performance Distributed Computing.

[13]  Shiyi Chen,et al.  Simulation of Cavity Flow by the Lattice Boltzmann Method , 1994, comp-gas/9401003.

[14]  郑楚光,et al.  Non-equilibrium extrapolation method for velocity and pressure boundary conditions in the lattice Boltzmann method , 2005 .