In this article, for the first time the numerical solution of Godunov method with HLLC Riemann solver is extended for a hyperbolic five equations two-fluid model. The flow field is considered for two-space dimensional case. So far, two main difficulties include non-monotonic behavior of mixture sound relation and inability of shock transition from interface was mentioned during working with this model. In this research these difficulties were overcome by selecting an appropriate mixture sound relation and appropriate discretization of non-conservative term. The mutual effect of shock wave impact with a droplet and two droplets with different diameters were simulated and studied. During the shock wave impact with 1.47 and 6 Mach with the droplet, a complicated shape of interface was formed with high pressure zone and low pressure zone of cavitations. The results obtained from the present attempt were compared with the experimental and related similar results of that obtained by the other numerical methods and models. The comparison of the results was good. It was also concluded that the numerical method used in the present work has enough accuracy with high capability in capturing two-phase flow interfacial instability and shock wave impact transmission from the droplet.