Multiphase flow in porous media is related to various natural and industrial processes, such as geological CO 2 sequestration (Middleton et al., 2012; Pentland et al., 2011), enhanced oil/gas recovery (Lake, 1989; Orr & Taber, 1984), and groundwater contamination by D(L)NAPL (Dawson & Roberts, 1997; Molnar et al., 2020). When a less viscous fluid invades into the porous media to displace another more viscous one, it always causes interfacial instability and generates different displacement patterns. In geological CO 2 sequestration, the stages of CO 2 injection and postinjection give rise to complex flow behaviors in geological media and result in fingering flow that significantly increases the efficiency of storage by improving capillary/residual trapping (Bachu, 2015; Y. Wang et al., 2012; Yamabe et al., 2014). In enhanced oil recovery, due to the unstable water-oil displacement interface, the sweep efficiency of water-flooding is greatly reduced and more than 80% native oil is left in reservoirs (Orr & Taber, 1984; Veld & Phillips, 2010). Understanding the fundamental mechanism and controlling displacement patterns of multiphase flow in porous media is therefore, critical for optimizing subsurface resources management. Fluid-fluid interfacial instability is a central issue for immiscible multiphase flow in permeable media. Numerous experimental, numerical and theoretical efforts have been devoted in the past semi-century (Arm