Buoyancy stabilised Kelvin-Helmholtz instability (lattice Boltzmann and finite volume)

Re = 1200 Pr = 9 Ri = 16 million nodes for the flow field simulated using the lattice Boltzmann method and 128 million cells for buoyancy simulated using the finite volume method. All calculations performed on GPU.