An Accurate Curved Boundary Treatment in the Lattice Boltzmann Method

Abstract

The lattice Boltzmann equation (LBE) is an alternative kinetic method capable of solving hydrodynamics for various systems. Major advantages of the method are due to the fact that the solution for the particle distribution functions is explicit, easy to implement, and natural to parallelize. Because the method often uses uniform regular Cartesian lattices in space, curved boundaries are often approximated by a series of stairs that leads to reduction in computational accuracy. In this work, a second-order accurate treatment of the boundary condition in the LBE method is developed for a curved boundary. The proposed treatment of the curved boundaries is an improvement of a scheme due to O. Filippova and D. Hänel (1998, J. Comput. Phys.147, 219). The proposed treatment for curved boundaries is tested against several flow problems: 2-D channel flows with constant and oscillating pressure gradients for which analytic solutions are known, flow due to an impulsively started wall, lid-driven square cavity flow, and uniform flow over a column of circular cylinders. The second-order accuracy is observed with a solid boundary arbitrarily placed between lattice nodes. The proposed boundary condition has well-behaved stability characteristics when the relaxation time is close to 1/2, the zero limit of viscosity. The improvement can make a substantial contribution toward simulating practical fluid flow problems using the lattice Boltzmann method.