Two approaches to compute unsteady compressible two-phase flow models with stiff relaxation terms

The paper deals with the numerical modeling of two-phase flows while using Baer-Nunziato type models. Focus is given here on the numerical treatment of source terms that involve three (or four) relaxation time scales. A new coupled approach relying on the continuous analysis of the system of ODEs is compared with a more widely used strategy grounded on the fractional step approach. Properties of schemes are given in both cases. Several numerical applications show that the coupled approach should be preferred for both stability and accuracy reasons