Multi-fluid flows involving shock-accelerated inhomogeneities and shock-induced instability play essential roles in a wide variety of problems including, but not limited to, supersonic combustion [1], inertial confinement fusion [2], and supernova explosion [3]. Numerical simulations of these complex flows prove to be challenging in the presence of moving and deformable material interfaces, especially for fluids with large differences in their densities or thermodynamic properties. Therefore, a discontinuity-capturing, mass-conserving, and positivity-preserving scheme is desirable for compressible multi-fluid simulations.