In this work, an improved multiphase computational fluid dynamics (CFD) model is developed for simulating medium segregation in a dense medium cyclone. This model approach uses the modified mixture model with the granular option and Reynolds stress model (RSM) to resolve the turbulent mixing of the particles. Rheological behavior of the medium is considered through various forms of viscosity models such as granular viscosity, Newtonian and non-Newtonian model corrected with particle loading and fraction of ultra-fines. Multiphase simulations using the granular viscosity option although predict the overall slurry volume split, product medium densities better than Newtonian models, but the mixture viscosity values are restricted to substantially lower values close to water viscosity levels. Simulations using the Newtonian viscosity model corrected with particle loading and fine fraction below 53. μm size are predicted the overall medium viscosity levels well above water viscosity levels but underpredicted the cyclone underflow medium density to the experimental data. Multiphase simulations with the non-Newtonian Herschel-Bulkley viscosity model are able to predict the medium segregation close to the gamma ray tomography data, and the predicted medium viscosity levels are increased up to 18 cp. Overall predicted product medium densities, underflow volume split and radial medium segregation levels are compared to the experimental data. Lagrangian particle tracking (LPT) is super imposed on converged medium simulations used for coal partitioning predictions. This model can be further used in developing new designs in the future. © 2015 Elsevier B.V.