"The importance of two-phase fluid flow modelling arises from many applications. However, the non- linearity of the system makes it a complicated task for the numerical methods. While a variety of numerical techniques to solve these problems exist, these strategies can lead to spurious oscillations of the solution near the interface. In this talk a problem of the shock-induced bubble collapse near a wall computed based on the explicit finite volume solver with underlying four-equation model will be discussed. The physical dynamics involved into this problem are characterised by high speeds and very small spacial-temporal scales. A very fine grid and fast converging and compact high-order numerical schemes are, thus, required. The mesh stretching methods coupled with modified numerical schemes implemented by using Open MP and MPI paradigm are used to reduce the CPU cost. Hence, the novelty of our work is a construction of the high-order numerical tool for solving a 3D problem of two-phase shock-interface interaction on non-uniform grid."