Two-phase flow in heterogeneous porous media: A multiscale digital model approach