Solving stiff mass transfers in compositional multiphase flow models: Numerical stability and spurious solutions

Author: Tardy Philippe   Quintard Michel  

Publisher: Springer Publishing Company

ISSN: 1420-0597

Source: Computational Geosciences, Vol.3, Iss.2, 1999-10, pp. : 161-183

Disclaimer: Any content in publications that violate the sovereignty, the constitution or regulations of the PRC is not accepted or approved by CNPIEC.

Previous Menu Next

Abstract

This paper points out two numerical problems linked to the resolution of compositional multiphase flow models for porous media with the finite-volume technique. In particular, we consider fluid mixtures featuring fast mass transfers between the phases, hence stiff. In this context, we show how the computation of mass exchange kinetics can be expensive and that erroneous saturation front locations arise. A numerical splitting method is developed which is proven to be stable with advection-type time steps, whatever the stiffness of the mass transfer. The erroneous front location problem is illustrated and shown to be intrinsically linked to the numerical diffusion. If we assume that the fluids are in thermodynamical equilibrium, we find that spurious solutions can be avoided by deriving and solving a new uncoupled hyperbolic equation for the saturation.