Relaxation process in an immiscible three-phase flow model

This paper deals with the numerical modeling of immiscible three-phase flows. The main focus here is on the numerical treatment of the source terms of the model. A new scheme based on a more coupled approach than the preexisting fractional step strategy is presented. Properties of this scheme are given. Numerical applications highlight the benefits of this scheme in terms of both accuracy and stability.