In this paper a general drift-flux model describing a subsonic and isentropic multi-phase fluid in connected pipes is considered. Each phase is assumed to be isentropic with its own sonic speed. The components are gamma-law gases with γ > 1. For such the computational challenge at a junction is the computation of rarefaction waves which do not have a readily available analytical form. Firstly, ...