Benchmark of different CFL conditions for IMPES

The IMplicit Pressure Explicit Saturation (IMPES) method is a prevalent way to simulate multiphase flows in porous media. The numerical stability of this sequential method implies limitations on the time step, which may depend on the flow regime studied. In this note, three stability criteria related to the IMPES method, that differ in their construction on the observed variables, are compared on homogeneous and heterogeneous configurations for different two-phase flow regimes (viscous/capillary/gravitational). This highlights that there is no single optimal criterion always ensuring stability and efficiency. For capillary dominated flows, the Todd's condition is the most efficient one, while the standard Coat condition should be preferred for viscous flows. When gravity effects are present, Coat's condition must be restricted, but remains more efficient than the Todd's condition.

RésuméL'IMplicit Pressure Explicit Saturation method (IMPES) est l'une des principales méthodes pour traiter les cas d'écoulements multiphasiques en milieu poreux. La stabilité numérique de cette méthode séquentielle implique des contraintes différentes sur le pas de temps selon le régime d'écoulement étudié. Dans cette note, les trois principaux critères de stabilité liés à l'IMPES sont testés sur des milieux homogènes et hétérogènes pour différents régimes (visqueux/capillaire/gravitaire). Cette étude montre qu'aucun critère optimal, réunissant stabilité et efficacité, ne se dégage. Pour les écoulements capillaires, la condition de Todd est la plus efficace, tandis que la condition standard de Coats est préférable pour les écoulements visqueux. Quand les effets gravitaires sont pris en compte, la condition de Coats doit être restreinte, mais demeure plus efficace que celle de Todd.