Numerical resolution for a compressible two-phase flow model based on the theory of thermodynamically compatible systems

This paper is concerned with the numerical solution of the equations governing two-phase gas-solid mixture in the framework of thermodynamically compatible systems theory. The equations constitute a non-homogeneous system of nonlinear hyperbolic conservation laws. A total variation diminishing (TVD) slope limiter centre (SLIC) numerical scheme, based on the splitting approach, is presented and applied for the solution of the initial-boundary value problem for the equations. The model equations and the numerical methods are systematically assessed through a series of numerical test cases. Strong numerical evidence shows that the model and the methods are accurate, robust and conservative. The model correctly describes the formations of shocks and rarefactions in two-phase gas-solid flow.