Numerical approach based on the combination of JFNK and SIMPLE methods are implemented for two-phase flow simulation. The equations are discretized through staggered grid by the finite volume method. Fully implicit scheme is implemented for the time discretization. Accuracy, CPU time and dynamic time step of several combination of JFNK and SIMPLE algorithm are analyzed.