A numerical code is presented for the numerical simulation of multi-fluid flows in pyrometallurgical unit operations. The code, MFVOF, is a finite difference code that allows transient solutions to immiscible multi-fluid flow problems to be generated in 2-D Cartesian and cylindrical geometries. The code is based on the use of an accurate PLIC (piecewise linear interface calculation) volume-tracking scheme to track the distortions of fluid bodies, with a redistribution procedure to ameliorate the formation of subgrid-scale fluid ligaments. Other recent enhancements to the code include swirl and expansion/compression modelling. The MFVOF code is not only able to model flows with low-curvature interface traditional applications of volume tracking - it is also able to generate robust and realistic transient representations of fragmentation and coalescence involving higher-curvature interfaces. This suggests that volume tracking can become increasingly useful for multi-phase flow modelling in chemical engineering unit operations, beyond traditional civil engineering and metallurgical applications