Numerical simulation of two phase flows