I would say take a look at the TwoPhaseVOF library by DLR-RY, specifically the cutCell, cutFace and surfaceIterator classes Edit: link https://github.com/DLR-RY/VoFLibrary