Numerical solution for hyperbolic conservative two-phase flow equations