Turbulent characteristics of Czochralski melt flow are presented using the unsteady Reynolds-averaged Navier–Stokes (URANS) turbulence modeling approach. Three-dimensional, transient computations were performed using the Launder and Sharma low-Re k-ε model and Menter shear stress transport (SST) k-ω model on an idealized Czochralski setup with counterrotating crystal and crucible. A comparative assessment is performed between these two Reynolds-averaged Navier–Stokes (RANS) models in capturing turbulent thermal and flow behaviors. It was observed that the SST k-ω model predicted a better resolution of the Czochralski melt flow capturing the near wall thermal gradients, resolving stronger convective flow at the melt free surface, deciphering more number of characteristics Czochralski recirculating cells along with predicting large number of coherent eddy structures and vortex cores distributed in the melt and hence a larger level of turbulent intensity in the Czochralski melt compared with that by Launder and Sharma low-Re k-ε model.