Numerical approximation for a Baer-Nunziato model of two-phase flows