A lattice kinetic scheme for bubble flows

A lattice kinetic scheme for two–phase immiscible fluids with large density ratios is presented. The difficulty in the treatment of large density ratio is resolved by using the projection method. The method can simulate two–phase fluid flows with a density ratio up to 1000. The method is applied to the simulation of many bubbles rising in a long square duct. The complicated unsteady structures of the interface and the velocity field can be simulated in a stable manner.