The simulation of multicomponent fluids at low Reynolds number and low capillary number is of interest in a variety of applications such as the modeling of venule scale blood flow and microfluidics; however, such simulations are computationally demanding. An improved multicomponent lattice Boltzmann scheme, designed to represent interfaces in the continuum approximation, is presented and shown (i) significantly to reduce common algorithmic artifacts and (ii) to recover full Galilean invariance. The method is used to model drop dynamics in shear flow in two dimensions where it recovers correct results over a range of Reynolds and capillary number greater than that which may be addressed with previous methods.

• Received 26 May 2005

DOI:https://doi.org/10.1103/PhysRevE.73.056708