A set of global heat transfer simulations in Czochralski (Cz) furnace for producing silicon single crystals have been performed to investigate the effect of crystal and crucible rotations on melt convection and crystal/melt interface shape. The 2D axisymmetric simulations are carried out taking into account radiative and conductive heat transfers between furnace components, melt convection including thermocapillary forces (Marangoni) and gas flow. Melt flow pattern and temperature distribution have been studied for several combinations of crystal and crucible rotations and for different crystal heights. The result shows that crystal/melt interface shape and melt flow regime are strongly sensitive to the rate of rotations of both crystal and crucible. Interfaces with low deflection can be achieved for certain combinations of crystal and crucible rotation rates. (C) 2010 Elsevier B.V. All rights reserved.