COMMIX-M, a three-dimensional transient and steady-state computer program written in Cartesian and cylindrical coordinates, has been developed by Argonne National Laboratory. This computer program is capable of analyzing multiphase flow and heat transfer and utilizes the separate phases model wherein each phase has its own mass, momentum, and energy equations. This computer program is in its early stages of development for application to test various interphase interaction models and to predict design and processing of dense fluid-solids suspension systems. COMMIX-M contains preliminary constitutive relationships for interfacial drag, solids viscosities and stresses to describe the solids rheology, and shear lift forces from the literature. Also included is a solids partial slip boundary condition to allow non-zero tangential velocity at the tube walls. Analyses of some of the steady-state, fully-developed isothermal carrier fluid velocity and solids concentration data of Altobelli et al. and Sinton and Chow are presented. These experimental data obtained using three-dimensional time-of-flight nuclear magnetic (NMR) imaging techniques were carefully performed and represent some of the best available open literature data of their kind. NMR imaging offers powerful techniques to non-intrusively determine three-dimensional time-dependent velocity and concentration fields to assist development and validation of the constitutive models and the computer programs describing concentrated suspensions. Analyses of these NMR data, together with comparisons of computed and measured concentration and velocity profiles provide some insights into the mechanisms governing the observed phenomena. Recommendations for future research are given. To the authors` knowledge, these are the first such comparisons of theory and experiment