Embedding sharp interfaces within the lattice Boltzmann method for fluids with arbitrary density ratios

Current lattice Boltzmann methods for simulating two fluids create a diffuse interface between the fluids. In this work, we develop a novel technique for embedding sharp interfaces between fluids with unbounded density ratios for the LB method. Distribution functions streamed across an interface are transformed so that the receiving node is passed information corresponding to its fluid phase. Two different methods are employed to determine the transformation. The first uses analytical distribution functions from steady Poiseuille flow to determine the jump in moments of the distribution functions across the interface. The second uses approximate expansions of distribution functions to determine jumps in distribution functions. The accuracy and stability of the methods are examined in simulations of Poiseuille-Couette flows with an interface parallel to the walls. Both methods show linear convergence to the analytical solution.