Incompressible flow in 2D of two fluids with different densities and viscosities. Buoyancy is included, but not surface tension (yet).

Source code and inputs here