Enhanced, mass-conserving closure scheme for lattice boltzmann equation hydrodynamics

We derive and demonstrate a practical and accurate algorithm designed to impose Dirichlet boundary conditions (specified boundary velocity) on the edge nodes of a lattice Boltzmann fluid simulation space. The current algorithm models the lattice fluid on boundary and bulk nodes to identical accuracy and in a demonstrably equivalent manner. Whilst this new method applies to rectangular geometries, it adds to our previous method [2] (i) the condition of mass conservation and (ii) the ability to treat with unimpaired accuracy both internal and external corners.