Numerics of the Lattice Boltzmann Method: Effects of Collision Models on the Lattice Boltzmann Simulations

We conduct a comparative study to evaluate several lattice Boltzmann (LB) models for solving the near incompressible Navier-Stokes equations, including the lattice Boltzmann equation with the multiple-relaxation-time (MRT), the two-relaxation-time (TRT), the single-relaxation-time (SRT) collision models, and the entropic lattice Boltzmann equation (ELBE). The lid-driven square cavity flow in two dimensions is used as a benchmark test. Our results demonstrate that the ELBE does not improve the numerical stability of the SRT or the lattice Bhatnagar-Gross-Krook (LBGK) model. Our results also show that the MRT and TRT LB models are superior to the ELBE and LBGK models in terms of accuracy, stability, and computational efficiency and that the ELBE scheme is the most inferior among the LB models tested in this study, thus is unfit for carrying out numerical simulations in practice. Our study suggests that, to optimize the accuracy, stability, and efficiency in the MRT model, it requires at least three independently adjustable relaxation rates: one for the shear viscosity ν (or the Reynolds number Re), one for the bulk viscosity ζ, and one to satisfy the criterion imposed by the Dirichlet boundary conditions which are realized by the bounce-back-type boundary conditions. © 2011 American Physical Society

Fluid-Structure Interaction Simulation of a Coriolis Mass Flowmeter using a Lattice Boltzmann Method

In this paper we use a fluid-structure interaction (FSI) approach to simulate a Coriolis mass flowmeter (CMF). The fluid dynamics are calculated by the open source framework OpenLB, based on the lattice Boltzmann method (LBM). For the structural dynamics we employ the open source software Elmer, an implementation of the finite element method (FEM). A staggered coupling approach between the two software packages is presented. The finite element mesh is created by the mesh generator Gmsh to ensure a complete open source workflow. The Eigenmodes of the CMF, which are calculated by modal analysis are compared with measurement data. Using the estimated excitation frequency, a fully coupled, partitioned, FSI simulation is applied to simulate the phase shift of the investigated CMF design. The calculated phaseshift values are in good agreement to the measurement data and verify the suitability of the model to numerically describe the working principle of a CMF