The capillary bundle model, wherein the flow dynamics of a porous medium is predicted from that independent cylindrical tubes/capillaries whose radii are distributed according to medium's pore size distribution, has been used extensively. However, as it lacks an interaction between channels, this model fails at predicting complex configuration, including those involving two-phase flow. We propo...