Straight velocity boundaries in the lattice Boltzmann method

Abstract

Various ways of implementing boundary conditions for the numerical solution of the Navier-Stokes equations by a lattice Boltzmann method are discussed. Five commonly adopted approaches are reviewed, analyzed, and compared, including local and nonlocal methods. The discussion is restricted to velocity Dirichlet boundary conditions, and to straight on-lattice boundaries which are aligned with the horizontal and vertical lattice directions. The boundary conditions are first inspected analytically by applying systematically the results of a multiscale analysis to boundary nodes. This procedure makes it possible to compare boundary conditions on an equal footing, although they were originally derived from very different principles. It is concluded that all five boundary conditions exhibit second-order accuracy, consistent with the accuracy of the lattice Boltzmann method. The five methods are then compared numerically for accuracy and stability through benchmarks of two-dimensional and three-dimensional flows. None of the methods is found to be throughout superior to the others. Instead, the choice of a best boundary condition depends on the flow geometry, and on the desired trade-off between accuracy and stability. From the findings of the benchmarks, the boundary conditions can be classified into two major groups. The first group comprehends boundary conditions that preserve the information streaming from the bulk into boundary nodes and complete the missing information through closure relations. Boundary conditions in this group are found to be exceptionally accurate at low Reynolds number. Boundary conditions of the second group replace all variables on boundary nodes by new values. They exhibit generally much better numerical stability and are therefore dedicated for use in high Reynolds number flows.

Feedback control of inertial microfluidics using axial control forces

“…The internal pressure follows an ideal gas law with p = c 2 s ρ, where c s = 1/ √ 3 is the speed of sound, and the kinematic viscosity is given by ν = c 2 s (τ − 1/2). To implement the no-slip boundary condition on the channel walls, we employ the regularized boundary condition introduced by Latt and Chopard 29 . It treats boundary nodes just like fluid nodes but modifies the distribution function before the collision such that the correct velocity is imposed.…”

Three-dimensional lattice Boltzmann method benchmarks between color-gradient and pseudo-potential immiscible multi-component models

“…Therefore, although several calculations reported in this work could have been performed with an SRT operator, we decided to use MRT for all of them. On the moving boundaries, regularized velocity boundary conditions 61 are applied to the blue°uid, while bounce-back boundary conditions are used for the red°u id. Here, we also adopt a fourth-order isotropic scheme for the calculation of the density gradients.…”

Progress in particle-based multiscale and hybrid methods for flow applications

“…The velocity component of u ⊥ is either prescribed via Dirichlet-or Neumann-type boundary conditions and u can be set for moving boundary problems. A comprehensive overview of various boundary conditions can be found in Latt et al (2008) and Chen et al (1996). Classical and recent applications of the LBM can be found in Chen and Doolen (1998) and Aidun and Clausen (2010).…”