Abstract: A method for tracking N fluid materials and their associated interfaces during simulated fluid flow is disclosed. A microgrid cell methodology is embedded on a regular macrogrid to subdivide and then tag fluid materials in a computational system preferably using a prime numbering algorithm. The motion of microgrid cells is tracked based on local velocity conditions, rectifying small anomalies by a coupled evaluation of local volume fraction fields and global mass conservation. Volume fractions can be calculated at any time step via an evaluation of the prime locations so that average cellular density and viscosity values can be regularly updated.