“…During displacement, the governing equations that describe the mass balance and macro DCP have been given as (Dahle et al, ; Tang et al, ; Yu et al, ) where is the Darcy velocity vector of phase α (the arrow above a variable denotes a vector hereinafter, and since we focus on two‐phase flow in this study, α is either “ n ” for nonwetting phase or “ w ” for wetting phase.) η α is the viscosity of phase α, ϕ is the porosity, k s is the absolute permeability, and ρ α and K ra are the corresponding density and the relative permeability (Blunt et al, ; Falode & Manuel, ), respectively, S α is the saturation and bounded with a constraint of S w + S n = 1, and P n and P w are the macro pressures of the nonwetting and wetting phases, respectively. denotes DCP, which is also critical to analyze the phase behavior in equations –.…”