Implicit penalty methods on eulerian grids for the simulation of incompressible multiphase flows