For understanding the monosodium aluminate hydrate crystallization from the supersaturated aluminate solution containing red mud as the leaching liquor of bauxite, the liquid-solid-solid dispersion of a simulant system, i.e. glycerite, red mud and sand, in a stirred reactor has been experimentally investigated as well as simulated using computational fluid dynamics model (CFD) for the first time. The computational model is based on the Eulerian multi-fluid model along with RNG k-epsilon turbulence model, where Syamlal-obrien-symmetric drag force model (Syamlal, 1987) of the interphase momentum transfer between two dispersed solid phases is taken into account. A good agreement is obtained between the experimental data of solid distributions and the simulation results in the flow fields of liquid-solid-solid as well as liquid-solid systems. The solid suspension qualities of both liquid-solid and liquid-solid-solid systems in the stirred reactors with and without draft tube were also studied in detail based on mixing time, the standard deviation of solid concentration proposed by Bohnet and Niesmak (1980), the flow pattern and power number. The influence of the interaction between two dispersed solid phases on the suspension of red mud is found significantly greater than that of sand. The holdup of sand below the impeller is considerably larger than that above the impeller and the red mud dispersion approaches homogeneous in the reactor. The mixing time of liquid-solid-solid suspension is longer than that of liquid-solid suspension under the same conditions, and the mixing times of both systems in the stirred reactor with draft tube are longer than that in the reactor without draft tube. Furthermore, the distributions of sand and red mud in the reactor with draft tube were found less homogeneous than those without draft tube in most cases. (C) 2010 Elsevier Ltd. All rights reserved.