##
**On iterative IMPES formulation for two phase flow with capillarity in heterogeneous porous media.**
*(English)*
Zbl 1253.76124

Summary: This work is a continuation of [J. Kou and S. Sun, Comput. Fluids 39, No. 10, 1923–1931 (2010; Zbl 1245.76147)] where we presented an efficient improvement on the implicit pressure explicit saturation (IMPES) method for two-phase immiscible fluid flow in porous media with different capillarity pressures. In the previous work, we present an implicit treatment of capillary pressure appearing in the pressure equation. A linear approximation of capillary function is used to couple the implicit saturation equation into the pressure equation that is solved implicitly. In this paper, we present an iterative version of this method. It is well-known that the fully implicit scheme has unconditional stability. The new method can be used for solving the coupled system of nonlinear equations arisen after the fully implicit scheme.

We follow the idea of the previous work, and use the linear approximation of capillary function at the current iteration. This is different from iterative IMPES that computes capillary pressure by the saturations at the previous iteration. From this approximation, we couple the saturation equation into the pressure equation, and establish the coupling relation between the pressure and saturation. We employ the relaxation technique to control the convergence of the new method, and we give a choice of relaxation factor. The convergence theorem of our method is established under the natural conditions. Numerical examples are provided to demonstrate the performance of our approach, and the results show that our method is efficient and stable.

We follow the idea of the previous work, and use the linear approximation of capillary function at the current iteration. This is different from iterative IMPES that computes capillary pressure by the saturations at the previous iteration. From this approximation, we couple the saturation equation into the pressure equation, and establish the coupling relation between the pressure and saturation. We employ the relaxation technique to control the convergence of the new method, and we give a choice of relaxation factor. The convergence theorem of our method is established under the natural conditions. Numerical examples are provided to demonstrate the performance of our approach, and the results show that our method is efficient and stable.

### MSC:

76S05 | Flows in porous media; filtration; seepage |

65M06 | Finite difference methods for initial value and initial-boundary value problems involving PDEs |

65M12 | Stability and convergence of numerical methods for initial value and initial-boundary value problems involving PDEs |

76M20 | Finite difference methods applied to problems in fluid mechanics |

76T20 | Suspensions |