<p>Interface tracking simulation of multifluid flow is challenging due to the unstable multiphase interface involves complex topological changes. In this paper, an ISPH-FVM coupling method combining the incompressible smooth particle hydrodynamics (ISPH) method and the finite volume method (FVM) is developed for interface tracking simulation of multifluid flow. In the present ISPH-FVM coupling method, the flow filed is calculated by FVM, while the tracking of phase interfaces is characterized by ISPH particles. The bidirectional information transfer between ISPH particles and FVM grids is implemented by the particle approximation interpolation technique. The continuous surface stress (CSS) model is introduced into the present ISPH-FVM coupling method to evaluate the effect of the interphase surface tension. Moreover, the surface tension between phase interfaces is identified and calculated by the volume fraction on the grid. In addition, the volume fraction on the grid can also be easily obtained through the particle approximation interpolation technique. To verify and test the ability of the present ISPH-FVM coupling method in interface tracking simulation of multifluid flow, several numerical examples are implemented. The results show that the present ISPH-FVM coupling method has good performance in interface tracking simulation of multifluid flow with complex interface topology changes.</p>