For the purpose of optimization by numerical analysis, a highly efficient flotation column named cyclonicstatic micro-bubble flotation apparatus, which is widely used in mineral processing, was simulated with computational fluid dynamics to study the distributions of velocity magnitude, pathlines, turbulent intensity, and turbulent kinetic energy under different discretization schemes, turbulence models, and rotational speeds of pump. This simulation was based on the volume of fluid multiphase and multiple reference frame models. Spatial discretization scheme had an obvious effect on the velocity distribution in the column and revolving flotation sections, especially in the former. The introduction of a turbulence model changed the velocity distribution, the turbulence state, and the pathlines in the column and revolving flotation sections. Moreover, the standard, renormalization group (RNG), and realizable k-epsilon turbulence models had numerous differences. Among these models, the RNG and realizable turbulence models had the most similar results. As the rotational speed of pump increased, velocity and turbulence field intensified. The flow field in the column flotation section differed from that in the revolving flotation section.