COUPLING MULTIPHASE PORE-SCALE MODELS TO ACCOUNT FOR BOUNDARY CONDITIONS: APPLICATION TO 2D QUASI-STATIC PORE NETWORKS