An efficient operator-split CICSAM scheme for three-dimensional multiphase-flow problems on Cartesian grids