Industrial Cz growth of Si crystal of 300 mm and higher diameter usually requires DC magnetic fields (MFs) to suppress turbulence in the melt. We present 3D unsteady analysis of melt turbulent convection in an industrial Cz system coupled with the effect of the transversal MF for different argon gas flow rates for the stage before crystal seeding. We have performed detailed 2D axisymmetric modeling of global heat transfer in the whole Cz furnace. Radiative heat fluxes obtained in 2D modeling have been used in detailed 3D steady and unsteady modeling of crystallization zone. LES method is applied as a predictive approach for modeling of turbulent flow of silicon melt. We have obtained flow structure and temperature distribution in the melt, which were different from previously reported data. We have observed a well-fixed dark spike which includes low temperature melt area on the melt free surface in MF cases. These results indicates that MF and argon flow rate conditions are important to achieve stable positioning of the dark spike on the melt free surface for optimized crystal seeding without uncontrollable meltdown and single crystal structure loss.