Numerical simulation of roll waves in pipelines using the two-fluid model

A finite volume discretization of the incompressible two-fluid model in four-equation form is proposed for simulating roll waves appearing in multiphase pipelines. The new formulation has two important advantages compared to existing roll wave simulators: (i) it is conservative by construction, meaning that the correct shock magnitude is obtained at the hydraulic jump, and (ii) it can be more easily extended with additional physics (e.g. Compressibility, axial diffusion, surface tension), without rederiving the model equations. A simple, robust, first-order upwind discretization of the four-equation model is able to capture the roll wave profiles, although a fine grid is needed to achieve converged results. The four-equation model leads to new roll wave solutions that differ from existing analytical and numerical results. Our solutions are believed to be physically more correct because the shock relations satisfy physically conserved quantities.