Trapping and mobilization of multiphase flow in porous media are important subsurface processes with broad applications such as subsurface CO 2 storage, oil and gas production, remediation of nonaqueous phase liquid contaminants in aquifers, and vadose zone hydrogeology. Water, assumed as wetting fluid in this paper for the convenience of analysis, mainly flows in small pores or the corners of large pores at the multiphase interface regions; oil, assumed as nonwetting fluid in this analysis, mainly occupies the central part of pores. When meeting the snap-off criteria (Deng et al., 2014(Deng et al., , 2015Roman et al., 2017;Roof, 1970), water film would swell at narrow pore throats, and oil would be broken up into smaller clusters and/or droplets due to snap-off events. These processes contribute to one of the main mechanisms of oil trapping. Another main mechanism for oil trapping is called by-pass cutoff which can be described by the pore doublet model (Chatzis & Dullien, 1983;Lenormand et al., 1983); as percolating oil becomes disconnected, the small oil clusters would further be fragmented into smaller droplets and trapped in the individual pores (Berg et al., 2015;Pak et al., 2015;Rücker et al., 2015). When the background pressure drop cannot overcome the capillary force exerted on these droplets, the fragmented droplets keep trapped in individual pores. When entrapped, the oil droplets stop flowing and start to contact with pore surface (Pak et al., 2015) by squeezing out water film and get pinned by their contact with mineral surface. These pinned droplets could be subjected to contact angle hysteresis (CAH), which can be attributed to geometric or chemical heterogeneity of solid surface (Rücker et al., 2020;Schmatz et al., 2015). In the ganglion dynamics, pinning and capillary trapping due to pore radius variation at pore constriction work together as the barriers for the Newtonian droplet mobilization.The equilibrium-based concept of multiphase flow by using Darcy's law becomes more and more insufficient to address the problems of trapping and mobilization processes in porous media. The conventional percolation model assumes that disconnected oil droplets and ganglia become immobile and remain stationary (Wilkinson, 1986;Wilkinson & Willemsen, 1983). However, the effect of ganglion dynamics on fluid topology and residual saturation was analyzed (Berg et al., 2015;Rücker et al., 2015;Youssef et al., 2014) and was found to be able to reduce another 10% of oil saturation after the oil get disconnected. In both two-dimensional (2D) micromodels (Avraam & Payatakes, 1995) and three-dimensional (3D) porous medium experiments (Rücker et al., 2015), it was observed that rich pore-scale events could cause local nonequilibrium