“…To simulate fluid flow in real pore structures, different pore‐scale simulation methods can be used, including the lattice Boltzmann method (LBM), traditional computational fluid dynamic (CFD) method based on Navier‐Stokes (N‐S) equations (Prodanović & Bryant, ; Raeini et al, ), dissipative particle dynamics (DPD) method (Groot & Warren, ), and smoothed particle hydrodynamics (SPH) method (Tartakovsky & Meakin, ). Because of its compelling advantages such as the simplicity in coding, dealing with the complex solid boundaries and parallelizing, LBM has become a very popular CFD method in many areas, such as multiphase and multicomponent flow (Chen et al, ; Liu et al, ; Martys & Chen, ; Zhang et al, ), reactive flow (Jiang & Tsuji, ; Kang et al, ), microscale and nanoscale flow (Zhao et al, ). There are different multiphase flow LB models (Huang et al, ), including the color‐gradient model (Ba et al, ; Gunstensen et al, ; Liu et al, ), pseudo‐potential model (Shan & Chen, ; Sukop, ), free‐energy model (Swift et al, ), and mean‐field model (He et al, ).…”