Nice simulation. Thanks for sharing! I've been working on similar CFD but I gave up on plain Navier-Stokes because it is hard to implement in the GPU shaders, and then it becomes too slow to compute. So I've been working with LBM (Lattice Boltzmann Method). Have you tried that approach?