Implicit linearization scheme for nonstandard two-phase flow in porous media

In this article, we consider a nonlocal (in time) two-phase flow model. The nonlocality is introduced through the wettability alteration induced dynamic capillary pressure function. We present a monotone fixed-point iterative linearization scheme for the resulting nonstandard model. The scheme treats the dynamic capillary pressure functions semiimplicitly and introduces anL-scheme type stabilization term in the pressure as well as the transport equations. We prove the convergence of the proposed scheme theoretically under physically acceptable assumptions, and verify the theoretical analysis with numerical simulations. The scheme is implemented and tested for a variety of reservoir heterogeneities in addition to the dynamic change of the capillary pressure function. The proposed scheme satisfies the predefined stopping criterion within a few number of iterations. We also compared the performance of the proposed scheme against the iterative implicit pressure explicit saturation scheme.