Numerical simulation of incompressible flow driven by density variations during phase change