A code based on finite element method was built and applied on the variable density incompressible Navier-Stokes equations for accurately simulating immiscible two phase flows. The algorithm simulates the interface ...