Two-dimensional (2D) transient numerical simulations are performed to investigate the evolution of the thermal and flow fields during the growth of multi-crystalline silicon ingots with two different silicon feedstock capacities, 800 kg and 1600 kg. The simulation results show that there are differences in the structure of the melt flow. In the 1600 kg case, there is a reduction in the concavity of the crystal-melt interface near the crucible wall and an increase in the convexity of the interface at higher solidification fractions. Moreover, the Voronkov ratios, which are indicative of the formation of defects, become lower during the solidification process.