Solving incompressible two-phase flows on multi-GPU clusters