Implicit temporal discretization is widely applied in reservoir simulation because of its
unconditional numerical stability with respect to time-step size. At each time-step, implicit …