Accurately capturing sharp interfaces on a fixed Eulerian grid is a major challenge for computational fluid dynamics. Attempts to minimise phase boundaries in multiphase LBMs cause them to become fixed to the underlying computational grid instead of advecting with the fluid velocity. Paul Dellar and I proposed lattice Boltzmann implementations of conservation equations with stiff source terms to better understand lattice pinning. By incorporating stochastic sharping we were able to eliminating pinning a facetting and correctly compute very narrow interfaces.

An overview of this work is given in the following slides, which were intended for a recent UKCOMES meeting.

The work above offers a theoretical explanation of lattice pinning, and some ways to eliminate (or at least reduce ) it, but the model does not readily lend itself to flows with significant local changes in interface topology. This is because of the variable, and re-calculated, sharpening threshold (but on the plus side, no numerical differentiation in needed so the sharpening is a purely local and algebraic source term!). I proposed a alternative and perhaps more usable (in terms of engineering applications) Conservative Interface Sharpening Lattice Boltzmann Model. The price to pay is some (nonlocal) numerical differentiation on the lattice, but this is required in all multiphase LBMs.

The resulting numerical scheme is nevertheless simple, accurate, does not leak mass, and can capture significant interface deformations, as can be seen in these two plots of the stringent Vortex Reversal benchmark test. This work was inspired by the Artificial Compression Model of Olsson, Kreiss,and Zahedi.