Numerical simulations of two-fluid flow models based on the full Navier-Stokes equations are presented. The models include six and seven partial differential equations, namely, six- and seven-equation models. The seven-equation model consists of a non-conservative equation for volume fraction evolution of one of the fluids and two sets of balance equations. Each set describes the motion of the corresponding fluid, which has its own pressure, velocity, and temperature. The closure is achieved by two stiffened gas equations of state. Instantaneous relaxation towards equilibrium is achieved by velocity and pressure relaxation terms. The six-equation model is deduced from the seven-equation model by assuming an infinite rate of velocity relaxation. In this model, a single velocity is used for both fluids. The numerical solutions are obtained by applying the Strang splitting technique. The numerical solutions are examined in a set of one, two, and three dimensions for both the six- and seven-equation models. The results indicate very good agreement with the experimental results. There is an insignificant difference between the results of the two models, but the six-equation model is much more economical compared to the seven-equation model. (C) 2015 Elsevier Ltd. All rights reserved.