In this talk an overview of the development and application of a two phase compressible Euler numerical solver will be given. Such a solver permits the study of the propagation of shock waves in both gases and liquids. The numerical solver uses finite difference methods in both the spatial and temporal dimensions whilst the interface between the two phases is captured using a level set function. The ghost fluid method is employed to impose accurate, thermodynamically correct, boundary conditions at the interface which also suppresses non-physical oscillations in the solution variables. A sample of results illustrating gas bubble/shock interactions in a number of different geometry configurations will be given and an overview of the intended future developments of the compressible solver to aid the understanding of the process of cavitation damage will be given.