Abstract- The aim of this paper is to develop a numerical procedure for simulating underwater explosion phenomena with a simplified mathematical and two fluid model. The two fluid Kapila five equation model is selected as the governing equations and the ideal gas and Stiffened gas equations of state (SG-EOS) are used to obtain pressure in the gas bubble and the surrounding water zone, respectively. The modified Schmidt EOS is used to simulate the cavitation regions with low pressure. A Godunov numerical method and HLLC Reiman solver is extended for Kapila two fluid model. The numerical results of the present method and comparing them with available experimental results, verify that the proposed method has good capablity of predicting complex physics involved in a spherical underwater explosion and its interaction with free surface. The method also shows a very good performance with no spurious oscillation in cavitation zone simulation in two-dimensional problems